Selfing is the safest sex for Caenorhabditis tropicalis

An animal toxin-antidote system kills cells by creating a novel cation channel

Toxin-antidote systems are selfish genetic elements composed of a linked toxin and antidote. The peel-1 zeel-1 toxin-antidote system in C. elegans consists of a transmembrane toxin protein PEEL-1 which acts cell autonomously to kill cells. Here we investigate the molecular mechanism of PEEL-1 toxicity. We find that PEEL-1 requires a small membrane protein, PMPL-1, for toxicity. Together, PEEL-1 and PMPL-1 are sufficient for toxicity in a heterologous system, HEK293T cells, and cause cell swelling and increased cell permeability to monovalent cations. Using purified proteins, we show that PEEL-1 and PMPL-1 allow ion flux through lipid bilayers and generate currents which resemble ion channel gating. Our work suggests that PEEL-1 kills cells by co-opting PMPL-1 and creating a cation channel.

The importance of using local species in ecotoxicological studies: nematodes of Amazonian occurrence vs. Caenorhabditis elegans in the analysis of lethal and sublethal effects of aluminium

It is recognized that in bioassays, especially those conducted for ecotoxicological purposes, preference should be given to the use of species that are adapted to the physical-chemical conditions of the environment to be monitored. However, to establish the use of alternative species instead of the standardized ones, it is recommended to carry out tests to assess/compare their sensitivity to contaminants. This study assessed the lethal and sublethal effects (growth, fertility, and reproduction) of different aluminium concentrations, including environmentally relevant concentrations recorded in the Amazon, on two nematode species (C. tropicalis and C. briggsae) with Amazonian occurrence and C. elegans. The species' responses to aluminium exposure were different. In tests to assess lethal effect, C. elegans was the most sensitive (LC50 = 3.32 ± 1.89 mg/L), while C. tropicalis was the least sensitive (LC50 = 6.98 ± 2.20 mg/L). The LC50 for C. briggsae could not be estimated due to the lack of a concentration-dependent response. On the other hand, when sublethal effects were assessed at low aluminium concentrations (environmentally relevant concentrations), C. tropicalis was the most sensitive with an inhibition rate in both reproduction and growth; C. elegans was the least sensitive, and C. briggsae showed an intermediate response. Therefore, C. tropicalis and C. elegans adopted opposite strategies in response to aluminium exposure. This study reinforces the use of local species in ecotoxicological tests and suggests the use of C. tropicalis as a test organism in future bioassays to evaluate the effects of contaminants, particularly in the tropical/Amazon region.

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

Independent mechanisms of benzimidazole resistance across Caenorhabditis nematodes

Benzimidazoles (BZs), a widely used class of anthelmintic drugs, target beta-tubulin proteins, disrupt microtubule formation, and cause nematode death. In parasitic nematode species, mutations in beta-tubulin genes (e.g., isotype-1 beta-tubulin) are predicted to inhibit BZ binding and are associated with BZ resistance. Similarly, in the free-living nematode Caenorhabditis elegans, mutations in an isotype-1 beta-tubulin ortholog, ben-1, are the primary drivers of BZ resistance. The recurrent association of BZ resistance with beta-tubulins suggests that BZ resistance is repeatedly caused by mutations in beta-tubulin genes, an example of repeated evolution of drug resistance across nematode species. To evaluate the hypothesis of repeated evolution of BZ resistance mediated by beta-tubulin, we identified predicted resistance alleles in beta-tubulin genes across wild strains from three Caenorhabditis species: C. elegans, Caenorhabditis briggsae, and Caenorhabditis tropicalis. We hypothesized that, if these species experienced similar selective pressures, they would evolve resistance to BZs by mutations in any of three beta-tubulin genes (ben-1, tbb-1, and tbb-2). Using high-throughput development assays, we tested the association of predicted beta-tubulin alleles with BZ resistance. We found that a heterogeneous set of variants identified in C. elegans ben-1 were associated with BZ resistance. In C. briggsae, only two variants in ben-1, predicted to encode a premature stop codon (W21stop) and a missense substitution (Q134H), were associated with BZ resistance. In C. tropicalis, two missense variants were identified in ben-1, but neither was associated with BZ resistance. C. briggsae and C. tropicalis might have evolved BZ resistance by mutations in other beta-tubulin genes, but we found that variants in tbb-1 or tbb-2 in these species were not associated with BZ resistance. Our findings reveal a lack of repeated evolution of BZ resistance across the three Caenorhabditis species and highlight the importance of defining BZ resistance mechanisms outside of beta-tubulins.

Pervasive Conservation of Intron Number and Other Genetic Elements Revealed by a Chromosome-level Genome Assembly of the Hyper-polymorphic Nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gamut of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and, notably, across Metazoa. Here, we present a high-quality, gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat-rich arms. A comparison of C. brenneri with other nematodes from the "Elegans" group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation of orthogroup size, indicative of high rates of gene turnover, consistent with previous studies. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. A comparison of gene structures revealed a strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

Parental-effect gene-drive elements under partial selfing, or why do Caenorhabditis genomes have hyperdivergent regions?

Self-fertile Caenorhabditis nematodes carry a surprising number of Medea elements, alleles that act in heterozygous mothers and cause death or developmental delay in offspring that do not inherit them. At some loci, both alleles in a cross operate as independent Medeas, affecting all the homozygous progeny of a selfing heterozygote. The genomic coincidence of Medea elements and ancient, deeply coalescing haplotypes, which pepper the otherwise homogeneous genomes of these animals, raises questions about how these apparent gene-drive elements persist for long periods of time. Here, I investigate how mating system affects the evolution of Medeas, and their paternal-effect counterparts, peels. Despite an intuition that antagonistic alleles should induce balancing selection by killing homozygotes, models show that, under partial selfing, antagonistic elements experience positive frequency dependence: the common allele drives the rare one extinct, even if the rare one is more penetrant. Analytical results for the threshold frequency required for one allele to invade a population show that a very weakly penetrant allele, one whose effects would escape laboratory detection, could nevertheless prevent a much more penetrant allele from invading under high rates of selfing. Ubiquitous weak antagonistic Medeas and peels could then act as localized barriers to gene flow between populations, generating genomic islands of deep coalescence. Analysis of gene expression data, however, suggests that this cannot be the whole story. A complementary explanation is that ordinary ecological balancing selection generates ancient haplotypes on which Medeas can evolve, while high homozygosity in these selfers minimizes the role of gene drive in their evolution.

Selection Can Favor a Recombination Landscape That Limits Polygenic Adaptation

Modifiers of recombination rates have been described but the selective pressures acting on them and their effect on adaptation to novel environments remain unclear. We performed experimental evolution in the nematode Caenorhabditis elegans using alternative rec-1 alleles modifying the position of meiotic crossovers along chromosomes without detectable direct fitness effects. We show that adaptation to a novel environment is impaired by the allele that decreases recombination rates in the genomic regions containing fitness variation. However, the allele that impairs adaptation is indirectly favored by selection, because it increases recombination rates and reduces the associations among beneficial and deleterious variation located in its chromosomal vicinity. These results validate theoretical expectations about the evolution of recombination but suggest that genome-wide polygenic adaptation is of little consequence to indirect selection on recombination rate modifiers.

Divergent C. elegans toxin alleles are suppressed by distinct mechanisms

Toxin-antidote elements (TAs) are selfish DNA sequences that bias their transmission to the next generation. TAs typically consist of two linked genes: a toxin and an antidote. The toxin kills progeny that do not inherit the TA, while the antidote counteracts the toxin in progeny that inherit the TA. We previously discovered two TAs in C. elegans that follow the canonical TA model of two linked genes: peel-1/zeel-1 and sup-35/pha-1. Here, we report a new TA that exists in three distinct states across the C. elegans population. The canonical TA, which is found in isolates from the Hawaiian islands, consists of two genes that encode a maternally deposited toxin (MLL-1) and a zygotically expressed antidote (SMLL-1). The toxin induces larval lethality in embryos that do not inherit the antidote gene. A second version of the TA has lost the toxin gene but retains a partially functional antidote. Most C. elegans isolates, including the standard laboratory strain N2, carry a highly divergent allele of the toxin that has retained its activity, but have lost the antidote through pseudogenization. We show that the N2 toxin allele has acquired mutations that enable piRNA binding to initiate MUT-16-dependent 22G small RNA amplification that targets the transcript for degradation. The N2 haplotype represents the first naturally occurring unlinked toxin-antidote system where the toxin is post-transcriptionally suppressed by endogenous small RNA pathways.

Life history in Caenorhabditis elegans: from molecular genetics to evolutionary ecology

Life history is defined by traits that reflect key components of fitness, especially those relating to reproduction and survival. Research in life history seeks to unravel the relationships among these traits and understand how life history strategies evolve to maximize fitness. As such, life history research integrates the study of the genetic and developmental mechanisms underlying trait determination with the evolutionary and ecological context of Darwinian fitness. As a leading model organism for molecular and developmental genetics, Caenorhabditis elegans is unmatched in the characterization of life history-related processes, including developmental timing and plasticity, reproductive behaviors, sex determination, stress tolerance, and aging. Building on recent studies of natural populations and ecology, the combination of C. elegans' historical research strengths with new insights into trait variation now positions it as a uniquely valuable model for life history research. In this review, we summarize the contributions of C. elegans and related species to life history and its evolution. We begin by reviewing the key characteristics of C. elegans life history, with an emphasis on its distinctive reproductive strategies and notable life cycle plasticity. Next, we explore intraspecific variation in life history traits and its underlying genetic architecture. Finally, we provide an overview of how C. elegans has guided research on major life history transitions both within the genus Caenorhabditis and across the broader phylum Nematoda. While C. elegans is relatively new to life history research, significant progress has been made by leveraging its distinctive biological traits, establishing it as a highly cross-disciplinary system for life history studies.

Parental-effect gene-drive elements under partial selfing, or why do Caenorhabditis genomes have hyperdivergent regions?

Self-fertile Caenorhabditis nematodes carry a surprising number of Medea elements, alleles that act in heterozygous mothers and cause death or developmental delay in offspring that don't inherit them. At some loci, both alleles in a cross operate as independent Medeas, affecting all the homozygous progeny of a selfing heterozygote. The genomic coincidence of Medea elements and ancient, deeply coalescing haplotypes, which pepper the otherwise homogeneous genomes of these animals, raises questions about how these apparent gene-drive elements persist for long periods of time. Here I investigate how mating system affects the evolution of Medeas, and their paternal-effect counterparts, peels. Despite an intuition that antagonistic alleles should induce balancing selection by killing homozygotes, models show that, under partial selfing, antagonistic elements experience positive frequency dependence: the common allele drives the rare one extinct, even if the rare one is more penetrant. Analytical results for the threshold frequency required for one allele to invade a population show that a very weakly penetrant allele, one whose effects would escape laboratory detection, could nevertheless prevent a much more penetrant allele from invading under high rates of selfing. Ubiquitous weak antagonistic Medeas and peels could then act as localized barriers to gene flow between populations, generating genomic islands of deep coalescence. Analysis of gene expression data, however, suggest that this cannot be the whole story. A complementary explanation is that ordinary ecological balancing selection generates ancient haplotypes on which Medeas can evolve, while high homozygosity in these selfers minimizes the role of gene drive in their evolution.

Selfing Promotes Spread and Introgression of Segregation Distorters in Hermaphroditic Plants

Segregation distorters (SDs) are genetic elements that distort the Mendelian segregation ratio to favor their own transmission and are able to spread even when they incur fitness costs on organisms carrying them. Depending on the biology of the host organisms and the genetic architecture of the SDs, the population dynamics of SDs can be highly variable. Inbreeding is considered an effective mechanism for inhibiting the spread of SDs in populations, and can evolve as a defense mechanism against SDs in some systems. However, we show that inbreeding in the form of selfing in fact promotes the spread of SDs acting as pollen killers in a toxin-antidote system in hermaphroditic plants by two mechanisms: (i) By reducing the effective recombination rate between killer and antidote loci in the two-locus system and (ii) by increasing the proportion of SD alleles in individual flowers, rather than in the general gene-pool. We also show that in rice (Oryza sativa L.), a typical hermaphroditic plant, all molecularly characterized SDs associated with pollen killing were involved in population hybridization and have introgressed across different species. Paradoxically, these loci, which are associated with hybrid incompatibility and can be thought of as Bateson-Dobzhansky-Muller incompatibility loci are expected to reduce gene-flow between species, in fact cross species boundaries more frequently than random loci, and may act as important drivers of introgression.

Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

Mechanism of an animal toxin-antidote system

Toxin-antidote systems are selfish genetic elements composed of a linked toxin and antidote. The peel-1 zeel-1 toxin-antidote system in C. elegans consists of a transmembrane toxin protein PEEL-1 which acts cell autonomously to kill cells. Here we investigate the molecular mechanism of PEEL-1 toxicity. We find that PEEL-1 requires a small membrane protein, PMPL-1, for toxicity. Together, PEEL-1 and PMPL-1 are sufficient for toxicity in a heterologous system, HEK293T cells, and cause cell swelling and increased cell permeability to monovalent cations. Using purified proteins, we show that PEEL-1 and PMPL-1 allow ion flux through lipid bilayers and generate currents which resemble ion channel gating. Our work suggests that PEEL-1 kills cells by co-opting PMPL-1 and creating a cation channel.

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.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

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.

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.

Meiotic drive, postzygotic isolation, and the Snowball Effect

As populations diverge, they accumulate incompatibilities which reduce gene flow and facilitate the formation of new species. Simple models suggest that the genes that cause Dobzhansky-Muller incompatibilities should accumulate at least as fast as the square of the number of substitutions between taxa, the so-called snowball effect. We show, however, that in the special- but possibly common- case in which hybrid sterility is due primarily to cryptic meiotic (gametic) drive, the number of genes that cause postzygotic isolation may increase nearly linearly with the number of substitutions between species.

A toxin-antidote selfish element increases fitness of its host

Selfish genetic elements can promote their transmission at the expense of individual survival, creating conflict between the element and the rest of the genome. Recently, a large number of toxin-antidote (TA) post-segregation distorters have been identified in non-obligate outcrossing nematodes. Their origin and the evolutionary forces that keep them at intermediate population frequencies are poorly understood. Here, we study a TA element in Caenorhabditis elegans called zeel-1;peel-1. Two major haplotypes of this locus, with and without the selfish element, segregate in C. elegans. We evaluate the fitness consequences of the zeel-1;peel-1 element outside of its role in gene drive in non-outcrossing animals and demonstrate that loss of the toxin peel-1 decreased fitness of hermaphrodites and resulted in reductions in fecundity and body size. These findings suggest a biological role for peel-1 beyond toxin lethality. This work demonstrates that a TA element can provide a fitness benefit to its hosts either during their initial evolution or by being co-opted by the animals following their selfish spread. These findings guide our understanding on how TA elements can remain in a population where gene drive is minimized, helping resolve the mystery of prevalent TA elements in selfing animals.

The parasitic nematode Strongyloides ratti exists predominantly as populations of long-lived asexual lineages

Nematodes are important parasites of people and animals, and in natural ecosystems they are a major ecological force. Strongyloides ratti is a common parasitic nematode of wild rats and we have investigated its population genetics using single-worm, whole-genome sequencing. We find that S. ratti populations in the UK consist of mixtures of mainly asexual lineages that are widely dispersed across a host population. These parasite lineages are likely very old and may have originated in Asia from where rats originated. Genes that underly the parasitic phase of the parasite's life cycle are hyperdiverse compared with the rest of the genome, and this may allow the parasites to maximise their fitness in a diverse host population. These patterns of parasitic nematode population genetics have not been found before and may also apply to Strongyloides spp. that infect people, which will affect how we should approach their control.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

BACKGROUND

The nematode Caenorhabditis briggsae has been used as a model in comparative genomics studies with Caenorhabditis elegans because of their striking morphological and behavioral similarities. However, the potential of C. briggsae for comparative studies is limited by the quality of its genome resources. The genome resources for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans. The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of software-derived gene predictions that contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 gene models and underlying transcriptomic data to repair software-derived errors.

RESULTS

We designed a detailed workflow to train a team of nine students to manually curate gene models using RNA read alignments. We manually inspected the gene models, proposed corrections to the coding sequences of over 8,000 genes, and modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality and showed that manual curation led to substantial improvements in the protein sequence length accuracy of QX1410 genes. Additionally, collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

CONCLUSIONS

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. Our manual curation efforts have brought the QX1410 gene models to a comparable level of quality as the extensively curated AF16 gene models. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

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.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

Background

The nematode Caenorhabditis briggsae has been used as a model for genomics studies compared to Caenorhabditis elegans because of its striking morphological and behavioral similarities. These studies yielded numerous findings that have expanded our understanding of nematode development and evolution. However, the potential of C. briggsae to study nematode biology is limited by the quality of its genome resources. The reference genome and gene models for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans . The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of protein-coding gene predictions generated from short- and long-read transcriptomic data. Because of the limitations of gene prediction software, the existing gene models for QX1410 contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 software-derived gene models and underlying transcriptomic data to improve the protein-coding gene models of the C. briggsae QX1410 genome.

Results

We designed a detailed workflow to train a team of nine students to manually curate genes using RNA read alignments and predicted gene models. We manually inspected the gene models using the genome annotation editor, Apollo, and proposed corrections to the coding sequences of over 8,000 genes. Additionally, we modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality before and after curation. Manual curation led to a substantial improvement in the protein sequence length accuracy of QX1410 genes. We also compared the curated QX1410 gene models against the existing AF16 gene models. The manual curation efforts yielded QX1410 gene models that are similar in quality to the extensively curated AF16 gene models in terms of protein-length accuracy and biological completeness scores. Collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

Conclusions

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. Comparative genomic analysis using a related species with high-quality reference genome(s) and gene models can be used to quantify improvements in gene model quality in a newly sequenced genome. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. The chromosome-level reference genome for the C. briggsae strain QX1410 far surpasses the quality of the genome of the laboratory strain AF16, and our manual curation efforts have brought the QX1410 gene models to a comparable level of quality to the previous reference, AF16. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

Higher-order epistasis shapes natural variation in germ stem cell niche activity

To study how natural allelic variation explains quantitative developmental system variation, we characterized natural differences in germ stem cell niche activity, measured as progenitor zone (PZ) size, between two Caenorhabditis elegans isolates. Linkage mapping yielded candidate loci on chromosomes II and V, and we found that the isolate with a smaller PZ size harbours a 148 bp promoter deletion in the Notch ligand, lag-2/Delta, a central signal promoting germ stem cell fate. As predicted, introducing this deletion into the isolate with a large PZ resulted in a smaller PZ size. Unexpectedly, restoring the deleted ancestral sequence in the isolate with a smaller PZ did not increase-but instead further reduced-PZ size. These seemingly contradictory phenotypic effects are explained by epistatic interactions between the lag-2/Delta promoter, the chromosome II locus, and additional background loci. These results provide first insights into the quantitative genetic architecture regulating an animal stem cell system.

Transcriptomic profiling of sex-specific olfactory neurons reveals subset-specific receptor expression in Caenorhabditis elegans

The nematode Caenorhabditis elegans utilizes chemosensation to navigate an ever-changing environment for its survival. A class of secreted small-molecule pheromones, termed ascarosides, play an important role in olfactory perception by affecting biological functions ranging from development to behavior. The ascaroside #8 (ascr#8) mediates sex-specific behaviors, driving avoidance in hermaphrodites and attraction in males. Males sense ascr#8 via the ciliated male-specific cephalic sensory (CEM) neurons, which exhibit radial symmetry along dorsal-ventral and left-right axes. Calcium imaging studies suggest a complex neural coding mechanism that translates stochastic physiological responses in these neurons to reliable behavioral outputs. To test the hypothesis that neurophysiological complexity arises from differential expression of genes, we performed cell-specific transcriptomic profiling; this revealed between 18 and 62 genes with at least twofold higher expression in a specific CEM neuron subtype vs both other CEM neurons and adult males. These included two G protein-coupled receptor (GPCR) genes, srw-97 and dmsr-12, that were specifically expressed in nonoverlapping subsets of CEM neurons and whose expression was confirmed by GFP reporter analysis. Single CRISPR-Cas9 knockouts of either srw-97 or dmsr-12 resulted in partial defects, while a double knockout of both srw-97 and dmsr-12 completely abolished the attractive response to ascr#8. Together, our results suggest that the evolutionarily distinct GPCRs SRW-97 and DMSR-12 act nonredundantly in discrete olfactory neurons to facilitate male-specific sensation of ascr#8.

Widespread sex ratio polymorphism in Caenorhabditis nematodes

Although equal sex ratio is ubiquitous and represents an equilibrium in evolutionary theory, biased sex ratios are predicted for certain local conditions. Cases of sex ratio bias have been mostly reported for single species, but little is known about its evolution above the species level. Here, we surveyed progeny sex ratios in 23 species of the nematode genus Caenorhabditis, including 19 for which we tested multiple strains. For the species with multiple strains, five species had female-biased and two had non-biased sex ratios in all strains, respectively. The other 12 species showed polymorphic sex ratios across strains. Female-biased sex ratios could be due to sperm competition whereby X-bearing sperm outcompete nullo-X sperm during fertilization. In this model, when sperm are limited allowing all sperm to be used, sex ratios are expected to be equal. However, in assays limiting mating to a few hours, most strains showed similarly biased sex ratios compared with unlimited mating experiments, except that one C. becei strain showed significantly reduced female bias compared with unlimited mating. Our study shows frequent polymorphism in sex ratios within Caenorhabditis species and that sperm competition alone cannot explain the sex ratio bias.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

The compact genome of Caenorhabditis niphades n. sp., isolated from a wood-boring weevil, Niphades variegatus

BACKGROUND

The first metazoan genome sequenced, that of Caenorhabditis elegans, has motivated animal genome evolution studies. To date > 50 species from the genus Caenorhabditis have been sequenced, allowing research on genome variation.

RESULTS

In the present study, we describe a new gonochoristic species, Caenorhabditis niphades n. sp., previously referred as C. sp. 36, isolated from adult weevils (Niphades variegatus), with whom they appear to be tightly associated during its life cycle. Along with a species description, we sequenced the genome of C. niphades n. sp. and produced a chromosome-level assembly. A genome comparison highlighted that C. niphades n. sp. has the smallest genome (59 Mbp) so far sequenced in the Elegans supergroup, despite being closely related to a species with an exceptionally large genome, C. japonica.

CONCLUSIONS

The compact genome of C. niphades n. sp. can serve as a key resource for comparative evolutionary studies of genome and gene number expansions in Caenorhabditis species.

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.

C. elegans as a model for inter-individual variation in metabolism

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors1-4. It is difficult to connect differences in metabolism to genomic variation and derive underlying molecular mechanisms in humans, owing to differences in diet and lifestyle, among others. Here we use the nematode Caenorhabditis elegans as a model to study inter-individual variation in metabolism. By comparing three wild strains and the commonly used N2 laboratory strain, we find differences in the abundances of both known metabolites and those that have not to our knowledge been previously described. The latter metabolites include conjugates between 3-hydroxypropionate (3HP) and several amino acids (3HP-AAs), which are much higher in abundance in one of the wild strains. 3HP is an intermediate in the propionate shunt pathway, which is activated when flux through the canonical, vitamin-B12-dependent propionate breakdown pathway is perturbed5. We show that increased accumulation of 3HP-AAs is caused by genetic variation in HPHD-1, for which 3HP is a substrate. Our results suggest that the production of 3HP-AAs represents a 'shunt-within-a-shunt' pathway to accommodate a reduction-of-function allele in hphd-1. This study provides a step towards the development of metabolic network models that capture individual-specific differences of metabolism and more closely represent the diversity that is found in entire species.

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.

Gene drives and population persistence vs elimination: The impact of spatial structure and inbreeding at low density

Synthetic gene drive constructs are being developed to control disease vectors, invasive species, and other pest species. In a well-mixed random mating population a sufficiently strong gene drive is expected to eliminate a target population, but it is not clear whether the same is true when spatial processes play a role. In species with an appropriate biology it is possible that drive-induced reductions in density might lead to increased inbreeding, reducing the efficacy of drive, eventually leading to suppression rather than elimination, regardless of how strong the drive is. To investigate this question we analyse a series of explicitly solvable stochastic models considering a range of scenarios for the relative timing of mating, reproduction, and dispersal and analyse the impact of two different types of gene drive, a Driving Y chromosome and a homing construct targeting an essential gene. We find in all cases a sufficiently strong Driving Y will go to fixation and the population will be eliminated, except in the one life history scenario (reproduction and mating in patches followed by dispersal) where low density leads to increased inbreeding, in which case the population persists indefinitely, tending to either a stable equilibrium or a limit cycle. These dynamics arise because Driving Y males have reduced mating success, particularly at low densities, due to having fewer sisters to mate with. Increased inbreeding at low densities can also prevent a homing construct from eliminating a population. For both types of drive, if there is strong inbreeding depression, then the population cannot be rescued by inbreeding and it is eliminated. These results highlight the potentially critical role that low-density-induced inbreeding and inbreeding depression (and, by extension, other sources of Allee effects) can have on the eventual impact of a gene drive on a target population.

Genetic diversity estimates for the Caenorhabditis Intervention Testing Program screening panel

The Caenorhabditis Intervention Testing Program (CITP) was founded on the principle that compounds with positive effects across a genetically diverse test-set should have an increased probability of engaging conserved biochemical pathways with mammalian translational potential. To fulfill its mandate, the CITP uses a genetic diversity panel of Caenorhabditis strains for assaying longevity effects of candidate compounds. The panel comprises 22 strains from three different species, collected globally, to achieve inter-population genetic diversity. The three represented species, C. elegans, C. briggsae, and C. tropicalis, are all sequential hermaphrodites, which simplifies experimental procedures while maximizing intra-population homogeneity. Here, we present estimates of the genetic diversity encapsulated by the constituent strains in the panel based on their most recently published and publicly available whole-genome sequences, as well as two newly generated genomic data sets. We observed average genome-wide nucleotide diversity (π) within the C. elegans (1.2e-3), C. briggsae (7.5e-3), and C. tropicalis strains (2.6e-3) greater than estimates for human populations, and comparable to that found in mouse populations. Our analysis supports the assumption that the CITP screening panel encompasses broad genetic diversity, suggesting that lifespan-extending chemicals with efficacy across the panel should be enriched for interventions that function on conserved processes that are shared across genetic backgrounds. While the diversity panel was established by the CITP for studying longevity interventions, the panel may prove useful for the broader research community when seeking broadly efficacious interventions for any phenotype with potential genetic background effects.

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.

Genetic exchange with an outcrossing sister species causes severe genome-wide dysregulation in a selfing Caenorhabditis nematode

Different modes of reproduction evolve rapidly, with important consequences for genome composition. Selfing species often occupy a similar niche as their outcrossing sister species with which they are able to mate and produce viable hybrid progeny, raising the question of how they maintain genomic identity. Here, we investigate this issue by using the nematode Caenorhabditis briggsae, which reproduces as a hermaphrodite, and its outcrossing sister species Caenorhabditis nigoni We hypothesize that selfing species might develop some barriers to prevent gene intrusions through gene regulation. We therefore examined gene regulation in the hybrid F2 embryos resulting from reciprocal backcrosses between F1 hybrid progeny and C. nigoni or C. briggsae F2 hybrid embryos with ∼75% of their genome derived from C. briggsae (termed as bB2) were inviable, whereas those with ∼75% of their genome derived from C. nigoni (termed as nB2) were viable. Misregulation of transposable elements, coding genes, and small regulatory RNAs was more widespread in the bB2 compared with the nB2 hybrids, which is a plausible explanation for the differential phenotypes between the two hybrids. Our results show that regulation of the C. briggsae genome is strongly affected by genetic exchanges with its outcrossing sister species, C. nigoni, whereas regulation of the C. nigoni genome is more robust on genetic exchange with C. briggsae The results provide new insights into how selfing species might maintain their identity despite genetic exchanges with closely related outcrossing species.

Evolution: Drivers in a Traffic Jam

Genetic loci coding for a toxin and its antidote behave like selfish elements. Two new studies find an accumulation of such elements in one species.