The Orphan Gene dauerless Regulates Dauer Development and Intraspecific Competition in Nematodes by Copy Number Variation

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.

Space radiation damage rescued by inhibition of key spaceflight associated miRNAs

Our previous research revealed a key microRNA signature that is associated with spaceflight that can be used as a biomarker and to develop countermeasure treatments to mitigate the damage caused by space radiation. Here, we expand on this work to determine the biological factors rescued by the countermeasure treatment. We performed RNA-sequencing and transcriptomic analysis on 3D microvessel cell cultures exposed to simulated deep space radiation (0.5 Gy of Galactic Cosmic Radiation) with and without the antagonists to three microRNAs: miR-16-5p, miR-125b-5p, and let-7a-5p (i.e., antagomirs). Significant reduction of inflammation and DNA double strand breaks (DSBs) activity and rescue of mitochondria functions are observed after antagomir treatment. Using data from astronaut participants in the NASA Twin Study, Inspiration4, and JAXA missions, we reveal the genes and pathways implicated in the action of these antagomirs are altered in humans. Our findings indicate a countermeasure strategy that can potentially be utilized by astronauts in spaceflight missions to mitigate space radiation damage.

Nematode genome announcement: A chromosome-scale genome assembly for the Pristionchus pacificus reference mapping strain PS1843

Pristionchus pacificus is a free-living nematode that shares many features with Caenorhabditis elegans, such as its short generation time and hermaphroditism, but also exhibits novel traits, i.e., a mouth-form dimorphism that enables predation. The availability of various genetic tools and genomic resources make it a powerful model organism for comparative studies. Here, we present an updated genome of the P. pacificus strain PS1843 (Washington) that is most widely used for genetic analysis. Assembly of PacBio reads together with reference-guided scaffolding resulted in a chromosome-scale genome spanning 171Mb for the PS1843 strain. Whole genome alignments between the P. pacificus PS1843 genome and the genome of the P. pacificus reference strain PS312 (California) revealed megabase-sized regions on chromosomes III, IV, and X that explain the majority of genome size difference between both strains. The improved PS1843 genome will be useful for future forward genetic studies and evolutionary genomic comparisons at the intra-species level.

The role of plasticity and stochasticity in coexistence

Species coexistence in ecological communities is a central feature of biodiversity. Different concepts, i.e., contemporary niche theory, modern coexistence theory, and the unified neutral theory, have identified many building blocks of such ecological assemblies. However, other factors, such as phenotypic plasticity and stochastic inter-individual variation, have received little attention, in particular in animals. For example, how resource polyphenisms resulting in predator-prey interactions affect coexistence is currently unknown. Here, we present an integrative theoretical-experimental framework using the nematode plasticity model Pristionchus pacificus with its well-studied mouth-form dimorphism resulting in cannibalism. We develop an individual-based model that relies upon synthetic data based on our empirical measurements of fecundity and polyphenism to preserve demographic heterogeneity. We demonstrate how the interplay between plasticity and individual stochasticity result in all-or-nothing outcomes at the local level. Coexistence is made possible when spatial structure is introduced.

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 Lost and Found: Unraveling the Functions of Orphan Genes

Orphan Genes (OGs) are a mysterious class of genes that have recently gained significant attention. Despite lacking a clear evolutionary history, they are found in nearly all living organisms, from bacteria to humans, and they play important roles in diverse biological processes. The discovery of OGs was first made through comparative genomics followed by the identification of unique genes across different species. OGs tend to be more prevalent in species with larger genomes, such as plants and animals, and their evolutionary origins remain unclear but potentially arise from gene duplication, horizontal gene transfer (HGT), or de novo origination. Although their precise function is not well understood, OGs have been implicated in crucial biological processes such as development, metabolism, and stress responses. To better understand their significance, researchers are using a variety of approaches, including transcriptomics, functional genomics, and molecular biology. This review offers a comprehensive overview of the current knowledge of OGs in all domains of life, highlighting the possible role of dark transcriptomics in their evolution. More research is needed to fully comprehend the role of OGs in biology and their impact on various biological processes.

Nematode Pheromones: Structures and Functions

Pheromones are chemical signals secreted by one individual that can affect the behaviors of other individuals within the same species. Ascaroside is an evolutionarily conserved family of nematode pheromones that play an integral role in the development, lifespan, propagation, and stress response of nematodes. Their general structure comprises the dideoxysugar ascarylose and fatty-acid-like side chains. Ascarosides can vary structurally and functionally according to the lengths of their side chains and how they are derivatized with different moieties. In this review, we mainly describe the chemical structures of ascarosides and their different effects on the development, mating, and aggregation of nematodes, as well as how they are synthesized and regulated. In addition, we discuss their influences on other species in various aspects. This review provides a reference for the functions and structures of ascarosides and enables their better application.

Dauer fate in a Caenorhabditis elegans Boolean network model

Cellular fates are determined by genes interacting across large, complex biological networks. A critical question is how to identify causal relationships spanning distinct signaling pathways and underlying organismal phenotypes. Here, we address this question by constructing a Boolean model of a well-studied developmental network and analyzing information flows through the system. Depending on environmental signals Caenorhabditis elegans develop normally to sexual maturity or enter a reproductively delayed, developmentally quiescent 'dauer' state, progressing to maturity when the environment changes. The developmental network that starts with environmental signal and ends in the dauer/no dauer fate involves genes across 4 signaling pathways including cyclic GMP, Insulin/IGF-1, TGF-β and steroid hormone synthesis. We identified three stable motifs leading to normal development, each composed of genes interacting across the Insulin/IGF-1, TGF-β and steroid hormone synthesis pathways. Three genes known to influence dauer fate, daf-2, daf-7 and hsf-1, acted as driver nodes in the system. Using causal logic analysis, we identified a five gene cyclic subgraph integrating the information flow from environmental signal to dauer fate. Perturbation analysis showed that a multifactorial insulin profile determined the stable motifs the system entered and interacted with daf-12 as the switchpoint driving the dauer/no dauer fate. Our results show that complex organismal systems can be distilled into abstract representations that permit full characterization of the causal relationships driving developmental fates. Analyzing organismal systems from this perspective of logic and function has important implications for studies examining the evolution and conservation of signaling pathways.

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.

Multidimensional competition of nematodes affects plastic traits in a beetle ecosystem

Resource competition has driven the evolution of novel polyphenisms in numerous organisms, enhancing fitness in constantly changing environmental conditions. In natural communities, the myriad interactions among diverse species are difficult to disentangle, but the multidimensional microscopic environment of a decaying insect teeming with bacteria and fighting nematodes provides pliable systems to investigate. Necromenic nematodes of the family Diplogastridae live on beetles worldwide, innocuously waiting for their hosts’ deaths to feast on the blooming bacteria. Often, more than one worm species either affiliates with the insect or joins the microbial meal; thus, competition over limited food ensues, and phenotypic plasticity provides perks for species capable of employing polyphenisms. The recently established system of cockchafer Gymnogaster bupthalma and its occasional co-infestation of Pristionchus mayeri and Acrostichus spp. has revealed that these worms will simultaneously utilize two polyphenisms to thrive in a competitive environment. While both genera maintain plastic capacities in mouth form (strictly bacterial-feeding and omnivorous predation) and developmental pathway (direct and arrested development, dauer), P. mayeri employs both when faced with competition from Acrostichus. Here, we took advantage of the malleable system and added a third competitor, model nematode Pristionchus pacificus. Intriguingly, with a third competitor, P. mayeri is quicker to exit dauer and devour available food, while Acrostichus hides in dauer, waiting for the two Pristionchus species to leave the immediate environment before resuming development. Thus, experimental manipulation of short-lived ecosystems can be used to study the roles of polyphenisms in organismal interactions and their potential significance for evolution.

Adaptation to environmental temperature in divergent clades of the nematode Pristionchus pacificus

Because of ongoing climate change, populations of organisms are being subjected to stressful temperatures more often. This is especially problematic for ectothermic organisms, which are likely to be more sensitive to changes in temperature. Therefore, we need to know if ectotherms have adapted to environmental temperature and, if so, what are the evolutionary mechanisms behind such adaptation. Here, we use the nematode Pristionchus pacificus as a case study to investigate thermal adaptation on the Indian Ocean island of La Réunion, which experiences a range of temperatures from coast to summit. We study the evolution of high-temperature tolerance by constructing a phylogenetic tree of strains collected from many different thermal niches. We show that populations of P. pacificus at low altitudes have higher fertility at warmer temperatures. Most likely, this phenotype has arisen recently and at least twice independently, consistent with parallel evolution. We also studied low-temperature tolerance and showed that populations from high altitudes have increased their fertility at cooler temperatures. Together, these data indicate that P. pacificus strains on La Réunion are subject to divergent selection, adapting to hot and cold niches at the coast and summit of the volcano. Precisely defining these thermal niches provides essential information for models that predict the impact of future climate change on these populations.

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.

Introduction to Pristionchus pacificus anatomy

Pristionchus pacificus has emerged as an important nematode species used to understand the evolution of development and behavior. While P. pacificus (Diplogasteridae) is only distantly related to Caenorhabditis elegans (Rhabditidae), both use an identical reproductive strategy, are easily reared on bacteria in Petri dishes and complete their life cycles within a few days. Over the past 25 years, several detailed light and electron microscopy studies have elucidated the anatomy of P. pacificus and have demonstrated clear homology to many cells in C. elegans. Despite this similarity, sufficient anatomical differences between C. elegans and P. pacificus have allowed the latter to be used in comparative evo-devo studies. For example, the stoma of P. pacificus contains a large dorsal tooth used during predation on other nematodes when supplementing its primarily bacterial diet. This review discusses the main anatomical features of P. pacificus with emphasis on comparison to C. elegans.

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.

Enhancers Facilitate the Birth of De Novo Genes and Gene Integration into Regulatory Networks

Regulatory networks control the spatiotemporal gene expression patterns that give rise to and define the individual cell types of multicellular organisms. In eumetazoa, distal regulatory elements called enhancers play a key role in determining the structure of such networks, particularly the wiring diagram of “who regulates whom.” Mutations that affect enhancer activity can therefore rewire regulatory networks, potentially causing adaptive changes in gene expression. Here, we use whole-tissue and single-cell transcriptomic and chromatin accessibility data from mouse to show that enhancers play an additional role in the evolution of regulatory networks: They facilitate network growth by creating transcriptionally active regions of open chromatin that are conducive to de novo gene evolution. Specifically, our comparative transcriptomic analysis with three other mammalian species shows that young, mouse-specific intergenic open reading frames are preferentially located near enhancers, whereas older open reading frames are not. Mouse-specific intergenic open reading frames that are proximal to enhancers are more highly and stably transcribed than those that are not proximal to enhancers or promoters, and they are transcribed in a limited diversity of cellular contexts. Furthermore, we report several instances of mouse-specific intergenic open reading frames proximal to promoters showing evidence of being repurposed enhancers. We also show that open reading frames gradually acquire interactions with enhancers over macroevolutionary timescales, helping integrate genes—those that have arisen de novo or by other means—into existing regulatory networks. Taken together, our results highlight a dual role of enhancers in expanding and rewiring gene regulatory networks.

Nematode biphasic 'boom and bust' dynamics are dependent on host bacterial load while linking dauer and mouth-form polyphenisms

Cross-kingdom interactions involve dynamic processes that shape terrestrial ecosystems and represent striking examples of co-evolution. The multifaceted relationships of entomopathogenic nematodes with their insect hosts and symbiotic bacteria are well-studied cases of co-evolution and pathogenicity. In contrast, microbial interactions in soil after the natural death of insects and other invertebrates are minimally understood. In particular, the turnover and succession of nematodes and bacteria during insect decay have not been well documented - although it represents a rich ecological niche with multiple species interactions. Here, we utilize developmentally plastic nematode Pristionchus pacificus and its associated scarab beetles as models. On La Réunion Island, we collected rhinoceros beetle Oryctes borbonicus, induced death, and placed carcasses in cages both on the island and in a mock-natural environment in the laboratory controlling for high spatial and temporal resolution. Investigating nematode population density and dispersal dynamics, we were able to connect two imperative plasticities, dauer and mouth form. We observed a biphasic 'boom and bust' dispersal dynamic of dauer larvae that corresponds to bacterial load on carcasses but not bacterial type. Strikingly, all post-dauer adults have the predatory mouth form, demonstrating novel intricate interactions on decaying insect hosts. Thus, ecologically relevant survival strategies incorporate critical plastic traits.

Spatial Transcriptomics of Nematodes Identifies Sperm Cells as a Source of Genomic Novelty and Rapid Evolution

Divergence of gene function and expression during development can give rise to phenotypic differences at the level of cells, tissues, organs, and ultimately whole organisms. To gain insights into the evolution of gene expression and novel genes at spatial resolution, we compared the spatially resolved transcriptomes of two distantly related nematodes, Caenorhabditis elegans and Pristionchus pacificus, that diverged 60-90 Ma. The spatial transcriptomes of adult worms show little evidence for strong conservation at the level of single genes. Instead, regional expression is largely driven by recent duplication and emergence of novel genes. Estimation of gene ages across anatomical structures revealed an enrichment of novel genes in sperm-related regions. This provides first evidence in nematodes for the "out of testis" hypothesis that has been previously postulated based on studies in Drosophila and mammals. "Out of testis" genes represent a mix of products of pervasive transcription as well as fast evolving members of ancient gene families. Strikingly, numerous novel genes have known functions during meiosis in Caenorhabditis elegans indicating that even universal processes such as meiosis may be targets of rapid evolution. Our study highlights the importance of novel genes in generating phenotypic diversity and explicitly characterizes gene origination in sperm-related regions. Furthermore, it proposes new functions for previously uncharacterized genes and establishes the spatial transcriptome of Pristionchus pacificus as a catalog for future studies on the evolution of gene expression and function.

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.

Screening for CRISPR/Cas9-induced mutations using a co-injection marker in the nematode Pristionchus pacificus

CRISPR/Cas9 genome-editing methods are used to reveal functions of genes and molecular mechanisms underlying biological processes in many species, including nematodes. In evolutionary biology, the nematode Pristionchus pacificus is a satellite model and has been used to understand interesting phenomena such as phenotypic plasticity and self-recognition. In P. pacificus, CRISPR/Cas9-mediated mutations are induced by microinjecting a guide RNA (gRNA) and Cas9 protein into the gonads. However, mutant screening is laborious and time-consuming due to the absence of visual markers. In this study, we established a Co-CRISPR strategy by using a dominant roller marker in P. pacificus. We found that heterozygous mutations in Ppa-prl-1 induced the roller phenotype, which can be used as an injection marker. After the co-injection of Ppa-prl-1 gRNA, target gRNA, and the Cas9 protein, roller progeny and their siblings were examined using the heteroduplex mobility assay and DNA sequencing. We found that some of the roller and non-roller siblings had mutations at the target site. We used varying Cas9 concentrations and found that a higher concentration of Cas9 did not increase genome-editing events. The Co-CRISPR strategy promotes the screening for genome-editing events and will facilitate the development of new genome-editing methods in P. pacificus.

Convergent evolution of small molecule pheromones in Pristionchus nematodes

The small molecules that mediate chemical communication between nematodes-so-called 'nematode-derived-modular-metabolites' (NDMMs)-are of major interest because of their ability to regulate development, behavior, and life-history. Pristionchus pacificus nematodes produce an impressive diversity of structurally complex NDMMs, some of which act as primer pheromones that are capable of triggering irreversible developmental switches. Many of these NDMMs have only ever been found in P. pacificus but no attempts have been made to study their evolution by profiling closely related species. This study brings a comparative perspective to the biochemical study of NDMMs through the systematic MS/MS- and NMR-based analysis of exo-metabolomes from over 30 Pristionchus species. We identified 36 novel compounds and found evidence for the convergent evolution of complex NDMMs in separate branches of the Pristionchus phylogeny. Our results demonstrate that biochemical innovation is a recurrent process in Pristionchus nematodes, a pattern that is probably typical across the animal kingdom.

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.

The Enigmatic Canal-Associated Neurons Regulate Caenorhabditis elegans Larval Development Through a cAMP Signaling Pathway

Caenorhabditis elegans larval development requires the function of the two Canal-Associated Neurons (CANs): killing the CANs by laser microsurgery or disrupting their development by mutating the gene ceh-10 results in early larval arrest. How these cells promote larval development, however, remains a mystery. In screens for mutations that bypass CAN function, we identified the gene kin-29, which encodes a member of the Salt-Inducible Kinase (SIK) family and a component of a conserved pathway that regulates various C. elegans phenotypes. Like kin-29 loss, gain-of-function mutations in genes that may act upstream of kin-29 or growth in cyclic-AMP analogs bypassed ceh-10 larval arrest, suggesting that a conserved adenylyl cyclase/PKA pathway inhibits KIN-29 to promote larval development, and that loss of CAN function results in dysregulation of KIN-29 and larval arrest. The adenylyl cyclase ACY-2 mediates CAN-dependent larval development: acy-2 mutant larvae arrested development with a similar phenotype to ceh-10 mutants, and the arrest phenotype was suppressed by mutations in kin-29 ACY-2 is expressed predominantly in the CANs, and we provide evidence that the acy-2 functions in the CANs to promote larval development. By contrast, cell-specific expression experiments suggest that kin-29 acts in both the hypodermis and neurons, but not in the CANs. Based on our findings, we propose two models for how ACY-2 activity in the CANs regulates KIN-29 in target cells.

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.

Young genes have distinct gene structure, epigenetic profiles, and transcriptional regulation

Species-specific, new, or "orphan" genes account for 10%-30% of eukaryotic genomes. Although initially considered to have limited function, an increasing number of orphan genes have been shown to provide important phenotypic innovation. How new genes acquire regulatory sequences for proper temporal and spatial expression is unknown. Orphan gene regulation may rely in part on origination in open chromatin adjacent to preexisting promoters, although this has not yet been assessed by genome-wide analysis of chromatin states. Here, we combine taxon-rich nematode phylogenies with Iso-Seq, RNA-seq, ChIP-seq, and ATAC-seq to identify the gene structure and epigenetic signature of orphan genes in the satellite model nematode Pristionchus pacificus Consistent with previous findings, we find young genes are shorter, contain fewer exons, and are on average less strongly expressed than older genes. However, the subset of orphan genes that are expressed exhibit distinct chromatin states from similarly expressed conserved genes. Orphan gene transcription is determined by a lack of repressive histone modifications, confirming long-held hypotheses that open chromatin is important for new gene formation. Yet orphan gene start sites more closely resemble enhancers defined by H3K4me1, H3K27ac, and ATAC-seq peaks, in contrast to conserved genes that exhibit traditional promoters defined by H3K4me3 and H3K27ac. Although the majority of orphan genes are located on chromosome arms that contain high recombination rates and repressive histone marks, strongly expressed orphan genes are more randomly distributed. Our results support a model of new gene origination by rare integration into open chromatin near enhancers.

Deep taxon sampling reveals the evolutionary dynamics of novel gene families in Pristionchus nematodes

The widespread identification of genes without detectable homology in related taxa is a hallmark of genome sequencing projects in animals, together with the abundance of gene duplications. Such genes have been called novel, young, taxon-restricted, or orphans, but little is known about the mechanisms accounting for their origin, age, and mode of evolution. Phylogenomic studies relying on deep and systematic taxon sampling and using the comparative method can provide insight into the evolutionary dynamics acting on novel genes. We used a phylogenomic approach for the nematode model organism Pristionchus pacificus and sequenced six additional Pristionchus and two outgroup species. This resulted in 10 genomes with a ladder-like phylogeny, sequenced in one laboratory using the same platform and analyzed by the same bioinformatic procedures. Our analysis revealed that 68%-81% of genes are assignable to orthologous gene families, the majority of which defined nine age classes with presence/absence patterns that can be explained by single evolutionary events. Contrasting different age classes, we find that older age classes are concentrated at chromosome centers, whereas novel gene families preferentially arise at the periphery, are weakly expressed, evolve rapidly, and have a high propensity of being lost. Over time, they increase in expression and become more constrained. Thus, the detailed phylogenetic resolution allowed a comprehensive characterization of the evolutionary dynamics of Pristionchus genomes indicating that distribution of age classes and their associated differences shape chromosomal divergence. This study establishes the Pristionchus system for future research on the mechanisms that drive the formation of novel genes.

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

The nematode Pristionchus pacificus is an established model for integrative evolutionary biology and comparative studies with Caenorhabditis elegans. While an existing genome draft facilitated the identification of several genes controlling various developmental processes, its high degree of fragmentation complicated virtually all genomic analyses. Here, we present a de novo genome assembly from single-molecule, long-read sequencing data consisting of 135 P. pacificus contigs. When combined with a genetic linkage map, 99% of the assembly could be ordered and oriented into six chromosomes. This allowed us to robustly characterize chromosomal patterns of gene density, repeat content, nucleotide diversity, linkage disequilibrium, and macrosynteny in P. pacificus. Despite widespread conservation of synteny between P. pacificus and C. elegans, we identified one major translocation from an autosome to the sex chromosome in the lineage leading to C. elegans. This highlights the potential of the chromosome-scale assembly for future genomic studies of P. pacificus.

Comparative Genomics of Gene Loss and Gain in Caenorhabditis and Other Nematodes

Nematodes, such as Caenorhabditis elegans, form one of the most species-rich animal phyla. By now more than 30 nematode genomes have been published allowing for comparative genomic analyses at various different time-scales. The majority of a nematode's gene repertoire is represented by either duplicated or so-called orphan genes of unknown origin. This indicates the importance of mechanisms that generate new genes during the course of evolution. While it is certain that nematodes have acquired genes by horizontal gene transfer from various donors, this process only explains a small portion of the nematode gene content. As evolutionary genomic analyses strongly support that most orphan genes are indeed protein-coding, future studies will have to decide, whether they are result from extreme divergence or evolved de novo from previously noncoding sequences. In this contribution, I summarize several studies investigating gene loss and gain in nematodes and discuss the strengths and weaknesses of individual approaches and datasets. These approaches can be used to ask nematode-specific questions such as associated with the evolution of parasitism or with switches in mating systems, but also can complement studies in other animal phyla like vertebrates and insects to broaden our general view on genome evolution.

The Natural Biotic Environment of Caenorhabditis elegans

Organisms evolve in response to their natural environment. Consideration of natural ecological parameters are thus of key importance for our understanding of an organism's biology. Curiously, the natural ecology of the model species Caenorhabditis elegans has long been neglected, even though this nematode has become one of the most intensively studied models in biological research. This lack of interest changed ∼10 yr ago. Since then, an increasing number of studies have focused on the nematode's natural ecology. Yet many unknowns still remain. Here, we provide an overview of the currently available information on the natural environment of C. elegans We focus on the biotic environment, which is usually less predictable and thus can create high selective constraints that are likely to have had a strong impact on C. elegans evolution. This nematode is particularly abundant in microbe-rich environments, especially rotting plant matter such as decomposing fruits and stems. In this environment, it is part of a complex interaction network, which is particularly shaped by a species-rich microbial community. These microbes can be food, part of a beneficial gut microbiome, parasites and pathogens, and possibly competitors. C. elegans is additionally confronted with predators; it interacts with vector organisms that facilitate dispersal to new habitats, and also with competitors for similar food environments, including competitors from congeneric and also the same species. Full appreciation of this nematode's biology warrants further exploration of its natural environment and subsequent integration of this information into the well-established laboratory-based research approaches.

Oxygen-induced social behaviours in Pristionchus pacificus have a distinct evolutionary history and genetic regulation from Caenorhabditis elegans

Wild isolates of the nematode Caenorhabditis elegans perform social behaviours, namely clumping and bordering, to avoid hyperoxia under laboratory conditions. In contrast, the laboratory reference strain N2 has acquired a solitary behaviour in the laboratory, related to a gain-of-function variant in the neuropeptide Y-like receptor NPR-1. Here, we study the evolution and natural variation of clumping and bordering behaviours in Pristionchus pacificus nematodes in a natural context, using strains collected from 22 to 2400 metres above sea level on La Réunion Island. Through the analysis of 106 wild isolates, we show that the majority of strains display a solitary behaviour similar to C. elegans N2, whereas social behaviours are predominantly seen in strains that inhabit high-altitude locations. We show experimentally that P. pacificus social strains perform clumping and bordering to avoid hyperoxic conditions in the laboratory, suggesting that social strains may have adapted to or evolved a preference for the lower relative oxygen levels available at high altitude in nature. In contrast to C. elegans, clumping and bordering in P. pacificus do not correlate with locomotive behaviours in response to changes in oxygen conditions. Furthermore, QTL analysis indicates clumping and bordering to represent complex quantitative traits. Thus, clumping and bordering behaviours represent an example of phenotypic convergence with a different evolutionary history and distinct genetic control in both nematode species.

Rapid Functional and Sequence Differentiation of a Tandemly Repeated Species-Specific Multigene Family in Drosophila

Gene clusters of recently duplicated genes are hotbeds for evolutionary change. However, our understanding of how mutational mechanisms and evolutionary forces shape the structural and functional evolution of these clusters is hindered by the high sequence identity among the copies, which typically results in their inaccurate representation in genome assemblies. The presumed testis-specific, chimeric gene Sdic originated, and tandemly expanded in Drosophila melanogaster, contributing to increased male-male competition. Using various types of massively parallel sequencing data, we studied the organization, sequence evolution, and functional attributes of the different Sdic copies. By leveraging long-read sequencing data, we uncovered both copy number and order differences from the currently accepted annotation for the Sdic region. Despite evidence for pervasive gene conversion affecting the Sdic copies, we also detected signatures of two episodes of diversifying selection, which have contributed to the evolution of a variety of C-termini and miRNA binding site compositions. Expression analyses involving RNA-seq datasets from 59 different biological conditions revealed distinctive expression breadths among the copies, with three copies being transcribed in females, opening the possibility to a sexually antagonistic effect. Phenotypic assays using Sdic knock-out strains indicated that should this antagonistic effect exist, it does not compromise female fertility. Our results strongly suggest that the genome consolidation of the Sdic gene cluster is more the result of a quick exploration of different paths of molecular tinkering by different copies than a mere dosage increase, which could be a recurrent evolutionary outcome in the presence of persistent sexual selection.

First insights into the nature and evolution of antisense transcription in nematodes

BACKGROUND

The development of multicellular organisms is coordinated by various gene regulatory mechanisms that ensure correct spatio-temporal patterns of gene expression. Recently, the role of antisense transcription in gene regulation has moved into focus of research. To characterize genome-wide patterns of antisense transcription and to study their evolutionary conservation, we sequenced a strand-specific RNA-seq library of the nematode Pristionchus pacificus.

RESULTS

We identified 1112 antisense configurations of which the largest group represents 465 antisense transcripts (ASTs) that are fully embedded in introns of their host genes. We find that most ASTs show homology to protein-coding genes and are overrepresented in proteomic data. Together with the finding, that expression levels of ASTs and host genes are uncorrelated, this indicates that most ASTs in P. pacificus do not represent non-coding RNAs and do not exhibit regulatory functions on their host genes. We studied the evolution of antisense gene pairs across 20 nematode genomes, showing that the majority of pairs is lineage-specific and even the highly conserved vps-4, ddx-27, and sel-2 loci show abundant structural changes including duplications, deletions, intron gains and loss of antisense transcription. In contrast, host genes in general, are remarkably conserved and encode exceptionally long introns leading to unusually large blocks of conserved synteny.

CONCLUSIONS

Our study has shown that in P. pacificus antisense transcription as such does not define non-coding RNAs but is rather a feature of highly conserved genes with long introns. We hypothesize that the presence of regulatory elements imposes evolutionary constraint on the intron length, but simultaneously, their large size makes them a likely target for translocation of genomic elements including protein-coding genes that eventually end up as ASTs.

Phenotypic plasticity and developmental innovations in nematodes

Developmental plasticity has been implicated as a facilitator for phenotypic diversification, but the molecular mechanisms controlling it are largely unknown. We review recent comparative analyses in non-Caenorhabditis nematodes that display polyphenisms in larval development, mouth morphology and reproductive mode. Some of the challenges ahead will be to connect how these phenotypic traits are linked to each other at the molecular level, and at the ecological level. This will require sampling of several nematode species, the characterization of their ecology and the employment of both classical genetics and recently developed technological advances, such as genome editing.

Are orphan genes protein-coding, prediction artifacts, or non-coding RNAs?

Background

Current genome sequencing projects reveal substantial numbers of taxonomically restricted, so called orphan genes that lack homology with genes from other evolutionary lineages. However, it is not clear to what extent orphan genes are real, genomic artifacts, or represent non-coding RNAs.

Results

Here, we use a simple set of assumptions to test the nature of orphan genes. First, a sequence that is transcribed is considered a real biological entity. Second, every sequence that is supported by proteome data or shows a depletion of non-synonymous substitutions is a protein-coding gene. Using genomic, transcriptomic and proteomic data for the nematode Pristionchus pacificus, we show that between 4129–7997 (42–81 %) of predicted orphan genes are expressed and 3818–7545 (39–76 %) of orphan genes are under negative selection. In three cases that exhibited strong evolutionary constraint but lacked expression evidence in 14 RNA-seq samples, we could experimentally validate the predicted gene structures. Comparing different data sets to infer selection on orphan gene clusters, we find that the presence of a closely related genome provides the most powerful resource to robustly identify evidence of negative selection. However, even in the absence of other genomic data, the availability of paralogous sequences was enough to show negative selection in 8–10 % of orphan genes.

Conclusions

Our study shows that the great majority of previously identified orphan genes in P. pacificus are indeed protein-coding genes. Even though this work represents a case study on a single species, our approach can be transferred to genomic data of other non-model organisms in order to ascertain the protein-coding nature of orphan genes.

Comparative transcriptomics of the nematode gut identifies global shifts in feeding mode and pathogen susceptibility

Background

The nematode Pristionchus pacificus has been established as a model for comparative studies using the well known Caenorhabditis elegans as a reference. Despite their relatedness, previous studies have revealed highly divergent development and a number of morphological differences including the lack of a pharyngal structure, the grinder, used to physically lyse the ingested bacteria in C. elegans.

Results

To complement current knowledge about developmental and ecological differences with a better understanding of their feeding modes, we have sequenced the intestinal transcriptomes of both nematodes. In total, we found 464 intestine-enriched genes in P. pacificus and 724 in C. elegans, of which the majority (66 %) has been identified by previous studies. Interestingly, only 15 genes could be identified with shared intestinal enrichment in both species, of which three genes are Hedgehog signaling molecules supporting a highly conserved role of this pathway for intestinal development across all metazoa. At the level of gene families, we find similar divergent trends with only five families displaying significant intestinal enrichment in both species. We compared our data with transcriptomic responses to various pathogens. Strikingly, C. elegans intestine-enriched genes showed highly significant overlaps with pathogen response genes whereas this was not the case for P. pacificus, indicating shifts in pathogen susceptibility that might be explained by altered feeding modes.

Conclusions

Our study reveals first insights into the evolution of feeding systems and the associated changes in intestinal gene expression that might have facilitated nematodes of the P. pacificus lineage to colonize new environments. These findings deepen our understanding about how morphological and genomic diversity is created during the course of evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-1886-9) contains supplementary material, which is available to authorized users.

Nematode orphan genes are adopted by conserved regulatory networks and find a home in ecology

Nematode dauer formation represents an essential survival and dispersal strategy and is one of a few ecologically relevant traits that can be studied in laboratory approaches. Under harsh environmental conditions, the nematode model organisms Caenorhabditis elegans and Pristionchus pacificus arrest their development and induce the formation of stress-resistant dauer larvae in response to dauer pheromones, representing a key example of phenotypic plasticity. Previous studies have indicated that in P. pacificus, many wild isolates show cross-preference of dauer pheromones and compete for access to a limited food source. When investigating the genetic mechanisms underlying this intraspecific competition, we recently discovered that the orphan gene dauerless (dau-1) controls dauer formation by copy number variation. Our results show that dau-1 acts in parallel to or downstream of steroid hormone signaling but upstream of the nuclear hormone receptor daf-12, suggesting that DAU-1 represents a novel inhibitor of DAF-12. Phylogenetic analysis reveals that the observed copy number variation is part of a complex series of gene duplication events that occurred over short evolutionary time scales. Here, we comment on the incorporation of novel or fast-evolving genes into conserved genetic networks as a common principle for the evolution of phenotypic plasticity and intraspecific competition. We discuss the possibility that orphan genes might often function in the regulation and execution of ecologically relevant traits. Given that only few ecological processes can be studied in model organisms, the function of such genes might often go unnoticed, explaining the large number of uncharacterized genes in model system genomes.

Microevolution of Duplications and Deletions and Their Impact on Gene Expression in the Nematode Pristionchus pacificus

The evolution of diversity across the animal kingdom has been accompanied by tremendous gene loss and gain. While comparative genomics has been fruitful to characterize differences in gene content across highly diverged species, little is known about the microevolution of structural variations that cause these differences in the first place. In order to investigate the genomic impact of structural variations, we made use of genomic and transcriptomic data from the nematode Pristionchus pacificus, which has been established as a satellite model to Caenorhabditis elegans for comparative biology. We exploit the fact that P. pacificus is a highly diverse species for which various genomic data including the draft genome of a sister species P. exspectatus is available. Based on resequencing coverage data for two natural isolates we identified large (> 2kb) deletions and duplications relative to the reference strain. By restriction to completely syntenic regions between P. pacificus and P. exspectatus, we were able to polarize the comparison and to assess the impact of structural variations on expression levels. We found that while loss of genes correlates with lack of expression, duplication of genes has virtually no effect on gene expression. Further investigating expression of individual copies at sites that segregate between the duplicates, we found in the majority of cases only one of the copies to be expressed. Nevertheless, we still find that certain gene classes are strongly depleted in deletions as well as duplications, suggesting evolutionary constraint acting on synteny. In summary, our results are consistent with a model, where most structural variations are either deleterious or neutral and provide first insights into the microevolution of structural variations in the P. pacificus genome.