| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|
| ||||||||||||||||||||||||||||||||
a Department of Human Genetics, University of Michigan School of Medicine,
b Institute of Gerontology, University of Michigan, Ann Arbor
c Ann Arbor VA Medical Center, University of Michigan, Ann Arbor
d Department of Pathology and Geriatrics Center, University of Michigan, Ann Arbor
Richard A. Miller, The Geriatrics Center, University of Michigan, 1500 East Medical Center Drive, Ann Arbor, MI 48109-0940 E-mail: millerr{at}umich.edu.
Decision Editor: Edward Masoro, PhD
| Abstract |
|---|
|
|
|---|
VERY little is known about the number or position of loci with detectably large effects on life span in mammals. Heritability estimates (1) have indicated that genetic differences may account for as much as 2530% of the variance in life span in populations of flies, rodents, and humans, but they give no inkling as to what portion of this variation is attributable to loci with major effects. Studies of recombinant inbred mouse lines (2) and of backcrosses between heterogeneous lines of mice selected for differences in antibody production (3) have suggested that longevity may be influenced by a fairly small number of loci. Quantitative trait locus (QTL) mapping in Drosophila (4) has provided evidence for loci whose effects on longevity are limited to male or female flies. To see if the QTL approach could also help enumerate and map mammalian loci with substantial (i.e., detectable) influence on longevity, we have carried out a genomewide scan for loci associated with differential survival in a four-way cross among four inbred mouse lines commonly used for studies of aging and disease. At the same time we have looked for evidence that loci with effects on life span might be sex specific in their action, and we tested for possible epistatic interactions among the loci with strongest effects.
| Methods |
|---|
|
|
|---|
Genotyping
Genotyping was performed by standard polymerase chain reaction (PCR) amplification of genomic DNA from each animal, using marker loci obtained from the Mouse Simple Sequence Length Polymorphism Database, Whitehead/MIT Center for Genome Research (carbon.wi.mit.edu:8000/ftp/distribution/mouse_sslp_releases/may99) (6). Polyacrylamide gels were scored by silver staining or by using the ALFExpress automated sequencer as described (7). Analyses at 78 marker loci were performed on 253 individuals (110 males and 143 females) with an average intermarker interval of 23 cM. Of the markers, 69 are fully informative for all four grandparental alleles; the other 9 are informative for only one of the two parents. Genotypes were 95% complete on average for each marker. A full listing of loci as well as the complete genotype and phenotype data sets used in this analysis are available at sitemaker.med.umich.edu/dtburke/files/253mice_78markers.xls.
Statistical Analysis
The initial genomewide search for QTL was performed by using a single-point locus scan. An analysis of variance (ANOVA) was performed for single loci by using Proc MIXED in SAS version 7.0 (SAS Institute, Cary, NC). For fully informative loci, the effects specified were maternal allele, paternal allele, and maternal allele by paternal allele interaction. For partially informative loci, only the informative allelic pairs were examined. The single-locus analysis for the male only, female only, and combined data used the following statistical model (8):
![]() |
mi and
pj are the additive effects of the maternal alleles Mmi and the paternal alleles Mpi, respectively, of marker M;
ij is the maternal by paternal allelic interaction (dominance effect); and
k is the error term with N(0,
'2). For reasons described in the text, the linkage analyses were conducted twice: once on the entire data set and once on a data set (EDE, for early deaths excluded) that removed the 20% of the mice dying prior to 657 days of age. For linkage analysis, the data were transformed to satisfy the assumptions of the classical regression model, that is, normality and constant variance of residuals. For this to be accomplished, the data were logarithmically transformed for the EDE populations and for the male-only population, and they were raised to the power of 1.8 for the female-only and male-plus-female populations.
Following the single-locus genome scan, a two-locus search was performed by examining all of the possible pairwise allelic interactions in the data. The computational difficulty of the two-locus search necessitated a change in the analytical procedure. Rather than viewing the genome scan data set as 69 four-way informative loci (plus nine biallelic informative loci), the two-locus genome scan defined 147 informative biallelic genotypes per animal (77 maternally informative and 70 paternally informative). Each pairwise combination of biallelic genotypes was examined by using an ANOVA with the following statistical model:
![]() |
k is the error term with N(0,
'2). As with the single-locus analysis, the two-locus genome scan was performed on the male only, female only, and combined data sets, sequentially, and it was conducted on both the entire data set and then on the EDE data set.
Assessment of Statistical Significance
Permutation testing was the principal method used to estimate experimentwide statistical significance. Permutation tests were performed essentially as described by Churchill and Doerge (9). For each round of permutation, phenotypes were shuffled and distributed at random among the individual animals, with each animal retaining its complete genomewide genotype. The resulting synthetic data set disrupts the genotypephenotype association. The reshuffling was performed 1000 times, and the genome search was performed by using the original statistical model on each synthetic data set. The maximum test statistic for each synthetic genome scan is then recorded, and the threshold test statistic to obtain
= 0.05 is derived from the resulting null distribution.
In addition, for each of the single-locus hypotheses (male only, female only, and combined), a Bonferroni correction was used that multiplies the pointwise probability by 216 = 69 x 3 + 9. This correction reflects the use of 69 markers that were informative about maternal, paternal, and maternal by paternal interaction effects (i.e., N = 3 hypotheses for each marker), and 9 other markers that were informative about maternally inherited or paternally inherited alleles but not both. For the two-locus search, the Bonferroni correction multiplies the pointwise probability by 10,731 = 147 x 146/2, because the genotype information was derived from 147 informative biallelic genotypings in the genome scan.
A post hoc analysis of epistatic interactions used an ANOVA in which life span was modeled as a function of allele at each of the two interacting loci plus an interaction term; reported p values refer to the significance of the interaction term.
| Results |
|---|
|
|
|---|
The genomewide search for life-span-associated QTL was performed with two statistical analyses: once using all of the data, and then a second time using the EDE data set that excluded the first 20% of the mice to die (i.e., excluding all deaths prior to 657 days of age). We chose to conduct a separate examination of this truncated EDE data set in part because of earlier observations that the heritability of life span in mice increases with age at death (3). In addition, we considered that deaths at early ages are unlikely to reflect the effects of aging per se, and their inclusion in life tables thus can confound the search for genes that modify aging processes. In fact, necropsy data showed that many of the early deaths in males were due to a urinary syndrome (10)(11) that is typically seen only in group-housed males and is thought to reflect stresses associated with adjustments in dominance hierarchy rather than an effect of the aging process. Indeed, this urinary syndrome was responsible for 58% of the male deaths prior to 657 days of age, but only 5% of the male deaths after this age.
The main analysis method was based on an ANOVA (8) to analyze single loci. The statistical strategy was also designed to detect loci whose effect on life span might be limited to males or to females, because loci with sex-specific effects on longevity have been documented in QTL studies in Drosophila (4). For these reasons we tested three hypotheses in each of two data sets (untruncated and EDE), first looking for loci that influence males alone, then for loci in females alone, and, last, loci in data pooled across genders. The combined data set provides improved power to detect QTL that are not sex specific. The combined data set was treated as independent for the purposes of Bonferroni correction, even though it pools data from the two single-sex data sets. Significance levels were based on permutation tests (9), using the same statistical model and experimental data to develop experimentwise confidence thresholds that adjust for the simultaneous evaluation of multiple nonindependent hypotheses.
No loci were found to be predictors of longevity in the entire population (using an experimentwise significance criterion of p < .05 to reduce the reporting of false positive associations; (12)), perhaps because the genetic factors that influence mortality risk at early stages of the life span differ from those that modulate aging and late-life disease.
In the EDE data set, however, three loci were found to be predictors of life-span differences with comparisonwise (i.e., unadjusted) probabilities ranging between .0004 and .00002. These are summarized in Table 1 . These pointwise probability estimates require adjustment for the simultaneous testing of multiple hypotheses, that is, simultaneous assessment of the 78 informative genetic polymorphisms used in the analysis. Our principal method of adjustment was to calculate an experimentwise probability by permutation analysis as suggested by Churchill and Doerge (7), and these values are included in Table 1 . By these criteria, both D12Mit167 and D9Mit110 were significant at the p < .01 level in the EDE data set, and D10Mit15 met the p < .05 criterion. We also calculated the Bonferroni-corrected significance threshold for multiple independent tests and noted that two of the three loci met this conservative test (p < .05).
|
|
|
Interaction Between Maternal and Paternal Genotypes
Fig. 3 presents the survival curves, among EDE mice, for the four genotypes discriminated by alleles segregating at D12Mit167. For the longest lived 80% of the population, D12Mit167 genotypes differ in life span (experimentwise p = .01), with a significant interaction between the maternal and paternal alleles (post hoc p = .0001 by two-factor ANOVA). Both males and females are equally affected (not shown). For D12Mit167, the mean longevity values (± SD) in the EDE data set were as follows: genotype B6/C3, 886 ± 111; genotype C/D2, 867 ± 110; genotype B6/D2, 828 ± 116; and genotype C/C3, 790 ± 105. It is noteworthy that superior longevity is not associated with any one allele at this locus, but instead with a particular combination of alleles: the combination B6 + C3 is associated with longer survival, though neither the B6 nor the C3 allele confers increased longevity in combination with the alternate counteralleles. Further work will be needed to determine whether this pattern represents interaction among alleles of a single effector locus, or interaction among alleles of linked but distinct loci on chromosome 12.
|
|
|
| Discussion |
|---|
|
|
|---|
We found no genes that had significant effects on longevity in males, or in both sexes combined, until we reanalyzed our data set by excluding the first 20% of mice to die. We considered a number of alternate strategies, including excluding larger proportions of the entire population, but we found that the statistical significance of our candidate loci diminished when the number of mice included in the study group became too small to provide adequate statistical power. The decision to exclude the first 20% is thus to some extent arbitrary, based in part on empirical considerations and in part on a desire to separate early-life deaths from those deaths, at later ages, which we think are more likely to be timed by aging. When in follow-up studies we accumulate sufficient numbers of mice, we will be able to test more thoroughly the idea that these (or other) QTL influence mortality risks in the oldest-surviving groups of mice.
An earlier report (14) summarized preliminary evidence for QTL associations with longevity in a subset (n = 129) of the same mouse population presented in this paper, using a slightly different set of genotypic markers. The earlier report listed seven candidate QTL with pointwise p < .01, but none of these met experimentwise significance criteria. Two of the seven candidate QTL are now supported, with stronger evidence, in the present report. The C allele of D16Mit182 was found, in the initial report, to be associated with increased survival in female mice, and it is now shown to reach experimentwise p = .06 when considered with D2Mit58. The D2 allele of D10Mit40 was, in the initial report, found to be associated with increased life span in male mice; D10Mit40, at position 21 cM, is closely linked to D10Mit10 (position 25 cM), which is in the present report shown to be associated with increased life span in males with experimentwise p < .05. (The initial paper excluded male mice dying with MUS; the current paper, instead, examined mice of either sex dying after 657 days.) The current report cannot be taken as an independent replication of the earlier report, because the mice examined are a superset of the initial population. It is noteworthy that six of the seven initially reported QTL were also apparently sex specific in their effects. An analysis of larger numbers of mice will be needed to determine whether any of the other markers reported in the initial publication will eventually prove to be associated with life span to a significant degree.
Our new data support the surprising conclusion that QTL may influence life span through sex-specific pathways. Of the three individual loci with significant effects on life span (Table 1 ), two influence life span in males but not in females, and the pair of loci illustrated in Fig. 4 influence survival only in females. Only one locus, D12Mit167, has equivalent effects in mice of both sexes. Most models of genetic influences on life span propose mechanisms in which these hypothetical polymorphisms act on pathways such as oxidant resistance, tumor suppression, mitochondrial function, or other processes thought likely to influence disease resistance equally in males and females. It is noteworthy that QTL analyses in Drosophila have also produced evidence for alleles with effects on life span limited to males or females only (4). Genome scans have also been successful in revealing QTL that affect life span in Caenorhabditis elegans (15)(16)(17)(18).
Our data also led to a second interesting finding: alleles whose effects on life span are conditional on the alleles inherited at another locus (see Fig. 5). Analyses of additional mice will be needed to address the question of whether the interaction between maternal and paternal alleles at D12Mit167 (Fig. 3) reflects epistatic effects between two distinct but closely linked loci or an interaction among alleles at a single locus. The paired locus scan that produced evidence for the additive effect illustrated in Fig. 4 also produced suggestive evidence (not shown) for epistatic interactions among other pairs of loci, but the numbers of mice tested are at this stage too small to allow detection of most paired-locus effects at high statistical power. It is striking that three (and possibly four, if one includes D12Mit167) of the five loci detected in the single and paired-locus genome scans have effects on life span that are conditional on inheritance at other loci. We therefore suspect that the significant interactions documented in Fig. 5 will eventually prove to be a minimal estimate of the importance of interlocus effects on life span. A recent study of QTL with effects on life span in C. elegans (17) also found evidence for two cases of epistatic interactions between pairs of loci, with suggestive evidence for several others, and epistatic interactions among life span QTL have also been noted in studies of Drosophila melanogaster (19).
It will be highly informative to determine whether the genetic differences that influence life span among these sibling mice can also influence the pattern of progression of age-sensitive traits, such as changes in immune function, muscle atrophy, protein glycation, and cataract progression, and to see whether the loci influence either the frequency or the rate of progression of specific neoplastic and nonneoplastic diseases. Higher resolution mapping of the loci that influence life span in this systemwhich will require the analysis of larger numbers of animalsmay be able to narrow the interval of interest to the point where it becomes feasible to identify the genes that influence life span by a candidate gene approach. It should be possible to use the QTL mapping data to identify, at birth, cohorts of mice that are very likely to live longer and age more slowly than their siblings. These genetic predispositions may then be exploited to test specific, mechanistic hypotheses about the connections that link genotype, aging, disease, and life span.
| Acknowledgments |
|---|
We thank Luann Linsalata, Gretchen Buehner, Emily Gray, and Jen Chisa for technical assistance. The necropsy diagnosis of mouse urinary syndrome was conducted by the late Dr. Clarence Chrisp.
Received April 25, 2001
Accepted August 2, 2001
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Flurkey, Y. Brandvain, S. Klebanov, S. N. Austad, R. A. Miller, R. Yuan, and D. E. Harrison PohnB6F1: A Cross of Wild and Domestic Mice That Is a New Model of Extended Female Reproductive Life Span J. Gerontol. A Biol. Sci. Med. Sci., November 1, 2007; 62(11): 1187 - 1198. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Kriete Biomarkers of Aging: Combinatorial or Systems Model? Sci. Aging Knowl. Environ., January 4, 2006; 2006(1): pe1 - pe1. [Abstract] [Full Text] |
||||
![]() |
J. K. Quarrie and K. T. Riabowol Murine Models of Life Span Extension Sci. Aging Knowl. Environ., August 4, 2004; 2004(31): re5 - re5. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Wolf, A. Galecki, R. Lipman, S. Chen, M. Smith-Wheelock, D. Burke, and R. Miller Quantitative Trait Locus Mapping for Age-Related Cataract Severity and Synechia Prevalence Using Four-Way Cross Mice Invest. Ophthalmol. Vis. Sci., June 1, 2004; 45(6): 1922 - 1929. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. N. Butler, S. N. Austad, N. Barzilai, A. Braun, S. Helfand, P. L. Larsen, A. M. McCormick, T. T. Perls, A. R. Shuldiner, R. L. Sprott, et al. Longevity Genes: From Primitive Organisms to Humans J. Gerontol. A Biol. Sci. Med. Sci., July 1, 2003; 58(7): B581 - 584. [Full Text] [PDF] |
||||
![]() |
J. W. Curtsinger Sex Specificity, Life-Span QTLs, and Statistical Power J. Gerontol. A Biol. Sci. Med. Sci., December 1, 2002; 57(12): B409 - 414. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Miller, C. Chrisp, A. U. Jackson, A. T. Galecki, and D. T. Burke Coordinated Genetic Control of Neoplastic and Nonneoplastic Diseases in Mice J. Gerontol. A Biol. Sci. Med. Sci., January 1, 2002; 57(1): B3 - 8. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||
| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|