Single worm transcriptomics identifies a developmental core network of oscillating genes with deep conservation across nematodes

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Characterization of the Pristionchus pacificus "epigenetic toolkit" reveals the evolutionary loss of the histone methyltransferase complex PRC2

Comparative approaches have revealed both divergent and convergent paths to achieving shared developmental outcomes. Thus, only through assembling multiple case studies can we understand biological principles. Yet, despite appreciating the conservation-or lack thereof-of developmental networks, the conservation of epigenetic mechanisms regulating these networks is poorly understood. The nematode Pristionchus pacificus has emerged as a model system of plasticity and epigenetic regulation as it exhibits a bacterivorous or omnivorous morph depending on its environment. Here, we determined the "epigenetic toolkit" available to P. pacificus as a resource for future functional work on plasticity, and as a comparison with Caenorhabditis elegans to investigate the conservation of epigenetic mechanisms. Broadly, we observed a similar cast of genes with putative epigenetic function between C. elegans and P. pacificus. However, we also found striking differences. Most notably, the histone methyltransferase complex PRC2 appears to be missing in P. pacificus. We described the deletion/pseudogenization of the PRC2 genes mes-2 and mes-6 and concluded that both were lost in the last common ancestor of P. pacificus and a related species P. arcanus. Interestingly, we observed the enzymatic product of PRC2 (H3K27me3) by mass spectrometry and immunofluorescence, suggesting that a currently unknown methyltransferase has been co-opted for heterochromatin silencing. Altogether, we have provided an inventory of epigenetic genes in P. pacificus to compare with C. elegans. This inventory will enable reverse-genetic experiments related to plasticity and has revealed the first loss of PRC2 in a multicellular organism.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data

SUMMARY

For model species, single-cell RNA-based cell atlases are available. A good cell atlas includes all major stages in a species' ontogeny, and soon, they will be standard even for nonmodel species. Here, we propose a Python package called oggmap, which allows for the easy extraction of an orthomap (gene ages per orthogroup) for any given query species from OrthoFinder and other gene family data resources, like homologous groups from eggNOG or PLAZA. oggmap provides extracted gene ages for more than thousand eukaryotic species which can be further used to calculate gene age-weighted expression data from scRNA sequencing objects using the Python Scanpy toolkit. Not limited to one transcriptome evolutionary index, oggmap can visualize the individual gene category (e.g. age class, nucleotide diversity bin) and their corresponding expression profiles to investigate scRNA-based cell type assignments in an evolutionary context.

AVAILABILITY AND IMPLEMENTATION

oggmap source code is available at https://github.com/kullrich/oggmap, documentation is available at https://oggmap.readthedocs.io/en/latest/. oggmap can be installed via PyPi or directly used via a docker container.

Thousands of Pristionchus pacificus orphan genes were integrated into developmental networks that respond to diverse environmental microbiota

Adaptation of organisms to environmental change may be facilitated by the creation of new genes. New genes without homologs in other lineages are known as taxonomically-restricted orphan genes and may result from divergence or de novo formation. Previously, we have extensively characterized the evolution and origin of such orphan genes in the nematode model organism Pristionchus pacificus. Here, we employ large-scale transcriptomics to establish potential functional associations and to measure the degree of transcriptional plasticity among orphan genes. Specifically, we analyzed 24 RNA-seq samples from adult P. pacificus worms raised on 24 different monoxenic bacterial cultures. Based on coexpression analysis, we identified 28 large modules that harbor 3,727 diplogastrid-specific orphan genes and that respond dynamically to different bacteria. These coexpression modules have distinct regulatory architecture and also exhibit differential expression patterns across development suggesting a link between bacterial response networks and development. Phylostratigraphy revealed a considerably high number of family- and even species-specific orphan genes in certain coexpression modules. This suggests that new genes are not attached randomly to existing cellular networks and that integration can happen very fast. Integrative analysis of protein domains, gene expression and ortholog data facilitated the assignments of biological labels for 22 coexpression modules with one of the largest, fast-evolving module being associated with spermatogenesis. In summary, this work presents the first functional annotation for thousands of P. pacificus orphan genes and reveals insights into their integration into environmentally responsive gene networks.

Chitin contributes to the formation of a feeding structure in a predatory nematode

Some nematode predators and parasites form teeth-like denticles that are histologically different from vertebrate teeth, but their biochemical composition remains elusive. Here, we show a role of chitin in the formation of teeth-like denticles in Pristionchus pacificus, a model system for studying predation and feeding structure plasticity. Pristionchus forms two alternative mouth morphs with one tooth or two teeth, respectively. The P. pacificus genome encodes two chitin synthases, with the highly conserved chs-2 gene being composed of 60 exons forming at least four isoforms. Generating CRISPR-Cas9-based gene knockouts, we found that Ppa-chs-2 mutations that eliminate the chitin-synthase domain are lethal. However, mutations in the C terminus result in viable but teethless worms, with severe malformation of the mouth. Similarly, treatment with the chitin-synthase inhibitor Nikkomycin Z also results in teethless animals. Teethless worms can feed on various bacterial food sources but are incapable of predation. High-resolution transcriptomics revealed that Ppa-chs-2 expression is controlled by the sulfatase-encoding developmental switch Ppa-eud-1. This study indicates a key role of chitin in the formation of teeth-like denticles and the complex feeding apparatus in nematodes.

Evolution and Diversity of TGF-β Pathways are Linked with Novel Developmental and Behavioral Traits

Transforming growth factor-β (TGF-β) signaling is essential for numerous biologic functions. It is a highly conserved pathway found in all metazoans including the nematode Caenorhabditis elegans, which has also been pivotal in identifying many components. Utilizing a comparative evolutionary approach, we explored TGF-β signaling in nine nematode species and revealed striking variability in TGF-β gene frequency across the lineage. Of the species analyzed, gene duplications in the DAF-7 pathway appear common with the greatest disparity observed in Pristionchus pacificus. Specifically, multiple paralogues of daf-3, daf-4 and daf-7 were detected. To investigate this additional diversity, we induced mutations in 22 TGF-β components and generated corresponding double, triple, and quadruple mutants revealing both conservation and diversification in function. Although the DBL-1 pathway regulating body morphology appears highly conserved, the DAF-7 pathway exhibits functional divergence, notably in some aspects of dauer formation. Furthermore, the formation of the phenotypically plastic mouth in P. pacificus is partially influenced through TGF-β with the strongest effect in Ppa-tag-68. This appears important for numerous processes in P. pacificus but has no known function in C. elegans. Finally, we observe behavioral differences in TGF-β mutants including in chemosensation and the establishment of the P. pacificus kin-recognition signal. Thus, TGF-β signaling in nematodes represents a stochastic genetic network capable of generating novel functions through the duplication and deletion of associated genes.

RNA-Sequencing of Heterorhabditis nematodes to identify factors involved in symbiosis with Photorhabdus bacteria

Background

Nematodes are a major group of soil inhabiting organisms. Heterorhabditis nematodes are insect-pathogenic nematodes and live in a close symbiotic association with Photorhabdus bacteria. Heterorhabditis-Photorhabdus pair offers a powerful and genetically tractable model to study animal-microbe symbiosis. It is possible to generate symbiont bacteria free (axenic) stages in Heterorhabditis. Here, we compared the transcriptome of symbiotic early-adult stage Heterorhabditis nematodes with axenic early-adult nematodes to determine the nematode genes and pathways involved in symbiosis with Photorhabdus bacteria.

Results

A de-novo reference transcriptome assembly of 95.7 Mb was created for H. bacteriophora by using all the reads. The assembly contained 46,599 transcripts with N50 value of 2,681 bp and the average transcript length was 2,054 bp. The differentially expressed transcripts were identified by mapping reads from symbiotic and axenic nematodes to the reference assembly. A total of 754 differentially expressed transcripts were identified in symbiotic nematodes as compared to the axenic nematodes. The ribosomal pathway was identified as the most affected among the differentially expressed transcripts. Additionally, 12,151 transcripts were unique to symbiotic nematodes. Endocytosis, cAMP signalling and focal adhesion were the top three enriched pathways in symbiotic nematodes, while a large number of transcripts coding for various responses against bacteria, such as bacterial recognition, canonical immune signalling pathways, and antimicrobial effectors could also be identified.

Conclusions

The symbiotic Heterorhabditis nematodes respond to the presence of symbiotic bacteria by expressing various transcripts involved in a multi-layered immune response which might represent non-systemic and evolved localized responses to maintain mutualistic bacteria at non-threatening levels. Subject to further functional validation of the identified transcripts, our findings suggest that Heterorhabditis nematode immune system plays a critical role in maintenance of symbiosis with Photorhabdus bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08952-4.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

The improved genome of the nematode Parapristionchus giblindavisi provides insights into lineage-specific gene family evolution

Nematodes such as Caenorhabditis elegans and Pristionchus pacificus are extremely successful model organisms for comparative biology. Several studies have shown that phenotypic novelty but also conserved processes are controlled by taxon-restricted genes. To trace back the evolution of such new or rapidly evolving genes, a robust phylogenomic framework is indispensable. Here, we present an improved version of the genome of Parapristionchus giblindavisi which is the only known member of the sister group of Pristionchus. Relative to the previous short-read assembly, the new genome is based on long reads and displays higher levels of contiguity, completeness, and correctness. Specifically, the number of contigs dropped from over 7,303 to 735 resulting in an N50 increase from 112 to 791 kb. We made use of the new genome to revisit the evolution of multiple gene families. This revealed Pristionchus-specific expansions of several environmentally responsive gene families and a Pristionchus-specific loss of the de novo purine biosynthesis pathway. Focusing on the evolution of sulfatases and sulfotransferases, which control the mouth form plasticity in P. pacificus, reveals differences in copy number and genomic configurations between the genera Pristionchus and Parapristionchus. Altogether, this demonstrates the utility of the P. giblindavisi genome to date and polarizes lineage-specific patterns.

Is Strongyloides stercoralis hyperinfection induced by glucocorticoids a result of both suppressed host immunity and altered parasite genetics?

The gastrointestinal (GI) nematode Strongyloides stercoralis (S.s.) causes human strongyloidiasis, a potentially life-threatening disease that currently affects over 600 million people globally. The uniquely pernicious aspect of S.s. infection, as compared to all other GI nematodes, is its autoinfective larval stage (L3a) that maintains a low-grade chronic infection, allowing undetectable persistence for decades. Infected individuals who are administered glucocorticoid therapy can develop a rapid and often lethal hyperinfection syndrome within days. Hyperinfection patients often present with dramatic increases in first- and second-stage larvae and L3a in their GI tract, with L3a widely disseminating throughout host organs leading to sepsis. How glucocorticoid administration drives hyperinfection remains a critical unanswered question; specifically, it is unknown whether these steroids promote hyperinfection through eliminating essential host protective mechanisms and/or through dysregulating parasite development. This current deficiency in understanding is largely due to the previous absence of a genetically defined mouse model that would support all S.s. life-cycle stages and the lack of successful approaches for S.s. genetic manipulation. However, there are currently new possibilities through the recent demonstration that immunodeficient NOD.Cg-Prkdc^scidIl2rg^tm1Wjl/SzJ (NSG) mice support sub-clinical infections that can be transformed to lethal hyperinfection syndrome following glucocorticoid administration. This is coupled with advances in transcriptomics, transgenesis, and gene inactivation strategies that now allow rigorous scientific inquiry into S.s. biology. We propose that combining in vivo manipulation of host immunity and deep immunoprofiling strategies with the latest advances in S.s. transcriptomics, piggyBac transposon-mediated transgene insertion, and CRISPR/Cas-9-mediated gene inactivation will facilitate new insights into the mechanisms that could be targeted to block lethality in humans with S.s. hyperinfection.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Synergistic interaction of gut microbiota enhances the growth of nematode through neuroendocrine signaling

Animals are associated with a diverse bacterial community that impacts host physiology. It is well known that nutrients and enzymes synthesized by bacteria largely expand host metabolic capacity. Bacteria also impact a wide range of animal physiology that solely depends on host genetics through direct interaction. However, studying the synergistic effects of the bacterial community remains challenging due to its complexity. The omnivorous nematode Pristionchus pacificus has limited digestive efficiency on bacteria. Therefore, we established a bacterial collection that represents the natural gut microbiota that are resistant to digestion. Using this collection, we show that the bacterium Lysinibacillus xylanilyticus by itself provides limited nutritional value, but in combination with Escherichia coli, it significantly promotes life-history traits of P. pacificus by regulating the neuroendocrine peptide in sensory neurons. This gut-to-brain communication depends on undigested L. xylanilyticus providing Pristionchus nematodes a specific fitness advantage to compete with nematodes that rupture bacteria efficiently. Using RNA-seq and CRISPR-induced mutants, we show that 1-h exposure to L. xylanilyticus is sufficient to stimulate the expression of daf-7-type TGF-β signaling ligands, which induce a global transcriptome change. In addition, several effects of L. xylanilyticus depend on TGF-β signaling, including olfaction, body size regulation, and a switch of energy allocation from lipid storage to reproduction. Our results reveal the beneficial effects of a gut bacterium to modify life-history traits and maximize nematode survival in natural habitats.

Vitamin B12 and predatory behavior in nematodes

The round worms or nematodes are the largest phylum of animals with an estimated species number of more than one million. Nematodes have invaded all ecosystems and are known from all continents including Antarctica. Parasitic species infest plants, animals and humans often with high host-specificity. Free-living species are known from marine, fresh water and soil systems, the latter of which contain many culturable species. This includes Caenorhabditis elegans, a species that was developed as one of the most prominent model systems in modern biology since the 1960ies. Pristionchus pacificus is a second nematode model organism that can easily be cultured in the laboratory. This species shows a number of complex traits including omnivorous feeding and the capability of predation on other nematodes. Predation depends on the formation of teeth-like denticles in the mouth of P. pacificus, structures unknown from C. elegans and most other nematodes. Here, we review the current knowledge about the role of vitamin B₁₂ for the predatory behavior in P. pacificus and correlate its role with that on the physiology and development in C. elegans.

The oscillating Mucin-type protein DPY-6 has a conserved role in nematode mouth and cuticle formation

Nematodes show an extraordinary diversity of mouth structures and strikingly different feeding strategies, which has enabled an invasion of all ecosystems. However, nearly nothing is known about the structural and molecular architecture of the nematode mouth (stoma). Pristionchus pacificus is an intensively studied nematode that exhibits unique life history traits, including predation, teeth-like denticle formation, and mouth-form plasticity. Here, we used a large-scale genetic screen to identify genes involved in mouth formation. We identified Ppa-dpy-6 to encode a Mucin-type hydrogel-forming protein that is macroscopically involved in the specification of the cheilostom, the anterior part of the mouth. We used a recently developed protocol for geometric morphometrics of miniature animals to characterize these defects further and found additional defects that affect mouth form, shape, and size resulting in an overall malformation of the mouth. Additionally, Ppa-dpy-6 is shorter than wild-type with a typical Dumpy phenotype, indicating a role in the formation of the external cuticle. This concomitant phenotype of the cheilostom and cuticle provides the first molecular support for the continuity of these structures and for the separation of the cheilostom from the rest of the stoma. In Caenorhabditis elegans, dpy-6 was an early mapping mutant but its molecular identity was only determined during genome-wide RNAi screens and not further investigated. Strikingly, geometric morphometric analysis revealed previously unrecognized cheilostom and gymnostom defects in Cel-dpy-6 mutants. Thus, the Mucin-type protein DPY-6 represents to the best of our knowledge, the first protein involved in nematode mouth formation with a conserved role in cuticle deposition. This study opens new research avenues to characterize the molecular composition of the nematode mouth, which is associated with extreme ecological diversification.