Single-Molecule Sequencing Reveals the Chromosome-Scale Genomic Architecture of the Nematode Model Organism Pristionchus pacificus

Microbiota succession influences nematode physiology in a beetle microcosm ecosystem

Unravelling the multifaceted and bidirectional interactions between microbiota and host physiology represents a major scientific challenge. Here, we utilise the nematode model, Pristionchus pacificus, coupled to a laboratory-simulated decay process of its insect host, to mimic natural microbiota succession and investigate associated tripartite interactions. Metagenomics reveal that during initial decay stages, the population of vitamin B-producing bacteria diminishes, potentially due to a preferential selection by nematodes. As decay progresses to nutrient-depleted stages, bacteria with smaller genomes producing less nutrients become more prevalent. Lipid utilisation and dauer formation, representing key nematode survival strategies, are influenced by microbiota changes. Additionally, horizontally acquired cellulases extend the nematodes' reproductive phase due to more efficient foraging. Lastly, the expressions of Pristionchus species-specific genes are more responsive to natural microbiota compared to conserved genes, suggesting their importance in the organisms' adaptation to its ecological niche. In summary, we show the importance of microbial successions and their reciprocal interaction with nematodes for insect decay in semi-artificial ecosystems.

Chromosome fusion and programmed DNA elimination shape karyotypes of nematodes

An increasing number of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and increased numbers of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the nematode Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among other nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these nematodes, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes. These karyotype changes may lead to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms are dynamic and may play a role during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

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.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Chromosome fusion and programmed DNA elimination shape karyotypes of parasitic nematodes

A growing list of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and an increased number of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the horse parasite Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among ascarid nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these parasites, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes, leading to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms play a dynamic role in the Parascaris germline chromosome during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

Kin-recognition and predation shape collective behaviors in the cannibalistic nematode Pristionchus pacificus

Kin-recognition is observed across diverse species forming an important behavioral adaptation influencing organismal interactions. In many species, the molecular mechanisms involved are difficult to characterize, but in the nematode Pristionchus pacificus molecular components regulating its kin-recognition system have been identified. These determine its predatory behaviors towards other con-specifics which prevents the killing and cannibalization of kin. Importantly, their impact on other interactions including collective behaviors is unknown. Here, we explored a high altitude adapted clade of this species which aggregates abundantly under laboratory conditions, to investigate the influence of the kin-recognition system on their group behaviours. By utilizing pairwise aggregation assays between distinct strains of P. pacificus with differing degrees of genetic relatedness, we observe aggregation between kin but not distantly related strains. In assays between distantly related strains, the aggregation ratio is frequently reduced. Furthermore, abolishing predation behaviors through CRISPR/Cas9 induced mutations in Ppa-nhr-40 result in rival strains successfully aggregating together. Finally, as Caenorhabditis elegans are found naturally occurring with P. pacificus, we also explored aggregation events between these species. Here, aggregates were dominated by P. pacificus with the presence of only a small number of predators proving sufficient to disrupt C. elegans aggregation dynamics. Thus, aggregating strains of P. pacificus preferentially group with kin, revealing competition and nepotism as previously unknown components influencing collective behaviors in nematodes.

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.

Divergent combinations of cis-regulatory elements control the evolution of phenotypic plasticity

The widespread occurrence of phenotypic plasticity across all domains of life demonstrates its evolutionary significance. However, how plasticity itself evolves and how it contributes to evolution is poorly understood. Here, we investigate the predatory nematode Pristionchus pacificus with its feeding structure plasticity using recombinant-inbred-line and quantitative-trait-locus (QTL) analyses between natural isolates. We show that a single QTL at a core developmental gene controls the expression of the cannibalistic morph. This QTL is composed of several cis-regulatory elements. Through CRISPR/Cas-9 engineering, we identify copy number variation of potential transcription factor binding sites that interacts with a single intronic nucleotide polymorphism. Another intronic element eliminates gene expression altogether, mimicking knockouts of the locus. Comparisons of additional isolates further support the rapid evolution of these cis-regulatory elements. Finally, an independent QTL study reveals evidence for parallel evolution at the same locus. Thus, combinations of cis-regulatory elements shape plastic trait expression and control nematode cannibalism.

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.

Influence of environmental temperature on mouth-form plasticity in Pristionchus pacificus acts through daf-11-dependent cGMP signaling

Mouth-form plasticity in the nematode Pristionchus pacificus has become a powerful system to identify the genetic and molecular mechanisms associated with developmental (phenotypic) plasticity. In particular, the identification of developmental switch genes that can sense environmental stimuli and reprogram developmental processes has confirmed long-standing evolutionary theory. However, how these genes are involved in the direct sensing of the environment, or if the switch genes act downstream of another, primary environmental sensing mechanism, remains currently unknown. Here, we study the influence of environmental temperature on mouth-form plasticity. We find that environmental temperature does influence mouth-form plasticity in most of the 10 wild isolates of P. pacificus tested in this study. We used one of these strains, P. pacificus RSA635, for detailed molecular analysis. Using forward and reverse genetic technology including CRISPR/Cas9, we show that mutations in the guanylyl cyclase Ppa-daf-11, the Ppa-daf-25/AnkMy2, and the cyclic nucleotide-gated channel Ppa-tax-2 eliminate the response to elevated temperatures. Together, our study indicates that DAF-11, DAF-25, and TAX-2 have been co-opted for environmental sensing during mouth-form plasticity regulation in P. pacificus.

Experimental and theoretical support for costs of plasticity and phenotype in a nematode cannibalistic trait

Developmental plasticity is the ability of a genotype to express multiple phenotypes under different environmental conditions and has been shown to facilitate the evolution of novel traits. However, while the associated cost of plasticity, i.e., the loss in fitness due to the ability to express plasticity in response to environmental change, and the cost of phenotype, i.e., the loss of fitness due to expressing a fixed phenotype across environments, have been theoretically predicted, empirically such costs remain poorly documented and little understood. Here, we use a plasticity model system, hermaphroditic nematode Pristionchus pacificus, to experimentally measure these costs in wild isolates under controlled laboratory conditions. P. pacificus can develop either a bacterial feeding or predatory mouth morph in response to different external stimuli, with natural variation of mouth-morph ratios between strains. We first demonstrated the cost of phenotype by analyzing fecundity and developmental speed in relation to mouth morphs across the P. pacificus phylogenetic tree. Then, we exposed P. pacificus strains to two distinct microbial diets that induce strain-specific mouth-form ratios. Our results indicate that the plastic strain does shoulder a cost of plasticity, i.e., the diet-induced predatory mouth morph is associated with reduced fecundity and slower developmental speed. In contrast, the non-plastic strain suffers from the cost of phenotype since its phenotype does not change to match the unfavorable bacterial diet but shows increased fitness and higher developmental speed on the favorable diet. Furthermore, using a stage-structured population model based on empirically derived life history parameters, we show how population structure can alleviate the cost of plasticity in P. pacificus. The results of the model illustrate the extent to which the costs associated with plasticity and its effect on competition depend on ecological factors. This study provides support for costs of plasticity and phenotype based on empirical and modeling approaches.

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.

A New Hope: A Hermaphroditic Nematode Enables Analysis of a Recent Whole Genome Duplication Event

Whole genome duplication (WGD) is often considered a major driver of evolution that leads to phenotypic novelties. However, the importance of WGD for evolution is still controversial because most documented WGD events occurred anciently and few experimental systems amenable to genetic analysis are available. Here, we report a recent WGD event in the hermaphroditic nematode Allodiplogaster sudhausi and present a comparison with a gonochoristic (male/female) sister species that did not undergo WGD. Self-fertilizing reproduction of A. sudhausi makes it amenable to functional analysis and an ideal system to study WGD events. We document WGD in A. sudhausi through karyotype analysis and whole genome sequencing, the latter of which allowed us to 1) identify functional bias in retention of protein domains and metabolic pathways, 2) show most duplicate genes are under evolutionary constraint, 3) show a link between sequence and expression divergence, and 4) characterize differentially expressed duplicates. We additionally show WGD is associated with increased body size and an abundance of repeat elements (36% of the genome), including a recent expansion of the DNA-hAT/Ac transposon family. Finally, we demonstrate the use of CRISPR/Cas9 to generate mutant knockouts, whereby two WGD-derived duplicate genes display functional redundancy in that they both need to be knocked out to generate a phenotype. Together, we present a novel experimental system that is convenient for examining and characterizing WGD-derived genes both computationally and functionally.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

Analysis of repeat elements in the Pristionchus pacificus genome reveals an ancient invasion by horizontally transferred transposons

BACKGROUND

Repetitive sequences and mobile elements make up considerable fractions of individual genomes. While transposition events can be detrimental for organismal fitness, repetitive sequences form an enormous reservoir for molecular innovation. In this study, we aim to add repetitive elements to the annotation of the Pristionchus pacificus genome and assess their impact on novel gene formation.

RESULTS

Different computational approaches define up to 24% of the P. pacificus genome as repetitive sequences. While retroelements are more frequently found at the chromosome arms, DNA transposons are distributed more evenly. We found multiple DNA transposons, as well as LTR and LINE elements with abundant evidence of expression as single-exon transcripts. When testing whether transposons disproportionately contribute towards new gene formation, we found that roughly 10-20% of genes across all age classes overlap transposable elements with the strongest trend being an enrichment of low complexity regions among the oldest genes. Finally, we characterized a horizontal gene transfer of Zisupton elements into diplogastrid nematodes. These DNA transposons invaded nematodes from eukaryotic donor species and experienced a recent burst of activity in the P. pacificus lineage.

CONCLUSIONS

The comprehensive annotation of repetitive elements in the P. pacificus genome builds a resource for future functional genomic analyses as well as for more detailed investigations of molecular innovations.

Robust, versatile DNA FISH probes for chromosome-specific repeats in Caenorhabditis elegans and Pristionchus pacificus

Repetitive DNA sequences are useful targets for chromosomal fluorescence in situ hybridization. We analyzed recent genome assemblies of Caenorhabditis elegans and Pristionchus pacificus to identify tandem repeats with a unique genomic localization. Based on these findings, we designed and validated sets of oligonucleotide probes for each species targeting at least 1 locus per chromosome. These probes yielded reliable fluorescent signals in different tissues and can easily be combined with the immunolocalization of cellular proteins. Synthesis and labeling of these probes are highly cost-effective and require no hands-on labor. The methods presented here can be easily applied in other model and nonmodel organisms with a sequenced genome.

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.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

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.

Nematode chromosomes

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

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.

Horizontally acquired cellulases assist the expansion of dietary range in Pristionchus nematodes

Horizontal gene transfer (HGT) enables the acquisition of novel traits via non-Mendelian inheritance of genetic material. HGT plays a prominent role in the evolution of prokaryotes, whereas in animals, HGT is rare and its functional significance is often uncertain. Here, we investigate horizontally acquired cellulase genes in the free-living nematode model organism Pristionchus pacificus. We show that these cellulase genes 1) are likely of eukaryotic origin, 2) are expressed, 3) have protein products that are secreted and functional, and 4) result in endo-cellulase activity. Using CRISPR/Cas9, we generated an octuple cellulase mutant, which lacks all eight cellulase genes and cellulase activity altogether. Nonetheless, this cellulase-null mutant is viable and therefore allows a detailed analysis of a gene family that was horizontally acquired. We show that the octuple cellulase mutant has associated fitness costs with reduced fecundity and slower developmental speed. Furthermore, by using various Escherichia coli K-12 strains as a model for cellulosic biofilms, we demonstrate that cellulases facilitate the procurement of nutrients from bacterial biofilms. Together, our analysis of cellulases in Pristionchus provides comprehensive evidence from biochemistry, genetics, and phylogeny, which supports the integration of horizontally acquired genes into the complex life history strategy of this soil nematode.

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.

hlh-12, a gene that is necessary and sufficient to promote migration of gonadal regulatory cells in Caenorhabditis elegans, evolved within the Caenorhabditis clade

Specialized cells of the somatic gonad primordium of nematodes play important roles in the final form and function of the mature gonad. Caenorhabditis elegans hermaphrodites are somatic females that have a two-armed, U-shaped gonad that connects to the vulva at the midbody. The outgrowth of each gonad arm from the somatic gonad primordium is led by two female distal tip cells (fDTCs), while the anchor cell (AC) remains stationary and central to coordinate uterine and vulval development. The bHLH protein HLH-2 and its dimerization partners LIN-32 and HLH-12 had previously been shown to be required for fDTC specification. Here, we show that ectopic expression of both HLH-12 and LIN-32 in cells with AC potential transiently transforms them into fDTC-like cells. Furthermore, hlh-12 was known to be required for the fDTCs to sustain gonad arm outgrowth. Here, we show that ectopic expression of HLH-12 in the normally stationary AC causes displacement from its normal position and that displacement likely results from activation of the leader program of fDTCs because it requires genes necessary for gonad arm outgrowth. Thus, HLH-12 is both necessary and sufficient to promote gonadal regulatory cell migration. As differences in female gonadal morphology of different nematode species reflect differences in the fate or migratory properties of the fDTCs or of the AC, we hypothesized that evolutionary changes in the expression of hlh-12 may underlie the evolution of such morphological diversity. However, we were unable to identify an hlh-12 ortholog outside of Caenorhabditis. Instead, by performing a comprehensive phylogenetic analysis of all Class II bHLH proteins in multiple nematode species, we found that hlh-12 evolved within the Caenorhabditis clade, possibly by duplicative transposition of hlh-10. Our analysis suggests that control of gene regulatory hierarchies for gonadogenesis can be remarkably plastic during evolution without adverse phenotypic consequence.

Frequent lineage-specific substitution rate changes support an episodic model for protein evolution

Since the inception of the molecular clock model for sequence evolution, the investigation of protein divergence has revolved around the question of a more or less constant change of amino acid sequences, with specific overall rates for each family. Although anomalies in clock-like divergence are well known, the assumption of a constant decay rate for a given protein family is usually taken as the null model for protein evolution. However, systematic tests of this null model at a genome-wide scale have lagged behind, despite the databases’ enormous growth. We focus here on divergence rate comparisons between very closely related lineages since this allows clear orthology assignments by synteny and reliable alignments, which are crucial for determining substitution rate changes. We generated a high-confidence dataset of syntenic orthologs from four ape species, including humans. We find that despite the appearance of an overall clock-like substitution pattern, several hundred protein families show lineage-specific acceleration and deceleration in divergence rates, or combinations of both in different lineages. Hence, our analysis uncovers a rather dynamic history of substitution rate changes, even between these closely related lineages, implying that one should expect that a large fraction of proteins will have had a history of episodic rate changes in deeper phylogenies. Furthermore, each of the lineages has a separate set of particularly fast diverging proteins. The genes with the highest percentage of branch-specific substitutions are ADCYAP1 in the human lineage (9.7%), CALU in chimpanzees (7.1%), SLC39A14 in the internal branch leading to humans and chimpanzees (4.1%), RNF128 in gorillas (9%), and S100Z in gibbons (15.2%). The mutational pattern in ADCYAP1 suggests a biased mutation process, possibly through asymmetric gene conversion effects. We conclude that a null model of constant change can be problematic for predicting the evolutionary trajectories of individual proteins.

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

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

Sex or cannibalism: Polyphenism and kin recognition control social action strategies in nematodes

Resource polyphenisms, where single genotypes produce alternative feeding strategies in response to changing environments, are thought to be facilitators of evolutionary novelty. However, understanding the interplay between environment, morphology, and behavior and its significance is complex. We explore a radiation of Pristionchus nematodes with discrete polyphenic mouth forms and associated microbivorous versus cannibalistic traits. Notably, comparing 29 Pristionchus species reveals that reproductive mode strongly correlates with mouth-form plasticity. Male-female species exhibit the microbivorous morph and avoid parent-offspring conflict as indicated by genetic hybrids. In contrast, hermaphroditic species display cannibalistic morphs encouraging competition. Testing predation between 36 co-occurring strains of the hermaphrodite P. pacificus showed that killing inversely correlates with genomic relatedness. These empirical data together with theory reveal that polyphenism (plasticity), kin recognition, and relatedness are three major factors that shape cannibalistic behaviors. Thus, developmental plasticity influences cooperative versus competitive social action strategies in diverse animals.

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

High-resolution spatial and temporal maps of gene expression have facilitated a comprehensive understanding of animal development and evolution. In nematodes, the small body size represented a major challenge for such studies, but recent advancements have helped overcome this limitation. Here, we have implemented single worm transcriptomics (SWT) in the nematode model organism Pristionchus pacificus to provide a high-resolution map of the developmental transcriptome. We selected 38 time points from hatching of the J2 larvae to young adults to perform transcriptome analysis over 60 h of postembryonic development. A mean sequencing depth of 4.5 million read pairs allowed the detection of more than 23,135 (80%) of all genes. Nearly 3000 (10%) genes showed oscillatory expression with discrete expression levels, phases, and amplitudes. Gene age analysis revealed an overrepresentation of ancient gene classes among oscillating genes, and around one-third of them have 1:1 orthologs in C. elegans. One important gene family overrepresented among oscillating genes is collagens. Several of these collagen genes are regulated by the developmental switch gene eud-1, indicating a potential function in the regulation of mouth-form plasticity, a key developmental process in this facultative predatory nematode. Together, our analysis provides (1) an updated protocol for SWT in nematodes that is applicable to many microscopic species, (2) a 1- to 2-h high-resolution catalog of P. pacificus gene expression throughout postembryonic development, and (3) a comparative analysis of oscillatory gene expression between the two model organisms P. pacificus and C. elegans and associated evolutionary dynamics.

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.

Steroid hormone pathways coordinate developmental diapause and olfactory remodeling in Pristionchus pacificus

Developmental and behavioral plasticity allow animals to prioritize alternative genetic programs during fluctuating environments. Behavioral remodeling may be acute in animals that interact with host organisms, since reproductive adults and the developmentally arrested larvae often have different ethological needs for chemical stimuli. To understand the genes that coordinate the development and host-seeking behavior, we used the entomophilic nematode Pristionchus pacificus to characterize dauer-constitutive mutants (Daf-c) that inappropriately enter developmental diapause to become dauer larvae. We found two Daf-c loci with dauer-constitutive and cuticle exsheathment phenotypes that can be rescued by the feeding of Δ7-dafachronic acid, and that are dependent on the conserved canonical steroid hormone receptor Ppa-DAF-12. Specifically at one locus, deletions in the sole hydroxysteroid dehydrogenase (HSD) in P. pacificus resulted in Daf-c phenotypes. Ppa-hsd-2 is expressed in the canal-associated neurons (CANs) and excretory cells whose homologous cells in Caenorhabditis elegans are not known to be involved in the dauer decision. While in wildtype only dauer larvae are attracted to host odors, hsd-2 mutant adults show enhanced attraction to the host beetle pheromone, along with ectopic activation of a marker for putative olfactory neurons, Ppa-odr-3. Surprisingly, this enhanced odor attraction acts independently of the Δ7-DA/DAF-12 module, suggesting that Ppa-HSD-2 may be responsible for several steroid hormone products involved in coordinating the dauer decision and host-seeking behavior in P. pacificus.

The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes

Background

The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans.

Results

Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated.

Conclusions

This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07529-x.

Different combinations of serotonin receptors regulate predatory and bacterial feeding behaviors in the nematode Pristionchus pacificus

Feeding behavior is one of the most fundamental behaviors in animals, and regulation of this behavior is critical for proper food intake. The nematode Pristionchus pacificus exhibits dimorphism in feeding behavior, bacterial feeding and predatory feeding on other nematodes, and the latter behavior is assumed to be an evolutionarily novel behavior. Both types of feeding behavior are modulated by serotonin; however, the downstream mechanism that modulates these behaviors is still to be clarified. Here, we focused on serotonin receptors and examined their expression patterns in P. pacificus. We also generated knockout mutants of the serotonin receptors using the CRISPR/Cas9 system and examined feeding behaviors. We found that Ppa-ser-5 mutants and the Ppa-ser-1; Ppa-ser-7 double mutant decreased predation. Detailed observation of the pharyngeal movement revealed that the Ppa-ser-1; Ppa-ser-7 double mutant reduces tooth movement, which is required for efficient predatory feeding. Conversely, Ppa-ser-7 and Ppa-mod-1 mutants decreased bacterial feeding. This study revealed that specific combinations of serotonin receptors are essential for the modulation of these distinct feeding behaviors, providing insight into the evolution of neural pathways to regulate novel feeding behavior.

Comparative genomics and community curation further improve gene annotations in the nematode Pristionchus pacificus

Background

Nematode model organisms such as Caenorhabditis elegans and Pristionchus pacificus are powerful systems for studying the evolution of gene function at a mechanistic level. However, the identification of P. pacificus orthologs of candidate genes known from C. elegans is complicated by the discrepancy in the quality of gene annotations, a common problem in nematode and invertebrate genomics.

Results

Here, we combine comparative genomic screens for suspicious gene models with community-based curation to further improve the quality of gene annotations in P. pacificus. We extend previous curations of one-to-one orthologs to larger gene families and also orphan genes. Cross-species comparisons of protein lengths, screens for atypical domain combinations and species-specific orphan genes resulted in 4311 candidate genes that were subject to community-based curation. Corrections for 2946 gene models were implemented in a new version of the P. pacificus gene annotations. The new set of gene annotations contains 28,896 genes and has a single copy ortholog completeness level of 97.6%.

Conclusions

Our work demonstrates the effectiveness of comparative genomic screens to identify suspicious gene models and the scalability of community-based approaches to improve the quality of thousands of gene models. Similar community-based approaches can help to improve the quality of gene annotations in other invertebrate species, including parasitic nematodes.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.

Genomic studies on Strongyloides stercoralis in northern and western Thailand

Background

Strongyloidiasis is a soil borne helminthiasis, which in most cases is caused by Strongyloides stercoralis. Human infections with S. fuelleborni fuelleborni and S. fuelleborni kellyi also occur. Although up to 370 million people are currently estimated to be infected with S. stercoralis, this parasite is frequently overlooked. Strongyloides stercoralis is prevalent among humans in Thailand; however, S. fuelleborni fuelleborni has also been reported. Three recent genomic studies of individual S. stercoralis worms found genetically diverse populations of S. stercoralis, with comparably low heterozygosity in Cambodia and Myanmar, and less diverse populations with high heterozygosity in Japan and southern China that presumably reproduce asexually.

Methods

We isolated individual Strongyloides spp. from different localities in northern and western Thailand and determined their nuclear small ribosomal subunit rDNA (18S rDNA, SSU), in particular the hypervariable regions I and IV (HVR-I and HVR-IV), mitochondrial cytochrome c oxidase subunit 1 (cox1) and for a subset whole genome sequences. These sequences were then compared with each other and with published sequences from different geographical locations.

Results

All 237 worms isolated from 16 different human hosts were S. stercoralis, no S. fuelleborni was found. All worms had the common S. stercoralis SSU HVR IV haplotype A. Two different SSU HVR I haplotypes (I and II), both previously described in S. stercoralis, were found. No animal heterozygous for the two haplotypes was identified. Among the twelve cox1 haplotypes found, five had not been previously described. Based upon the mitochondrial cox1 and the nuclear whole genome sequences, S. stercoralis in Thailand was phylogenetically intermixed with the samples from other Southeast Asian countries and did not form its own branch. The genomic heterozygosity was even slightly lower than in the samples from the neighboring countries.

Conclusions

In our sample from humans, all Strongyloides spp. were S. stercoralis. The S. stercoralis from northern and western Thailand appear to be part of a diverse, intermixing continental Southeast Asian population. No obvious indication for genetic sub-structuring of S. stercoralis within Thailand or within the Southeast Asian peninsula was detected.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

Conserved nuclear hormone receptors controlling a novel plastic trait target fast-evolving genes expressed in a single cell

Environment shapes development through a phenomenon called developmental plasticity. Deciphering its genetic basis has potential to shed light on the origin of novel traits and adaptation to environmental change. However, molecular studies are scarce, and little is known about molecular mechanisms associated with plasticity. We investigated the gene regulatory network controlling predatory vs. non-predatory dimorphism in the nematode Pristionchus pacificus and found that it consists of genes of extremely different age classes. We isolated mutants in the conserved nuclear hormone receptor nhr-1 with previously unseen phenotypic effects. They disrupt mouth-form determination and result in animals combining features of both wild-type morphs. In contrast, mutants in another conserved nuclear hormone receptor nhr-40 display altered morph ratios, but no intermediate morphology. Despite divergent modes of control, NHR-1 and NHR-40 share transcriptional targets, which encode extracellular proteins that have no orthologs in Caenorhabditis elegans and result from lineage-specific expansions. An array of transcriptional reporters revealed co-expression of all tested targets in the same pharyngeal gland cell. Major morphological changes in this gland cell accompanied the evolution of teeth and predation, linking rapid gene turnover with morphological innovations. Thus, the origin of feeding plasticity involved novelty at the level of genes, cells and behavior.

Crowdsourcing and the feasibility of manual gene annotation: A pilot study in the nematode Pristionchus pacificus

Nematodes such as Caenorhabditis elegans are powerful systems to study basically all aspects of biology. Their species richness together with tremendous genetic knowledge from C. elegans facilitate the evolutionary study of biological functions using reverse genetics. However, the ability to identify orthologs of candidate genes in other species can be hampered by erroneous gene annotations. To improve gene annotation in the nematode model organism Pristionchus pacificus, we performed a genome-wide screen for C. elegans genes with potentially incorrectly annotated P. pacificus orthologs. We initiated a community-based project to manually inspect more than two thousand candidate loci and to propose new gene models based on recently generated Iso-seq and RNA-seq data. In most cases, misannotation of C. elegans orthologs was due to artificially fused gene predictions and completely missing gene models. The community-based curation raised the gene count from 25,517 to 28,036 and increased the single copy ortholog completeness level from 86% to 97%. This pilot study demonstrates how even small-scale crowdsourcing can drastically improve gene annotations. In future, similar approaches can be used for other species, gene sets, and even larger communities thus making manual annotation of large parts of the genome feasible.

A Developmental Switch Generating Phenotypic Plasticity Is Part of a Conserved Multi-gene Locus

Switching between alternative complex phenotypes is often regulated by "supergenes," polymorphic clusters of linked genes such as in butterfly mimicry. In contrast, phenotypic plasticity results in alternative complex phenotypes controlled by environmental influences rather than polymorphisms. Here, we show that the developmental switch gene regulating predatory versus non-predatory mouth-form plasticity in the nematode Pristionchus pacificus is part of a multi-gene locus containing two sulfatases and two α-N-acetylglucosaminidases (nag). We provide functional characterization of all four genes, using CRISPR-Cas9-based reverse genetics, and show that nag genes and the previously identified eud-1/sulfatase have opposing influences. Members of the multi-gene locus show non-overlapping neuronal expression and epistatic relationships. The locus architecture is conserved in the entire genus Pristionchus. Interestingly, divergence between paralogs is counteracted by gene conversion, as inferred from phylogenies and genotypes of CRISPR-Cas9-induced mutants. Thus, we found that physical linkage accompanies regulatory linkage between switch genes controlling plasticity in P. pacificus.

Evolution of neuronal anatomy and circuitry in two highly divergent nematode species

The nematodes C. elegans and P. pacificus populate diverse habitats and display distinct patterns of behavior. To understand how their nervous systems have diverged, we undertook a detailed examination of the neuroanatomy of the chemosensory system of P. pacificus. Using independent features such as cell body position, axon projections and lipophilic dye uptake, we have assigned homologies between the amphid neurons, their first-layer interneurons, and several internal receptor neurons of P. pacificus and C. elegans. We found that neuronal number and soma position are highly conserved. However, the morphological elaborations of several amphid cilia are different between them, most notably in the absence of ‘winged’ cilia morphology in P. pacificus. We established a synaptic wiring diagram of amphid sensory neurons and amphid interneurons in P. pacificus and found striking patterns of conservation and divergence in connectivity relative to C. elegans, but very little changes in relative neighborhood of neuronal processes. These findings demonstrate the existence of several constraints in patterning the nervous system and suggest that major substrates for evolutionary novelty lie in the alterations of dendritic structures and synaptic connectivity.

Nerve cells, also called neurons, are responsible both for sensing signals from the environment and for determining how organisms react. This means that the unique features of an animal’s nervous system underpin its characteristic behaviors. Comparing the anatomy of the nervous systems in different animals could therefore yield valuable insights into how structural and behavioral differences emerge over time.

Behavioral variation often occurs even in similar-looking animals. One example is a group of microscopic worms, called nematodes. Although many nematode species exist, their overall body plans are the same, and the worms of each species contain a fixed number of cells. Despite these apparent similarities, different species of nematodes inhabit a variety of environments and may respond differently to the same signals.

The main sensory organs in nematodes are called the amphid sensilla. They are used to detect chemicals, as well as other inputs from the environment such as temperature and pheromones from other nematodes. Although researchers have often speculated that the number of cells in these organs and their arrangement are broadly the same across species, their anatomy had not been studied in detail.

Hong, Riebesell et al. compared the detailed structure and genetic features of the sensory systems in two distantly related species of nematode worms, Pristionchus pacificus and Caenorhabditis elegans. These two species behave in different ways, for example, P. pacificus is usually found in association with different species of beetles, while C. elegans is free-living and usually found on rotting fruit. By comparing the two, Hong, Riebesell et al. wanted to determine whether the diverse behaviors observed in the two species could be determined by differences between their sensory systems.

Experiments using electron microscopy yielded several thousand high resolution images spanning the entire sensory organ. These images were then used to create detailed reconstructions of the sensory nervous system in each worm species, demonstrating that both species had the same number of sensory nerve cells, allowing one-to-one comparisons between them.

Further analysis showed that while the overall structure of the neuronal connections remains the same between the two species, the neurons in P. pacificus made more diverse connections than those in C. elegans. Detailed studies of gene activity also revealed that neurons in each species switched on a slightly different group of genes, possibly indicating that each type of worm processes sensory signals in different ways.

These results shed new light on how nervous systems in related species can change over time without any change in neuron count. In the future, a better understanding of these changes could link the evolution of the nervous system to the emergence of different behaviors, in both simple and more complex organisms.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

Multiple plasticity regulators reveal targets specifying an induced predatory form in nematodes

The ability to translate a single genome into multiple phenotypes, or developmental plasticity, defines how phenotype derives from more than just genes. However, to study the evolutionary targets of plasticity and their evolutionary fates, we need to understand how genetic regulators of plasticity control downstream gene expression. Here, we have identified a transcriptional response specific to polyphenism (i.e., discrete plasticity) in the nematode Pristionchus pacificus. This species produces alternative resource-use morphs - microbivorous and predatory forms, differing in the form of their teeth, a morphological novelty - as influenced by resource availability. Transcriptional profiles common to multiple polyphenism-controlling genes in P. pacificus reveal a suite of environmentally sensitive loci, or ultimate target genes, that make up an induced developmental response. Additionally, in vitro assays show that one polyphenism regulator, the nuclear receptor (NR) NHR-40, physically binds to promoters with putative HNF4⍺ (the NR class including NHR-40) binding sites, suggesting this receptor may directly regulate genes that describe alternative morphs. Among differentially expressed genes were morph-limited genes, highlighting factors with putative "on-off" function in plasticity regulation. Further, predatory morph-biased genes included candidates - namely, all four P. pacificus homologs of Hsp70, which have HNF4⍺ motifs - whose natural variation in expression matches phenotypic differences among P. pacificus wild isolates. In summary, our study links polyphenism regulatory loci to the transcription producing alternative forms of a morphological novelty. Consequently, our findings establish a platform for determining how specific regulators of morph-biased genes may influence selection on plastic phenotypes.

De Novo, Divergence, and Mixed Origin Contribute to the Emergence of Orphan Genes in Pristionchus Nematodes

Homology is a fundamental concept in comparative biology. It is extensively used at the sequence level to make phylogenetic hypotheses and functional inferences. Nonetheless, the majority of eukaryotic genomes contain large numbers of orphan genes lacking homologs in other taxa. Generally, the fraction of orphan genes is higher in genomically undersampled clades, and in the absence of closely related genomes any hypothesis about their origin and evolution remains untestable. Previously, we sequenced ten genomes with an underlying ladder-like phylogeny to establish a phylogenomic framework for studying genome evolution in diplogastrid nematodes. Here, we use this deeply sampled data set to understand the processes that generate orphan genes in our focal species Pristionchus pacificus Based on phylostratigraphic analysis and additional bioinformatic filters, we obtained 29 high-confidence candidate genes for which mechanisms of orphan origin were proposed based on manual inspection. This revealed diverse mechanisms including annotation artifacts, chimeric origin, alternative reading frame usage, and gene splitting with subsequent gain of de novo exons. In addition, we present two cases of complete de novo origination from non-coding regions, which represents one of the first reports of de novo genes in nematodes. Thus, we conclude that de novo emergence, divergence, and mixed mechanisms contribute to novel gene formation in Pristionchus nematodes.

Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

Characterization of a non-sexual population of Strongyloides stercoralis with hybrid 18S rDNA haplotypes in Guangxi, Southern China

Strongyloidiasis is a much-neglected but sometimes fatal soil born helminthiasis. The causing agent, the small intestinal parasitic nematode Strongyloides stercoralis can reproduce sexually through the indirect/heterogonic life cycle, or asexually through the auto-infective or the direct/homogonic life cycles. Usually, among the progeny of the parasitic females both, parthenogenetic parasitic (females only) and sexual free-living (females and males) individuals, are present simultaneously. We isolated S. stercoralis from people living in a village with a high incidence of parasitic helminths, in particular liver flukes (Clonorchis sinensis) and hookworms, in the southern Chinese province Guangxi. We determined nuclear and mitochondrial DNA sequences of individual S. stercoralis isolated from this village and from close by hospitals and we compared these S. stercoralis among themselves and with selected published S. stercoralis from other geographic locations. For comparison, we also analyzed the hookworms present in the same location. We found that, compared to earlier studies of S. stercoralis populations in South East Asia, all S. stercoralis sampled in our study area were very closely related, suggesting a recent common source of infection for all patients. In contrast, the hookworms from the same location, while all belonging to the species Necator americanus, showed rather extensive genetic diversity even within host individuals. Different from earlier studies conducted in other geographic locations, almost all S. stercoralis in this study appeared heterozygous for different sequence variants of the 18S rDNA hypervariable regions (HVR) I and IV. In contrast to earlier investigations, except for three males, all S. stercoralis we isolated in this study were infective larvae, suggesting that the sampled population reproduces predominantly, if not exclusively through the clonal life cycles. Consistently, whole genome sequencing of individual worms revealed higher heterozygosity than reported earlier for likely sexual populations of S. stercoralis. Elevated heterozygosity is frequently associated with asexual clonal reproduction.

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.

Small peptide–mediated self-recognition prevents cannibalism in predatory nematodes

Self-recognition is observed abundantly throughout the natural world, regulating diverse biological processes. Although ubiquitous, often little is known of the associated molecular machinery, and so far, organismal self-recognition has never been described in nematodes. We investigated the predatory nematode Pristionchus pacificus and, through interactions with its prey, revealed a self-recognition mechanism acting on the nematode surface, capable of distinguishing self-progeny from closely related strains. We identified the small peptide SELF-1, which is composed of an invariant domain and a hypervariable C terminus, as a key component of self-recognition. Modifications to the hypervariable region, including single–amino acid substitutions, are sufficient to eliminate self-recognition. Thus, the P. pacificus self-recognition system enables this nematode to avoid cannibalism while promoting the killing of competing nematodes.

Common workflow language (CWL)-based software pipeline for de novo genome assembly from long- and short-read data

BACKGROUND

Here, we created an automated pipeline for the de novoassembly of genomes from Pacific Biosciences long-read and Illumina short-read data using common workflow language (CWL). To evaluate the performance of this pipeline, we assembled the nuclear genomes of the eukaryotes Caenorhabditis elegans (∼100 Mb), Drosophila melanogaster (∼138 Mb), and Plasmodium falciparum (∼23 Mb) directly from publicly accessible nucleotide sequence datasets and assessed the quality of the assemblies against curated reference genomes.

FINDINGS

We showed a dependency of the accuracy of assembly on sequencing technology and GC content and repeatedly achieved assemblies that meet the high standards set by the National Human Genome Research Institute, being applicable to gene prediction and subsequent genomic analyses.

CONCLUSIONS

This CWL pipeline overcomes current challenges of achieving repeatability and reproducibility of assembly results and offers a platform for the re-use of the workflow and the integration of diverse datasets. This workflow is publicly available via GitHub (https://github.com/vetscience/Assemblosis) and is currently applicable to the assembly of haploid and diploid genomes of eukaryotes.

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