Beyond induced mutants: using worms to study natural variation in genetic pathways

REVIEW OF THE DAUER HYPOTHESIS: WHAT NON-PARASITIC SPECIES CAN TELL US ABOUT THE EVOLUTION OF PARASITISM

Parasitic lineages have acquired suites of new traits compared to their nearest free-living relatives. When and why did these traits arise? We can envision lineages evolving through multiple stable intermediate steps such as a series of increasingly exploitative species interactions. This view allows us to use non-parasitic species that approximate those intermediate steps to uncover the timing and original function of parasitic traits, knowledge critical to understanding the evolution of parasitism. The dauer hypothesis proposes that free-living nematode lineages evolved into parasites through two intermediate steps, phoresy and necromeny. Here we delve into the proposed steps of the dauer hypothesis by collecting and organizing data from genetic, behavioral, and ecological studies in a range of nematode species. We argue that hypotheses on the evolution of parasites will be strengthened by complementing comparative genomic studies with ecological studies on non-parasites that approximate intermediate steps.

Neuro-genetic plasticity of Caenorhabditis elegans behavioral thermal tolerance

BACKGROUND

Animal responses to thermal stimuli involve intricate contributions of genetics, neurobiology and physiology, with temperature variation providing a pervasive environmental factor for natural selection. Thermal behavior thus exemplifies a dynamic trait that requires non-trivial phenotypic summaries to appropriately capture the trait in response to a changing environment. To characterize the deterministic and plastic components of thermal responses, we developed a novel micro-droplet assay of nematode behavior that permits information-dense summaries of dynamic behavioral phenotypes as reaction norms in response to increasing temperature (thermal tolerance curves, TTC).

RESULTS

We found that C. elegans TTCs shift predictably with rearing conditions and developmental stage, with significant differences between distinct wildtype genetic backgrounds. Moreover, after screening TTCs for 58 C. elegans genetic mutant strains, we determined that genes affecting thermosensation, including cmk-1 and tax-4, potentially play important roles in the behavioral control of locomotion at high temperature, implicating neural decision-making in TTC shape rather than just generalized physiological limits. However, expression of the transient receptor potential ion channel TRPA-1 in the nervous system is not sufficient to rescue rearing-dependent plasticity in TTCs conferred by normal expression of this gene, indicating instead a role for intestinal signaling involving TRPA-1 in the adaptive plasticity of thermal performance.

CONCLUSIONS

These results implicate nervous system and non-nervous system contributions to behavior, in addition to basic cellular physiology, as key mediators of evolutionary responses to selection from temperature variation in nature.

Genetic background modifies phenotypic and transcriptional responses in a C. elegans model of α-synuclein toxicity

Background

Accumulation of protein aggregates are a major hallmark of progressive neurodegenerative disorders such as Parkinson’s disease and Alzheimer’s disease. Transgenic Caenorhabditis elegans nematodes expressing the human synaptic protein α-synuclein in body wall muscle show inclusions of aggregated protein, which affects similar genetic pathways as in humans. It is not however known how the effects of α-synuclein expression in C. elegans differs among genetic backgrounds. Here, we compared gene expression patterns and investigated the phenotypic consequences of transgenic α-synuclein expression in five different C. elegans genetic backgrounds.

Results

Transcriptome analysis indicates that α-synuclein expression effects pathways associated with nutrient storage, lipid transportation and ion exchange and that effects vary depending on the genetic background. These gene expression changes predict that a range of phenotypes will be affected by α-synuclein expression. We confirm this, showing that α-synuclein expression delayed development, reduced lifespan, increased rate of matricidal hatching, and slows pharyngeal pumping. Critically, these phenotypic effects depend on the genetic background and coincide with the core changes in gene expression.

Conclusions

Together, our results show genotype-specific effects and core alterations in both gene expression and in phenotype in response to α-synuclein expression. We conclude that the effects of α-synuclein expression are substantially modified by the genetic background, illustrating that genetic background needs to be considered in C. elegans models of neurodegenerative disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5597-1) contains supplementary material, which is available to authorized users.

Forward and reverse genetics approaches to uncover metabolic aging pathways in Caenorhabditis elegans

The biological mechanisms of aging have been studied in depth and prominent findings in this field promote the development of new therapies for age-associated disorders. Various model organisms are used for research on aging; among these, the nematode Caenorhabditis elegans has been widely used and has provided valuable knowledge in determining the regulatory mechanisms driving the aging process. Many genes involved in lifespan regulation are associated with metabolic pathways and are influenced by genetic and environmental factors. In line with this, C. elegans provides a promising platform to study such gene by environment interactions, in either a reverse or forward genetics approach. In this review, we discuss longevity mechanisms related to metabolic networks that have been discovered in C. elegans. We also highlight the use of wild populations to study the complex genetic basis of natural variation for quantitative traits that mediate longevity.

Ras/MAPK Modifier Loci Revealed by eQTL in Caenorhabditis elegans

The oncogenic Ras/MAPK pathway is evolutionarily conserved across metazoans. Yet, almost all our knowledge on this pathway comes from studies using single genetic backgrounds, whereas mutational effects can be highly background dependent. Therefore, we lack insight in the interplay between genetic backgrounds and the Ras/MAPK-signaling pathway. Here, we used a Caenorhabditis elegans RIL population containing a gain-of-function mutation in the Ras/MAPK-pathway gene let-60 and measured how gene expression regulation is affected by this mutation. We mapped eQTL and found that the majority (∼73%) of the 1516 detected cis-eQTL were not specific for the let-60 mutation, whereas most (∼76%) of the 898 detected trans-eQTL were associated with the let-60 mutation. We detected six eQTL trans-bands specific for the interaction between the genetic background and the mutation, one of which colocalized with the polymorphic Ras/MAPK modifier amx-2. Comparison between transgenic lines expressing allelic variants of amx-2 showed the involvement of amx-2 in 79% of the trans-eQTL for genes mapping to this trans-band. Together, our results have revealed hidden loci affecting Ras/MAPK signaling using sensitized backgrounds in C. elegans. These loci harbor putative polymorphic modifier genes that would not have been detected using mutant screens in single genetic backgrounds.

Genetic variation in neurodegenerative diseases and its accessibility in the model organism Caenorhabditis elegans

BACKGROUND

Neurodegenerative diseases (NGDs) such as Alzheimer's and Parkinson's are debilitating and largely untreatable conditions strongly linked to age. The clinical, neuropathological, and genetic components of NGDs indicate that neurodegeneration is a complex trait determined by multiple genes and by the environment.

MAIN BODY

The symptoms of NGDs differ among individuals due to their genetic background, and this variation affects the onset and progression of NGD and NGD-like states. Such genetic variation affects the molecular and cellular processes underlying NGDs, leading to differential clinical phenotypes. So far, we have a limited understanding of the mechanisms of individual background variation. Here, we consider how variation between genetic backgrounds affects the mechanisms of aging and proteostasis in NGD phenotypes. We discuss how the nematode Caenorhabditis elegans can be used to identify the role of variation between genetic backgrounds. Additionally, we review advances in C. elegans methods that can facilitate the identification of NGD regulators and/or networks.

CONCLUSION

Genetic variation both in disease genes and in regulatory factors that modulate onset and progression of NGDs are incompletely understood. The nematode C. elegans represents a valuable system in which to address such questions.

Systemic Regulation of RAS/MAPK Signaling by the Serotonin Metabolite 5-HIAA

Human cancer is caused by the interplay of mutations in oncogenes and tumor suppressor genes and inherited variations in cancer susceptibility genes. While many of the tumor initiating mutations are well characterized, the effect of genetic background variation on disease onset and progression is less understood. We have used C. elegans genetics to identify genetic modifiers of the oncogenic RAS/MAPK signaling pathway. Quantitative trait locus analysis of two highly diverged C. elegans isolates combined with allele swapping experiments identified the polymorphic monoamine oxidase A (MAOA) gene amx-2 as a negative regulator of RAS/MAPK signaling. We further show that the serotonin metabolite 5-hydroxyindoleacetic acid (5-HIAA), which is a product of MAOA catalysis, systemically inhibits RAS/MAPK signaling in different organs of C. elegans. Thus, MAOA activity sets a global threshold for MAPK activation by controlling 5-HIAA levels. To our knowledge, 5-HIAA is the first endogenous small molecule that acts as a systemic inhibitor of RAS/MAPK signaling.

Highly polygenic variation in environmental perception determines dauer larvae formation in growing populations of Caenorhabditis elegans

Background

Determining how complex traits are genetically controlled is a requirement if we are to predict how they evolve and how they might respond to selection. This requires understanding how distinct, and often more simple, life history traits interact and change in response to environmental conditions. In order to begin addressing such issues, we have been analyzing the formation of the developmentally arrested dauer larvae of Caenorhabditis elegans under different conditions.

Results

We find that 18 of 22 previously identified quantitative trait loci (QTLs) affecting dauer larvae formation in growing populations, assayed by determining the number of dauer larvae present at food patch exhaustion, can be recovered under various environmental conditions. We also show that food patch size affects both the ability to detect QTLs and estimates of effect size, and demonstrate that an allele of nath-10 affects dauer larvae formation in growing populations. To investigate the component traits that affect dauer larvae formation in growing populations we map, using the same introgression lines, QTLs that affect dauer larvae formation in response to defined amounts of pheromone. This identifies 36 QTLs, again demonstrating the highly polygenic nature of the genetic variation underlying dauer larvae formation.

Conclusions

These data indicate that QTLs affecting the number of dauer larvae at food exhaustion in growing populations of C. elegans are highly reproducible, and that nearly all can be explained by variation affecting dauer larvae formation in response to defined amounts of pheromone. This suggests that most variation in dauer larvae formation in growing populations is a consequence of variation in the perception of the food and pheromone environment (i.e. chemosensory variation) and in the integration of these cues.

Widespread genomic incompatibilities in Caenorhabditis elegans

In the Bateson-Dobzhansky-Muller (BDM) model of speciation, incompatibilities emerge from the deleterious interactions between alleles that are neutral or advantageous in the original genetic backgrounds, i.e., negative epistatic effects. Within species such interactions are responsible for outbreeding depression and F2 (hybrid) breakdown. We sought to identify BDM incompatibilities in the nematode Caenorhabditis elegans by looking for genomic regions that disrupt egg laying; a complex, highly regulated, and coordinated phenotype. Investigation of introgression lines and recombinant inbred lines derived from the isolates CB4856 and N2 uncovered multiple incompatibility quantitative trait loci (QTL). These QTL produce a synthetic egg-laying defective phenotype not seen in CB4856 and N2 nor in other wild isolates. For two of the QTL regions, results are inconsistent with a model of pairwise interaction between two loci, suggesting that the incompatibilities are a consequence of complex interactions between multiple loci. Analysis of additional life history traits indicates that the QTL regions identified in these screens are associated with effects on other traits such as lifespan and reproduction, suggesting that the incompatibilities are likely to be deleterious. Taken together, these results indicate that numerous BDM incompatibilities that could contribute to reproductive isolation can be detected and mapped within C. elegans.

Genome-wide variations in a natural isolate of the nematode Caenorhabditis elegans

Background

Increasing genetic and phenotypic differences found among natural isolates of C. elegans have encouraged researchers to explore the natural variation of this nematode species.

Results

Here we report on the identification of genomic differences between the reference strain N2 and the Hawaiian strain CB4856, one of the most genetically distant strains from N2. To identify both small- and large-scale genomic variations (GVs), we have sequenced the CB4856 genome using both Roche 454 (~400 bps single reads) and Illumina GA DNA sequencing methods (101 bps paired-end reads). Compared to previously described variants (available in WormBase), our effort uncovered twice as many single nucleotide variants (SNVs) and increased the number of small InDels almost 20-fold. Moreover, we identified and validated large insertions, most of which range from 150 bps to 1.2 kb in length in the CB4856 strain. Identified GVs had a widespread impact on protein-coding sequences, including 585 single-copy genes that have associated severe phenotypes of reduced viability in RNAi and genetics studies. Sixty of these genes are homologs of human genes associated with diseases. Furthermore, our work confirms previously identified GVs associated with differences in behavioural and biological traits between the N2 and CB4856 strains.

Conclusions

The identified GVs provide a rich resource for future studies that aim to explain the genetic basis for other trait differences between the N2 and CB4856 strains.

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.

A heritable antiviral RNAi response limits Orsay virus infection in Caenorhabditis elegans N2

Orsay virus (OrV) is the first virus known to be able to complete a full infection cycle in the model nematode species Caenorhabditis elegans. OrV is transmitted horizontally and its infection is limited by antiviral RNA interference (RNAi). However, we have no insight into the kinetics of OrV replication in C. elegans. We developed an assay that infects worms in liquid, allowing precise monitoring of the infection. The assay revealed a dual role for the RNAi response in limiting Orsay virus infection in C. elegans. Firstly, it limits the progression of the initial infection at the step of recognition of dsRNA. Secondly, it provides an inherited protection against infection in the offspring. This establishes the heritable RNAi response as anti-viral mechanism during OrV infections in C. elegans. Our results further illustrate that the inheritance of the anti-viral response is important in controlling the infection in the canonical wild type Bristol N2. The OrV replication kinetics were established throughout the worm life-cycle, setting a standard for further quantitative assays with the OrV-C. elegans infection model.

A rapid and massive gene expression shift marking adolescent transition in C. elegans

Organismal development is the most dynamic period of the life cycle, yet we have only a rough understanding of the dynamics of gene expression during adolescent transition. Here we show that adolescence in Caenorhabditis elegans is characterized by a spectacular expression shift of conserved and highly polymorphic genes. Using a high resolution time series we found that in adolescent worms over 10,000 genes changed their expression. These genes were clustered according to their expression patterns. One cluster involved in chromatin remodelling showed a brief up-regulation around 50 h post-hatch. At the same time a spectacular shift in expression was observed. Sequence comparisons for this cluster across many genotypes revealed diversifying selection. Strongly up-regulated genes showed signs of purifying selection in non-coding regions, indicating that adolescence-active genes are constrained on their regulatory properties. Our findings improve our understanding of adolescent transition and help to eliminate experimental artefacts due to incorrect developmental timing.

Specialist versus generalist life histories and nucleotide diversity in Caenorhabditis nematodes

Species with broad ecological amplitudes with respect to a key focal resource, niche generalists, should maintain larger and more connected populations than niche specialists, leading to the prediction that nucleotide diversity will be lower and more subdivided in specialists relative to their generalist relatives. This logic describes the specialist-generalist variation hypothesis (SGVH). Some outbreeding species of Caenorhabditis nematodes use a variety of invertebrate dispersal vectors and have high molecular diversity. By contrast, Caenorhabditis japonica lives in a strict association and synchronized life cycle with its dispersal host, the shield bug Parastrachia japonensis, itself a diet specialist. Here, we characterize sequence variation for 20 nuclear loci to investigate how C. japonica's life history shapes nucleotide diversity. We find that C. japonica has more than threefold lower polymorphism than other outbreeding Caenorhabditis species, but that local populations are not genetically disconnected. Coupled with its restricted range, we propose that its specialist host association contributes to a smaller effective population size and lower genetic variation than host generalist Caenorhabditis species with outbreeding reproductive modes. A literature survey of diverse organisms provides broader support for the SGVH. These findings encourage further testing of ecological and evolutionary hypotheses with comparative population genetics in Caenorhabditis and other taxa.

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.

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.

Natural genetic variation determines susceptibility to aggregation or toxicity in a C. elegans model for polyglutamine disease

Background

Monogenic gain-of-function protein aggregation diseases, including Huntington’s disease, exhibit substantial variability in age of onset, penetrance, and clinical symptoms, even between individuals with similar or identical mutations. This difference in phenotypic expression of proteotoxic mutations is proposed to be due, at least in part, to the variability in genetic background. To address this, we examined the role of natural variation in defining the susceptibility of genetically diverse individuals to protein aggregation and toxicity, using the Caenorhabditis elegans polyglutamine model.

Results

Introgression of polyQ40 into three wild genetic backgrounds uncovered wide variation in onset of aggregation and corresponding toxicity, as well as alteration in the cell-specific susceptibility to aggregation. To further dissect these relationships, we established a panel of 21 recombinant inbred lines that showed a broad range of aggregation phenotypes, independent of differences in expression levels. We found that aggregation is a transgressive trait, and does not always correlate with measures of toxicity, such as early onset of muscle dysfunction, egg-laying deficits, or reduced lifespan. Moreover, distinct measures of proteotoxicity were independently modified by the genetic background.

Conclusions

Resistance to protein aggregation and the ability to restrict its associated cellular dysfunction are independently controlled by the natural variation in genetic background, revealing important new considerations in the search for targets for therapeutic intervention in conformational diseases. Thus, our C. elegans model can serve as a powerful tool to dissect the contribution of natural variation to individual susceptibility to proteotoxicity.

Please see related commentary by Kaeberlein, http://www.biomedcentral.com/1741-7007/11/102.

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.

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.

Population dynamics and habitat sharing of natural populations of Caenorhabditis elegans and C. briggsae

BACKGROUND

The nematode Caenorhabditis elegans is a major model organism in laboratory biology. Very little is known, however, about its ecology, including where it proliferates. In the past, C. elegans was mainly isolated from human-made compost heaps, where it was overwhelmingly found in the non-feeding dauer diapause stage.

RESULTS

C. elegans and C. briggsae were found in large, proliferating populations in rotting plant material (fruits and stems) in several locations in mainland France. Both species were found to co-occur in samples isolated from a given plant species. Population counts spanned a range from one to more than 10,000 Caenorhabditis individuals on a single fruit or stem. Some populations with an intermediate census size (10 to 1,000) contained no dauer larvae at all, whereas larger populations always included some larvae in the pre-dauer or dauer stages. We report on associated micro-organisms, including pathogens. We systematically sampled a spatio-temporally structured set of rotting apples in an apple orchard in Orsay over four years. C. elegans and C. briggsae were abundantly found every year, but their temporal distributions did not coincide. C. briggsae was found alone in summer, whereas both species co-occurred in early fall and C. elegans was found alone in late fall. Competition experiments in the laboratory at different temperatures show that C. briggsae out-competes C. elegans at high temperatures, whereas C. elegans out-competes C. briggsae at lower temperatures.

CONCLUSIONS

C. elegans and C. briggsae proliferate in the same rotting vegetal substrates. In contrast to previous surveys of populations in compost heaps, we found fully proliferating populations with no dauer larvae. The temporal sharing of the habitat by the two species coincides with their temperature preference in the laboratory, with C. briggsae populations growing faster than C. elegans at higher temperatures, and vice at lower temperatures.

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 phylogeny and molecular barcodes for Caenorhabditis, with numerous new species from rotting fruits

BACKGROUND

The nematode Caenorhabditis elegans is a major laboratory model in biology. Only ten Caenorhabditis species were available in culture at the onset of this study. Many of them, like C. elegans, were mostly isolated from artificial compost heaps, and their more natural habitat was unknown.

RESULTS

Caenorhabditis nematodes were found to be proliferating in rotten fruits, flowers and stems. By collecting a large worldwide set of such samples, 16 new Caenorhabditis species were discovered. We performed mating tests to establish biological species status and found some instances of semi-fertile or sterile hybrid progeny. We established barcodes for all species using ITS2 rDNA sequences. By obtaining sequence data for two rRNA and nine protein-coding genes, we determined the likely phylogenetic relationships among the 26 species in culture. The new species are part of two well-resolved sister clades that we call the Elegans super-group and the Drosophilae super-group. We further scored phenotypic characters such as reproductive mode, mating behavior and male tail morphology, and discuss their congruence with the phylogeny. A small space between rays 2 and 3 evolved once in the stem species of the Elegans super-group; a narrow fan and spiral copulation evolved once in the stem species of C. angaria, C. sp. 8 and C. sp. 12. Several other character changes occurred convergently. For example, hermaphroditism evolved three times independently in C. elegans, C. briggsae and C. sp. 11. Several species can co-occur in the same location or even the same fruit. At the global level, some species have a cosmopolitan distribution: C. briggsae is particularly widespread, while C. elegans and C. remanei are found mostly or exclusively in temperate regions, and C. brenneri and C. sp. 11 exclusively in tropical zones. Other species have limited distributions, for example C. sp. 5 appears to be restricted to China, C. sp. 7 to West Africa and C. sp. 8 to the Eastern United States.

CONCLUSIONS

Caenorhabditis are "fruit worms", not soil nematodes. The 16 new species provide a resource and their phylogeny offers a framework for further studies into the evolution of genomic and phenotypic characters.

MicroRNA sequence variation potentially contributes to within-species functional divergence in the nematode Caenorhabditis briggsae

Mounting evidence points to differences in gene regulation as a major source of phenotypic variation. MicroRNA-mediated post-transcriptional regulation has emerged recently as a key factor controlling gene activity during development. MicroRNA genes are abundant in genomes, acting as managers of gene expression by directing translational repression. Thus, understanding the role of microRNA sequence variation within populations is essential for fully dissecting the origin and maintenance of phenotypic diversity in nature. In this study, we investigate allelic variation at microRNA loci in the nematode Caenorhabditis briggsae, a close relative of C. elegans. Phylogeographic structure in C. briggsae partitions most strains from around the globe into a "temperate" or a "tropical" clade, with a few strains having divergent, geographically restricted genotypes. Remarkably, strains that follow this latitudinal dichotomy also differ in temperature-associated fitness. With this phylogeographic pattern in mind, we examined polymorphisms in 18 miRNAs in a global sample of C. briggsae isolates and tested whether newly isolated strains conform to this phylogeography. Surprisingly, nucleotide diversity is relatively high in this class of gene that generally experiences strong purifying selection. In particular, we find that miRNAs in C. briggsae are substantially more polymorphic than in Arabidopsis thaliana, despite similar background levels of neutral site diversity between the two species. We find that some mutations suggest functional divergence on the basis of requirements for target site recognition and computational prediction of the effects of the polymorphisms on RNA folding. These findings demonstrate the potential for miRNA polymorphisms to contribute to phenotypic variation within a species. Sequences were deposited in GenBank under accession nos. JN251323-JN251744.

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.

Heritable determinants of male fertilization success in the nematode Caenorhabditis elegans

BACKGROUND

Sperm competition is a driving force in the evolution of male sperm characteristics in many species. In the nematode Caenorhabditis elegans, larger male sperm evolve under experimentally increased sperm competition and larger male sperm outcompete smaller hermaphrodite sperm for fertilization within the hermaphrodite reproductive tract. To further elucidate the relative importance of sperm-related traits that contribute to differential reproductive success among males, we quantified within- and among-strain variation in sperm traits (size, rate of production, number transferred, competitive ability) for seven male genetic backgrounds known previously to differ with respect to some sperm traits. We also quantified male mating ability in assays for rates of courtship and successful copulation, and then assessed the roles of these pre- and post-mating traits in first- and second-male fertilization success.

RESULTS

We document significant variation in courtship ability, mating ability, sperm size and sperm production rate. Sperm size and production rate were strong indicators of early fertilization success for males that mated second, but male genetic backgrounds conferring faster sperm production make smaller sperm, despite virgin males of all genetic backgrounds transferring indistinguishable numbers of sperm to mating partners.

CONCLUSIONS

We have demonstrated that sperm size and the rate of sperm production represent dominant factors in determining male fertilization success and that C. elegans harbors substantial heritable variation for traits contributing to male reproductive success. C. elegans provides a powerful, tractable system for studying sexual selection and for dissecting the genetic basis and evolution of reproduction-related traits.

Heat avoidance is regulated by transient receptor potential (TRP) channels and a neuropeptide signaling pathway in Caenorhabditis elegans

The ability to avoid noxious extremes of hot and cold is critical for survival and depends on thermal nociception. The TRPV subset of transient receptor potential (TRP) channels is heat activated and proposed to be responsible for heat detection in vertebrates and fruit flies. To gain insight into the genetic and neural basis of thermal nociception, we developed assays that quantify noxious heat avoidance in the nematode Caenorhabditis elegans and used them to investigate the genetic basis of this behavior. First, we screened mutants for 18 TRP channel genes (including all TRPV orthologs) and found only minor defects in heat avoidance in single and selected double and triple mutants, indicating that other genes are involved. Next, we compared two wild isolates of C. elegans that diverge in their threshold for heat avoidance and linked this phenotypic variation to a polymorphism in the neuropeptide receptor gene npr-1. Further analysis revealed that loss of either the NPR-1 receptor or its ligand, FLP-21, increases the threshold for heat avoidance. Cell-specific rescue of npr-1 implicates the interneuron RMG in the circuit regulating heat avoidance. This neuropeptide signaling pathway operates independently of the TRPV genes, osm-9 and ocr-2, since mutants lacking npr-1 and both TRPV channels had more severe defects in heat avoidance than mutants lacking only npr-1 or both osm-9 and ocr-2. Our results show that TRPV channels and the FLP-21/NPR-1 neuropeptide signaling pathway determine the threshold for heat avoidance in C. elegans.

Extremely high molecular diversity within the East Asian nematode Caenorhabditis sp. 5

Most relatives of the self-fertilizing hermaphroditic nematode model organism Caenorhabditis elegans reproduce via obligate outbreeding between males and females, which also represents the ancestral mode of reproduction within the genus. However, little is known about the scope of genetic diversity and differentiation within such gonochoristic species, especially those found outside of temperate Europe and North America. It is critical to understand the evolutionary processes operating in these species to provide a framework for deciphering the evolution of hermaphroditism and a baseline for the application of outcrossing Caenorhabditis to problems in evolutionary genetics. Here, we investigate for the first time molecular sequence variation for Caenorhabditis sp. 5, a species found commonly in eastern Asia. We identify enormous levels of standing genetic variation that approach the levels observed in the marine broadcast-spawning sea squirt, Ciona savignyi. Although we document significant isolation by distance, we demonstrate that the high polymorphism within C. sp. 5 is not because of strong differentiation among populations or to the presence of cryptic species. These findings illustrate that molecular population genetic approaches to studying obligately outbreeding species of Caenorhabditis will prove powerful in identifying and characterizing functionally and evolutionarily important features of the genome.

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.

Molecular population genetics and phenotypic sensitivity to ethanol for a globally diverse sample of the nematode Caenorhabditis briggsae

New genomic resources and genetic tools of the past few years have advanced the nematode genus Caenorhabditis as a model for comparative biology. However, understanding of natural genetic variation at molecular and phenotypic levels remains rudimentary for most species in this genus, and for C. briggsae in particular. Here we characterize phenotypic variation in C. briggsae's sensitivity to the potentially important and variable environmental toxin, ethanol, for globally diverse strains. We also quantify nucleotide variation in a new sample of 32 strains from four continents, including small islands, and for the closest-known relative of this species (C. sp. 9). We demonstrate that C. briggsae exhibits little heritable variation for the effects of ethanol on the norm of reaction for survival and reproduction. Moreover, C. briggsae does not differ significantly from C. elegans in our assays of its response to this substance that both species likely encounter regularly in habitats of rotting fruit and vegetation. However, we uncover drastically more molecular genetic variation than was known previously for this species, despite most strains, including all island strains, conforming to the broad biogeographic patterns described previously. Using patterns of sequence divergence between populations and between species, we estimate that the self-fertilizing mode of reproduction by hermaphrodites in C. briggsae likely evolved sometime between 0.9 and 10 million generations ago. These insights into C. briggsae's natural history and natural genetic variation greatly expand the potential of this organism as an emerging model for studies in molecular and quantitative genetics, the evolution of development, and ecological genetics.

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

A sensitized genetic background reveals evolution near the terminus of the Caenorhabditis germline sex determination pathway

Caenorhabditis elegans and Caenorhabditis briggsae are both self-fertile hermaphroditic nematodes that evolved independently from male/female ancestors. In C. elegans, FEM-1, FEM-2, and FEM-3 specify male fates by promoting proteolysis of the male-repressing transcription factor, TRA-1. Phenotypes of tra-1 and fem mutants are consistent with this simple linear model in the soma, but not in the germline. While both XX and XO tra-1(lf) mutants have functional male somas, they produce both sperm and oocytes. Further, all three tra-1; fem double mutants retain the expected male soma, but make only oocytes (the germline fem phenotype). Thus, a poorly characterized tra-1 activity is important for sustained male spermatogenesis, and the fem genes affect germline sexual fate independently of their role in regulating TRA-1. C. briggsae tra-1 mutants are phenotypically identical to their C. elegans counterparts, while the fem mutants differ in the germline: XX and XO C. elegans fem mutants are true females, but in C. briggsae they are self-fertile hermaphrodites. To further explore how C. briggsae hermaphrodites regulate germline sex, we analyzed Cb-tra-1/Cb-fem interactions. Cb-tra-1 is fully epistatic to Cb-fem-2 in the germline, unlike the orthologous C. elegans combination. In contrast, Cb-fem-3 shifts the Cb-tra-1(lf) germline phenotype to that of a nearly normal hermaphrodite in the context of a male somatic gonad. This suggests that Cb-fem-3 is epistatic to Cb-tra-1(lf) (as in C. elegans), and that the normal control of C. briggsae XX spermatogenesis targets Cb-tra-1-independent factors downstream of Cb-fem-3. The effect of Cb-fem-3(lf) on Cb-tra-1(lf) is not mediated by change in the expression of Cb-fog-3, a likely direct germline target of Cb-tra-1. As Cb-fem-2 and Cb-fem-3 have identical single mutant phenotypes, Cb-tra-1 provides a sensitized background that reveals differences in how they promote male germline development. These results represent another way in which C. briggsae germline sex determination is incongruent with that of the outwardly similar C. elegans.

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.

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.

Rapid high resolution single nucleotide polymorphism-comparative genome hybridization mapping in Caenorhabditis elegans

We have developed a significantly improved and simplified method for high-resolution mapping of phenotypic traits in Caenorhabditis elegans using a combination of single nucleotide polymorphisms (SNPs) and oligo array comparative genome hybridization (array CGH). We designed a custom oligonucleotide array using a subset of confirmed SNPs between the canonical wild-type Bristol strain N2 and the Hawaiian isolate CB4856, populated with densely overlapping 50-mer probes corresponding to both N2 and CB4856 SNP sequences. Using this method a mutation can be mapped to a resolution of approximately 200 kb in a single genetic cross. Six mutations representing each of the C. elegans chromosomes were detected unambiguously and at high resolution using genomic DNA from populations derived from as few as 100 homozygous mutant segregants of mutant N2/CB4856 heterozygotes. Our method completely dispenses with the PCR, restriction digest, and gel analysis of standard SNP mapping and should be easy to extend to any organism with interbreeding strains. This method will be particularly powerful when applied to difficult or hard-to-map low-penetrance phenotypes. It should also be possible to map polygenic traits using this method.