| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2002;105:1824.)
© 2002 American Heart Association, Inc.
Basic Science Reports |
From the Department of Medicine and Genetics, Duke University Medical Center, Durham, NC.
Correspondence to Howard A. Rockman, MD, DUMC 3104, Durham, NC 27710. E-mail h.rockman{at}duke.edu
| Abstract |
|---|
|
|
|---|
Methods and Results In this study, we used transgenic mice that develop severe dilated cardiomyopathy due to the cardiac-specific overexpression of calsequestrin. Transgenic mice showed marked strain-specific variation of cardiac function and survival, independent of transgene expression. A reciprocal backcross strategy was employed using two inbred strains showing distinct differences in survival and cardiac function. To map the genes that modified the heart failure phenotype, progeny from the 2 reciprocal backcrosses were used in a genome-wide scan for linkage. We identified two loci significantly linked to survival with a maximum likelihood ratio statistic of 36.2 (LOD score
7.8) on chromosome 2 and of 26.5 (LOD score
5.7) on chromosome 3. The chromosome 3 locus was also significantly linked to cardiac function with a maximum likelihood ratio statistic of 42.9 (LOD score
9.3). Because only a single strong modifier locus was found in each backcross, we applied a haplotype analysis to map crossovers and successfully narrowed the critical intervals for each locus.
Conclusion Using a sensitized mouse model, we identified major modifier loci that affect the genetically complex disease of heart failure. This approach should allow the rapid identification of candidate genes involved in disease susceptibility in human populations and new insights into the pathogenesis of heart failure.
Key Words: heart failure genes linkage mapping transgenic mice genetics
| Introduction |
|---|
|
|
|---|
Although rodent models have been successfully used to map disease susceptibility loci for complex disease, these studies have often resulted in the identification of multiple loci of weak to moderate effect, making it difficult to identify the genes involved.5 In order to circumvent these difficulties, we sought to simplify the genetic contribution to an otherwise complex phenotype by using a mouse model of heart failure as a "disease-sensitized" strain. The mouse model was generated by the cardiac overexpression of the sarcoplasmic reticulum Ca2+-binding protein calsequestrin (CSQ),6,7 which sequesters the contractile-dependent pool of Ca2+.8 CSQ transgenic mice recapitulate many of the phenotypic features of human dilated cardiomyopathy, including cardiac enlargement, depressed contractile function, abnormal ß-adrenergic receptor signaling, and premature death.6,7,9 However, the effects of the CSQ transgene are highly strain-specific, showing wide variability in phenotype when crossed into different inbred mouse strains. The goal of this study was to use a genetically sensitized mouse model to map modifier genes that affect the heart failure phenotype.
| Methods |
|---|
|
|
|---|
Echocardiography
Transthoracic M-mode echocardiography was performed on all mice at the age of 4 weeks in the conscious state using a HDI 5000CV echocardiograph machine with a 10.5-MHz frequency probe (ATL Ultrasound).10
Western Blotting
Whole hearts from wild-type DBA/B6, DBA CSQ, F1 (B6xDBA) CSQ, N2 (F1xB6) CSQ, and N2 (F1xDBA) CSQ mice were homogenized in lysis buffer, and 15 µg of cytosolic extracts were used for immunoblotting with a monoclonal CSQ antibody (kindly provided by Dr Larry Jones, Krannert Institute, Indiana University, Indianapolis).
Genome-Wide Scan
Seventy progeny from the N2 (F1xB6) and 56 progeny from the N2 (F1xDBA), chosen because their phenotypic values were greater than 0.5 standard deviations from the mean, were screened using 100 microsatellite markers (Research Genetics) spaced at approximately 15-cM intervals through the 19 mouse autosomes (http://www-genome.wi.mit.edu). An additional 39 progeny from the N2 (F1xB6) and 21 progeny from the N2 (F1xDBA) chosen at random were included in the analysis to follow-up on regions with evidence of a QTL. DNA was visualized by one of the following protocols: fluorescent labeling using an ABI Prism 377 sequencer (Applied Biosystems), ethidium bromide staining on agarose gels, or radioactive labeling on acrylamide gels. Data were collected with GeneScan 2.1 software (Applied Biosystems).
Statistical Analysis
MapManager QTX11 was used to perform all the statistical functions used in QTL mapping, including the estimation of the phenotypic variability due to each locus identified. In order to establish linkage thresholds, permutation tests12 were performed on our data at 1-cM intervals for 10 000 permutations. Suggestive (P=0.67), significant (P=0.05), and highly significant (P=0.001) thresholds were established based on the guidelines suggested by Lander and Kruglyak.13 Single locus association tests14 were performed between each marker and the phenotype to identify regions of interest in our initial genome scan. Chromosomes containing markers that reached or exceeded the established genome-wide significant linkage threshold (P<0.05) were analyzed further using simple interval mapping at 1-cM intervals.14 The logarithm of the odds (LOD) score equivalent was calculated by dividing the likelihood ratio statistic (LRS) value by 4.6. Genetic distances were based on those reported in the Mouse Genome Database (MGD) (http://www.informatics.jax.org/).
Haplotype Analysis
We selected animals that had a recombinant haplotype in our regions of interest (chromosome 2 between markers D2Mit91 andD2Mit526 for the F1 backcross to B6 and chromosome 3 between markers D3Mit42 and D3Mit116 for the F1 backcross to DBA). Animals that were included in the haplotype analysis had a phenotypic value for survival beyond a 0.5 standard deviation from the mean (survival for the N2 (F1xDBA) mice was 137.6±46.7 and the N2 (F1xB6) mice was 33.5±5.3 days).
| Results |
|---|
|
|
|---|
|
To further assess the effect of genetic background on cardiac function, we backcrossed F1 CSQ males to either wild-type DBA or B6 females to generate N2 progeny. Because the F1 CSQ mice do not survive longer than 8 weeks, breeding F1 CSQ females to obtain the next generation was not feasible. Backcrossing the F1 CSQ mice with B6 led to a further shortening of survival and worsening in cardiac function in the N2 (F1xB6) progeny (Figures 2A through 2C). In contrast, the range in survival of the N2 (F1xDBA) CSQ progeny was broad and spread over 431 days with a mean of 137.6±93.4 days (Figure 2A).
|
Although we were unable to obtain echocardiographic data on 33 out of 94 (39.4%) N2 (F1xB6) CSQ progeny because they died before their 4-week echocardiogram, echocardiograms on the remaining 61 N2 (F1xB6) CSQ progeny showed a significant increase in left ventricular (LV) dimension and reduction in cardiac function (measured by FS) compared with the DBA CSQ and the F1 CSQ mice (Figures 2B and 2C). These data indicate that the heart failure phenotype was accelerated in the N2 progeny of the backcross to B6 compared with the F1 CSQ mice. When the reciprocal backcross to DBA was performed (F1xDBA), we observed a lessening of the cardiac phenotype and an increase in the variance of phenotypic measures of cardiac function (Figures 2B and 2C). Although some N2 (F1xDBA) CSQ progeny showed marked enlargement of the LV and poor contraction compared with the F1 mice, other N2 (F1xDBA) progeny had relatively preserved cardiac function with similar values of LV end-diastolic diameter (LVEDD) as wild-type (Figures 2B and 2C). There were no significant differences in echocardiographic parameters between sexes in any of the genetic backgrounds (data not shown).
To determine whether the difference in the phenotypic expression of the N2 progeny might be due to differences in expression level of the CSQ transgene, we performed Western blotting analysis on heart extracts from the parental DBA CSQ, the F1, and the N2 progeny. Importantly, the level of CSQ protein expression was equal in all progeny tested, indicating that the variation in phenotype of the CSQ mouse cannot be accounted for by variation in the level of CSQ protein in the heart (Figure 2D).
To map the genes that modify the heart failure phenotype, progeny from the 2 reciprocal backcrosses were used in a genome-wide scan for linkage. One hundred microsatellite markers were chosen with an average spacing of 15 cM to detect quantitative trait loci (QTLs) linked to LVEDD and survival in the N2 progeny. For the N2 (F1xB6) backcross, the initial genome screen included 70 progeny to detect QTLs linked to survival. We identified one marker on chromosome 2 that showed linkage to survival with a LRS 31.7 (P<0.001) (Figure 3A). For interval mapping of this locus at 1-cM intervals, 39 more N2 progeny were genotyped with additional markers providing a denser map of the interval with a maximum LRS of 36.2 (equivalent to a LOD score of 7.8) (Figure 3B). We demonstrate the existence of Hfm1 (for heart failure modifier 1) on chromosome 2 and estimate that it accounts for approximately 29% of the variance in survival in this cross. Because 39% of the N2 (F1xB6) progeny died before their 4-week echocardiogram, we could not screen for QTLs affecting cardiac function in this backcross. To show a relationship between phenotype and genotype, survival was plotted against genotype at the chromosome 2 locus that showed the highest LRS (marker D2Mit327) (Figure 3C). As shown, a heterozygous genotype (DBA/B6) in the Hfm1 region is more often associated with increased survival than a homozygous (B6/B6) genotype. These data support our hypothesis that the DBA genetic strain is contributing a dominant protective allele to the phenotype.
|
Fifty-six progeny from the reciprocal backcross N2 (F1xDBA) were used to map modifier genes that alter survival and cardiac size. One marker on chromosome 3 showed linkage to survival with a LRS 23.0 (P<0.001) (Figure 4A) using single locus association tests. With the addition of 21 N2 progeny and a denser map of the interval, this linkage was confirmed with a LRS of 26.5 (equivalent to a LOD score of 5.7) (Figure 4B). We demonstrate the existence of Hfm2 on chromosome 3 and show that this modifier locus accounts for approximately 28% of the variance in survival in this backcross. There was also one marker mapping near Hfm2 that shows linkage to LVEDD with a peak LRS of 33.4 (P<0.001) (Figure 5A). Linkage to this region was also confirmed by interval mapping with a maximum LRS of 42.9 (equivalent to a LOD score of 9.3) and accounts for approximately 42% of the variance in LVEDD in this backcross (Figure 5B). Both traits (survival and LVEDD) appear to be linked to the same QTL (Hfm2) on chromosome 3, which is supported by the multipoint linkage analysis for both regions. We focused on the markerD3Mit86, which showed one of the highest LRS for both phenotypes, to determine whether there was an association between LVEDD and survival at that specific locus (Figure 5C). Importantly, a heterozygous genotype (DBA/B6) in the Hfm2 region is more often associated with a larger ventricle and shortened survival than a homozygous (DBA/DBA) genotype.
|
|
Because only one significant QTL was identified in each of the backcrosses, we used a haplotype analysis to further narrow the 2 candidate intervals. We identified 15 N2 progeny from the F1 backcross to DBA with both a recombinant haplotype for Hfm2 and a survival value greater than 0.5 standard deviation from the mean (Figure 6A). Ten mice were grouped into an early death category, and genotyping revealed that these animals were heterozygous (B6/DBA) in the region of the QTL. Five mice were grouped into a late death category and showed a homozygous DBA genotype for the QTL. The minimal candidate interval is defined by 11 independent crossovers (animal Nos. 25, 91, 12, 411, 501, and 590 for the proximal limit; animal Nos. 153, 64, 123, 72, and 316 for the distal limit), which narrow the Hfm2 candidate interval to a 1.1-cM region between markers D3Mit320 and D3Mit260.
|
The same approach was used to refine the Hfm1 interval; we identified 14 N2 progeny from our F1 backcross to B6 that have both a recombinant haplotype in the QTL region and a survival value greater than 0.5 standard deviation from the mean (Figure 6B). Six of these animals were grouped into an early death category, and genotyping revealed that these mice were homozygous for B6 within the region of the QTL. Eight animals were grouped into a late death category and showed a heterozygous (B6/DBA) genotype for the QTL. The minimal candidate interval is defined by 6 independent crossovers (animal Nos. 31, 237, 1025, and 470 for the proximal limit; animal Nos. 478 and 1350 for the distal limit), which narrow the Hfm1 interval to 5.2 cM between markers D2Mit325 and D2Mit418. Interestingly, 1 animal (No. 1181) appears discordant, with a survival just above the upper limit of mean±1/2SD exclusion.
| Discussion |
|---|
|
|
|---|
The phenotype-genotype relationship (Figures 3C and 5C) supports our hypothesis that the B6 genetic strain is contributing a dominant susceptibility allele to the cardiomyopathic phenotype in this model and is consistent with the clinical data showing the relationship between cardiac size and outcome (survival).15 However, our estimates of the percentage of genetic contribution to the phenotypes defined by Hfm1 and Hfm2 are overestimates, because we excluded animals clustered around the mean values for each. Regardless, these loci collectively do not account for 100% of the phenotypic variance in LVEDD and survival. We note that there is reproducible variation in the phenotype even in genetically identical F1 animals (Figure 2), suggesting that experimental error and stochastic events contribute to the variance. Furthermore, because the measured traits are continuous, we assume that the interaction of multiple modifiers with weaker effects may be at work but each may fall below the threshold for detection. Nonetheless, when we hold constant the Hfm1 and Hfm2 loci in the respective backcross mapping analyses, we find that only a single locus on chromosome 9 reaches the P=0.05 level, in only the N2 progeny from the DBA backcross.
Because only one major modifier locus was identified in each of the backcrosses, we were able to apply a haplotype analysis to map individual crossovers and further narrow the Hfm1 and Hfm2 intervals. The large number of independent crossovers that lie within the critical interval confirms that the modifiers are acting in a dominant manner and validates the experimental approach of using haplotype analysis. In the N2 (F1xB6) backcross, the haplotype of one animal (No. 1181) was found to be discordant. This is not surprising given the severity of the phenotype in this backcross and the clustering of survival around the mean. Indeed, animal No. 1181 died within 2 days of the upper threshold for exclusion and underscores the need for multiple independent crossover events to narrow the Hfm1 interval.
Within the Hfm1 region, a large number of genes map, including the gene titin, where it has recently been shown that mutations can cause autosomal-dominant dilated cardiomyopathy.16 Future studies will determine whether DNA sequence alterations in titin are responsible for the observed differences in the heart failure phenotype. Importantly, for both the Hfm1 and Hfm2 intervals, additional recombinant events within these regions will further refine the loci and lead to the identification of potential candidate genes.
Population studies have long suggested that common diseases (such as heart failure) have a strong genetic component. Although progress has been made in identifying disease loci causing inheritable forms of cardiomyopathy, such as that caused by mutations in cardiac actin,17 it is apparent that modifying genes define susceptibility to progressive cardiac dysfunction and sudden death in the heart failure population and are likely to be particularly important in mediating outcome. However, detecting which genes are relevant is difficult in human studies because environmental factors cannot be standardized and because human individuals differ in many genes in addition to the genes of interest. In this regard, the mouse provides a solution to both these difficulties and provides the opportunity to dissect the genetic factors involved in the disease process. Although not all mouse modifier loci will lead to the identification of corresponding loci in humans,18 we demonstrate here how the use of a sensitized mouse strain bred into different genetic backgrounds can significantly enhance the likelihood of successful mapping of modifier QTLs. In this regard, modifier genes have been successfully mapped in other mouse models of human disease, such as the secretory phospholipase gene Pla2g2a, identified as a modifier of adenomatous polyps,19,20 and a mouse gene that modifies cystic fibrosis.21,22 Whether modifier genes identified using this approach will be applicable to human heart failure will ultimately be determined through association studies in patient populations. Nonetheless, mapping of QTLs in mice may provide novel mechanistic insights into the pathogenesis and progression of heart failure.
| Acknowledgments |
|---|
| Footnotes |
|---|
This article originally appeared Online on March 25, 2002 (Circulation. 2002;105:r47r52).
Received January 28, 2002; revision received February 21, 2002; accepted February 21, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. J. Barrick, M. Rojas, R. Schoonhoven, S. S. Smyth, and D. W. Threadgill Cardiac response to pressure overload in 129S1/SvImJ and C57BL/6J mice: temporal- and background-dependent development of concentric left ventricular hypertrophy Am J Physiol Heart Circ Physiol, May 1, 2007; 292(5): H2119 - H2130. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Li, R. E. Waters, S. I. Redfern, M. Zhang, L. Mao, B. H. Annex, and Z. Yan Oxidative Phenotype Protects Myofibers from Pathological Insults Induced by Chronic Heart Failure in Mice Am. J. Pathol., February 1, 2007; 170(2): 599 - 608. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Donahue, D. A. Marchuk, and H. A. Rockman Redefining Heart Failure: The Utility of Genomics J. Am. Coll. Cardiol., October 3, 2006; 48(7): 1289 - 1298. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. P. Maddatu, S. M. Garvey, D. G. Schroeder, W. Zhang, S.-Y. Kim, A. I. Nicholson, C. J. Davis, and G. A. Cox Dilated cardiomyopathy in the nmd mouse: transgenic rescue and QTLs that improve cardiac function and survival Hum. Mol. Genet., November 1, 2005; 14(21): 3179 - 3189. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. I. Caron, L. R. James, H.-S. Kim, J. Knowles, R. Uhlir, L. Mao, J. R. Hagaman, W. Cascio, H. Rockman, and O. Smithies Cardiac hypertrophy and sudden death in mice with a genetically clamped renin transgene PNAS, March 2, 2004; 101(9): 3106 - 3111. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Le Corvoisier, H.-Y. Park, K. M. Carlson, D. A. Marchuk, and H. A. Rockman Multiple quantitative trait loci modify the heart failure phenotype in murine cardiomyopathy Hum. Mol. Genet., December 1, 2003; 12(23): 3097 - 3107. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. A. Doevendans Genetic Polymorphisms and Cardiac Failure Seminars in Cardiothoracic and Vascular Anesthesia, March 1, 2003; 7(1): 23 - 29. [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2002 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |