Widespread genomic incompatibilities in Caenorhabditis elegans

A role for worm cutl-24 in background- and parent-of-origin-dependent ER stress resistance

BACKGROUND

Organisms in the wild can acquire disease- and stress-resistance traits that outstrip the programs endogenous to humans. Finding the molecular basis of such natural resistance characters is a key goal of evolutionary genetics. Standard statistical-genetic methods toward this end can perform poorly in organismal systems that lack high rates of meiotic recombination, like Caenorhabditis worms.

RESULTS

Here we discovered unique ER stress resistance in a wild Kenyan C. elegans isolate, which in inter-strain crosses was passed by hermaphrodite mothers to hybrid offspring. We developed an unbiased version of the reciprocal hemizygosity test, RH-seq, to explore the genetics of this parent-of-origin-dependent phenotype. Among top-scoring gene candidates from a partial-coverage RH-seq screen, we focused on the neuronally-expressed, cuticlin-like gene cutl-24 for validation. In gene-disruption and controlled crossing experiments, we found that cutl-24 was required in Kenyan hermaphrodite mothers for ER stress tolerance in their inter-strain hybrid offspring; cutl-24 was also a contributor to the trait in purebred backgrounds.

CONCLUSIONS

These data establish the Kenyan strain allele of cutl-24 as a determinant of a natural stress-resistant state, and they set a precedent for the dissection of natural trait diversity in invertebrate animals without the need for a panel of meiotic recombinants.

Mating system and speciation I: Accumulation of genetic incompatibilities in allopatry

Self-fertilisation is widespread among hermaphroditic species across the tree of life. Selfing has many consequences on the genetic diversity and the evolutionary dynamics of populations, which may in turn affect macroevolutionary processes such as speciation. On the one hand, because selfing increases genetic drift and reduces migration rate among populations, it may be expected to promote speciation. On the other hand, because selfing reduces the efficacy of selection, it may be expected to hamper ecological speciation. To better understand under which conditions and in which direction selfing affects the build-up of reproductive isolation, an explicit population genetics model is required. Here, we focus on the interplay between genetic drift, selection and genetic linkage by studying speciation without gene flow. We test how fast populations with different rates of selfing accumulate mutations leading to genetic incompatibilities. When speciation requires populations to pass through a fitness valley caused by underdominant and compensatory mutations, selfing reduces the depth and/or breadth of the valley, and thus overall facilitates the fixation of incompatibilities. When speciation does not require populations to pass through a fitness valley, as for Bateson-Dobzhanzky-Muller incompatibilities (BDMi), the lower effective population size and higher genetic linkage in selfing populations both facilitate the fixation of incompatibilities. Interestingly, and contrary to intuitive expectations, local adaptation does not always accelerate the fixation of incompatibilities in outcrossing relative to selfing populations. Our work helps to clarify how incompatibilities accumulate in selfing vs. outcrossing lineages, and has repercussions on the pace of speciation as well as on the genetic architecture of reproductive isolation.

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.

Hybrid fitness effects modify fixation probabilities of introgressed alleles

Hybridization is a common occurrence in natural populations, and introgression is a major source of genetic variation. Despite the evolutionary importance of adaptive introgression, classical population genetics theory does not take into account hybrid fitness effects. Specifically, heterosis (i.e. hybrid vigor) and Dobzhansky–Muller incompatibilities influence the fates of introgressed alleles. Here, we explicitly account for polygenic, unlinked hybrid fitness effects when tracking a rare introgressed marker allele. These hybrid fitness effects quickly decay over time due to repeated backcrossing, enabling a separation-of-timescales approach. Using diffusion and branching process theory in combination with computer simulations, we formalize the intuition behind how hybrid fitness effects affect introgressed alleles. We find that hybrid fitness effects can significantly hinder or boost the fixation probability of introgressed alleles, depending on the relative strength of heterosis and Dobzhansky–Muller incompatibilities effects. We show that the inclusion of a correction factor (α, representing the compounded effects of hybrid fitness effects over time) into classic population genetics theory yields accurate fixation probabilities. Despite having a strong impact on the probability of fixation, hybrid fitness effects only subtly change the distribution of fitness effects of introgressed alleles that reach fixation. Although strong Dobzhansky–Muller incompatibility effects may expedite the loss of introgressed alleles, fixation times are largely unchanged by hybrid fitness effects.

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

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

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.

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

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

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

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

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 Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

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

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

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

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.

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.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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.

Patterns of Transmission Ratio Distortion in Interspecific Lettuce Hybrids Reveal a Sex-Independent Gametophytic Barrier

Interspecific crosses can result in progeny with reduced vitality or fertility due to genetic incompatibilities between species, a phenomenon known as hybrid incompatibility (HI). HI is often caused by a bias against deleterious allele combinations, which results in transmission ratio distortion (TRD). Here, we determined the genome-wide distribution of HI between wild lettuce, Lactuca saligna, and cultivated lettuce, L. sativa, in a set of backcross inbred lines (BILs) with single introgression segments from L. saligna introgressed into a L. sativa genetic background. Almost all BILs contained an introgression segment in a homozygous state except a few BILs, for which we were able to obtain only a single heterozygous introgression. Their inbred progenies displayed severe TRD with a bias toward the L. sativa allele and complete nontransmission of the homozygous L. saligna introgression, i.e., absolute HI. These HI might be caused by deleterious heterospecific allele combinations at two loci. We used an multilocus segregating interspecific F2 population to identify candidate conspecific loci that can nullify the HI in BILs. Segregation analysis of developed double-introgression progenies showed nullification of three HI and proved that these HI are explained by nuclear pairwise incompatibilities. One of these digenic HI showed 29% reduced seed set and its pattern of TRD pointed to a sex-independent gametophytic barrier. Namely, this HI was caused by complete nontransmission of one heterospecific allele combination at the haploid stage, surprisingly in both male and female gametophytes. Our study shows that two-locus incompatibility systems contribute to reproductive barriers among Lactuca species.

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.

Dramatic evolution of body length due to postembryonic changes in cell size in a newly discovered close relative of Caenorhabditis elegans

Understanding morphological diversity-and morphological constraint-has been a central question in evolutionary biology since its inception. Nematodes of the genus Caenorhabditis, which contains the well-studied model organism C. elegans, display remarkable morphological consistency in the face of extensive genetic divergence. Here, we provide a description of the broad developmental patterns of a newly discovered species, C. sp. 34, which was isolated from fresh figs in Okinawa and which is among the closest known relatives of C. elegans. C. sp. 34 displays an extremely large body size; it can grow to be nearly twice as long as C. elegans and all other known members of the genus. Observations of the timing of developmental milestones reveal that C. sp. 34 develops about twice as slowly as C. elegans. Measurements of embryonic and larval size show that the size difference between C. sp. 34 and C. elegans is largely due to postembryonic events, particularly during the transition from larval to adult stages. This difference in size is not attributable to differences in germ line chromosome number or the number of somatic cells. The overall difference in body size is therefore largely attributable to changes in cell size via increased cytoplasmic volume. Because of its close relationship to C. elegans, the distinctness of C. sp. 34 provides an ideal system for the detailed analysis of evolutionary diversification. The context of over 40 years of C. elegans developmental genetics also reveals clues into how natural selection and developmental constraint act jointly to promote patterns of morphological stasis and divergence in this group.

Use of a Sibling Subtraction Method for Identifying Causal Mutations in Caenorhabditis elegans by Whole-Genome Sequencing

Whole-genome sequencing (WGS) is an indispensable tool for identifying causal mutations obtained from genetic screens. To reduce the number of causal mutation candidates typically uncovered by WGS, Caenorhabditis elegans researchers have developed several strategies. One involves crossing N2-background mutants to the polymorphic Hawaiian (HA) strain, which can be used to simultaneously identify mutant strain variants and obtain high-density mapping information. This approach, however, is not well suited for uncovering mutations in complex genetic backgrounds, and HA polymorphisms can alter phenotypes. Other approaches make use of DNA variants present in the initial background or introduced by mutagenesis. This information is used to implicate genomic regions with high densities of DNA lesions that persist after backcrossing, but these methods can provide lower resolution than HA mapping. To identify suppressor mutations using WGS, we developed an approach termed the sibling subtraction method (SSM). This method works by eliminating variants present in both mutants and their nonmutant siblings, thus greatly reducing the number of candidates. We used this method with two members of the C. elegans NimA-related kinase family, nekl-2 and nekl-3. Combining weak aphenotypic alleles of nekl-2 and nekl-3 leads to penetrant molting defects and larval arrest. We isolated ∼50 suppressors of nekl-2; nekl-3 synthetic lethality using F1 clonal screening methods and a peel-1–based counterselection strategy. When applied to five of the suppressors, SSM led to only one to four suppressor candidates per strain. Thus SSM is a powerful approach for identifying causal mutations in any genetic background and provides an alternative to current methods.

X exceptionalism in Caenorhabditis speciation

Speciation genetics research in diverse organisms shows the X-chromosome to be exceptional in how it contributes to "rules" of speciation. Until recently, however, the nematode phylum has been nearly silent on this issue, despite the model organism Caenorhabditis elegans having touched most other topics in biology. Studies of speciation with Caenorhabditis accelerated with the recent discovery of species pairs showing partial interfertility. The resulting genetic analyses of reproductive isolation in nematodes demonstrate key roles for the X-chromosome in hybrid male sterility and inviability, opening up new understanding of the genetic causes of Haldane's rule, Darwin's corollary to Haldane's rule, and enabling tests of the large-X effect hypothesis. Studies to date implicate improper chromatin regulation of the X-chromosome by small RNA pathways as integral to hybrid male dysfunction. Sexual transitions in reproductive mode to self-fertilizing hermaphroditism inject distinctive molecular evolutionary features into the speciation process for some species. Caenorhabditis also provides unique opportunities for analysis in a system with XO sex determination that lacks a Y-chromosome, sex chromosome-dependent sperm competition differences and mechanisms of gametic isolation, exceptional accessibility to the development process and rapid experimental evolution. As genetic analysis of reproductive isolation matures with investigation of multiple pairs of Caenorhabditis species and new species discovery, nematodes will provide a powerful complement to more established study organisms for deciphering the genetic basis of and rules to speciation.

Variation in reproductive isolation across a species range

Reproductive isolation is often variable within species, a phenomenon that while largely ignored by speciation studies, can be leveraged to gain insight into the potential mechanisms driving the evolution of genetic incompatibilities. We used experimental greenhouse crosses to characterize patterns of reproductive isolation among three divergent genetic lineages of Campanulastrum americanum that occur in close geographic proximity in the Appalachian Mountains. Substantial, asymmetrical reproductive isolation for survival due to cytonuclear incompatibility was found among the lineages (up to 94% reduction). Moderate reductions in pollen viability, as well as cytoplasmic male sterility, were also found between some Mountain populations. We then compared these results to previously established patterns of reproductive isolation between these Mountain lineages and a fourth, widespread Western lineage to fully characterize reproductive isolation across the complete geographic and genetic range of C. americanum. Reproductive isolation for survival and pollen viability was consistent across studies, indicating the evolution of the underlying genetic incompatibilities is primarily determined by intrinsic factors. In contrast, reproductive isolation for germination was only found when crossing Mountain populations with the Western lineage, suggesting the underlying genetic incompatibility is likely influenced by environmental or demographic differences between the two lineages. Cytoplasmic male sterility was also limited in occurrence, being restricted to a handful of Mountain populations in a narrow geographic range. These findings illustrate the complexity of speciation by demonstrating multiple, independent genetic incompatibilities that lead to a mosaic of genetic divergence and reproductive isolation across a species range.

Genetic Dissection of Morphometric Traits Reveals That Phytochrome B Affects Nucleus Size and Heterochromatin Organization in Arabidopsis thaliana

Microscopically visible chromatin is partitioned into two major components in Arabidopsis thaliana nuclei. On one hand, chromocenters are conspicuous foci of highly condensed "heterochromatic" domains that contain mostly repeated sequences. On the other hand, less condensed and gene-rich "euchromatin" emanates from these chromocenters. This differentiation, together with the dynamic nature of chromatin compaction in response to developmental and environmental stimuli, makes Arabidopsis a powerful system for studying chromatin organization and dynamics. Heterochromatin dynamics can be monitored by measuring the Heterochromatin Index, i.e., the proportion of nuclei displaying well-defined chromocenters, or the DNA fraction of chromocenters (relative heterochromatin fraction). Both measures are composite traits, thus their values represent the sum of effects of various underlying morphometric properties. We exploited genetic variation between natural occurring accessions to determine the genetic basis of individual nucleus and chromocenter morphometric parameters (area, perimeter, density, roundness, and heterogeneity) that together determine chromatin compaction. Our novel reductionist genetic approach revealed quantitative trait loci (QTL) for all measured traits. Genomic colocalization among QTL was limited, which suggests a complex genetic regulation of chromatin compaction. Yet genomic intervals of QTL for nucleus size (area and perimeter) both overlap with a known QTL for heterochromatin compaction that is explained by natural polymorphism in the red/far-red light and temperature receptor Phytochrome B. Mutant analyses and genetic complementation assays show that Phytochrome B is a negative regulator of nucleus size, revealing that perception of climatic conditions by a Phytochrome-mediated hub is a major determinant for coordinating nucleus size and heterochromatin compaction.

Species-wide survey reveals the various flavors of intraspecific reproductive isolation in yeast

Exploring the origin and extent of reproductive isolation within the same species is valuable to capture early events to the onset of speciation. In multiple genetic models, reproductive isolation was recently observed at the intraspecific scale, indicating that the raw potential for speciation segregates readily within populations, which could be a rule rather than an exception in a broad context. We briefly recapitulate the molecular evidence of intrinsic post-zygotic isolation in major model organisms including Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster and their close relatives. We then focus on recent advances in yeast and review the genetic basis of post-zygotic isolation within and between multiple members of the Saccharomyces genus, especially in Saccharomyces cerevisiae We discuss the role of various mechanisms involved in the onset of reproductive isolation including DNA sequence divergence, chromosomal rearrangement, cytonuclear as well as nuclear-nuclear genetic incompatibilities and provide a comparative view along a continuum of genetic differentiation, which encompasses intraspecific populations, recent delineating nascent species as well as closely related sister species in the same subphylum.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Natural Genetic Variation Influences Protein Abundances in C. elegans Developmental Signalling Pathways

Complex traits, including common disease-related traits, are affected by many different genes that function in multiple pathways and networks. The apoptosis, MAPK, Notch, and Wnt signalling pathways play important roles in development and disease progression. At the moment we have a poor understanding of how allelic variation affects gene expression in these pathways at the level of translation. Here we report the effect of natural genetic variation on transcript and protein abundance involved in developmental signalling pathways in Caenorhabditis elegans. We used selected reaction monitoring to analyse proteins from the abovementioned four pathways in a set of recombinant inbred lines (RILs) generated from the wild-type strains N2 (Bristol) and CB4856 (Hawaii) to enable quantitative trait locus (QTL) mapping. About half of the cases from the 44 genes tested showed a statistically significant change in protein abundance between various strains, most of these were however very weak (below 1.3-fold change). We detected a distant QTL on the left arm of chromosome II that affected protein abundance of the phosphatidylserine receptor protein PSR-1, and two separate QTLs that influenced embryonic and ionizing radiation-induced apoptosis on chromosome IV. Our results demonstrate that natural variation in C. elegans is sufficient to cause significant changes in signalling pathways both at the gene expression (transcript and protein abundance) and phenotypic levels.

Collective effects of common SNPs in foraging decisions in Caenorhabditis elegans and an integrative method of identification of candidate genes

Optimal foraging decision is a quantitative flexible behavior, which describes the time at which animals choose to abandon a depleting food supply. The total minor allele content (MAC) in an individual has been shown to correlate with quantitative variations in complex traits. We have studied the role of MAC in the decision to leave a food lawn in recombinant inbred advanced intercross lines (RIAILs) of Caenorhabditis elegans. We found a strong link between MAC and the food lawn leaving rates (Spearman r = 0.4, P = 0.005). We identified 28 genes of unknown functions whose expression levels correlated with both MAC and leaving rates. When examined by RNAi experiments, 8 of 10 tested among the 28 affected leaving rates, whereas only 2 of 9 did among genes that were only associated with leaving rates but not MAC (8/10 vs 2/9, P < 0.05). The results establish a link between MAC and the foraging behavior and identify 8 genes that may play a role in linking MAC with the quantitative nature of the trait. The method of correlations with both MAC and traits may find broad applications in high efficiency identification of target genes for other complex traits in model organisms and humans.

Mitochondrial-Nuclear Epistasis Impacts Fitness and Mitochondrial Physiology of Interpopulation Caenorhabditis briggsae Hybrids

In order to identify the earliest genetic changes that precipitate species formation, it is useful to study genetic incompatibilities that cause only mild dysfunction when incompatible alleles are combined in an interpopulation hybrid. Such hybridization within the nematode species Caenorhabditis briggsae has been suggested to result in selection against certain combinations of nuclear and mitochondrial alleles, raising the possibility that mitochondrial-nuclear (mitonuclear) epistasis reduces hybrid fitness. To test this hypothesis, cytoplasmic-nuclear hybrids (cybrids) were created to purposefully disrupt any epistatic interactions. Experimental analysis of the cybrids suggests that mitonuclear discord can result in decreased fecundity, increased lipid content, and increased mitochondrial reactive oxygen species levels. Many of these effects were asymmetric with respect to cross direction, as expected if cytoplasmic-nuclear Dobzhansky-Muller incompatibilities exist. One such effect is consistent with the interpretation that disrupting coevolved mitochondrial and nuclear loci impacts mitochondrial function and organismal fitness. These findings enhance efforts to study the genesis, identity, and maintenance of genetic incompatibilities that precipitate the speciation process.

Genotype-dependent lifespan effects in peptone deprived Caenorhabditis elegans

Dietary restriction appears to act as a general non-genetic mechanism that can robustly prolong lifespan. There have however been reports in many systems of cases where restricted food intake either shortens, or does not affect, lifespan. Here we analyze lifespan and the effect of food restriction via deprived peptone levels on lifespan in wild isolates and introgression lines (ILs) of the nematode Caenorhabditis elegans. These analyses identify genetic variation in lifespan, in the effect of this variation in diet on lifespan and also in the likelihood of maternal, matricidal, hatching. Importantly, in the wild isolates and the ILs, we identify genotypes in which peptone deprivation mediated dietary restriction reduces lifespan. We also identify, in recombinant inbred lines, a locus that affects maternal hatching, a phenotype closely linked to dietary restriction in C. elegans. These results indicate that peptone deprivation mediated dietary restriction affects lifespan in C. elegans in a genotype-dependent manner, reducing lifespan in some genotypes. This may operate by a mechanism similar to dietary restriction.

The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations

I review theoretical models for the evolution of supergenes in the cases of Batesian mimicry in butterflies, distylous plants and sex chromosomes. For each of these systems, I outline the genetic evidence that led to the proposal that they involve multiple genes that interact during ‘complex adaptations’, and at which the mutations involved are not unconditionally advantageous, but show advantages that trade‐off against some disadvantages. I describe recent molecular genetic studies of these systems and questions they raise about the evolution of suppressed recombination. Nonrecombining regions of sex chromosomes have long been known, but it is not yet fully understood why recombination suppression repeatedly evolved in systems in distantly related taxa, but does not always evolve. Recent studies of distylous plants are tending to support the existence of recombination‐suppressed genome regions, which may include modest numbers of genes and resemble recently evolved sex‐linked regions. For Batesian mimicry, however, molecular genetic work in two butterfly species suggests a new supergene scenario, with a single gene mutating to produce initial adaptive phenotypes, perhaps followed by modifiers specifically refining and perfecting the new phenotype.

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.

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

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

Mutations that activate a RAS oncogene are found in a large proportion of human cancers. In this study, we have used the roundworm Caenorhabditis elegans (C. elegans) as a model to investigate how the genetic composition of the animal affects the outcome of oncogenic RAS mutations that activate the MAPK pathway. By comparing the effects of activated RAS/MAPK signaling in two genetically different C. elegans strains, we have identified the monoamine oxidase A (MAOA) gene amx-2 as a negative regulator of RAS/MAPK signaling. MAOA enzymes are primarily known to catalyze the degradation of the neurotransmitters dopamine and serotonin. Here, we show that a specific serotonin degradation product that is produced by MAOA (5-HIAA) inhibits RAS signaling in different organs of C. elegans. Thus, by producing the inhibitory serotonin metabolite 5-HIAA the MAOA enzyme systemically controls the activation of the RAS/MAPK pathway.

The laboratory domestication of Caenorhabditis elegans

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

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

Background

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

Results

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

Conclusions

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