Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Free Will as Defined by the Constrained Disorder Principle: a Restricted, Mandatory, Personalized, Regulated Process for Decision-Making

The concept of free will has challenged physicists, biologists, philosophers, and other professionals for decades. The constrained disorder principle (CDP) is a fundamental law that defines systems according to their inherent variability. It provides mechanisms for adapting to dynamic environments. This work examines the CDP's perspective of free will concerning various free will theories. Per the CDP, systems lack intentions, and the "freedom" to select and act is built into their design. The "freedom" is embedded within the response range determined by the boundaries of the systems' variability. This built-in and self-generating mechanism enables systems to cope with perturbations. According to the CDP, neither dualism nor an unknown metaphysical apparatus dictates choices. Brain variability facilitates cognitive adaptation to complex, unpredictable situations across various environments. Human behaviors and decisions reflect an underlying physical variability in the brain and other organs for dealing with unpredictable noises. Choices are not predetermined but reflect the ongoing adaptation processes to dynamic prssu½res. Malfunctions and disease states are characterized by inappropriate variability, reflecting an inability to respond adequately to perturbations. Incorporating CDP-based interventions can overcome malfunctions and disease states and improve decision processes. CDP-based second-generation artificial intelligence platforms improve interventions and are being evaluated to augment personal development, wellness, and health.

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

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

Caenorhabditis elegans immune responses to microsporidia and viruses

The model organism Caenorhabditis elegans is susceptible to infection by obligate intracellular pathogens, specifically microsporidia and viruses. These intracellular pathogens infect intestinal cells, or, for some microsporidia, epidermal cells. Strikingly, intestinal cell infections by viruses or microsporidia trigger a common transcriptional response, activated in part by the ZIP-1 transcription factor. Among the strongest activated genes in this response are ubiquitin-pathway members and members of the pals family, an intriguing gene family with cross-regulations of different members of genomic clusters. Some of the induced genes participate in host defense against the pathogens, for example through ubiquitin-mediated inhibition. Other mechanisms defend the host specifically against viral infections, including antiviral RNA interference and uridylation. These various immune responses are altered by environmental factors and by intraspecific genetic variation of the host. These pathogens were first isolated 15 years ago and much remains to be discovered using C. elegans genetics; also, other intracellular pathogens of C. elegans may yet to be discovered.

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

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

Getting around the roundworms: Identifying knowledge gaps and research priorities for the ascarids

The ascarids are a large group of parasitic nematodes that infect a wide range of animal species. In humans, they cause neglected diseases of poverty; many animal parasites also cause zoonotic infections in people. Control measures include hygiene and anthelmintic treatments, but they are not always appropriate or effective and this creates a continuing need to search for better ways to reduce the human, welfare and economic costs of these infections. To this end, Le Studium Institute of Advanced Studies organized a two-day conference to identify major gaps in our understanding of ascarid parasites with a view to setting research priorities that would allow for improved control. The participants identified several key areas for future focus, comprising of advances in genomic analysis and the use of model organisms, especially Caenorhabditis elegans, a more thorough appreciation of the complexity of host-parasite (and parasite-parasite) communications, a search for novel anthelmintic drugs and the development of effective vaccines. The participants agreed to try and maintain informal links in the future that could form the basis for collaborative projects, and to co-operate to organize future meetings and workshops to promote ascarid research.

The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

Expansion of the split hygromycin toolkit for transgene insertion in Caenorhabditis elegans

Engineered sites for genetic transformation have simplified transgene insertion in Caenorhabditis elegans . These strategies include our split hygromycin system ​(Stevenson et al. 2020)​ which allows for integration-specific selection of transgenes. Here we have expanded the split hygromycin selection system to include two additional chromosomal locations, both of which are permissive for germline expression, as well as engineered landing pads in three additional natural isolates. Corresponding guide and empty repair template plasmids are also available for each of these sites.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

CaeNDR, the Caenorhabditis Natural Diversity Resource

Studies of model organisms have provided important insights into how natural genetic differences shape trait variation. These discoveries are driven by the growing availability of genomes and the expansive experimental toolkits afforded to researchers using these species. For example, Caenorhabditis elegans is increasingly being used to identify and measure the effects of natural genetic variants on traits using quantitative genetics. Since 2016, the C. elegans Natural Diversity Resource (CeNDR) has facilitated many of these studies by providing an archive of wild strains, genome-wide sequence and variant data for each strain, and a genome-wide association (GWA) mapping portal for the C. elegans community. Here, we present an updated platform, the Caenorhabditis Natural Diversity Resource (CaeNDR), that enables quantitative genetics and genomics studies across the three Caenorhabditis species: C. elegans, C. briggsae and C. tropicalis. The CaeNDR platform hosts several databases that are continually updated by the addition of new strains, whole-genome sequence data and annotated variants. Additionally, CaeNDR provides new interactive tools to explore natural variation and enable GWA mappings. All CaeNDR data and tools are accessible through a freely available web portal located at caendr.org.

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

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

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

A dominant dpy-10 co-transformation marker using CRISPR/Cas9 and a linear repair template in Caenorhabditis tropicalis

Caenorhabditis elegans is an excellent genetic model system with a large arsenal of forward and reverse genetic techniques. However, not all approaches are easily ported to related Caenorhabditis species (which are useful for gene conservation and gene pathway evolution studies). For CRISPR/Cas9 genetic editing, an easily screenable and dominant co-transformation marker is required - a secondary mutation that won't impact the phenotype of a desired mutation but is capable of being screened for in heterozygous mutants. We describe here the adaptation of a dominant dumpy/roller CRISPR/Cas9-induced mutation in the C. tropicalis dpy-10 orthologue.

Identification of a spontaneously arising variant affecting thermotaxis behavior in a recombinant inbred Caenorhabditis elegans line

Analyses of the contributions of genetic variants in wild strains to phenotypic differences have led to a more complete description of the pathways underlying cellular functions. Causal loci are typically identified via interbreeding of strains with distinct phenotypes in order to establish recombinant inbred lines (RILs). Since the generation of RILs requires growth for multiple generations, their genomes may contain not only different combinations of parental alleles but also genetic changes that arose de novo during the establishment of these lines. Here, we report that in the course of generating RILs between Caenorhabditis elegans strains that exhibit distinct thermotaxis behavioral phenotypes, we identified spontaneously arising variants in the ttx-1 locus. ttx-1 encodes the terminal selector factor for the AFD thermosensory neurons, and loss-of-function mutations in ttx-1 abolish thermotaxis behaviors. The identified genetic changes in ttx-1 in the RIL are predicted to decrease ttx-1 function in part via specifically affecting a subset of AFD-expressed ttx-1 isoforms. Introduction of the relevant missense mutation in the laboratory C. elegans strain via gene editing recapitulates the thermotaxis behavioral defects of the RIL. Our results suggest that spontaneously occurring genomic changes in RILs may complicate identification of loci contributing to phenotypic variation, but that these mutations may nevertheless lead to the identification of important causal molecules and mechanisms.

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

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Repurposing degradation pathways for modular metabolite biosynthesis in nematodes

Recent studies have revealed that Caenorhabditis elegans and other nematodes repurpose products from biochemical degradation pathways for the combinatorial assembly of complex modular structures that serve diverse signaling functions. Building blocks from neurotransmitter, amino acid, nucleoside and fatty acid metabolism are attached to scaffolds based on the dideoxyhexose ascarylose or glucose, resulting in hundreds of modular ascarosides and glucosides. Genome-wide association studies have identified carboxylesterases as the key enzymes mediating modular assembly, enabling rapid compound discovery via untargeted metabolomics and suggesting that modular metabolite biosynthesis originates from the 'hijacking' of conserved detoxification mechanisms. Modular metabolites thus represent a distinct biosynthetic strategy for generating structural and functional diversity in nematodes, complementing the primarily polyketide synthase- and nonribosomal peptide synthetase-derived universe of microbial natural products. Although many aspects of modular metabolite biosynthesis and function remain to be elucidated, their identification demonstrates how phenotype-driven compound discovery, untargeted metabolomics and genomic approaches can synergize to facilitate the annotation of metabolic dark matter.

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

Diet can alter the cost of resistance to a natural parasite in Caenorhabditis elegans

Resistance to parasites confers a fitness advantage, yet hosts show substantial variation in resistance in natural populations. Evolutionary theory indicates that resistant and susceptible genotypes can coexist if resistance is costly, but there is mixed evidence that resistant individuals have lower fitness in the absence of parasites. One explanation for this discrepancy is that the cost of resistance varies with environmental context. We tested this hypothesis using Caenorhabditis elegans and its natural microsporidian parasite, Nematocida ironsii. We used multiple metrics to compare the fitness of two near-isogenic host genotypes differing at regions associated with resistance to N. ironsii. To quantify the effect of the environment on the cost associated with these known resistance regions, we measured fitness on three microbial diets. We found that the cost of resistance varied with both diet and the measure of fitness. We detected no cost to resistance, irrespective of diet, when fitness was measured as fecundity. However, we detected a cost when fitness was measured in terms of population growth, and the magnitude of this cost varied with diet. These results provide a proof of concept that, by mediating the cost of resistance, environmental context may govern the rate and nature of resistance evolution in heterogeneous environments.

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.

Caenorhabditis nematodes colonize ephemeral resource patches in neotropical forests

Factors shaping the distribution and abundance of species include life-history traits, population structure, and stochastic colonization-extinction dynamics. Field studies of model species groups help reveal the roles of these factors. Species of Caenorhabditis nematodes are highly divergent at the sequence level but exhibit highly conserved morphology, and many of these species live in sympatry on microbe-rich patches of rotten material. Here, we use field experiments and large-scale opportunistic collections to investigate species composition, abundance, and colonization efficiency of Caenorhabditis species in two of the world's best-studied lowland tropical field sites: Barro Colorado Island in Panamá and La Selva in Sarapiquí, Costa Rica. We observed seven species of Caenorhabditis, four of them known only from these collections. We formally describe two species and place them within the Caenorhabditis phylogeny. While these localities contain species from many parts of the phylogeny, both localities were dominated by globally distributed androdiecious species. We found that Caenorhabditis individuals were able to colonize baits accessible only through phoresy and preferentially colonized baits that were in direct contact with the ground. We estimate the number of colonization events per patch to be low.