Environmental influence on the genetic correlations between life-history traits in Caenorhabditis elegans

Rapid evolutionary change in trait correlations of single proteins

Many organismal traits are genetically determined and covary in evolving populations. The resulting trait correlations can either help or hinder evolvability - the ability to bring forth new and adaptive phenotypes. The evolution of evolvability requires that trait correlations themselves must be able to evolve, but we know little about this ability. To learn more about it, we here study two evolvable systems, a yellow fluorescent protein and the antibiotic resistance protein VIM-2 metallo beta-lactamase. We consider two traits in the fluorescent protein, namely the ability to emit yellow and green light, and three traits in our enzyme, namely the resistance against ampicillin, cefotaxime, and meropenem. We show that correlations between these traits can evolve rapidly through both mutation and selection on short evolutionary time scales. In addition, we show that these correlations are driven by a protein's ability to fold, because single mutations that alter foldability can dramatically change trait correlations. Since foldability is important for most proteins and their traits, mutations affecting protein folding may alter trait correlations mediated by many other proteins. Thus, mutations that affect protein foldability may also help shape the correlations of complex traits that are affected by hundreds of proteins.

Microbes are potential key players in the evolution of life histories and aging in Caenorhabditis elegans

Microbes can have profound effects on host fitness and health and the appearance of late-onset diseases. Host-microbe interactions thus represent a major environmental context for healthy aging of the host and might also mediate trade-offs between life-history traits in the evolution of host senescence. Here, we have used the nematode Caenorhabditis elegans to study how host-microbe interactions may modulate the evolution of life histories and aging. We first characterized the effects of two non-pathogenic and one pathogenic Escherichia coli strains, together with the pathogenic Serratia marcescens DB11 strain, on population growth rates and survival of C. elegans from five different genetic backgrounds. We then focused on an outbred C. elegans population, to understand if microbe-specific effects on the reproductive schedule and in traits such as developmental rate and survival were also expressed in the presence of males and standing genetic variation, which could be relevant for the evolution of C. elegans and other nematode species in nature. Our results show that host-microbe interactions have a substantial host-genotype-dependent impact on the reproductive aging and survival of the nematode host. Although both pathogenic bacteria reduced host survival in comparison with benign strains, they differed in how they affected other host traits. Host fertility and population growth rate were affected by S. marcescens DB11 only during early adulthood, whereas this occurred at later ages with the pathogenic E. coli IAI1. In both cases, these effects were largely dependent on the host genotypes. Given such microbe-specific genotypic differences in host life history, we predict that the evolution of reproductive schedules and senescence might be critically contingent on host-microbe interactions in nature.

The genetic architecture underlying body-size traits plasticity over different temperatures and developmental stages in Caenorhabditis elegans

Most ectotherms obey the temperature-size rule, meaning they grow larger in a colder environment. This raises the question of how the interplay between genes and temperature affects the body size of ectotherms. Despite the growing body of literature on the physiological life-history and molecular genetic mechanism underlying the temperature-size rule, the overall genetic architecture orchestrating this complex phenotype is not yet fully understood. One approach to identify genetic regulators of complex phenotypes is quantitative trait locus (QTL) mapping. Here, we explore the genetic architecture of body-size phenotypes, and plasticity of body-size phenotypes at different temperatures using Caenorhabditis elegans as a model ectotherm. We used 40 recombinant inbred lines (RILs) derived from N2 and CB4856, which were reared at four different temperatures (16, 20, 24, and 26 °C) and measured at two developmental stages (L4 and adult). The animals were measured for body length, width at vulva, body volume, length/width ratio, and seven other body-size traits. The genetically diverse RILs varied in their body-size phenotypes with heritabilities ranging from 0.0 to 0.99. We detected 18 QTL underlying the body-size traits across all treatment combinations, with the majority clustering on Chromosome X. We hypothesize that the Chromosome X QTL could result from a known pleiotropic regulator-npr-1-known to affect the body size of C. elegans through behavioral changes. We also found five plasticity QTL of body-size traits where three colocalized with body-size QTL. In conclusion, our findings shed more light on multiple loci affecting body-size plasticity and the possibility of co-regulation of traits and traits plasticity by the same loci under different environments.

Genetic coupling of life-history and aerobic performance in Atlantic salmon

A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions (vgll3 and six6) makes Atlantic salmon (Salmo salar) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes-standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS)-and the life-history genomic regions, and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Furthermore, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve.

From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans

Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.

Unidirectional response to bidirectional selection on body size II. Quantitative genetics

Anticipating the genetic and phenotypic changes induced by natural or artificial selection requires reliable estimates of trait evolvabilities (genetic variances and covariances). However, whether or not multivariate quantitative genetics models are able to predict precisely the evolution of traits of interest, especially fitness-related, life history traits, remains an open empirical question. Here, we assessed to what extent the response to bivariate artificial selection on both body size and maturity in the medaka Oryzias latipes, a model fish species, fits the theoretical predictions. Three lines (Large, Small, and Control lines) were differentially selected for body length at 75 days of age, conditional on maturity. As maturity and body size were phenotypically correlated, this selection procedure generated a bi-dimensional selection pattern on two life history traits. After removal of nonheritable trends and noise with a random effect ("animal") model, the observed selection response did not match the expected bidirectional response. For body size, Large and Control lines responded along selection gradients (larger body size and stasis, respectively), but, surprisingly, the Small did not evolve a smaller body length and remained identical to the Control line throughout the experiment. The magnitude of the empirical response was smaller than the theoretical prediction in both selected directions. For maturity, the response was opposite to the expectation (the Large line evolved late maturity compared to the Control line, while the Small line evolved early maturity, while the opposite pattern was predicted due to the strong positive genetic correlation between both traits). The mismatch between predicted and observed response was substantial and could not be explained by usual sources of uncertainties (including sampling effects, genetic drift, and error in G matrix estimates).

The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

Pleiotropy, the concept that a single gene controls multiple distinct traits, is prevalent in most organisms and has broad implications for medicine and agriculture. The identification of the molecular mechanisms underlying pleiotropy has the power to reveal previously unknown biological connections between seemingly unrelated traits. Additionally, the discovery of pleiotropic genes increases our understanding of both genetic and phenotypic complexity by characterizing novel gene functions. Quantitative trait locus (QTL) mapping has been used to identify several pleiotropic regions in many organisms. However, gene knockout studies are needed to eliminate the possibility of tightly linked, non-pleiotropic loci. Here, we use a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to identify a single large-effect QTL on the center of chromosome V associated with variation in responses to eight chemotherapeutics. We validate this QTL with near-isogenic lines and pair genome-wide gene expression data with drug response traits to perform mediation analysis, leading to the identification of a pleiotropic candidate gene, scb-1, for some of the eight chemotherapeutics. Using deletion strains created by genome editing, we show that scb-1, which was previously implicated in response to bleomycin, also underlies responses to other double-strand DNA break-inducing chemotherapeutics. This finding provides new evidence for the role of scb-1 in the nematode drug response and highlights the power of mediation analysis to identify causal genes.

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

Quantitative genetics provides the tools for linking polymorphic loci to trait variation. Linkage analysis of gene expression is an established and widely applied method, leading to the identification of expression quantitative trait loci (eQTLs). (e)QTL detection facilitates the identification and understanding of the underlying molecular components and pathways, yet (e)QTL data access and mining often is a bottleneck. Here, we present WormQTL2, a database and platform for comparative investigations and meta-analyses of published (e)QTL data sets in the model nematode worm C. elegans. WormQTL2 integrates six eQTL studies spanning 11 conditions as well as over 1000 traits from 32 studies and allows experimental results to be compared, reused and extended upon to guide further experiments and conduct systems-genetic analyses. For example, one can easily screen a locus for specific cis-eQTLs that could be linked to variation in other traits, detect gene-by-environment interactions by comparing eQTLs under different conditions, or find correlations between QTL profiles of classical traits and gene expression. WormQTL2 makes data on natural variation in C. elegans and the identified QTLs interactively accessible, allowing studies beyond the original publications. Database URL: www.bioinformatics.nl/WormQTL2/.

Carryover effects of larval environment on individual variation in a facultatively diadromous fish

Intraspecific trait variation may result from "carryover effects" of variability of environments experienced at an earlier life stage. This phenomenon is particularly relevant in partially migrating populations composed of individuals with divergent early life histories. While many studies have addressed the causes of partial migration, few have investigated the consequences for between-individual variability later in life.We studied carryover effects of larval environment in a facultatively diadromous New Zealand fish, Gobiomorphus cotidianus, along an estuarine salinity gradient. We investigated the implications of varying environmental conditions during this critical stage of ontogeny for adult phenotype.We inferred past environmental history of wild-caught adult fish using otolith microchemistry (Sr/Ca) as a proxy for salinity. We tested for main and interactive effects of larval and adult environment on a suite of traits, including growth rates, behavior (exploration and activity), parasite load, and diet (stable isotopes and gut contents).We found a Sr/Ca consistent with a continuum from freshwater to brackish environments, and with different trajectories from juvenile to adult habitat. Fish with Sr/Ca indicating upstream migration were more vulnerable to trematode infection, suggesting a mismatch to freshwater habitat. Diet analysis suggested an interactive effect of larval and adult environments on trophic position and diet preference, while behavioral traits were unrelated to environment at any life stage. Growth rates did not seem to be affected by past environment.Overall, we show that early life environment can have multiple effects on adult performance and ecology, with the potential for lifetime fitness trade-offs associated with life history. Our study highlights that even relatively minor variation in rearing conditions may be enough to generate individual variation in natural populations.

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits

BACKGROUND

The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale.

RESULTS

To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping.

CONCLUSION

Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.

Genetic parameters for faecal worm egg count at different ages in Australian sheep under natural challenge

The data used in the present study consisted of 24535 worm egg count records on sheep observed from 63 to 560 days of age under conditions of the natural challenge of trichostrongylid species. Records were extracted from the Information Nucleus Flock database of the Australia Sheep Cooperative Research Centre program from 2007 to 2011. Records were observed at various ages and subdivided into weaning (W, ~3 months), post-weaning (P, ~4 months), yearling (Y, ~12 months) and hogget (H, ~18 months) age stages and were used to investigate genetic variation at different age stages in univariate analyses and estimate genetic correlations between age stages in multi-trait analyses. The full data were also analysed by random regression models to study how heritability and genetic correlations varied with age. Heritability estimates from univariate analyses were 0.20 ± 0.05, 0.15 ± 0.02, 0.36 ± 0.09, 0.22 ± 0.06 for W, P, Y and H age stages respectively. A similar trend of heritability over ages was found from random regression analyses, which decreased from 0.16 at 90 days to 0.09 at 120 days, following a steady increase to 0.32 at ~410 days, and then decreased afterwards to 0.24 at 520 days. Strong genetic correlations (>0.8) were found between W and P age stages, along with Y and H age stages. Sire by flock interaction effects were significant, and accounted for the reduced estimates of heritability and increased genetic correlations between age stages. The results indicated that a multiple-trait approach is required for genetic evaluation of worm egg count when measurements are at different ages, and the accuracy of evaluations would benefit from recording at least two separate age stages.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation.

Polygenicity and Epistasis Underlie Fitness-Proximal Traits in the Caenorhabditis elegans Multiparental Experimental Evolution (CeMEE) Panel

Understanding the genetic basis of complex traits remains a major challenge in biology. Polygenicity, phenotypic plasticity, and epistasis contribute to phenotypic variance in ways that are rarely clear. This uncertainty can be problematic for estimating heritability, for predicting individual phenotypes from genomic data, and for parameterizing models of phenotypic evolution. Here, we report an advanced recombinant inbred line (RIL) quantitative trait locus mapping panel for the hermaphroditic nematode Caenorhabditis elegans, the C. elegans multiparental experimental evolution (CeMEE) panel. The CeMEE panel, comprising 507 RILs at present, was created by hybridization of 16 wild isolates, experimental evolution for 140-190 generations, and inbreeding by selfing for 13-16 generations. The panel contains 22% of single-nucleotide polymorphisms known to segregate in natural populations, and complements existing C. elegans mapping resources by providing fine resolution and high nucleotide diversity across > 95% of the genome. We apply it to study the genetic basis of two fitness components, fertility and hermaphrodite body size at time of reproduction, with high broad-sense heritability in the CeMEE. While simulations show that we should detect common alleles with additive effects as small as 5%, at gene-level resolution, the genetic architectures of these traits do not feature such alleles. We instead find that a significant fraction of trait variance, approaching 40% for fertility, can be explained by sign epistasis with main effects below the detection limit. In congruence, phenotype prediction from genomic similarity, while generally poor ([Formula: see text]), requires modeling epistasis for optimal accuracy, with most variance attributed to the rapidly evolving chromosome arms.

Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans

Background

Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.

Results

We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.

Conclusions

These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3899-8) contains supplementary material, which is available to authorized users.

Correlations of Genotype with Climate Parameters Suggest Caenorhabditis elegans Niche Adaptations

Species inhabit a variety of environmental niches, and the adaptation to a particular niche is often controlled by genetic factors, including gene-by-environment interactions. The genes that vary in order to regulate the ability to colonize a niche are often difficult to identify, especially in the context of complex ecological systems and in experimentally uncontrolled natural environments. Quantitative genetic approaches provide an opportunity to investigate correlations between genetic factors and environmental parameters that might define a niche. Previously, we have shown how a collection of 208 whole-genome sequenced wild Caenorhabditis elegans can facilitate association mapping approaches. To correlate climate parameters with the variation found in this collection of wild strains, we used geographic data to exhaustively curate daily weather measurements in short-term (3 month), middle-term (one year), and long-term (three year) durations surrounding the date of strain isolation. These climate parameters were used as quantitative traits in association mapping approaches, where we identified 11 quantitative trait loci (QTL) for three climatic variables: elevation, relative humidity, and average temperature. We then narrowed the genomic interval of interest to identify gene candidates with variants potentially underlying phenotypic differences. Additionally, we performed two-strain competition assays at high and low temperatures to validate a QTL that could underlie adaptation to temperature and found suggestive evidence supporting that hypothesis.

Contrasting invertebrate immune defense behaviors caused by a single gene, the Caenorhabditis elegans neuropeptide receptor gene npr-1

Background

The invertebrate immune system comprises physiological mechanisms, physical barriers and also behavioral responses. It is generally related to the vertebrate innate immune system and widely believed to provide nonspecific defense against pathogens, whereby the response to different pathogen types is usually mediated by distinct signalling cascades. Recent work suggests that invertebrate immune defense can be more specific at least at the phenotypic level. The underlying genetic mechanisms are as yet poorly understood.

Results

We demonstrate in the model invertebrate Caenorhabditis elegans that a single gene, a homolog of the mammalian neuropeptide Y receptor gene, npr-1, mediates contrasting defense phenotypes towards two distinct pathogens, the Gram-positive Bacillus thuringiensis and the Gram-negative Pseudomonas aeruginosa. Our findings are based on combining quantitative trait loci (QTLs) analysis with functional genetic analysis and RNAseq-based transcriptomics. The QTL analysis focused on behavioral immune defense against B. thuringiensis, using recombinant inbred lines (RILs) and introgression lines (ILs). It revealed several defense QTLs, including one on chromosome X comprising the npr-1 gene. The wildtype N2 allele for the latter QTL was associated with reduced defense against B. thuringiensis and thus produced an opposite phenotype to that previously reported for the N2 npr-1 allele against P. aeruginosa. Analysis of npr-1 mutants confirmed these contrasting immune phenotypes for both avoidance behavior and nematode survival. Subsequent transcriptional profiling of C. elegans wildtype and npr-1 mutant suggested that npr-1 mediates defense against both pathogens through p38 MAPK signaling, insulin-like signaling, and C-type lectins. Importantly, increased defense towards P. aeruginosa seems to be additionally influenced through the induction of oxidative stress genes and activation of GATA transcription factors, while the repression of oxidative stress genes combined with activation of Ebox transcription factors appears to enhance susceptibility to B. thuringiensis.

Conclusions

Our findings highlight the role of a single gene, npr-1, in fine-tuning nematode immune defense, showing the ability of the invertebrate immune system to produce highly specialized and potentially opposing immune responses via single regulatory genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2603-8) contains supplementary material, which is available to authorized users.

The laboratory domestication of Caenorhabditis elegans

Model organisms are of great importance to our understanding of basic biology and to making advances in biomedical research. However, the influence of laboratory cultivation on these organisms is underappreciated, and especially how that environment can affect research outcomes. Recent experiments led to insights into how the widely used laboratory reference strain of the nematode Caenorhabditis elegans compares with natural strains. Here we describe potential selective pressures that led to the fixation of laboratory-derived alleles for the genes npr-1, glb-5, and nath-10. These alleles influence a large number of traits, resulting in behaviors that affect experimental interpretations. Furthermore, strong phenotypic effects caused by these laboratory-derived alleles hinder the discovery of natural alleles. We highlight strategies to reduce the influence of laboratory-derived alleles and to harness the full power of C. elegans.

Pollution breaks down the genetic architecture of life history traits in Caenorhabditis elegans

When pollution occurs in an environment, populations present suffer numerous negative and immediate effects on their life history traits. Their evolutionary potential to live in a highly stressful environment will depend on the selection pressure strengths and on the genetic structure, the trait heritability, and the genetic correlations between them. If expression of this structure changes in a stressful environment, it becomes necessary to quantify these changes to estimate the evolutionary potential of the population in this new environment. We studied the genetic structure for survival, fecundity, and early and late growth in isogenic lines of a Caenorhabditis elegans population subject to three different environments: a control environment, an environment polluted with uranium, and a high salt concentration environment. We found a heritability decrease in the polluted environments for fecundity and early growth, two traits that were the most heritable in the control environment. The genetic structure of the traits was particularly affected in the uranium polluted environment, probably due to generally low heritability in this environment. This could prevent selection from acting on traits despite the strong selection pressures exerted on them. Moreover, phenotypic traits were more strongly affected in the salt than in the uranium environment and the heritabilities were also lower in the latter environment. Consequently the decrease in heritability was not proportional to the population fitness reduction in the polluted environments. Our results suggest that pollution can alter the genetic structure of a C. elegans population, and thus modify its evolutionary potential.

A streamlined system for species diagnosis in Caenorhabditis (Nematoda: Rhabditidae) with name designations for 15 distinct biological species

The rapid pace of species discovery outstrips the rate of species description in many taxa. This problem is especially acute for Caenorhabditis nematodes, where the naming of distinct species would greatly improve their visibility and usage for biological research, given the thousands of scientists studying Caenorhabditis. Species description and naming has been hampered in Caenorhabditis, in part due to the presence of morphologically cryptic species despite complete biological reproductive isolation and often enormous molecular divergence. With the aim of expediting species designations, here we propose and apply a revised framework for species diagnosis and description in this group. Our solution prioritizes reproductive isolation over traditional morphological characters as the key feature in delineating and diagnosing new species, reflecting both practical considerations and conceptual justifications. DNA sequence divergence criteria help prioritize crosses for establishing patterns of reproductive isolation among the many species of Caenorhabditis known to science, such as with the ribosomal internal transcribed spacer-2 (ITS2) DNA barcode. By adopting this approach, we provide new species name designations for 15 distinct biological species, thus increasing the number of named Caenorhabditis species in laboratory culture by nearly 3-fold. We anticipate that the improved accessibility of these species to the research community will expand the opportunities for study and accelerate our understanding of diverse biological phenomena.

Worm variation made accessible: Take your shopping cart to store, link, and investigate!

In Caenorhabditis elegans, the recent advances in high-throughput quantitative analyses of natural genetic and phenotypic variation have led to a wealth of data on genotype phenotype relations. This data has resulted in the discovery of genes with major allelic effects and insights in the effect of natural genetic variation on a whole range of complex traits as well as how this variation is distributed across the genome. Regardless of the advances presented in specific studies, the majority of the data generated in these studies had yet to be made easily accessible, allowing for meta-analysis. Not only data in figures or tables but meta-data should be accessible for further investigation and comparison between studies. A platform was created where all the data, phenotypic measurements, genotypes, and mappings can be stored, compared, and new linkages within and between published studies can be discovered. WormQTL focuses on quantitative genetics in Caenorhabditis and other nematode species, whereas WormQTL^HD quantitatively links gene expression quantitative trait loci (eQTL) in C. elegans to gene–disease associations in humans.

Quantitative trait loci × environment interactions for plant morphology vary over ontogeny in Brassica rapa

Growth in plants occurs via the addition of repeating modules, suggesting that the genetic architecture of similar subunits may vary between earlier- and later-developing modules. These complex environment × ontogeny interactions are not well elucidated, as studies examining quantitative trait loci (QTLs) expression over ontogeny have not included multiple environments. Here, we characterized the genetic architecture of vegetative traits and onset of reproduction over ontogeny in recombinant inbred lines of Brassica rapa in the field and glasshouse. The magnitude of genetic variation in plasticity of seedling internodes was greater than in those produced later in ontogeny. We correspondingly detected that QTLs for seedling internode length were environment-specific, whereas later in ontogeny the majority of QTLs affected internode lengths in all treatments. The relationship between internode traits and onset of reproduction varied with environment and ontogenetic stage. This relationship was observed only in the glasshouse environment and was largely attributable to one environment-specific QTL. Our results provide the first evidence of a QTL × environment × ontogeny interaction, and provide QTL resolution for differences between early- and later-stage plasticity for stem elongation. These results also suggest potential constraints on morphological evolution in early vs later modules as a result of associations with reproductive timing.

WormQTLHD--a web database for linking human disease to natural variation data in C. elegans

Interactions between proteins are highly conserved across species. As a result, the molecular basis of multiple diseases affecting humans can be studied in model organisms that offer many alternative experimental opportunities. One such organism—Caenorhabditis elegans—has been used to produce much molecular quantitative genetics and systems biology data over the past decade. We present WormQTL^HD (Human Disease), a database that quantitatively and systematically links expression Quantitative Trait Loci (eQTL) findings in C. elegans to gene–disease associations in man. WormQTL^HD, available online at http://www.wormqtl-hd.org, is a user-friendly set of tools to reveal functionally coherent, evolutionary conserved gene networks. These can be used to predict novel gene-to-gene associations and the functions of genes underlying the disease of interest. We created a new database that links C. elegans eQTL data sets to human diseases (34 337 gene–disease associations from OMIM, DGA, GWAS Central and NHGRI GWAS Catalogue) based on overlapping sets of orthologous genes associated to phenotypes in these two species. We utilized QTL results, high-throughput molecular phenotypes, classical phenotypes and genotype data covering different developmental stages and environments from WormQTL database. All software is available as open source, built on MOLGENIS and xQTL workbench.

Gene-environment and protein-degradation signatures characterize genomic and phenotypic diversity in wild Caenorhabditis elegans populations

Background

Analyzing and understanding the relationship between genotypes and phenotypes is at the heart of genetics. Research on the nematode Caenorhabditis elegans has been instrumental for unraveling genotype-phenotype relations, and has important implications for understanding the biology of mammals, but almost all studies, including forward and reverse genetic screens, are limited by investigations in only one canonical genotype. This hampers the detection and functional analysis of allelic variants, which play a key role in controlling many complex traits. It is therefore essential to explore the full potential of the natural genetic variation and evolutionary context of the genotype-phenotype map in wild C. elegans populations.

Results

We used multiple wild C. elegans populations freshly isolated from local sites to investigate gene sequence polymorphisms and a multitude of phenotypes including the transcriptome, fitness, and behavioral traits. The genotype, transcriptome, and a number of fitness traits showed a direct link with the original site of the strains. The separation between the isolation sites was prevalent on all chromosomes, but chromosome V was the largest contributor to this variation. These results were supported by a differential food preference of the wild isolates for naturally co-existing bacterial species. Comparing polymorphic genes between the populations with a set of genes extracted from 19 different studies on gene expression in C. elegans exposed to biotic and abiotic factors, such as bacteria, osmotic pressure, and temperature, revealed a significant enrichment for genes involved in gene-environment interactions and protein degradation.

Conclusions

We found that wild C. elegans populations are characterized by gene-environment signatures, and we have unlocked a wealth of genotype-phenotype relations for the first time. Studying natural isolates provides a treasure trove of evidence compared with that unearthed by the current research in C. elegans, which covers only a diminutive part of the myriad of genotype-phenotype relations that are present in the wild.

Genotypic and phenotypic variation in transmission traits of a complex life cycle parasite

Characterizing genetic variation in parasite transmission traits and its contribution to parasite vigor is essential for understanding the evolution of parasite life-history traits. We measured genetic variation in output, activity, survival, and infection success of clonal transmission stages (cercaria larvae) of a complex life cycle parasite (Diplostomum pseudospathaceum). We further tested if variation in host nutritional stage had an effect on these traits by keeping hosts on limited or ad libitum diet. The traits we measured were highly variable among parasite genotypes indicating significant genetic variation in these life-history traits. Traits were also phenotypically variable, for example, there was significant variation in the measured traits over time within each genotype. However, host nutritional stage had no effect on the parasite traits suggesting that a short-term reduction in host resources was not limiting the cercarial output or performance. Overall, these results suggest significant interclonal and phenotypic variation in parasite transmission traits that are not affected by host nutritional status.

Worms under stress: C. elegans stress response and its relevance to complex human disease and aging

Many organisms have stress response pathways, components of which share homology with players in complex human disease pathways. Research on stress response in the nematode worm Caenorhabditis elegans has provided detailed insights into the genetic and molecular mechanisms underlying complex human diseases. In this review we focus on four different types of environmental stress responses - heat shock, oxidative stress, hypoxia, and osmotic stress - and on how these can be used to study the genetics of complex human diseases. All four types of responses involve the genetic machineries that underlie a number of complex human diseases such as cancer and neurodegenerative diseases, including Alzheimer's and Parkinson's. We highlight the types of stress response experiments required to detect the genes and pathways underlying human disease and suggest that studying stress biology in worms can be translated to understanding human disease and provide potential targets for drug discovery.

WormQTL--public archive and analysis web portal for natural variation data in Caenorhabditis spp

Here, we present WormQTL (http://www.wormqtl.org), an easily accessible database enabling search, comparative analysis and meta-analysis of all data on variation in Caenorhabditis spp. Over the past decade, Caenorhabditis elegans has become instrumental for molecular quantitative genetics and the systems biology of natural variation. These efforts have resulted in a valuable amount of phenotypic, high-throughput molecular and genotypic data across different developmental worm stages and environments in hundreds of C. elegans strains. WormQTL provides a workbench of analysis tools for genotype-phenotype linkage and association mapping based on but not limited to R/qtl (http://www.rqtl.org). All data can be uploaded and downloaded using simple delimited text or Excel formats and are accessible via a public web user interface for biologists and R statistic and web service interfaces for bioinformaticians, based on open source MOLGENIS and xQTL workbench software. WormQTL welcomes data submissions from other worm researchers.

The genetics of primary and secondary sexual character trade-offs in a horned beetle

When structures compete for shared resources, this may lead to acquisition and allocation trade-offs so that the enlargement of one structure occurs at the expense of another. Among the studies of morphological trade-offs, their importance has been demonstrated primarily through experimental manipulations and comparative analyses. Relatively, a few studies have investigated the underlying genetic basis of phenotypic patterns. Here, we use a half-sibling breeding design to determine the genetic underpinnings of the phenotypic trade-off between head horns and the male copulatory organ or aedeagus that has been found in the dung beetle Onthophagus taurus. Instead of the predicted negative genetic covariance among characters that trade-off, we find positive genetic covariance between absolute horn and aedeagus length and zero genetic covariance between relative horn and aedeagus length. Therefore, although the genetic covariance between absolute horn and aedeagus length would constrain the independent evolution of primary and secondary sexual characters in this population, there was no evidence of a trade-off. We discuss alternative hypotheses for the observed patterns of genetic correlation between traits that compete for resources and the implications that these have for selection and the evolution of such traits.

Genetic variation for stress-response hormesis in C. elegans lifespan

Increased lifespan can be associated with greater resistance to many different stressors, most notably thermal stress. Such hormetic effects have also been found in C. elegans where short-term exposure to heat lengthens the lifespan. Genetic investigations have been carried out using mutation perturbations in a single genotype, the wild type Bristol N2. Yet, induced mutations do not yield insight regarding the natural genetic variation of thermal tolerance and lifespan. We investigated the genetic variation of heat-shock recovery, i.e. hormetic effects on lifespan and associated quantitative trait loci (QTL) in C. elegans. Heat-shock resulted in an 18% lifespan increase in wild type CB4856 whereas N2 did not show a lifespan elongation. Using recombinant inbred lines (RILs) derived from a cross between wild types N2 and CB4856 we found natural variation in stress-response hormesis in lifespan. Approx. 28% of the RILs displayed a hormesis effect in lifespan. We did not find any hormesis effects for total offspring. Across the RILs there was no relation between lifespan and offspring. The ability to recover from heat-shock mapped to a significant QTL on chromosome II which overlapped with a QTL for offspring under heat-shock conditions. The QTL was confirmed by introgressing relatively small CB4856 regions into chromosome II of N2. Our observations show that there is natural variation in hormetic effects on C. elegans lifespan for heat-shock and that this variation is genetically determined.

Aging Uncouples Heritability and Expression-QTL in Caenorhabditis elegans

The number and distribution of gene expression QTL (eQTL) represent the genetic architecture of many complex traits, including common human diseases. We previously reported that the heritable eQTL patterns are highly dynamic with age in an N2 × CB4856 recombinant inbred population of the nematode Caenorhabditis elegans. In particular, we showed that the number of eQTL decreased with age. Here, we investigated the reason for this decrease by combining gene expression profiles at three ages in the wild types N2 and CB4856 with the reported expression profiles of the RIL population. We determined heritability and transgression (when gene expression levels in the RILs are more extreme than the parents) and investigated their relation with eQTL changes with age. Transgressive segregation was widespread but depended on physiological age. The percentage of genes with an eQTL increased with a higher heritability in young worms. However, for old worms this percentage hardly increased. Using a single marker approach, we found that almost 20% of genes with heritability >0.9 had an eQTL in developing worms. Surprisingly, only 10% was found in old worms. Using a multimarker approach, this percentage increased to almost 30% for both age groups. Comparison of the single marker to a multiple marker eQTL mapping indicated that heritable regulation of gene expression becomes more polygenic in aging worms due to multiple loci and possible epistatic interactions. We conclude that linkage studies should account for the relation between increased polygenic regulation and diminished effects at older ages.

A fitness assay for comparing RNAi effects across multiple C. elegans genotypes

Background

RNAi technology by feeding of E. coli containing dsRNA in C. elegans has significantly contributed to further our understanding of many different fields, including genetics, molecular biology, developmental biology and functional genomics. Most of this research has been carried out in a single genotype or genetic background. However, RNAi effects in one genotype do not reveal the allelic effects that segregate in natural populations and contribute to phenotypic variation.

Results

Here we present a method that allows for rapidly comparing RNAi effects among diverse genotypes at an improved high throughput rate. It is based on assessing the fitness of a population of worms by measuring the rate at which E. coli is consumed. Critically, we demonstrate the analytical power of this method by QTL mapping the loss of RNAi sensitivity (in the germline) in a recombinant inbred population derived from a cross between Bristol and a natural isolate from Hawaii. Hawaii has lost RNAi sensitivity in the germline. We found that polymorphisms in ppw-1 contribute to this loss of RNAi sensitivity, but that other loci are also likely to be important.

Conclusions

In summary, we have established a fast method that improves the throughput of RNAi in liquid, that generates quantitative data, that is easy to implement in most laboratories, and importantly that enables QTL mapping using RNAi.

Genetic architecture of life history traits and environment-specific trade-offs

Life history theory predicts the evolution of trait combinations that enhance fitness, and the occurrence of trade-offs depends in part on the magnitude of variation in growth rate or acquisition. Using recombinant inbred lines, we examined the genetic architecture of age and size at reproduction across abiotic conditions encountered by cultivars and naturalized populations of Brassica rapa. We found that genotypes are plastic to seasonal setting, such that reproduction was accelerated under conditions encountered by summer annual populations and genetic variances for age at reproduction varied across simulated seasonal settings. Using an acquisition-allocation model, we predicted the likelihood of trade-offs. Consistent with predicted relationships, we observed a trade-off where early maturity is associated with small size at maturity under simulated summer and fall annual conditions but not under winter annual conditions. The trade-off in the summer annual setting was observed despite significant genotypic variation in growth rate, which is often expected to decouple age and size at reproduction because rapidly growing genotypes could mature early and attain a larger size relative to slowly growing genotypes that mature later. The absence of a trade-off in the winter setting is presumably attributable to the absence of genotypic differences in age at reproduction. We observed QTL for age at reproduction that jointly regulated size at reproduction in both the summer and fall annual settings, but these QTL were environment-specific (i.e. different QTL contributed to the trade-off in the fall vs. summer annual settings). Thus, at least some of the genetic mechanisms underlying observed trade-offs differed across environments.

Caenorhabditis elegans as a platform for molecular quantitative genetics and the systems biology of natural variation

Over the past 30 years, the characteristics that have made the nematode Caenorhabditis elegans one of the premier animal model systems have also allowed it to emerge as a powerful model system for determining the genetic basis of quantitative traits, particularly for the identification of naturally segregating and/or lab-adapted alleles with large phenotypic effects. To better understand the genetic underpinnings of natural variation in other complex phenotypes, C. elegans is uniquely poised in the emerging field of quantitative systems biology because of the extensive knowledge of cellular and neural bases to such traits. However, perturbations in standing genetic variation and patterns of linkage disequilibrium among loci are likely to limit our ability to tie understanding of molecular function to a broader evolutionary context. Coupling the experimental strengths of the C. elegans system with the ecological advantages of closely related nematodes should provide a powerful means of understanding both the molecular and evolutionary genetics of quantitative traits.

Does thermoregulatory behavior maximize reproductive fitness of natural isolates of Caenorhabditis elegans?

BACKGROUND

A central premise of physiological ecology is that an animal's preferred body temperature should correspond closely with the temperature maximizing performance and Darwinian fitness. Testing this co-adaptational hypothesis has been problematic for several reasons. First, reproductive fitness is the appropriate measure, but is difficult to measure in most animals. Second, no single fitness measure applies to all demographic situations, complicating interpretations. Here we test the co-adaptation hypothesis by studying an organism (Caenorhabditis elegans) in which both fitness and thermal preference can be reliably measured.

RESULTS

We find that natural isolates of C. elegans display a range of mean thermal preferences and also vary in their thermal sensitivities for fitness. Hot-seeking isolates CB4854 and CB4857 prefer temperatures that favor population growth rate (r), whereas the cold-seeking isolate CB4856 prefers temperatures that favor Lifetime Reproductive Success (LRS).

CONCLUSIONS

Correlations between fitness and thermal preference in natural isolates of C. elegans are driven primarily by isolate-specific differences in thermal preference. If these differences are the result of natural selection, then this suggests that the appropriate measure of fitness for use in evolutionary ecology studies might differ even within species, depending on the unique ecological and evolutionary history of each population.

The stress of protein misfolding: from single cells to multicellular organisms

Organisms survive changes in the environment by altering their rates of metabolism, growth, and reproduction. At the same time, the system must ensure the stability and functionality of its macromolecules. Fluctuations in the environment are sensed by highly conserved stress responses and homeostatic mechanisms, and of these, the heat shock response (HSR) represents an essential response to acute and chronic proteotoxic damage. However, unlike the strategies employed to maintain the integrity of the genome, protection of the proteome must be tailored to accommodate the normal flux of nonnative proteins and the differences in protein composition between cells, and among individuals. Moreover, adult cells are likely to have significant differences in the rates of synthesis and clearance that are influenced by intrinsic errors in protein expression, genetic polymorphisms, and fluctuations in physiological and environmental conditions. Here, we will address how protein homeostasis (proteostasis) is achieved at the level of the cell and organism, and how the threshold of the stress response is set to detect and combat protein misfolding. For metazoans, the requirement for coordinated function and growth imposes additional constraints on the detection, signaling, and response to misfolding, and requires that the HSR is integrated into various aspects of organismal physiology, such as lifespan. This is achieved by hierarchical regulation of heat shock factor 1 (HSF1) by the metabolic state of the cell and centralized neuronal control that could allow optimal resource allocation between cells and tissues. We will examine how protein folding quality control mechanisms in individual cells may be integrated into a multicellular level of control, and further, even custom-designed to support individual variability and impose additional constraints on evolutionary adaptation.

Switch between life history strategies due to changes in glycolytic enzyme gene dosage in Saccharomyces cerevisiae

Adaptation is the process whereby a population or species becomes better fitted to its habitat through modifications of various life history traits which can be positively or negatively correlated. The molecular factors underlying these covariations remain to be elucidated. Using Saccharomyces cerevisiae as a model system, we have investigated the effects on life history traits of varying the dosage of genes involved in the transformation of resources into energy. Changing gene dosage for each of three glycolytic enzyme genes (hexokinase 2, phosphoglucose isomerase, and fructose-1,6-bisphosphate aldolase) resulted in variation in enzyme activities, glucose consumption rate, and life history traits (growth rate, carrying capacity, and cell size). However, the range of effects depended on which enzyme was expressed differently. Most interestingly, these changes revealed a genetic trade-off between carrying capacity and cell size, supporting the discovery of two extreme life history strategies already described in yeast populations: the "ants," which have lower glycolytic gene dosage, take up glucose slowly, and have a small cell size but reach a high carrying capacity, and the "grasshoppers," which have higher glycolytic gene dosage, consume glucose more rapidly, and allocate it to a larger cell size but reach a lower carrying capacity. These results demonstrate antagonist pleiotropy for glycolytic genes and show that altered dosage of a single gene drives a switch between two life history strategies in yeast.

Temperature-dependent fecundity associates with latitude in Caenorhabditis briggsae

Populations of organisms separated by latitude provide striking examples of local adaptation, by virtue of ecological gradients that correlate with latitudinal position on the globe. Ambient temperature forms one key ecological variable that varies with latitude, and here we investigate its effects on the fecundity of self-fertilizing nematodes of the species Caenorhabditis briggsae that exhibits strong genetically based differentiation in association with latitude. We find that isogenic strains from a Tropical phylogeographic clade have greater lifetime fecundity when reared at extreme high temperatures and lower lifetime fecundity at extreme low temperatures than do strains from a Temperate phylogeographic clade, consistent with adaptation to local temperature regimes. Further, we determine experimentally that the mechanism underlying reduced fecundity at extreme temperatures differs for low versus high temperature extremes, but that the total number of sperm produced by the gonad is unaffected by rearing temperature. Low rearing temperatures result in facultatively reduced oocyte production by hermaphrodites, whereas extreme high temperatures experienced during development induce permanent defects in sperm fertility. Available and emerging genetic tools for this organism will permit the characterization of the evolutionary genetic basis to this putative example of adaptation in latitudinally separated populations.

A genome-wide library of CB4856/N2 introgression lines of Caenorhabditis elegans

Recombinant inbred lines (RILs) derived from Caenorhabditis elegans wild-type N2 and CB4856 are increasingly being used for mapping genes underlying complex traits. To speed up mapping and gene discovery, introgression lines (ILs) offer a powerful tool for more efficient QTL identification. We constructed a library of 90 ILs, each carrying a single homozygous CB4856 genomic segment introgressed into the genetic background of N2. The ILs were genotyped by 123 single-nucleotide polymorphism (SNP) markers. The proportion of the CB4856 segments in most lines does not exceed 3%, and together the introgressions cover 96% of the CB4856 genome. The value of the IL library was demonstrated by identifying novel loci underlying natural variation in two ageing-related traits, i.e. lifespan and pharyngeal pumping rate. Bin mapping of lifespan resulted in six QTLs, which all have a lifespan-shortening effect on the CB4856 allele. We found five QTLs for the decrease in pumping rate, of which four colocated with QTLs found for average lifespan. This suggests pleiotropic or closely linked QTL associated with lifespan and pumping rate. Overall, the presented IL library provides a versatile resource toward easier and efficient fine mapping and functional analyses of loci and genes underlying complex traits in C. elegans.

Thermotaxis is a robust mechanism for thermoregulation in Caenorhabditis elegans nematodes

Many biochemical networks are robust to variations in network or stimulus parameters. Although robustness is considered an important design principle of such networks, it is not known whether this principle also applies to higher-level biological processes such as animal behavior. In thermal gradients, Caenorhabditis elegans uses thermotaxis to bias its movement along the direction of the gradient. Here we develop a detailed, quantitative map of C. elegans thermotaxis and use these data to derive a computational model of thermotaxis in the soil, a natural environment of C. elegans. This computational analysis indicates that thermotaxis enables animals to avoid temperatures at which they cannot reproduce, to limit excursions from their adapted temperature, and to remain relatively close to the surface of the soil, where oxygen is abundant. Furthermore, our analysis reveals that this mechanism is robust to large variations in the parameters governing both worm locomotion and temperature fluctuations in the soil. We suggest that, similar to biochemical networks, animals evolve behavioral strategies that are robust, rather than strategies that rely on fine tuning of specific behavioral parameters.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

Recombinational landscape and population genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.

C. elegans is a model system for diverse fields of biology, but its ability to serve as a model for quantitative trait gene mapping depends on its recombination rate in the laboratory and in nature. The latter is a function of how worms mate and migrate in the wild. We examined the patterns of recombination in a population that we put through thousands of meioses in the laboratory and in a collection of strains isolated from nature. The data suggest that meiotic recombination rate is highly regular in worms, with discrete domains whose boundaries we identify. The pattern in natural strains suggests that population structure, population size, outcrossing rate, and selection combine to suppress the overall effects of recombination. Moreover, some “wild” strains appear to be laboratory contaminants. Nevertheless, the history of recombination in wild worms is sufficient to permit correlations between genotype and phenotype to pinpoint the loci responsible for phenotypic variation.

Integrating the stress response: lessons for neurodegenerative diseases from C. elegans

All cells possess surveillance and homeostatic mechanisms to adjust protein biogenesis to the demands of growth, differentiation, ageing and environmental stress. However, under certain circumstances, these mechanisms fail to adequately respond to proteotoxic imbalances and result in the accumulation of misfolded proteins. In humans, this can lead to neurodegeneration and other protein conformational diseases. To protect itself, the cell employs highly conserved stress responses and chaperone networks to maintain protein-folding homeostasis (proteostasis). Although the regulation of stress responses, such as the heat-shock response, and of proteostasis have been widely considered to be cell autonomous, recent studies using Caenorhabditis elegans have shown that these processes are regulated by neuronal signaling and endocrine pathways and integrated into other functions of the organism. The hierarchical control of the cellular proteostasis machinery affords insight into the organization of stress regulatory networks in multicellular organisms and offers novel targets for the treatment of human protein conformational diseases.