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

Natural variation in infection specificity of Caenorhabditis briggsae isolates by two RNA viruses

Antagonistic relationships such as host-virus interactions potentially lead to rapid evolution and specificity in interactions. The Orsay virus is so far the only horizontal virus naturally infecting the nematode C. elegans. In contrast, several related RNA viruses infect its congener C. briggsae, including Santeuil (SANTV) and Le Blanc (LEBV) viruses. Here we focus on the host's intraspecific variation in sensitivity to these two intestinal viruses. Many temperate-origin C. briggsae strains, including JU1264 and JU1498, are sensitive to both, while many tropical strains, such as AF16, are resistant to both. Interestingly, some C. briggsae strains exhibit a specific resistance, such as the HK104 strain, specifically resistant to LEBV. The viral sensitivity pattern matches the strains' geographic and genomic relationships. The heavily infected strains mount a seemingly normal small RNA response that is insufficient to suppress viral infection, while the resistant strains show no small RNA response, suggesting an early block in viral entry or replication. We use a genetic approach from the host side to map genomic regions participating in viral resistance polymorphisms. Using Advanced Intercrossed Recombinant Inbred Lines (RILs) between virus-resistant AF16 and SANTV-sensitive HK104, we detect Quantitative Trait Loci (QTLs) on chromosomes IV and III. Building RILs between virus-sensitive JU1498 and LEBV-resistant HK104 followed by bulk segregant analysis, we identify a chromosome II QTL. In both cases, further introgressions of the regions confirmed the QTLs. This diversity provides an avenue for studying virus entry, replication, and exit mechanisms, as well as host-virus specificity and the host response to a specific virus infection.

Environmental radiation exposure at Chornobyl has not systematically affected the genomes or chemical mutagen tolerance phenotypes of local worms

The 1986 disaster at the Chornobyl Nuclear Power Plant transformed the surrounding region into the most radioactive landscape known on the planet. Whether or not this sudden environmental shift selected for species, or even individuals within a species, that are naturally more resistant to mutagen exposure remains an open question. In this study, we collected, cultured, and cryopreserved 298 wild nematode isolates from areas varying in radioactivity within the Chornobyl Exclusion Zone. We sequenced and assembled genomes de novo for 20 Oscheius tipulae strains, analyzed their genomes for evidence of recent mutation acquisition in the field, and observed no evidence of an association between mutation and radioactivity at the sites of collection. Multigenerational exposure of each of these strains to several chemical mutagens in the lab revealed that strains vary heritably in tolerance to each mutagen, but mutagen tolerance cannot be predicted based on the radiation levels at collection sites, and Chornobyl isolates were not systematically more resistant than strains from undisturbed habitats. In sum, the absence of mutational signatures does not reflect unique capacity for tolerating DNA damage.

rec-1 loss of function increases recombination in the central gene clusters at the expense of autosomal pairing centers

Meiotic control of crossover (CO) number and position is critical for homologous chromosome segregation and organismal fertility, recombination of parental genotypes, and the generation of novel genetic combinations. We here characterize the recombination rate landscape of a rec-1 loss of function modifier of CO position in Caenorhabditis elegans, one of the first ever modifiers discovered. By averaging CO position across hermaphrodite and male meioses and by genotyping 203 single-nucleotide variants covering about 95% of the genome, we find that the characteristic chromosomal arm-center recombination rate domain structure is lost in the loss of function rec-1 mutant. The rec-1 loss of function mutant smooths the recombination rate landscape but is insufficient to eliminate the nonuniform position of CO. Lower recombination rates in the rec-1 mutant are particularly found in the autosomal arm domains containing the pairing centers. We further find that the rec-1 mutant is of little consequence for organismal fertility and egg viability and thus for rates of autosomal nondisjunction. It nonetheless increases X chromosome nondisjunction rates and thus male appearance. Our findings question the maintenance of recombination rate heritability and genetic diversity among C. elegans natural populations, and they further suggest that manipulating genetic modifiers of CO position will help find quantitative trait loci located in low-recombining genomic regions normally refractory to discovery.

Heritable epigenetic variation facilitates long-term maintenance of epigenetic and genetic variation

How genetic and phenotypic variation are maintained has long been one of the fundamental questions in population and quantitative genetics. A variety of factors have been implicated to explain the maintenance of genetic variation in some contexts (e.g. balancing selection), but the potential role of epigenetic regulation to influence population dynamics has been understudied. It is well recognized that epigenetic regulation, including histone methylation, small RNA expression, and DNA methylation, helps to define differences between cell types and facilitate phenotypic plasticity. In recent years, empirical studies have shown the potential for epigenetic regulation to also be heritable for at least a few generations without selection, raising the possibility that differences in epigenetic regulation can act alongside genetic variation to shape evolutionary trajectories. Heritable differences in epigenetic regulation that arise spontaneously are termed "epimutations." Epimutations differ from genetic mutations in 2 key ways-they occur at a higher rate and the loci at which they occur often revert back to their original state within a few generations. Here, we present an extension of the standard population genetic model with selection to incorporate epigenetic variation arising via epimutation. Our model assumes a diploid, sexually reproducing population with random mating. In addition to spontaneous genetic mutation, we included parameters for spontaneous epimutation and back-epimutation, allowing for 4 potential epialleles at a single locus (2 genetic alleles, each with 2 epigenetic states), each of which affect fitness. We then analyzed the conditions under which stable epialleles were maintained. Our results show that highly reversible epialleles can be maintained in long-term equilibrium under neutral conditions in a manner that depends on the epimutation and back-epimutation rates, which we term epimutation-back-epimutation equilibrium. On the other hand, epialleles that compensate for deleterious mutations cause deviations from the expectations of mutation-selection balance by a simple factor that depends on the epimutation and back-epimutation rates. We also numerically analyze several sets of fitness parameters for which large deviations from mutation-selection balance occur. Together, these results demonstrate that transient epigenetic regulation may be an important factor in the maintenance of both epigenetic and genetic variation in populations.

Genome-wide association and environmental suppression of the mortal germline phenotype of wild C. elegans

The animal germline lineage needs to be maintained along generations. However, some Caenorhabditis elegans wild isolates display a mortal germline phenotype, leading to sterility after several generations at 25°C. Using a genome-wide association approach, we detect a significant peak on chromosome III around 5 Mb, confirmed by introgressions. Thus, a seemingly deleterious genotype is maintained at intermediate frequency in the species. Environmental rescue is a likely explanation, and indeed associated bacteria and microsporidia suppress the phenotype of wild isolates as well as mutants in small RNA inheritance (nrde-2) and histone modifications (set-2). Escherichia coli strains of the K-12 lineage suppress the phenotype compared to B strains. By shifting a wild strain from E. coli K-12 to E. coli B, we find that memory of the suppressing condition is maintained over several generations. Thus, the mortal germline phenotype of wild C. elegans is in part revealed by laboratory conditions and may represent variation in epigenetic inheritance and environmental interactions. This study also points to the importance of non-genetic memory in the face of environmental variation.

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.

Selection and the direction of phenotypic evolution

Predicting adaptive phenotypic evolution depends on invariable selection gradients and on the stability of the genetic covariances between the component traits of the multivariate phenotype. We describe the evolution of six traits of locomotion behavior and body size in the nematode Caenorhabditis elegans for 50 generations of adaptation to a novel environment. We show that the direction of adaptive multivariate phenotypic evolution can be predicted from the ancestral selection differentials, particularly when the traits were measured in the new environment. Interestingly, the evolution of individual traits does not always occur in the direction of selection, nor are trait responses to selection always homogeneous among replicate populations. These observations are explained because the phenotypic dimension with most of the ancestral standing genetic variation only partially aligns with the phenotypic dimension under directional selection. These findings validate selection theory and suggest that the direction of multivariate adaptive phenotypic evolution is predictable for tens of generations.

Embracing Complexity: Yeast Evolution Experiments Featuring Standing Genetic Variation

The yeast Saccharomyces cerevisiae has a long and esteemed history as a model system for laboratory selection experiments. The majority of yeast evolution experiments begin with an isogenic ancestor, impose selection as cells divide asexually, and track mutations that arise and accumulate over time. Within the last decade, the popularity of S. cerevisiae as a model system for exploring the evolution of standing genetic variation has grown considerably. As a facultatively sexual microbe, it is possible to initiate experiments with populations that harbor diversity and also to maintain that diversity by promoting sexual recombination as the experiment progresses. These experimental choices expand the scope of evolutionary hypotheses that can be tested with yeast. And, in this review, I argue that yeast is one of the best model systems for testing such hypotheses relevant to eukaryotic species. Here, I compile a list of yeast evolution experiments that involve standing genetic variation, initially and/or by implementing protocols that induce sexual recombination in evolving populations. I also provide an overview of experimental methods required to set up such an experiment and discuss the unique challenges that arise in this type of research. Throughout the article, I emphasize the best practices emerging from this small but growing niche of the literature.

Environmental radiation exposure at Chornobyl has not systematically affected the genomes or mutagen tolerance phenotypes of local worms

The 1986 disaster at the Chornobyl Nuclear Power Plant transformed the surrounding region into the most radioactive landscape known on the planet. Questions remain regarding whether this sudden environmental shift selected for species, or even individuals within a species, that are naturally more resistant to radiation exposure. We collected, cultured, and cryopreserved 298 wild nematodes isolates from areas varying in radioactivity within the Chornobyl Exclusion Zone. We sequenced and assembled genomes de novo for 20 Oschieus tipulae strains, analyzed their genomes for evidence of recent mutation acquisition in the field and saw no evidence of an association between mutation and radiation level at the sites of collection. Multigenerational exposure of each of these strains to several mutagens in the lab revealed that strains vary heritably in tolerance to each mutagen, but mutagen tolerance cannot be predicted based on the radiation levels at collection sites.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

Variation in mutational (co)variances

Because of pleiotropy, mutations affect the expression and inheritance of multiple traits and, together with selection, are expected to shape standing genetic covariances between traits and eventual phenotypic divergence between populations. It is therefore important to find if the M matrix, describing mutational variances of each trait and covariances between traits, varies between genotypes. We here estimate the M matrix for six locomotion behavior traits in lines of two genotypes of the nematode Caenorhabditis elegans that accumulated mutations in a nearly neutral manner for 250 generations. We find significant mutational variance along at least one phenotypic dimension of the M matrices, but neither their size nor their orientation had detectable differences between genotypes. The number of generations of mutation accumulation, or the number of MA lines measured, was likely insufficient to sample enough mutations and detect potentially small differences between the two M matrices. We then tested if the M matrices were similar to one G matrix describing the standing genetic (co)variances of a population derived by the hybridization of several genotypes, including the two measured for M, and domesticated to a lab-defined environment for 140 generations. We found that the M and G were different because the genetic covariances caused by mutational pleiotropy in the two genotypes are smaller than those caused by linkage disequilibrium in the lab population. We further show that M matrices differed in their alignment with the lab population G matrix. If generalized to other founder genotypes of the lab population, these observations indicate that selection does not shape the evolution of the M matrix for locomotion behavior in the short-term of a few tens to hundreds of generations and suggests that the hybridization of C. elegans genotypes allows selection on new phenotypic dimensions of locomotion behavior.

Ionizing radiation affects the demography and the evolution of Caenorhabditis elegans populations

Ionizing radiation can reduce survival, reproduction and affect development, and lead to the extinction of populations if their evolutionary response is insufficient. However, demographic and evolutionary studies on the effects of ionizing radiation are still scarce. Using an experimental evolution approach, we analyzed population growth rate and associated change in life history traits across generations in Caenorhabditis elegans populations exposed to 0, 1.4, and 50.0 mGy.h^-1 of ionizing radiation (gamma external irradiation). We found a higher population growth rate in the 1.4 mGy.h^-1 treatment and a lower in the 50.0 mGy.h^-1 treatment compared to the control. Realized fecundity was lower in both 1.4 and 50.0 mGy.h^-1 than control treatment. High irradiation levels decreased brood size from self-fertilized hermaphrodites, specifically early brood size. Finally, high irradiation levels decreased hatching success compared to the control condition. In reciprocal-transplant experiments, we found that life in low irradiation conditions led to the evolution of higher hatching success and late brood size. These changes could provide better tolerance against ionizing radiation, investing more in self-maintenance than in reproduction. These evolutionary changes were with some costs of adaptation. This study shows that ionizing radiation has both demographic and evolutionary consequences on populations.

Linkage mapping reveals loci that underlie differences in Caenorhabditis elegans growth

Growth rate and body size are complex traits that contribute to the fitness of organisms. The identification of loci that underlie differences in these traits provides insights into the genetic contributions to development. Leveraging Caenorhabditis elegans as a tractable metazoan model for quantitative genetics, we can identify genomic regions that underlie differences in growth. We measured postembryonic growth of the laboratory-adapted wild-type strain (N2) and a wild strain from Hawaii (CB4856) and found differences in body size. Using linkage mapping, we identified three distinct quantitative trait loci (QTL) on chromosomes IV, V, and X that are associated with variation in body growth. We further examined these growth-associated quantitative trait loci using chromosome substitution strains and near-isogenic lines and validated the chromosome X quantitative trait loci. In addition, we generated a list of candidate genes for the chromosome X quantitative trait loci. These genes could potentially contribute to differences in animal growth and should be evaluated in subsequent studies. Our work reveals the genetic architecture underlying animal growth variation and highlights the genetic complexity of growth in Caenorhabditis elegans natural populations.

Selection on modifiers of genetic architecture under migration load

Gene flow between populations adapting to differing local environmental conditions might be costly because individuals can disperse to habitats where their survival is low or because they can reproduce with locally maladapted individuals. The amount by which the mean relative population fitness is kept below one creates an opportunity for modifiers of the genetic architecture to spread due to selection. Prior work that separately considered modifiers changing dispersal, recombination rates, or altering dominance or epistasis, has typically focused on the direction of selection rather than its absolute magnitude. We here develop methods to determine the strength of selection on modifiers of the genetic architecture, including modifiers of the dispersal rate, in populations that have previously evolved local adaptation. We consider scenarios with up to five loci contributing to local adaptation and derive a new model for the deterministic spread of modifiers. We find that selection for modifiers of epistasis and dominance is stronger than selection for decreased recombination, and that selection for partial reductions in recombination are extremely weak, regardless of the number of loci contributing to local adaptation. The spread of modifiers that reduce dispersal depends on the number of loci, epistasis and extent of local adaptation in the ancestral population. We identify a novel effect, that modifiers of dominance are more strongly selected when they are unlinked to the locus that they modify. These findings help explain population differentiation and reproductive isolation and provide a benchmark to compare selection on modifiers under finite population sizes and demographic stochasticity.

When populations of a species are spread over different habitats the populations can adapt to their local conditions, provided dispersal between habitats is low enough. Natural selection allows the populations to maintain local adaptation, but dispersal and gene flow create a cost called the migration load. The migration load measures how much fitness is lost because of dispersal between different habitats, and also creates an opportunity for selection to act on the arrangement and interaction between genes that are involved in local adaptation. Modifier genes can spread in these linked populations and cause functional, local adaptation genes, to become more closely linked on a chromosome, or change the way that these genes are expressed so that the locally adapted gene copy becomes dominant. We modeled this process and found that selection on modifiers that create tighter linkage between locally adapted genes is generally weak, and modifiers that cause gene interactions are more strongly selected. Even after these gene interactions have begun to evolve, further selection for increased gene interaction is still strong. Our results show that populations are more likely to adapt to local conditions by evolving new gene interactions than by evolving tightly linked gene clusters.

Male frequency in Caenorhabditis elegans increases in response to chronic irradiation

Outcrossing can be advantageous in a changing environment because it promotes the purge of deleterious mutations and increases the genetic diversity within a population, which may improve population persistence and evolutionary potential. Some species may, therefore, switch their reproductive mode from inbreeding to outcrossing when under environmental stress. This switch may have consequences on the demographic dynamics and evolutionary trajectory of populations. For example, it may directly influence the sex ratio of a population. However, much remains to be discovered about the mechanisms and evolutionary implications of sex ratio changes in a population in response to environmental stress. Populations of the androdioecious nematode Caenorhabditis elegans, are composed of selfing hermaphrodites and rare males. Here, we investigate the changes in the sex ratio of C. elegans populations exposed to radioactive pollution for 60 days or around 20 generations. We experimentally exposed populations to three levels of ionizing radiation (i.e., 0, 1.4, and 50 mGy.h^-1). We then performed reciprocal transplant experiments to evaluate genetic divergence between populations submitted to different treatments. Finally, we used a mathematical model to examine the evolutionary mechanisms that could be responsible for the change in sex ratio. Our results showed an increase in male frequency in irradiated populations, and this effect increased with the dose rate. The model showed that an increase in male fertilization success or a decrease in hermaphrodite self-fertilization could explain this increase in the frequency of males. Moreover, males persisted in populations after transplant back into the control conditions. These results suggested selection favoring outcrossing under irradiation conditions. This study shows that ionizing radiation can sustainably alter the reproductive strategy of a population, likely impacting its long-term evolutionary history. This study highlights the need to evaluate the impact of pollutants on the reproductive strategies of populations when assessing the ecological risks.

Population bottleneck has only marginal effect on fitness evolution and its repeatability in dioecious Caenorhabditis elegans

The predictability of evolution is expected to depend on the relative contribution of deterministic and stochastic processes. This ratio is modulated by effective population size. Smaller effective populations harbor less genetic diversity and stochastic processes are generally expected to play a larger role, leading to less repeatable evolutionary trajectories. Empirical insight into the relationship between effective population size and repeatability is limited and focused mostly on asexual organisms. Here, we tested whether fitness evolution was less repeatable after a population bottleneck in obligately outcrossing populations of Caenorhabditis elegans. Replicated populations founded by 500, 50, or five individuals (no/moderate/strong bottleneck) were exposed to a novel environment with a different bacterial prey. As a proxy for fitness, population size was measured after one week of growth before and after 15 weeks of evolution. Surprisingly, we found no significant differences among treatments in their fitness evolution. Even though the strong bottleneck reduced the relative contribution of selection to fitness variation, this did not translate to a significant reduction in the repeatability of fitness evolution. Thus, although a bottleneck reduced the contribution of deterministic processes, we conclude that the predictability of evolution may not universally depend on effective population size, especially in sexual organisms.

Evaluating the power and limitations of genome-wide association studies in Caenorhabditis elegans

Quantitative genetics in Caenorhabditis elegans seeks to identify naturally segregating genetic variants that underlie complex traits. Genome-wide association studies scan the genome for individual genetic variants that are significantly correlated with phenotypic variation in a population, or quantitative trait loci. Genome-wide association studies are a popular choice for quantitative genetic analyses because the quantitative trait loci that are discovered segregate in natural populations. Despite numerous successful mapping experiments, the empirical performance of genome-wide association study has not, to date, been formally evaluated in C. elegans. We developed an open-source genome-wide association study pipeline called NemaScan and used a simulation-based approach to provide benchmarks of mapping performance in collections of wild C. elegans strains. Simulated trait heritability and complexity determined the spectrum of quantitative trait loci detected by genome-wide association studies. Power to detect smaller-effect quantitative trait loci increased with the number of strains sampled from the C. elegans Natural Diversity Resource. Population structure was a major driver of variation in mapping performance, with populations shaped by recent selection exhibiting significantly lower false discovery rates than populations composed of more divergent strains. We also recapitulated previous genome-wide association studies of experimentally validated quantitative trait variants. Our simulation-based evaluation of performance provides the community with critical context to pursue quantitative genetic studies using the C. elegans Natural Diversity Resource to elucidate the genetic basis of complex traits in C. elegans natural populations.

HGGA: hierarchical guided genome assembler

BACKGROUND

De novo genome assembly typically produces a set of contigs instead of the complete genome. Thus additional data such as genetic linkage maps, optical maps, or Hi-C data is needed to resolve the complete structure of the genome. Most of the previous work uses the additional data to order and orient contigs.

RESULTS

Here we introduce a framework to guide genome assembly with additional data. Our approach is based on clustering the reads, such that each read in each cluster originates from nearby positions in the genome according to the additional data. These sets are then assembled independently and the resulting contigs are further assembled in a hierarchical manner. We implemented our approach for genetic linkage maps in a tool called HGGA.

CONCLUSIONS

Our experiments on simulated and real Pacific Biosciences long reads and genetic linkage maps show that HGGA produces a more contiguous assembly with less contigs and from 1.2 to 9.8 times higher NGA50 or N50 than a plain assembly of the reads and 1.03 to 6.5 times higher NGA50 or N50 than a previous approach integrating genetic linkage maps with contig assembly. Furthermore, also the correctness of the assembly remains similar or improves as compared to an assembly using only the read data.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5–46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype–phenotype relationship.

The Genetic Architecture of a Congenital Heart Defect Is Related to Its Fitness Cost

In newborns, severe congenital heart defects are rarer than mild ones. This epidemiological relationship between heart defect severity and incidence lacks explanation. Here, an analysis of ~10,000 Nkx2-5^+/− mice from two inbred strain crosses illustrates the fundamental role of epistasis. Modifier genes raise or lower the risk of specific defects via pairwise (G×G_Nkx) and higher-order (G×G×G_Nkx) interactions with Nkx2-5. Their effect sizes correlate with the severity of a defect. The risk loci for mild, atrial septal defects exert predominantly small G×G_Nkx effects, while the loci for severe, atrioventricular septal defects exert large G×G_Nkx and G×G×G_Nkx effects. The loci for moderately severe ventricular septal defects have intermediate effects. Interestingly, G×G×G_Nkx effects are three times more likely to suppress risk when the genotypes at the first two loci are from the same rather than different parental inbred strains. This suggests the genetic coadaptation of interacting G×G×G_Nkx loci, a phenomenon that Dobzhansky first described in Drosophila. Thus, epistasis plays dual roles in the pathogenesis of congenital heart disease and the robustness of cardiac development. The empirical results suggest a relationship between the fitness cost and genetic architecture of a disease phenotype and a means for phenotypic robustness to have evolved.

Gene-level quantitative trait mapping in Caenorhabditis elegans

The Caenorhabditis elegans multiparental experimental evolution (CeMEE) panel is a collection of genome-sequenced, cryopreserved recombinant inbred lines useful for mapping the evolution and genetic basis of quantitative traits. We have expanded the resource with new lines and new populations, and here report the genotype and haplotype composition of CeMEE version 2, including a large set of putative de novo mutations, and updated additive and epistatic mapping simulations. Additive quantitative trait loci explaining 4% of trait variance are detected with >80% power, and the median detection interval approaches single-gene resolution on the highly recombinant chromosome arms. Although CeMEE populations are derived from a long-term evolution experiment, genetic structure is dominated by variation present in the ancestral population.

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.

Multiple haplotype reconstruction from allele frequency data

Because haplotype information is of widespread interest in biomedical applications, effort has been put into their reconstruction. Here, we propose an efficient method, called haploSep, that is able to accurately infer major haplotypes and their frequencies just from multiple samples of allele frequency data. Even the accuracy of experimentally obtained allele frequencies can be improved by re-estimating them from our reconstructed haplotypes. From a methodological point of view, we model our problem as a multivariate regression problem where both the design matrix and the coefficient matrix are unknown. Compared to other methods, haploSep is very fast, with linear computational complexity in the haplotype length. We illustrate our method on simulated and real data focusing on experimental evolution and microbial data.

Evolution: Drivers in a Traffic Jam

Genetic loci coding for a toxin and its antidote behave like selfish elements. Two new studies find an accumulation of such elements in one species.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

The Ancestral Caenorhabditis elegans Cuticle Suppresses rol-1

Genetic background commonly modifies the effects of mutations. We discovered that worms mutant for the canonical rol-1 gene, identified by Brenner in 1974, do not roll in the genetic background of the wild strain CB4856. Using linkage mapping, association analysis and gene editing, we determined that N2 carries an insertion in the collagen gene col-182 that acts as a recessive enhancer of rol-1 rolling. From population and comparative genomics, we infer the insertion is derived in N2 and related laboratory lines, likely arising during the domestication of Caenorhabditis elegans, and breaking a conserved protein. The ancestral version of col-182 also modifies the phenotypes of four other classical cuticle mutant alleles, and the effects of natural genetic variation on worm shape and locomotion. These results underscore the importance of genetic background and the serendipity of Brenner's choice of strain.

Phenotypic Responses to and Genetic Architecture of Sterility Following Exposure to Sub-Lethal Temperature During Development

Thermal tolerance range, based on temperatures that result in incapacitating effects, influences species’ distributions and has been used to predict species’ response to increasing temperature. Reproductive performance may also be negatively affected at less extreme temperatures, but such sublethal heat-induced sterility has been relatively ignored in studies addressing the potential effects of, and ability of species’ to respond to, predicted climate warming. The few studies examining the link between increased temperature and reproductive performance typically focus on adults, although effects can vary between life history stages. Here we assessed how sublethal heat stress during development impacted subsequent adult fertility and its plasticity, both of which can provide the raw material for evolutionary responses to increased temperature. We quantified phenotypic and genetic variation in fertility of Drosophila melanogaster reared at standardized densities in three temperatures (25, 27, and 29°C) from a set of lines of the Drosophila Genetic Reference Panel (DGRP). We found little phenotypic variation at the two lower temperatures with more variation at the highest temperature and for plasticity. Males were more affected than females. Despite reasonably large broad-sense heritabilities, a genome-wide association study found little evidence for additive genetic variance and no genetic variants were robustly linked with reproductive performance at specific temperatures or for phenotypic plasticity. We compared results on heat-induced male sterility with other DGRP results on relevant fitness traits measured after abiotic stress and found an association between male susceptibility to sterility and male lifespan reduction following oxidative stress. Our results suggest that sublethal stress during development has profound negative consequences on male adult reproduction, but despite phenotypic variation in a population for this response, there is limited evolutionary potential, either through adaptation to a specific developmental temperature or plasticity in response to developmental heat-induced sterility.

Genetic Variation in Caenorhabditis elegans Responses to Pathogenic Microbiota

The bacterivorous nematode Caenorhabditis elegans is an important model species for understanding genetic variation of complex traits. So far, most studies involve axenic laboratory settings using Escherichia coli as the sole bacterial species. Over the past decade, however, investigations into the genetic variation of responses to pathogenic microbiota have increasingly received attention. Quantitative genetic analyses have revealed detailed insight into loci, genetic variants, and pathways in C. elegans underlying interactions with bacteria, microsporidia, and viruses. As various quantitative genetic platforms and resources like C. elegans Natural Diversity Resource (CeNDR) and Worm Quantitative Trait Loci (WormQTL) have been developed, we anticipate that expanding C. elegans research along the lines of genetic variation will be a treasure trove for opening up new insights into genetic pathways and gene functionality of microbiota interactions.

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/.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

Population Selection and Sequencing of Caenorhabditis elegans Wild Isolates Identifies a Region on Chromosome III Affecting Starvation Resistance

To understand the genetic basis of complex traits, it is important to be able to efficiently phenotype many genetically distinct individuals. In the nematode Caenorhabditis elegans, individuals have been isolated from diverse populations around the globe and whole-genome sequenced. As a result, hundreds of wild strains with known genome sequences can be used for genome-wide association studies (GWAS). However, phenotypic analysis of these strains can be laborious, particularly for quantitative traits requiring multiple measurements per strain. Starvation resistance is likely a fitness-proximal trait for nematodes, and it is related to metabolic disease risk in humans. However, natural variation in C. elegans starvation resistance has not been systematically characterized, and precise measurement of the trait is time-intensive. Here, we developed a population-selection-and-sequencing-based approach to phenotype starvation resistance in a pool of 96 wild strains. We used restriction site-associated DNA sequencing (RAD-seq) to infer the frequency of each strain among survivors in a mixed culture over time during starvation. We used manual starvation survival assays to validate the trait data, confirming that strains that increased in frequency over time are starvation-resistant relative to strains that decreased in frequency. Further, we found that variation in starvation resistance is significantly associated with variation at a region on chromosome III. Using a near-isogenic line (NIL), we showed the importance of this genomic interval for starvation resistance. This study demonstrates the feasibility of using population selection and sequencing in an animal model for phenotypic analysis of quantitative traits, documents natural variation of starvation resistance in C. elegans, and identifies a genomic region that contributes to such variation.

Tightly linked antagonistic-effect loci underlie polygenic phenotypic variation in C. elegans

Recent work has provided strong empirical support for the classic polygenic model for trait variation. Population-based findings suggest that most regions of genome harbor variation affecting most traits. Here, we use the approach of experimental genetics to show that, indeed, most genomic regions carry variants with detectable effects on growth and reproduction in Caenorhabditis elegans populations sensitized by nickel stress. Nine of 15 adjacent intervals on the X chromosome, each encompassing ∼0.001 of the genome, have significant effects when tested individually in near-isogenic lines (NILs). These intervals have effects that are similar in magnitude to those of genome-wide significant loci that we mapped in a panel of recombinant inbred advanced intercross lines (RIAILs). If NIL-like effects were randomly distributed across the genome, the RIAILs would exhibit phenotypic variance that far exceeds the observed variance. However, the NIL intervals are arranged in a pattern that significantly reduces phenotypic variance relative to a random arrangement; adjacent intervals antagonize one another, cancelling each other's effects. Contrary to the expectation of small additive effects, our findings point to large-effect variants whose effects are masked by epistasis or linkage disequilibrium between alleles of opposing effect.

Partial Selfing Can Reduce Genetic Loads While Maintaining Diversity During Experimental Evolution

Partial selfing, whereby self- and cross- fertilization occur in populations at intermediate frequencies, is generally thought to be evolutionarily unstable. Yet, it is found in natural populations. This could be explained if populations with partial selfing are able to reduce genetic loads and the possibility for inbreeding depression while keeping genetic diversity that may be important for future adaptation. To address this hypothesis, we compare the experimental evolution of Caenorhabditis elegans populations under partial selfing, exclusive selfing or predominant outcrossing, while they adapt to osmotically challenging conditions. We find that the ancestral genetic load, as measured by the risk of extinction upon inbreeding by selfing, is maintained as long as outcrossing is the main reproductive mode, but becomes reduced otherwise. Analysis of genome-wide single-nucleotide polymorphisms (SNPs) during experimental evolution and among the inbred lines that survived enforced inbreeding indicates that populations with predominant outcrossing or partial selfing maintained more genetic diversity than expected with neutrality or purifying selection. We discuss the conditions under which this could be explained by the presence of recessive deleterious alleles and/or overdominant loci. Taken together, our observations suggest that populations evolving under partial selfing can gain some of the benefits of eliminating unlinked deleterious recessive alleles and also the benefits of maintaining genetic diversity at partially dominant or overdominant loci that become associated due to variance of inbreeding levels.

How differing modes of non-genetic inheritance affect population viability in fluctuating environments

Different modes of non-genetic inheritance are expected to affect population persistence in fluctuating environments. We here analyse Caenorhabditis elegans density-independent per capita growth rate time series on 36 populations experiencing six controlled sequences of challenging oxygen level fluctuations across 60 generations, and parameterise competing models of non-genetic inheritance in order to explain observed dynamics. Our analysis shows that phenotypic plasticity and anticipatory maternal effects are sufficient to explain growth rate dynamics, but that a carryover model where 'epigenetic' memory is imperfectly transmitted and might be reset at each generation is a better fit to the data. We further find that this epigenetic memory is asymmetric since it is kept for longer when populations are exposed to the more challenging environment. Our analysis suggests that population persistence in fluctuating environments depends on the non-genetic inheritance of phenotypes whose expression is regulated across multiple generations.

Fast genetic mapping of complex traits in C. elegans using millions of individuals in bulk

Genetic studies of complex traits in animals have been hindered by the need to generate, maintain, and phenotype large panels of recombinant lines. We developed a new method, C. elegans eXtreme Quantitative Trait Locus (ceX-QTL) mapping, that overcomes this obstacle via bulk selection on millions of unique recombinant individuals. We use ceX-QTL to map a drug resistance locus with high resolution. We also map differences in gene expression in live worms and discovered that mutations in the co-chaperone sti-1 upregulate the transcription of HSP-90. Lastly, we use ceX-QTL to map loci that influence fitness genome-wide confirming previously reported causal variants and uncovering new fitness loci. ceX-QTL is fast, powerful and cost-effective, and will accelerate the study of complex traits in animals.

Determinants of QTL Mapping Power in the Realized Collaborative Cross

The Collaborative Cross (CC) is a mouse genetic reference population whose range of applications includes quantitative trait loci (QTL) mapping. The design of a CC QTL mapping study involves multiple decisions, including which and how many strains to use, and how many replicates per strain to phenotype, all viewed within the context of hypothesized QTL architecture. Until now, these decisions have been informed largely by early power analyses that were based on simulated, hypothetical CC genomes. Now that more than 50 CC strains are available and more than 70 CC genomes have been observed, it is possible to characterize power based on realized CC genomes. We report power analyses from extensive simulations and examine several key considerations: 1) the number of strains and biological replicates, 2) the QTL effect size, 3) the presence of population structure, and 4) the distribution of functionally distinct alleles among the founder strains at the QTL. We also provide general power estimates to aide in the design of future experiments. All analyses were conducted with our R package, SPARCC (Simulated Power Analysis in the Realized Collaborative Cross), developed for performing either large scale power analyses or those tailored to particular CC experiments.

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.

The genomic basis of Red Queen dynamics during rapid reciprocal host-pathogen coevolution

Pathogens are omnipresent and by definition detrimental to their hosts. Pathogens thus exert high selection on their hosts, which, if adapting, can exert similar levels of selection on the pathogen, resulting in ongoing cycles of reciprocal adaptation between the antagonists. Such coevolutionary interactions have a central influence on the evolution of organisms. Surprisingly, we still know little about the exact selection dynamics and the genome regions involved. Our study uses a controlled experimental approach with an animal host to dissect coevolutionary selection. We find that distinct selective processes underlie rapid coadaptation in the two antagonists, including antagonistic frequency-dependent selection on toxin gene copy number in the pathogen, while the host response is likely influenced by changes in multiple genome regions.

Red Queen dynamics, involving coevolutionary interactions between species, are ubiquitous, shaping the evolution of diverse biological systems. To date, information on the underlying selection dynamics and the involved genome regions is mainly available for bacteria–phage systems or only one of the antagonists of a eukaryotic host–pathogen interaction. We add to our understanding of these important coevolutionary interactions using an experimental host–pathogen model, which includes the nematode Caenorhabditis elegans and its pathogen Bacillus thuringiensis. We combined experimental evolution with time-shift experiments, in which a focal host or pathogen is tested against a coevolved antagonist from the past, present, or future, followed by genomic analysis. We show that (i) coevolution occurs rapidly within few generations, (ii) temporal coadaptation at the phenotypic level is found in parallel across replicate populations, consistent with antagonistic frequency-dependent selection, (iii) genomic changes in the pathogen match the phenotypic pattern and include copy number variations of a toxin-encoding plasmid, and (iv) host genomic changes do not match the phenotypic pattern and likely involve selective responses at more than one locus. By exploring the dynamics of coevolution at the phenotypic and genomic level for both host and pathogen simultaneously, our findings demonstrate a more complex model of the Red Queen, consisting of distinct selective processes acting on the two antagonists during rapid and reciprocal coadaptation.

R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely-used R/qtl software package to include multiparental populations, better handles modern high-dimensional data....

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely used R/qtl software package to include multiparent populations derived from more than two founder strains, such as the Collaborative Cross and Diversity Outbred mice, heterogeneous stocks, and MAGIC plant populations. R/qtl2 is designed to handle modern high-density genotyping data and high-dimensional molecular phenotypes, including gene expression and proteomics. R/qtl2 includes the ability to perform genome scans using a linear mixed model to account for population structure, and also includes features to impute SNPs based on founder strain genomes and to carry out association mapping. The R/qtl2 software provides all of the basic features needed for QTL mapping, including graphical displays and summary reports, and it can be extended through the creation of add-on packages. R/qtl2, which is free and open source software written in the R and C++ programming languages, comes with a test framework.

Slower environmental change hinders adaptation from standing genetic variation

Evolutionary responses to environmental change depend on the time available for adaptation before environmental degradation leads to extinction. Explicit tests of this relationship are limited to microbes where adaptation usually depends on the sequential fixation of de novo mutations, excluding standing variation for genotype-by-environment fitness interactions that should be key for most natural species. For natural species evolving from standing genetic variation, adaptation at slower rates of environmental change may be impeded since the best genotypes at the most extreme environments can be lost during evolution due to genetic drift or founder effects. To address this hypothesis, we perform experimental evolution with self-fertilizing populations of the nematode Caenorhabditis elegans and develop an inference model to describe natural selection on extant genotypes under environmental change. Under a sudden environmental change, we find that selection rapidly increases the frequency of genotypes with high fitness in the most extreme environment. In contrast, under a gradual environmental change selection first favors genotypes that are worse at the most extreme environment. We demonstrate with a second set of evolution experiments that, as a consequence of slower environmental change and thus longer periods to reach the most extreme environments, genetic drift and founder effects can lead to the loss of the most beneficial genotypes. We further find that maintenance of standing genetic variation can retard the fixation of the best genotypes in the most extreme environment because of interference between them. Taken together, these results show that slower environmental change can hamper adaptation from standing genetic variation and they support theoretical models indicating that standing variation for genotype-by-environment fitness interactions critically alters the pace and outcome of adaptation under environmental change.

Analysis of Epistasis in Natural Traits Using Model Organisms

The ability to detect and understand epistasis in natural populations is important for understanding how biological traits are influenced by genetic variation. However, identification and characterization of epistasis in natural populations remains difficult due to statistical issues that arise as a result of multiple comparisons, and the fact that most genetic variants segregate at low allele frequencies. In this review, we discuss how model organisms may be used to manipulate genotypic combinations to power the detection of epistasis as well as test interactions between specific genes. Findings from a number of species indicate that statistical epistasis is pervasive between natural genetic variants. However, the properties of experimental systems that enable analysis of epistasis also constrain extrapolation of these results back into natural populations.

Changes to social feeding behaviors are not sufficient for fitness gains of the Caenorhabditis elegans N2 reference strain

The standard reference Caenorhabditis elegans strain, N2, has evolved marked behavioral changes in social feeding behavior since its isolation from the wild. We show that the causal, laboratory-derived mutations in two genes, npr-1 and glb-5, confer large fitness advantages in standard laboratory conditions. Using environmental manipulations that suppress social/solitary behavior differences, we show the fitness advantages of the derived alleles remained unchanged, suggesting selection on these alleles acted through pleiotropic traits. Transcriptomics, developmental timing, and food consumption assays showed that N2 animals mature faster, produce more sperm, and consume more food than a strain containing ancestral alleles of these genes regardless of behavioral strategies. Our data suggest that the pleiotropic effects of glb-5 and npr-1 are a consequence of changes to O₂ -sensing neurons that regulate both aerotaxis and energy homeostasis. Our results demonstrate how pleiotropy can lead to profound behavioral changes in a popular laboratory model.

Why do humans walk on two feet? And what makes us smarter than our ape ancestors? The answers to these questions, and countless others about the particular traits of any number of species, is often said to be natural selection – a process where genes that ensure the survival of a species are favored of others. But it is not always the answer. Other evolutionary forces, such as random changes to the frequency of certain gene variants, restrictions on the development of a certain trait and pleiotropy (where one gene influences other, seemingly unrelated traits) can also cause differences between species.

Designing experiments to test whether a trait difference is due to natural selection or other factors is notoriously difficult. However, the humble nematode worm, Caenorhabditis elegans, has proven to be particularly useful in this respect. One subtype or strain of C. elegans with certain changes to its genes is used internationally as a ‘reference strain’, to ensure results between labs are comparable. This strain, N2, has been bred in the laboratory for hundreds of generations, isolated from its wild counterparts. N2 shows several differences in behavior from the wildtype, including its feeding habits. Wild C. elegans tend to feed together socially, whereas N2 prefers to feed alone.

In 1998 and 2009, researchers – including some involved in the current study – have identified the genetic modifications responsible for this change in behavior. Now, Zhao et al. set out to determine whether this was due to natural selection, and if so, was there a benefit to solitary feeding in laboratory conditions that was driving this genetic change? Zhao et al. found that the genetic changes in the N2 strain gave the worms a considerable advantage in the artificial environment. However, experiments to modify the conditions the animals grew in revealed that the solitary feeding habits were not necessary for the fitness advantage. In other words, the changes in feeding habits were a symptom of the genetic changes that gave N2 a selective advantage, but they were not the cause. In other words, the changes in feeding behavior were not a result of natural selection, but rather of pleiotropy.

The findings highlight that not every change in a trait is down to natural selection and must therefore be put to the test. With declining costs of DNA sequencing, researchers can now easily identify genes and regions of DNA that are likely to be under selection. However, they must be careful before leaping to the conclusion that behavioral differences linked to genetic changes are adaptive. In addition, the findings show that the laboratories relying on N2 as a model organism should be aware that the strain has evolved fundamental differences in its brain connections compared with the wildtype.