Long-range regulatory polymorphisms affecting a GABA receptor constitute a quantitative trait locus (QTL) for social behavior in Caenorhabditis elegans

Neuronal enhancers fine-tune adaptive circuit plasticity

Neuronal activity-regulated gene expression plays a crucial role in sculpting neural circuits that underpin adaptive brain function. Transcriptional enhancers are now recognized as key components of gene regulation that orchestrate spatiotemporally precise patterns of gene transcription. We propose that the dynamics of enhancer activation uniquely position these genomic elements to finely tune activity-dependent cellular plasticity. Enhancer specificity and modularity can be exploited to gain selective genetic access to specific cell states, and the precise modulation of target gene expression within restricted cellular contexts enabled by targeted enhancer manipulation allows for fine-grained evaluation of gene function. Mounting evidence also suggests that enduring stimulus-induced changes in enhancer states can modify target gene activation upon restimulation, thereby contributing to a form of cell-wide metaplasticity. We advocate for focused exploration of activity-dependent enhancer function to gain new insight into the mechanisms underlying brain plasticity and cognitive dysfunction.

Genetic variants that modify neuroendocrine gene expression and foraging behavior of C. elegans

The molecular mechanisms underlying diversity in animal behavior are not well understood. A major experimental challenge is determining the contribution of genetic variants that affect neuronal gene expression to differences in behavioral traits. In Caenorhabditis elegans, the neuroendocrine transforming growth factor-β ligand, DAF-7, regulates diverse behavioral responses to bacterial food and pathogens. The dynamic neuron-specific expression of daf-7 is modulated by environmental and endogenous bacteria-derived cues. Here, we investigated natural variation in the expression of daf-7 from the ASJ pair of chemosensory neurons. We identified common genetic variants in gap-2, encoding a Ras guanosine triphosphatase (GTPase)-activating protein homologous to mammalian synaptic Ras GTPase-activating protein, which modify daf-7 expression cell nonautonomously and promote exploratory foraging behavior in a partially DAF-7-dependent manner. Our data connect natural variation in neuron-specific gene expression to differences in behavior and suggest that genetic variation in neuroendocrine signaling pathways mediating host-microbe interactions may give rise to diversity in animal behavior.

Genetic Variants That Modify the Neuroendocrine Regulation of Foraging Behavior in C. elegans

The molecular mechanisms underlying diversity in animal behavior are not well understood. A major experimental challenge is determining the contribution of genetic variants that affect neuronal gene expression to differences in behavioral traits. The neuroendocrine TGF-beta ligand, DAF-7, regulates diverse behavioral responses of Caenorhabditis elegans to bacterial food and pathogens. The dynamic neuron-specific expression of daf-7 is modulated by environmental and endogenous bacteria-derived cues. Here, we investigated natural variation in the expression of daf-7 from the ASJ pair of chemosensory neurons and identified common variants in gap-2, encoding a GTPase-Activating Protein homologous to mammalian SynGAP proteins, which modify daf-7 expression cell-non-autonomously and promote exploratory foraging behavior in a DAF-7-dependent manner. Our data connect natural variation in neuron-specific gene expression to differences in behavior and suggest that genetic variation in neuroendocrine signaling pathways mediating host-microbe interactions may give rise to diversity in animal behavior.

Beyond the reference: gene expression variation and transcriptional response to RNA interference in Caenorhabditis elegans

Though natural systems harbor genetic and phenotypic variation, research in model organisms is often restricted to a reference strain. Focusing on a reference strain yields a great depth of knowledge but potentially at the cost of breadth of understanding. Furthermore, tools developed in the reference context may introduce bias when applied to other strains, posing challenges to defining the scope of variation within model systems. Here, we evaluate how genetic differences among 5 wild Caenorhabditis elegans strains affect gene expression and its quantification, in general and after induction of the RNA interference (RNAi) response. Across strains, 34% of genes were differentially expressed in the control condition, including 411 genes that were not expressed at all in at least 1 strain; 49 of these were unexpressed in reference strain N2. Reference genome mapping bias caused limited concern: despite hyperdiverse hotspots throughout the genome, 92% of variably expressed genes were robust to mapping issues. The transcriptional response to RNAi was highly strain- and target-gene-specific and did not correlate with RNAi efficiency, as the 2 RNAi-insensitive strains showed more differentially expressed genes following RNAi treatment than the RNAi-sensitive reference strain. We conclude that gene expression, generally and in response to RNAi, differs across C. elegans strains such that the choice of strain may meaningfully influence scientific inferences. Finally, we introduce a resource for querying gene expression variation in this dataset at https://wildworm.biosci.gatech.edu/rnai/.

Beyond the reference: gene expression variation and transcriptional response to RNAi in C. elegans

A universal feature of living systems is that natural variation in genotype underpins variation in phenotype. Yet, research in model organisms is often constrained to a single genetic background, the reference strain. Further, genomic studies that do evaluate wild strains typically rely on the reference strain genome for read alignment, leading to the possibility of biased inferences based on incomplete or inaccurate mapping; the extent of reference bias can be difficult to quantify. As an intermediary between genome and organismal traits, gene expression is well positioned to describe natural variability across genotypes generally and in the context of environmental responses, which can represent complex adaptive phenotypes. C. elegans sits at the forefront of investigation into small-RNA gene regulatory mechanisms, or RNA interference (RNAi), and wild strains exhibit natural variation in RNAi competency following environmental triggers. Here, we examine how genetic differences among five wild strains affect the C. elegans transcriptome in general and after inducing RNAi responses to two germline target genes. Approximately 34% of genes were differentially expressed across strains; 411 genes were not expressed at all in at least one strain despite robust expression in others, including 49 genes not expressed in reference strain N2. Despite the presence of hyper-diverse hotspots throughout the C. elegans genome, reference mapping bias was of limited concern: over 92% of variably expressed genes were robust to mapping issues. Overall, the transcriptional response to RNAi was strongly strain-specific and highly specific to the target gene, and the laboratory strain N2 was not representative of the other strains. Moreover, the transcriptional response to RNAi was not correlated with RNAi phenotypic penetrance; the two germline RNAi incompetent strains exhibited substantial differential gene expression following RNAi treatment, indicating an RNAi response despite failure to reduce expression of the target gene. We conclude that gene expression, both generally and in response to RNAi, differs across C. elegans strains such that choice of strain may meaningfully influence scientific conclusions. To provide a public, easily accessible resource for querying gene expression variation in this dataset, we introduce an interactive website at https://wildworm.biosci.gatech.edu/rnai/ .

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

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.

Evolution and Diversity of TGF-β Pathways are Linked with Novel Developmental and Behavioral Traits

Transforming growth factor-β (TGF-β) signaling is essential for numerous biologic functions. It is a highly conserved pathway found in all metazoans including the nematode Caenorhabditis elegans, which has also been pivotal in identifying many components. Utilizing a comparative evolutionary approach, we explored TGF-β signaling in nine nematode species and revealed striking variability in TGF-β gene frequency across the lineage. Of the species analyzed, gene duplications in the DAF-7 pathway appear common with the greatest disparity observed in Pristionchus pacificus. Specifically, multiple paralogues of daf-3, daf-4 and daf-7 were detected. To investigate this additional diversity, we induced mutations in 22 TGF-β components and generated corresponding double, triple, and quadruple mutants revealing both conservation and diversification in function. Although the DBL-1 pathway regulating body morphology appears highly conserved, the DAF-7 pathway exhibits functional divergence, notably in some aspects of dauer formation. Furthermore, the formation of the phenotypically plastic mouth in P. pacificus is partially influenced through TGF-β with the strongest effect in Ppa-tag-68. This appears important for numerous processes in P. pacificus but has no known function in C. elegans. Finally, we observe behavioral differences in TGF-β mutants including in chemosensation and the establishment of the P. pacificus kin-recognition signal. Thus, TGF-β signaling in nematodes represents a stochastic genetic network capable of generating novel functions through the duplication and deletion of associated genes.

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.

Dauer Formation in C. elegans Is Modulated through AWC and ASI-Dependent Chemosensation

The perception of our surrounding environment is an amalgamation of stimuli detected by sensory neurons. In Caenorhabditis elegans, olfaction is an essential behavior that determines various behavioral functions such as locomotion, feeding and development. Sensory olfactory cues also initiate downstream neuroendocrine signaling that controls aging, learning, development and reproduction. Innate sensory preferences toward odors (food, pathogens) and reproductive pheromones are modulated by 11 pairs of amphid chemosensory neurons in the head region of C. elegans. Amongst these sensory neurons, the ASI neuron has neuroendocrine functions and secretes neuropeptides, insulin-like peptide (DAF-28) and the TGF-β protein, DAF-7. Its expression levels are modulated by the presence of food (increased levels) and population density (decreased levels). A recent study has shown that EXP-1, an excitatory GABA receptor regulates DAF-7/TGF-β levels and participates in DAF-7/TGF-β-mediated behaviors such as aggregation and bordering. Here, we show that exp-1 mutants show defective responses toward AWC-sensed attractive odors in a non-autonomous manner through ASI neurons. Our dauer experiments reveal that in daf-7 mutants, ASI expressed EXP-1 and STR-2 (a G-protein-coupled receptor; GPCR) that partially maintained reproductive growth of animals. Further, studies suggest that neuronal connections between ASI and AWC neurons are allowed at least partially through ASI secreted DAF-7 or through alternate TGF- β pathway/s regulated by EXP-1 and STR-2. Together, our behavioral, genetic and imaging experiments propose that EXP-1 and STR-2 integrate food cues and allow the animals to display DAF-7/TGF-β neuroendocrine dependent or independent behavioral responses contributing to chemosensensory and developmental plasticity.

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.

Social and sexual behaviors in C. elegans: the first fifty years

For the first 25 years after the landmark 1974 paper that launched the field, most C. elegans biologists were content to think of their subjects as solitary creatures. C. elegans presented no shortage of fascinating biological problems, but some of the features that led Brenner to settle on this species-in particular, its free-living, self-fertilizing lifestyle-also seemed to reduce its potential for interesting social behavior. That perspective soon changed, with the last two decades bringing remarkable progress in identifying and understanding the complex interactions between worms. The growing appreciation that C. elegans behavior can only be meaningfully understood in the context of its ecology and evolution ensures that the coming years will see similarly exciting progress.

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.

The Ancestral Caenorhabditis elegans Cuticle Suppresses rol-1

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

Whole-organism behavioral profiling reveals a role for dopamine in state-dependent motor program coupling in C. elegans

Animal behaviors are commonly organized into long-lasting states that coordinately impact the generation of diverse motor outputs such as feeding, locomotion, and grooming. However, the neural mechanisms that coordinate these distinct motor programs remain poorly understood. Here, we examine how the distinct motor programs of the nematode C. elegans are coupled together across behavioral states. We describe a new imaging platform that permits automated, simultaneous quantification of each of the main C. elegans motor programs over hours or days. Analysis of these whole-organism behavioral profiles shows that the motor programs coordinately change as animals switch behavioral states. Utilizing genetics, optogenetics, and calcium imaging, we identify a new role for dopamine in coupling locomotion and egg-laying together across states. These results provide new insights into how the diverse motor programs throughout an organism are coordinated and suggest that neuromodulators like dopamine can couple motor circuits together in a state-dependent manner.

How Natural Genetic Variation Shapes Behavior

Nervous systems allow animals to acutely respond and behaviorally adapt to changes and recurring patterns in their environment at multiple timescales-from milliseconds to years. Behavior is further shaped at intergenerational timescales by genetic variation, drift, and selection. This sophistication and flexibility of behavior makes it challenging to measure behavior consistently in individual subjects and to compare it across individuals. In spite of these challenges, careful behavioral observations in nature and controlled measurements in the laboratory, combined with modern technologies and powerful genetic approaches, have led to important discoveries about the way genetic variation shapes behavior. A critical mass of genes whose variation is known to modulate behavior in nature is finally accumulating, allowing us to recognize emerging patterns. In this review, we first discuss genetic mapping approaches useful for studying behavior. We then survey how variation acts at different levels-in environmental sensation, in internal neuronal circuits, and outside the nervous system altogether-and then discuss the sources and types of molecular variation linked to behavior and the mechanisms that shape such variation. We end by discussing remaining questions in the field.

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

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.

Shared behavioral mechanisms underlie C. elegans aggregation and swarming

In complex biological systems, simple individual-level behavioral rules can give rise to emergent group-level behavior. While collective behavior has been well studied in cells and larger organisms, the mesoscopic scale is less understood, as it is unclear which sensory inputs and physical processes matter a priori. Here, we investigate collective feeding in the roundworm C. elegans at this intermediate scale, using quantitative phenotyping and agent-based modeling to identify behavioral rules underlying both aggregation and swarming-a dynamic phenotype only observed at longer timescales. Using fluorescence multi-worm tracking, we quantify aggregation in terms of individual dynamics and population-level statistics. Then we use agent-based simulations and approximate Bayesian inference to identify three key behavioral rules for aggregation: cluster-edge reversals, a density-dependent switch between crawling speeds, and taxis towards neighboring worms. Our simulations suggest that swarming is simply driven by local food depletion but otherwise employs the same behavioral mechanisms as the initial aggregation.

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.

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.

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

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

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

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

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

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

Decoding the intensity of sensory input by two glutamate receptors in one C. elegans interneuron

How neurons are capable of decoding stimulus intensity and translate this information into complex behavioral outputs is poorly defined. Here, we demonstrate that the C. elegans interneuron AIB regulates two types of behaviors: reversal initiation and feeding suppression in response to different concentrations of quinine. Low concentrations of quinine are decoded in AIB by a low-threshold, fast-inactivation glutamate receptor GLR-1 and translated into reversal initiation. In contrast, high concentrations of quinine are decoded by a high-threshold, slow-inactivation glutamate receptor GLR-5 in AIB. After activation, GLR-5 evokes sustained Ca²⁺ release from the inositol 1,4,5-trisphosphate (IP₃)-sensitive Ca²⁺ stores and triggers neuropeptide secretion, which in turn activates the downstream neuron RIM and inhibits feeding. Our results reveal that distinct signal patterns in a single interneuron AIB can encode differential behavioral outputs depending on the stimulus intensity, thus highlighting the importance of functional mapping of information propagation at the single-neuron level during connectome construction.

Genetics and Genomics of Social Behavior in a Chicken Model

The identification of genes affecting sociality can give insights into the maintenance and development of sociality and personality. In this study, we used the combination of an advanced intercross between wild and domestic chickens with a combined QTL and eQTL genetical genomics approach to identify genes for social reinstatement, a social and anxiety-related behavior. A total of 24 social reinstatement QTL were identified and overlaid with over 600 eQTL obtained from the same birds using hypothalamic tissue. Correlations between overlapping QTL and eQTL indicated five strong candidate genes, with the gene TTRAP being strongly significantly correlated with multiple aspects of social reinstatement behavior, as well as possessing a highly significant eQTL.

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.

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

BACKGROUND

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

MAIN BODY

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

CONCLUSION

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

Dopamine signaling tunes spatial pattern selectivity in C. elegans

Animals with complex brains can discriminate the spatial arrangement of physical features in the environment. It is unknown whether such sensitivity to spatial patterns can be accomplished in simpler nervous systems that lack long-range sensory modalities such as vision and hearing. Here we show that the nematode Caenorhabditis elegans can discriminate spatial patterns in its surroundings, despite having a nervous system of only 302 neurons. This spatial pattern selectivity requires touch-dependent dopamine signaling, including the mechanosensory TRP-4 channel in dopaminergic neurons and the D2-like dopamine receptor DOP-3. We find that spatial pattern selectivity varies significantly among C. elegans wild isolates. Electrophysiological recordings show that natural variations in TRP-4 reduce the mechanosensitivity of dopaminergic neurons. Polymorphic substitutions in either TRP-4 or DOP-3 alter the selectivity of spatial patterns. Together, these results demonstrate an ancestral role for dopamine signaling in tuning spatial pattern preferences in a simple nervous system.

Natural courtship song variation caused by an intronic retroelement in an ion channel gene

Animal species display enormous variation for innate behaviours, but little is known about how this diversity arose. Here, using an unbiased genetic approach, we map a courtship song difference between wild isolates of Drosophila simulans and Drosophila mauritiana to a 966 base pair region within the slowpoke (slo) locus, which encodes a calcium-activated potassium channel. Using the reciprocal hemizygosity test, we confirm that slo is the causal locus and resolve the causal mutation to the evolutionarily recent insertion of a retroelement in a slo intron within D. simulans. Targeted deletion of this retroelement reverts the song phenotype and alters slo splicing. Like many ion channel genes, slo is expressed widely in the nervous system and influences a variety of behaviours; slo-null males sing little song with severely disrupted features. By contrast, the natural variant of slo alters a specific component of courtship song, illustrating that regulatory evolution of a highly pleiotropic ion channel gene can cause modular changes in behaviour.

Natural Genetic Variation Differentially Affects the Proteome and Transcriptome in Caenorhabditis elegans

Natural genetic variation is the raw material of evolution and influences disease development and progression. An important question is how this genetic variation translates into variation in protein abundance. To analyze the effects of the genetic background on gene and protein expression in the nematode Caenorhabditis elegans, we quantitatively compared the two genetically highly divergent wild-type strains N2 and CB4856. Gene expression was analyzed by microarray assays, and proteins were quantified using stable isotope labeling by amino acids in cell culture. Among all transcribed genes, we found 1,532 genes to be differentially transcribed between the two wild types. Of the total 3,238 quantified proteins, 129 proteins were significantly differentially expressed between N2 and CB4856. The differentially expressed proteins were enriched for genes that function in insulin-signaling and stress-response pathways, underlining strong divergence of these pathways in nematodes. The protein abundance of the two wild-type strains correlates more strongly than protein abundance versus transcript abundance within each wild type. Our findings indicate that in C. elegans only a fraction of the changes in protein abundance can be explained by the changes in mRNA abundance. These findings corroborate with the observations made across species.

Genetics of Intraspecies Variation in Avoidance Behavior Induced by a Thermal Stimulus in Caenorhabditis elegans

Individuals within a species vary in their responses to a wide range of stimuli, partly as a result of differences in their genetic makeup. Relatively little is known about the genetic and neuronal mechanisms contributing to diversity of behavior in natural populations. By studying intraspecies variation in innate avoidance behavior to thermal stimuli in the nematode Caenorhabditis elegans, we uncovered genetic principles of how different components of a behavioral response can be altered in nature to generate behavioral diversity. Using a thermal pulse assay, we uncovered heritable variation in responses to a transient temperature increase. Quantitative trait locus mapping revealed that separate components of this response were controlled by distinct genomic loci. The loci we identified contributed to variation in components of thermal pulse avoidance behavior in an additive fashion. Our results show that the escape behavior induced by thermal stimuli is composed of simpler behavioral components that are influenced by at least six distinct genetic loci. The loci that decouple components of the escape behavior reveal a genetic system that allows independent modification of behavioral parameters. Our work sets the foundation for future studies of evolution of innate behaviors at the molecular and neuronal level.

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.

A Powerful New Quantitative Genetics Platform, Combining Caenorhabditis elegans High-Throughput Fitness Assays with a Large Collection of Recombinant Strains

The genetic variants underlying complex traits are often elusive even in powerful model organisms such as Caenorhabditis elegans with controlled genetic backgrounds and environmental conditions. Two major contributing factors are: (1) the lack of statistical power from measuring the phenotypes of small numbers of individuals, and (2) the use of phenotyping platforms that do not scale to hundreds of individuals and are prone to noisy measurements. Here, we generated a new resource of 359 recombinant inbred strains that augments the existing C. elegans N2xCB4856 recombinant inbred advanced intercross line population. This new strain collection removes variation in the neuropeptide receptor gene npr-1, known to have large physiological and behavioral effects on C. elegans and mitigates the hybrid strain incompatibility caused by zeel-1 and peel-1, allowing for identification of quantitative trait loci that otherwise would have been masked by those effects. Additionally, we optimized highly scalable and accurate high-throughput assays of fecundity and body size using the COPAS BIOSORT large particle nematode sorter. Using these assays, we identified quantitative trait loci involved in fecundity and growth under normal growth conditions and after exposure to the herbicide paraquat, including independent genetic loci that regulate different stages of larval growth. Our results offer a powerful platform for the discovery of the genetic variants that control differences in responses to drugs, other aqueous compounds, bacterial foods, and pathogenic stresses.

A wild C. elegans strain has enhanced epithelial immunity to a natural microsporidian parasite

Microbial pathogens impose selective pressures on their hosts, and combatting these pathogens is fundamental to the propagation of a species. Innate immunity is an ancient system that provides the foundation for pathogen resistance, with epithelial cells in humans increasingly appreciated to play key roles in innate defense. Here, we show that the nematode C. elegans displays genetic variation in epithelial immunity against intestinal infection by its natural pathogen, Nematocida parisii. This pathogen belongs to the microsporidia phylum, which comprises a large phylum of over 1400 species of fungal-related parasites that can infect all animals, including humans, but are poorly understood. Strikingly, we find that a wild C. elegans strain from Hawaii is able to clear intracellular infection by N. parisii, with this ability restricted to young larval animals. Notably, infection of older larvae does not impair progeny production, while infection of younger larvae does. The early-life immunity of Hawaiian larvae enables them to produce more progeny later in life, providing a selective advantage in a laboratory setting—in the presence of parasite it is able to out-compete a susceptible strain in just a few generations. We show that enhanced immunity is dominant to susceptibility, and we use quantitative trait locus mapping to identify four genomic loci associated with resistance. Furthermore, we generate near-isogenic strains to directly demonstrate that two of these loci influence resistance. Thus, our findings show that early-life immunity of C. elegans against microsporidia is a complex trait that enables the host to produce more progeny later in life, likely improving its evolutionary success.

Infectious diseases caused by microbes create some of the strongest forces in evolution, by killing their hosts, and impairing their ability to produce progeny. Microsporidia are very common microbes that cause disease in all animals, including roundworms, insects, fish and people. We investigated microsporidia infection in the roundworm C. elegans, and found that strains from diverse parts of the world have differing levels of resistance against infection. Interestingly, a C. elegans strain from Hawaii can clear infection but only during the earliest stage of life. This resistance appears to be evolutionarily important, because it is during this early stage of life when infection can greatly reduce the number of progeny produced by the host. Consistent with this idea, if the Hawaiian strain is infected when young, it will ultimately produce more progeny than a susceptible strain of C. elegans. We find that this early life resistance of Hawaiian animals is due to a combination of genetic regions, which together provide enhanced immunity against a natural pathogen, thus enabling this strain to have more offspring.

Why we need more ecology for genetic models such as C. elegans

Functional information about the large majority of the genes is still lacking in the classical eukaryotic model species Drosophila melanogaster, Caenorhabditis elegans, and Mus musculus. Because many of these genes are likely to be important in natural settings, considering explicit ecological information should increase our knowledge of gene function. Using C. elegans as an example, we discuss the importance of biotic factors as a driving force in shaping the composition and structure of the nematode genome. We highlight examples for which consideration of ecological information and natural variation have been key to the identification of novel, unexpected gene functions, and use these examples to define future research avenues for the classical genetic model taxa.

Behavioral avoidance of pathogenic bacteria by Caenorhabditis elegans

The simple animal host Caenorhabditis elegans utilizes its nervous system to respond to diverse microbial cues, and can engage in a protective behavioral avoidance response to environmental pathogens. This behavior bears hallmarks of an immune response, with sensors and recognition systems that trigger a protective response following a learning experience. Neuronal circuits required for aversive learning have been defined, revealing conserved signaling modules with dual roles in immunity and neuronal responses to pathogenic bacteria. Identification of natural polymorphisms that modulate avoidance behavior has enabled an improved understanding of host-microbe interactions at the molecular level. We review here these findings and discuss how the microbial cues and host responses defined in C. elegans may provide insight into evolutionarily diverse host-microbe interactions.

Genetic influences on brain gene expression in rats selected for tameness and aggression

Interindividual differences in many behaviors are partly due to genetic differences, but the identification of the genes and variants that influence behavior remains challenging. Here, we studied an F2 intercross of two outbred lines of rats selected for tame and aggressive behavior toward humans for >64 generations. By using a mapping approach that is able to identify genetic loci segregating within the lines, we identified four times more loci influencing tameness and aggression than by an approach that assumes fixation of causative alleles, suggesting that many causative loci were not driven to fixation by the selection. We used RNA sequencing in 150 F2 animals to identify hundreds of loci that influence brain gene expression. Several of these loci colocalize with tameness loci and may reflect the same genetic variants. Through analyses of correlations between allele effects on behavior and gene expression, differential expression between the tame and aggressive rat selection lines, and correlations between gene expression and tameness in F2 animals, we identify the genes Gltscr2, Lgi4, Zfp40, and Slc17a7 as candidate contributors to the strikingly different behavior of the tame and aggressive animals.

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

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

A variant in the neuropeptide receptor npr-1 is a major determinant of Caenorhabditis elegans growth and physiology

The mechanistic basis for how genetic variants cause differences in phenotypic traits is often elusive. We identified a quantitative trait locus in Caenorhabditis elegans that affects three seemingly unrelated phenotypic traits: lifetime fecundity, adult body size, and susceptibility to the human pathogen Staphyloccus aureus. We found a QTL for all three traits arises from variation in the neuropeptide receptor gene npr-1. Moreover, we found that variation in npr-1 is also responsible for differences in 247 gene expression traits. Variation in npr-1 is known to determine whether animals disperse throughout a bacterial lawn or aggregate at the edges of the lawn. We found that the allele that leads to aggregation is associated with reduced growth and reproductive output. The altered gene expression pattern caused by this allele suggests that the aggregation behavior might cause a weak starvation state, which is known to reduce growth rate and fecundity. Importantly, we show that variation in npr-1 causes each of these phenotypic differences through behavioral avoidance of ambient oxygen concentrations. These results suggest that variation in npr-1 has broad pleiotropic effects mediated by altered exposure to bacterial food.

Using the nematode roundworm Caenorhabditis elegans, we identified differences in lifetime fecundity, adult body size, and susceptibility to the human pathogen Staphyloccus aureus between the laboratory strain (N2) from Bristol, England and a wild strain (CB4856) from Hawaii, USA. Using linkage mapping and other genetic tests, we found a QTL for all three traits arises from variation in the neuropeptide receptor gene npr-1. Moreover, we found that variation in npr-1 is also responsive for differences in 247 gene expression traits. Variation in npr-1 is known to determine whether animals disperse throughout a bacterial lawn or aggregate at the edges of the lawn. We found that the allele that leads to aggregation is associated with reduced growth and reproductive output likely caused by a weak chronic starvation state. These results suggest that variation in npr-1 has broad effects on the phenotype of an organism mediated by altered exposure to bacterial food.