Skip to main content

Genome-wide transcriptome analysis of hypothalamus in rats with inherited stress-induced arterial hypertension

Abstract

Background

The hypothalamus has an important role in the onset and maintenance of hypertension and stress responses. Rats with inherited stress-induced arterial hypertension (ISIAH), reproducing the human stress-sensitive hypertensive state with predominant involvement of the neuroendocrine hypothalamic-pituitary-adrenal and sympathoadrenal axes, were used for analysis of the hypothalamus transcriptome.

Results

RNA-seq analysis revealed 139 genes differentially expressed in the hypothalami of hypertensive ISIAH and normotensive Wistar Albino Glaxo (WAG) rats. According to the annotation in databases, 18 of the differentially expressed genes (DEGs) were associated with arterial hypertension. The Gene Ontology (GO) functional annotation showed that these genes were related to different biological processes that may contribute to the hypertension development in the ISIAH rats. The most significantly affected processes were the following: regulation of hormone levels, immune system process, regulation of response to stimulus, blood circulation, response to stress, response to hormone stimulus, transport, metabolic processes, and endocrine system development. The most significantly affected metabolic pathways were those associated with the function of the immune system and cell adhesion molecules and the metabolism of retinol and arachidonic acid. Of the top 40 DEGs making the greatest contribution to the interstrain differences, there were 3 genes (Ephx2, Cst3 and Ltbp2) associated with hypertension that were considered to be suitable for further studies as potential targets for the stress-sensitive hypertension therapy. Seven DEGs were found to be common between hypothalamic transcriptomes of ISIAH rats and Schlager mice with established neurogenic hypertension.

Conclusions

The results of this study revealed multiple DEGs and possible mechanisms specifying the hypothalamic function in the hypertensive ISIAH rats. These results provide a basis for further investigation of the signalling mechanisms that affect hypothalamic output related to stress-sensitive hypertension development.

Background

Brain regulatory mechanisms play a major role in the control of sympathetic outflow and blood pressure (BP) regulation [1]. Based on the review of many studies, it was noted that most forms of hypertension are associated with a wide variety of functional changes in the hypothalamus [2]. In recent studies, progress has been made in understanding neurogenic aspects by the determination of global alterations in gene expression in key brain regions of animal models of neurogenic hypertension [3]. However, ongoing novel approaches and studies on different animal models are expected to lead to the identification of genes that serve as a common link between different forms of the disease [4] and a better understanding of the processes responsible for the increased sympathetic outflow and hypertension development in humans.

The aetiology of essential hypertension is multifactorial. As the studies on experimental animal models provide valuable information to elucidate the nature of polygenic traits [5], a number of animal models for essential hypertension are widely used. One of these is the ISIAH (Inherited Stress-Induced Arterial Hypertension) rat strain, which has been developed to study the mechanisms of stress-sensitive hypertension and its complications. The ISIAH rat strain was selected from outbred normotensive Wistar rats for an increased BP response to emotional stress caused by 30 min restraint in a cylindrical wire-mesh cell. The current population of ISIAH rats is characterized by elevation of both the basal arterial BP, which reaches up to 175.0 ± 3.5 mmHg in males and 165.0 ± 3.0 mmHg in females, and the stress-induced BP, which increases just after the mild emotional stress up to 195.0 ± 2.4 mmHg in males and 174.0 ± 3.2 mmHg in females [6], [7]. The high genetic homogeneity of the ISIAH strain was demonstrated by a DNA fingerprinting procedure using a multilocus microsatellite (CAC)5 probe [8]. Earlier studies showed that the ISIAH rats may be regarded as a model for the human stress sensitive hypertensive disease with predominant involvement of the neuroendocrine hypothalamic-pituitary-adrenal (HPA) and sympathoadrenal systems in the pathogenesis of the hypertensive state [9].

The recently developed next-generation sequencing technologies, widely used for transcriptome analysis (RNA-Seq approach), are making substantial contributions to the understanding of genome-wide expression and regulation [10]. In the current work, RNA-Seq technology was used for comparative analysis of the hypothalamic transcriptome in hypertensive ISIAH and inbred normotensive WAG rats. The goal of the study was to identify the differentially expressed genes (DEGs) and pathways involved in the differences of hypothalamic functions in the hypertensive ISIAH and normotensive WAG rats.

The study revealed multiple DEGs in the hypothalami of hypertensive ISIAH and normotensive WAG rats, including 18 DEGs known to be related to hypertension and the regulation of BP. These DEGs were associated with a diverse group of biological processes and pathways. Three DEGs were in the list of the top 40 DEGs defined as genes making the largest contribution to inter-strain differences. These DEGs may be considered potential candidates for further studies to better understand the mechanisms of hypertension development in ISIAH rats.

Results

Altogether, 11,369 annotated genes were defined as being expressed in the hypothalami of ISIAH and WAG rats and were used in comparative expression analysis, which revealed 139 DEGs. For a complete listing of these results, see Additional file 1. The hierarchical clustering using Euclidean distance is shown in Fig. 1. About half of the DEGs (71 genes, i.e., 51.1 %) were down-regulated in ISIAH rats. Two of the DEGs, Ahsp (alpha haemoglobin stabilizing protein, 1.9 FPKM) and RGD1359290 (ribosomal_L22 domain containing protein, 46.5 FPKM), were detected in the hypothalamus of ISIAH rats but not in WAG rats.

Fig. 1
figure1

Heatmap of the differentially expressed genes in the hypothalami of the ISIAH and WAG rats. Dendrograms were constructed using hierarchical ‘complete linkage’ clustering by Euclidean distance. Yellow depicts genes that were up-regulated in a sample, and red depicts those that were down-regulated

Functional annotation of DEGs

The most significant GO terms for biological processes, which may be related to the hypertensive state of ISIAH rats, and genes in related groups of DEGs are represented in Additional file 2. The major DEG groups were associated with the regulation of hormone levels, immune system processes, blood circulation, response to stimulus and regulation of the response to stimulus, transport, metabolic process, gland development, and regulation of the immune response. Some of these groups consisted of several subgroups describing the specificity of the processes. For instance, the group of DEGs associated with blood circulation included 5 genes that participate in the regulation of BP. The group labelled, ‘response to stimulus’, consisted of several subgroups based on the specificity of the stimuli. These subgroups predominantly consisted of ‘response to stress’, ‘response to hormone stimulus’, ‘response to inorganic substance’, and ‘response to external stimulus’. Many DEGs were associated with different metabolic processes including lipid metabolism, hormone metabolism, and vitamin A and catecholamine metabolic processes. These and some other DEG groups related to immune system processes, oxidation reduction, endocrine system development, cell activation and proliferation, indicated that there were multiple impairments in hypothalamic functions in hypertensive rats.

Two DEGs known to be key players in corticosteroid biosynthesis and BP regulation (Cyp11b1 and Cyp11b2) were estimated to be up-regulated in the hypothalamus of hypertensive ISIAH rats. However, these genes were merged during the Cufflinks assembly and thus were tested as a single unit by the Cuffdiff program. To distinguish between their levels of expression, we performed quantitative real time polymerase chain reaction (qPCR). The results of qPCR confirmed the enhanced expression of both the Cyp11b1 and Cyp11b2 genes in the hypothalamus of ISIAH rats (Fig. 2).

Fig. 2
figure2

The relative mRNA abundance measured by qPCR. The relative mRNA abundance was calculated as the ratio of normalized (mRNA/Rpl30) mRNA level in experimental ISIAH samples to normalized (mRNA/Rpl30) mRNA level in control WAG samples. The normalized mRNA level in control samples of the WAG rats was assigned a value of 1. Vertical bars show the standard error of the mean, and significance of inter-strain difference is indicated by *p < 0.05, **p < 0.01, ***p < 0.001

The analysis in the Kyoto Encyclopedia of Genes and Genomes Pathway Database (KEGG) showed that the inter-strain differences in gene expression are related to six pathways that were significantly enriched (p < 0.05; Additional file 3). Four of these pathways were associated with the function of the immune system. All of the pathways (except for retinol metabolism) contained genes associated with hypertension and with central nervous system (CNS) diseases.

Genes associated with hypertension and CNS diseases

Of the genes found to be differentially expressed in the hypothalami of ISIAH and WAG rats, 16 genes are annotated in the Rat Genome Database (RGD) as associated with hypertension, and two additional DEGs (P2rx4, and Sult1a1) were related to the regulation of BP according to the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool. These 18 genes might be essential for hypertension development in ISIAH rats. Most of these genes are annotated in RGD as associated with diseases concomitant to hypertension, such as insulin resistance and CNS diseases (Table 1). Altogether, we found 26 DEGs referred to in RGD as associated with CNS diseases, including cerebrovascular disorders, brain infarction, and brain ischemia. Many of these genes were also related to immune system diseases (Table 2).

Table 1 Genes differentially expressed in ISIAH and WAG hypothalami and referred to in Databases as associated with hypertension and blood pressure regulation
Table 2 Genes differentially expressed in ISIAH and WAG hypothalami and referred to in Rat Genome Database as associated with central nervous system diseases

Transcription factors

Five transcription factor genes were differentially expressed in ISIAH and WAG hypothalami (Table 3). One of these (Ercc2) is annotated in RGD as being associated with hypertension. According to the GO functional annotation, two genes (Ercc2 and Grhl3) were associated with response to stress, and three genes (Foxg1, Mcm7, and Pitx1) were related to metabolic processes. Moreover, Pitx1 was related to the GO term ‘endocrine system development’ (Additional file 2).

Table 3 Transcription factor genes differentially expressed in ISIAH and WAG hypothalami

Genes making the greatest contribution to the inter-strain differences

The partial-least squares discriminant analysis (PLS-DA) was used to detect the genes that made the greatest contribution to the inter-strain differences. The PLS-DA Axes maximizing the distances between ISIAH and WAG rats were constructed (Additional file 4), and the correlation between gene expression and PLS-DA Axis 1 was defined. The genes that were designated as differentially expressed are shown in red in Figure in Additional file 5, and their polar position along the first PLS-DA Axis confirms their contribution to the existing inter-strain variations.

The top 40 DEGs that make the most significant contribution to the inter-strain differences possess different biological potencies (Table 4). According to RGD annotations, there were several genes related to CNS diseases (Abcg2, Ephx2, RGD1563482, Tgfb2, Mal, and Cst3), and 3 genes associated with hypertension (Ephx2, Cst3, and Ltbp2) are found on this list. These 3 genes may be considered to be potential candidates for further studies to better understand the mechanisms of hypertension development in ISIAH rats. Their differential transcription was validated by qPCR (Fig. 2). The comparison of the relative mRNA abundance is shown in Additional file 6. The correlation coefficient between the results from the two methods was 0.98.

Table 4 The top 40 differentially expressed genes making the greatest contribution to the inter-strain differences in hypothalamic function of ISIAH and WAG rats

Discussion

Mapping and quantifying mammalian transcriptomes by deep sequencing (RNA-Seq) provides a digital measure of the presence and prevalence of transcripts [11]. In the current study, an RNA-Seq approach was performed to detect genes with altered transcriptional activity in the hypothalamus of hypertensive ISIAH rats compared to normotensive WAG rats. However, only 18 out of the 139 DEGs found in the current study were characterized as being associated with the hypertension. This could be explained by genetic drift during the selection process, which could result in randomly fixed alleles contributing to differences in gene expression but not related to disease itself. The functional annotation of DEGs helped to define the main biological processes that might contribute to inter-strain differences. This analysis showed that almost all GO groups contained the DEGs associated with hypertension (see Additional file 2). These results highlight once again that many genes and biological processes may contribute to the stress-sensitive hypertension development.

As expected, GO annotation revealed multiple DEGs related to the hormonal nature of hypothalamic function. The group of eight DEGs participating in the regulation of hormone levels (see Additional file 2) might play an important role in the orchestration of the changes in physiological and metabolic processes in ISIAH hypothalamus. This group contained several genes associated with hypertension. Two genes in this group, cytochrome P450, family 11, subfamily b, polypeptide 1 (Cyp11b1) and cytochrome P450, family 11, subfamily b, polypeptide 2 (Cyp11b2), are known to be key players in steroidogenic pathways [12]. These genes are particularly known to be involved in corticosterone and aldosterone biosynthesis and are widely known to be regulators of BP. It was shown that the small amounts of aldosterone synthesized in the brain of Dahl salt-sensitive hypertensive rats could provide a local ligand for autocrine or paracrine activation of the mineralocorticoid receptor [13], and this was important in the genesis of salt-sensitive hypertension [14]. Aldosterone synthesized in the brain contributes to salt-sensitive hypertension through an increased salt appetite, sympathetic drive and vasopressin release [15]. The up-regulation of the Cyp11b2 gene found in the current study may result in aldosterone excess in the hypothalamus of ISIAH rats and may also contribute to the development of stress-sensitive hypertension.

Several other DEGs participated in the regulation of hormone levels were related to retinol metabolism: aldehyde dehydrogenase 1 family, member A2 (Aldh1a2), cytochrome P450, family 26, subfamily b, polypeptide 1 (Cyp26b1), retinol binding protein 4, plasma (Rbp4), and retinol saturase (Retsat). Retinoic acid (RA) is an essential signalling molecule. In the adult brain, RA has been shown to be involved in synaptic plasticity, the induction of neural differentiation, motor axon outgrowth and nerve regeneration [16], [17]. In the current study, Rbp4 and RetSat were in a list of the top 40 DEGs that contribute the most to inter-strain differences (Table 4). Additionally, in this list, there was one additional gene, stimulated by retinoic acid 6 (Stra6), related to retinol transport across the cell membrane [18]. These data suggest that retinol metabolism is likely to make a substantial contribution to inter-strain differences in hypothalamic functions; however, the direct influence of this pathway to BP regulation has not been defined so far.

One additional gene associated with the regulation of hormone level was Cga (glycoprotein hormones, alpha polypeptide). It encodes the alpha subunit of several glycoprotein hormones: chorionic gonadotropin, luteinizing hormone, follicle stimulating hormone, and thyroid stimulating hormone. Its up-regulation in the ISIAH hypothalamus may contribute to a variety of biological processes related to endocrine system functioning in hypertensive rats.

The Sult1a1 (sulfotransferase family, cytosolic, 1A, phenol-preferring, member 1) gene was associated with several biological processes including the regulation of hormone levels, the regulation of BP, and the response to hormone stimulus (Additional file 2). Rat sulfotransferase Sult1a1 catalyzes the sulfation of neurotransmitters and xenobiotic phenolic compounds [19]. Sulfation usually leads to the inactivation of biological signalling molecules including catecholamines [20]. It was shown that the experimental reduction of norepinephrine content in the anterior hypothalamus induced an increase in arterial blood pressure and heart rate [21]. As Sult1a1 sulfates catecholamines, the down-regulation of this gene in the ISIAH hypothalamus may be adaptive and directed against the excessive elevation of BP.

Two other genes related to the regulation of the BP level that are differentially expressed in ISIAH and WAG hypothalami were the purinergic receptor P2X, the ligand-gated ion channel 4 (P2rx4) and epoxide hydrolase 2, cytoplasmic (Ephx2). The product of the P2rx4 gene belongs to the family of purinoceptors for adenosine triphosphate. This receptor functions as a ligand-gated cation-selective channel. Endothelial P2X4 channels regulate BP and vascular remodelling. P2rx4 has a key role in the response of endothelial cells to changes in blood flow. P2rx4(−/−) mice have higher BP levels and excrete smaller amounts of NO products in their urine than wild-type mice [22]. Based on this, we suggest that down-regulation of the P2rx4 gene in the ISIAH hypothalamus might contribute to hypertension development.

Ephx2 encodes the soluble epoxide hydrolase (sEH) that metabolizes the epoxyeicosatrienoic acids (EETs). EETs modulate ion transport and gene expression and produce vasorelaxation, as well as anti-inflammatory and pro-fibrinolytic effects [23]. sEH inhibition lowers arterial BP in angiotensin II-induced hypertension [24]. sEH overexpression was linked to hypertension in the spontaneously hypertensive rats [25]. Ephx2 was also considered to be one of the gatekeeper genes that contributes to programmed hypertension [4]. In our study, Ephx2 was one of several DEGs related to arachidonic acid (AA) metabolism. The two other DEGs, cytochrome P450, family 2, subfamily j, polypeptide 10 (Cyp2j10) and cytochrome P450, family 4, subfamily x, polypeptide 1 (Cyp4x1), are also known to encode enzymes that regulate EET production. The down-regulation of these two genes and the enhanced expression of the Ephx2 gene in the ISIAH hypothalamus may lead to a reduced amount of EETs and contribute to hypertension development in ISIAH rats. In the current study, Ephx2 and Cyp2j10 were in a list of the top 40 DEGs that contribute the most to inter-strain differences (Table 4).

One more DEG related to AA metabolism was prostaglandin D2 synthase (brain) (Ptgds). The decreased expression of Ptgds found in the current study may lead to a reduced amount of prostaglandin D2, which is a major brain prostaglandin that functions as a neuromodulator and/or trophic factor in the CNS [26]. It is involved in smooth muscle contraction/relaxation and is a potent inhibitor of platelet aggregation [27].

Both GO annotation and KEGG analysis revealed highly significant differences in the expression of genes belonging to immune system processes in the hypothalami of ISIAH and WAG rats. It is long known that immune system changes play a role in hypertension and that extensive bidirectional interactions between the sympathetic nervous system and the immune system exist [28], [29]. Recent studies have shown that both innate and adaptive immunity contribute to hypertension [30]. Major histocompatibility complex (MHC) class I molecules are ligands for the killer-cell immunoglobulin-like receptors that are expressed by natural killer (NK) cells and T cells. The interactions between these molecules contribute to both innate and adaptive immunity [31]. MHC class-II molecules are key participants in immune activation events in autoimmunity [32]. It was shown that mice lacking adaptive immune cells, including recombinase-activating gene-deficient mice and rats and mice with severe combined immunodeficiency, have blunted hypertension responses to stimuli such as ANG II, high salt, and norepinephrine [33]. MHC molecules are rare in healthy brain tissue, but are heavily expressed on microglial cells after inflammatory or neurodegenerative processes [34]. The transcription of most of the DEGs that encode MHC molecules was decreased in the hypothalamus of ISIAH rats. However, two genes, RT1-A1 (RT1 class Ia, locus A1) and RT1-M3-1 (RT1 class Ib, locus M3, gene 1), were up-regulated. Earlier, it was shown that the expression of the RT1-A1 molecule protects targets from NK lysis [35]. In another study it was found that the expression of RT1-A1 became up-regulated in glial cells after mechanical nerve injury [36]. We can suppose that the elevated transcription of RT1-A1 in the ISIAH hypothalamus may play a protective role.

Altogether, there were three genes from the list of the top 40 DEGs, RT1 class Ia, locus A2 (RT1-A2), RT1 class Ib, locus S3 (RT1-S3), and transforming growth factor, beta 2 (Tgfb2), that made the most significant contribution to the inter-strain differences in the group ‘immune system process’ (Table 4). The possible role of the MHC molecules in pathology development was discussed above. TGFb2 is known to promote the survival of dopaminergic neurons [37]. The hypothalamic changes in neurotransmitter and neuromodulator functions via the sympathetic nervous system appear to play a predominant role in most forms of hypertension [2]. Tgfb2 was one of several DEGs identified in the current study that were related to endocrine system development (see Additional file 2). Two other genes related to this group and not discussed above were solute carrier family 6 (neurotransmitter transporter), member 3 (Slc6a3) and paired-like homeodomain 1 (Pitx1). The Slc6a3 gene encodes the dopamine transporter. Its primary physiological role is the regulation of extracellular dopamine by means of rapid transport (reuptake) of dopamine back into the dopaminergic presynaptic terminals. Increased expression of the dopamine transporter leads to a loss of dopamine neurons and oxidative stress [38]. Therefore, the up-regulation of Tgfb2 and the decreased expression of Slc6a3 found in the hypothalamus of ISIAH rats might have a protective role for dopaminergic neurons.

Pitx1 encodes a homeobox transcription factor involved in transcriptional activation of the pro-opiomelanocortin (Pomc) gene [39], [40]. The increased hypothalamic POMC presumably attenuates weight gain [41], activates sympathetic neurotransmission to adipose tissues, accelerates energy expenditure [42], and could be, at least in part, responsible for the increased sensitivity to insulin [43]. However, Pomc transcription was not significantly increased in the ISIAH hypothalamus in our study. Therefore, the up-regulation of Pitx1, may be essential for the other signalling pathways in the hypothalamus of hypertensive rats.

The hypothalamus is known as a brain region involved in the regulation of body homeostasis and stress responses. In the current study, multiple DEGs were associated with the response to different stimuli such as stress, hormonal stimulus, inorganic substance, and external stimuli (Additional file 2). As the study was conducted on rats not exposed to the stress, the DEGs in the group ‘response to stress’ might be essential for stress-sensitive hypertension development. Altogether, we found nine genes (Cst3, Cyp11b2, Ephx2, Fn1, Igfbp2, P2rx4, Prkcd, RT1-Ba, and Sult1a1) annotated in the RGD as associated with hypertension in this group. Two of these genes cystatin C (Cst3) and epoxide hydrolase 2, cytoplasmic (Ephx2) were also present in the list of top 40 DEGs defined as making the largest contribution to inter-strain differences (Table 4). Cystatin C (Cst3 gene) is an inhibitor of cysteine proteinases. Expression of cystatin C has protective effects against various oxidative stresses that induce cell death [44]. The studies of hypothalamic mechanisms underlying the neurogenic hypertension in mice showed that inflammation and oxidative stress via reactive oxygen species generation could contribute to oxidative damage and activation of the sympathetic nervous system leading to an increase in BP [45]. Therefore, the decreased transcription of the Cst3 gene in the ISIAH hypothalamus may contribute to hypertension development as well as excessive transcription of the Ephx2 gene discussed above.

Altogether, the group ‘response to stimulus’ contained 13 genes (Abcg2, Cst3, Ephx2, Grhl3, Grm2, Gstm4, Lcat, Lyz2, Mx2, Rbp4, RT1-A2, RT1-S3, and Tgfb2) included in the list of top 40 DEGs that contribute the most to inter-strain differences in hypothalamic function in the current study. The known activities of genes described in Table 4 or discussed above underscore the diversity of the mechanisms that contribute to the inter-strain differences. These genes might be useful in further studies to elucidate the mechanisms accompanying the stress-sensitive hypertension development.

In the current study, only 3 genes associated with hypertension (Ephx2, Cst3 and Ltbp2) appeared on the list of the top 40 DEGs that make the largest contribution to inter-strain differences. We suggest that these 3 genes might be of interest for their further studies as potential therapeutic targets for stress-sensitive hypertension therapy. One of these genes, Ephx2, was already proposed as a target for the treatment of hypertension [46]. Ltbp2 had previously been proposed to be a potential biomarker for atrial stress [47] and was patented as a biomarker for cardiovascular, hematological, neurological, endocrinological, and urological diseases, as well as cancer [48]. As discussed above, cystatin C has protective effects against various oxidative stresses. Therefore, its contribution to BP regulation may be expected.

The analysis of different animal models of hypertension helps to identify the genes which may serve as a common link between different forms of the disease [4]. In the current study, we found several common genes with changed transcriptional activity in the hypothalami of hypertensive ISIAH rats and hypertensive Schlager mice that represent a genetic model of neurogenic hypertension (Table 5) [45]. One of these genes, chitinase 3-like 1 (cartilage glycoprotein-39) (Chi3l1), is known to be associated with hypertension. The protein encoded by this gene is thought to play a role in the process of inflammation [49], tissue remodelling [50] and reactive oxygen species production [51]. The other genes listed in Table 5 were not previously related to the regulation of blood pressure, but, as they are common between two models of the disease, these genes may be considered promising targets for further studies directed to understanding the processes contributing to increased sympathetic outflow and hypertension development. During the last decade, the results of high-throughput biological experiments made a significant contribution to understanding the molecular bases of complex diseases. RNA-Seq is one of the recently developed approaches to transcriptome profiling which provides a far more precise measurement of levels of transcripts than other methods [52]. However, it has some limitations. The main limitation of RNA-Seq analysis is its expensiveness. The other limitation is the relative impossibility of technical validation of all the multiple DEGs found in the experiment. Therefore, the current study was conducted with the use of only three biological replicates in both experimental and control groups, and the expression of only several DEGs was validated technically. However, it was shown earlier that the results from RNA-Seq analysis showed high levels of reproducibility for both technical and biological replicates [53], [54]. In the current study, we found a strong correlation (r = 0.98) between the qPCR and RNA-Seq data. These results are similar to those reported by other groups [53].

Table 5 The list of the common genes with changed transcriptional activity in hypothalami of hypertensive ISIAH rats and hypertensive Schlager mice (26 weeks old), which represent a genetic model of neurogenic hypertension

Conclusion

The results of comparative transcriptional profiling of the hypothalamus in hypertensive ISIAH and normotensive WAG rats helped to identify multiple DEGs, including 18 genes associated with hypertension and regulation of BP. These DEGs were related to the diversity of biological processes and pathways contributing significantly to the inter-strain differences related to hypothalamic function. The discussion demonstrated that many genes found to be differentially expressed in ISIAH and WAG hypothalami might be responsible for pathology development, and that the changed expression of other genes may exert compensatory effects and be directed toward the restoration of homeostasis. As the number of hypertensive genes was considerable, we were not able to discuss all of them in details. Moreover, the discussion of some genes was complicated by uncertainty about their roles, especially within the hypothalamus. In the current study, we defined 3 genes associated with hypertension (Ephx2, Cst3 and Ltbp2) that might be of interest for further studies as potential therapeutic targets for stress-sensitive hypertension therapies. Our findings provide a basis for identification of potential biomarkers of stress-sensitive hypertension and for further investigation of the signalling mechanisms that affect hypothalamic output related to hypertension development and the resulting sympathetic nervous system effects on other tissues and the organism as a whole.

Methods

Animals

Two inbred rat strains were used in the current study. Both hypertensive ISIAH (Inherited Stress Induced Arterial Hypertension) and normotensive Wistar Albino Glaxo (WAG/GSto-Icgn) rats were bred in the Center for Genetic Resources of Laboratory Animals at the Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences (Novosibirsk, Russia, RFMEFI61914X0005 and RFMEFI61914X0010). All rats were maintained under standard conditions with free access to food and water. One week before the systolic arterial BP measurement, animals were placed in individual cages. The systolic arterial BP was measured indirectly by the tail-cuff method under short-term ether anaesthesia to exclude the effect of psychological stress induced by the measuring procedure. In RNA-seq experiments, the 3-month-old ISIAH (n = 3) and WAG (n = 3) males were used. Their systolic arterial BP was 171.7 ± 1.22 mmHg in ISIAH and 116.33 ± 1.86 mmHg in WAG males. The rats were decapitated 6 days after BP measurement, and their hypothalami were isolated according to known anatomical location [55]. All samples were stored in RNALater (Qiagen, Chatsworth, CA) at −70 °C until use. All animal experiments were approved by the Institute’s Animal Care and Use Committee.

RNA-Seq analysis

The collected samples were sent to JSC Genoanalytica (Moscow, Russia) where mRNA was extracted using Dynabeads mRNA Purification Kit (Ambion, USA). cDNA libraries were constructed using NEBNext mRNA Library Prep Reagent Set for Illumina (NEB, USA) following the manufacturer’s protocol and subjected to single-end sequencing on Illumina HiSeq1500 with read length of 50 bases. Three hypothalami from ISIAH and three hypothalami from WAG rats were run as experimental replicates. The resulting fastq-formatted sequencing data were mapped to the RGSC Rnor_5.0\rn5 reference genome using Tophat2 aligner [56] and ENSEMBL/RefSeq gene annotation. A quality assessment of the mapped data was performed using the module ‘CollectRnaSeqMetrics’ from Picard tools suite (http://broadinstitute.github.io/picard/). The summary statistics for each sequenced library is given in Additional file 7. The Cufflinks/Cuffdiff programs were then used to estimate gene expression levels in FPKM (fragments per kilobase of transcript per million mapped reads) and to perform differential expression analysis [57]. Genes were defined as being expressed if they were assigned to test status “OK” (i.e., test successful) in the Cufflinks program. The false discovery rate (FDR) adjusted p-value (q value) was calculated automatically in Cuffdiff program. Genes were considered to be differentially expressed at an FDR q value < 0.05. Although the used Cufflinks pipeline included the reference-guided assembly stage, only known annotations were used in the following analysis. The RNA-Seq data were deposited in the NCBI Short Read Archive database with Accession number: PRJNA299102.

Functional annotation

The functional analysis of DEGs (see Additional file 1) was performed using the DAVID (The Database for Annotation, Visualization and Integrated Discovery) tool (http://david.abcc.ncifcrf.gov/) [58], [59]. The Rattus norvegicus genome was used as the background list for the over-representation analysis in DAVID. The GO option was utilized, and the significantly (p < 0.05) enriched biological processes and groups of genes possibly contributing to hypertensive phenotype in ISIAH rats were identified. The Kyoto Encyclopedia of Genes and Genomes Pathway Database (KEGG, http://www.genome.jp/kegg/) was used to identify pathways that were most significant to the data set. The genes related to hypertension and CNS diseases were detected according to the DEGs annotation in Rat Genome Database (RGD, http://rgd.mcw.edu/). The detection of transcription factors among DEGs was performed using gene annotations from GenBank (http://www.ncbi.nlm.nih.gov/gene/), an atlas of combinatorial transcriptional regulation in mouse and man [60] and the Panther classification system (http://www.pantherdb.org/).

Quantitative real-time PCR (qPCR)

The relative amount of target mRNA was measured by qPCR. Hypothalami were analyzed in 3-month old ISIAH and WAG rats (13 animals in each group). Total RNA was extracted using the TRI reagent (Molecular research center, USA). Remaining traces of genomic DNA were removed from the RNA samples using DNase I (Promega, USA) treatment, according to the manufacturer’s instructions.

cDNA was synthesized using 3 μg of total RNA, 0.25 nmol of random nonanucleotide primers (Biosan, Russia), 36 μl of reverse transcription buffer, 40 units of MoMLV (Vektor-Best, RF), and 0.4 mM dNTP. The cDNA was synthesized at 37 °C (1 h), 42 °C (30 min), and 50 °C (10 min). The enzyme was inactivated by heating the mixture at 75 °C for 5 min.

qPCR was performed in a final volume of 20 μl. The reaction volume contained a master mix with SYBR Green, forward and reverse primers (0.15 mM each), 1 unit of HotStart Taq polymerase (Vektor-Best, RF), and the cDNA template. The housekeeping gene Rpl30 encoding ribosomal protein L30 was used as a reference gene. Primer sequences are given in Additional file 8.

qPCR was carried out in an iCycler iQ4 Real-Time PCR Detection System (Bio-Rad Laboratories, USA) with an initial denaturation of 1 min at 95 °C followed by 35 cycles of 15 s at 95 °C, 20 s at each primer’s annealing temperatures (see Additional file 8), 20 s at 72 °C, fluorescence signal acquisition for 10 s at 83 °C and generation of the melting curve from 65 °C to 94 °C. The standard-curve quantitation method was applied [61]. Standard cDNA solution for plotting calibration curves was obtained by mixing aliquots from each of the synthesized cDNA samples. In each experiment, cDNA samples with primers for the target gene (three replicates per cDNA sample), the same samples with primers for the reference gene (three replicates), and the standard cDNA dilutions (1 : 4, 1 : 8, 1 : 16, and 1 : 32) with the primers for the target gene (two replicates) and with the primers for the reference gene (two replicates) were placed on the same plate. The relative amount of the tested cDNA was determined using calibration curves derived from the dilutions of the standard cDNA. Calibration curves were built using iCycler iQ4 Real-Time PCR Detection System software. The value for the target gene was further normalized against the qPCR level of the reference gene. The relative mRNA abundance was calculated as a ratio of the normalized (mRNA/Rpl30) mRNA level in the experimental ISIAH samples to the normalized (mRNA/Rpl30) mRNA level in the control WAG samples. The normalized mRNA level in the control samples of the WAG rats was assigned a value of 1. The data points that did not fit a normal distribution (assessed by 3-sigma limits) were removed from the analysis.

Statistical methods

Statistical calculations for qPCR data were performed with the software package Statistica v.6.0 (Statsoft, USA) using Student’s t-test. Differences were considered statistically significant at p < 0.05. The data were presented as the means and their standard errors (M ± S.E.M.).

The acquired RNA-seq data (FPKM values) were log transformed, centered and normalized. The principal coordinates method based on Euclidean metric distances was used for scaling the data sets. Then, the partial-least squares discriminant analysis (PLS-DA) approach was used to explore the pattern of co-variation for linear combinations between two blocks of variables [62], [63]. In the current experiment, the PLS-DA method and Pearson correlation were used to find a set of variables (expressed genes) that maximize the covariance between gene expression in two strains and fixed matrix with contrast values 1 and −1 for ISIAH rats and WAG rats, correspondingly. These procedures helped to construct the PLS-DA Axes maximizing the distances between ISIAH and WAG rats and to define the correlation between gene expression and PLS-DA Axis 1. The genes showing the most deviation along the first functionally meaningful synthetic PLS-DA Axis were considered to be genes that contribute the most to inter-strain differences.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files.

Additional files

Abbreviations

AA:

arachidonic acid

BP:

blood pressure

CNS:

central nervous system

DAVID:

Database for Annotation, Visualization and Integrated Discovery

DEG:

differentially expressed genes

EETs:

epoxyeicosatrienoic acids

FPKM:

fragments per kilobase of transcript per million mapped reads

GO:

Gene Ontology

HPA:

hypothalamic-pituitary-adrenal

ISIAH:

Inherited Stress-Induced Arterial Hypertension

KEGG:

Kyoto Encyclopedia of Genes and Genomes Pathway Database

MHC:

major histocompatibility complex

NK:

natural killer

PLS-DA:

partial-least squares discriminant analysis

qPCR:

quantitative real time polymerase chain reaction

RA:

retinoic acid

RGD:

Rat Genome Database

RNA-seq:

RNA sequencing

sEH:

soluble epoxide hydrolase

WAG:

Wistar Albino Glaxo

References

  1. 1.

    Dampney RA, Coleman MJ, Fontes MA, Hirooka Y, Horiuchi J, Li YW, et al: Central mechanisms underlying short- and long-term regulation of the cardiovascular system. Clin Exp Pharmacol Physiol. 2002, 29: 261-8. 10.1046/j.1440-1681.2002.03640.x.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    de Wardener HE: The hypothalamus and hypertension. Physiol Rev. 2001, 81: 1599-658.

    CAS  PubMed  Google Scholar 

  3. 3.

    Marques FZ, Morris BJ: Neurogenic hypertension: revelations from genome-wide gene expression profiling. Curr Hypertens Rep. 2012, 14: 485-91. 10.1007/s11906-012-0282-7.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Tain YL, Huang LT, Chan JY, Lee CT: Transcriptome analysis in rat kidneys: importance of genes involved in programmed hypertension. Int J Mol Sci. 2015, 16: 4744-58. 10.3390/ijms16034744.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Dornas WC, Silva ME: Animal models for the study of arterial hypertension. J Biosci. 2011, 36: 731-7. 10.1007/s12038-011-9097-y.

    Article  PubMed  Google Scholar 

  6. 6.

    Markel AL: Development of a new strain of rats with inherited stress-induced arterial hypertension. Genetic hypertension. Edited by: Sassard J. 1992, Colloque INSERM, Paris, 405-7.

    Google Scholar 

  7. 7.

    Markel AL, Maslova LN, Shishkina GT, Bulygina VV, Machanova NA, Jacobson GS: Developmental influences on blood pressure regulation in ISIAH rats. Development of the hypertensive phenotype: basic and clinical studies. Edited by: McCarty R, Blizard DA, Chevalier RL. 1999, Elsevier, Amsterdam- Lausanne- NewYork- Oxford- Shannon- Singapore- Tokyo, 493-526.

    Google Scholar 

  8. 8.

    Adarichev VA, Korokhov NP, Ostapchuk Ia V, Dymshits GM, Markel AL: Characterization of rat lines with normotensive and hypertensive status using genomic fingerprinting. Genetika. 1996, 32: 1669-72.

    CAS  PubMed  Google Scholar 

  9. 9.

    Markel AL, Redina OE, Gilinsky MA, Dymshits GM, Kalashnikova EV, Khvorostova YV, et al: Neuroendocrine profiling in inherited stress-induced arterial hypertension rat strain with stress-sensitive arterial hypertension. J Endocrinol. 2007, 195: 439-50. 10.1677/JOE-07-0254.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Marguerat S, Bahler J: RNA-seq: from technology to biology. Cell Mol Life Sci. 2010, 67: 569-79. 10.1007/s00018-009-0180-6.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-8. 10.1038/nmeth.1226.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Hu MC, Chiang EF, Tong SK, Lai W, Hsu NC, Wang LC, et al: Regulation of steroidogenesis in transgenic mice and zebrafish. Mol Cell Endocrinol. 2001, 171: 9-14. 10.1016/S0303-7207(00)00385-3.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Gomez-Sanchez EP, Ahmad N, Romero DG, Gomez-Sanchez CE: Is aldosterone synthesized within the rat brain?. Am J Physiol Endocrinol Metab. 2005, 288: E342-6. 10.1152/ajpendo.00355.2004.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Gomez-Sanchez EP, Gomez-Sanchez CM, Plonczynski M, Gomez-Sanchez C: Aldosterone synthesis in the brain contributes to Dahl salt-sensitive rat hypertension. Exp Physiol. 2010, 95: 120-30. 10.1113/expphysiol.2009.048900.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Oki K, Gomez-Sanchez EP, Gomez-Sanchez CE: Role of mineralocorticoid action in the brain in salt-sensitive hypertension. Clin Exp Pharmacol Physiol. 2012, 39: 90-5. 10.1111/j.1440-1681.2011.05538.x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Blomhoff R, Blomhoff HK: Overview of retinoid metabolism and function. J Neurobiol. 2006, 66: 606-30. 10.1002/neu.20242.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Maden M: Retinoic acid in the development, regeneration and maintenance of the nervous system. Nat Rev Neurosci. 2007, 8: 755-65. 10.1038/nrn2212.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Zhong M, Kawaguchi R, Ter-Stepanian M, Kassai M, Sun H: Vitamin A transport and the transmembrane pore in the cell-surface receptor for plasma retinol binding protein. PLoS One. 2013, 8: e73838-10.1371/journal.pone.0073838.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Zhou T, Huang C, Chen Y, Xu J, Shanbhag PD, Chen G: Methamphetamine regulation of sulfotransferase 1A1 and 2A1 expression in rat brain sections. Neurotoxicology. 2013, 34: 212-8. 10.1016/j.neuro.2012.09.010.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Borchardt RT, Schasteen CS: Phenol sulfotransferase. I. Purification of a rat liver enzyme by affinity chromatography. Biochim Biophys Acta. 1982, 708: 272-9. 10.1016/0167-4838(82)90437-X.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Benarroch EE, Balda MS, Finkielman S, Nahmod VE: Neurogenic hypertension after depletion of norepinephrine in anterior hypothalamus induced by 6-hydroxydopamine administration into the ventral pons: role of serotonin. Neuropharmacology. 1983, 22: 29-34. 10.1016/0028-3908(83)90257-5.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Yamamoto K, Sokabe T, Matsumoto T, Yoshimura K, Shibata M, Ohura N, et al: Impaired flow-dependent control of vascular tone and remodeling in P2X4-deficient mice. Nat Med. 2006, 12: 133-7. 10.1038/nm1338.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Spector AA, Fang X, Snyder GD, Weintraub NL: Epoxyeicosatrienoic acids (EETs): metabolism and biochemical function. Prog Lipid Res. 2004, 43: 55-90. 10.1016/S0163-7827(03)00049-3.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Imig JD, Zhao X, Capdevila JH, Morisseau C, Hammock BD: Soluble epoxide hydrolase inhibition lowers arterial blood pressure in angiotensin II hypertension. Hypertension. 2002, 39: 690-4. 10.1161/hy0202.103788.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Sellers KW, Sun C, Diez-Freire C, Waki H, Morisseau C, Falck JR, et al: Novel mechanism of brain soluble epoxide hydrolase-mediated blood pressure regulation in the spontaneously hypertensive rat. FASEB J. 2005, 19: 626-8.

    CAS  PubMed  Google Scholar 

  26. 26.

    Nagata A, Suzuki Y, Igarashi M, Eguchi N, Toh H, Urade Y, et al: Human brain prostaglandin D synthase has been evolutionarily differentiated from lipophilic-ligand carrier proteins. Proc Natl Acad Sci U S A. 1991, 88: 4020-4. 10.1073/pnas.88.9.4020.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Giles H, Leff P: The biology and pharmacology of PGD2. Prostaglandins. 1988, 35: 277-300. 10.1016/0090-6980(88)90093-7.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Fu ML: Do immune system changes have a role in hypertension?. J Hypertens. 1995, 13: 1259-65. 10.1097/00004872-199511000-00007.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Purcell ES, Gattone VH: Immune system of the spontaneously hypertensive rat. I. Sympathetic innervation. Exp Neurol. 1992, 117: 44-50. 10.1016/0014-4886(92)90109-4.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Harrison DG, Vinh A, Lob H, Madhur MS: Role of the adaptive immune system in hypertension. Curr Opin Pharmacol. 2010, 10: 203-7. 10.1016/j.coph.2010.01.006.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Parham P: MHC class I molecules and KIRs in human history, health and survival. Nat Rev Immunol. 2005, 5: 201-14. 10.1038/nri1570.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Nepom GT, Erlich H: MHC class-II molecules and autoimmunity. Annu Rev Immunol. 1991, 9: 493-525. 10.1146/annurev.iy.09.040191.002425.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Trott DW, Harrison DG: The immune system in hypertension. Adv Physiol Educ. 2014, 38: 20-4. 10.1152/advan.00063.2013.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Neumann H, Misgeld T, Matsumuro K, Wekerle H: Neurotrophins inhibit major histocompatibility class II inducibility of microglia: involvement of the p75 neurotrophin receptor. Proc Natl Acad Sci U S A. 1998, 95: 5779-84. 10.1073/pnas.95.10.5779.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Backman-Petersson E, Butcher GW, Hedlund G: Self-MHC class Ia (RT1-A(n)) protects cells co-expressing the activatory allogeneic MHC class Ib molecule (RT1-E(u)) from NK lysis. Int Immunol. 2000, 12: 843-50. 10.1093/intimm/12.6.843.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Lidman O, Olsson T, Piehl F: Expression of nonclassical MHC class I (RT1-U) in certain neuronal populations of the central nervous system. Eur J Neurosci. 1999, 11: 4468-72. 10.1046/j.1460-9568.1999.00904.x.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Poulsen KT, Armanini MP, Klein RD, Hynes MA, Phillips HS, Rosenthal A: TGFb2 and TGFb3 are potent survival factors for midbrain dopaminergic neurons. Neuron. 1994, 13: 1245-52. 10.1016/0896-6273(94)90062-0.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Masoud ST, Vecchio LM, Bergeron Y, Hossain MM, Nguyen LT, Bermejo MK, et al: Increased expression of the dopamine transporter leads to loss of dopamine neurons, oxidative stress and l-DOPA reversible motor deficits. Neurobiol Dis. 2015, 74: 66-75. 10.1016/j.nbd.2014.10.016.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Lamonerie T, Tremblay JJ, Lanctot C, Therrien M, Gauthier Y, Drouin J: Ptx1, a bicoid-related homeo box transcription factor involved in transcription of the pro-opiomelanocortin gene. Genes Dev. 1996, 10: 1284-95. 10.1101/gad.10.10.1284.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Lamolet B, Pulichino AM, Lamonerie T, Gauthier Y, Brue T, Enjalbert A, et al: A pituitary cell-restricted T box factor, Tpit, activates POMC transcription in cooperation with Pitx homeoproteins. Cell. 2001, 104: 849-59. 10.1016/S0092-8674(01)00282-3.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Bergen HT, Mizuno T, Taylor J, Mobbs CV: Resistance to diet-induced obesity is associated with increased proopiomelanocortin mRNA and decreased neuropeptide Y mRNA in the hypothalamus. Brain Res. 1999, 851: 198-203. 10.1016/S0006-8993(99)02186-1.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Konishi Y, Koosaka Y, Maruyama R, Imanishi K, Kasahara K, Matsuda A, et al: L-Ornithine intake affects sympathetic nerve outflows and reduces body weight and food intake in rats. Brain Res Bull. 2015, 111: 48-52. 10.1016/j.brainresbull.2014.11.004.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Knauf C, Rieusset J, Foretz M, Cani PD, Uldry M, Hosokawa M, et al: Peroxisome proliferator-activated receptor-alpha-null mice have increased white adipose tissue glucose utilization, GLUT4, and fat mass: Role in liver and brain. Endocrinology. 2006, 147: 4067-78. 10.1210/en.2005-1536.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Nishiyama K, Konishi A, Nishio C, Araki-Yoshida K, Hatanaka H, Kojima M, et al: Expression of cystatin C prevents oxidative stress-induced death in PC12 cells. Brain Res Bull. 2005, 67: 94-9. 10.1016/j.brainresbull.2005.05.020.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Marques FZ, Campain AE, Davern PJ, Yang YH, Head GA, Morris BJ: Global identification of the genes and pathways differentially expressed in hypothalamus in early and established neurogenic hypertension. Physiol Genomics. 2011, 43: 766-71. 10.1152/physiolgenomics.00009.2011.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Fang X: Soluble epoxide hydrolase: a novel target for the treatment of hypertension. Recent Pat Cardiovasc Drug Discov. 2006, 1: 67-72. 10.2174/157489006775244227.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Maass AH, De Jong AM, Smit MD, Gouweleeuw L, de Boer RA, Van Gilst WH, et al: Cardiac gene expression profiling - the quest for an atrium-specific biomarker. Neth Heart J. 2010, 18: 610-4. 10.1007/s12471-010-0844-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Golz S, Summer H, Geerts A, Bruggemeier U, Albrecht-Kupper B, Klein M, Steppan S, Ellinghaus P, D′Urso D, Seewald M, Milting H. LTBP2 as a biomarker and diagnostic target. 2012. Patent US 8088584 B2.

  49. 49.

    Recklies AD, Ling H, White C, Bernier SM: Inflammatory cytokines induce production of CHI3L1 by articular chondrocytes. J Biol Chem. 2005, 280: 41213-21. 10.1074/jbc.M510146200.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Francescone RA, Scully S, Faibish M, Taylor SL, Oh D, Moral L, et al: Role of YKL-40 in the angiogenesis, radioresistance, and progression of glioblastoma. J Biol Chem. 2011, 286: 15332-43. 10.1074/jbc.M110.212514.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Ma JY, Li RH, Huang K, Tan G, Li C, Zhi FC: Increased expression and possible role of chitinase 3-like-1 in a colitis-associated carcinoma model. World J Gastroenterol. 2014, 20: 15736-44. 10.3748/wjg.v20.i42.15736.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, et al: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320: 1344-9. 10.1126/science.1158441.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Cloonan N, Forrest AR, Kolle G, Gardiner BB, Faulkner GJ, Brown MK, et al: Stem cell transcriptome profiling via massive-scale mRNA sequencing. Nat Methods. 2008, 5: 613-9. 10.1038/nmeth.1223.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Paxinos G, Watson C: The rat brain in stereotaxic coordinates. 2006, Academic Press, Elsevier, San Diego, California

    Google Scholar 

  56. 56.

    Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14: R36-10.1186/gb-2013-14-4-r36.

    Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L: Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013, 31: 46-53. 10.1038/nbt.2450.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57. 10.1038/nprot.2008.211.

    CAS  Article  Google Scholar 

  59. 59.

    Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.

    Article  Google Scholar 

  60. 60.

    Ravasi T, Suzuki H, Cannistraci CV, Katayama S, Bajic VB, Tan K, et al: An atlas of combinatorial transcriptional regulation in mouse and man. Cell. 2010, 140: 744-52. 10.1016/j.cell.2010.01.044.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Ginzinger DG: Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp Hematol. 2002, 30: 503-12. 10.1016/S0301-472X(02)00806-8.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Barker M, Rayens W: Partial least squares for discrimination. J Chemom. 2003, 17: 166-73. 10.1002/cem.785.

    CAS  Article  Google Scholar 

  63. 63.

    Szymanska E, Saccenti E, Smilde AK, Westerhuis JA: Double-check: validation of diagnostic statistics for PLS-DA models in metabolomics studies. Metabolomics. 2012, 8 (Suppl 1): 3-16. 10.1007/s11306-011-0330-3.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Moise AR, Kuksa V, Imanishi Y, Palczewski K: Identification of all-trans-retinol:all-trans-13,14-dihydroretinol saturase. J Biol Chem. 2004, 279: 50230-42. 10.1074/jbc.M409130200.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Takakuwa Y, Nishino H, Ishibe Y, Ishibashi T: Properties and kinetics of membrane-bound enzymes when both the enzyme and substrate are components of the same microsomal membrane. Studies on lathosterol 5-desaturase. J Biol Chem. 1994, 269: 27889-93.

    CAS  PubMed  Google Scholar 

  66. 66.

    Saether T, Tran TN, Rootwelt H, Christophersen BO, Haugen TB: Expression and regulation of delta5-desaturase, delta6-desaturase, stearoyl-coenzyme A (CoA) desaturase 1, and stearoyl-CoA desaturase 2 in rat testis. Biol Reprod. 2003, 69: 117-24. 10.1095/biolreprod.102.014035.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Hori S, Ohtsuki S, Tachikawa M, Kimura N, Kondo T, Watanabe M, et al: Functional expression of rat ABCG2 on the luminal side of brain capillaries and its enhancement by astrocyte-derived soluble factor(s). J Neurochem. 2004, 90: 526-36. 10.1111/j.1471-4159.2004.02537.x.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Alabert C, Bukowski-Wills JC, Lee SB, Kustatscher G, Nakamura K, de Lima AF, et al: Nascent chromatin capture proteomics determines chromatin dynamics during DNA replication and identifies unknown fork components. Nat Cell Biol. 2014, 16: 281-93. 10.1038/ncb2918.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Barnes LD, Garrison PN, Siprashvili Z, Guranowski A, Robinson AK, Ingram SW, et al: Fhit, a putative tumor suppressor in humans, is a dinucleoside 5′,5‴-P1, P3-triphosphate hydrolase. Biochemistry. 1996, 35: 11529-35. 10.1021/bi961415t.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Yu H, Tardivo L, Tam S, Weiner E, Gebreab F, Fan C, et al: Next-generation sequencing to generate interactome datasets. Nat Methods. 2011, 8: 478-80. 10.1038/nmeth.1597.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Guardiola-Serrano F, Haendeler J, Lukosz M, Sturm K, Melchner H, Altschmied J: Gene trapping identifies a putative tumor suppressor and a new inducer of cell migration. Biochem Biophys Res Commun. 2008, 376: 748-52. 10.1016/j.bbrc.2008.09.070.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Kremmidiotis G, Gardner AE, Settasatian C, Savoia A, Sutherland GR, Callen DF: Molecular and functional analyses of the human and mouse genes encoding AFG3L1, a mitochondrial metalloprotease homologous to the human spastic paraplegia protein. Genomics. 2001, 76: 58-65. 10.1006/geno.2001.6560.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Casari G, De Fusco M, Ciarmatori S, Zeviani M, Mora M, Fernandez P, et al: Spastic paraplegia and OXPHOS impairment caused by mutations in paraplegin, a nuclear-encoded mitochondrial metalloprotease. Cell. 1998, 93: 973-83. 10.1016/S0092-8674(00)81203-9.

    CAS  Article  PubMed  Google Scholar 

  74. 74.

    Schaeren-Wiemers N, Valenzuela DM, Frank M, Schwab ME: Characterization of a rat gene, rMAL, encoding a protein with four hydrophobic domains in central and peripheral myelin. J Neurosci. 1995, 15: 5753-64.

    CAS  PubMed  Google Scholar 

  75. 75.

    Erne B, Sansano S, Frank M, Schaeren-Wiemers N: Rafts in adult peripheral nerve myelin contain major structural myelin proteins and myelin and lymphocyte protein (MAL) and CD59 as specific markers. J Neurochem. 2002, 82: 550-62. 10.1046/j.1471-4159.2002.00987.x.

    CAS  Article  PubMed  Google Scholar 

  76. 76.

    Wang J, Gebre AK, Anderson RA, Parks JS: Cloning and in vitro expression of rat lecithin:cholesterol acyltransferase. Biochim Biophys Acta. 1997, 1346: 207-11. 10.1016/S0005-2760(97)00055-6.

    CAS  Article  PubMed  Google Scholar 

  77. 77.

    Bojesen E: Reversible inhibition of LCAT by penicillins. Kinetic and mechanistic investigations. Scand J Clin Lab Invest Suppl. 1978, 150: 26-31. 10.3109/00365517809104896.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Steiglitz BM, Keene DR, Greenspan DS: PCOLCE2 encodes a functional procollagen C-proteinase enhancer (PCPE2) that is a collagen-binding protein differing in distribution of expression and post-translational modification from the previously described PCPE1. J Biol Chem. 2002, 277: 49820-30. 10.1074/jbc.M209891200.

    CAS  Article  PubMed  Google Scholar 

  79. 79.

    Zeldin DC: Epoxygenase pathways of arachidonic acid metabolism. J Biol Chem. 2001, 276: 36059-62. 10.1074/jbc.R100030200.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Poisik O, Raju DV, Verreault M, Rodriguez A, Abeniyi OA, Conn PJ, et al: Metabotropic glutamate receptor 2 modulates excitatory synaptic transmission in the rat globus pallidus. Neuropharmacology. 2005, 49 (Suppl 1): 57-69. 10.1016/j.neuropharm.2005.03.006.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Lee SH, Vidal SM: Functional diversity of Mx proteins: variations on a theme of host resistance to infection. Genome Res. 2002, 12: 527-30. 10.1101/gr.20102.

    CAS  Article  PubMed  Google Scholar 

  82. 82.

    Brotto MA, Biesiadecki BJ, Brotto LS, Nosek TM, Jin JP: Coupled expression of troponin T and troponin I isoforms in single skeletal muscle fibers correlates with contractility. Am J Physiol Cell Physiol. 2006, 290: C567-76. 10.1152/ajpcell.00422.2005.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Yamamoto Y, Yoshizawa T, Kamio S, Aoki O, Kawamata Y, Masushige S, et al: Interactions of transthyretin (TTR) and retinol-binding protein (RBP) in the uptake of retinol by primary rat hepatocytes. Exp Cell Res. 1997, 234: 373-8. 10.1006/excr.1997.3642.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Hyytiainen M, Keski-Oja J: Latent TGF-beta binding protein LTBP-2 decreases fibroblast adhesion to fibronectin. J Cell Biol. 2003, 163: 1363-74. 10.1083/jcb.200309105.

    Article  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Tsuruga E, Oka K, Hatakeyama Y, Isokawa K, Sawa Y: Latent transforming growth factor-beta binding protein 2 negatively regulates coalescence of oxytalan fibers induced by stretching stress. Connect Tissue Res. 2012, 53: 521-7. 10.3109/03008207.2012.702816.

    CAS  Article  PubMed  Google Scholar 

  86. 86.

    Akizu N, Cantagrel V, Zaki MS, Al-Gazali L, Wang X, Rosti RO, et al: Biallelic mutations in SNX14 cause a syndromic form of cerebellar atrophy and lysosome-autophagosome dysfunction. Nat Genet. 2015, 47: 528-34. 10.1038/ng.3256.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are grateful to JSC Genoanalytica (Moscow, Russia) for conducting the technological part of the experiment and the primary statistical analysis. The RNA-Seq experiments, analysis of transcriptomes, and validation of results by qPCR were supported by the Russian Science Foundation, grant No. 14-15-00118. The development of conception and use of multivariate statistical analysis was supported by RFBR, grant No. 13-07-00315.

Declarations

Publication costs for this article were funded by the Russian Science Foundation, grant No. 14-15-00118.

This article has been published as part of BMC Genetics Volume 17 Supplement 1, 2016: Selected articles from the 7th International Young Scientists School “Systems Biology and Bioinformatics” (SBB’2015): Genetics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenetics/supplements/17/S1.

Author information

Affiliations

Authors

Corresponding author

Correspondence to Olga E. Redina.

Additional information

Competing interests

The authors declare that they have no competing interests. The authors alone are responsible for the content and writing of the paper.

Authors’ contributions

LK performed quantitative real time PCR, participated in interpretation of data, and helped to draft the manuscript; NE has made a substantial contribution to bioinformatics analysis; VE has made substantial contributions to conception and the use of multivariate statistical analysis; AM has made substantial contributions to the conception and design of the study and has participated in the interpretation of data; OR performed the multivariate statistical analysis, participated in the interpretation of data and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Genes differentially expressed in hypothalamus of hypertensive ISIAH and normotensive control WAG rats. (XLS 40 kb) (XLS 40 KB)

Additional file 2: Differentially expressed genes related to Gene Ontlogy (GO) terms for biological processes. (XLS 58 kb) (XLS 58 KB)

Additional file 3: Metabolic pathways and genes differentially expressed in ISIAH and WAG hypothalami. (DOC 86 kb) (DOC 86 KB)

Additional file 4: PLS-DA Axis maximizing the distance between ISIAH and WAG rats. (JPG 57 kb) (JPG 57 KB)

Additional file 5: The distribution of expressed genes along the PLS-DA Axis 1. (JPG 117 kb) (JPG 118 KB)

Additional file 6: Comparison of gene expression level measurements obtain by RNA-seq and qPCR. (JPG 62 kb) (JPG 63 KB)

Additional file 7: The summary statistics for the sequenced libraries. (XLS 20 kb) (XLS 20 KB)

Additional file 8: Primers used in real-time PCR. (DOC 35 kb) (DOC 35 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Klimov, L.O., Ershov, N.I., Efimov, V.M. et al. Genome-wide transcriptome analysis of hypothalamus in rats with inherited stress-induced arterial hypertension. BMC Genet 17, S13 (2016). https://doi.org/10.1186/s12863-015-0307-8

Download citation

Keywords

  • Stress-sensitive hypertension
  • Hypothalamus
  • Transcriptional profiling
  • RNA-Seq
  • ISIAH rats