Toward a synthesis of developmental biology with evolutionary theory and ecology.

The evolutionary conservation of developmental mechanisms is a truism in biology, but few attempts have been made to integrate development with evolutionary theory and ecology. To work toward such a synthesis, we summarize studies in the nematode model Pristionchus pacificus, focusing on the development of the dauer, a stress-resistant, alternative larval stage. Integrative approaches combining molecular and genetic principles of development with natural variation and ecological studies in wild populations have identified a key role for a developmental switch mechanism in dauer development and evolution, one that involves the nuclear hormone receptor DAF-12. DAF-12 is a crucial regulator and convergence point for different signaling inputs, and its function is conserved among free-living and parasitic nematodes. Furthermore, DAF-12 is the target of regulatory loops that rely on novel or fast-evolving components to control the intraspecific competition of dauer larvae. We propose developmental switches as paradigms for understanding the integration of development, evolution, and ecology at the molecular level.

System-wide rewiring underlies behavioral differences in predatory and bacterial-feeding nematodes.

The relationship between neural circuit function and patterns of synaptic connectivity is poorly understood, in part due to a lack of comparative data for larger complete systems. We compare system-wide maps of synaptic connectivity generated from serial transmission electron microscopy for the pharyngeal nervous systems of two nematodes with divergent feeding behavior: the microbivore Caenorhabditis elegans and the predatory nematode Pristionchus pacificus. We uncover a massive rewiring in a complex system of identified neurons, all of which are homologous based on neurite anatomy and cell body position. Comparative graph theoretical analysis reveals a striking pattern of neuronal wiring with increased connectional complexity in the anterior pharynx correlating with tooth-like denticles, a morphological feature in the mouth of P. pacificus. We apply focused centrality methods to identify neurons I1 and I2 as candidates for regulating predatory feeding and predict substantial divergence in the function of pharyngeal glands.

A developmental switch coupled to the evolution of plasticity acts through a sulfatase.

Developmental plasticity has been suggested to facilitate phenotypic diversity, but the molecular mechanisms underlying this relationship are little understood. We analyzed a feeding dimorphism in Pristionchus nematodes whereby one of two alternative adult mouth forms is executed after an irreversible developmental decision. By integrating developmental genetics with functional tests in phenotypically divergent populations and species, we identified a regulator of plasticity, eud-1, that acts in a developmental switch. eud-1 mutations eliminate one mouth form, whereas overexpression of eud-1 fixes it. EUD-1 is a sulfatase that acts dosage dependently, is necessary and sufficient to control the sexual dimorphism of feeding forms, and has a conserved function in Pristionchus evolution. It is epistatic to known signaling cascades and results from lineage-specific gene duplications. EUD-1 thus executes a developmental switch for morphological plasticity in the adult stage, showing that regulatory pathways can evolve by terminal addition of new genes.

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

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

A safety mechanism enables tissue-specific resistance to protein aggregation during aging in C. elegans.

During aging, proteostasis capacity declines and distinct proteins become unstable and can accumulate as protein aggregates inside and outside of cells. Both in disease and during aging, proteins selectively aggregate in certain tissues and not others. Yet, tissue-specific regulation of cytoplasmic protein aggregation remains poorly understood. Surprisingly, we found that the inhibition of 3 core protein quality control systems, namely chaperones, the proteasome, and macroautophagy, leads to lower levels of age-dependent protein aggregation in Caenorhabditis elegans pharyngeal muscles, but higher levels in body-wall muscles. We describe a novel safety mechanism that selectively targets newly synthesized proteins to suppress their aggregation and associated proteotoxicity. The safety mechanism relies on macroautophagy-independent lysosomal degradation and involves several previously uncharacterized components of the intracellular pathogen response (IPR). We propose that this protective mechanism engages an anti-aggregation machinery targeting aggregating proteins for lysosomal degradation.

Extracellular proteostasis prevents aggregation during pathogenic attack.

In metazoans, the secreted proteome participates in intercellular signalling and innate immunity, and builds the extracellular matrix scaffold around cells. Compared with the relatively constant intracellular environment, conditions for proteins in the extracellular space are harsher, and low concentrations of ATP prevent the activity of intracellular components of the protein quality-control machinery. Until now, only a few bona fide extracellular chaperones and proteases have been shown to limit the aggregation of extracellular proteins1-5. Here we performed a systematic analysis of the extracellular proteostasis network in Caenorhabditis elegans with an RNA interference screen that targets genes that encode the secreted proteome. We discovered 57 regulators of extracellular protein aggregation, including several proteins related to innate immunity. Because intracellular proteostasis is upregulated in response to pathogens6-9, we investigated whether pathogens also stimulate extracellular proteostasis. Using a pore-forming toxin to mimic a pathogenic attack, we found that C. elegans responded by increasing the expression of components of extracellular proteostasis and by limiting aggregation of extracellular proteins. The activation of extracellular proteostasis was dependent on stress-activated MAP kinase signalling. Notably, the overexpression of components of extracellular proteostasis delayed ageing and rendered worms resistant to intoxication. We propose that enhanced extracellular proteostasis contributes to systemic host defence by maintaining a functional secreted proteome and avoiding proteotoxicity.

Spatial and temporal heterogeneity alter the cost of plasticity in Pristionchus pacificus.

Phenotypic plasticity, the ability of a single genotype to produce distinct phenotypes under different environmental conditions, has become a leading concept in ecology and evolutionary biology, with the most extreme examples being the formation of alternative phenotypes (polyphenisms). However, several aspects associated with phenotypic plasticity remain controversial, such as the existence of associated costs. While already predicted by some of the pioneers of plasticity research, i.e. Schmalhausen and Bradshaw, experimental and theoretical approaches have provided limited support for the costs of plasticity. In experimental studies, one common restriction is the measurement of all relevant parameters over long time periods. Similarly, theoretical studies rarely use modelling approaches that incorporate specific experimentally-derived fitness parameters. Therefore, the existence of the costs of plasticity remains disputed. Here, we provide an integrative approach to understand the cost of adaptive plasticity and its ecological ramifications, by combining laboratory data from the nematode plasticity model system Pristionchus pacificus with a stage-structured population model. Taking advantage of measurements of two isogenic strains grown on two distinct diets, we illustrate how spatial and temporal heterogeneity with regard to the distribution of resources on a metapopulation can alter the outcome of the competition and alleviate the realized cost of plasticity.

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.

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.

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.

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.

Horizontally Acquired Cellulases Assist the Expansion of Dietary Range in Pristionchus Nematodes.

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

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

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

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.

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

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

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.

New Gene Origin and Deep Taxon Phylogenomics: Opportunities and Challenges.

The immense morphological and phenotypic diversity within eukaryotes coincides with large-scale differences in genic repertoires, including the presence of thousands of new genes in every genome. New genes arise through duplication and divergence of existing coding sequences or de novo from noncoding sequences. These processes together cause individual genomes to contain up to one-third of orphan genes without any detectable homology in other lineages. Recently, deep taxon phylogenomics, the genome comparisons of extremely closely related species, provided novel insight into the evolutionary dynamics of such rapidly evolving genes. This review focuses on deep taxon phylogenomics and its importance in studying the evolution of new genes and discusses challenges and opportunities.

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.

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.

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.