Chromosome-scale selective sweeps shape Caenorhabditis elegans genomic diversity

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains.

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.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

A Novel Gene Underlies Bleomycin-Response Variation in Caenorhabditis elegans

Bleomycin is a powerful chemotherapeutic drug used to treat a variety of cancers. However, individual patients vary in their responses to bleomycin. The identification of genetic differences that underlie this response variation could improve treatment outcomes by tailoring bleomycin dosages to each patient. We used the model organism Caenorhabditis elegans to identify genetic determinants of bleomycin-response differences by performing linkage mapping on recombinants derived from a cross between the laboratory strain (N2) and a wild strain (CB4856). This approach identified a small genomic region on chromosome V that underlies bleomycin-response variation. Using near-isogenic lines, and strains with CRISPR-Cas9 mediated deletions and allele replacements, we discovered that a novel nematode-specific gene (scb-1) is required for bleomycin resistance. Although the mechanism by which this gene causes variation in bleomycin responses is unknown, we suggest that a rare variant present in the CB4856 strain might cause differences in the potential stress-response function of scb-1 between the N2 and CB4856 strains, thereby leading to differences in bleomycin resistance.

Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

Natural variation in C. elegans arsenic toxicity is explained by differences in branched chain amino acid metabolism

We find that variation in the dbt-1 gene underlies natural differences in Caenorhabditis elegans responses to the toxin arsenic. This gene encodes the E2 subunit of the branched-chain α-keto acid dehydrogenase (BCKDH) complex, a core component of branched-chain amino acid (BCAA) metabolism. We causally linked a non-synonymous variant in the conserved lipoyl domain of DBT-1 to differential arsenic responses. Using targeted metabolomics and chemical supplementation, we demonstrate that differences in responses to arsenic are caused by variation in iso-branched chain fatty acids. Additionally, we show that levels of branched chain fatty acids in human cells are perturbed by arsenic treatment. This finding has broad implications for arsenic toxicity and for arsenic-focused chemotherapeutics across human populations. Our study implicates the BCKDH complex and BCAA metabolism in arsenic responses, demonstrating the power of C. elegans natural genetic diversity to identify novel mechanisms by which environmental toxins affect organismal physiology.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

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.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

The population genetics of parasitic nematodes of wild animals

Parasitic nematodes are highly diverse and common, infecting virtually all animal species, and the importance of their roles in natural ecosystems is increasingly becoming apparent. How genes flow within and among populations of these parasites - their population genetics - has profound implications for the epidemiology of host infection and disease, and for the response of parasite populations to selection pressures. The population genetics of nematode parasites of wild animals may have consequences for host conservation, or influence the risk of zoonotic disease. Host movement has long been recognised as an important determinant of parasitic nematode population genetic structure, and recent research has also highlighted the importance of nematode life histories, environmental conditions, and other aspects of host ecology. Commonly, factors influencing parasitic nematode population genetics have been studied in isolation, such that an integrated view of the drivers of population genetic structure of parasitic nematodes is still lacking. Here, we seek to provide a comprehensive, broad, and integrative picture of these factors in parasitic nematodes of wild animals that will be a useful resource for investigators studying non-model parasitic nematodes in natural ecosystems. Increasingly, new methods of analysing the population genetics of nematodes are becoming available, and we consider the opportunities that these afford in resolving hitherto inaccessible questions of the population genetics of these important animals.

Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles

Benzimidazoles (BZ) are essential components of the limited chemotherapeutic arsenal available to control the global burden of parasitic nematodes. The emerging threat of BZ resistance among multiple nematode species necessitates the development of novel strategies to identify genetic and molecular mechanisms underlying this resistance. All detection of parasitic helminth resistance to BZ is focused on the genotyping of three variant sites in the orthologs of the β-tubulin gene found to confer resistance in the free-living nematode Caenorhabditis elegans. Because of the limitations of laboratory and field experiments in parasitic nematodes, it is difficult to look beyond these three sites to identify additional mechanisms that might contribute to BZ resistance in the field. Here, we took an unbiased genome-wide mapping approach in the free-living nematode species C. elegans to identify the genetic underpinnings of natural resistance to the commonly used BZ, albendazole (ABZ). We found a wide range of natural variation in ABZ resistance in natural C. elegans populations. In agreement with known mechanisms of BZ resistance in parasites, we found that a majority of the variation in ABZ resistance among wild C. elegans strains is caused by variation in the β-tubulin gene ben-1. This result shows empirically that resistance to ABZ naturally exists and segregates within the C. elegans population, suggesting that selection in natural niches could enrich for resistant alleles. We identified 25 distinct ben-1 alleles that are segregating at low frequencies within the C. elegans population, including many novel molecular variants. Population genetic analyses indicate that ben-1 variation arose multiple times during the evolutionary history of C. elegans and provide evidence that these alleles likely occurred recently because of local selective pressures. Additionally, we find purifying selection at all five β-tubulin genes, despite predicted loss-of-function variants in ben-1, indicating that BZ resistance in natural niches is a stronger selective pressure than loss of one β-tubulin gene. Furthermore, we used genome-editing to show that the most common parasitic nematode β-tubulin allele that confers BZ resistance, F200Y, confers resistance in C. elegans. Importantly, we identified a novel genomic region that is correlated with ABZ resistance in the C. elegans population but independent of ben-1 and the other β-tubulin loci, suggesting that there are multiple mechanisms underlying BZ resistance. Taken together, our results establish a population-level resource of nematode natural diversity as an important model for the study of mechanisms that give rise to BZ resistance.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

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

Acute Effects of Drugs on Caenorhabditis elegans Movement Reveal Complex Responses and Plasticity

Many drugs act very rapidly - they can turn on or off their targets within minutes in a whole animal. What are the acute effects of drug treatment and how does an animal respond to these? We developed a simple assay to measure the acute effects of drugs on C. elegans movement and examined the effects of a range of compounds including neuroactive drugs, toxins, environmental stresses and novel compounds on worm movement over a time period of 3 hr. We found a wide variety of acute responses. Many compounds cause rapid paralysis which may be permanent or followed by one or more recovery phases. The recoveries are not the result of some generic stress response but are specific to the drug e.g., recovery from paralysis due to a neuroactive drug requires neurotransmitter pathways whereas recovery from a metabolic inhibitor requires metabolic changes. Finally, we also find that acute responses can vary greatly across development and that there is extensive natural variation in acute responses. In summary, acute responses are sensitive probes of the ability of biological networks to respond to drug treatment and these responses can reveal the action of unexplored pathways.

Field studies reveal a close relative of C. elegans thrives in the fresh figs of Ficus septica and disperses on its Ceratosolen pollinating wasps

BACKGROUND

Biotic interactions are ubiquitous and require information from ecology, evolutionary biology, and functional genetics in order to be understood. However, study systems that are amenable to investigations across such disparate fields are rare. Figs and fig wasps are a classic system for ecology and evolutionary biology with poor functional genetics; Caenorhabditis elegans is a classic system for functional genetics with poor ecology. In order to help bridge these disciplines, here we describe the natural history of a close relative of C. elegans, Caenorhabditis inopinata, that is associated with the fig Ficus septica and its pollinating Ceratosolen wasps.

RESULTS

To understand the natural context of fig-associated Caenorhabditis, fresh F. septica figs from four Okinawan islands were sampled, dissected, and observed under microscopy. C. inopinata was found in all islands where F. septica figs were found. C.i nopinata was routinely found in the fig interior and almost never observed on the outside surface. C. inopinata was only found in pollinated figs, and C. inopinata was more likely to be observed in figs with more foundress pollinating wasps. Actively reproducing C. inopinata dominated early phase figs, whereas late phase figs with emerging wasp progeny harbored C. inopinata dauer larvae. Additionally, C. inopinata was observed dismounting from Ceratosolen pollinating wasps that were placed on agar plates. C. inopinata was not found on non-pollinating, parasitic Philotrypesis wasps. Finally, C. inopinata was only observed in F. septica figs among five Okinawan Ficus species sampled.

CONCLUSION

These are the first detailed field observations of C. inopinata, and they suggest a natural history where this species proliferates in early phase F. septica figs and disperses from late phase figs on Ceratosolen pollinating fig wasps. While consistent with other examples of nematode diversification in the fig microcosm, the fig and wasp host specificity of C. inopinata is highly divergent from the life histories of its close relatives and frames hypotheses for future investigations. This natural co-occurrence of the fig/fig wasp and C. inopinata study systems sets the stage for an integrated research program that can help to explain the evolution of interspecific interactions.

Natural Genetic Variation in a Multigenerational Phenotype in C. elegans

Although heredity mostly relies on the transmission of DNA sequence, additional molecular and cellular features are heritable across several generations. In the nematode Caenorhabditis elegans, insights into such unconventional inheritance result from two lines of work. First, the mortal germline (Mrt) phenotype was defined as a multigenerational phenotype whereby a selfing lineage becomes sterile after several generations, implying multigenerational memory [1, 2]. Second, certain RNAi effects are heritable over several generations in the absence of the initial trigger [3-5]. Both lines of work converged when the subset of Mrt mutants that are heat sensitive were found to closely correspond to mutants defective in the RNAi-inheritance machinery, including histone modifiers [6-9]. Here, we report the surprising finding that several C. elegans wild isolates display a heat-sensitive mortal germline phenotype in laboratory conditions: upon chronic exposure to higher temperatures, such as 25°C, lines reproducibly become sterile after several generations. This phenomenon is reversible, as it can be suppressed by temperature alternations at each generation, suggesting a non-genetic basis for the sterility. We tested whether natural variation in the temperature-induced Mrt phenotype was of genetic nature by building recombinant inbred lines between the isolates MY10 (Mrt) and JU1395 (non-Mrt). Using bulk segregant analysis, we detected two quantitative trait loci. After further recombinant mapping and genome editing, we identified the major causal locus as a polymorphism in the set-24 gene, encoding a SET- and SPK-domain protein. We conclude that C. elegans natural populations may harbor natural genetic variation in epigenetic inheritance phenomena.

Natural genetic variation in C. elegans identified genomic loci controlling metabolite levels

Metabolic homeostasis is sustained by complex biological networks that respond to nutrient availability. Genetic and environmental factors may disrupt this equilibrium, leading to metabolic disorders, including obesity and type 2 diabetes. To identify the genetic factors controlling metabolism, we performed quantitative genetic analysis using a population of 199 recombinant inbred lines (RILs) in the nematode Caenorhabditis elegans We focused on the genomic regions that control metabolite levels by measuring fatty acid (FA) and amino acid (AA) composition in the RILs using targeted metabolomics. The genetically diverse RILs showed a large variation in their FA and AA levels with a heritability ranging from 32% to 82%. We detected strongly co-correlated metabolite clusters and 36 significant metabolite quantitative trait loci (mQTL). We focused on mQTL displaying highly significant linkage and heritability, including an mQTL for the FA C14:1 on Chromosome I, and another mQTL for the FA C18:2 on Chromosome IV. Using introgression lines (ILs), we were able to narrow down both mQTL to a 1.4-Mbp and a 3.6-Mbp region, respectively. RNAi-based screening focusing on the Chromosome I mQTL identified several candidate genes for the C14:1 mQTL, including lagr-1, Y87G2A.2, nhr-265, nhr-276, and nhr-81 Overall, this systems approach provides us with a powerful platform to study the genetic basis of C. elegans metabolism. Furthermore, it allows us to investigate interventions such as nutrients and stresses that maintain or disturb the regulatory network controlling metabolic homeostasis, and identify gene-by-environment interactions.

Human demographic history has amplified the effects of background selection across the genome

Natural populations often grow, shrink, and migrate over time. Such demographic processes can affect genome-wide levels of genetic diversity. Additionally, genetic variation in functional regions of the genome can be altered by natural selection, which drives adaptive mutations to higher frequencies or purges deleterious ones. Such selective processes affect not only the sites directly under selection but also nearby neutral variation through genetic linkage via processes referred to as genetic hitchhiking in the context of positive selection and background selection (BGS) in the context of purifying selection. While there is extensive literature examining the consequences of selection at linked sites at demographic equilibrium, less is known about how non-equilibrium demographic processes influence the effects of hitchhiking and BGS. Utilizing a global sample of human whole-genome sequences from the Thousand Genomes Project and extensive simulations, we investigate how non-equilibrium demographic processes magnify and dampen the consequences of selection at linked sites across the human genome. When binning the genome by inferred strength of BGS, we observe that, compared to Africans, non-African populations have experienced larger proportional decreases in neutral genetic diversity in strong BGS regions. We replicate these findings in admixed populations by showing that non-African ancestral components of the genome have also been affected more severely in these regions. We attribute these differences to the strong, sustained/recurrent population bottlenecks that non-Africans experienced as they migrated out of Africa and throughout the globe. Furthermore, we observe a strong correlation between FST and the inferred strength of BGS, suggesting a stronger rate of genetic drift. Forward simulations of human demographic history with a model of BGS support these observations. Our results show that non-equilibrium demography significantly alters the consequences of selection at linked sites and support the need for more work investigating the dynamic process of multiple evolutionary forces operating in concert.

Single-Molecule Sequencing Reveals the Chromosome-Scale Genomic Architecture of the Nematode Model Organism Pristionchus pacificus

The nematode Pristionchus pacificus is an established model for integrative evolutionary biology and comparative studies with Caenorhabditis elegans. While an existing genome draft facilitated the identification of several genes controlling various developmental processes, its high degree of fragmentation complicated virtually all genomic analyses. Here, we present a de novo genome assembly from single-molecule, long-read sequencing data consisting of 135 P. pacificus contigs. When combined with a genetic linkage map, 99% of the assembly could be ordered and oriented into six chromosomes. This allowed us to robustly characterize chromosomal patterns of gene density, repeat content, nucleotide diversity, linkage disequilibrium, and macrosynteny in P. pacificus. Despite widespread conservation of synteny between P. pacificus and C. elegans, we identified one major translocation from an autosome to the sex chromosome in the lineage leading to C. elegans. This highlights the potential of the chromosome-scale assembly for future genomic studies of P. pacificus.

Linking Genomic and Metabolomic Natural Variation Uncovers Nematode Pheromone Biosynthesis

In the nematodes Caenorhabditis elegans and Pristionchus pacificus, a modular library of small molecules control behavior, lifespan, and development. However, little is known about the final steps of their biosynthesis, in which diverse building blocks from primary metabolism are attached to glycosides of the dideoxysugar ascarylose, the ascarosides. We combine metabolomic analysis of natural isolates of P. pacificus with genome-wide association mapping to identify a putative carboxylesterase, Ppa-uar-1, that is required for attachment of a pyrimidine-derived moiety in the biosynthesis of ubas#1, a major dauer pheromone component. Comparative metabolomic analysis of wild-type and Ppa-uar-1 mutants showed that Ppa-uar-1 is required specifically for the biosynthesis of ubas#1 and related metabolites. Heterologous expression of Ppa-UAR-1 in C. elegans yielded a non-endogenous ascaroside, whose structure confirmed that Ppa-uar-1 is involved in modification of a specific position in ascarosides. Our study demonstrates the utility of natural variation-based approaches for uncovering biosynthetic pathways.

Early experiences mediate distinct adult gene expression and reproductive programs in Caenorhabditis elegans

Environmental stress during early development in animals can have profound effects on adult phenotypes via programmed changes in gene expression. Using the nematode C. elegans, we demonstrated previously that adults retain a cellular memory of their developmental experience that is manifested by differences in gene expression and life history traits; however, the sophistication of this system in response to different environmental stresses, and how it dictates phenotypic plasticity in adults that contribute to increased fitness in response to distinct environmental challenges, was unknown. Using transcriptional profiling, we show here that C. elegans adults indeed retain distinct cellular memories of different environmental conditions. We identified approximately 500 genes in adults that entered dauer due to starvation that exhibit significant opposite (“seesaw”) transcriptional phenotypes compared to adults that entered dauer due to crowding, and are distinct from animals that bypassed dauer. Moreover, we show that two-thirds of the genes in the genome experience a 2-fold or greater seesaw trend in gene expression, and based upon the direction of change, are enriched in large, tightly linked regions on different chromosomes. Importantly, these transcriptional programs correspond to significant changes in brood size depending on the experienced stress. In addition, we demonstrate that while the observed seesaw gene expression changes occur in both somatic and germline tissue, only starvation-induced changes require a functional GLP-4 protein necessary for germline development, and both programs require the Argonaute CSR-1. Thus, our results suggest that signaling between the soma and the germ line can generate phenotypic plasticity as a result of early environmental experience, and likely contribute to increased fitness in adverse conditions and the evolution of the C. elegans genome.

Environmental stress during early development in animals can have profound effects on adult behavior and physiology due to programmed changes in gene expression. However, whether different stresses result in distinct changes in traits that allow stressed animals to better survive and reproduce in future adverse conditions is largely unknown. Using the animal model system, C. elegans, we show that adults that experienced starvation exhibit opposite (“seesaw”) genome-wide gene expression changes compared to adults that experienced crowding, and are distinct from animals that experienced favorable conditions. Genes that are similarly up- or downregulated due to either starvation or crowding are located in clusters on the same chromosomes. Importantly, these gene expression changes of differently-stressed animals result in corresponding changes in progeny number, a life history trait of evolutionary significance. These distinct gene expression programs require different signaling pathways that communicate across somatic and germline tissue types. Thus, different environmental stresses experienced early in development induce distinct signaling mechanisms to result in changes in gene expression and reproduction in adults, and likely contribute to increased survival in future adverse conditions.

Rictor/TORC2 mediates gut-to-brain signaling in the regulation of phenotypic plasticity in C. elegans

Animals integrate external cues with information about internal conditions such as metabolic state to execute the appropriate behavioral and developmental decisions. Information about food quality and quantity is assessed by the intestine and transmitted to modulate neuronal functions via mechanisms that are not fully understood. The conserved Target of Rapamycin complex 2 (TORC2) controls multiple processes in response to cellular stressors and growth factors. Here we show that TORC2 coordinates larval development and adult behaviors in response to environmental cues and feeding state in the bacterivorous nematode C. elegans. During development, pheromone, bacterial food, and temperature regulate expression of the daf-7 TGF-β and daf-28 insulin-like peptide in sensory neurons to promote a binary decision between reproductive growth and entry into the alternate dauer larval stage. We find that TORC2 acts in the intestine to regulate neuronal expression of both daf-7 and daf-28, which together reflect bacterial-diet dependent feeding status, thus providing a mechanism for integration of food signals with external cues in the regulation of neuroendocrine gene expression. In the adult, TORC2 similarly acts in the intestine to modulate food-regulated foraging behaviors via a PDF-2/PDFR-1 neuropeptide signaling-dependent pathway. We also demonstrate that genetic variation affects food-dependent larval and adult phenotypes, and identify quantitative trait loci (QTL) associated with these traits. Together, these results suggest that TORC2 acts as a hub for communication of feeding state information from the gut to the brain, thereby contributing to modulation of neuronal function by internal state.

Decision-making in all animals, including humans, involves weighing available information about the external environment as well as the animals’ internal conditions. Information about the environment is obtained via the sensory nervous system, whereas internal state can be assessed via cues such as levels of hormones or nutrients. How multiple external and internal inputs are processed in the nervous system to drive behavior or development is not fully understood. In this study, we examine how the nematode C. elegans integrates dietary information received by the gut with environmental signals to alter nervous system function. We have found that a signaling complex, called TORC2, acts in the gut to relay nutrition signals to alter hormonal signaling by the nervous system in C. elegans. Altered neuronal signaling in turn affects a food-dependent binary developmental decision in larvae, as well as food-dependent foraging behaviors in adults. Our results provide a mechanism by which animals prioritize specific signals such as feeding status to appropriately alter their development and/or behavior.

The Local Coexistence Pattern of Selfing Genotypes in Caenorhabditis elegans Natural Metapopulations

To study the interplay of rare outcrossing and metapopulation structure, we focus on the nematode Caenorhabditis elegans Its remarkably low outcrossing rate is at the extreme end of the spectrum for facultative selfing organisms. At the demographic level, C. elegans natural populations undergo boom and bust dynamics on ephemeral resources, with the dauer diapause larva acting as the dispersal form. Here we investigate the small-scale genetic structure of C. elegans populations in two localities over several years, using 2b restriction-associated DNA sequencing of nearly 1000 individuals. We find a remarkably small number of genome-wide haplotypes, almost exclusively in the homozygous state, confirming the low effective outcrossing rate. Most strikingly, the major haplotypes in a locality remain intact and do not effectively recombine over several years. From the spatial pattern of diversity, we estimate that each subpopulation or deme is seeded by a mean of 3-10 immigrating individuals. Populations are thus formed by clones that compete at two levels, within a subpopulation and at the metapopulation level. We test for the presence of local phenotypic variation in pathogen resistance and dauer larva nictation, which could possibly explain the maintenance of different genotypes by heterogeneous selection in different local environments or lifecycles. This study is the first to address the local spatiotemporal genetic structure of C. elegans on feeding substrates. We conclude that these animals coexist as competing homozygous clones at the smallest population scale as well as in the metapopulation.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

Evolution of mitotic spindle behavior during the first asymmetric embryonic division of nematodes

Asymmetric cell division is essential to generate cellular diversity. In many animal cells, the cleavage plane lies perpendicular to the mitotic spindle, and it is the spindle positioning that dictates the size of the daughter cells. Although some properties of spindle positioning are conserved between distantly related model species and different cell types, little is known of the evolutionary robustness of the mechanisms underlying this event. We recorded the first embryonic division of 42 species of nematodes closely related to Caenorhabditis elegans, which is an excellent model system to study the biophysical properties of asymmetric spindle positioning. Our recordings, corresponding to 128 strains from 27 Caenorhabditis and 15 non-Caenorhabditis species (accessible at http://www.ens-lyon.fr/LBMC/NematodeCell/videos/), constitute a powerful collection of subcellular phenotypes to study the evolution of various cellular processes across species. In the present work, we analyzed our collection to the study of asymmetric spindle positioning. Although all the strains underwent an asymmetric first cell division, they exhibited large intra- and inter-species variations in the degree of cell asymmetry and in several parameters controlling spindle movement, including spindle oscillation, elongation, and displacement. Notably, these parameters changed frequently during evolution with no apparent directionality in the species phylogeny, with the exception of spindle transverse oscillations, which were an evolutionary innovation at the base of the Caenorhabditis genus. These changes were also unrelated to evolutionary variations in embryo size. Importantly, spindle elongation, displacement, and oscillation each evolved independently. This finding contrasts starkly with expectations based on C. elegans studies and reveals previously unrecognized evolutionary changes in spindle mechanics. Collectively, these data demonstrate that, while the essential process of asymmetric cell division has been conserved over the course of nematode evolution, the underlying spindle movement parameters can combine in various ways. Like other developmental processes, asymmetric cell division is subject to system drift.

Methods and Tools in Genome-wide Association Studies

Many traits, such as height, the response to a given drug, or the susceptibility to certain diseases are presumably co-determined by genetics. Especially in the field of medicine, it is of major interest to identify genetic aberrations that alter an individual's risk to develop a certain phenotypic trait. Addressing this question requires the availability of comprehensive, high-quality genetic datasets. The technological advancements and the decreasing cost of genotyping in the last decade led to an increase in such datasets. Parallel to and in line with this technological progress, an analysis framework under the name of genome-wide association studies was developed to properly collect and analyze these data. Genome-wide association studies aim at finding statistical dependencies-or associations-between a trait of interest and point-mutations in the DNA. The statistical models used to detect such associations are diverse, spanning the whole range from the frequentist to the Bayesian setting.Since genetic datasets are inherently high-dimensional, the search for associations poses not only a statistical but also a computational challenge. As a result, a variety of toolboxes and software packages have been developed, each implementing different statistical methods while using various optimizations and mathematical techniques to enhance the computations.This chapter is devoted to the discussion of widely used methods and tools in genome-wide association studies. We present the different statistical models and the assumptions on which they are based, explain peculiarities of the data that have to be accounted for and, most importantly, introduce commonly used tools and software packages for the different tasks in a genome-wide association study, complemented with examples for their application.

Ephemeral-habitat colonization and neotropical species richness of Caenorhabditis nematodes

Background

The drivers of species co-existence in local communities are especially enigmatic for assemblages of morphologically cryptic species. Here we characterize the colonization dynamics and abundance of nine species of Caenorhabditis nematodes in neotropical French Guiana, the most speciose known assemblage of this genus, with resource use overlap and notoriously similar external morphology despite deep genomic divergence.

Methods

To characterize the dynamics and specificity of colonization and exploitation of ephemeral resource patches, we conducted manipulative field experiments and the largest sampling effort to date for Caenorhabditis outside of Europe. This effort provides the first in-depth quantitative analysis of substrate specificity for Caenorhabditis in natural, unperturbed habitats.

Results

We amassed a total of 626 strain isolates from nine species of Caenorhabditis among 2865 substrate samples. With the two new species described here (C. astrocarya and C. dolens), we estimate that our sampling procedures will discover few additional species of these microbivorous animals in this tropical rainforest system. We demonstrate experimentally that the two most prevalent species (C. nouraguensis and C. tropicalis) rapidly colonize fresh resource patches, whereas at least one rarer species shows specialist micro-habitat fidelity.

Conclusion

Despite the potential to colonize rapidly, these ephemeral patchy resources of rotting fruits and flowers are likely to often remain uncolonized by Caenorhabditis prior to their complete decay, implying dispersal-limited resource exploitation. We hypothesize that a combination of rapid colonization, high ephemerality of resource patches, and species heterogeneity in degree of specialization on micro-habitats and life histories enables a dynamic co-existence of so many morphologically cryptic species of Caenorhabditis.

Electronic supplementary material

The online version of this article (10.1186/s12898-017-0150-z) contains supplementary material, which is available to authorized users.

The extraordinary AFD thermosensor of C. elegans

The nematode C. elegans exhibits complex thermal experience-dependent navigation behaviors in response to environmental temperature changes of as little as 0.01°C over a > 10°C temperature range. The remarkable thermosensory abilities of this animal are mediated primarily via the single pair of AFD sensory neurons in its head. In this review, we describe the contributions of AFD to thermosensory behaviors and temperature-dependent regulation of organismal physiology. We also discuss the mechanisms that enable this neuron type to adapt to recent temperature experience and to exhibit extraordinary thermosensitivity over a wide dynamic range.

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

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

Inverse correlation between longevity and developmental rate among wild C. elegans strains

Genetic studies using model organisms have shown that many long-lived mutants display impaired fitness, such as reduced fecundity and delayed development. However, in several wild animals, the association between longevity and fitness does not seem to be inevitable. Thus, the relationship between longevity and fitness in wild organisms remains inconclusive. Here, we determined the correlation between lifespan and fitness, developmental rate and brood size, by using 16 wild-derived C. elegans strains originated from various geographic areas. We found a negative correlation between lifespan and developmental rate. In contrast, we did not find such negative correlation between longevity and developmental rate among the individuals of C. elegans strains. These data imply that polymorphic genetic variants among wild isolates determine resource allocation to longevity and developmental rate.

Natural Variation in the Distribution and Abundance of Transposable Elements Across the Caenorhabditis elegans Species

Transposons are mobile DNA elements that generate both adaptive and deleterious phenotypic variation thereby driving genome evolution. For these reasons, genomes have mechanisms to regulate transposable element (TE) activity. Approximately 12–16% of the Caenorhabditis elegans genome is composed of TEs, of which the majority are likely inactive. However, most studies of TE activity have been conducted in the laboratory strain N2, which limits our knowledge of the effects of these mobile elements across natural populations. We analyzed the distribution and abundance of TEs in 208 wild C. elegans strains to better understand how transposons contribute to variation in natural populations. We identified 3,397 TEs as compared with the reference strain, of which 2,771 are novel insertions and 241 are TEs that have been excised in at least one wild strain. Likely because of their hypothesized deleterious effects, we find that TEs are found at low allele frequencies throughout the population, and we predict functional effects of TE insertions. The abundances of TEs reflect their activities, and these data allowed us to perform both genome-wide association mappings and rare variant correlations to reveal several candidate genes that impact TE regulation, including small regulatory piwi-interacting RNAs and chromatin factors. Because TE variation in natural populations could underlie phenotypic variation for organismal and behavioral traits, the transposons that we identified and their regulatory mechanisms can be used in future studies to explore the genomics of complex traits and evolutionary changes.

The mutational decay of male-male and hermaphrodite-hermaphrodite competitive fitness in the androdioecious nematode C. elegans

Androdioecious Caenorhabditis have a high frequency of self-compatible hermaphrodites and a low frequency of males. The effects of mutations on male fitness are of interest for two reasons. First, when males are rare, selection on male-specific mutations is less efficient than in hermaphrodites. Second, males may present a larger mutational target than hermaphrodites because of the different ways in which fitness accrues in the two sexes. We report the first estimates of male-specific mutational effects in an androdioecious organism. The rate of male-specific inviable or sterile mutations is ⩽5 × 10^-4/generation, below the rate at which males would be lost solely due to those kinds of mutations. The rate of mutational decay of male competitive fitness is ~ 0.17%/generation; that of hermaphrodite competitive fitness is ~ 0.11%/generation. The point estimate of ~ 1.5X faster rate of mutational decay of male fitness is nearly identical to the same ratio in Drosophila. Estimates of mutational variance (V_M) for male mating success and competitive fitness are not significantly different from zero, whereas V_M for hermaphrodite competitive fitness is similar to that of non-competitive fitness. Two independent estimates of the average selection coefficient against mutations affecting hermaphrodite competitive fitness agree to within two-fold, 0.33-0.5%.

Genetic variation in glia-neuron signalling modulates ageing rate

The rate of behavioural decline in the ageing population is remarkably variable among individuals. Despite the considerable interest in studying natural variation in ageing rate to identify factors that control healthy ageing, no such factor has yet been found. Here we report a genetic basis for variation in ageing rates in Caenorhabditis elegans. We find that C. elegans isolates show diverse lifespan and age-related declines in virility, pharyngeal pumping, and locomotion. DNA polymorphisms in a novel peptide-coding gene, named regulatory-gene-for-behavioural-ageing-1 (rgba-1), and the neuropeptide receptor gene npr-28 influence the rate of age-related decline of worm mating behaviour; these two genes might have been subjected to recent selective sweeps. Glia-derived RGBA-1 activates NPR-28 signalling, which acts in serotonergic and dopaminergic neurons to accelerate behavioural deterioration. This signalling involves the SIR-2.1-dependent activation of the mitochondrial unfolded protein response, a pathway that modulates ageing. Thus, natural variation in neuropeptide-mediated glia-neuron signalling modulates the rate of ageing in C. elegans.

Hierarchical compression of Caenorhabditis elegans locomotion reveals phenotypic differences in the organization of behaviour

Regularities in animal behaviour offer insights into the underlying organizational and functional principles of nervous systems and automated tracking provides the opportunity to extract features of behaviour directly from large-scale video data. Yet how to effectively analyse such behavioural data remains an open question. Here, we explore whether a minimum description length principle can be exploited to identify meaningful behaviours and phenotypes. We apply a dictionary compression algorithm to behavioural sequences from the nematode worm Caenorhabditis elegans freely crawling on an agar plate both with and without food and during chemotaxis. We find that the motifs identified by the compression algorithm are rare but relevant for comparisons between worms in different environments, suggesting that hierarchical compression can be a useful step in behaviour analysis. We also use compressibility as a new quantitative phenotype and find that the behaviour of wild-isolated strains of C. elegans is more compressible than that of the laboratory strain N2 as well as the majority of mutant strains examined. Importantly, in distinction to more conventional phenotypes such as overall motor activity or aggregation behaviour, the increased compressibility of wild isolates is not explained by the loss of function of the gene npr-1, which suggests that erratic locomotion is a laboratory-derived trait with a novel genetic basis. Because hierarchical compression can be applied to any sequence, we anticipate that compressibility can offer insights into the organization of behaviour in other animals including humans.

Non-Mendelian assortment of homologous autosomes of different sizes in males is the ancestral state in the Caenorhabditis lineage

Organismal genome sizes vary by six orders of magnitude and appear positively correlated with organismal size and complexity. Neutral models have been proposed to explain the broad patterns of genome size variation based on organism population sizes. In the Caenorhabditis genus, hermaphrodite genomes are smaller than those of gonochoristic species. One possible driving force for this genome size difference could be non-random chromosome segregation. In Caenorhabditis elegans, chromosome assortment is non-independent and violates Mendel's second law. In males, the shorter homologue of a heterozygous autosome pair preferentially co-segregates with the X chromosome while the longer one preferentially co-segregates with the nullo-X (O) chromosome in a process we call "skew". Since hermaphrodites preferentially receive the shorter chromosomes and can start populations independently, their genome size would be predicted to decrease over evolutionary time. If skew is an important driver for genome size reduction in hermaphroditic Caenorhabditis species, then it should be present in all congeneric species. In this study, we tested this hypothesis and found that skew is present in all eight examined species. Our results suggest that skew is likely the ancestral state in this genus. More speculatively, skew may drive genome size patterns in hermaphroditic species in other nematodes.

Evolution: A Parthenogenetic Nematode Shows How Animals Become Sexless

Most animals have male and female sexes, implying that sex is ancient and beneficial; yet some have survived for millions of years without sex. The genome of the parthenogenetic nematode Diploscapter pachys gives clues as to how 'ancient asexual' animals can exist.

Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.

Why do Individuals Differ in Viral Susceptibility? A Story Told by Model Organisms

Viral susceptibility and disease progression is determined by host genetic variation that underlies individual differences. Genetic polymorphisms that affect the phenotype upon infection have been well-studied for only a few viruses, such as HIV-1 and Hepatitis C virus. However, even for well-studied viruses the genetic basis of individual susceptibility differences remains elusive. Investigating the effect of causal polymorphisms in humans is complicated, because genetic methods to detect rare or small-effect polymorphisms are limited and genetic manipulation is not possible in human populations. Model organisms have proven a powerful experimental platform to identify and characterize polymorphisms that underlie natural variations in viral susceptibility using quantitative genetic tools. We summarize and compare the genetic tools available in three main model organisms, Mus musculus, Drosophila melanogaster, and Caenorhabditis elegans, and illustrate how these tools can be applied to detect polymorphisms that determine the viral susceptibility. Finally, we analyse how candidate polymorphisms from model organisms can be used to shed light on the underlying mechanism of individual variation. Insights in causal polymorphisms and mechanisms underlying individual differences in viral susceptibility in model organisms likely provide a better understanding in humans.

The genetic basis of natural variation in a phoretic behavior

Phoresy is a widespread form of commensalism that facilitates dispersal of one species through an association with a more mobile second species. Dauer larvae of the nematode Caenorhabditis elegans exhibit a phoretic behavior called nictation, which could enable interactions with animals such as isopods or snails. Here, we show that natural C. elegans isolates differ in nictation. We use quantitative behavioral assays and linkage mapping to identify a genetic locus (nict-1) that mediates the phoretic interaction with terrestrial isopods. The nict-1 locus contains a Piwi-interacting small RNA (piRNA) cluster; we observe that the Piwi Argonaute PRG-1 is involved in the regulation of nictation. Additionally, this locus underlies a trade-off between offspring production and dispersal. Variation in the nict-1 locus contributes directly to differences in association between nematodes and terrestrial isopods in a laboratory assay. In summary, the piRNA-rich nict-1 locus could define a novel mechanism underlying phoretic interactions.

Nematodes use a characteristic set of movements, called nictation, to hitchhike on more mobile animals. Here, Lee et al. identify a genetic locus in the nematode Caenorhabditis elegans that underlies nictation and contributes to successful hitchhiking, but at expense of reduced offspring production.

Natural variation in a single amino acid substitution underlies physiological responses to topoisomerase II poisons

Many chemotherapeutic drugs are differentially effective from one patient to the next. Understanding the causes of this variability is a critical step towards the development of personalized treatments and improvements to existing medications. Here, we investigate sensitivity to a group of anti-neoplastic drugs that target topoisomerase II using the model organism Caenorhabditis elegans. We show that wild strains of C. elegans vary in their sensitivity to these drugs, and we use an unbiased genetic approach to demonstrate that this natural variation is explained by a methionine-to-glutamine substitution in topoisomerase II (TOP-2). The presence of a non-polar methionine at this residue increases hydrophobic interactions between TOP-2 and its poison etoposide, as compared to a polar glutamine. We hypothesize that this stabilizing interaction results in increased genomic instability in strains that contain a methionine residue. The residue affected by this substitution is conserved from yeast to humans and is one of the few differences between the two human topoisomerase II isoforms (methionine in hTOPIIα and glutamine in hTOPIIβ). We go on to show that this amino acid difference between the two human topoisomerase isoforms influences cytotoxicity of topoisomerase II poisons in human cell lines. These results explain why hTOPIIα and hTOPIIβ are differentially affected by various poisons and demonstrate the utility of C. elegans in understanding the genetics of drug responses.

The severe cytotoxic effects associated with anti-neoplastic treatment regimens make it difficult to assess the contributions of genetic variation on treatment responses in clinical settings. Therefore, we leveraged genetic diversity present in the metazoan model nematode Caenorhabditis elegans to identify genetic variants that contribute to differential susceptibility to a broadly administered class of anti-neoplastic compounds that poison the activity of topoisomerase II enzymes. We show that wild C. elegans isolates contain either glutamine or methionine at a highly conserved residue of the topoisomerase II (TOP-2) protein and that this substitution is predictive of animal responses to the topoisomerase II poisons etoposide, teniposide, dactinomycin, and XK469. Interestingly, the two human versions of this protein, hTOPIIα and hTOPIIβ, contain a methionine or glutamine at the corresponding residue, respectively. We show that this difference between the two human topoisomerase II isoforms contributes to the differential cytotoxicity induced by these drugs. Taken together, our results highlight the power of studying the effects of natural genetic variation on drug responses in a model organism and propose methods to develop new drugs that have increased affinity for the desired hTOPIIα isoform expressed in tumor cells.

Experimental Evolution with Caenorhabditis Nematodes

The hermaphroditic nematode Caenorhabditis elegans has been one of the primary model systems in biology since the 1970s, but only within the last two decades has this nematode also become a useful model for experimental evolution. Here, we outline the goals and major foci of experimental evolution with C. elegans and related species, such as C. briggsae and C. remanei, by discussing the principles of experimental design, and highlighting the strengths and limitations of Caenorhabditis as model systems. We then review three exemplars of Caenorhabditis experimental evolution studies, underlining representative evolution experiments that have addressed the: (1) maintenance of genetic variation; (2) role of natural selection during transitions from outcrossing to selfing, as well as the maintenance of mixed breeding modes during evolution; and (3) evolution of phenotypic plasticity and its role in adaptation to variable environments, including host-pathogen coevolution. We conclude by suggesting some future directions for which experimental evolution with Caenorhabditis would be particularly informative.

A maternal-effect selfish genetic element in Caenorhabditis elegans

Selfish genetic elements spread in natural populations and have an important role in genome evolution. We discovered a selfish element causing embryonic lethality in crosses between wild strains of the nematode Caenorhabditis elegans The element is made up of sup-35, a maternal-effect toxin that kills developing embryos, and pha-1, its zygotically expressed antidote. pha-1 has long been considered essential for pharynx development on the basis of its mutant phenotype, but this phenotype arises from a loss of suppression of sup-35 toxicity. Inactive copies of the sup-35/pha-1 element show high sequence divergence from active copies, and phylogenetic reconstruction suggests that they represent ancestral stages in the evolution of the element. Our results suggest that other essential genes identified by genetic screens may turn out to be components of selfish elements.

The Evolutionary Interplay between Adaptation and Self-Fertilization

Genome-wide surveys of nucleotide polymorphisms, obtained from next-generation sequencing, have uncovered numerous examples of adaptation in self-fertilizing organisms, especially regarding changes to climate, geography, and reproductive systems. Yet existing models for inferring attributes of adaptive mutations often assume idealized outcrossing populations, which risks mischaracterizing properties of these variants. Recent theoretical work is emphasizing how various aspects of self-fertilization affects adaptation, yet empirical data on these properties are lacking. We review theoretical and empirical studies demonstrating how self-fertilization alters the process of adaptation, illustrated using examples from current sequencing projects. We propose ideas for how future research can more accurately quantify aspects of adaptation in self-fertilizers, including incorporating the effects of standing variation, demographic history, and polygenic adaptation.

Analysis of large-scale next-generation sequencing datasets are finding more examples of adaptive evolution at the genomic level.

Advances in theoretical work has demonstrated how self-fertilisation affects different aspects of adaptation in these organisms, compared to outcrossers.

Current software and statistical methods do not take different mating systems into account, which risks mischaracterising the presence or strength of adaptive mutations from genome scans.

Development of new mathematical and statistical methods that explicitly consider self-fertilization and associated demographic effects will enable researchers to more accurately quantify adaptation in these organisms.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

The Natural Biotic Environment of Caenorhabditis elegans

Organisms evolve in response to their natural environment. Consideration of natural ecological parameters are thus of key importance for our understanding of an organism's biology. Curiously, the natural ecology of the model species Caenorhabditis elegans has long been neglected, even though this nematode has become one of the most intensively studied models in biological research. This lack of interest changed ∼10 yr ago. Since then, an increasing number of studies have focused on the nematode's natural ecology. Yet many unknowns still remain. Here, we provide an overview of the currently available information on the natural environment of C. elegans We focus on the biotic environment, which is usually less predictable and thus can create high selective constraints that are likely to have had a strong impact on C. elegans evolution. This nematode is particularly abundant in microbe-rich environments, especially rotting plant matter such as decomposing fruits and stems. In this environment, it is part of a complex interaction network, which is particularly shaped by a species-rich microbial community. These microbes can be food, part of a beneficial gut microbiome, parasites and pathogens, and possibly competitors. C. elegans is additionally confronted with predators; it interacts with vector organisms that facilitate dispersal to new habitats, and also with competitors for similar food environments, including competitors from congeneric and also the same species. Full appreciation of this nematode's biology warrants further exploration of its natural environment and subsequent integration of this information into the well-established laboratory-based research approaches.

Genetic relationships among Eriobotrya species revealed by genome-wide RAD sequence data

Restriction site‐associated DNA sequencing (RAD‐seq) was used to illuminate the genetic relationships among Eriobotrya species. The raw data were filtered, and 221 million clean reads were used for further analysis. A total of 1,983,332 SNPs were obtained from 23 Eriobotrya species and two relative genera. We obtained similar results by neighbor‐joining and maximum likelihood phylogenetic trees. All Eriobotrya plants grouped together into a big clade, and two out‐groups clustered together into a single or separate clade. Chinese and Vietnam accessions were distributed throughout the dendrogram. There was nonsignificant correlation between genotype and geographical distance. However, clustering results were correlated with leaf size to some extent. The Eriobotrya species could be divided into following three groups based on leaf size and phylogenetic analysis: group A and group B comprised of small leaves with <10 cm length except E. stipularis (16.76 cm), and group C can be further divided into two subgroups, which contained medium‐size leaves with a leaf length ranged from 10 to 20 cm and a leaf length bigger than 20 cm.

Fitness Effects of Thermal Stress Differ Between Outcrossing and Selfing Populations in Caenorhabditis elegans

The maintenance of males and outcrossing is widespread, despite considerable costs of males. By enabling recombination between distinct genotypes, outcrossing may be advantageous during adaptation to novel environments and if so, it should be selected for under environmental challenge. However, a given environmental change may influence fitness of male, female, and hermaphrodite or asexual individuals differently, and hence the relationship between reproductive system and dynamics of adaptation to novel conditions may not be driven solely by the level of outcrossing and recombination. This has important implications for studies investigating the evolution of reproductive modes in the context of environmental changes, and for the extent to which their findings can be generalized. Here, we use Caenorhabditis elegans-a free-living nematode species in which hermaphrodites (capable of selfing but not cross-fertilizing each other) coexist with males (capable of fertilizing hermaphrodites)-to investigate the response of wild type as well as obligatorily outcrossing and obligatorily selfing lines to stressfully increased ambient temperature. We found that thermal stress affects fitness of outcrossers much more drastically than that of selfers. This shows that apart from the potential for recombination, the selective pressures imposed by the same environmental change can differ between populations expressing different reproductive systems and affect their adaptive potential.

Variation in rates of spontaneous male production within the nematode species Pristionchus pacificus supports an adaptive role for males and outcrossing

Background

The nematode species Pristionchus pacificus has an androdioecious mating system in which populations consist of self-fertilizing hermaphrodites and relatively few males. The prevalence of males in such a system is likely to depend on the relative pros and cons of outcrossing. While outcrossing generates novel allelic combinations and can therefore increase adaptive potential, it may also disrupt the potentially beneficial consequences of repeated generations of selfing. These include purging of deleterious alleles, inheritance of co-adapted allele complexes, improved hermaphrodite fitness and increased population growth. Here we use experimental and population genetic approaches to test hypotheses relating to male production and outcrossing in laboratory and natural populations of P. pacificus sampled from the volcanic island of La Réunion.

Results

We find a significant interaction between sampling locality and temperature treatment influencing rates of spontaneous male production in the laboratory. While strains isolated at higher altitude, cooler localities produce a higher proportion of male offspring at 25 °C relative to 20 or 15 °C, the reverse pattern is seen in strains isolated from warmer, low altitude localities. Linkage disequilibrium extends across long physical distances, but fails to approach levels reported for the partially selfing nematode species Caenorhabditis elegans. Finally, we find evidence for admixture between divergent genetic lineages.

Conclusions

Elevated rates of laboratory male generation appear to occur under environmental conditions which differ from those experienced by populations in nature. Such elevated male generation may result in higher outcrossing rates, hence driving increased effective recombination and the creation of potentially adaptive novel allelic combinations. Patterns of linkage disequilibrium decay support selfing as the predominant reproductive strategy in P. pacificus. Finally, despite the potential for outcrossing depression, our results suggest admixture has occurred between distinct genetic lineages since their independent colonization of the island, suggesting outcrossing depression may not be uniform in this species.

Electronic supplementary material

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