Recombinational landscape and population genomics of Caenorhabditis elegans

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three-year period, measured various environmental parameters at each sampling site, and whole-genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500-1,500 m) where ambient air temperatures are cool (15-20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper-divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

Chromosome-level reference genomes for two strains of Caenorhabditis briggsae: an improved platform for comparative genomics

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines (RILs) generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing (RNA-seq) data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Post-insemination selection dominates pre-insemination selection in driving rapid evolution of male competitive ability

Sexual reproduction is a complex process that contributes to differences between the sexes and divergence between species. From a male’s perspective, sexual selection can optimize reproductive success by acting on the variance in mating success (pre-insemination selection) as well as the variance in fertilization success (post-insemination selection). The balance between pre- and post-insemination selection has not yet been investigated using a strong hypothesis-testing framework that directly quantifies the effects of post-insemination selection on the evolution of reproductive success. Here we use experimental evolution of a uniquely engineered genetic system that allows sperm production to be turned off and on in obligate male-female populations of Caenorhabditis elegans. We show that enhanced post-insemination competition increases the efficacy of selection and surpasses pre-insemination sexual selection in driving a polygenic response in male reproductive success. We find that after 10 selective events occurring over 30 generations post-insemination selection increased male reproductive success by an average of 5- to 7-fold. Contrary to expectation, enhanced pre-insemination competition hindered selection and slowed the rate of evolution. Furthermore, we found that post-insemination selection resulted in a strong polygenic response at the whole-genome level. Our results demonstrate that post-insemination sexual selection plays a critical role in the rapid optimization of male reproductive fitness. Therefore, explicit consideration should be given to post-insemination dynamics when considering the population effects of sexual selection.

Some of the most dramatic and diverse phenotypes observed in nature––such as head-butting in wild sheep and the elaborate tails of peacocks––are sexually dimorphic. These remarkable phenotypes are a result of sexual selection optimizing reproductive success in females and males independently. For males, total reproductive success is comprised of winning a mating event and then translating that mating event into a fertilization event. Therefore, to understand not only how male reproductive success is comprised, but also how it evolves, we must examine the interaction between pre- and post-insemination sexual selection. We combine environmentally-inducible control of sperm production within a highly reproducible factorial experimental evolution design to directly quantify the contribution of post-insemination selection to male reproductive evolution. We demonstrate that enhanced sperm competition increases the efficacy of selection and enhances the rate of male evolution. Alternatively, we show that enhanced pre-insemination competition slows the evolutionary rate. Using whole-genome approaches, we identify over 60 genes that contribute to male fertilization success. Brought together, our new approaches and results demonstrate that the unseen world of molecular interactions occurring during post-insemination are as fundamentally important as pre-mating factors.

Microbial protection favors parasite tolerance and alters host-parasite coevolutionary dynamics

Coevolution between hosts and parasites is a major driver of rapid evolutionary change¹ and diversification.^2,3 However, direct antagonistic interactions between hosts and parasites could be disrupted⁴ when host microbiota form a line of defense, a phenomenon widespread across animal and plant species.^5,6 By suppressing parasite infection, protective microbiota could reduce the need for host-based defenses and favor host support for microbiota colonization,⁶ raising the possibility that the microbiota can alter host-parasite coevolutionary patterns and processes.⁷ Here, using an experimental evolution approach, we co-passaged populations of nematode host (Caenorhabditis elegans) and parasites (Staphylococcus aureus) when hosts were colonized (or not) by protective bacteria (Enterococcus faecalis). We found that microbial protection during coevolution resulted in the evolution of host mortality tolerance—higher survival following parasite infection—and in parasites adapting to microbial defenses. Compared to unprotected host-parasite coevolution, the protected treatment was associated with reduced dominance of fluctuating selection dynamics in host populations. No differences in host recombination rate or genetic diversity were detected. Genomic divergence was observed between parasite populations coevolved in protected and unprotected hosts. These findings indicate that protective host microbiota can determine the evolution of host defense strategies and shape host-parasite coevolutionary dynamics.

•

Microbial protection resulted in the evolution of host mortality tolerance

•

Parasites adapted to counter microbial defenses within hosts

•

Protective microbes reduced fluctuating selection dynamics

•

Microbial protection did not impact host genetic diversity or recombination rates

Rapid Isolation of Wild Nematodes by Baermann Funnel

Beyond being robust experimental model organisms, Caenorhabditis elegans and its relatives are also real animals that live in nature. Studies of wild nematodes in their natural environments are valuable for understanding many aspects of biology, including the selective regimes in which distinctive genomic and phenotypic characters evolve, the genetic basis for complex trait variation, and the natural genetic diversity fundamental to all animal populations. This manuscript describes a simple and efficient method for extracting nematodes from their natural substrates, including rotting fruits, flowers, fungi, leaf litter, and soil. The Baermann funnel method, a classical nematology technique, selectively isolates active nematodes from their substrates. Because it recovers nearly all active worms from the sample, the Baermann funnel technique allows for the recovery of rare and slow-growing genotypes that co-occur with abundant and fast-growing genotypes, which might be missed in extraction methods that involve multiple generations of reproduction. The technique is also well suited to addressing metagenetic, population-genetic, and ecological questions. It captures the entire population in a sample simultaneously, allowing an unbiased view of the natural distribution of ages, sexes, and genotypes. The protocol allows for deployment at scale in the field, rapidly converting substrates into worm plates, and the authors have validated it through fieldwork on multiple continents.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

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.

Mutation rate and spectrum in obligately outcrossing Caenorhabditis elegans mutation accumulation lines subjected to RNAi-induced knockdown of the mismatch repair gene msh-2

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutSα homolog, msh-2, in Caenorhabditis elegans male–female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male–female background substantially increased the frequency of nuclear base substitutions (∼23×) and small indels (∼328×) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

The distribution of mutational effects on fitness in Caenorhabditis elegans inferred from standing genetic variation

The distribution of fitness effects (DFE) for new mutations is one of the most theoretically important but difficult to estimate properties in population genetics. A crucial challenge to inferring the DFE from natural genetic variation is the sensitivity of the site frequency spectrum to factors like population size change, population substructure, genome structure, and nonrandom mating. Although inference methods aim to control for population size changes, the influence of nonrandom mating remains incompletely understood, despite being a common feature of many species. We report the DFE estimated from 326 genomes of Caenorhabditis elegans, a nematode roundworm with a high rate of self-fertilization. We evaluate the robustness of DFE inferences using simulated data that mimics the genomic structure and reproductive life history of C. elegans. Our observations demonstrate how the combined influence of self-fertilization, genome structure, and natural selection on linked sites can conspire to compromise estimates of the DFE from extant polymorphisms with existing methods. These factors together tend to bias inferences toward weakly deleterious mutations, making it challenging to have full confidence in the inferred DFE of new mutations as deduced from standing genetic variation in species like C. elegans. Improved methods for inferring the DFE are needed to appropriately handle strong linked selection and selfing. These results highlight the importance of understanding the combined effects of processes that can bias our interpretations of evolution in natural populations.

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.

Analysis of meiosis in Pristionchus pacificus reveals plasticity in homolog pairing and synapsis in the nematode lineage

Meiosis is conserved across eukaryotes yet varies in the details of its execution. Here we describe a new comparative model system for molecular analysis of meiosis, the nematode Pristionchus pacificus, a distant relative of the widely studied model organism Caenorhabditis elegans. P. pacificus shares many anatomical and other features that facilitate analysis of meiosis in C. elegans. However, while C. elegans has lost the meiosis-specific recombinase Dmc1 and evolved a recombination-independent mechanism to synapse its chromosomes, P. pacificus expresses both DMC-1 and RAD-51. We find that SPO-11 and DMC-1 are required for stable homolog pairing, synapsis, and crossover formation, while RAD-51 is dispensable for these key meiotic processes. RAD-51 and DMC-1 localize sequentially to chromosomes during meiotic prophase and show nonoverlapping functions. We also present a new genetic map for P. pacificus that reveals a crossover landscape very similar to that of C. elegans, despite marked divergence in the regulation of synapsis and crossing-over between these lineages.

Natural variation in fecundity is correlated with species-wide levels of divergence in Caenorhabditis elegans

Life history traits underlie the fitness of organisms and are under strong natural selection. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g., a selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of a species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fecundity of 121 wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions have significantly higher fecundity than strains without evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fecundity variation. In addition, we mapped previous fecundity data from wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines that were grown in various conditions and detected eight QTL using GWA and linkage mappings. These QTL show the genetic complexity of fecundity across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fecundity than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that caused increased fecundity could underlie the selective sweeps during the worldwide expansion of C. elegans.

Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

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.

Balancing selection maintains hyper-divergent haplotypes in Caenorhabditis elegans

Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing decreases the effective population size, effective recombination rate and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis, selfing has evolved at least three times, and all three species, including the model organism Caenorhabditis elegans, show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long-read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example of how species can avoid the evolutionary dead end associated with selfing.

Let's get physical - mechanisms of crossover interference

The formation of crossovers between homologous chromosomes is key to sexual reproduction. In most species, crossovers are spaced further apart than would be expected if they formed independently, a phenomenon termed crossover interference. Despite more than a century of study, the molecular mechanisms implementing crossover interference remain a subject of active debate. Recent findings of how signaling proteins control the formation of crossovers and about the interchromosomal interface in which crossovers form offer new insights into this process. In this Review, we present a cell biological and biophysical perspective on crossover interference, summarizing the evidence that links interference to the spatial, dynamic, mechanical and molecular properties of meiotic chromosomes. We synthesize this physical understanding in the context of prevailing mechanistic models that aim to explain how crossover interference is implemented.

Punctuated Loci on Chromosome IV Determine Natural Variation in Orsay Virus Susceptibility of Caenorhabditis elegans Strains Bristol N2 and Hawaiian CB4856

Host-pathogen interactions play a major role in evolutionary selection and shape natural genetic variation. The genetically distinct Caenorhabditis elegans strains, Bristol N2 and Hawaiian CB4856, are differentially susceptible to the Orsay virus (OrV). Here, we report the dissection of the genetic architecture of susceptibility to OrV infection. We compare OrV infection in the relatively resistant wild-type CB4856 strain to the more susceptible canonical N2 strain. To gain insight into the genetic architecture of viral susceptibility, 52 fully sequenced recombinant inbred lines (CB4856 × N2 RILs) were exposed to OrV. This led to the identification of two loci on chromosome IV associated with OrV resistance. To verify the two loci and gain additional insight into the genetic architecture controlling virus infection, introgression lines (ILs) that together cover chromosome IV, were exposed to OrV. Of the 27 ILs used, 17 had an CB4856 introgression in an N2 background, and 10 had an N2 introgression in a CB4856 background. Infection of the ILs confirmed and fine-mapped the locus underlying variation in OrV susceptibility, and we found that a single nucleotide polymorphism in cul-6 may contribute to the difference in OrV susceptibility between N2 and CB4856. An allele swap experiment showed the strain CB4856 became as susceptible as the N2 strain by having an N2 cul-6 allele, although having the CB4856 cul-6 allele did not increase resistance in N2. In addition, we found that multiple strains with nonoverlapping introgressions showed a distinct infection phenotype from the parental strain, indicating that there are punctuated locations on chromosome IV determining OrV susceptibility. Thus, our findings reveal the genetic complexity of OrV susceptibility in C. elegans and suggest that viral susceptibility is governed by multiple genes.IMPORTANCE Genetic variation determines the viral susceptibility of hosts. Yet, pinpointing which genetic variants determine viral susceptibility remains challenging. Here, we have exploited the genetic tractability of the model organism Caenorhabditis elegans to dissect the genetic architecture of Orsay virus infection. Our results provide novel insight into natural determinants of Orsay virus infection.

HIM-17 regulates the position of recombination events and GSP-1/2 localization to establish short arm identity on bivalents in meiosis

The position of recombination events established along chromosomes in early prophase I and the chromosome remodeling that takes place in late prophase I are intrinsically linked steps of meiosis that need to be tightly regulated to ensure accurate chromosome segregation and haploid gamete formation. Here, we show that RAD-51 foci, which form at the sites of programmed meiotic DNA double-strand breaks (DSBs), exhibit a biased distribution toward off-centered positions along the chromosomes in wild-type Caenorhabditis elegans, and we identify two meiotic roles for chromatin-associated protein HIM-17 that ensure normal chromosome remodeling in late prophase I. During early prophase I, HIM-17 regulates the distribution of DSB-dependent RAD-51 foci and crossovers on chromosomes, which is critical for the formation of distinct chromosome subdomains (short and long arms of the bivalents) later during chromosome remodeling. During late prophase I, HIM-17 promotes the normal expression and localization of protein phosphatases GSP-1/2 to the surface of the bivalent chromosomes and may promote GSP-1 phosphorylation, thereby antagonizing Aurora B kinase AIR-2 loading on the long arms and preventing premature loss of sister chromatid cohesion. We propose that HIM-17 plays distinct roles at different stages during meiotic progression that converge to promote normal chromosome remodeling and accurate chromosome segregation.

Whole-organism eQTL mapping at cellular resolution with single-cell sequencing

Genetic regulation of gene expression underlies variation in disease risk and other complex traits. The effect of expression quantitative trait loci (eQTLs) varies across cell types; however, the complexity of mammalian tissues makes studying cell-type eQTLs highly challenging. We developed a novel approach in the model nematode Caenorhabditis elegans that uses single-cell RNA sequencing to map eQTLs at cellular resolution in a single one-pot experiment. We mapped eQTLs across cell types in an extremely large population of genetically distinct C. elegans individuals. We found cell-type-specific trans eQTL hotspots that affect the expression of core pathways in the relevant cell types. Finally, we found single-cell-specific eQTL effects in the nervous system, including an eQTL with opposite effects in two individual neurons. Our results show that eQTL effects can be specific down to the level of single cells.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Synaptonemal Complex dimerization regulates chromosome alignment and crossover patterning in meiosis

During sexual reproduction the parental homologous chromosomes find each other (pair) and align along their lengths by integrating local sequence homology with large-scale contiguity, thereby allowing for precise exchange of genetic information. The Synaptonemal Complex (SC) is a conserved zipper-like structure that assembles between the homologous chromosomes, bringing them together and regulating exchanges between them. However, the molecular mechanisms by which the SC carries out these functions remain poorly understood. Here we isolated and characterized two mutations in the dimerization interface in the middle of the SC zipper in C. elegans. The mutations perturb both chromosome alignment and the regulation of genetic exchanges. Underlying the chromosome-scale phenotypes are distinct alterations to the way SC subunits interact with one another. We propose a model whereby the SC brings homologous chromosomes together through two activities: obligate zipping that prevents assembly on unpaired chromosomes; and a tendency to extend pairing interactions along the entire length of the chromosomes.

Gene-level quantitative trait mapping in Caenorhabditis elegans

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

A single-nucleotide change underlies the genetic assimilation of a plastic trait

Genetic assimilation-the evolutionary process by which an environmentally induced phenotype is made constitutive-represents a fundamental concept in evolutionary biology. Thought to reflect adaptive phenotypic plasticity, matricidal hatching in nematodes is triggered by maternal nutrient deprivation to allow for protection or resource provisioning of offspring. Here, we report natural Caenorhabditis elegans populations harboring genetic variants expressing a derived state of near-constitutive matricidal hatching. These variants exhibit a single amino acid change (V530L) in KCNL-1, a small-conductance calcium-activated potassium channel subunit. This gain-of-function mutation causes matricidal hatching by strongly reducing the sensitivity to environmental stimuli triggering egg-laying. We show that reestablishing the canonical KCNL-1 protein in matricidal isolates is sufficient to restore canonical egg-laying. While highly deleterious in constant food environments, KCNL-1 V530L is maintained under fluctuating resource availability. A single point mutation can therefore underlie the genetic assimilation-by either genetic drift or selection-of an ancestrally plastic trait.

Ubiquitous Selfish Toxin-Antidote Elements in Caenorhabditis Species

Toxin-antidote elements (TAs) are selfish genetic dyads that spread in populations by selectively killing non-carriers. TAs are common in prokaryotes, but very few examples are known in animals. Here, we report the discovery of maternal-effect TAs in both C. tropicalis and C. briggsae, two distant relatives of C. elegans. In C. tropicalis, multiple TAs combine to cause a striking degree of intraspecific incompatibility: five elements reduce the fitness of >70% of the F₂ hybrid progeny of two Caribbean isolates. We identified the genes underlying one of the novel TAs, slow-1/grow-1, and found that its toxin, slow-1, is homologous to nuclear hormone receptors. Remarkably, although previously known TAs act during embryonic development, maternal loading of slow-1 in oocytes specifically slows down larval development, delaying the onset of reproduction by several days. Finally, we found that balancing selection acting on linked, conflicting TAs hampers their ability to spread in populations, leading to more stable genetic incompatibilities. Our findings indicate that TAs are widespread in Caenorhabditis species and target a wide range of developmental processes and that antagonism between them may cause lasting incompatibilities in natural populations. We expect that similar phenomena exist in other animal species.

Meiosis Progression and Recombination in Holocentric Plants: What Is Known?

Differently from the common monocentric organization of eukaryotic chromosomes, the so-called holocentric chromosomes present many centromeric regions along their length. This chromosomal organization can be found in animal and plant lineages, whose distribution suggests that it has evolved independently several times. Holocentric chromosomes present an advantage: even broken chromosome parts can be correctly segregated upon cell division. However, the evolution of holocentricity brought about consequences to nuclear processes and several adaptations are necessary to cope with this new organization. Centromeres of monocentric chromosomes are involved in a two-step cohesion release during meiosis. To deal with that holocentric lineages developed different adaptations, like the chromosome remodeling strategy in Caenorhabditis elegans or the inverted meiosis in plants. Furthermore, the frequency of recombination at or around centromeres is normally very low and the presence of centromeric regions throughout the entire length of the chromosomes could potentially pose a problem for recombination in holocentric organisms. However, meiotic recombination happens, with exceptions, in those lineages in spite of their holocentric organization suggesting that the role of centromere as recombination suppressor might be altered in these lineages. Most of the available information about adaptations to meiosis in holocentric organisms is derived from the animal model C. elegans. As holocentricity evolved independently in different lineages, adaptations observed in C. elegans probably do not apply to other lineages and very limited research is available for holocentric plants. Currently, we still lack a holocentric model for plants, but good candidates may be found among Cyperaceae, a large angiosperm family. Besides holocentricity, chiasmatic and achiasmatic inverted meiosis are found in the family. Here, we introduce the main concepts of meiotic constraints and adaptations with special focus in meiosis progression and recombination in holocentric plants. Finally, we present the main challenges and perspectives for future research in the field of chromosome biology and meiosis in holocentric plants.

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.

Toxin-Antidote Elements Across the Tree of Life

In life's constant battle for survival, it takes one to kill but two to conquer. Toxin-antitoxin or toxin-antidote (TA) elements are genetic dyads that cheat the laws of inheritance to guarantee their transmission to the next generation. This seemingly simple genetic arrangement—a toxin linked to its antidote—is capable of quickly spreading and persisting in natural populations. TA elements were first discovered in bacterial plasmids in the 1980s and have recently been characterized in fungi, plants, and animals, where they underlie genetic incompatibilities and sterility in crosses between wild isolates. In this review, we provide a unified view of TA elements in both prokaryotic and eukaryotic organisms and highlight their similarities and differences at the evolutionary, genetic, and molecular levels. Finally, we propose several scenarios that could explain the paradox of the evolutionary origin of TA elements and argue that these elements may be key evolutionary players and that the full scope of their roles is only beginning to be uncovered.

Phosphoregulation of HORMA domain protein HIM-3 promotes asymmetric synaptonemal complex disassembly in meiotic prophase in Caenorhabditis elegans

In the two cell divisions of meiosis, diploid genomes are reduced into complementary haploid sets through the discrete, two-step removal of chromosome cohesion, a task carried out in most eukaryotes by protecting cohesion at the centromere until the second division. In eukaryotes without defined centromeres, however, alternative strategies have been innovated. The best-understood of these is found in the nematode Caenorhabditis elegans: after the single off-center crossover divides the chromosome into two segments, or arms, several chromosome-associated proteins or post-translational modifications become specifically partitioned to either the shorter or longer arm, where they promote the correct timing of cohesion loss through as-yet unknown mechanisms. Here, we investigate the meiotic axis HORMA-domain protein HIM-3 and show that it becomes phosphorylated at its C-terminus, within the conserved “closure motif” region bound by the related HORMA-domain proteins HTP-1 and HTP-2. Binding of HTP-2 is abrogated by phosphorylation of the closure motif in in vitro assays, strongly suggesting that in vivo phosphorylation of HIM-3 likely modulates the hierarchical structure of the chromosome axis. Phosphorylation of HIM-3 only occurs on synapsed chromosomes, and similarly to other previously-described phosphorylated proteins of the synaptonemal complex, becomes restricted to the short arm after designation of crossover sites. Regulation of HIM-3 phosphorylation status is required for timely disassembly of synaptonemal complex central elements from the long arm, and is also required for proper timing of HTP-1 and HTP-2 dissociation from the short arm. Phosphorylation of HIM-3 thus plays a role in establishing the identity of short and long arms, thereby contributing to the robustness of the two-step chromosome segregation.

To segregate properly in meiosis, cohesion between replicated chromosomes must remain after the first meiotic cell division, so chromosomes can be held together until they finally separate in the second division. While the majority of organisms use centromeres to protect chromosome cohesion in the first division, the nematode worm C. elegans, which lacks single centromeres, instead protects cohesion only on a segment of the chromosome known as the “long arm”. The long arm (and its complement, the short arm) are known to accumulate specific proteins and protein modifications, but it is not known how the short and long arms are first distinguished, nor how their separate functions are carried out. We report here that the chromosome axis protein HIM-3 and its modification by phosphorylation is important for ensuring the robust establishment of short and long arm functions. We show that phosphorylated HIM-3 partitions to the short arms after crossover recombination sites are designated, and HIM-3 mutants that mimic constitutive phosphorylation delay the normal establishment of the two complementary arm domains. Our findings reveal another layer of regulation to an outstanding mystery in chromosome biology.

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.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Meiotic Double-Strand Break Processing and Crossover Patterning Are Regulated in a Sex-Specific Manner by BRCA1-BARD1 in Caenorhabditis elegans

Meiosis is regulated in a sex-specific manner to produce two distinct gametes, sperm and oocytes, for sexual reproduction. To determine how meiotic recombination is regulated in spermatogenesis, we analyzed the meiotic phenotypes of mutants in the tumor suppressor E3 ubiquitin ligase BRC-1-BRD-1 complex in Caenorhabditis elegans male meiosis. Unlike in mammals, this complex is not required for meiotic sex chromosome inactivation, the process whereby hemizygous sex chromosomes are transcriptionally silenced. Interestingly, brc-1 and brd-1 mutants show meiotic recombination phenotypes that are largely opposing to those previously reported for female meiosis. Fewer meiotic recombination intermediates marked by the recombinase RAD-51 were observed in brc-1 and brd-1 mutants, and the reduction in RAD-51 foci could be suppressed by mutation of nonhomologous-end-joining proteins. Analysis of GFP::RPA-1 revealed fewer foci in the brc-1brd-1 mutant and concentration of BRC-1-BRD-1 to sites of meiotic recombination was dependent on DNA end resection, suggesting that the complex regulates the processing of meiotic double-strand breaks to promote repair by homologous recombination. Further, BRC-1-BRD-1 is important to promote progeny viability when male meiosis is perturbed by mutations that block the pairing and synapsis of different chromosome pairs, although the complex is not required to stabilize the RAD-51 filament as in female meiosis under the same conditions. Analyses of crossover designation and formation revealed that BRC-1-BRD-1 inhibits supernumerary COs when meiosis is perturbed. Together, our findings suggest that BRC-1-BRD-1 regulates different aspects of meiotic recombination in male and female meiosis.

Stoichiometric interactions explain spindle dynamics and scaling across 100 million years of nematode evolution

The spindle shows remarkable diversity, and changes in an integrated fashion, as cells vary over evolution. Here, we provide a mechanistic explanation for variations in the first mitotic spindle in nematodes. We used a combination of quantitative genetics and biophysics to rule out broad classes of models of the regulation of spindle length and dynamics, and to establish the importance of a balance of cortical pulling forces acting in different directions. These experiments led us to construct a model of cortical pulling forces in which the stoichiometric interactions of microtubules and force generators (each force generator can bind only one microtubule), is key to explaining the dynamics of spindle positioning and elongation, and spindle final length and scaling with cell size. This model accounts for variations in all the spindle traits we studied here, both within species and across nematode species spanning over 100 million years of evolution.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

Natural variation in a glucuronosyltransferase modulates propionate sensitivity in a C. elegans propionic acidemia model

Mutations in human metabolic genes can lead to rare diseases known as inborn errors of human metabolism. For instance, patients with loss-of-function mutations in either subunit of propionyl-CoA carboxylase suffer from propionic acidemia because they cannot catabolize propionate, leading to its harmful accumulation. Both the penetrance and expressivity of metabolic disorders can be modulated by genetic background. However, modifiers of these diseases are difficult to identify because of the lack of statistical power for rare diseases in human genetics. Here, we use a model of propionic acidemia in the nematode Caenorhabditis elegans to identify genetic modifiers of propionate sensitivity. Using genome-wide association (GWA) mapping across wild strains, we identify several genomic regions correlated with reduced propionate sensitivity. We find that natural variation in the putative glucuronosyltransferase GLCT-3, a homolog of human B3GAT, partly explains differences in propionate sensitivity in one of these genomic intervals. We demonstrate that loss-of-function alleles in glct-3 render the animals less sensitive to propionate. Additionally, we find that C. elegans has an expansion of the glct gene family, suggesting that the number of members of this family could influence sensitivity to excess propionate. Our findings demonstrate that natural variation in genes that are not directly associated with propionate breakdown can modulate propionate sensitivity. Our study provides a framework for using C. elegans to characterize the contributions of genetic background in models of human inborn errors in metabolism.

Caenorhabditis elegans Deficient in DOT-1.1 Exhibit Increases in H3K9me2 at Enhancer and Certain RNAi-Regulated Regions

The methylation of histone H3 at lysine 79 is a feature of open chromatin. It is deposited by the conserved histone methyltransferase DOT1. Recently, DOT1 localization and H3K79 methylation (H3K79me) have been correlated with enhancers in C. elegans and mammalian cells. Since earlier research implicated H3K79me in preventing heterochromatin formation both in yeast and leukemic cells, we sought to inquire whether a H3K79me deficiency would lead to higher levels of heterochromatic histone modifications, specifically H3K9me2, at developmental enhancers in C. elegans. Therefore, we used H3K9me2 ChIP-seq to compare its abundance in control and dot-1.1 loss-of-function mutant worms, as well as in rde-4; dot-1.1 and rde-1; dot-1.1 double mutants. The rde-1 and rde-4 genes are components of the RNAi pathway in C. elegans, and RNAi is known to initiate H3K9 methylation in many organisms, including C. elegans. We have previously shown that dot-1.1(-) lethality is rescued by rde-1 and rde-4 loss-of-function. Here we found that H3K9me2 was elevated in enhancer, but not promoter, regions bound by the DOT-1.1/ZFP-1 complex in dot-1.1(-) worms. We also found increased H3K9me2 at genes targeted by the ALG-3/4-dependent small RNAs and repeat regions. Our results suggest that ectopic H3K9me2 in dot-1.1(-) could, in some cases, be induced by small RNAs.

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.

The nematode Caenorhabditis elegans and the terrestrial isopod Porcellio scaber likely interact opportunistically

Phoresy is a behavior in which an organism, the phoront, travels from one location to another by 'hitching a ride' on the body of a host as it disperses. Some phoronts are generalists, taking advantage of any available host. Others are specialists and travel only when specific hosts are located using chemical cues to identify and move (chemotax) toward the preferred host. Free-living nematodes, like Caenorhabditis elegans, are often found in natural environments that contain terrestrial isopods and other invertebrates. Additionally, the C. elegans wild strain PB306 was isolated associated with the isopod Porcellio scaber. However, it is currently unclear if C. elegans is a phoront of terrestrial isopods, and if so, whether it is a specialist, generalist, or developmental stage-specific combination of both strategies. Because the relevant chemical stimuli might be secreted compounds or volatile odorants, we used different types of chemotaxis assays across diverse extractions of compounds or odorants to test whether C. elegans is attracted to P. scaber. We show that two different strains-the wild isolate PB306 and the laboratory-adapted strain N2 -are not attracted to P. scaber during either the dauer or adult life stages. Our results indicate that C. elegans was not attracted to chemical compounds or volatile odorants from P. scaber, providing valuable empirical evidence to suggest that any associations between these two species are likely opportunistic rather than specific phoresy.

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.

Crossover Position Drives Chromosome Remodeling for Accurate Meiotic Chromosome Segregation

Interhomolog crossovers (COs) are a prerequisite for achieving accurate chromosome segregation during meiosis [1, 2]. COs are not randomly positioned, occurring at distinct genomic intervals during meiosis in all species examined [3-10]. The role of CO position as a major determinant of accurate chromosome segregation has not been previously directly analyzed in a metazoan. Here, we use spo-11 mutants, which lack endogenous DNA double-strand breaks (DSBs), to induce a single DSB by Mos1 transposon excision at defined chromosomal locations in the C. elegans germline and show that the position of the resulting CO directly affects the formation of distinct chromosome subdomains during meiotic chromosome remodeling. CO formation in the typically CO-deprived center region of autosomes leads to premature loss of sister chromatid cohesion and chromosome missegregation, whereas COs at an off-centered position, as in wild type, can result in normal remodeling and accurate segregation. Ionizing radiation (IR)-induced DSBs lead to the same outcomes, and modeling of IR dose-response reveals that the CO-unfavorable center region encompasses up to 6% of the total chromosome length. DSBs proximal to telomeres rarely form COs, likely because of formation of unstable recombination intermediates that cannot be sustained as chiasmata until late prophase. Our work supports a model in which regulation of CO position early in meiotic prophase is required for proper designation of chromosome subdomains and normal chromosome remodeling in late meiotic prophase I, resulting in accurate chromosome segregation and providing a mechanism to prevent aneuploid gamete formation.

Chromosome-Level Assembly of the Caenorhabditis remanei Genome Reveals Conserved Patterns of Nematode Genome Organization

The nematode Caenorhabditis elegans is one of the key model systems in biology, including possessing the first fully assembled animal genome. Whereas C. elegans is a self-reproducing hermaphrodite with fairly limited within-population variation, its relative C. remanei is an outcrossing species with much more extensive genetic variation, making it an ideal parallel model system for evolutionary genetic investigations. Here, we greatly improve on previous assemblies by generating a chromosome-level assembly of the entire C. remanei genome (124.8 Mb of total size) using long-read sequencing and chromatin conformation capture data. Like other fully assembled genomes in the genus, we find that the C. remanei genome displays a high degree of synteny with C. elegans despite multiple within-chromosome rearrangements. Both genomes have high gene density in central regions of chromosomes relative to chromosome ends and the opposite pattern for the accumulation of repetitive elements. C. elegans and C. remanei also show similar patterns of interchromosome interactions, with the central regions of chromosomes appearing to interact with one another more than the distal ends. The new C. remanei genome presented here greatly augments the use of the Caenorhabditis as a platform for comparative genomics and serves as a basis for molecular population genetics within this highly diverse species.

The Genome of Caenorhabditis bovis

The free-living nematode Caenorhabditis elegans is a key laboratory model for metazoan biology. C. elegans has also become a model for parasitic nematodes despite being only distantly related to most parasitic species. All of the ∼65 Caenorhabditis species currently in culture are free-living, with most having been isolated from decaying plant or fungal matter. Caenorhabditis bovis is a particularly unusual species that has been isolated several times from the inflamed ears of Zebu cattle in Eastern Africa, where it is associated with the disease bovine parasitic otitis. C. bovis is therefore of particular interest to researchers interested in the evolution of nematode parasitism. However, as C. bovis is not in laboratory culture, it remains little studied. Here, by sampling livestock markets and slaughterhouses in Western Kenya, we successfully reisolated C. bovis from the ear of adult female Zebu. We sequenced the genome of C. bovis using the Oxford Nanopore MinION platform in a nearby field laboratory and used the data to generate a chromosome-scale draft genome sequence. We exploited this draft genome sequence to reconstruct the phylogenetic relationships of C. bovis to other Caenorhabditis species and reveal the changes in genome size and content that have occurred during its evolution. We also identified expansions in several gene families that have been implicated in parasitism in other nematode species. The high-quality draft genome and our analyses thereof represent a significant advancement in our understanding of this unusual Caenorhabditis species.

The whole transcriptome regulation as a function of mitochondrial polymorphisms and aging in Caenorhabditis elegans

Recently, mitochondrial-nuclear interaction in aging has been widely studied. However, the nuclear genome controlled by natural mitochondrial variations that influence aging has not been comprehensively understood so far. We hypothesized that mitochondrial polymorphisms could play critical roles in the aging process, probably by regulation of the whole-transcriptome expression. Our results showed that mitochondria polymorphisms not only decreased the mitochondrial mass but also miRNA, lncRNA, mRNA, circRNA and metabolite profiles. Furthermore, most genes that are associated with mitochondria show age-related expression features (P = 3.58E-35). We also constructed a differentially expressed circRNA-lncRNA-miRNA-mRNA regulatory network and a ceRNA network affected by the mitochondrial variations. In addition, Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the genes affected by the mitochondrial variation were enriched in metabolic activity. We finally constructed a multi-level regulatory network with aging which affected by the mitochondrial variation in Caenorhabditis elegans. The interactions between these genes and metabolites have great values for further aging research. In sum, our findings provide new evidence for understanding the molecular mechanisms of how mitochondria influence aging.

Transcription-independent TFIIIC-bound sites cluster near heterochromatin boundaries within lamina-associated domains in C. elegans

BACKGROUND

Chromatin organization is central to precise control of gene expression. In various eukaryotic species, domains of pervasive cis-chromatin interactions demarcate functional domains of the genomes. In nematode Caenorhabditis elegans, however, pervasive chromatin contact domains are limited to the dosage-compensated sex chromosome, leaving the principle of C. elegans chromatin organization unclear. Transcription factor III C (TFIIIC) is a basal transcription factor complex for RNA polymerase III, and is implicated in chromatin organization. TFIIIC binding without RNA polymerase III co-occupancy, referred to as extra-TFIIIC binding, has been implicated in insulating active and inactive chromatin domains in yeasts, flies, and mammalian cells. Whether extra-TFIIIC sites are present and contribute to chromatin organization in C. elegans remains unknown.

RESULTS

We identified 504 TFIIIC-bound sites absent of RNA polymerase III and TATA-binding protein co-occupancy characteristic of extra-TFIIIC sites in C. elegans embryos. Extra-TFIIIC sites constituted half of all identified TFIIIC binding sites in the genome. Extra-TFIIIC sites formed dense clusters in cis. The clusters of extra-TFIIIC sites were highly over-represented within the distal arm domains of the autosomes that presented a high level of heterochromatin-associated histone H3K9 trimethylation (H3K9me3). Furthermore, extra-TFIIIC clusters were embedded in the lamina-associated domains. Despite the heterochromatin environment of extra-TFIIIC sites, the individual clusters of extra-TFIIIC sites were devoid of and resided near the individual H3K9me3-marked regions.

CONCLUSION

Clusters of extra-TFIIIC sites were pervasive in the arm domains of C. elegans autosomes, near the outer boundaries of H3K9me3-marked regions. Given the reported activity of extra-TFIIIC sites in heterochromatin insulation in yeasts, our observation raised the possibility that TFIIIC may also demarcate heterochromatin in C. elegans.

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

Deep sampling of Hawaiian Caenorhabditis elegans reveals high genetic diversity and admixture with global populations

Hawaiian isolates of the nematode species Caenorhabditis elegans have long been known to harbor genetic diversity greater than the rest of the worldwide population, but this observation was supported by only a small number of wild strains. To better characterize the niche and genetic diversity of Hawaiian C. elegans and other Caenorhabditis species, we sampled different substrates and niches across the Hawaiian islands. We identified hundreds of new Caenorhabditis strains from known species and a new species, Caenorhabditis oiwi. Hawaiian C. elegans are found in cooler climates at high elevations but are not associated with any specific substrate, as compared to other Caenorhabditis species. Surprisingly, admixture analysis revealed evidence of shared ancestry between some Hawaiian and non-Hawaiian C. elegans strains. We suggest that the deep diversity we observed in Hawaii might represent patterns of ancestral genetic diversity in the C. elegans species before human influence.

Hybridization promotes asexual reproduction in Caenorhabditis nematodes

Although most unicellular organisms reproduce asexually, most multicellular eukaryotes are obligately sexual. This implies that there are strong barriers that prevent the origin or maintenance of asexuality arising from an obligately sexual ancestor. By studying rare asexual animal species we can gain a better understanding of the circumstances that facilitate their evolution from a sexual ancestor. Of the known asexual animal species, many originated by hybridization between two ancestral sexual species. The balance hypothesis predicts that genetic incompatibilities between the divergent genomes in hybrids can modify meiosis and facilitate asexual reproduction, but there are few instances where this has been shown. Here we report that hybridizing two sexual Caenorhabditis nematode species (C. nouraguensis females and C. becei males) alters the normal inheritance of the maternal and paternal genomes during the formation of hybrid zygotes. Most offspring of this interspecies cross die during embryogenesis, exhibiting inheritance of a diploid C. nouraguensis maternal genome and incomplete inheritance of C. becei paternal DNA. However, a small fraction of offspring develop into viable adults that can be either fertile or sterile. Fertile offspring are produced asexually by sperm-dependent parthenogenesis (also called gynogenesis or pseudogamy); these progeny inherit a diploid maternal genome but fail to inherit a paternal genome. Sterile offspring are hybrids that inherit both a diploid maternal genome and a haploid paternal genome. Whole-genome sequencing of individual viable worms shows that diploid maternal inheritance in both fertile and sterile offspring results from an altered meiosis in C. nouraguensis oocytes and the inheritance of two randomly selected homologous chromatids. We hypothesize that hybrid incompatibility between C. nouraguensis and C. becei modifies maternal and paternal genome inheritance and indirectly induces gynogenetic reproduction. This system can be used to dissect the molecular mechanisms by which hybrid incompatibilities can facilitate the emergence of asexual reproduction.