Ubiquitous Selfish Toxin-Antidote Elements in Caenorhabditis Species

Selfish conflict underlies RNA-mediated parent-of-origin effects

Genomic imprinting-the non-equivalence of maternal and paternal genomes-is a critical process that has evolved independently in many plant and mammalian species1,2. According to kinship theory, imprinting is the inevitable consequence of conflictive selective forces acting on differentially expressed parental alleles3,4. Yet, how these epigenetic differences evolve in the first place is poorly understood3,5,6. Here we report the identification and molecular dissection of a parent-of-origin effect on gene expression that might help to clarify this fundamental question. Toxin-antidote elements (TAs) are selfish elements that spread in populations by poisoning non-carrier individuals7-9. In reciprocal crosses between two Caenorhabditis tropicalis wild isolates, we found that the slow-1/grow-1 TA is specifically inactive when paternally inherited. This parent-of-origin effect stems from transcriptional repression of the slow-1 toxin by the PIWI-interacting RNA (piRNA) host defence pathway. The repression requires PIWI Argonaute and SET-32 histone methyltransferase activities and is transgenerationally inherited via small RNAs. Remarkably, when slow-1/grow-1 is maternally inherited, slow-1 repression is halted by a translation-independent role of its maternal mRNA. That is, slow-1 transcripts loaded into eggs-but not SLOW-1 protein-are necessary and sufficient to counteract piRNA-mediated repression. Our findings show that parent-of-origin effects can evolve by co-option of the piRNA pathway and hinder the spread of selfish genes that require sex for their propagation.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

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.

A toxin-antidote system contributes to interspecific reproductive isolation in rice

Breakdown of reproductive isolation facilitates flow of useful trait genes into crop plants from their wild relatives. Hybrid sterility, a major form of reproductive isolation exists between cultivated rice (Oryza sativa) and wild rice (O. meridionalis, Mer). Here, we report the cloning of qHMS1, a quantitative trait locus controlling hybrid male sterility between these two species. Like qHMS7, another locus we cloned previously, qHMS1 encodes a toxin-antidote system, but differs in the encoded proteins, their evolutionary origin, and action time point during pollen development. In plants heterozygous at qHMS1, ~ 50% of pollens carrying qHMS1-D (an allele from cultivated rice) are selectively killed. In plants heterozygous at both qHMS1 and qHMS7, ~ 75% pollens without co-presence of qHMS1-Mer and qHMS7-D are selectively killed, indicating that the antidotes function in a toxin-dependent manner. Our results indicate that different toxin-antidote systems provide stacked reproductive isolation for maintaining species identity and shed light on breakdown of hybrid male sterility.

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 dominant dpy-10 co-transformation marker using CRISPR/Cas9 and a linear repair template in Caenorhabditis tropicalis

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

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.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

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.

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.

Genetic architecture and temporal analysis of Caenorhabditis briggsae hybrid developmental delay

Identifying the alleles that reduce hybrid fitness is a major goal in the study of speciation genetics. It is rare to identify systems in which hybrid incompatibilities with minor phenotypic effects are segregating in genetically diverse populations of the same biological species. Such traits do not themselves cause reproductive isolation but might initiate the process. In the nematode Caenorhabditis briggsae, a small percent of F2 generation hybrids between two natural populations suffer from developmental delay, in which adulthood is reached after approximately 33% more time than their wild-type siblings. Prior efforts to identify the genetic basis for this hybrid incompatibility assessed linkage using one or two genetic markers on chromosome III and suggested that delay is caused by a toxin-antidote element. Here, we have genotyped F2 hybrids using multiple chromosome III markers to refine the developmental delay locus. Also, to better define the developmental delay phenotype, we measured the development rate of 66 F2 hybrids and found that delay is not restricted to a particular larval developmental stage. Deviation of the developmental delay frequency from hypothetical expectations for a toxin-antidote element adds support to the assertion that the epistatic interaction is not fully penetrant. Our mapping and refinement of the delay phenotype motivates future efforts to study the genetic architecture of hybrid dysfunction between genetically distinct populations of one species by identifying the underlying loci.

Unbalanced selection: the challenge of maintaining a social polymorphism when a supergene is selfish

Supergenes often have multiple phenotypic effects, including unexpected detrimental ones, because recombination suppression maintains associations among co-adapted alleles but also allows the accumulation of recessive deleterious mutations and selfish genetic elements. Yet, supergenes often persist over long evolutionary periods. How are such polymorphisms maintained in the face of selection, drive and drift? We present a population genetic model that investigates the conditions necessary for a stable polymorphic equilibrium when one of the supergene haplotypes is a selfish genetic element. The model fits the characteristics of the Alpine silver ant, Formica selysi, in which a large supergene underlies colony social organization, and one haplotype distorts Mendelian transmission by killing progeny that did not inherit it. The model shows that such maternal-effect killing strongly limits the maintenance of social polymorphism. Under random mating, transmission ratio distortion prevents rare single-queen colonies from invading populations of multiple-queen colonies, regardless of the fitness of each genotype. A stable polymorphic equilibrium can, however, be reached when high rates of assortative mating are combined with large fitness differences among supergene genotypes. The model reveals that the persistence of the social polymorphism is non-trivial and expected to occur only under restrictive conditions that deserve further empirical investigation. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

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.

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.

Selfing is the safest sex for Caenorhabditis tropicalis

Mating systems have profound effects on genetic diversity and compatibility. The convergent evolution of self-fertilization in three Caenorhabditis species provides a powerful lens to examine causes and consequences of mating system transitions. Among the selfers, Caenorhabditis tropicalis is the least genetically diverse and most afflicted by outbreeding depression. We generated a chromosomal-scale genome for C. tropicalis and surveyed global diversity. Population structure is very strong, and islands of extreme divergence punctuate a genomic background that is highly homogeneous around the globe. Outbreeding depression in the laboratory is caused largely by multiple Medea-like elements, genetically consistent with maternal toxin/zygotic antidote systems. Loci with Medea activity harbor novel and duplicated genes, and their activity is modified by mito-nuclear background. Segregating Medea elements dramatically reduce fitness, and simulations show that selfing limits their spread. Frequent selfing in C. tropicalis may therefore be a strategy to avoid Medea-mediated outbreeding depression.