Polyphenism of a Novel Trait Integrated Rapidly Evolving Genes into Ancestrally Plastic Networks

Pristionchus - Beetle associations: Towards a new natural history

The free-living nematode Pristionchus pacificus has been established as a model system in integrative evolutionary biology by combining laboratory studies with field work and evolutionary biology. Multiple genetic, molecular and experimental tools and a collection of more than 2,500 P. pacificus strains and more than 50 Pristionchus species, which are available as living cultures or frozen stock collections, support research on various life history traits. Species of Pristionchus exhibit a number of complex traits unknown from Caenorhabditis elegans and most other free-living nematodes. First, P. pacificus can form two alternative mouth forms, an example of developmental plasticity that is increasingly studied to investigate the role of plasticity as a facilitator of evolutionary novelty. More than a decade of work has identified associated genetic and epigenetic mechanisms and revealed the evolutionary and ecological significance of feeding structure plasticity. Second, one of the two mouth morphs results in predatory behavior against other nematodes and is currently used to investigate the neurobiology of predation. Third, potential predation results in the risk of cannibalism among conspecifics. Strikingly, Pristionchus nematodes have developed a self-recognition system that allows the distinction of self (kin) and non-self. Given all these organismal features, this nematode has recently been considered a key example for research towards a new natural history (West-Eberhard, 2024). Here, I summarize recent work on Pristionchus with a focus on a 'new natural history'. In addition, I review some recent studies that indicate an interaction of Pristionchus with EPNs that was suggested based on various surveys in different ecological habitats.

The Role of Epigenetic Switches in Polyphenism Control: Implications from a Nematode Model for the Developmental Regulation of Alternative Phenotypes

Polyphenisms, the capability of organisms to form two or more alternative phenotypes in response to environmental variation, are prevalent in nature. However, associated molecular mechanisms and potential general principles of polyphenisms among major organismal groups remain currently unknown. This review focuses on an emerging model system for developmental plasticity and polyphenism research, the nematode Pristionchus pacificus and explores mechanistic insight obtained through unbiased genetic, experimental and natural variation studies. Resulting findings identify a central role for epigenetic switches in the environmental control of alternative phenotypes and their micro-and macroevolution. Several features observed in P. pacificus are shared with insects and plants and might become general principles for the control of polyphenisms during development.

Developmental transcriptomics in Pristionchus reveals the logic of a plasticity gene regulatory network

Developmental plasticity enables the production of alternative phenotypes in response to different environmental conditions. While significant advances in understanding the ecological and evolutionary implications of plasticity have been made, understanding its genetic basis has lagged. However, a decade of genetic screens in the model nematode Pristionchus pacificus has culminated in 30 genes which affect mouth-form plasticity. We also recently reported the critical window of environmental sensitivity, and therefore have clear expectations for when differential gene expression should matter. Here, we collated previous data into a gene-regulatory network (GRN), and performed developmental transcriptomics across different environmental conditions, genetic backgrounds, and mouth-form mutants to assess the regulatory logic of plasticity. We found that only two genes in the GRN (eud-1 and seud-1/sult-1) are sensitive to the environment during the critical window. Interestingly, the time points of their sensitivity differ, suggesting that they act as sequential checkpoints. We also observed temporal constraint upon the transcriptional effects of mutating the GRN and revealed unexpected feedback between mouth-form genes. Surprisingly, expression of seud-1/sult-1, but not eud-1, correlated with mouth form biases across different strains and species. Finally, a comprehensive analysis of all samples identified metabolism as a shared pathway for regulating mouth-form plasticity. These data are presented in a Shiny app to facilitate gene-expression comparisons across development in up to 14 different conditions. Collectively, our results suggest that mouth-form plasticity evolved a constrained, two-tiered logic to integrate environmental information leading up to the final developmental decision.

Feeding-structure morphogenesis in "rhabditid" and diplogastrid nematodes is not controlled by a conserved genetic module

Disentangling the evolution of the molecular processes and genetic networks that facilitate the emergence of morphological novelties is one of the main objectives in evolutionary developmental biology. Here, we investigated the evolutionary history of a gene regulatory network controlling the development of novel tooth-like feeding structures in diplogastrid nematodes. Focusing on NHR-1 and NHR-40, the two transcription factors that regulate the morphogenesis of these feeding structures in Pristionchus pacificus, we sought to determine whether they have a similar function in Caenorhabditis elegans, an outgroup species to the Diplogastridae which has typical "rhabditid" flaps instead of teeth. Contrary to our initial expectations, we found that they do not have a similar function. While both receptors are co-expressed in the tissues that produce the feeding structures in the two nematodes, genetic inactivation of either receptor had no impact on feeding-structure morphogenesis in C. elegans. Transcriptomic experiments revealed that NHR-1 and NHR-40 have highly species-specific regulatory targets. These results suggest two possible evolutionary scenarios: either the genetic module responsible for feeding-structure morphogenesis in Diplogastridae already existed in the last common ancestor of C. elegans and P. pacificus, and subsequently disintegrated in the former as NHR-1 and NHR-40 acquired new targets, or it evolved in conjunction with teeth in Diplogastridae. These findings indicate that feeding-structure morphogenesis is regulated by different genetic programs in P. pacificus and C. elegans, hinting at developmental systems drift during the flap-to-tooth transformation. Further research in other "rhabditid" species is needed to fully reconstruct the developmental genetic changes which facilitated the evolution of novel feeding structures in Diplogastridae.

A histone demethylase links the loss of plasticity to nongenetic inheritance and morphological change

Plasticity is a widespread feature of development, enabling phenotypic change based on the environment. Although the evolutionary loss of plasticity has been linked both theoretically and empirically to increased rates of phenotypic diversification, molecular insights into how this process might unfold are generally lacking. Here, we show that a regulator of nongenetic inheritance links evolutionary loss of plasticity in nature to changes in plasticity and morphology as selected in the laboratory. Across nematodes of Diplogastridae, which ancestrally had a polyphenism, or discrete plasticity, in their feeding morphology, we use molecular evolutionary analyses to screen for change associated with independent losses of plasticity. Having inferred a set of ancestrally polyphenism-biased genes from phylogenetically informed gene-knockouts and gene-expression comparisons, selection signatures associated with plasticity's loss identify the histone H3K4 di/monodemethylase gene spr-5/LSD1/KDM1A. Manipulations of this gene affect both sensitivity and variation in plastic morphologies, and artificial selection of manipulated lines drive multigenerational shifts in these phenotypes. Our findings thus give mechanistic insight into how traits are modified as they traverse the continuum of greater to lesser environmental sensitivity.

Comparative transcriptomics reveals that a novel form of phenotypic plasticity evolved via lineage-specific changes in gene expression

Novel forms of phenotypic plasticity may evolve by lineage-specific changes or by co-opting mechanisms from more general forms of plasticity. Here, we evaluated whether a novel resource polyphenism in New World spadefoot toads (genus Spea) evolved by co-opting mechanisms from an ancestral form of plasticity common in anurans-accelerating larval development rate in response to pond drying. We compared overlap in differentially expressed genes between alternative trophic morphs constituting the polyphenism in Spea versus those found between tadpoles of Old World spadefoot toads (genus Pelobates) when experiencing different pond-drying regimes. Specifically, we (1) generated a de novo transcriptome and conducted differential gene expression analysis in Spea multiplicata, (2) utilized existing gene expression data and a recently published transcriptome for Pelobates cultripes when exposed to different drying regimes, and (3) identified unique and overlapping differentially expressed transcripts. We found thousands of differentially expressed genes between S. multiplicata morphs that were involved in major developmental reorganization, but the vast majority of these were not differentially expressed in P. cultripes. Thus, S. multiplicata's novel polyphenism appears to have arisen primarily through lineage-specific changes in gene expression and not by co-opting existing patterns of gene expression involved in pond-drying plasticity. Therefore, although ancestral stress responses might jump-start evolutionary innovation, substantial lineage-specific modification might be needed to refine these responses into more complex forms of plasticity.

A Mediator subunit imparts robustness to a polyphenism decision

Polyphenism is a type of developmental plasticity that translates continuous environmental variability into discontinuous phenotypes. Such discontinuity likely requires a switch between alternative gene-regulatory networks, a principle that has been borne out by mechanisms found to promote morph-specific gene expression. However, whether robustness is required to execute a polyphenism decision has awaited testing at the molecular level. Here, we used a nematode model for polyphenism, Pristionchus pacificus, to identify the molecular regulatory factors that ensure the development of alternative forms. This species has a dimorphism in its adult feeding structures, specifically teeth, which are a morphological novelty that allows predation on other nematodes. Through a forward genetic screen, we determined that a duplicate homolog of the Mediator subunit MDT-15/MED15, P. pacificus MDT-15.1, is necessary for the polyphenism and the robustness of the resulting phenotypes. This transcriptional coregulator, which has a conserved role in metabolic responses to nutritional stress, coordinates these processes with its effects on this diet-induced polyphenism. Moreover, this MED15 homolog genetically interacts with two nuclear receptors, NHR-1 and NHR-40, to achieve dimorphism: Single and double mutants for these three factors result in morphologies that together produce a continuum of forms between the extremes of the polyphenism. In summary, we have identified a molecular regulator that confers discontinuity to a morphological polyphenism, while also identifying a role for MED15 as a plasticity effector.

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.

The GRN concept as a guide for evolutionary developmental biology

Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.

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.

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.

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.

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.

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.

Genomic integration of transgenes using UV irradiation in Pristionchus pacificus

Transgenes are widely used throughout molecular biology for numerous applications. In Caenorhabditis elegans, stable transgenes are usually generated by microinjection into the germline establishing extrachromosomal arrays. Furthermore, numerous technologies exist to integrate transgenes into the C. elegans genome. In the nematode Pristionchus pacificus, transgenes are possible, however, their establishment is less efficient and dependent on the formation of complex arrays containing the transgene of interest and host carrier DNA. Additionally, genomic integration has only been reported via biolistic methods. Here we describe a simple technique using UV irradiation to facilitate the integration of transgenes into the P. pacificus genome.

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

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.

Linking Molecular Mechanisms and Evolutionary Consequences of Resource Polyphenism

Resource polyphenism-the occurrence of environmentally induced, discrete, and intraspecific morphs showing differential niche use-is taxonomically widespread and fundamental to the evolution of ecological function where it has arisen. Despite longstanding appreciation for the ecological and evolutionary significance of resource polyphenism, only recently have its proximate mechanisms begun to be uncovered. Polyphenism switches, especially those influencing and influenced by trophic interactions, offer a route to integrating proximate and ultimate causation in studies of plasticity, and its potential influence on evolution more generally. Here, we use the major events in generalized polyphenic development as a scaffold for linking the molecular mechanisms of polyphenic switching with potential evolutionary outcomes of polyphenism and for discussing challenges and opportunities at each step in this process. Not only does the study of resource polyphenism uncover interesting details of discrete plasticity, it also illuminates and informs general principles at the intersection of development, ecology, and evolution.

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.

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.

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.