Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Model nematodes as a practical innovation to promote high throughput screening of natural products for anthelmintics discovery in South Asia: Current challenges, proposed practical and conceptual solutions

With the increasing prevalence of anthelmintic resistance in animals recorded globally, and the threat of resistance in human helminths, the need for novel anthelmintic drugs is greater than ever. Most research aimed at discovering novel anthelmintic leads relies on high throughput screening (HTS) of large libraries of synthetic small molecules in industrial and academic settings in developed countries, even though it is the tropical countries that are most plagued by helminth infections. Tropical countries, however, have the advantage of possessing a rich flora that may yield natural products (NP) with promising anthelmintic activity. Focusing on South Asia, which produces one of the world's highest research outputs in NP and NP-based anthelmintic discovery, we find that limited basic research and funding, a lack of awareness of the utility of model organisms, poor industry-academia partnerships and lack of technological innovations greatly limit anthelmintics research in the region. Here we propose that utilizing model organisms including the free-living nematode Caenorhabditis elegans, that can potentially allow rapid target identification of novel anthelmintics, and Oscheius tipulae, a closely related, free-living nematode which is found abundantly in soil in hotter temperatures, could be a much-needed innovation that can enable cost-effective and efficient HTS of NPs for discovering compounds with anthelmintic/antiparasitic potential in South Asia and other tropical regions that historically have devoted limited funding for such research. Additionally, increased collaborations at the national, regional and international level between parasitologists and pharmacologists/ethnobotanists, setting up government-industry-academia partnerships to fund academic research, creating a centralized, regional collection of plant extracts or purified NPs as a dereplication strategy and HTS library, and holding regional C. elegans/O. tipulae-based anthelmintics workshops and conferences to share knowledge and resources regarding model organisms may collectively promote and foster a NP-based anthelmintics landscape in South Asia and beyond.

A QTL on chromosome IV explains a natural variation of QR.pap final position in Caenorhabditis elegans

In Caenorhabditis elegans , the QR neuroblast and its progeny migrate from the posterior to the anterior part of the animal during the L1 stage. We previously showed that the final position of QR.pa daughters varies among C. elegans wild isolates, with CB4932 displaying a particularly anterior QR.pap position (Dubois et al., 2021). Here, we study the genetic basis of the variation between isolates CB4932 and JU1242. We show that JU1242 alleles behave in a mostly dominant fashion. Using a Bulk Segregant Analysis, we detect a quantitative trait locus (QTL) region on chromosome IV. This QTL was confirmed using reciprocal chromosome IV introgressions.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

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.

A broad mutational target explains a fast rate of phenotypic evolution

The rapid evolution of a trait in a clade of organisms can be explained by the sustained action of natural selection or by a high mutational variance, that is the propensity to change under spontaneous mutation. The causes for a high mutational variance are still elusive. In some cases, fast evolution depends on the high mutation rate of one or few loci with short tandem repeats. Here, we report on the fastest evolving cell fate among vulva precursor cells in Caenorhabditis nematodes, that of P3.p. We identify and validate causal mutations underlying P3.p's high mutational variance. We find that these positions do not present any characteristics of a high mutation rate, are scattered across the genome and the corresponding genes belong to distinct biological pathways. Our data indicate that a broad mutational target size is the cause of the high mutational variance and of the corresponding fast phenotypic evolutionary rate.

Heritable characteristics or traits of a group of organisms, for example the large brain size of primates or the hooves of a horse, are determined by genes, the environment, and by the interactions between them. Traits can change over time and generations when enough mutations in these genes have spread in a species to result in visible differences.

However, some traits, such as the large brain of primates, evolve faster than others, but why this is the case has been unclear. It could be that a few specific genes important for that trait in question mutate at a high rate, or, that many genes affect the trait, creating a lot of variation for natural selection to choose from.

Here, Besnard, Picao-Osorio et al. studied the roundworm Caenorhabditis elegans to better understand the causes underlying the different rates of trait evolution. These worms have a short life cycle and evolve quickly over many generations, making them an ideal candidate for studying mutation rates in different traits.

Previous studies have shown that one of C. elegans’ six cells of the reproductive system evolves faster than the others. To investigate this further, Besnard, Picao-Osorio et al. analysed the genetic mutations driving change in this cell in 250 worm generations. The results showed that five mutations in five different genes – all responsible for different processes in the cells – were behind the supercharged evolution of this particular cell. This suggests that fast evolution results from natural selection acting upon a collection of genes, rather than one gene, and that many genes and pathways shape this trait.

In conclusion, these results demonstrate that how traits are coded at the molecular level, in one gene or many, can influence the rate at which they evolve.

Interspecific Variation in Nematode Responses to Metals

Performing toxicity testing on multiple species with differing degrees of evolutionary relatedness can provide important information on how chemical sensitivity varies among species and can help pinpoint the biological drivers of species sensitivity. Such knowledge could ultimately be used to design better multispecies predictive ecological risk assessment models and identify particularly sensitive species. However, laboratory toxicity tests involving multiple species can also be resource intensive, especially when each species has unique husbandry conditions. We performed lethality tests with 2 metals, copper chloride and zinc chloride, on 5 different nematode species, which are nested in their degree of evolutionary relatedness: Caenorhabditis briggsae, Caenorhabditis elegans, Oscheius myriophila, Oscheius tipulae, and Pristionchus pacificus. All species were successfully cultured and tested concurrently with limited resources, demonstrating that inexpensive, multispecies nematode toxicity testing systems are achievable. The results indicate that P. pacificus is the most sensitive to both metals. Conversely, C. elegans is the least sensitive species to copper, but the second most sensitive to zinc, indicating that species relationships do not necessarily predict species sensitivity. Toxicity testing with additional nematode species and types of chemicals is feasible and will help form more generalizable conclusions about relative species sensitivity. Environ Toxicol Chem 2020;39:1006-1016. © 2020 SETAC.

Chromosome-Wide Evolution and Sex Determination in the Three-Sexed Nematode Auanema rhodensis

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis, males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM (doublesex/mab-3) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis, instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

Comparative Epigenomics Reveals that RNA Polymerase II Pausing and Chromatin Domain Organization Control Nematode piRNA Biogenesis

Piwi-interacting RNAs (piRNAs) are important for genome regulation across metazoans, but their biogenesis evolves rapidly. In Caenorhabditis elegans, piRNA loci are clustered within two 3-Mb regions on chromosome IV. Each piRNA locus possesses an upstream motif that recruits RNA polymerase II to produce an ∼28 nt primary transcript. We used comparative epigenomics across nematodes to gain insight into the origin, evolution, and mechanism of nematode piRNA biogenesis. We show that the piRNA upstream motif is derived from core promoter elements controlling snRNA transcription. We describe two alternative modes of piRNA organization in nematodes: in C. elegans and closely related nematodes, piRNAs are clustered within repressive H3K27me3 chromatin, while in other species, typified by Pristionchus pacificus, piRNAs are found within introns of active genes. Additionally, we discover that piRNA production depends on sequence signals associated with RNA polymerase II pausing. We show that pausing signals synergize with chromatin to control piRNA transcription.

•

Nematode piRNA transcription evolved from small nuclear RNA biogenesis

•

Clustered piRNAs are produced from regulated (H3K27me3) chromatin domains

•

Dispersed piRNAs are produced from active (H3K36me3) chromatin domains

•

RNA polymerase II pausing determines the short (∼28 nt) length of piRNA precursors

Necessity and Contingency in Developmental Genetic Screens: EGF, Wnt, and Semaphorin Pathways in Vulval Induction of the Nematode Oscheius tipulae

Genetic screens in the nematode Caenorhabditis elegans have identified EGF and Notch pathways as key for vulval precursor cell fate patterning. Here, Vargas-Velazquez, Besnard, and Félix report on the molecular identification of...

Genetic screens in the nematode Caenorhabditis elegans identified the EGF/Ras and Notch pathways as central for vulval precursor cell fate patterning. Schematically, the anchor cell secretes EGF, inducing the P6.p cell to a primary (1°) vulval fate; P6.p in turn induces its neighbors to a secondary (2°) fate through Delta-Notch signaling and represses Ras signaling. In the nematode Oscheius tipulae, the anchor cell successively induces 2° then 1° vulval fates. Here, we report on the molecular identification of mutations affecting vulval induction in O. tipulae. A single Induction Vulvaless mutation was found, which we identify as a cis-regulatory deletion in a tissue-specific enhancer of the O. tipulae lin-3 homolog, confirmed by clustered regularly interspaced short palindromic repeats/Cas9 mutation. In contrast to this predictable Vulvaless mutation, mutations resulting in an excess of 2° fates unexpectedly correspond to the plexin/semaphorin pathway. Hyperinduction of P4.p and P8.p in these mutants likely results from mispositioning of these cells due to a lack of contact inhibition. The third signaling pathway found by forward genetics in O. tipulae is the Wnt pathway; a decrease in Wnt pathway activity results in loss of vulval precursor competence and induction, and 1° fate miscentering on P5.p. Our results suggest that the EGF and Wnt pathways have qualitatively similar activities in vulval induction in C. elegans and O. tipulae, albeit with quantitative differences in the effects of mutation. Thus, the derived induction process in C. elegans with an early induction of the 1° fate appeared during evolution, after the recruitment of the EGF pathway for vulval induction.

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Natural Genetic Variation in a Multigenerational Phenotype in C. elegans

Although heredity mostly relies on the transmission of DNA sequence, additional molecular and cellular features are heritable across several generations. In the nematode Caenorhabditis elegans, insights into such unconventional inheritance result from two lines of work. First, the mortal germline (Mrt) phenotype was defined as a multigenerational phenotype whereby a selfing lineage becomes sterile after several generations, implying multigenerational memory [1, 2]. Second, certain RNAi effects are heritable over several generations in the absence of the initial trigger [3-5]. Both lines of work converged when the subset of Mrt mutants that are heat sensitive were found to closely correspond to mutants defective in the RNAi-inheritance machinery, including histone modifiers [6-9]. Here, we report the surprising finding that several C. elegans wild isolates display a heat-sensitive mortal germline phenotype in laboratory conditions: upon chronic exposure to higher temperatures, such as 25°C, lines reproducibly become sterile after several generations. This phenomenon is reversible, as it can be suppressed by temperature alternations at each generation, suggesting a non-genetic basis for the sterility. We tested whether natural variation in the temperature-induced Mrt phenotype was of genetic nature by building recombinant inbred lines between the isolates MY10 (Mrt) and JU1395 (non-Mrt). Using bulk segregant analysis, we detected two quantitative trait loci. After further recombinant mapping and genome editing, we identified the major causal locus as a polymorphism in the set-24 gene, encoding a SET- and SPK-domain protein. We conclude that C. elegans natural populations may harbor natural genetic variation in epigenetic inheritance phenomena.