(Hypertension. 2003;41:847.)
© 2003 American Heart Association, Inc.
Scientific Contributions |
From the BHF Glasgow Cardiovascular Research Centre, Division of Cardiovascular and Medical Sciences, University of Glasgow, Western Infirmary, Glasgow, Scotland.
Correspondence to Prof Anna F. Dominiczak, BHF Glasgow Cardiovascular Research Centre, Division of Cardiovascular and Medical Sciences, University of Glasgow, Western Infirmary, Glasgow, G11 6NT, Scotland. E-mail ad7e{at}clinmed.gla.ac.uk
| Abstract |
|---|
|
|
|---|
Key Words: genetics hypertension, genetic gene expression rats, stroke-prone SHR
| Introduction |
|---|
|
|
|---|
The stroke-prone spontaneously hypertensive rat (SHRSP) is a well-characterized experimental model for human essential hypertension. Similar to human disease, the genetic determination of blood pressure variation in this model is complex and due to multiple gene-gene and gene-environment interactions. The SHRSP develops a number of vascular complications, including cardiac hypertrophy, cardiac failure, and stroke.24 Genome-wide linkage studies have proved successful in the localization of large chromosomal regions containing quantitative trait loci (QTLs) for blood regulation in the SHRSP.5 In particular, previous work in our laboratory has identified at least 2 blood pressure QTLs mapping to rat chromosome 2.6 This region of chromosome 2 is a classic example of a common or overlapping QTL, as it has been implicated in several crosses.7 Moreover, the QTL contained in the telomeric region of rat chromosome 2 has been confirmed in 2 different congenic strains by introgressing the relevant region from Milan normotensive rat or Wistar Kyoto rat (WKY) into the Dahl salt-sensitive background.810
Congenic strains derived from hypertensive rat models have been used to confirm the existence of several other QTLs, and 1 locus has already been fine-mapped to a 0.54-centimorgan (cM) interval.11 However, the requirement for the construction of minimal congenic strains may be circumvented, at least in some cases, by use of new technology such as high-throughput microarray expression profiling. Microarrays are available in 2 principle forms: complementary DNA (cDNA)12 and oligonucleotides.13 Microarray technology has been used to accelerate gene identification by Aitman et al,14 who used cDNA microarrays and rat chromosome 4 congenic strains to identify Cd36, a gene responsible for defective fatty acid metabolism in the spontaneously hypertensive rat (SHR).
Several mapping studies investigating genome conservation at the genetic and physical level have determined that gene order is relatively invariant showing conserved synteny between mammals.15 Comparative genome analysis of rat, mouse, and human can be used to identify regions of the human genome likely to harbor genes involved in blood pressure regulation.7,16 Julier et al17 and Baima et al18 provided first successful examples of comparative mapping in which a region involved in blood pressure regulation on rat chromosome 10 indicated a susceptibility locus for human hypertension on human chromosome 17q.
Aims of the current study were to generate and phenotype a congenic strain containing a region of rat chromosome 2 from the normotensive WKY strain introgressed into the genetic background of the SHRSP. Genome-wide expression profiling was undertaken to identify differentially expressed genes among the parental SHRSP, WKY, and congenic strain. Furthermore, a comparative genome analysis was used to identify regions of conserved synteny on mouse and human chromosomes.
| Methods |
|---|
|
|
|---|
The Dataquest IV telemetry system (Data Sciences International) was used for the direct measurement of systolic and diastolic blood pressure.19,20 Briefly, rats were implanted at 12 weeks of age with 1-week recovery followed by 3 weeks of baseline measurements. Rats were killed at the end of the 15th week, and kidneys, spleen, and liver were snap-frozen in liquid nitrogen and stored at -70°C for either RNA or DNA extraction. Genomic DNA was isolated from a 4-mm tip from the tail of congenic animals and genotyping as we described previously.20 These studies were approved by the Home Office according to regulations regarding experiments with animals in the United Kingdom.
Genetic Linkage and Radiation Hybrid Maps
Genotyping was performed by polymerase chain reaction (PCR) amplification of microsatellite markers, and the genotypic results obtained were mapped relative to each other using the MAPMAKER/EXP 3.0 computer package with an error-detection procedure.21 We constructed a genetic linkage map of rat chromosome 2 consisting of 74 markers polymorphic between the WKY and SHRSP strains. Genetic distances were calculated with the Haldane mapping function. PCR was performed using the rat genome T55 radiation hybrid (RH) panel obtained from Research Genetics. Scored consensus data were submitted to the RHMAPPER program (http://rgd.mcw.edu/RHMAPSERVER/).
Gene Expression Profiling
Affymetrix GeneChip expression analysis was used to identify differentially expressed probe sets (representing unique gene or expressed sequence tag sequence on the Affymetrix GeneChip) between SHRSP, SP.WKYGla2c*, and WKY. Whole kidneys were homogenized and total RNA extracted from 3 rats from each strain by using the maxi RNeasy kit according to manufacturers protocol (Qiagen). Biotinylated amplified target cRNA was prepared and hybridized to the Affymetrix Rat U34 array set (chips U34A, U34B, and U34C; a total of 26 379 probe sets) as described by Affymetrix.22 After hybridization, microarray chips were washed, stained, and scanned. Using MAS 5.0 (Affymetrix) and each of the raw data probe sets, a perfect match versus mismatch value was calculated, and the data were compiled into one value for each probe set. Statistical tools available from Genespring 4.1 (Silicon Genetics) were used to normalize the chip data and identify differentially expressed probe sets among SHRSP, SP.WKYGla2c*, and WKY. Probe sets were identified as being differentially expressed by using the global error model, a variance components model. The final step in the filtering process was to remove genes and expressed sequence tags not expressed in either of the strains being compared.
Quantitative Reverse TranscriptionPCR
RNA was accurately quantified using Ribogreen (Molecular Probes) and normalized to 285 ng/µL. Normalization was confirmed by performing real-time reverse transcription (RT)PCR on the LightCylcer (Roche) of ß-actin (Promega) with comparable threshold cycles. RT-PCR was performed on 3 individual animals from each strain with 0.5 µmol/L each primer for glutathione S-transferase µ-type 2 (Gstm2) F 5'-TTT GAG CCC AAG TGC CTG GA-3' R 5'-GCAGGATCCAATGTGGACAG-3' and then expressed relative to the ß-actin standard curve generated by using serial dilutions (100, 1x10-1, 5x10-2, 1x10-2). Melting curve analysis was performed to confirm the presence of a single species. Concentrations of gene transcripts were calculated based on this standard curve.
Comparative Genome Analysis
Comparative genome analysis used publically available rat, mouse, and human data to identify genes with interspecies sequence identity (homology) for comparative map construction. The anchors for the comparative map are genes that have been mapped in all 3 species. Only genes that could be confirmed as mapping correctly in at least 2 of the species were added to the map.
Statistical Analysis
Blood pressure comparisons between the congenic and the parental strains were performed as previously described.9 Microarray data used the global error model, which is a variance components method that estimates both measurement and sample-to-sample variation and standard errors to compare mean expression levels between groups of samples. A nominal significance level of P=0.05 was taken to be the threshold for differential expression. Statistical analysis of RT-PCR by LightCycler was undertaken using the Student t test.
| Results |
|---|
|
|
|---|
|
|
Construction of Genetic and RH Maps
To create the new genetic linkage map, a total of 74 microsatellite markers, spanning rat chromosome 2, were genotyped in F2 animals. Of the 74 markers scored, we positioned 55 on the RH map of rat chromosome 2. Nineteen of the markers added to the map did not position within the region of rat chromosome 2, as they were multilinked and the reduction of stringency to accommodate them would result in a less accurate map. The markers covered a distance of 1224 centiRays and were all placed with a logarithm of odds ratio (LOD) >10. This equates to 1010 odds for linkage.
Gene Expression Profiling
Total numbers of differentially expressed genes between the SHRSP and SP.WKYGla2c*, as well as the SHRSP and WKY, are shown in Table 1. Of the 45 differentially expressed probe sets in the SHRSP and SP.WKYGla2c* comparison, 12 were also identified in the SHRSP and WKY comparison (Table 2). Four of these probe sets have a known rat chromosomal location, with 3 mapping to the congenic segment. All 3 probe sets that mapped to the congenic segment represented the same gene, glutathione S-transferase µ-type 2 (Gstm2). Of these, two were gene probe sets, and one was an expressed sequence tag DNA sequence. Quantitative RT-PCR relative to the same ß-actin standard curve further confirmed the microarray results (Figure 3). The LightCycler showed Gstm2 expression was significantly lower in SHRSP (8.56x10-4±1.6x10-4) versus SP.WKYGla2c* (3.67x10-3±2.8x10-4); (95% CI -3.9x10-3 to -1.8x10-3; P=0.0034) and versus WKY (4.03x10-3 ±5.1x10-4); (95% CI -5.4x10-3 to -8.9x10-4; P=0.027).
|
|
|
Comparative Genome Analysis
We used publically available databases to demonstrate that the SP.WKYGla2c* congenic segment showed highly conserved synteny to mouse chromosome 3 and human chromosome 1. Orthologous genes between the 3 species are shown in Figure 4. Gstm2 is one of the orthologs, as it maps to the region of conserved gene synteny in all 3 species.
|
| Discussion |
|---|
|
|
|---|
Other efforts to analyze the QTLs on rat chromosome 2 have focused on the construction of overlapping congenic strains to narrow down the region(s) under study.810 However, these studies have not resulted in the identification of positional candidate genes. Although theoretical considerations have suggested that it should be possible to obtain minimal congenic strains as small as 1 cM,5 this has been proven very difficult in practice. In the current study, we produced dense genetic and RH maps of rat chromosome 2 in an effort to narrow down the congenic region. Similar work has been performed for the remaining rat chromosomes by other groups.30,31 However, it seems unlikely that minimal congenic strains encompassing <3 cM will be frequently produced. Therefore, ancillary strategies such as microarray gene expression profiling might contribute to causative gene discovery. Aitman et al14 and Eaves et al32 used rat and mouse congenic strains, respectively, in combination with gene chip technology. The former study14 identified a strong positional and physiological candidate gene, the Cd36, although the congenic interval of the SHR.BN strain used was relatively large at 36 cM. However, Eaves et al32 failed to identify any strong candidates, despite a true minimal congenic strain at 0.35 cM. These contrasting results, together with the microarray data from the current study, confirm that although not every congenic/microarray experiment will identify causative genes, the size of congenic interval should not be seen as a limitation for a well-designed microarray experiment.
The choice of tissue for gene expression profiling has been a matter of recent debate.33 It seems that for gene hunting experiments in hypertension, the kidney is an ideal choice as several elegant transplantation experiments showed that hypertension always "travels with the kidney."3436 The other important issue, already stressed by Aitman et al14 but also clearly seen in the current study, is the relative ease of interpretation of microarray data. Although one would expect hundreds of differentially expressed genes between the 2 parental strains, this number is significantly reduced in the congenic and parental strain comparisons. Moreover, initially one would focus on genes that are differentially expressed and mapped to the congenic region, with a possibility of later analysis of genes mapping to other chromosomes but perhaps still of functional importance through the relevant physiological pathways.37
Perspectives
We have identified Gstm2 as a positional and physiological candidate gene for blood pressure regulation and oxidative stress in the SHRSP. Using comparative genome analysis between rat, mouse, and human, we successfully transferred the relevant QTL across species and showed that the rat Gstm2 has orthologs on mouse chromosome 3 and human chromosome 1. Despite highly significant changes in expression between the congenic and the SHRSP strain, we have not excluded the possibility that these could be secondary to blood pressure differences. Further studies in young animals during the development of hypertension may help elucidate primary mechanisms.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received October 4, 2002; first decision October 24, 2002; accepted November 4, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. H. Lee, B. J. Haas, N. E. Letwin, B. C. Frank, T. V. Luu, Q. Sun, C. D. House, S. Yerga-Woolwine, P. Farms, E. Manickavasagam, et al. Cross-Talk of Expression Quantitative Trait Loci Within 2 Interacting Blood Pressure Quantitative Trait Loci Hypertension, December 1, 2007; 50(6): 1126 - 1133. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Graham, M. W. McBride, M. Gaasenbeek, K. Gilday, E. Beattie, W. H. Miller, J. D. McClure, J. M. Polke, A. Montezano, R. M. Touyz, et al. Candidate Genes That Determine Response to Salt in the Stroke-Prone Spontaneously Hypertensive Rat: Congenic Analysis Hypertension, December 1, 2007; 50(6): 1134 - 1141. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. H. Lee Physiogenomic strategies and resources to associate genes with rat models of heart, lung and blood disorders Exp Physiol, November 1, 2007; 92(6): 992 - 1002. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Y. Deng Genetic basis of polygenic hypertension Hum. Mol. Genet., October 15, 2007; 16(R2): R195 - R202. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Sancho-Shimizu, R. Khan, S. Mostowy, L. Lariviere, R. Wilkinson, N. Riendeau, M. Behr, and D. Malo Molecular Genetic Analysis of Two Loci (Ity2 and Ity3) Involved in the Host Response to Infection With Salmonella Typhimurium Using Congenic Mice and Expression Profiling Genetics, October 1, 2007; 177(2): 1125 - 1139. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. R. Garrett, W. T. Gunning, T. Radecki, and A. Richard Dissection of a genetic locus influencing renal function in the rat and its concordance with kidney disease loci on human chromosome 1q21 Physiol Genomics, August 20, 2007; 30(3): 322 - 334. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Y. Deng Positional Cloning of Quantitative Trait Loci for Blood Pressure: How Close Are We?: A Critical Perspective Hypertension, April 1, 2007; 49(4): 740 - 747. [Full Text] [PDF] |
||||
![]() |
L. J. Mullins, M. A. Bailey, and J. J. Mullins Hypertension, Kidney, and Transgenics: A Fresh Perspective Physiol Rev, April 1, 2006; 86(2): 709 - 746. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Joe, N. E. Letwin, M. R. Garrett, S. Dhindaw, B. Frank, R. Sultana, K. Verratti, J. P. Rapp, and N. H. Lee Transcriptional profiling with a blood pressure QTL interval-specific oligonucleotide array Physiol Genomics, November 17, 2005; 23(3): 318 - 326. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Seubert, F. Xu, J. P. Graves, J. B. Collins, S. O. Sieber, R. S. Paules, D. L. Kroetz, and D. C. Zeldin Differential renal gene expression in prehypertensive and hypertensive spontaneously hypertensive rats Am J Physiol Renal Physiol, September 1, 2005; 289(3): F552 - F561. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. W. McBride, M. J. Brosnan, J. Mathers, L. I. McLellan, W. H. Miller, D. Graham, N. Hanlon, C. A. Hamilton, J. M. Polke, W. K. Lee, et al. Reduction of Gstm1 Expression in the Stroke-Prone Spontaneously Hypertension Rat Contributes to Increased Oxidative Stress Hypertension, April 1, 2005; 45(4): 786 - 792. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. F. Dominiczak, D. Graham, M. W. McBride, N. J.R. Brain, W. K. Lee, F. J. Charchar, M. Tomaszewski, C. Delles, and C. A. Hamilton Cardiovascular Genomics and Oxidative Stress Hypertension, April 1, 2005; 45(4): 636 - 642. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Hinojos, E. Boerwinkle, M. Fornage, and P. A. Doris Combined Genealogical, Mapping, and Expression Approaches to Identify Spontaneously Hypertensive Rat Hypertension Candidate Genes Hypertension, April 1, 2005; 45(4): 698 - 704. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. J. Wallace, R. H. Wallis, S. C. Collins, K. Argoud, P. J. Kaisaki, A. Ktorza, P. Y. Woon, M.-T. Bihoreau, and D. Gauguier Quantitative trait locus dissection in congenic strains of the Goto-Kakizaki rat identifies a region conserved with diabetes loci in human chromosome 1q Physiol Genomics, September 16, 2004; 19(1): 1 - 10. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Fries, P. Mahboubi, N. R. Mahapatra, S. K. Mahata, N. J. Schork, G. W. Schmid-Schoenbein, and D. T. O'Connor Neuroendocrine Transcriptome in Genetic Hypertension: Multiple Changes in Diverse Adrenal Physiological Systems Hypertension, June 1, 2004; 43(6): 1301 - 1311. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Mein, M. J. Caulfield, R. J. Dobson, and P. B. Munroe Genetics of essential hypertension Hum. Mol. Genet., April 1, 2004; 13(90001): R169 - 175. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. W. McBride, F. J. Charchar, D. Graham, W. H. Miller, P. Strahorn, F. J. Carr, and A. F. Dominiczak Functional genomics in rodent models of hypertension J. Physiol., January 1, 2004; 554(1): 56 - 63. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||