Genotype-Environment Interaction for Quantitative Trait Loci Affecting Life Span in Drosophila melanogaster

Corresponding author: Trudy F. C. Mackay, Department of Genetics, Box 7614, North Carolina State University, Raleigh, NC 27695., trudy_mackay{at}ncsu.edu (E-mail)

Communicating editor: P. D. KEIGHTLEY

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

The nature of genetic variation for Drosophila longevity in a population of recombinant inbred lines was investigated by estimating quantitative genetic parameters and mapping quantitative trait loci (QTL) for adult life span in five environments: standard culture conditions, high and low temperature, and heat-shock and starvation stress. There was highly significant genetic variation for life span within each sex and environment. In the analysis of variance of life span pooled over sexes and environments, however, the significant genetic variation appeared in the genotype x sex and genotype x environment interaction terms. The genetic correlation of longevity across the sexes and environments was not significantly different from zero in these lines. We estimated map positions and effects of QTL affecting life span by linkage to highly polymorphic roo transposable element markers, using a multiple-trait composite interval mapping procedure. A minimum of 17 QTL were detected; all were sex and/or environment-specific. Ten of the QTL had sexually antagonistic or antagonistic pleiotropic effects in different environments. These data provide support for the pleiotropy theory of senescence and the hypothesis that variation for longevity might be maintained by opposing selection pressures in males and females and variable environments. Further work is necessary to assess the generality of these results, using different strains, to determine heterozygous effects and to map the life span QTL to the level of genetic loci.

LIMITED life span and senescence, the progressive decline in survivorship and fertility with advancing age, are near-universal characteristics of eukaryotic organisms. However, there is great variation in total life span and rate of aging between and within species (FINCH 1990 ). Variation in life span within natural populations is partly attributable to both genetic and environmental effects (TOWER 1996 ; MCCLEARN et al. 1997 ). Determining the nature of, and interactions between, the genetic and environmental factors affecting aging is important from two rather different perspectives: human health and evolutionary theory. As the human population rapidly grows older (MARTIN and PRESTON 1994 ), population aging is becoming one of the most important social and health problems facing the next half century (ROWE 1997 ). The evolutionary enigma is why organisms age at all, and why genetic variation for life span segregates within populations.

The classical evolutionary explanations for aging are based on the reduction of the strength of natural selection as organisms grow older (MEDAWAR 1952 ; WILLIAMS 1957 ; CHARLESWORTH 1980 ). The mutation-accumulation hypothesis (MEDAWAR 1952 ) assumes that mutations that are neutral early in life but that have variable late-age-specific deleterious effects on viability will escape the vigilance of natural selection and cause both aging and variation in aging. The antagonistic pleiotropy hypothesis proposes that mutations with beneficial effects on survival and/or fertility early in life have deleterious effects later on, and vice versa, leading to an equilibrium situation in which both classes of alleles remain segregating (WILLIAMS 1957 ; ROSE and CHARLESWORTH 1980 ; KIRKWOOD and ROSE 1991 ). Additionally, hypotheses for the maintenance of genetic variation for quantitative traits in general apply to longevity. These include unconditionally deleterious mutation-selection balance (HOULE et al. 1996 ) and genotype x environment interaction (GEI). GEI can maintain quantitative genetic variation if heterozygotes are less sensitive to environmental variation than homozygotes (GILLESPIE and TURELLI 1989 ) and if there are reversals in effects of alleles affecting the life history trait in different environments (LEVENE 1953 ). None of these hypotheses are mutually exclusive.

Attempts to discriminate among the various evolutionary hypotheses using Drosophila melanogaster as a model system, by testing whether genetic variation in fitness components increases with age, as predicted by the mutation-accumulation hypothesis (ROSE and CHARLESWORTH 1980 ; ENGSTROM et al. 1989 ; HUGHES and CHARLESWORTH 1994 ; HUGHES 1995 ; PROMISLOW et al. 1996 ; TATAR et al. 1996 ), and testing for negative genetic correlations between longevity and other life history traits, as predicted by the antagonistic pleiotropy hypothesis (ROSE and CHARLESWORTH 1981A , ROSE and CHARLESWORTH 1981B ; ROSE 1984 ; LUCKINBILL et al. 1984 ; ENGSTROM et al. 1992 ; PARTRIDGE and FOWLER 1992 ; ROPER et al. 1993 ; SERVICE 1993 ; STEARNS and KAISER 1993 ; FORCE et al. 1995 ; ZWAAN et al. 1995 ), have produced mixed results. Given that the various studies did not have the same power to detect differences in variance, the potential for real genetic heterogeneity among the base populations used, and different population sizes and environments, such inconsistency is perhaps not surprising. Further, estimates of genetic correlation are only informative with respect to the net effects of all segregating loci and can conceal considerable heterogeneity in effects across loci. The evolutionary question is not whether each of the potential mechanisms operates, but rather at what quantitative trait loci (QTL) do alleles affecting variation in life span segregate, and at what fraction of the QTL is variation attributable to the segregation of rare, unconditionally deleterious alleles, alleles with late-age effects, alleles with opposing pleiotropic effects on early fitness components, and alleles that exhibit GEI.

What genes are likely to affect senescence and life span? Traditional screens for mutations altering longevity have not been considered feasible for a variety of reasons:

1. Mutations in almost any gene will have a deleterious pleiotropic effect on life span, but this does not imply the genes so identified directly affect life span.

2. Screens for mutations increasing life span are logistically difficult, because each potential mutation must be preserved prior to scoring the phenotype (date of death).

3. Chemical mutagenesis induces multiple mutations simultaneously; therefore, long-lived mutations will be difficult to discern given the background noise of mutations decreasing life span.

4. Single P-element mutagenesis in Drosophila is not a solution to the problem of deleterious background mutations, since the transgenic P element carries a dominant wild-type selectable marker gene and the host strain is homozygous for a mutant allele for the marker. The host strain thus has reduced fitness (and longevity) that can be "rescued" by many P-element insertions, which do not affect loci with direct effects on life span (LYMAN et al. 1996 ; but see LIN et al. 1998 ).

In Caenorhabditis elegans, however, mutations in genes affecting dauer-larva formation, fertility, and biological rhythm have been found to have extended longevity phenotypes (LITHGOW 1996 ). Several of these genes have been cloned and encode phosphatidyl inositol 3-kinase (age-1/daf-23; MORRIS et al. 1996 ), an insulin-like receptor (daf-2; KIMURA et al. 1997 ), a member of the HNK-3/forkhead transcription factor family (daf-16; LIN et al. 1997 ), and a short protein with an 82-amino-acid tandem repeat that is homologous to yeast CAT5/COQ7 (clk-1; HEKIMI et al. 1998 ). These genes strongly implicate the glucose metabolism pathway in organismal aging. These observations are consistent with observations of increased glycogen and lipid content of Drosophila lines selected for postponed senescence relative to the unselected controls (SERVICE 1987 ; GRAVES et al. 1992 ); an increase in frequency of the most active electromorph of the glycolytic enzyme phosphoglucomutase in the same long-lived Drosophila strains compared to the base population controls (DECKERT-CRUZ et al. 1997 ); and the effect of severe caloric restriction on extended life span in mice (HARRISON and ARCHER 1987 ).

Long-lived strains of flies and worms also exhibit increased resistance to starvation, desiccation (SERVICE et al. 1988 ; GRAVES et al. 1992 ; LIN et al. 1998 ), and thermal and UV-irradiation stress (LITHGOW 1996 ; LIN et al. 1998 ), implicating oxidative stress and heat-shock response genes as candidates. Numerous observations support this hypothesis. The more active superoxide dismutase (Sod) electromorph allele is at a higher frequency in the long- than in the short-lived replicates of one set of Drosophila selection lines (DECKERT-CRUZ et al. 1997 ), while in an independent set of lines, Sod, catalase (Cat), and xanthine dehydrogenase (Xdh) have higher mRNA and enzyme activity levels in the long-lived strains (FORCE et al. 1995 ). Long-lived mutant strains of C. elegans have elevated Sod and Cat enzyme activities (LITHGOW 1996 ). Transgenic Drosophila strains that jointly overexpress both Sod and Cat have extended life spans (ORR and SOHAL 1994 ). Heat-shock proteins (Hsps) are produced in response to thermal and other stress factors (BURTON et al. 1988 ). Not only is the ability to induce Hsps reduced in senescent animals (NIEDZWIECKI et al. 1991 ), but a brief, sublethal heat shock to young animals has been reported to extend life span (LITHGOW et al. 1995 ; KHAZAELI et al. 1997 ). Additional candidate genes affecting longevity are those associated with DNA repair and replication (WOODHEAD et al. 1985 ; HARRISON 1990 ; YU et al. 1996 ) and telomere length (BODNAR et al. 1998 ).

Analysis of effects of mutations at candidate genes and their expression patterns in young and old animals does not, however, directly address the question of whether allelic variation at these loci in natural populations causes quantitative variation in longevity. For this, allelic association studies are necessary, in which polymorphism(s) at the candidate gene(s) are identified, and the association between the allelic state of the polymorphism and longevity in a natural population is evaluated statistically (FINCH and TANZI 1997 ; see also MACKAY and LANGLEY 1990 ; LAI et al. 1994 ; LONG et al. 1998 ). This approach is potentially exceedingly powerful but is limited by the sheer number of putative candidate genes to be tested, our restricted knowledge of all relevant candidates, and the difficulty in assessing effects of the candidate gene alleles in a range of environments if individuals are sampled from nature.

QTL analysis of naturally occurring modifiers of life span directly addresses the issue of what genetic regions are associated with variation in longevity between the two parental strains used to generate the mapping population and what are their effects. This approach complements association studies at candidate loci, because it can refine the list of candidate genes to be investigated to those that map to the same location as the QTL and can point to genomic regions containing no known candidate genes as worthy of further study. Previously, we mapped five autosomal QTL affecting D. melanogaster life span that segregated between two inbred strains that had not been selected for life span, using a panel of 98 recombinant inbred (RI) lines derived from the parental strains and a dense cytogenetic marker map (NUZHDIN et al. 1997 ). The effects of longevity QTL, estimated in a single environment, were highly sex-specific. Here, we have reared the same set of RI lines under standard and stressful environmental conditions to determine whether longevity QTL exhibit GEI. The stressful conditions—high and low temperature, a brief, sublethal heat shock early in adult life, and starvation—were chosen on the basis of reports in the literature of direct and correlated effects of these treatments on life span. Surprisingly, all significant genetic variation among these RI lines is associated with GEI: longevity QTL exhibit sex- and environment-specific effects, as well as reversal in the signs of effects across sexes and environments.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Drosophila stocks:
The parental lines used were isogenic derivatives of two unrelated strains, Oregon-R (LINDSLEY and ZIMM 1992 ) and 2b (PASYUKOVA and NUZHDIN 1993 ). These strains were not selected for life span. F₁ progeny were backcrossed to 2b, and the backcross progeny randomly mated for four generations. A total of 98 RI lines were constructed by 25 generations of full-sib mating (NUZHDIN et al. 1997 ; GURGANUS et al. 1998 ).

Life span phenotypes:
Adult longevity was measured for all 98 RI and the two parental lines in five environments. All flies were reared from egg to adult on 10 ml standard cornmeal-agar-molasses medium at 25° in shell vials. The density of the stocks was controlled for three generations prior to the start of the longevity assays by restricting egg laying to 3 days and initiating the cultures with 10 pairs of flies. For each assay performed, 20 virgin males and females per line were collected in a 24-hr period. Assays were begun with 2-day-old flies housed in replicate vials with 10 same-sex individuals per vial. The five environments were as follows: control (C), 5 ml cornmeal-agar-molasses medium at 25°; high temperature (HT), 5 ml cornmeal-agar-molasses medium at 29°; low temperature (LT), 5 ml cornmeal-agar-molasses medium at 14°; heat shock (HS), 37° heat shock for 30 min followed by maintenance in the C environment; and starvation (S), 5 ml of 1.5% agar medium at 25°. To provide nutrition and humidity during the 30-min HS treatment, flies were held in vials containing a layer of filter paper soaked in a 25% sucrose solution. Numbers of dead flies were scored daily for the first four treatments, until all animals were dead, and vials were replaced weekly. The numbers of dead flies were scored every 8 hr for the starvation treatment. The five assays were conducted sequentially over a 2-mo period.

Phenotypic analyses:
Variation in male and female life span in each of the five environments was partitioned into sources attributable to line (L), replicate (R) within line, and error by random effects analysis of variance (ANOVA), according to the model y = µ + L + R(L) + error. µ is the overall mean and nested effects are in parentheses. Variation in life span for the full data set was partitioned into sources attributable to environment (E), sex (S), and L according to the mixed-model ANOVA: y = µ + E + S + L + E x S + E x L + S x L + E x S x L + R(E x S x L) + Error. E, S, and E x S are fixed effects; the rest are random. Tests of significance of F-ratios and estimates of variance components of the random effects were obtained using SAS procedures GLM and VARCOMP, respectively (SAS INSTITUTE 1988). Type III sums of squares were used in these analyses, because the design was not completely balanced at the lowest level. Genetic correlations between the sexes (r_GS) and across pairs of environments for each sex (r_GE) were computed as cov₁₂/_L1L2, where cov₁₂ is the covariance among line means between males and females or pairs of environments, and _L1 and _L2 are the square roots of the respective variance components among lines from the reduced-model analyses by sex and environment (ROBERTSON 1959 ). All phenotypic analyses were performed on untransformed life span records as well as on ln-transformed data.

Molecular marker map:
Highly polymorphic cytological insertion sites of high copy number roo transposable elements were used as molecular markers. roo element insertion sites in the RI lines were determined as described by NUZHDIN et al. 1997 . There were a total of 76 informative markers, including the recessive mutation spa^pol, which also segregated among the RI lines. The average spacing between the markers on the standard Drosophila map (LINDSLEY and ZIMM 1992 ) was 3.2 cM (NUZHDIN et al. 1997 ), but the additional generations of recombination during the construction of the RI lines resulted in an expanded genetic map. The map positions of the markers (d cM) on the expanded map were estimated from the observed recombination frequencies (r) using the Kosambi mapping function: 100d = 0.25 ln[(1 + 2r)/(1 - 2r)]. The insertion sites of the informative markers and their estimated map positions were as follows: 1B(1-0.00), 3E (1-14.85), 4F (1-33.16), 5D (1-39.36), 6E (1-42.43), 7D (1-57.28), 7E (1-58.30), 9A (1-76.81), 10D (1-96.29), 11C (1-99.36), 11D (1-102.43), 12E (1-106.52), 14C (1-126.01), 15A (1-127.03), 17C (1-131.16), 19A (1-134.27), 21E (2-0.00), 22F (2-8.59), 27B (2-61.19), 29F (2-85.53), 30AB (2-87.65), 30D (2-89.77), 33E (2-126.77), 34EF (2-135.75), 35B (2-140.15), 38A (2-146.00), 38E (2-148.49), 43A (2-152.16), 43E (2-153.38), 46C (2-154.52), 48D (2-163.49), 49D (2-169.92), 50B (2-172.04), 50D (2-175.15), 50F (2-178.47), 57C (3-0.00), 57F (3-7.65), 60E (3-39.23), 61A (4-0.00), 63A (4-16.72), 65A (4-43.80), 65D (4-67.70), 67D (4-83.06), 68B (4-103.00), 68C (4-114.92), 69D (4-124.71), 70C (4-142.04), 71E (4-146.59), 72A (4-148.85), 73D (4-153.26), 76A (4-156.76), 76B (4-160.18), 77A (4-162.34), 82D (4-164.65), 85F (4-169.26), 87B (4-177.28), 87E (4-180.67), 87F (4-181.81), 88E (4-192.05), 89B (4-196.56), 91A (4-207.31), 91D (4-215.33), 92A (4-218.65), 93A (4-230.09), 93B (4-231.15), 94D (4-243.34), 96A (4-260.05), 96F (4-273.56), 97D (4-284.86), 97E (4-287.06), 98A (4-295.83), 99A (4-313.98), 99B (4-326.17), 99E (4-337.22), 100A (4-338.28), spa (5-0.00). Cytological divisions 1-20, 21-60, and 61-100 correspond to Drosophila chromosomes 1, 2, and 3, respectively, on the standard map (LINDSLEY and ZIMM 1992 ). Note that chromosome 2 was formally split into two linkage groups due to free recombination between markers 50F and 57C. All markers were scored as homozygous Oregon, homozygous 2b, or heterozygous (NUZHDIN et al. 1997 ).

QTL mapping:
Multiple-trait composite interval mapping (JIANG and ZENG 1995 ) was used to test the hypotheses that an interval flanked by two adjacent markers contains QTL affecting the trait, and that it accounts for a QTL x trait interaction, while statistically accounting for the effects of additional segregating QTL by multiple regression on markers outside the tested interval. The analysis of one trait in multiple environments and the cross-environment genetic correlation is formally equivalent to the analysis of multiple genetically correlated traits in a single environment (FALCONER 1952 ). The likelihood-ratio test statistic, LR, is -2 ln (L₀/L₁), where L₀ is the maximum likelihood under the null hypothesis that a₁ = a₂ = ... = a_i = 0, where a_i is the additive effect of the QTL in environment i; and L₁ is the maximum likelihood under the alternative hypothesis that at least one a_i 0. In other words, this is a test for the presence of QTL in one or more environments. The test statistic at a genomic location is asymptotically distributed as ² with i + 1 d.f. (The additional degree of freedom is because one tests for a location parameter in addition to the i QTL effects.) The likelihood-ratio test for QTL x environment interaction evaluates the ratio of the maximum likelihood under the null hypothesis that a₁ = a₂ = ... = a_i, to that under the alternative hypothesis of a₁ a₂ ... a_i. If this test is performed only in the regions where QTL were detected in the joint mapping analysis, the test statistic is asympotically distributed as ² with i - 1 d.f. under the null hypothesis. If a genome-wide scan for interaction effects is performed, the asymptotic distribution is ² with i d.f.

The appropriate threshold for significance of each test must be adjusted for the number of independent tests in the genome-wide scan. Although there are 76 markers, there are not 76 independent tests, because the markers on the same chromosome are correlated. Under composite interval mapping (ZENG 1994 ), the effective number of independent tests in the search for QTL is determined by the size of the genome region to either side of the test interval in which marker cofactors are not fitted (the "window size"). In these analyses, we used a window size of 30 cM and thus conservatively estimate that adjacent regions of 50 cM may be regarded as independent. We estimate the number of independent tests as _i [(T_i/50) + 1], where T_i is the total estimated map length in centimorgans of the ith linkage group. For these data, there are 3.7, 4.6, 1.8, and 7.8 independent tests for each of the four linkage groups, respectively; or ~20 independent tests in total. Consequently, a type I error rate of = 0.0025 was used in the joint mapping analyses and genome-wide screens for interaction. For tests of QTL x environment interaction that were conditional on the presence of significant QTL, a conventional 5% significance level is appropriate. To be conservative, we used a type I error rate of = 0.01 in conditional tests for interaction. Seven joint composite interval mapping analyses were performed: one for each of the five environments, estimating male and female QTL effects separately, joint QTL, and QTL x sex effects; and one for each sex, estimating joint QTL and QTL x environment effects. All analyses used a window size of 30 cM and 15 background markers, selected by stepwise multiple regression. Heterozygous genotypes were treated as missing data. Analyses were performed on life span and on ln(life span).

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Life span phenotypes:
Mean life spans of the two parental lines, Oregon and 2b, and of the 98 RI lines derived from them, are given for each of the five environments in Table 1. The results of tests of significance of differences in life span between environments and sexes are presented in Table 2. Mean life spans of the RI lines ranged from a low of 2.3 days in the starvation treatment to a high of 80 days at low temperature. Relative to the control environment, the heat-shock treatment prolonged average life span by 6 days in males only, and the high temperature environment reduced average life span by 3.6 days in males and females. These treatment effects are consistent with previous reports documenting an inverse relationship between life span and temperature (e.g., ZWAAN et al. 1995 ) and demonstrating protective effects of early, sublethal heat shock on Drosophila longevity (LITHGOW et al. 1995 ; KHAZAELI et al. 1997 ). Mean differences in life span between all pairs of treatments were highly significant as was the overall effect of environment (Table 2). Females had significantly longer mean life spans than males in the control, high temperature, and starvation environments, but not under heat-shock and high temperature conditions; over all environments, sexual dimorphism for life span was only significant at P < 0.05 (Table 2). Differences in the magnitude of sexual dimorphism for life span across environments are reflected in highly significant (P < 0.0001) sex x environment interaction terms between the heat-shock and control, heat-shock and high temperature, and heat-shock and starvation environments, as well as between the high temperature and starvation environments (data not shown). Consequently, the sex x environment interaction term approached formal significance in the analysis over all five environments (Table 2).

Genetic variance in life span:
The among-line and residual variance components estimated from ANOVA for each sex in each of the five environments are given in Table 1. In all but one case, there is highly significant variation among the RI lines. (The exception is for female longevity at low temperature, for which the among-line variance component is not quite formally significant.) All among-line variance components were significant for the ln-transformed data (not shown). If one assumes that gene frequencies at all segregating loci affecting the trait are 0.5, then the variance among the RI lines is an estimate of the genetic variance, V_G, between the parental inbred lines (FALCONER and MACKAY 1996 ). Because the mean life span () of the RI lines varies across treatments, the coefficient of genetic variation CV_G = 100 is an appropriate statistic to compare relative magnitudes of genetic variance in the different environments (HOULE 1992 ). Estimates of CV_G ranged from 14–28% in males and 13–20% in females.

The ratio of V_G to (V_G + V_R), where V_R is the sum of the variance between and within replicates, indicates the extent to which variation in phenotypes in the population of RI lines is due to variation in genotypes. This ratio is 0.17, averaged over males and females in all but the starvation environment. The relatively large effect of uncontrolled environmental factors on individual phenotype is expected for life history traits such as life span. However, V_G/(V_G + V_R) was much higher for the starvation treatment: ~0.40 averaged over the sexes. Ratios of V_G to V_G + V_R were reduced on average by 70% in the analyses of ln(life span) (not shown).

Genotype x environment interaction for life span:
We assessed the extent to which alleles at loci affecting variation in life span among the RI lines cumulatively exhibited genotype x sex (GSI) and GEI interaction by determining the significance of these interaction effects in mixed model ANOVAs of life span among the RI lines within each environment, pooled across sexes, as well as pooled across environments. Highly significant GSI was observed for the control, heat-shock, and high temperature treatments; GSI was nominally significant in the starvation environment but not the low temperature environment (Table 1). The genotype x sex interaction effect was also significant in the analysis considering all five environments (Table 2). There was highly significant GEI for life span among the five environments considered here (Table 2). Although the three-way genotype x sex x environment term was not significantly different from zero in the overall analysis, this term was highly significant in the pairwise control-starvation, heat shock-starvation and high temperature-starvation analysis (data not shown). Thus, the genetic variation in sexual dimorphism for life span is itself sensitive to environmental conditions. The significant genotype x sex and genotype x environment interaction variances observed on the natural scale are also significant in the analyses on ln(life span) and are thus not an artifact of scale (data not shown).

Measures of the importance of GSI and GEI are r_GS and r_GE: the cross-sex and cross-environment genetic correlations, respectively. These correlations indicate the extent to which the same genes affect male and female life span within each environment, or the life span of males or females in different environments. Correlations approaching 0 indicate different constellations of genes affect the trait in the sexes or environments, and correlations approaching |1| indicate the same genes are responsible for variation in the trait (FALCONER 1952 ) in males and females or treatments. The correlation in mean life span of RI lines with the same genotype across both sexes and all environments is r_GSE = V_G/(V_G + V_GSI + V_GEI), where V_G is estimated by the among-line variance component and V_GSI and V_GEI are, respectively, the GSI and GEI variance components from the mixed-model ANOVA pooled over environments and sexes (BULMER 1985 ). For the untransformed data, r_GSE = 0.10 (Table 2) and for the ln-transformed data, r_GSE = 0.13 (data not shown). All of the significant genetic variation in life span among these RI lines is tied up in the interaction terms!

We computed cross-sex genetic correlations and cross-environment genetic correlations for pairs of environments from variance components using the method of ROBERTSON 1959 . Approximate confidence limits for these correlations were estimated using Fisher's z transformation (SOKAL and ROHLF 1981 ), which can be used to infer whether a given correlation is significantly different from zero or unity. The ANOVA results can also be used in this regard: cross-sex or cross-environment genetic correlations are significantly different from 1 if the GSI or GEI interaction terms from the pooled analysis are significant, and are significantly different from 0 if the among-line variance from the same analysis is significant. Estimates of r_GS within each environment are presented in Table 1. The cross-sex genetic correlations in the control and heat-shock environments are not significantly different from 0 and are highly significantly different from 1, indicating that different QTL affect the life span of the two sexes in these environments. In the low temperature environment, the cross-sex genetic correlation is not significantly different from unity, suggesting that the same genes affect variation in life span in this environment. Intermediate situations are represented by the high temperature and starvation treatments, where the genetic correlations of life span between the sexes are significantly different from both 0 and 1. Cross-sex genetic correlations can depart from unity if there are changes in the rank order of line mean longevities in males and females (crossing of reaction norms) and if there are changes in the among-line variance of life span in the two sexes. In all cases for which the sex x line interaction is significantly different from 0, >90% of the interaction variance is attributable to changes in rank order of life span between males and females. Consistent with this analysis, estimates of cross-sex genetic correlations of ln-transformed data are very similar to the estimates obtained on the natural scale (not shown).

Estimates of r_GE between pairs of environments and confidence limits of the estimates are given in Table 3. Cross-environment genetic correlations fall into four categories. The first group are high and positive; these are highly significantly different from zero and either not significantly different from unity or only nominally so. In this group are the control-high temperature and control-heat-shock environment pairs. Largely the same genes affect variation in life span in these environments. The second group has a single representative: the cross-environment genetic correlation between the high temperature and heat-shock environments is intermediate, and significantly different from both zero and one. Some common as well as different loci affect variation in life span in these environments. Most of the correlations fall into the third category: highly significantly different from one and not significantly different from zero. This group includes all correlations with the starvation environment and all male correlations with the low temperature environment. Completely different sets of loci affect longevity in these environment pairs. The fourth group is somewhat anomalous and includes the correlations of low temperature with the control, high temperature and heat-shock environments in females. These correlations are not significantly different from one or from zero. We regard these correlations as poor estimates; not unexpected given that the among-line variance component for females in the low temperature treatment was not formally significantly different from zero. The differences between cross-environment genetic correlations of ln-transformed compared to untransformed data are quantitative, not qualitative. Control, heat-shock, and high temperature treatments are more highly correlated with each other than with the low temperature treatment (particularly for males), and the correlation between starvation and the other treatments remains low (data not shown). Scale effects (changes in variance) are not the main cause of the observed departures of cross-environment genetic correlations from unity.

Life span QTL:
Cross-sex and cross-environment genetic correlations that are significantly less than one could arise if some loci are conditionally expressed in different sex or treatment environments, or if the effects change sign across environments (a kind of antagonistic pleiotropy). To determine the extent to which loci contributing to GSI and GEI exhibit these properties, we mapped the QTL affecting adult longevity and GSI and GEI for adult longevity in this set of RI lines, using multiple trait composite interval mapping (JIANG and ZENG 1995 ). Several analyses were done. First, we mapped QTL affecting life span and GSI for life span within each treatment environment. The LOD thresholds for significance of the main QTL effect, the conditional test for GSI, and the genome-wide scan for GSI were 3.11, 1.44, and 2.60, respectively. Second, we mapped QTL affecting life span and GEI for life span within each sex across treatment environments. It was necessary to exclude the starvation treatment from these analyses: there were convergence problems, most likely because this trait was not correlated with the others. The LOD thresholds for significance of the main QTL effect, the conditional test for GEI, and the genome-wide scan for GEI were 3.99, 2.46, and 3.56, respectively. A third analysis was performed for females, including only the control, heat-shock, and high temperature treatments, because the genetic variance among lines in the low temperature treatment was not significant in females. The LOD score threshold for the main QTL effect in this analysis was 3.56 and for the genome-wide scan for GEI was 3.11.

All composite interval mapping analyses were done on life span and on ln(life span). The results of these analyses on the untransformed scale are summarized in Table 4 and are depicted graphically in Fig 1 Fig 2 Fig 3 Fig 4 Fig 5 Fig 6. We have presented the results of the analyses on the natural scale for several reasons: (1) the ln-transformation did not improve the fit of the line means to a normal distribution; (2) the genetic variance was greater, relative to the total variance, for the untransformed data, thus giving greater power to detect QTL; and (3) there were few differences between the two sets of analyses. These differences are also summarized below.

View larger version (22K): In this window In a new window Download PPT slide   Figure 1. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of life span in the control environment. The marker positions are indicated by triangles along the abscissa. Results are shown for the joint analysis pooled across sexes and the QTL x sex interaction.

View larger version (22K): In this window In a new window Download PPT slide   Figure 2. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of life span in the heat-shock environment.

View larger version (22K): In this window In a new window Download PPT slide   Figure 3. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of life span in the high temperature environment.

View larger version (19K): In this window In a new window Download PPT slide   Figure 4. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of life span in the low temperature environment.

View larger version (19K): In this window In a new window Download PPT slide   Figure 5. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of life span in the starvation environment.

View larger version (31K): In this window In a new window Download PPT slide   Figure 6. Plot of LOD scores against estimated map position from multiple-trait composite interval mapping of male and female life span. The joint analysis considers the control, heat shock, high temperature, and low temperature environments. Total QTL and QTL x environment interaction LOD scores are shown.

We mapped 29 QTL affecting life span in the separate analyses; however, several QTL map to the same location. Conservatively, we estimate that a total of 17 QTL affect variation in life span in these lines, and have named them Ls1–Ls17 for the purposes of discussion and future reference. QTL detected in different analyses were considered to be the same if the map positions were the same or the interval in which the LOD score exceeded the threshold value overlapped, and the estimated effects were similar. Two QTL were not clearly the same or different based on these criteria. To be conservative, they were classified as being the same as others detected in the gene region, but with the suffix "a" to indicate the ambiguity. A brief description of each QTL follows.

Ls1 has a female-specific main effect, and female-specific GEI considering the control, heat-shock, and high temperature environments. The GEI is partly attributable to conditional expression (there is a highly significant positive effect in the control environment) and antagonistic pleiotropy (the effect in the high temperature environment is negative).

Ls2 is specific to the high temperature environment. The significant GSI effect is caused by antagonistic pleiotropy.

Ls3 is a starvation-specific QTL that is expressed in both sexes.

Ls4 is significant in both the control environment and in females over the four treatments. There is highly significant GSI and GEI due to conditional expression only in females in the control environment.

Ls5 is significant in the control environment and in males over the four treatments. The significant GSI effect is attributable to antagonistic pleiotropy, and the GEI effect is due to conditional expression in males of the control environment.

Ls6 is a female-specific QTL with antagonistic pleiotropic effects in the control and heat-shock environments.

Ls7 has opposite effects in males and females in the heat-shock environment. The GSI effect is significant at the level of a genome wide-scan, but the main QTL effect is not significant. Ls7a, which may be the same QTL, is significant in the starvation environment and has a significant GSI effect attributable to conditional expression in males. The effects in the heat-shock and starvation environments are of opposite sign.

Ls8 is significant in males in the four environment analyses, and in the control environment alone. The significant GEI effect is due to the large conditional effect in the control (and to some extent, the heat-shock) environment. The significant GSI effect is because the effects change sign in males and females. Ls8a maps to the same region and may be the same QTL. It is specific to the starvation environment; this effect is in the opposite direction to the control and heat-shock effects in males.

Ls9 is a QTL with a very large effect in the starvation treatment only, expressed in both sexes.

Ls10 has a main and a GSI effect in the high temperature environment. This QTL has opposite effects in males and females.

Ls11 is expressed in both sexes in the low temperature environment. This large-effect QTL contributes to the main and GEI effect in the female four-treatment analysis, and a GEI effect in the male four-treatment analysis. The GEI effects are attributable to the conditional expression in a single environment.

Ls12 has a large female-specific effect in the control environment, contributing to conditional GSI and conditional GEI in the four-treatment analysis. It also has a significant GSI effect in the high temperature environment, also due to conditional expression in females.

Ls13 is a female-specific QTL with GEI in the four-treatment analysis. The effects of this QTL change sign between the control and high temperature environments.

Ls14 is a male-specific QTL in the four-treatment analysis, with significant GEI due to conditional expression in the control and heat-shock environments. The heat-shock effect is also significant pooled over both sexes. This QTL also has a male-specific effect in the starvation treatment.

Ls15 is significant in the female four-treatment analysis and in the control. It contributes to GEI by antagonistic pleiotropic effects in the control and high temperature environments, and to GSI through sexually antagonistic effects in the control environment.

Ls16 is a female-specific QTL with GEI due to conditional expression in the control and heat-shock environments.

Ls17 is a heat-shock-specific QTL with GSI due to conditional expression in females.

We detected 27 QTL mapping to 18 locations in the analyses of ln(life span). Twenty-one QTL mapping to 14 separate locations were the same in the two analyses. Of the 21 QTL common to both analyses, 17 exhibited the same QTL x environment interactions on both scales. The major differences between the two analyses were as follows:

1. Three highly significant QTL were detected in the analysis of the untransformed data that were not detected in the analysis of ln-transformed data: Ls6, Ls10, and Ls13.

2. The QTL detected between markers 38E–43A in the analysis of the untransformed control data appear between markers 34EF–35B in the analysis of the ln-transformed data. Such shifts in inferred QTL location are common using composite interval mapping analysis because the inferred locations are sensitive to the marker cofactors used.

3. Although QTL mapping to the locations of Ls12 and Ls15 were detected in both sets of analyses, the exact analyses in which the QTL were significant were not the same. The Ls12 QTL detected in the C and female four-treatment (F, 4T) analyses of untransformed data were not detected in the ln-transformed data, and the Ls15 QTL detected in the F, 4T analysis of untransformed data was not detected in the ln-transformed data. All remaining differences between the two analyses were attributable to appearance and disappearance in one or the other analysis of QTL significant for the interaction effect only; these could be false-positive results.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

We have examined the genetic architecture of Drosophila life span by estimating quantitative genetic parameters and mapping QTL for adult life span in a population of RI lines reared in five environments—standard culture conditions, high and low temperature, and heat-shock and starvation stress. Although there was highly significant genetic variation for life span among the RI lines within each sex and environment (with the exception of the female, low temperature analysis on the natural scale), the only significant genetic variation in the analysis pooled over all environments and the two sexes appeared in the GSI and GEI interaction terms. The genetic correlation of life span across sexes and environments was not significantly different from zero.

The expression of QTL affecting longevity in Drosophila is thus highly sensitive to environmental conditions. No QTL were expressed in all environments. Two QTL, Ls4 and Ls17, were conditionally expressed in only one sex and environment. Four sex-specific QTL (Ls1, Ls6, Ls13, and Ls16) were expressed in two environments, and one sex-specific QTL (Ls14) was expressed in three environments. Ls1, Ls6, and Ls13 exhibited antagonistic pleiotropic effects in the different environments, whereas the effects of Ls14 and Ls16 were in the same direction in the different treatments. Five QTL affected both sexes in only one environment: Ls2 and Ls10 had sexually antagonistic effects, and Ls3, Ls9, and Ls11 had similar effects in both sexes and did not contribute to GSI. The remaining five QTL, Ls5, Ls7, Ls8, Ls12, and Ls15, had more complicated patterns of expression in both sexes and two or more environments. Each of these QTL had sexually antagonistic effects in one environment. In addition, Ls5 and Ls12 had sex-specific effects in the same direction in a second environment; and Ls7, Ls8, and Ls15 had antagonistic pleiotropic effects in two or more environments that were expressed in only one sex.

NUZHDIN et al. 1997 mapped QTL for adult longevity using the same set of RI lines. In this experiment, animals were housed individually in small tubes rather than in groups of 10 in standard culture vials; otherwise the environment was equivalent to the control conditions reported here. Five sex-specific QTL were detected, four of which correspond in location to Ls7, Ls8, Ls11, and Ls15.

The use of RI lines is the best experimental design for investigating the genetic basis of variation for traits whose phenotypic expression is highly sensitive to uncontrollable environmental variation, such as life history traits (SOLLER and BECKMANN 1990 ), and is also the optimal method for evaluating QTL x environment interactions. However, like all QTL mapping methods, these analyses are of necessity restricted to assessing the nature of the genetic variation that segregates between the two parental strains, and results from a single experiment beg the question of generality. Will the same QTL be found to segregate between different parental strains? To what extent are sex- and environment-specific effects of longevity QTL a peculiarity of the genetic constitution of the two parental lines used to establish the RI line mapping population studied here? These questions can only be addressed by extensive replication of the experimental design using different starting material or by introgressing gene regions containing the life span QTL from a large natural sample of chromosomes into a common genetic background (LYMAN and MACKAY 1998 ). If the same QTL are found in many parental lines, one would infer that alleles at these loci are at intermediate frequency in the natural population from which the strains were sampled; if not, a rare alleles model (BARTON and TURELLI 1989 ) of segregating variation is supported. There is some evidence that extensive GSI and GEI are general properties of longevity QTL. MAYNARD SMITH 1958 undertook a quantitative genetic analysis of longevity in D. subobscura, and found that genes with sex-limited effects were responsible for variation in life span within and between populations of this species, and that the sex differences were temperature-dependent. Highly environment-specific effects of QTL affecting longevity pose a practical problem for the genetic dissection of this trait, for it is clearly critical to evaluate life span in a wide range of environments to obtain a complete picture of the factors affecting variation in this trait.

The common occurrence of QTL with opposite effects in males and females, and between environments, immediately suggests that genetic variation for life span could be maintained by the balancing selection mechanism proposed by LEVENE 1953 , whereby the fitness effects of two alleles at a single locus change rank in alternate environments. Although the original formulation of the antagonistic pleiotropy hypothesis for aging was in terms of negative genetic correlations between fitness effects of genes expressed early and late in life (i.e., negative correlations between alleles affecting life span and other life history traits), this hypothesis can be readily extended to opposite effects of alleles affecting life span in males and females or in different environments. These observations thus provide support for WILLIAMS' (1957) pleiotropy theory of senescence. An important caveat regarding the generality of inferences drawn from analyses of life history traits in RI lines is that these are nearly totally homozygous genotypes, a contrived genetic situation that does not exist in natural populations of outbreeding species like Drosophila. It is critical in the future to determine homozygous and heterozygous effects of QTL affecting life span, and QTL x environment interaction effects, in a largely heterozygous background genotype. However, if homozygous genotypes exhibit more GSI and GEI than heterozygotes, genetic variation would still be maintained at these loci via the mechanism proposed by GILLESPIE and TURELLI 1989 .

It is interesting to note that opposing QTL effects in males and females and across environments did not lead to significantly negative estimates of genetic correlations between the sexes or environments, although many of the correlations were not significantly different from zero. The genetic correlations represent the summation of effects over all contributing loci and can conceal considerable heterogeneity of individual QTL properties. Such heterogeneity can confound efforts to dissect physiological and genetic mechanisms responsible for variation in longevity based on correlated responses to selection for postponed senescence and may provide a partial explanation for the variable results obtained in different experiments (in addition to real genetic differences among the starting base populations). The long-lived lines of ROSE 1984 had decreased early fertility and increased resistance to starvation and desiccation stress relative to the controls (SERVICE et al. 1988 ). However, lines selected for postponed senescence by PARTRIDGE and FOWLER 1992 did not differ in early female fertility from the control lines, and the long-lived line of LUCKINBILL et al. 1984 exhibited decreased early fertility and no difference in resistance to starvation and desiccation relative to the unselected line (ARKING 1987 ; FORCE et al. 1995 ). These conflicting results can only be resolved by analyses at the level of individual QTL. For example, the genetic correlation between starvation resistance and life span in the other environments in the population of RI lines studied here was not significantly different from zero, and most QTL affecting variation in starvation resistance were correspondingly specific to this treatment. However, one QTL, Ls14, exhibited a positive correlation between starvation resistance and life span in the control and heat-shock environments.

Based on the results reported here, different environmental conditions might be expected to yield variable correlated responses to selection, given that the effects of the same genotype on life span will vary according to the environment. Indeed, BUCK et al. 1993 report that their long-lived phenotype is conditional on high larval density, and LEROI et al. 1994 note that the genetic correlation between life span and early fecundity in Rose's long-lived and control lines can vary according to the conditions under which the phenotypes are assessed.

Many candidate genes affecting life span have been postulated: genes involved in resistance to heat shock and oxidative stress, DNA repair, cellular aging, metabolic energy storage, and loci with sex-specific effects on fertility and reproduction. With so many candidate genes and such a large fraction of the genome covered by QTL affecting life span, it is expected that many candidate gene and QTL map positions will overlap by chance. While fine-scale mapping is necessary before positional candidates can be proposed for further analysis, it is nevertheless interesting to note that life span QTL map to the same locations as genes encoding enzymes of carbohydrate metabolism (Phosphogluconate dehydrogenase, Ls1; Phosphoglucose isomerase, Ls8a); oxidative and heat-shock stress (Sod and the small heat shock proteins, Ls11; Heat shock protein 70A; Ls16); sex determination (Sex-lethal, Ls4; sisterless-a, Ls5); and sex-specific peptides (accessory gland proteins B (AcpB), Ls7; AcpK, Ls8).

	FOOTNOTES

¹ Present address: UMR 5558, Laboratoire de Biométrie, Génétique et Biologie des Populations, Université Claude Bernard Lyon1, 43 bd. 11 novembre, 69622 Villeurbanne Cedex, France.

	ACKNOWLEDGMENTS

We thank R. Anholt, J. Leips, and M. Wayne for comments on the manuscript. This work was supported by National Institutes of Health (NIH) grant GM 45146 to Z.-B.Z. and T.F.C.M., NIH GM 45344 to T.F.C.M., a Portuguese postdoctoral fellowship to C.V., grant 97-04-48101 from the Russian Fund of Basic Research to E.G.P., and NIH TW00997 to T.F.C.M. and E.G.P. This is a publication of the W. M. Keck Program for Behavioral Biology.

Manuscript received January 11, 1999; Accepted for publication September 7, 1999.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

ARKING, R., 1987  Successful selection for increased longevity in Drosophila: analysis of the survival data and presentation of a hypothesis on the genetic regulation of longevity. Exp. Gerontol. 22:199-220[Medline].

BARTON, N. H. and M. TURELLI, 1989  Evolutionary quantitative genetics: how little do we know? Annu. Rev. Genet. 23:337-370[Medline].

BODNAR, A. G., M. OUELLETTE, M. FROLKIS, S. E. HOLT, and C.-P. CHIU et al., 1998  Extension of life-span by introduction of telomerase into normal human cells. Science 279:349-352[Abstract/Free Full Text].

BUCK, S., M. NICHOLSON, S. DUDAS, R. WELLS, and A. FORCE et al., 1993  Larval regulation of adult longevity in a genetically-selected long-lived strain of Drosophila.. Heredity 71:23-32.

BULMER, M. G., 1985 The Mathematical Theory of Quantitative Genetics. Clarendon Press, Oxford.

BURTON, V., H. K. MITCHELL, P. YOUNG, and N. S. PETERSEN, 1988  Heat shock protection against cold stress of Drosophila melanogaster.. Mol. Cell. Biol. 8:3550-3552[Abstract/Free Full Text].

CHARLESWORTH, B., 1980 Evolution in Age-Structured Populations. Cambridge University Press, London.

DECKERT-CRUZ, D. J., R. H. TYLER, J. E. LANDMESSER, and M. R. ROSE, 1997  Allozymic differentiation in response to laboratory demographic selection of Drosophila melanogaster.. Evolution 51:865-872.

ENGSTRÖM, G., L.-E. LILJEDAHL, M. RASMUSON, and T. BJÖRKLUND, 1989  Expression of genetic and environmental variation during ageing. 1. Estimation of variance components for number of adult offspring in Drosophila melanogaster. Theor. Appl. Genet. 77:119-122.

ENGSTRÖM, G., L.-E. LILJEDAHL, and T. BJÖRKLUND, 1992  Expression of genetic and environmental variation during ageing. 2. Selection for increased life span in Drosophila melanogaster. Theor. Appl. Genet. 85:26-32.

FALCONER, D. S., 1952  The problem of environment and selection. Am. Nat. 86:293-298.

FALCONER, D. S., and T. F. C. MACKAY, 1996 Introduction to Quantitative Genetics, Ed. 4. Longman, Harlow, Essex, United Kingdom.

FINCH, C. E., 1990 Longevity, Senescence and the Genome. University of Chicago Press, Chicago.

FINCH, C. E. and R. E. TANZI, 1997  Genetics of aging. Science 278:407-411[Abstract/Free Full Text].

FORCE, A. G., T. STAPLES, S. SOLIMAN, and R. ARKING, 1995  Comparative biochemical and stress analysis of genetically selected Drosophila strains with different longevities. Dev. Genet. 17:340-351[Medline].

GILLESPIE, J. H. and M. TURELLI, 1989  Genotype-environment interactions and the maintenance of polygenic variation. Genetics 121:129-138[Abstract/Free Full Text].

GRAVES, J. L., E. C. TOOLSON, C. JEONG, L. N. VU, and M. R. ROSE, 1992  Desiccation, flight, glycogen, and postponed senescence in Drosophila melanogaster.. Physiol. Zool. 65:268-286.

GURGANUS, M. C., J. D. FRY, S. V. NUZHDIN, E. G. PASYUKOVA, and R. F. LYMAN et al., 1998  Genotype-environment interaction for quantitative trait loci affecting sensory bristle number in Drosophila melanogaster.. Genetics 149:1883-1898[Abstract/Free Full Text].

HARRISON, D. E., 1990 Genetic Effects on Aging II. The Telford Press, Caldwell, NJ.

HARRISON, D. E. and J. R. ARCHER, 1987  Genetic differences in effects of food restriction on aging in mice. J. Nutr. 117:376-382.

HEKIMI, S., B. LAKOWSKI, T. M. BARNES, and J. J. EWBANK, 1998  Molecular genetics of life span in C. elegans: how much does it teach us? Trends Genet. 14:14-19[Medline].

HOULE, D., 1992  Comparing evolvability and variability of quantitative traits. Genetics 130:195-204[Abstract].

HOULE, D., B. MORIKAWA, and M. LYNCH, 1996  Comparing mutational variabilities. Genetics 143:1467-1483[Abstract].

HUGHES, K. A., 1995  The evolutionary genetics of male life-history characters in Drosophila melanogaster.. Evolution 49:521-537.

HUGHES, K. A. and B. CHARLESWORTH, 1994  A genetic analysis of senescence in Drosophila.. Nature 367:64-66[Medline].

JIANG, C. and Z.-B. ZENG, 1995  Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics 140:1111-1127[Abstract].

KHAZAELI, A. A., M. TATAR, S. D. PLETCHER, and J. W. CURTSINGER, 1997  Heat-induced longevity extension in Drosophila. I. Heat treatment, mortality, and thermotolerance. J. Gerontol. 52:B48-B52.

KIMURA, K. D., H. A. TISSENBAUM, Y. LIU, and G. RUVKUN, 1997  daf-2, an insulin receptor-like gene that regulates longevity and diapause in Caenorhabditis elegans.. Science 277:942-946[Abstract/Free Full Text].

KIRKWOOD, T. B. L. and M. R. ROSE, 1991  Evolution of senescence: late survival sacrificed for reproduction. Philos. Trans. R. Soc. Lond. B Biol. Sci. 332:15-24[Medline].

LAI, C., R. F. LYMAN, A. D. LONG, C. H. LANGLEY, and T. F. C. MACKAY, 1994  Naturally occurring variation in bristle number and DNA polymorphisms at the scabrous locus of Drosophila melanogaster.. Science 266:1697-1702[Abstract/Free Full Text].

LEROI, A. M., A. K. CHIPPENDALE, and M. R. ROSE, 1994  Long-term laboratory evolution of a genetic life-history trade-off in Drosophila melanogaster. I. The role of genotype-by-environment interaction. Evolution 48:1244-1257.

LEVENE, H., 1953  Genetic equilibrium when more than one ecological niche is available. Am. Nat. 87:331-333.

LIN, K., J. B. DORMAN, A. RODAN, and C. KENYON, 1997  daf-16: an HNF-3/forkhead family member that can function to double the life-span of Caenorhabditis elegans. Science 278:1319-1322[Abstract/Free Full Text].

LIN, Y.-J., L. SEROUDE, and S. BENZER, 1998  Extended life-span and stress resistance in the Drosophila mutant methuselah.. Science 282:943-946[Abstract/Free Full Text].

LINDSLEY, D. L., and G. G. ZIMM, 1992 The Genome of Drosophila melanogaster. Academic Press, San Diego.

LITHGOW, G. J., 1996  Invertebrate gerontology: the age of mutations of Caenorhabditis elegans.. BioEssays 18:809-815[Medline].

LITHGOW, G. J., T. M. WHITE, S. MELOV, and T. E. JOHNSON, 1995  Thermotolerance and extended life-span conferred by single-gene mutations and induced by thermal stress. Proc. Natl. Acad. Sci. USA 92