Quantitative genetic analysis of life-history traits of Caenorhabditis elegans in stressful environments
BMC Evolutionary Biology volume 8, Article number: 15 (2008)
Organisms live in environments that vary. For life-history traits that vary across environments, fitness will be maximised when the phenotype is appropriately matched to the environmental conditions. For the free-living nematode Caenorhabditis elegans, we have investigated how two major life-history traits, (i) the development of environmentally resistant dauer larvae and (ii) reproduction, respond to environmental stress (high population density and low food availability), and how these traits vary between lines and the genetic basis of this variation.
We found that lines of C. elegans vary in their phenotypic plasticity of dauer larva development, i.e. there is variation in the likelihood of developing into a dauer larva for the same environmental change. There was also variation in how lifetime fecundity and the rate of reproduction changed under conditions of environmental stress. These traits were related, such that lines that are highly plastic for dauer larva development also maintain a high population growth rate when stressed. We identified quantitative trait loci (QTL) on two chromosomes that control the dauer larva development and population size phenotypes. The QTLs affecting the dauer larva development and population size phenotypes on chromosome II are closely linked, but are genetically separable. This chromosome II QTL controlling dauer larva development does not encompass any loci previously identified to control dauer larva development. This chromosome II region contains many predicted 7-transmembrane receptors. Such proteins are often involved in information transduction, which is clearly relevant to the control of dauer larva development.
C. elegans alters both its larval development and adult reproductive strategy in response to environmental stress. Together the phenotypic and genotypic data suggest that these two major life-history traits are co-ordinated responses to environmental stress and that they are, at least in part, controlled by the same genomic regions.
Organisms live in environments that vary both spatially and temporally. In such variable environments there are different ways to maximise fitness. Life-history traits can either be robust to environmental change (a buffered or canalised trait) or they can be variable in an environmentally-dependent manner (a phenotypically plastic trait). Phenotypic plasticity of a trait can be manifest as a continuous phenotypic range across an environmental gradient, such as the variation in Drosophila melanogaster body size metrics across temperature ranges . Alternatively, phenotypic plasticity may appear as a threshold trait, with apparently distinct phenotypes developing in different environments. An example of this is the switch between winged and wingless aphid morphs in response to host plant quality and, or aphid population density . These different phenotypic responses have been termed phenotypic modulation and developmental conversion, respectively .
A priori, fitness could be maximised by all traits being fully phenotypically plastic. However, phenotypically plastic traits vary both within and between populations, particularly in the magnitude and sensitivity of their response to environmental change: in the language of phenotypic plasticity, there may be different reaction norms. The existence of this variation suggests that there are limits or costs to the evolution of phenotypically plastic traits and of the reaction norms of traits, and therefore that fitness is maximised by not all traits being fully phenotypically plastic. These costs are likely to be (i) having sufficiently accurate and robust processes for environmental sensation, (ii) maintaining the genetic and cellular machinery for the development of alternate phenotypes and (iii) co-ordination between different phenotypically plastic traits [4–6]. Therefore, all traits can be considered on a continuum from those that are fully plastic, via those with a low level plasticity, to non-plastic, invariant traits. The molecular basis of the phenotypic plasticity of most traits is not clear, but progress in identifying genes involved in such environmental interactions is being made (e.g [7–10]). For many organisms, including intensively studied 'model' species, the role and function of the majority of genes remains unknown [11, 12]. It is probable that genes involved in phenotypically plastic traits, especially the means by which the phenotype is modulated in response to the environment, are among these genes with, as yet, unidentified functions. Given this, it is crucial to move towards integrating an understanding of the molecular basis of phenotypic plasticity with the ecology of the species in question .
The model free-living nematode Caenorhabditis elegans has a phenotypically plastic developmental switch in its life-cycle. In the 'normal' C. elegans life-cycle, progeny moult through four larval stages (L1 – L4) into adult worms. Under conditions that are unsuitable for growth and reproduction an alternative, developmentally arrested third larval stage, the dauer larva, is formed [14, 15]. Dauer larvae are environmentally resistant, have a greatly extended lifespan compared with non-dauer larvae, and will only resume normal development as fourth stage larvae (L4) when exposed to more favourable conditions. This developmental conversion (dauer vs. non-dauer) is made based on the worms' perception of environmental quality, which is determined by the concentrations of both food and dauer pheromone and temperature. Dauer pheromone is a cue produced by all worms that acts as a measure of con-specific population density  and appears to consist of three related molecules [17, 18]. The relative concentrations of food and dauer pheromone are used to assess environmental quality: dauer larvae develop under low food and high dauer pheromone concentration conditions, i.e. conditions of environmental stress; 'normal', non-dauer, larvae develop under high food and low dauer pheromone concentration conditions, i.e. conditions of plenty.
Geographically distinct isolates of C. elegans vary extensively in the phenotypic plasticity of dauer larva formation . That is, for the same change in environmental conditions, there is variation in the proportion of larvae that develop into dauer larvae. The adaptive value of these different plasticities of dauer larva formation is not known. There has been extensive investigation into the genetic and molecular genetic control of the development of dauer larvae [14, 15]. However, the genetic basis of the variation between isolates in the phenotypic plasticity of dauer larva development is not known. Furthermore, the fitness consequences of different phenotypic plasticities of dauer larva development are not known.
In the C. elegans life-cycle there are other phenotypic responses to conditions of environmental stress. For example adult lifespan, lifetime fecundity, the schedule of reproduction and body size all vary as a consequence of food quantity and quality [20, 21]. At present little is known about variation between C. elegans isolates in these genotype by environment (GxE) interactions and it is not known if they represent separate trait-specific responses to environmental change or are part of a co-ordinated response. Hence, it is possible that the fitness of a C. elegans genotype may be maximised by a co-ordinated response to environmental stress of both larval and adult traits. Further, it can be envisaged that different strategies may have evolved for how such major larval and adult life-history traits co-ordinately respond to environmental change.
To investigate this, we have analysed recombinant inbred lines (RILs) produced from C. elegans isolates N2 and DR1350, which differ significantly in their plasticity of dauer larvae formation . Using these, we determined how the development of dauer larvae and the population growth, as assayed by determining population size through time, varied between the RILs. Population size at any given time is a consequence of both survival and all the factors that contribute to reproduction (e.g. time to first reproduction, total lifetime fecundity, rate of reproduction) this broad analysis was most likely to identify any relationship with dauer larvae development. We found that the C. elegans RILs vary in both their phenotypic plasticity of dauer larva development and their population size and that these traits are related. Specifically, we found that lines of C. elegans that are highly plastic for dauer larva development also maintain a high population growth rate when stressed (i.e. their rate of reproduction is comparatively unchanged under conditions of environmental stress). This suggests that C. elegans may co-ordinate its larval development and adult reproductive strategy to respond to environmental stress. We have also mapped the quantitative trait loci (QTL) that control the phenotypic plasticity of dauer larva development and of population size. These analyses have identified QTL that affect dauer larva development and population size.
C. elegans varies in its phenotypic plasticity of dauer larva formation
The proportion of dauer larvae that developed in 163 RILs at two different (1% and 2% w/v) food concentrations, i.e. the reaction norm of dauer larvae development, is shown in Figure 1A; for clarity these data from one representative assay are shown in Figure 1B. RILs that are highly plastic for dauer larva formation have the greatest positive difference in the proportion of dauer larvae formed at low food concentrations, compared with high food concentrations (i.e. the steepest, positive, reaction norms); N2 is an example of such a line (Figure 1B). Other lines are less plastic; DR1350 is an example of such a line (Figure 1B). Approximately a quarter of the RILs (45/163) had a negatively plastic dauer larva development phenotype (i.e. a reaction norm with a negative slope); that is, a lower proportion of dauer larvae developed at a low food concentration, compared with a high food concentration (Figure 1A). In many of these 45 lines, the difference in dauer larva formation between the two food concentrations was very small and, given the inherent variation in this trait and its assay, it is likely that these lines have only a small, or no, plasticity of dauer larva development. However, ten lines had a negative plasticity ≥ 0.2, and were found to consistently show this response in additional assays (data not shown).
As expected dauer larva formation was significantly affected by food concentration (FOOD: F1,542 = 4.64, p = 0.048). There was significant variation between the RILs (LINE [ASSAY]: F153,542 = 6.43, p < 0.001) and the RILs also varied significantly in their dauer larva formation phenotype in response to changes in the concentration of food (FOOD*LINE [ASSAY]: F153,542 = 1.55, p < 0.001), i.e. there is genetic variation between the RILs for both overall propensity to form dauer larvae and in their plasticity of dauer larva development with respect to food concentration. There were significant differences between the assays, and in the effect of food concentration between the assays (ASSAY: F9,542 = 6.43, p < 0.001; FOOD*ASSAY: F9,542 = 3.28, p < 0.001).
The dauer larva development phenotype with respect to dauer pheromone concentration had previously been determined for 35 of these RILs . We used these data to compare the dauer larva development of these lines with respect to changes in the concentration of food and of pheromone. This showed that there was a significant positive correlation between these plasticities of dauer larvae development (r2 = 0.40, p = 0.02). This, therefore, suggests that the plasticity of dauer larva development is, at least partially, independent of the particular environmental cue (i.e. food or pheromone concentration) used for the initiation of dauer larva development.
C. elegans varies in its population growth
We wished to determine how C. elegans varied for another important life-history trait, namely adult reproduction. To do this we determined the population size of 45 RILs as they grew. This showed that the population sizes of the lines increased over time, but that the rate of this increase differed between the lines (Figure 2A; for clarity these data from one representative assay are shown in Figure 2B). There were significant differences between the population sizes of the lines and in how they grew over time (TIME: F2,1039 = 23431.01, p = 0.048; LINE [ASSAY]: F37,1039 = 9.98, p < 0.001; TIME*LINE [ASSAY]: F74,1039 = 4.14, p < 0.001). There were also significant differences between the assays, and in how the populations grew between assays (ASSAY: F7,1039 = 130.84, p < 0.001; TIME*ASSAY: F14,1039 = 102.49, p < 0.001). Hence, the RILs vary in their population sizes and in how the populations grew over time. At the start of these assays (day 5) there are no limitations to population growth; as the size of the population increases (day 7), food conditions become limiting and by day 9, the food source is exhausted. We therefore posited that the variation between the RILs was due to both differences in the reproductive capacity of the lines (the effect of LINE), possibly manifest either by variation in their lifetime fecundities or their temporal schedule of reproduction, and by differences in how the reproductive capacity of the lines altered in response to density (the effect of TIME*LINE [ASSAY]). These possible explanations were investigated further.
C. elegans varies in its fecundity
One of the principal features of high density environments is a reduction in the per capita food availability. To investigate how the fecundity of C. elegans was affected by this we firstly determined the lifetime fecundity of 12 RILs under ad libitum food conditions. This showed that these RILs had significantly different lifetime fecundities (LINE [BLOCK]: F10,222 = 10.35, p < 0.001); fecundity did not differ between experimental blocks (BLOCK: F1,222 = 3.14, p = 0.078). We calculated the rate of reproduction, q, a measure of the rate of increase in the proportion of lifetime fecundity, of the RILs and found that this also differed significantly between RILs (LINE [BLOCK]: F10,222 = 7.94, p < 0.001) and between experimental blocks (BLOCK: F1,222 = 131.46, p < 0.001).
These data were used to calculate the expected population size of the lines in the absence of any limits on population growth on days 5, 7 and 9. In all cases, the population size was greater than that which was observed. These predicted population sizes correlated with the observed population size on day 7 (r2 = 0.548, p = 0.002), but not on days 5 and 9 (r2 = 0.097, p = 0.26 and r2 = 0.12, p = 0.207, respectively). This may suggest that at day 7, growing populations of C. elegans are not constrained by density-dependent effects, but at day 9 that they are. The absence of a significant correlation on day 5 may be due to inter-RIL variation in the rate of reproduction.
We next considered how lifetime fecundity and the rate of reproduction, and hence the expected population sizes of the lines, was affected by low concentrations of food, a likely feature of high-density populations. Analogously, as above, the lifetime fecundity differed between the lines (LINE: F5,223 = 10.85, p < 0.001) and was significantly negatively affected by food concentration (FOOD: F3,223 = 35.63, p < 0.001). But, the lines did not differ in how their lifetime fecundities were affected by different food concentrations (LINE*FOOD: F15,223 = 1.30, p = 0.20). The rate of reproduction, q, did not differ between lines (LINE: F5,223 = 1.33, p = 0.25), but it was significantly affected by the food concentration (FOOD: F3,223 = 8.80, p < 0.001), such that at lower food concentrations q was lower. The lines also differed significantly in how q was affected by food concentration (FOOD*LINE: F15,223 = 2.41, p = 0.003). Together, these results therefore show that the rate of reproduction was affected by food concentration independently of effects on lifetime fecundity.
Why does low food concentration reduce fecundity?
The fecundity of wild-type C. elegans is reduced when there is limited food available. We sought to investigate whether this occurs due to (i) direct effects of caloric restriction per se or (ii) indirect effects on reproduction due to the perception of food availability, i.e. altered resource allocation in response to perceived food limitation. To do this we compared reproduction under a range of food concentrations in mutants defective in pharyngeal pumping with reproduction in wild-type worms. This was done with strains with mild pharyngeal defects (DA573 eat-14(ad573) and DA602 eat-15(ad602)), with severe pharyngeal defects (DA522 eat-13(ad522) and DA606 eat-10(ad606))  and N2 as a control. Here, as far as it is known, the perception of food by all strains will be the same, but the actual acquisition of food will be compromised in the strains with pharyngeal defects, and more strongly so in those with the more severe defects. i.e. the comparison allowed the separation of the acquisition of food from the perception of food. We hypothesised that if caloric restriction (i, above) controlled the food concentration-dependent change in the rate of reproduction, q, the rank order of the food-dependent reduction in q would be eat-10 and eat-13; eat-14 and eat-15; N2. In contrast, if food perception (ii, above) controlled the food concentration-dependent reduction in the rate of reproduction, q, then no relationship between strength of the pharyngeal defects and food-dependent reduction in fecundity would be observed.
The mutants with the strongest pharyngeal defects (eat-10, eat-13) had the lowest rate of reproduction, q, and this was reduced most under lower food concentrations, compared with mutants with mild pharyngeal defects (eat-14 and eat-15) and compared with N2 (Table 1). This therefore suggests that the food concentration-dependent reduction in the fecundity of C. elegans is directly due to caloric restriction.
To identify markers that would differentiate N2 and DR1350, 144 genomic regions of DR1350 were sequenced and a further 49 candidate polymorphisms were screened by restriction fragment length polymorphism (RFLP) analysis (Table 2 and additional file 1). Overall, 45 polymorphisms were detected that differentiated N2 and DR1350; this included 20 polymorphisms previously identified . 42 of these 45 markers were used to genotype the RILs (the remaining three markers were developed at a later date and used only in the fine-mapping, see below). These markers cover chromosomes I, II, III and X (Table 2). In total, 15.45 Kbp of chromosome IV and 12.60 Kbp of chromosome V of DR1350 were sequenced, but no polymorphisms between N2 and DR1350 were detected (see additional file 1 for details).
188 RILs were genotyped at 42 loci. Ten of these RILs were found to have an N2 genotype at all loci and these were excluded from all analyses. For the remaining 178 RILs, the measured genetic distance between markers was greater than the F2-derived genetic maps , with a map expansion of 3.45, 2.73, 2.43 and 3.21 for chromosomes I, II, III and X, respectively. On the X chromosome a group of three markers (2, 3 and 4) and a pair of markers (5 and 6) (Table 2) could not be separated genetically and so the non-informative markers were removed from all further analyses. Analysis of the 39 remaining markers revealed extensive segregation distortion among the RILs, with the overall allele frequency of the 39 markers in the 178 RILs biased to N2 (65%) and a variable rate of segregation distortion across chromosomes. This implies that there was selection for the N2 alleles at many of the loci during the construction of the RILs, a finding common in many studies analysing C. elegans RILs (e.g. [9, 24, 25]).
Dauer larva development QTLs
Comparative interval mapping (CIM) of dauer larvae formation in response to different food concentrations in 153 RILs identified six QTLs: two each for the plasticity of dauer larva formation, and for dauer larva formation at high and low food concentrations (Figure 3, Table 3). Five of these QTL occurred on the X chromosome, one, (for dauer larva formation at low food concentration) occurred on chromosome II. Multiple trait CIM (MT-CIM) analysis of the dauer larva formation at high (2% w/v) and at low (1% w/v) food concentrations supported the QTLs identified by CIM for dauer formation at low food concentration and one of the QTLs for dauer formation at high food concentration (marker 2 on the X chromosome from Table 3). Similarly, one of the QTLs identified by CIM for the plasticity of dauer larvae formation (marker 2 on the X chromosome from Table 3) was supported the MT-CIM analysis (significant GxE interaction detected in that area of the X chromosome). CIM and MT-CIM analyses of dauer larvae formation in response to different pheromone concentrations in 35 RILs identified no significant QTLs. However, there were significant single marker effects predominantly on the X chromosome (Table 3).
The chromosome II QTL affecting dauer larva formation at low food concentration had only a relatively small phenotypic effect (Table 3). However, analysis of dauer development in nearly isogenic lines (NILs), in which the chromosome II QTL region of the DR1350 genome was introgressed into N2, suggests that this QTL has a much greater effect and that it affects the dauer larva formation phenotype across a range of food concentrations. This is shown most clearly in the comparison of NIL-1 to N2 and DR1350 over a range of food concentrations (Figure 4).
Population size QTLs
CIM analysis of the population size at 5, 7 and 9 in 35 RILs identified two QTL each for population size on day 5 and on day 9, three of which were on chromosome II and one on chromosome X (Figure 3, Table 3). MT-CIM analysis of population size supported one of the day 5 population size QTLs on chromosome II (marker 2 on chromosome 2 from Table 3). No significant QTLs were detected by CIM for population size on day 7 nor were any significant single marker associations identified.
Combined, the analyses of dauer development and population size have therefore identified two regions of the C. elegans genome that affect both traits. Firstly, there is a region of the X chromosome near marker 2 that significantly affects both the population size on day 9 and dauer larva formation in response to changes in both the food and pheromone concentrations, by CIM and single marker analyses respectively. Comparison of these phenotypes among the 21 RILs for which both phenotypes had been measured showed that the population size on day 9 was positively correlated with the plasticity of dauer larva formation in response to pheromone concentration (r2 = 0.55, p = 0.01), but not with the plasticity in response to food concentration (r2 = 0.10, p = 0.57). Secondly, there is a region on chromosome II near marker 11 (see below) that significantly affects both the population size on day 5 and dauer larva formation at low food concentrations. However, these traits do not correlate (r2 = 0.10, p = 0.57).
We undertook further mapping to investigate the extent of the co-incidence of the chromosome II QTLs affecting both population size and dauer larva formation.
Dauer larva development: chromosome II further mapping
Among the RILs, those with a DR1350 genotype at chromosome II markers 11 and 13 (i.e. RILs with a DR1350 genotype at the chromosome II QTL) formed a significantly lower proportion of dauer larvae compared with those lines either with an N2 genotype at both markers and recombinants (i.e. RILs with a DR1350 genotype at marker 11 and an N2 genotype at marker 13) (F2,160 = 11.25, p < 0.001) (Figure 5A). This therefore indicates that the loci underlying the chromosome II QTL are not located between markers 9 and 11 and lie to the right of marker 11 on chromosome II (Figure 6A).
Analysis of dauer larva formation in NILs indicated that at a low food concentration (1% w/v), a significantly lower proportion of dauer larvae developed in NILs-9, -68 and -110 and DR1350 than in NIL-10 and N2 (F5,87 = 8.42, p < 0.001) (Figure 5B), i.e. the presence of the introgressed DR1350 DNA is sufficient to replicate the DR1350 dauer development phenotype. The genotype data for these NILs therefore defines the DR1350 allele of the chromosome II QTL affecting dauer development to lie between markers 9 and 12 on chromosome II (Figure 6B). This combined with the previous definition of the left extremity of the QTL (Figure 5A) therefore resolves the QTL to a 208,662 bp region between markers 11 and 12 (3.66 – 3.87 Mbp) (Figure 6).
This chromosome II QTL region in N2 is predicted to contain 72 genes and 10 pseudogenes (see additional file 2A for details) . The RNAi phenotypes, gene ontology (GO) terms, and both the spatial and temporal patterns of expression of these genes do not present any obvious candidates for the control of the dauer larva development phenotype. However, thirteen of these 72 predicted genes and 7 of the 10 predicted pseudogenes are predicted to be 7-transmembrane receptors . The role of these receptors in information transduction, may suggest that a polymorphism in one of these genes or pseudogenes, such that an N2 pseudogene is active in DR1350, may underlie the differences in dauer larva development of these lines.
The 5 QTL identified on the X chromosome could not be defined to small genomic regions. The combined confidence limits for these QTL encompass approximately a third of the X chromosome (5.5 – 12.3 Mbp). Within this region there are a number of genes known to be involved in dauer larva development, e.g. daf-9, daf-12, dyf-6 and dyf-12 [11, 23]. The phenotypic effect of the X chromosome QTL may therefore be due to one of these genes and, or other loci within this region.
Population size: chromosome II further mapping
Analysis of the population size trait indicated that on day 5, the population sizes of N2 and NIL-110 were significantly greater than those of NILs-1 and -9, which were in turn greater than that of DR1350 (F4,96 = 17.86, p < 0.001) (Figure 5C). The genotype data for these NILs therefore resolves the DR1350 allele affecting population size on day 5 to a 657,795 bp region between markers 12 and 14 (3.87 – 4.53 Mbp). In N2, this region is predicted to contain 116 genes and 2 pseudogenes (see additional file 2B for details) . The RNAi phenotypes, GO terms and both the spatial and temporal patterns of expression of these genes do not suggest any obvious candidates for the control of the population size phenotype.
Extensive variation in both the phenotypic plasticity of dauer larva development and population size phenotypes among the RILs was observed. For both phenotypes, there were RILs with phenotypes more extreme than the parental lines (Figures 1 and 2). Indeed, such transgressive segregation appears to be common in C. elegans and has been observed for a number of other life-history traits (e.g. [9, 24, 25]). The variation in dauer larva development between RILs was similar in magnitude to that previously observed among a small number of RILs . However, we observed that approximately 5% of the lines had substantial negative plasticity of dauer larva formation; that is, dauer larva formation was greatest at the higher food concentration. A similar phenomenon was observed in the phenotypic analysis of synthesised daumone, where across a range of food concentrations, negative phenotypic plasticity of dauer larva formation was observed at very low food concentrations, but beyond a mid-concentration of food, there was positive phenotypic plasticity . We hypothesise that in the RILs studied here this previously reported negative plasticity phenotype  has become manifest for some the RILs because of a shift in these lines of their sensitivity to the concentration of food.
The adaptive value of the C. elegans dauer larva and the variation in its formation is unknown, largely due to our limited knowledge of the natural history of this species. In compost-rich soils, C. elegans is mostly found as dauer larvae rather than as reproducing adults . The comparative abundance of this life-cycle morph in the wild therefore suggests its central importance for C. elegans and that the control of the formation of dauer larvae is likely to be under strong natural selection. This suggests that variation in the plasticity of dauer larva formation may also be under similarly strong selection and hence that this variation may be adaptive.
The observed population sizes of the RILs were smaller than those predicted from measures of lifetime fecundity and reproductive schedule. In these assays on day 9 the food supply is exhausted which suggests that by this time the growing population is subject to density-dependent effects. The C. elegans rate of reproduction is lower in low food concentration settings. Comparison of the rate of reproduction of eat mutants suggests that the rate of reproduction is determined directly by calorific restriction, thus, that density-dependent effects on reproduction come about by such caloric restriction. This trait is therefore different from the dauer larva developmental choice, in which a developmental decision is made based primarily on the multi-factor perceived environmental conditions. That is, the density-dependent effects on reproduction are a consequence of the internal state of the worm while the effects on dauer development are a consequence of the external environment.
In summary, under conditions of reduced food availability, a likely significant environmental stress, larval stages of C. elegans form dauer larvae and the reproduction of the adult hermaphrodites is altered, which reduces the rate of reproduction. Therefore, these different life-cycle stages both have strategies for responding to this environmental stress and that these larval and adult hermaphrodite responses to this stress vary between lines. There are consistent differences in C. elegans life-history traits when grown on two different food sources. On chemically defined liquid media, both development and lifetime fecundity are reduced and lifespan increased compared with worms grown on an NGM/OP50 food source [27–29]. These results suggest the existence of at least two distinct adult life histories: one maximizing the intrinsic rate of population increase (NGM/OP50 food), and the other the efficiency of exploitation of the carrying capacity of the environment (chemically defined liquid media) . Our observations that C. elegans varies in its reproductive response to altered food conditions, is therefore in general concordance with these observations.
We detected substantial genetic polymorphism between N2 and DR1350 for much of the genome. However, we did not detect any genetic differences between these isolates on chromosomes IV and V nor in ~5 Mbp of the X chromosome. This supports previous observations  that also did not detect any polymorphisms between N2 and DR1350 in several regions of the genome, including those noted here.
Our genetic mapping analyses identified 10 QTL that affect either dauer larva development or population size (Figure 3 and Table 3). There are two principal regions of the genome, on chromosomes II and X, that affect both dauer larva development and population size. The dauer larva development and population size QTLs identified on the X chromosome could not be separated and, as would be expected given the position and effects of the day 9 population size and dauer development QTLs on the X chromosome (Table 3), these traits are positively correlated in the RILs. However, given the limited resolution of these QTLs little can be inferred about the genes underlying these QTLs. In contrast, detailed analysis of the chromosome II region separates the day 5 population size and dauer larva development phenotypes (Figure 5) showing that they are not controlled by the same loci. Crucially, the chromosome II region that can control dauer larva formation does not include any genes previously identified to affect this phenotype; therefore a new locus or loci that affects dauer larva development has been discovered. This locus or loci may not have been identified previously by more usual mutagenesis analyses, perhaps because there may not be a severe (or any) dauer formation (daf) loss of function phenotype.
The QTLs for population size appeared, on average, to have greater phenotypic effects compared with those affecting dauer larva development (Table 3). However, for the chromosome II dauer larvae development QTL, the phenotypic effect size is greater in the NIL analysis compared with the QTL analysis (Figure 4 and 5C) and this QTL affects dauer development across a range of food concentrations (Figure 4). The difference between the effect size observed in the RILs and the NILs could have a number of explanations. It is unlikely that there are epistatic interactions between the DR1350 introgressed region and the N2 genome, since no such interactions were detected between QTLs. Similarly, we consider it unlikely that the difference between the estimated phenotypic effect size of the QTL in the RILs and the observed effect size in the NILs is due to the separation in the NILs of loci with opposing effects on dauer larva formation given the relatively large size (10.0 – 11.3 Mbp of chromosome II, data not shown) of the DR1350 region introgressed in NIL-1. The greater than expected phenotypic effect of the QTL in the NILs may therefore be an underestimate due to variation in trait measurement in these analyses (Figure 1). Alternatively, the increased effect of the QTL in the NILs may be a consequence of the use of daumone, rather than pheromone extract in the assay of the NILs. Because daumone is only one of three active compounds in dauer pheromone , a larger observed effect would be expected in the NIL analyses if the chromosome II QTL represented a polymorphism in a gene specifically involved in the response to daumone. Indeed, this would be consistent with the suggestion that no candidate dauer pheromone receptors have been isolated because different receptors respond to the different active compounds in the pheromone . Hence, the disruption of the receptor for one component of dauer pheromone would not be expected to produce a dauer defective phenotype. Within this region there are many predicted 7-transmembrane receptors [11, 23], and such receptors have been implicated in the response to dauer pheromone by the finding that the G-protein α-subunits GPA-1 and -2 are involved in the response to pheromone . One possibility is that one of these genes is the basis of the QTL.
This work has shown that the larval dauer/non-dauer developmental decision and adult reproduction are affected by environmental conditions and that this varies among C. elegans RILs. These analyses have been undertaken on natural variation in these traits, thereby reflecting some aspect of how these traits have evolved and, or are maintained. These results further suggest the hypothesis that the same, or closely linked genes, may co-ordinate the response of these different life-history traits to environmental stress and that there are different evolved strategies by which these traits can be deployed. Further, detailed genetic analyses are required to fully resolve this.
1. C. elegans alters both its propensity to form dauer larvae and its adult rate of reproduction and lifetime fecundity, in response to environmental stress.
2. C. elegans isolates vary in their dauer larva formation and adult reproductive phenotypes and how these vary in response to stress. Changes in dauer larva development and adult reproductive strategy are related, shown by the correlation of these traits among RILs.
3. QTL analyses have identified regions on chromosome II and X that control dauer larva development and population size. The chromosome II QTL was further resolved to genetically separate the control of these two phenotypes.
4. Genetic analysis of the chromosome II QTL affecting dauer larva formation, resolved this to a c. 210 Kbp region containing no genes known to be involved in dauer development.
5. C. elegans has a co-ordinated larval and adult response to environmental stress controlled, at least in part, by physically close genomic regions.
C. elegans isolates N2 and DR1350 were used for the construction of recombinant inbred lines (RILs) and nearly isogenic lines (NILs). All C. elegans lines were maintained on standard NGM plates with an Escherichia coli OP50 food source . Strains were obtained from the Caenorhabditis Genetics Center. RILs were generated as previously described . Briefly, this involved maintaining N2 × DR1350 F1 cross progeny, produced from both N2 hermaphrodite by DR1350 male crosses and reciprocal matings, individually and allowing them to self-fertilise for at least 30 generations, at which point they were cryopreserved. In total 200 RILs were produced. NILs were produced by backcrossing N2 × DR1350 F1 cross-progeny to N2 for nine generations and using polymorphic markers (see below) to genotype the progeny at every generation to identify lines with the desired genotype. All resultant NILs were also genotyped for six markers on other chromosomes to confirm homozygosity for N2 alleles.
Dauer larva formation assays
Dauer pheromone extract was prepared from N2 liquid culture media as previously described . A single batch of pheromone was used for all assays of the RILs. This batch of pheromone was different to that used previously  and so the quantities reported here are not comparable to this previous report . Artificial dauer pheromone, so-called daumone, was synthesized according to the published method . This artificial pheromone, daumone, was used for all assays of the NILs. Subsequently, it has been reported that daumone is one of three molecules that comprise C. elegans dauer pheromone and that daumone induces the lowest level of dauer larva formation, compared with the other components of dauer pheromone . All dauer formation assays were carried out as previously described , with the quantities of pheromone and, or food varied as described below. All dauer assays were performed at 25°C with worms maintained on 3.5 cm diameter plates containing 2 mL of dauer larva formation assay agar .
For the RILs, the dauer larva formation phenotype was measured at both high and low food concentrations against a standard concentration of dauer pheromone (15 μL/mL agar). To prepare the high and low food concentrations, liquid cultures of E. coli OP50 were grown overnight to saturation, centrifuged and the media removed, and resuspended, in water, at 2 and 1% w/v, respectively. RILs were assayed in 10 separate assays, with N2 and DR1350 controls included in each assay. Within each assay, there were three or four replicate plates for each RIL, and between five and ten replicate plates for N2 and DR1350 each at both food concentrations. Plates on which there were fewer than 20 worms and plates on which larvae could not be classified as dauer or non-dauer larvae with confidence, were excluded from all further analyses. Further, when data were only available for a single plate at one combination of conditions, these data were also excluded. The dauer formation phenotype of the NILs was measured at low food concentrations as described above.
To measure the population growth of 48 randomly selected RILs, L1 stages were obtained by allowing eggs liberated from hypochlorite treated gravid hermaphrodites , to develop on NGM plates in the absence of food for 24 hours at 19°C. These arrested L1s were then individually placed on fresh 5.5 cm diameter NGM plates seeded with 100 μL of a stock OP50 food source and maintained at 19°C. The population that grew from these individual L1s were assayed 5, 7 or 9 days later. This was done on day 5 by directly counting the worms (all larval stages and adults) visible on each plate; on days 7 and 9, worms were washed from the plates with M9 buffer, and the number of worms (all larval stages and adults) determined by dilution. Plates on which the worms failed to grow and plates on which the worms had burrowed into the agar were excluded from the analysis. In this manner, the population growth of 48 RILs was measured in 8 experimental blocks, each of which contained 15 replicate plates for 6 RILs. The parental lines, N2 and DR1350, were also included within each experimental block.
The lifetime fecundity and the rate of reproduction, q (see Statistical analyses, below), were determined for 12 RILs chosen to be representative of the range of dauer development phenotypes observed among the 48 RILs whose population growth had been determined. Late-stage L4s were transferred individually to NGM plates that had been previously seeded with an excess of OP50 food source, and the worms were then transferred to fresh plates every 8 hours, until egg laying had ceased. The plates from which the hermaphrodites had been removed were then maintained for 24 hours to allow the eggs to hatch and larvae to develop, which were then counted. This was done in two experimental blocks, with the parental lines, N2 and DR1350, included in each block, with 15 replicate plates for each for line at 19°C.
To determine the effect of different concentrations of food on lifetime fecundity and reproductive timing, these analyses were repeated except that the NGM plates were seeded with a food source consisting of 100 μL of 100, 50, 25, 12.5 or 6.25% w/v OP50 diluted in LB medium. Prior to the addition of these food sources, all plates were supplemented with streptomycin to a final concentration of 50 mg/mL, to prevent the growth of OP50. At low food concentrations, worms leave the site of the food and can be lost from an experiment; to prevent this, heated, sterile, 1 cm diameter brass rings were individually heated and melted into the agar prior to the application of the food and the worms, which were placed, and contained, within the rings . This was done for 6 RILs, which were representative of the range of dauer development phenotypes observed among the 12 RILs analysed above and the parental lines, N2 and DR1350, with 25 replicate plates for each line at each food concentration at 19°C. Plates on which worms had desiccated on the walls of the brass rings or in which worms had burrowed into the agar, were excluded from the analysis.
The lifetime fecundity and reproductive timing of five eat mutants (DA606 eat-10(ad606), DA522 eat-13(ad522), DA573 eat-14(ad573), DA602 eat-15(ad602))  and N2 were determined, as described above, but without the use of brass rings, for individual worms with a food source consisting of 100 μL of 100, 50, 25 or 12.5% w/v OP50 diluted in LB medium. At 100% w/v OP50, 25 replicate plates of each strain were used; 60 replicate plates were used for 50 and 25% w/v OP50 and 70 replicate plates used for 12.5% w/v OP50 at 19°C.
Dauer development and lifetime fecundity were analysed using generalized linear models (GLMs), with LINE and, as appropriate, FOOD and, or PHEROMONE concentration, ASSAY and TIME fitted as factors. When included in analyses, ASSAY was fitted such that LINE was nested within ASSAY. All analyses were performed in Minitab and all proportion data were arcsine transformed prior to analysis. All data presented in the text, figures and tables has been back transformed.
To determine the rate of reproduction, q, the daily cumulative fecundity, expressed as a proportion of lifetime fecundity, was determined for each worm. This gives a relationship where the proportion of eggs laid increases from 0 on day 0 and converges to an asymptote of 1 on the final day of egg laying. The reproductive timing, expressed as the cumulative proportion of eggs laid, can then be described by y = 1 - e(-qx), where q is the rate of reproduction and x is the day of egg laying. The reproductive rate, q, was then estimated individually for each worm by iteratively determining the value that gave the best fit to the data, which was determined by maximum likelihood using negative log likelihoods. A large value of q means that a greater proportion of eggs are laid earlier compared with a low value of q. These rates of reproduction, q, of different lines and under different conditions were analysed by ANOVA.
The lifetime fecundity and reproductive timing of individual worms were also used to calculate the expected population size after 5, 7 or 9 days and these predicted values were correlated to the observed population sizes. Expected population size was calculated using a model in which the lifetime fecundity and reproductive schedule were used to calculate the number of eggs laid in a particular 8 hour period (see additional file 3). The delay between an egg being laid and it, in turn, reproducing was then estimated using data on the timing of embryogenesis and the time from hatch to first egg laying. This process was iterated over the desired period of time, to determine the total population size.
Candidate polymorphic markers that would differentiate N2 and DR1350 were screened by RFLP or direct sequence analysis. Three groups of candidate markers were screened: single nucleotide polymorphisms (SNPs) that differentiated N2 and DR1350 ; SNPs that differentiate N2 and CB4856 [23, 34] and a number of intergenic and intronic regions that were within regions not covered by the two other groups of candidate markers (see additional file 1 for all regions screened). RILs were genotyped either by RFLP analysis or by allele-specific PCR.
Marker order was based on physical position in the N2 genome . Genetic distance between neighbouring markers was calculated from the proportion of recombinants between the markers in the RILs using Kosambi's mapping function . Segregation distortion from an expected 1:1 DR1350: N2 ratio was tested by χ2 analysis for each marker with Bonferroni-corrected significance levels.
QTL mapping was performed using QTL cartographer . In total, nine traits from three experiments were analysed. Dauer development in response to different food concentrations, from which dauer larva formation at high (2% w/v) and at low (1% w/v) food concentrations, and the plasticity of dauer larvae formation, which was calculated as the difference in the proportion of dauer larvae formed between low and high food concentrations, were analysed. Dauer development in response to different pheromone concentrations from previously available data , from which a directly analogous set of traits concerning dauer development in response to different pheromone concentrations were analysed. The population size at day 5, 7 and 9 was also analysed. Each trait was initially investigated by single marker analysis. Data were then analysed by composite interval mapping (CIM), using model 6, forward and backward stepwise regression, a 10 cM window size and a 2 cM walking speed. Genome-wide significance levels, p < 0.05, for each trait were determined by performing 1000 permutations on the data. Effects size and R2 values for significant QTLs were determined from these CIM analyses. To better control for the analysis of multiple traits, each of the three experiments, were also analysed by multiple trait CIM (MT-CIM) [36, 37]. Analysis settings were as described above; genome-wide significance levels, p < 0.05, were determined by performing 1000 permutations on the data. GxE interactions were investigated in two ways. Firstly, the plasticity of dauer larvae formation, calculated as described above, was analysed by CIM as a separate trait. Secondly, GxE interactions were analysed as part of the multiple trait CIM analyses. Epistasis between identified QTLs was assessed by multiple interval mapping (MIM).
James AC, Azevedo RB, Partridge L: Genetic and environmental responses to temperature of Drosophila melanogaster from a latitudinal cline. Genetics. 1997, 146: 881-890.
Müller CB, Williams IS, Hardie J: The role of nutrition, crowding and interspecific interactions in the development of winged aphids. Ecological Entomology. 2001, 26: 330-340. 10.1046/j.1365-2311.2001.00321.x.
Smith-Gill SJ: Developmental plasticity: developmental conversion versus phenotypic modulation. American Zoologist. 1983, 23: 47-55.
Scheiner SM, Lyman RF: The genetics of phenotypic plasticity I: Heritability. Journal of Evolutionary Biology. 1989, 2: 95-107. 10.1046/j.1420-9101.1989.2020095.x.
Scheiner SM: Selection experiments and the study of phenotypic plasticity. Journal of Evolutionary Biology. 2002, 15: 889-898. 10.1046/j.1420-9101.2002.00468.x.
Roff DA: Contributions of genomics to life-history theory. Nature Reviews Genetics. 2007, 8: 116-125. 10.1038/nrg2040.
Li Y, Alvarez OA, Gutteling EW, Tijsterman M, Fu J, Riksen JA, Hazendonk E, Prins P, Plasterk RH, Jansen RC, Breitling R, Kammenga JE: Mapping determinants of gene expression plasticity by genetical genomics in C. elegans. PLoS Genet. 2006, 2: e222-10.1371/journal.pgen.0020222.
Morgan TJ, Mackay TF: Quantitative trait loci for thermotolerance phenotypes in Drosophila melanogaster. Heredity. 2006, 96: 232-242. 10.1038/sj.hdy.6800786.
Gutteling EW, Riksen JA, Bakker J, Kammenga JE: Mapping phenotypic plasticity and genotype-environment interactions affecting life-history traits in Caenorhabditis elegans. Heredity. 2007, 98: 28-37. 10.1038/sj.hdy.6800894.
Mackay TF, Anholt RR: Ain't misbehavin' ? Genotype-environment interactions and the genetics of behavior. Trends in Genetics. 2007, 23: 311-314. 10.1016/j.tig.2007.03.013.
The C. elegans Sequencing Consortium: Genome sequence of the nematode C. elegans: A platform for investigating biology. Science. 1998, 282: 2012-2018. 10.1126/science.282.5396.2012.
Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chinwalla A, Clarke L, Clee C, Coghlan A, Coulson A, D'Eustachio P, Fitch DH, Fulton LA, Fulton RE, Griffiths-Jones S, Harris TW, Hillier LW, Kamath R, Kuwabara PE, Mardis ER, Marra MA, Miner TL, Minx P, Mullikin JC, Plumb RW, Rogers J, Schein JE, Sohrmann M, Spieth J, Stajich JE, Wei C, Willey D, Wilson RK, Durbin R, Waterston RH: The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics. PLoS Biology. 2003, 1: E45-10.1371/journal.pbio.0000045.
Pigliucci M: Evolution of phenotypic plasticity: where are we going now?. Trends in Ecology and Evolution. 2005, 20: 481-486. 10.1016/j.tree.2005.06.001.
Riddle DLA: Genetic and environmental regulation of dauer larvae development. C elegans II. Edited by: Riddle DS, Blumenthal T, Meyer BJ, Priess JR. 1997, NY, Cold Spring Harbour Laboratory Press., 739-768.
Hu PJ: Dauer. Wormbook. Edited by: Community TCR. 2007, Wormbook, [http://www.wormbook.org]
Golden JW, Riddle DL: A pheromone influences larval development in the nematode Caenorhabditis elegans. Science. 1982, 218: 578-580. 10.1126/science.6896933.
Jeong PY, Jung M, Yim YH, Kim H, Park M, Hong E, Lee W, Kim YH, Kim K, Paik YK: Chemical structure and biological activity of the Caenorhabditis elegans dauer-inducing pheromone. Nature. 2005, 433: 541-545. 10.1038/nature03201.
Butcher RA, Fujita M, Schroeder FC, Clardy J: Small-molecule pheromones that control dauer development in Caenorhabditis elegans. Nature Chemical Biology. 2007, 3: 420-422. 10.1038/nchembio.2007.3.
Viney ME, Gardner MP, Jackson JA: Variation in Caenorhabditis elegans dauer larva formation. Development Growth and Differentiation. 2003, 45: 389-396. 10.1046/j.1440-169X.2003.00703.x.
Klass MR: Aging in the nematode Caenorhabditis elegans: major biological and environmental factors influencing the lifespan. Mechanisms of Aging and Development. 1977, 6: 413-429. 10.1016/0047-6374(77)90043-4.
Avery L: The genetics of feeding in Caenorhabditis elegans. Genetics. 1993, 133: 897-917.
Koch R, van Luenen HG, van der Horst M, Thijssen KL, Plasterk RH: Single nucleotide polymorphisms in wild isolates of Caenorhabditis elegans. Genome Research. 2000, 10: 1690-1696. 10.1101/gr.GR-1471R.
WormBase: Wormbase. [http://www.wormbase.org]
Ebert RH, Cherkasova VA, Dennis RA, Wu JH, Ruggles S, Perrin TE, Shmookler Reis RJ: Longevity-determining genes in Caenorhabditis elegans: chromosomal mapping of multiple noninteractive loci. Genetics. 1993, 135: 1003-1010.
Shook DR, Brooks A, Johnson TE: Mapping quantitative trait loci affecting life history traits in the nematode Caenorhabditis elegans. Genetics. 1996, 142: 801-817.
Barrière A, Félix MA: Natural variation and population genetics of Caenorhabditis elegans. WormBook. Edited by: Community TCR. 2005, WormBook, [http://www.wormbook.org]
Decuyper C, Vanfleteran JR: Nutritional Alteration of Life-Span in the Nematode Caenorhabditis elegans. Age. 1982, 5: 42-45. 10.1007/BF02431722.
Houthoofd K, Braeckman BP, Lenaerts I, Brys K, De Vreese A, Van Eygen S, Vanfleteren JR: Axenic growth up-regulates mass-specific metabolic rate, stress resistance, and extends life span in Caenorhabditis elegans. Experimental Gerontology. 2002, 37: 1371-1378. 10.1016/S0531-5565(02)00173-0.
Szewczyk NJ, Udranszky IA, Kozak E, Sunga J, Kim SK, Jacobson LA, Conley CA: Delayed development and lifespan extension as features of metabolic lifestyle alteration in C. elegans under dietary restriction. Journal of Experimental Biology. 2006, 209: 4129-4139. 10.1242/jeb.02492.
Zwaal RR, Mendel JE, Sternberg PW, Plasterk RH: Two neuronal G proteins are involved in chemosensation of the Caenorhabditis elegans Dauer-inducing pheromone. Genetics. 1997, 145: 715-727.
Stiernagle T: Methods. C elegans: a practical approach. Edited by: Hope IA. 1999, , Oxford University Press, Oxford, 51-67.
Golden JW, Riddle DL: The Caenorhabditis elegans dauer larva: developmental effects of pheromone, food, and temperature. Developmental Biology. 1984, 102: 368-378. 10.1016/0012-1606(84)90201-X.
Davies AG, McIntire SL: Using C. elegans to screen for targets of ethanol and behavior-altering drugs. Biological Procedures Online. 2004, 6: 113-119. 10.1251/bpo79.
Wicks SR, Yeh RT, Gish WR, Waterston RH, Plasterk RH: Rapid gene mapping in Caenorhabditis elegans using a high density polymorphism map. Nature Genetics. 2001, 28: 160-164. 10.1038/88878.
Kosambi DD: The estimation of map distances from recombination values. Annals of Eugenics. 1944, 12: 172-175.
Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. 2007, [http://statgen.ncsu.edu/qtlcart/WQTLCart.htm]
Jiang C, Zeng ZB: Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics. 1995, 140: 1111-1127.
We thank Sue Coe for technical assistance and Mark Steer for helpful discussions. Daumone was synthesised by Philip Wall and Guy Lloyd-Jones. C. elegans isolates were supplied by the Caenorhabditis Genetics Center which is supported by the National Institutes of Health, National Center for Research Resources, USA. This work was supported by a NERC Environmental Genomics thematic programme grant.
MEV conceived and directed the study. SCH devised and executed the dauer larva phenotyping, genotyping, all molecular genetic and genetic analyses and data analyses; AS undertook the population growth and fecundity assays. The manuscript was written by SCH and MEV. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Genomic regions and markers screened in DR1350. (A) Genomic regions sequenced in DR1350, showing the physical location, the size of the amplified region and the PCR oligonucleotide primers (5'-3') and, where appropriate, the candidate SNP contained within the amplified region . (B) N2/CB4856 RFLP polymorphisms screened in DR1350, showing the polymorphic marker , the chromosome and physical location (Mbp), the restriction endonuclease used and the PCR oligonucleotide primers (5'-3'). No N2/DR1350 polymorphisms were detected. (XLS 39 KB)
Additional file 2: Candidate genes underlying the chromosome II QTLs. Genes and pseudogenes predicted to occur in N2 within the defined regions of the chromosome II QTLs for (A) dauer larva development and (B) for population size on day 5. For each gene, the Wormbase ID number, gene name, gene description and position and, where available, any RNAi phenotypes, the expression pattern and gene ontology (GO) terms are given . 7TM, 7-transmembrane; Adl, adult lethal; Bli, blistered, Bmb, organism morphology abnormal; Ced, cell death abnormal; Clr, clear; Dpy, dumpy; Egl, egg laying abnormal; Emb, embryonic lethal; Gro, slow growth; Let, lethal; Lva, larval arrest; Lvl, larval lethal; Muv, multivulva; Pch, patchy colouration; Prl, paralysed; Pvl, protruding vulva, Ric, aldicarb resistant; Rup, exploded through vulva; Sck, sick; Sma, small; Ste, maternal sterile; Stp, sterile progeny; Unc, locomotion abnormal. (XLS 62 KB)
Additional file 3: Model used in the calculation of predicted population size. The equation used in the calculation of the predicted population size and the definition of the parameters. (DOC 24 KB)
Authors’ original submitted files for images
About this article
Cite this article
Harvey, S.C., Shorto, A. & Viney, M.E. Quantitative genetic analysis of life-history traits of Caenorhabditis elegans in stressful environments. BMC Evol Biol 8, 15 (2008). https://doi.org/10.1186/1471-2148-8-15