- Research article
- Open Access
Vascular microarray profiling in two models of hypertension identifies caveolin-1, Rgs2 and Rgs5 as antihypertensive targets
BMC Genomics volume 8, Article number: 404 (2007)
Hypertension is a complex disease with many contributory genetic and environmental factors. We aimed to identify common targets for therapy by gene expression profiling of a resistance artery taken from animals representing two different models of hypertension. We studied gene expression and morphology of a saphenous artery branch in normotensive WKY rats, spontaneously hypertensive rats (SHR) and adrenocorticotropic hormone (ACTH)-induced hypertensive rats.
Differential remodeling of arteries occurred in SHR and ACTH-treated rats, involving changes in both smooth muscle and endothelium. Increased expression of smooth muscle cell growth promoters and decreased expression of growth suppressors confirmed smooth muscle cell proliferation in SHR but not in ACTH. Differential gene expression between arteries from the two hypertensive models extended to the renin-angiotensin system, MAP kinase pathways, mitochondrial activity, lipid metabolism, extracellular matrix and calcium handling. In contrast, arteries from both hypertensive models exhibited significant increases in caveolin-1 expression and decreases in the regulators of G-protein signalling, Rgs2 and Rgs5. Increased protein expression of caveolin-1 and increased incidence of caveolae was found in both smooth muscle and endothelial cells of arteries from both hypertensive models.
We conclude that the majority of differences in gene expression found in the saphenous artery taken from rats with two different forms of hypertension reflect distinctive morphological and physiological alterations. However, changes in common to caveolin-1 expression and G protein signalling, through attenuation of Rgs2 and Rgs5, may contribute to hypertension through augmentation of vasoconstrictor pathways and provide potential targets for common drug development.
Hypertension is a complex condition with many contributory genetic and environmental factors , which have constrained identification of the genes, biochemical pathways and pathophysiological processes that accompany its development. Hypertension is a critical risk factor for stroke, renal failure and cardiovascular disease and has been associated with a range of well recognized complications, including congestive heart failure, left ventricular hypertrophy, and atherosclerosis, and with more subtle structural and functional changes to target organs and tissues occurring as a consequence of enhanced salt sensitivity , enhanced angiotensin II sensitivity  and structural modifications to the arterial wall [4, 5]. Genetic factors are thought to partly account for structural changes by directly influencing the quantity, density and organization of the extracellular matrix and also indirectly through association with risk factors such as age, gender, blood pressure, total and high-density lipoprotein cholesterol, smoking and diabetes .
A variety of both genetic, such as the spontaneously hypertensive rat (SHR), and nongenetic, such as glucocorticoid-induced hypertension, animal models of hypertension have been developed . Administration of adrenocorticotropic hormone (ACTH) to rats provides a quick and reproducible means of inducing hypertension by stimulating the adrenal production of corticosterone . ACTH-induced hypertension was originally thought to involve altered kidney function and loss of control of blood volume  but is now known to involve disruption of the nitric oxide synthase (NOS) pathway at various points, affecting NO bioavailability, oxidative stress and function of the endothelium [10, 11]. The genetically hypertensive SHR is a renin-angiotensin system-dependent model of hypertension that also exhibits oxidative stress and disrupted endothelial NOS (eNOS) activity [12, 13]. Such convergence in effectors holds the promise of some overlap in causative factors.
Hypertension is generally characterized by remodeling of arteries, such as thickening of the media and intima, reduced lumen diameter and reduced distensibility . Endothelial dysfunction is also a characteristic of hypertension and alterations to endothelial morphology and density have been described [15–17]. Resistance vessels have a primary role in autoregulation of the blood supply to the vascular bed and changes to their structure and function are crucial to the onset of target organ damage . Angiotensin II has a prominent role in affecting resistance vessels as a pro-inflammatory stimulator of signaling pathways promoting vascular growth, remodeling and oxidative stress [18, 19]. Both hypertrophic and eutrophic remodeling have been documented in arteries of SHR, perhaps related to the size of the vessel and intraluminal pressures experienced [20–23]. Whether similar changes to resistance vessels accompany ACTH-induced hypertension is not known. Microarray expression profiling offers a means of studying the modifications to global gene expression that accompany these events  and hence may provide novel insights into potential targets for common therapeutic intervention. The aim of the present study was therefore to compare gene expression and morphology of a resistance artery from adult normotensive Wistar-Kyoto (WKY) rats, SHR and ACTH-induced hypertensive rats; the latter two as representatives of renin-angiotensin-dependent and steroid-dependent hypertension, respectively. Although the majority of genetic changes differed between the two hypertensive models, changes in common were found to caveolin-1 expression and the regulators of G protein signalling, Rgs2 and Rgs5.
Blood pressure of SHR and ACTH-treated rats was significantly elevated over that of age-matched WKY rats (Table 1).
Anatomy of saphenous arteries from normotensive and hypertensive rats
The internal vessel circumference of arteries taken from both ACTH-treated rats and SHR showed a trend towards a decrease, compared with WKY, although only saphenous arteries from SHR exhibited hypertrophy of the smooth muscle cell layers (Figure 1, Table 1). In the endothelium, there was a significant reduction in the area of endothelial cells in SHR, while area was significantly increased during ACTH-treatment (Table 1). There was no significant difference between experimental groups in the density of myoendothelial gap junctions per endothelial cell, which were rarely encountered in either WKY, SHR or ACTH-treated rats (Table 1).
Assessment of RNA degradation
The extent of RNA degradation was assessed by comparing the MAS5-derived 3'/5' and 3'/M' (middle) ratios of the housekeeping genes, GAPDH and hexokinase sample (Table 2). The 3'/5' ratios for GAPDH for all 12 samples were in the range 3–7, which is acceptable for doubly-amplified RNA, while the 3'/M' ratios for GAPDH were less than 2. Similar data were found for the 3'/5' and 3'/M ratios for hexokinase although some of the ratios tended to be greater than those for GAPDH, perhaps due to the longer hexokinase transcripts (3653 bp, versus 1268 bp for GAPDH).
Statistical analysis of microarrays
Of the 15,878 probesets on each array, 9098 were present in all replicates of at least one experimental group (Affymetrix MAS5 Statistical Algorithm) and considered for further analysis. One-way ANOVA on each of the 9098 probesets showed significant differences in expression of 870 genes amongst the 3 experimental groups (F-test p value cutoff of 0.02) or 136 genes (p value cutoff of 0.001). The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO) at the following website  and are accessible through GEO Series accession number GSE8051.
The 50 most significant probesets found by ANOVA were superimposed on a GE-biplot of the 9098 probesets present in all arrays of at least one experimental group (Figure 2A). The locations of the probesets are shown as red asterisks with the 50 most significant genes indicated. The locations of the arrays are connected to the centre of the plot and labelled by their respective names. Figure 2A shows that the replicates within each experimental group are clustered together relative to the other groups, indicating consistency in measurement. In this genome-wide view the three experimental groups, WKY, ACTH and SHR groups are well separated, providing evidence of overall transcriptional differences between the groups.
The genes separating from the dense group in the centre of Figure 2A have relatively high variance. The labelled genes have a statistically significant difference in their means in the three groups. Generally the statistically significant genes on the left of the plot have relatively low expressions in SHR but high and approximately equal expression in WKY and ACTH. Those statistically significant genes on the right of the centre have relatively high expression values in SHR and low and approximately equal expression in WKY and ACTH. These patterns in gene expression are the dominant genome-wide patterns. However there are a number of genes which are relatively highly expressed in ACTH (eg. 9, 31, 48) compared to WKY and SHR.
The heat map of the 50 most significant genes in Figure 2B shows that the arrays cluster according to experimental groups but the experimental groups differ substantially, in agreement with clustering in the GE-biplot.
Validation of microarray results by quantitative PCR
Eleven out of 12 genes (6 in SHR, 6 in ACTH) which showed significant variation from WKY in the microarrays were validated by quantitative PCR (Table 3). The one variant gene (adenylate cyclase3) was significantly different in the SHR by microarray, but found not to differ from WKY by quantitative PCR. In contrast, four of 14 genes (7 in SHR, 7 in ACTH) which did not show any significant difference in expression in the microarrays relative to WKY, were found to be significantly different to WKY when assayed using quantitative PCR (Table 3). Three of these genes were found to be decreased in expression relative to WKY and one gene increased in expression (Table 3). These data concord with a recent analysis of microarray data which showed a high false-negative rate, or low sensitivity, especially in the case of genes showing small fold-changes .
Differentially and co-ordinately expressed genes between the two models of hypertension
Genes that were identified as significantly altered in expression by ANOVA in eight multiple comparisons of SHR and ACTH relative to WKY are given in Additional Files 1, 2, 3, 4, 5, 6, 7, 8: Tables S1–8. Genes identified in six dual comparisons of either SHR or ACTH with WKY are given in Additional Files 9, 10, 11, 12, 13, 14: Tables S9–14. Gene names, gene descriptions, biological function, Affymetrix IDs, accession numbers, expression levels and p-values are provided. The biological association networks of selected genes for each model of hypertension are shown in Figure 3A,B.
Table 4 summarises the changes in genes which relate to a range of processes or pathways which could impact on blood pressure regulation. The majority of changes reflect the specific morphological and physiological alterations occurring in the two different types of hypertension. Saphenous arteries from SHR showed upregulation of genes associated with renin-angiotensin signalling through MAP kinase ERK1/2 pathway and were characterised by upregulation of growth promoters and down regulation of growth suppressors (Table 4). In contrast, arteries from ACTH-treated animals showed upregulation of the p38MAP kinase pathway, diverse effects on growth promoters and suppressors (Table 4) and prominent involvement of the transcription factor Jun (Figure 3B). Differences were also found in expression of genes involved in mitochondrial activity, lipid metabolism, calcium handling, adenylate cyclase, ion transport and protein kinase C. In contrast, coordinate upregulation in both forms of hypertension was found for caveolin-1 and its interaction and modification by the membrane-associated tyrosine kinase Fyn (Figure 3, blue shape), and coordinate downregulation of Rgs2 and Rgs5, the regulators of G protein signalling (Tables 3, 4).
Distribution of caveolin-1 and caveolae in saphenous arteries of normotensive and hypertensive rats
Expression of caveolin-1 was found in the membranes of both the endothelial and smooth muscle cells of saphenous arteries of WKY rats (Figure 4A, D, arrows). Expression varied both within and between cells in the endothelium (Figure 4A) and also throughout the muscle layers (Figure 4D). An increase in caveolin-1 expression was found in both endothelial cells and smooth muscle cells of arteries taken from both SHR and ACTH, compared to WKY rats (Figure 4B, C, compare 4A; Figure 4E, F compare 4D). Paired confocal images of either whole mounts or sections were taken with the same voltage and pinhole settings so that comparisons could be made of the extent of staining. Preincubation of the antibody with the immunogen abolished staining in endothelial and smooth muscle cells (Figure 4G).
At the ultrastructural level, caveolae were found along the membranes of both endothelial and smooth muscle cells of arteries from WKY rats (Figure 5A, B), although there were fewer on the abluminal than luminal endothelial surface. While the incidence of caveolae was increased in both cell types in the SHR (Figure 5C, D), the differential distribution in the endothelium was still apparent. In arteries from ACTH rats, a differential distribution in the endothelium was less obvious as caveolae were found in abundance on both sides of the endothelial cells (Figure 5E). Caveolae were also more abundant throughout the muscle layers of ACTH-treated rats, than of WKY rats (Figure 5F, compare 5B).
Comparative gene expression profiling of a resistance artery taken from rats with either renin-angiotensin-dependent or glucocorticoid-induced hypertension has identified over 800 genes showing significant changes. Statistical analysis of the microarrays showed consistency in gene expression amongst the replicates within each experimental group but evidence of substantial transcriptional differences between the two hypertensive models and the normotensive controls. The majority of these differences reflected differential vascular remodeling, as confirmed by our anatomical studies, as well as activation of specific intracellular pathways. Nevertheless, arteries from both models exhibited significantly increased expression of caveolin-1, which binds a number of important signalling molecules including endothelial nitric oxide synthase , and significantly decreased expression of Rgs2 and Rgs5, which terminate the signalling of the G proteins Gq and Gi/o . Together these common effects would be expected to elevate blood pressure through augmentation of vasoconstrictor mechanisms.
Saphenous arteries from SHRs were characterised by hypertrophic remodeling, in contrast to the arteries from the ACTH rats in which no alterations were found in the number of smooth muscle cell layers. Similar changes have been described in other arteries of SHR [20, 22, 23], typical of remodeling of similarly sized resistance vessels during angiotensin II-induced hypertension . The absence of hypertrophy in the saphenous artery of ACTH-treated rats suggests that inward eutrophic remodeling  may be more common to this form of hypertension, perhaps related to the lower blood pressures experienced in comparison with those of the SHR. In the endothelium of SHR arteries, there was a significant reduction in the area of endothelial cells, consistent with changes observed in other muscular arteries of SHR [29–31]. Surprisingly, the area of endothelial cells was increased during ACTH-treatment resulting in a decrease in overall endothelial cell numbers. Changes in gene expression were consistent with cell proliferation in angiotensin II-dependent hypertension but not in glucocorticoid-induced hypertension. Thus, expression of growth promoters was upregulated and growth inhibitors suppressed in SHR arteries, while effects on both growth promoters and inhibitors were more diverse in the ACTH arteries. A lack of change in myoendothelial gap junctions suggests that these structures are not involved in the remodeling processes in saphenous arteries in either form of hypertension.
A large number of other genetic differences in the arteries from the two hypertensive models emerged from the study. These related to lipid metabolism, mitochondrial activity, calcium handling, ion transport, the adenylate cyclase pathway and protein kinase C. Not surprisingly, saphenous arteries from SHR showed upregulation of genes consistent with renin-angiotensin hypertension, and angiotensin II signaling pathways through Raf and MEK1/2 . Upregulation of growth promoters and down regulation of growth suppressors as found in the SHR samples would be expected to stimulate cell growth and proliferation, mobilise intracellular Ca2+ and reactive oxygen species [18, 32]. In contrast, the saphenous arteries from ACTH-treated animals showed upregulation of the p38MAP kinase pathway which is involved in inflammation and apoptosis, as well as growth and differentiation , accentuating mechanistic differences between the two models of hypertension. While alterations to the extracellular matrix were found in both hypertensive models, substantial changes involving multiple collagen genes were found predominantly in the ACTH model. Such changes have been linked with hypertension  and could impact on cytoskeletal signalling through integrin receptors.
Downstream mechanisms underlying smooth muscle contraction and calcium sensitization also varied between the two hypertensive models. Calcium/calmodulin-dependent protein kinase II (CaMkII) has an important role in regulating smooth muscle contractility by inhibiting myosin light chain kinase (MLCK) and the phosphorylation of myosin light chain . In SHR, upregulation of CaMkII may shift smooth muscle contraction toward Rho-dependent mechanisms and decreased phosphatase activity. In contrast, in ACTH, CaMkII was downregulated suggesting that smooth muscle contraction could be more influenced by MLCK, in addition to Rho-dependent mechanisms.
In the face of these substantial differences in gene expression between the two hypertensive models, we found coordinate upregulation in gene expression of caveolin-1 and increased caveolin-1 protein in both endothelial and smooth muscle cells. Our ultrastructural observations of increased incidence of caveolae in these sites provides strong confirmation of the immunohistochemical findings since caveolae are absent in caveolin-1 knockout mice . While some of the genetic changes in SHR relative to WKY might have arisen through decades of isolation from the parent WKY strain following inbreeding, coordinate changes in gene expression between the SHR and the ACTH-treated WKY rats cannot be attributed to such a cause. Caveolae are now recognized as important sites for colocalisation of signaling molecules, including receptor tyrosine kinases, G-protein-coupled receptors and G-proteins, and hence are critical for vascular function . The inactivation of eNOS through binding to caveolin-1 is the best characterized function for caveolae in endothelial cells , highlighting their crucial role in vascular homeostasis. Alterations to caveolae in smooth muscle cells will also impact on cell signalling and recent studies of idiopathic pulmonary arterial hypertension in humans have concluded that increased smooth muscle expression of caveolin-1 contributes to the pathophysiology of the disease through increased capacitative calcium entry . Whether upregulation of caveolin-1 in smooth muscle cells and in endothelial cells could serve to induce hypertension warrants further investigation.
In contrast to caveolin-1, the G protein regulators, Rgs2 and Rgs5, were coordinately downregulated in both SHR and ACTH saphenous arteries. These molecules terminate G protein signalling through Gq and Gi/o by activating the GTPase activity of the Gα subunit. Downregulation thus potentiates the effect of vasoconstrictors like endothelin and angiotensin II . A role in blood pressure regulation has been demonstrated for Rgs2, since knockout mice are severely hypertensive . Activation of Rgs2 occurs in resistance and conduit arteries through binding of protein kinase G, the downstream mediator of nitric oxide action [40, 41]. Interestingly, Rgs2 was also decreased in aortic smooth muscle during L-NAME-induced hypertension . Upregulation of caveolin-1 and decreased eNOS function may therefore play a role in the downregulation of Rgs2.
Less is known about a role in blood pressure regulation for Rgs5, which has been described as a marker of arterial but not venous smooth muscle cells  and is downregulated during proliferation to form atherosclerotic plaques . In contrast, upregulation of Rgs5 has been reported during neovascularisation of tumours . Rgs5 has also been identified in endothelial cells of cerebral capillaries where it was downregulated in stroke-prone SHR compared to SHR . Our results demonstrating downregulation in arteries from both SHR and ACTH rats would suggest a role in hypertension rather than proliferation since smooth muscle hypertrophy was not a feature of the ACTH model.
Study of changes in gene expression in two different models of hypertension have highlighted numerous differences which can be linked to variation in the form of vascular remodeling and the unique signalling pathways affecting contractility and cell growth. In contrast, we have identified common changes to upstream targets crucial for the alteration in the balance between vasoconstriction and vasodilation as occurs in hypertension. Common increases in expression of caveolin-1 are likely to increase vasoconstriction through direct effects on calcium influx in smooth muscle cells and indirect inhibitory effects on eNOS in endothelial cells. Decreases in Rgs2 and Rgs5 will augment these effects by prolonging the actions of Gq and Gi/o and their downstream vasoconstrictor pathways. Further studies examining the functional implications of these data and aimed at reversing these gene changes during the different forms of hypertension are needed to determine whether these molecules could be considered as therapeutic drug targets.
Sixteen-week-old male WKY and age-matched SHR were used. These rats were originally from stock derived from the National Institutes of Health, Bethesda, MD. A group of 12-week-old male WKY rats was administered ACTH (0.2 mg/kg/day) subcutaneously for 4 weeks. Systolic blood pressure was measured using tail-cuff plethysmography. All experiments were approved by the Australian National University Animal Experimentation Ethics Committee under guidelines of the National Health and Medical Research Council of Australia.
All analyses were made on the first lateral branch of the saphenous artery. This vessel was defined as a resistance vessel due to its small diameter (<250 μm), since feed arteries and the microcirculation contribute to blood pressure regulation in rats .
Ultrastructural studies were conducted on a minimum of 3 arteries each from a different rat, as previously described [48, 49]. WKY, ACTH and SHR rats were anaesthetized (i.p., 44 and 8 mg/kg, ketamine and xylazine, respectively) and perfused via the left ventricle with 0.1% BSA, 0.1% NaNO3 and 10 U/ml heparin and fixed with 1% paraformaldehyde, 3% glutaraldehyde in 0.1 M sodium cacodylate buffer, with 10 mM betaine, 0.2 mM CaCl and 0.15 M sucrose, pH 7.4, using standard procedures . The first lateral branch of the saphenous artery diverging towards the knee was removed and processed for electron microscopy [48, 49]. Serial transverse sections (~100 nm thick) totalling ~5 μ m of vessel length were cut from each vessel using 3–4 animals/experimental group. Myoendothelial gap junctions connecting the endothelium to the inner layer of smooth muscle were counted and photographed at ×3,500 to ×60,000 on a Phillips 7100 transmission electron microscope. Vessel circumference at the internal elastic lamina, the number of smooth muscle cell layers and endothelial cell areas were determined for each vessel as previously described [48, 49]. Statistical significance was tested using one-way ANOVA followed by Student's t-test with Bonferroni correction for multiple groups. A p-value < 0.05 was taken as significant. The incidence of caveolae was assessed in montages of the vessel wall using a minimum of 3 arteries, each from a different rat, for each rat group.
Rats were anaesthetized as above and perfused with 2% paraformaldehyde in 0.1 M sodium phosphate buffer. The lateral branch of the saphenous artery was removed from both hind legs. Arteries were processed either as whole mounts or were sectioned at 10 μm. Preparations were incubated overnight at room temperature with rabbit antibodies against human caveolin-1 (Santa Cruz; 1:200 for whole mounts, 1:1000 for sections) followed by 1 hour in Cy3 conjugated anti-rabbit IgG (1:100, Jackson Immunoresearch Laboratories). For whole mounts, optical confocal series were taken throughout the thickness of the endothelium and recombined to a single image. For tissue sections, optical confocal series were taken throughout the section thickness and recombined to a single image. All images were collected using the same pinhole and voltage settings so that comparisons could be made of the intensity and extent of staining in arteries from the different rat groups (3 animals/group). Specificity of the staining was tested by preincubation of the antibody with a ten-fold excess of the immunogen.
Extraction of total RNA
Rats were anesthetized with ether and killed by decapitation. The saphenous artery branch was removed, stripped of extraneous material and stored in RNAlater (Ambion). Three different samples were prepared for each experimental group, i.e. WKY, SHR and ACTH. The choice of three biological replicates was based on the constraints of limited availability of mRNA and the high costs of microarray experiments. This is a common choice in the microarray literature . Each sample, which comprised arteries from 2–3 animals, was snap frozen in liquid nitrogen and processed using the RNeasy Micro Kit (Qiagen), including 20 min DNAse treatment. Concentrations of RNA were measured using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies). RNA quality was monitored using RNA 6000 Pico LabChip kit (Agilent 2100 Bioanalyzer). RNA integrity numbers for the 9 samples ranged from 6.4 to 8.1 .
Preparation of cRNA and hybridization on microarrays
Sample preparation followed the Affymetrix GeneChip Eukaryotic Small Sample Target Labeling Assay Version II available at the website . After first round amplification, the yields of unlabeled cRNA ranged from 2.8–7.8 μg. Following the second round amplification, 600 ng of unlabeled cRNA was used for the second round of cDNA synthesis. Double-stranded cDNA was synthesized using oligo(dT) 24-anchored T7 primer, Superscript III (Invitrogen) and 100 ng total RNA for the first strand followed by DNA polymerase I and T4 DNA polymerase for the second strand. The double-stranded cDNA was precipitated with ethanol, air-dried and used as template for the first round of in vitro transcription of cRNA with MEGAscript enzyme mix (Ambion) at 37°C for 7 h. The cRNA was purified (RNeasy, Qiagen) and used as template for the second round of double-stranded cDNA synthesis with random primers (Invitrogen). The double-stranded cDNA product of second round synthesis was purified by ethanol precipitation and used as template for the second round of in vitro transcription and labeling of cRNA, using the ENZO BioArray HighYield RNA Transcript Labeling kit (Affymetrix) according to the manufacturer's protocol. The yields of labelled cRNA ranged from 58–193.8 μg. All components were mixed, incubated at 37°C for 7 h, and the biotin-labeled cRNA (1 μg for each sample) was purified (RNeasy, Qiagen) and fragmented prior to hybridisation to Affymetrix GeneChip RAE230A arrays. Hybridization, washing and staining followed the Affymetrix GeneChip Expression Analysis protocols.
Microarray data analysis
Probeset signal values were calculated using the Statistical Algorithm in the Affymetrix MAS5 software. Raw data were normalised by scaling the global mean method to a TGT value of 150. Quality of all arrays was assessed using MAS5 quality control metrics and software procedures available in R/Bioconductor (Bioconductor version 1.8.0 ) and at the following website . Based on these metrics, the quality of all arrays was considered to be satisfactory (see Additional file 15).
One-way ANOVA was carried out on the log2 transformed MAS5 signal values for each of the 9098 probesets which were present in all replicates of at least one experimental group (Affymetrix MAS5 Statistical Algorithm). Differences between experimental groups were examined using t-tests (post-hoc contrasts).
GE-biplot and Heat Map
An ordination method, the GE-biplot , based on the singular value decomposition of the data matrix of scaled log2 signal values, was obtained using an R function available at the website . The GE-biplot is an unsupervised visualisation tool that facilitates the detection of subsets of genes associated with subsets of microarrays. It provides a genome-wide view of the data, showing genes and microarrays simultaneously. Co-regulated (correlated) genes are located close to each other, with genes having reversed response profiles (negatively correlated) being located on the opposite side of the plot. Genes lying further from the centre have higher variance. Genes which lie close to, or lie in the direction of, a group of samples are relatively upregulated in that group compared with other samples, while those lying in the opposite direction are relatively down-regulated in that group. In a GE-biplot, arrays that have approximately the same relative values of gene expression across probesets cluster together. Further details are described elsewhere .
A heat map of the MAS 5.0 signal values of the 50 most significant probesets as evaluated by one way ANOVA was constructed . Hierarchical clustering with complete linkage was used to order the arrays in the columns of the heat map and the rank order of the p-values from the ANOVA was used to order the columns.
Validation of microarray analysis by quantitative PCR
Quantitative PCR was used to determine the copy numbers of mRNA for 13 genes in the WKY, SHR, and ACTH samples (n = 3 for each). Details of the genes, primers and product sizes are given in Table 5. Nine of the genes were selected on the basis that their expression was significantly altered in either SHR or ACTH and that they were of biological interest to the processes involved in hypertension. The three vascular connexins and eNOS, although unaltered in the microarrays, were chosen because of their significance to vascular function and potential involvement in hypertension . Validation of the microarray data was considered to be a significant (p < 0.05) change in normalised PCR copy number in the same direction, or lack of a change in normalised PCR copy number which corresponded to that found in the microarray data.
Messenger RNA of all 9 saphenous artery samples was reverse transcribed to cDNA (42°C 1 h, 50°C 1 h, 90°C 10 min) using oligo dT (500 ng.μl-1, Invitrogen) primers and Superscript II (200 U.μl-1, Invitrogen). For every sample, control reactions from which either reverse transcriptase or RNA was omitted were run in parallel. Each 25 μl PCR reaction mixture contained 1 μl of each primer (final concentration 800 nM), 2 μl of 12.5 mM dNTPs, 3 μl of 25 mM MgCl2, 0.125 μl of AmpliTaq Gold (5 U. μl-1), 2.5 μl of 10× SYBR Green buffer (SYBR Green Core Reagents kit, Applied Biosystems), 10.375 μl H2O and 5 μl of template containing either 10 ng of cDNA or diluted plasmid DNA containing the relevant cloned PCR product. All samples were diluted in water containing tRNA (2 ng. μl-1). All reactions were performed in duplicate (ABI Prism 7700 Sequence Detection System, ACRF Biomolecular Resource Facility) at 95°C for 10 min, and 40 cycles of 20 sec at 95°C, 20 sec at 65°C, 45 sec at 72°C. Every experiment included duplicate control reactions and the expression of all genes was determined in every sample. The integrity of each PCR reaction was checked using dissociation curve analysis after every reaction.
For each gene, the mRNA concentration was determined by comparing the fluorescent signal at threshold (equal to 10 × S.D. of baseline fluorescence) to that generated by a standard curve using serial dilutions of a known concentration of purified plasmid DNA containing the relevant cloned and sequenced PCR product. Sequences were obtained using an ABI 3730 sequencer (ACRF Biomolecular Resource Facility, JCSMR, ANU). The concentration of each mRNA was normalized by determining 18S rRNA concentration, using primers complementary to a region of the 5' end of the molecule and a standard curve using serial dilutions of plasmid DNA containing the cloned fragment of the 18S rRNA gene . Results were analyzed for statistical significance (p < 0.05) using ANOVA and Student's t-test.
Direct interactions between selected genes
To examine differences in the major pathways affected in SHR and ACTH-induced hypertension, genes that were identified by t-tests as significantly altered in expression in SHR or ACTH compared with WKY (at the 0.015 level) and that were significant in ANOVA at the 0.02 level, were imported into a demonstration version of PathwayStudio 4.0 (Ariadne Genomics, Rockville, MD). This software contains the proprietary ResNet 3.0 database of functional relationships between genes, proteins and small molecules automatically extracted from PubMed and scientific journals using the MedScan natural language processing software. Functional links between genes and proteins were also investigated and cross checked using publicly available gene and genome databases: Online Mendelian Inheritance in Man (OMIM) and Genomic Biology and software provided by NCBI , and the protein database Human Protein Reference Database (HPRD; ) and software .
Hopkins PN, Hunt SC: Genetics of hypertension. Genet Med. 2003, 5 (6): 413-429.
Frohlich ED, Varagic J: Sodium directly impairs target organ function in hypertension. Curr Opin Cardiol. 2005, 20 (5): 424-429. 10.1097/01.hco.0000175519.34933.a5.
Kubo T: [Mechanisms of hypertension in the central nervous system]. Yakugaku Zasshi. 2006, 126 (9): 695-709. 10.1248/yakushi.126.695.
Izzard AS, Rizzoni D, Agabiti-Rosei E, Heagerty AM: Small artery structure and hypertension: adaptive changes and target organ damage. J Hypertens. 2005, 23 (2): 247-250. 10.1097/00004872-200502000-00002.
Laurent S, Boutouyrie P, Lacolley P: Structural and genetic bases of arterial stiffness. Hypertension. 2005, 45 (6): 1050-1055. 10.1161/01.HYP.0000164580.39991.3d.
North KE, MacCluer JW, Devereux RB, Howard BV, Welty TK, Best LG, Lee ET, Fabsitz RR, Roman MJ: Heritability of carotid artery structure and function: the Strong Heart Family Study. Arterioscler Thromb Vasc Biol. 2002, 22 (10): 1698-1703. 10.1161/01.ATV.0000032656.91352.5E.
Lerman LO, Chade AR, Sica V, Napoli C: Animal models of hypertension: an overview. J Lab Clin Med. 2005, 146 (3): 160-173. 10.1016/j.lab.2005.05.005.
Mangos GJ, Turner SW, Fraser TB, Whitworth JA: The role of corticosterone in corticotrophin (ACTH)-induced hypertension in the rat. J Hypertens. 2000, 18 (12): 1849-1855. 10.1097/00004872-200018120-00020.
Connell JM, Whitworth JA, Davies DL, Lever AF, Richards AM, Fraser R: Effects of ACTH and cortisol administration on blood pressure, electrolyte metabolism, atrial natriuretic peptide and renal function in normal man. J Hypertens. 1987, 5 (4): 425-433. 10.1097/00004872-198708000-00007.
Whitworth JA, Schyvens CG, Zhang Y, Andrews MC, Mangos GJ, Kelly JJ: The nitric oxide system in glucocorticoid-induced hypertension. J Hypertens. 2002, 20 (6): 1035-1043. 10.1097/00004872-200206000-00003.
Zhang Y, Chan MM, Andrews MC, Mori TA, Croft KD, McKenzie KU, Schyvens CG, Whitworth JA: Apocynin but not allopurinol prevents and reverses adrenocorticotropic hormone-induced hypertension in the rat. Am J Hypertens. 2005, 18 (7): 910-916. 10.1016/j.amjhyper.2005.02.017.
Ulker S, McKeown PP, Bayraktutan U: Vitamins reverse endothelial dysfunction through regulation of eNOS and NAD(P)H oxidase activities. Hypertension. 2003, 41 (3): 534-539. 10.1161/01.HYP.0000057421.28533.37.
Demirci B, McKeown PP, Bayraktutan U: Blockade of angiotensin II provides additional benefits in hypertension- and ageing-related cardiac and vascular dysfunctions beyond its blood pressure-lowering effects. J Hypertens. 2005, 23 (12): 2219-2227. 10.1097/01.hjh.0000191906.03983.ee.
Mulvany MJ, Baumbach GL, Aalkjaer C, Heagerty AM, Korsgaard N, Schiffrin EL, Heistad DD: Vascular remodeling. Hypertension. 1996, 28 (3): 505-506.
Arribas SM, Hillier C, Gonzalez C, McGrory S, Dominiczak AF, McGrath JC: Cellular aspects of vascular remodeling in hypertension revealed by confocal microscopy. Hypertension. 1997, 30 (6): 1455-1464.
Huttner I, Costabella PM, De Chastonay C, Gabbiani G: Volume, surface, and junctions of rat aortic endothelium during experimental hypertension: a morphometric and freeze fracture study. Lab Invest. 1982, 46 (5): 489-504.
Rummery NM, McKenzie KU, Whitworth JA, Hill CE: Decreased endothelial size and connexin expression in rat caudal arteries during hypertension. J Hypertens. 2002, 20 (2): 247-253. 10.1097/00004872-200202000-00014.
Schiffrin EL, Touyz RM: From bedside to bench to bedside: role of renin-angiotensin-aldosterone system in remodeling of resistance arteries in hypertension. Am J Physiol Heart Circ Physiol. 2004, 287 (2): H435-46. 10.1152/ajpheart.00262.2004.
Touyz RM: Intracellular mechanisms involved in vascular remodelling of resistance arteries in hypertension: role of angiotensin II. Exp Physiol. 2005, 90 (4): 449-455. 10.1113/expphysiol.2005.030080.
Inoue T, Masuda T, Kishi K: Structural and functional alterations of mesenteric vascular beds in spontaneously hypertensive rats. Jpn Heart J. 1990, 31 (3): 393-403.
Intengan HD, Thibault G, Li JS, Schiffrin EL: Resistance artery mechanics, structure, and extracellular components in spontaneously hypertensive rats : effects of angiotensin receptor antagonism and converting enzyme inhibition. Circulation. 1999, 100 (22): 2267-2275.
Sandow SL, Bramich NJ, Bandi HP, Rummery NM, Hill CE: Structure, function, and endothelium-derived hyperpolarizing factor in the caudal artery of the SHR and WKY rat. Arterioscler Thromb Vasc Biol. 2003, 23 (5): 822-828. 10.1161/01.ATV.0000067425.06019.D7.
Yang L, Gao YJ, Lee RM: Quinapril effects on resistance artery structure and function in hypertension. Naunyn Schmiedebergs Arch Pharmacol. 2004, 370 (6): 444-451. 10.1007/s00210-004-0990-x.
Westhoff TH, Scheid S, Tolle M, Kaynak B, Schmidt S, Zidek W, Sperling S, van der Giet M: A physiogenomic approach to study the regulation of blood pressure. Physiol Genomics. 2005, 23 (1): 46-53. 10.1152/physiolgenomics.00077.2005.
NCBI Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
Choe SE, Boutros M, Michelson AM, Church GM, Halfon MS: Preferred analysis methods for Affymetrix GeneChips revealed by a wholly defined control dataset. Genome Biol. 2005, 6 (2): R16-10.1186/gb-2005-6-2-r16.
Hnasko R, Lisanti MP: The biology of caveolae: lessons from caveolin knockout mice and implications for human disease. Mol Interv. 2003, 3 (8): 445-464. 10.1124/mi.3.8.445.
Wieland T, Mittmann C: Regulators of G-protein signalling: multifunctional proteins with impact on signalling in the cardiovascular system. Pharmacol Ther. 2003, 97 (2): 95-115. 10.1016/S0163-7258(02)00326-1.
Rummery NM, Hickey H, McGurk G, Hill CE: Connexin37 is the major connexin expressed in the media of caudal artery. Arterioscler Thromb Vasc Biol. 2002, 22 (9): 1427-1432. 10.1161/01.ATV.0000028814.45706.E5.
Goto K, Rummery NM, Grayson TH, Hill CE: Attenuation of conducted vasodilatation in rat mesenteric arteries during hypertension: role of inwardly rectifying potassium channels. J Physiol. 2004, 561 (Pt 1): 215-231. 10.1113/jphysiol.2004.070458.
Rummery NM, Grayson TH, Hill CE: Angiotensin-converting enzyme inhibition restores endothelial but not medial connexin expression in hypertensive rats. J Hypertens. 2005, 23 (2): 317-328. 10.1097/00004872-200502000-00014.
Papparella I, Ceolotto G, Lenzini L, Mazzoni M, Franco L, Sartori M, Ciccariello L, Semplicini A: Angiotensin II-induced over-activation of p47phox in fibroblasts from hypertensives: which role in the enhanced ERK1/2 responsiveness to angiotensin II?. J Hypertens. 2005, 23 (4): 793-800. 10.1097/01.hjh.0000163148.97459.9d.
Ono K, Han J: The p38 signal transduction pathway: activation and function. Cell Signal. 2000, 12 (1): 1-13. 10.1016/S0898-6568(99)00071-6.
Safar ME: Systolic hypertension in the elderly: arterial wall mechanical properties and the renin-angiotensin-aldosterone system. J Hypertens. 2005, 23 (4): 673-681. 10.1097/01.hjh.0000163130.39149.fe.
Kamm KE, Stull JT: Dedicated myosin light chain kinases with diverse cellular functions. J Biol Chem. 2001, 276 (7): 4527-4530. 10.1074/jbc.R000028200.
Drab M, Verkade P, Elger M, Kasper M, Lohn M, Lauterbach B, Menne J, Lindschau C, Mende F, Luft FC, Schedl A, Haller H, Kurzchalia TV: Loss of caveolae, vascular dysfunction, and pulmonary defects in caveolin-1 gene-disrupted mice. Science. 2001, 293 (5539): 2449-2452. 10.1126/science.1062688.
Gratton JP, Bernatchez P, Sessa WC: Caveolae and caveolins in the cardiovascular system. Circ Res. 2004, 94 (11): 1408-1417. 10.1161/01.RES.0000129178.56294.17.
Patel HH, Zhang S, Murray F, Suda RY, Head BP, Yokoyama U, Swaney JS, Niesman IR, Schermuly RT, Pullamsetti SS, Thistlethwaite PA, Miyanohara A, Farquhar MG, Yuan JX, Insel PA: Increased smooth muscle cell expression of caveolin-1 and caveolae contribute to the pathophysiology of idiopathic pulmonary arterial hypertension. Faseb J. 2007, 21 (11): 2970-2979. 10.1096/fj.07-8424com.
Heximer SP, Knutsen RH, Sun X, Kaltenbronn KM, Rhee MH, Peng N, Oliveira-dos-Santos A, Penninger JM, Muslin AJ, Steinberg TH, Wyss JM, Mecham RP, Blumer KJ: Hypertension and prolonged vasoconstrictor signaling in RGS2-deficient mice. J Clin Invest. 2003, 111 (4): 445-452. 10.1172/JCI200315598.
Tang KM, Wang GR, Lu P, Karas RH, Aronovitz M, Heximer SP, Kaltenbronn KM, Blumer KJ, Siderovski DP, Zhu Y, Mendelsohn ME: Regulator of G-protein signaling-2 mediates vascular smooth muscle relaxation and blood pressure. Nat Med. 2003, 9 (12): 1506-1512. 10.1038/nm958.
Sun X, Kaltenbronn KM, Steinberg TH, Blumer KJ: RGS2 is a mediator of nitric oxide action on blood pressure and vasoconstrictor signaling. Mol Pharmacol. 2005, 67 (3): 631-639. 10.1124/mol.104.007724.
Dupuis M, Soubrier F, Brocheriou I, Raoux S, Haloui M, Louedec L, Michel JB, Nadaud S: Profiling of aortic smooth muscle cell gene expression in response to chronic inhibition of nitric oxide synthase in rats. Circulation. 2004, 110 (7): 867-873. 10.1161/01.CIR.0000138850.72900.FE.
Adams LD, Geary RL, McManus B, Schwartz SM: A comparison of aorta and vena cava medial message expression by cDNA array analysis identifies a set of 68 consistently differentially expressed genes, all in aortic media. Circ Res. 2000, 87 (7): 623-631.
Li J, Adams LD, Wang X, Pabon L, Schwartz SM, Sane DC, Geary RL: Regulator of G protein signaling 5 marks peripheral arterial smooth muscle cells and is downregulated in atherosclerotic plaque. J Vasc Surg. 2004, 40 (3): 519-528. 10.1016/j.jvs.2004.06.021.
Berger M, Bergers G, Arnold B, Hammerling GJ, Ganss R: Regulator of G-protein signaling-5 induction in pericytes coincides with active vessel remodeling during neovascularization. Blood. 2005, 105 (3): 1094-1101. 10.1182/blood-2004-06-2315.
Kirsch T, Wellner M, Luft FC, Haller H, Lippoldt A: Altered gene expression in cerebral capillaries of stroke-prone spontaneously hypertensive rats. Brain Res. 2001, 910 (1-2): 106-115. 10.1016/S0006-8993(01)02670-1.
Christensen KL, Mulvany MJ: Location of resistance arteries. J Vasc Res. 2001, 38 (1): 1-12. 10.1159/000051024.
Sandow SL, Tare M, Coleman HA, Hill CE, Parkington HC: Involvement of myoendothelial gap junctions in the actions of endothelium-derived hyperpolarizing factor. Circ Res. 2002, 90 (10): 1108-1113. 10.1161/01.RES.0000019756.88731.83.
Sandow SL, Goto K, Rummery NM, Hill CE: Developmental changes in myoendothelial gap junction mediated vasodilator activity in the rat saphenous artery. J Physiol. 2004, 556 (Pt 3): 875-886. 10.1113/jphysiol.2003.058669.
Wit E, McClure J: Statistics for Microarrays: Design, Analysis and Inference. 2004, John Wiley and Sons, Ltd.
Imbeaud S, Graudens E, Boulanger V, Barlet X, Zaborski P, Eveno E, Mueller O, Schroeder A, Auffray C: Towards standardization of RNA quality assessment using user-independent classifiers of microcapillary electrophoresis traces. Nucleic Acids Res. 2005, 33 (6): e56-10.1093/nar/gni054.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
The Centre for Bioinformation Science. [http://dayhoff.anu.edu.au/software.html]
Pittelkow YE, Wilson SR: Visualisation of gene expression data - the GE-biplot, the Chip-plot and the Gene-plot. Stat Appl Genet Mol Biol. 2003, 2: Article6-
Dudoit S, Gentleman RC, Quackenbush J: Open source software for the analysis of microarray data. Biotechniques. 2003, Suppl: 45-51.
Rummery NM, Hill CE: Vascular gap junctions and implications for hypertension. Clin Exp Pharmacol Physiol. 2004, 31 (10): 659-667. 10.1111/j.1440-1681.2004.04071.x.
National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
Human Protein Reference Database. [http://www.hprd.org]
Peri S, Navarro JD, Amanchy R, Kristiansen TZ, Jonnalagadda CK, Surendranath V, Niranjan V, Muthusamy B, Gandhi TK, Gronborg M, Ibarrola N, Deshpande N, Shanker K, Shivashankar HN, Rashmi BP, Ramya MA, Zhao Z, Chandrika KN, Padma N, Harsha HC, Yatish AJ, Kavitha MP, Menezes M, Choudhury DR, Suresh S, Ghosh N, Saravana R, Chandran S, Krishna S, Joy M, Anand SK, Madavan V, Joseph A, Wong GW, Schiemann WP, Constantinescu SN, Huang L, Khosravi-Far R, Steen H, Tewari M, Ghaffari S, Blobe GC, Dang CV, Garcia JG, Pevsner J, Jensen ON, Roepstorff P, Deshpande KS, Chinnaiyan AM, Hamosh A, Chakravarti A, Pandey A: Development of human protein reference database as an initial platform for approaching systems biology in humans. Genome Res. 2003, 13 (10): 2363-2371. 10.1101/gr.1680803.
The study was funded by the National Heart Foundation and National Health and Medical Research Council of Australia (NHMRC) Project Grant 224219. SLS was supported by NHMRC RD Wright Fellowship, 455342.
THG extracted RNA, performed quantitative PCR, contributed to the analysis and interpretation of microarray data, analysed biological association networks and interactions between genes and proteins, and drafted the manuscript. SJO performed quality assessment of the microarrays, contributed to the design and statistical analyses of the microarray data, prepared data for the figures, participated in interpretation of the results and contributed to writing the manuscript. TDB administered ACTH, measured blood pressure of experimental rat groups and performed serial section electron microscopy and morphometric measurements. KRM performed the immunohistochemistry and measured blood pressure of experimental animals. KP amplified and labeled mRNA and conducted all of the microarray hybridizations. YEP and SRW designed and performed statistical analyses of data and contributed to the interpretation of results and preparation of the manuscript. SLS performed serial section electron microscopy and morphometric measurements, assisted with experimental design and contributed to preparation of the manuscript. CEH conceived of the study and coordinated the experimental design, dissected arteries, contributed to the interpretation of data and made a substantial contribution to the writing of the manuscript. All authors have read and approved submission of the manuscript.
Electronic supplementary material
Additional file 4: Genes with decreased expression in SHR vs WKY and unchanged in ACTH vs WKY. Description of dataset: Gene lists of statistical analysis of microarray data including gene names, gene descriptions, biological function, Affymetrix IDs, accession numbers, expression levels and p-values. (XLS 73 KB)
Additional file 7: Genes with increased expression in ACTH vs WKY and decreased expression in SHR vs WKY. Gene lists of statistical analysis of microarray data including gene names, gene descriptions, biological function, Affymetrix IDs, accession numbers, expression levels and p-values. (XLS 24 KB)
Additional file 8: Genes with decreased expression in ACTH vs WKY and increased expression in SHR vs WKY. Gene lists of statistical analysis of microarray data including gene names, gene descriptions, biological function, Affymetrix IDs, accession numbers, expression levels and p-values. (XLS 24 KB)
Authors’ original submitted files for images
About this article
Cite this article
Grayson, T.H., Ohms, S.J., Brackenbury, T.D. et al. Vascular microarray profiling in two models of hypertension identifies caveolin-1, Rgs2 and Rgs5 as antihypertensive targets. BMC Genomics 8, 404 (2007). https://doi.org/10.1186/1471-2164-8-404
- Myosin Light Chain Kinase
- Hypertensive Model
- Smooth Muscle Cell Layer
- Saphenous Artery
- p38MAP Kinase Pathway