Characterization of genetic diversity in the nematode Pristionchus pacificus from population-scale resequencing data

Global analysis of neuropeptide receptor conservation across phylum Nematoda

BACKGROUND

The phylum Nematoda is incredibly diverse and includes many parasites of humans, livestock, and plants. Peptide-activated G protein-coupled receptors (GPCRs) are central to the regulation of physiology and numerous behaviors, and they represent appealing pharmacological targets for parasite control. Efforts are ongoing to characterize the functions and define the ligands of nematode GPCRs, with already most peptide GPCRs known or predicted in Caenorhabditis elegans. However, comparative analyses of peptide GPCR conservation between C. elegans and other nematode species are limited, and many nematode GPCRs remain orphan. A phylum-wide perspective on peptide GPCR profiles will benefit functional and applied studies of nematode peptide GPCRs.

RESULTS

We constructed a pan-phylum resource of C. elegans peptide GPCR orthologs in 125 nematode species using a semi-automated pipeline for analysis of predicted proteome datasets. The peptide GPCR profile varies between nematode species of different phylogenetic clades and multiple C. elegans peptide GPCRs have orthologs across the phylum Nematoda. We identified peptide ligands for two highly conserved orphan receptors, NPR-9 and NPR-16, that belong to the bilaterian galanin/allatostatin A (Gal/AstA) and somatostatin/allatostatin C (SST/AstC) receptor families. The AstA-like NLP-1 peptides activate NPR-9 in cultured cells and are cognate ligands of this receptor in vivo. In addition, we discovered an AstC-type peptide, NLP-99, that activates the AstC-type receptor NPR-16. In our pan-phylum resource, the phylum-wide representation of NPR-9 and NPR-16 resembles that of their cognate ligands more than those of allatostatin-like peptides that do not activate these receptors.

CONCLUSIONS

The repertoire of C. elegans peptide GPCR orthologs varies across phylogenetic clades and several peptide GPCRs show broad conservation in the phylum Nematoda. Our work functionally characterizes the conserved receptors NPR-9 and NPR-16 as the respective GPCRs for the AstA-like NLP-1 peptides and the AstC-related peptide NLP-99. NLP-1 and NLP-99 are widely conserved in nematodes and their representation matches that of their receptor in most species. These findings demonstrate the conservation of a functional Gal/AstA and SST/AstC signaling system in nematodes. Our dataset of C. elegans peptide GPCR orthologs also lays a foundation for further functional studies of peptide GPCRs in the widely diverse nematode phylum.

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.

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.

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.

Interpreting the pervasive observation of U-shaped Site Frequency Spectra

The standard neutral model of molecular evolution has traditionally been used as the null model for population genomics. We gathered a collection of 45 genome-wide site frequency spectra from a diverse set of species, most of which display an excess of low and high frequency variants compared to the expectation of the standard neutral model, resulting in U-shaped spectra. We show that multiple merger coalescent models often provide a better fit to these observations than the standard Kingman coalescent. Hence, in many circumstances these under-utilized models may serve as the more appropriate reference for genomic analyses. We further discuss the underlying evolutionary processes that may result in the widespread U-shape of frequency spectra.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

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

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

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

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

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.

Genealogical structure changes as range expansions transition from pushed to pulled

Range expansions accelerate evolution through multiple mechanisms, including gene surfing and genetic drift. The inference and control of these evolutionary processes ultimately rely on the information contained in genealogical trees. Currently, there are two opposing views on how range expansions shape genealogies. In invasion biology, expansions are typically approximated by a series of population bottlenecks producing genealogies with only pairwise mergers between lineages-a process known as the Kingman coalescent. Conversely, traveling wave models predict a coalescent with multiple mergers, known as the Bolthausen-Sznitman coalescent. Here, we unify these two approaches and show that expansions can generate an entire spectrum of coalescent topologies. Specifically, we show that tree topology is controlled by growth dynamics at the front and exhibits large differences between pulled and pushed expansions. These differences are explained by the fluctuations in the total number of descendants left by the early founders. High growth cooperativity leads to a narrow distribution of reproductive values and the Kingman coalescent. Conversely, low growth cooperativity results in a broad distribution, whose exponent controls the merger sizes in the genealogies. These broad distribution and non-Kingman tree topologies emerge due to the fluctuations in the front shape and position and do not occur in quasi-deterministic simulations. Overall, our results show that range expansions provide a robust mechanism for generating different types of multiple mergers, which could be similar to those observed in populations with strong selection or high fecundity. Thus, caution should be exercised in making inferences about the origin of non-Kingman genealogies.

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.

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

Developmental polyphenism, the ability to switch between phenotypes in response to environmental variation, involves the alternating activation of environmentally sensitive genes. Consequently, to understand how a polyphenic response evolves requires a comparative analysis of the components that make up environmentally sensitive networks. Here, we inferred coexpression networks for a morphological polyphenism, the feeding-structure dimorphism of the nematode Pristionchus pacificus. In this species, individuals produce alternative forms of a novel trait—moveable teeth, which in one morph enable predatory feeding—in response to environmental cues. To identify the origins of polyphenism network components, we independently inferred coexpression modules for more conserved transcriptional responses, including in an ancestrally nonpolyphenic nematode species. Further, through genome-wide analyses of these components across the nematode family (Diplogastridae) in which the polyphenism arose, we reconstructed how network components have changed. To achieve this, we assembled and resolved the phylogenetic context for five genomes of species representing the breadth of Diplogastridae and a hypothesized outgroup. We found that gene networks instructing alternative forms arose from ancestral plastic responses to environment, specifically starvation-induced metabolism and the formation of a conserved diapause (dauer) stage. Moreover, loci from rapidly evolving gene families were integrated into these networks with higher connectivity than throughout the rest of the P. pacificus transcriptome. In summary, we show that the modular regulatory outputs of a polyphenic response evolved through the integration of conserved plastic responses into networks with genes of high evolutionary turnover.

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.

Four Pristionchus species associated with two mass-occurring Parafontaria laminata populations

Phoretic nematodes associated with two mass-occurring populations of the millipede Parafontaria laminata were examined, focusing on Pristionchus spp. The nematodes that propagated on dissected millipedes were genotyped using the D2-D3 expansion segments of the 28S ribosomal RNA gene. Four Pristionchus spp. were detected: P. degawai, P. laevicollis, P. fukushimae, and P. entomophagus. Of the four, P. degawai dominated and it was isolated from more than 90% of the millipedes examined. The haplotypes of partial sequences of mitochondrial cytochrome oxidase subunit I examined for Pristionchus spp. and P. degawai showed high haplotype diversity.

A model of autowave self-organization as a hierarchy of active media in the biological evolution

Within the framework of the active media concept, we develop a biophysical model of autowave self-organization which is treated as a hierarchy of active media in the evolution of the biosphere. We also propose a mathematical model of the autowave process of speciation in a flow of mutations for the three main taxonometric groups (prokaryotes, unicellular and multicellular eukaryotes) with a naturally determined lower boundary of living matter (the appearance of prokaryotes) and an open upper boundary for the formation of new species. It is shown that the fluctuation-bifurcation description of the evolution for the formation of new taxonometric groups as a trajectory of transformation of small fluctuations into giant ones adequately reflects the process of self-organization during the formation of taxa. The major concepts of biological evolution, conditions of hierarchy formation as a fundamental manifestation of self-organization and complexity in the evolution of biological systems are considered.

Genomic studies on Strongyloides stercoralis in northern and western Thailand

Background

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

Methods

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

Results

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

Conclusions

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

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

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.

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.

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

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

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

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

Comparative Epigenomics Reveals that RNA Polymerase II Pausing and Chromatin Domain Organization Control Nematode piRNA Biogenesis

Piwi-interacting RNAs (piRNAs) are important for genome regulation across metazoans, but their biogenesis evolves rapidly. In Caenorhabditis elegans, piRNA loci are clustered within two 3-Mb regions on chromosome IV. Each piRNA locus possesses an upstream motif that recruits RNA polymerase II to produce an ∼28 nt primary transcript. We used comparative epigenomics across nematodes to gain insight into the origin, evolution, and mechanism of nematode piRNA biogenesis. We show that the piRNA upstream motif is derived from core promoter elements controlling snRNA transcription. We describe two alternative modes of piRNA organization in nematodes: in C. elegans and closely related nematodes, piRNAs are clustered within repressive H3K27me3 chromatin, while in other species, typified by Pristionchus pacificus, piRNAs are found within introns of active genes. Additionally, we discover that piRNA production depends on sequence signals associated with RNA polymerase II pausing. We show that pausing signals synergize with chromatin to control piRNA transcription.

•

Nematode piRNA transcription evolved from small nuclear RNA biogenesis

•

Clustered piRNAs are produced from regulated (H3K27me3) chromatin domains

•

Dispersed piRNAs are produced from active (H3K36me3) chromatin domains

•

RNA polymerase II pausing determines the short (∼28 nt) length of piRNA precursors

A sulfotransferase dosage-dependently regulates mouthpart polyphenism in the nematode Pristionchus pacificus

Polyphenism, the extreme form of developmental plasticity, is the ability of a genotype to produce discrete morphologies matched to alternative environments. Because polyphenism is likely to be under switch-like molecular control, a comparative genetic approach could reveal the molecular targets of plasticity evolution. Here we report that the lineage-specific sulfotransferase SEUD-1, which responds to environmental cues, dosage-dependently regulates polyphenism of mouthparts in the nematode Pristionchus pacificus. SEUD-1 is expressed in cells producing dimorphic morphologies, thereby integrating an intercellular signalling mechanism at its ultimate target. Additionally, multiple alterations of seud-1 support it as a potential target for plasticity evolution. First, a recent duplication of seud-1 in a sister species reveals a direct correlation between genomic dosage and polyphenism threshold. Second, inbreeding to produce divergent polyphenism thresholds resulted in changes in transcriptional dosage of seud-1. Our study thus offers a genetic explanation for how plastic responses evolve.

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.

The genetics of phenotypic plasticity in nematode feeding structures

Phenotypic plasticity has been proposed as an ecological and evolutionary concept. Ecologically, it can help study how genes and the environment interact to produce robust phenotypes. Evolutionarily, as a facilitator it might contribute to phenotypic novelty and diversification. However, the discussion of phenotypic plasticity remains contentious in parts due to the absence of model systems and rigorous genetic studies. Here, we summarize recent work on the nematode Pristionchus pacificus, which exhibits a feeding plasticity allowing predatory or bacteriovorous feeding. We show feeding plasticity to be controlled by developmental switch genes that are themselves under epigenetic control. Phylogenetic and comparative studies support phenotypic plasticity and its role as a facilitator of morphological novelty and diversity.

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.

The genomic features of parasitism, Polyembryony and immune evasion in the endoparasitic wasp Macrocentrus cingulum

BACKGROUND

Parasitoid wasps are well-known natural enemies of major agricultural pests and arthropod borne diseases. The parasitoid wasp Macrocentrus cingulum (Hymenoptera: Braconidae) has been widely used to control the notorious insect pests Ostrinia furnacalis (Asian Corn Borer) and O. nubilalis (European corn borer). One striking phenomenon exhibited by M. cingulum is polyembryony, the formation of multiple genetically identical offspring from a single zygote. Moreover, M. cingulum employs a passive parasitic strategy by preventing the host's immune system from recognizing the embryo as a foreign body. Thus, the embryos evade the host's immune system and are not encapsulated by host hemocytes. Unfortunately, the mechanism of both polyembryony and immune evasion remains largely unknown.

RESULTS

We report the genome of the parasitoid wasp M. cingulum. Comparative genomics analysis of M. cingulum and other 11 insects were conducted, finding some gene families with apparent expansion or contraction which might be linked to the parasitic behaviors or polyembryony of M. cingulum. Moreover, we present the evidence that the microRNA miR-14b regulates the polyembryonic development of M. cingulum by targeting the c-Myc Promoter-binding Protein 1 (MBP-1), histone-lysine N-methyltransferase 2E (KMT2E) and segmentation protein Runt. In addition, Hemomucin, an O-glycosylated transmembrane protein, protects the endoparasitoid wasp larvae from being encapsulated by host hemocytes. Motif and domain analysis showed that only the hemomucin in two endoparasitoids, M. cingulum and Venturia canescens, possessing the ability of passive immune evasion has intact mucin domain and similar O-glycosylation patterns, indicating that the hemomucin is a key factor modulating the immune evasion.

CONCLUSIONS

The microRNA miR-14b participates in the regulation of polyembryonic development, and the O-glycosylation of the mucin domain in the hemomucin confers the passive immune evasion in this wasp. These key findings provide new insights into the polyembryony and immune evasion.

Linking Genomic and Metabolomic Natural Variation Uncovers Nematode Pheromone Biosynthesis

In the nematodes Caenorhabditis elegans and Pristionchus pacificus, a modular library of small molecules control behavior, lifespan, and development. However, little is known about the final steps of their biosynthesis, in which diverse building blocks from primary metabolism are attached to glycosides of the dideoxysugar ascarylose, the ascarosides. We combine metabolomic analysis of natural isolates of P. pacificus with genome-wide association mapping to identify a putative carboxylesterase, Ppa-uar-1, that is required for attachment of a pyrimidine-derived moiety in the biosynthesis of ubas#1, a major dauer pheromone component. Comparative metabolomic analysis of wild-type and Ppa-uar-1 mutants showed that Ppa-uar-1 is required specifically for the biosynthesis of ubas#1 and related metabolites. Heterologous expression of Ppa-UAR-1 in C. elegans yielded a non-endogenous ascaroside, whose structure confirmed that Ppa-uar-1 is involved in modification of a specific position in ascarosides. Our study demonstrates the utility of natural variation-based approaches for uncovering biosynthetic pathways.

Full mitochondrial and nuclear genome comparison confirms that Onchocerca sp. "Siisa" is Onchocerca ochengi

Onchocerca ochengi is a nodule-forming filarial nematode parasite of cattle. It is the closest known relative of the human parasite Onchocerca volvulus, with which it shares the black fly vector Simulium damnosum. Onchocerca sp. "Siisa" was described in black flies and in cattle and, based on limited mitochondrial sequence information, appeared to be about equally phylogenetically distant from O. ochengi and O. volvulus. Based on molecular genetic markers and apparent interbreeding, we later proposed that O. sp. "Siisa" belongs to the species O. ochengi. However, we did not demonstrate directly that the hybrids were fertile, and we were still unable to resolve the phylogenetic relationship of O. ochengi, O. sp. "Siisa," and O. volvulus, leaving some concerns with the conclusion mentioned above. Here, we present fully assembled, manually curated mitochondrial genomes of O. ochengi and O. sp. "Siisa," and we compare multiple individuals of these two taxa with respect to their whole mitochondrial and nuclear genomes. Based on the mitochondrial genomes, O. ochengi and O. sp. "Siisa" are phylogenetically much closer to each other than to O. volvulus. The differences between them are well within the range of what is expected for within-species variation. The nuclear genome comparison provided no indication of genetic separation of O. ochengi and O. sp. "Siisa." From this, in combination with the earlier literature, we conclude that O. ochengi and O. sp. "Siisa" should be considered one species.

RNA-Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates

Conflict between organisms can lead to a reciprocal adaptation that manifests as an increased evolutionary rate in genes mediating the conflict. This adaptive signature has been observed in RNA-interference (RNAi) pathway genes involved in the suppression of viruses and transposable elements in Drosophila melanogaster, suggesting that a subset of Drosophila RNAi genes may be locked in an arms race with these parasites. However, it is not known whether rapid evolution of RNAi genes is a general phenomenon across invertebrates, or which RNAi genes generally evolve adaptively. Here we use population genomic data from eight invertebrate species to infer rates of adaptive sequence evolution, and to test for past and ongoing selective sweeps in RNAi genes. We assess rates of adaptive protein evolution across species using a formal meta-analytic framework to combine data across species and by implementing a multispecies generalized linear mixed model of mutation counts. Across species, we find that RNAi genes display a greater rate of adaptive protein substitution than other genes, and that this is primarily mediated by positive selection acting on the genes most likely to defend against viruses and transposable elements. In contrast, evidence for recent selective sweeps is broadly spread across functional classes of RNAi genes and differs substantially among species. Finally, we identify genes that exhibit elevated adaptive evolution across the analyzed insect species, perhaps due to concurrent parasite-mediated arms races.

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.

Regulation of hyperoxia-induced social behaviour in Pristionchus pacificus nematodes requires a novel cilia-mediated environmental input

Social behaviours are frequently utilised for defence and stress avoidance in nature. Both Caenorhabditis elegans and Pristionchus pacificus nematodes display social behaviours including clumping and bordering, to avoid hyperoxic stress conditions. Additionally, both species show natural variation in social behaviours with “social” and “solitary” strains. While the single solitary C. elegans N2 strain has evolved under laboratory domestication due to a gain-of-function mutation in the neuropeptide receptor gene npr-1, P. pacificus solitary strains are commonplace and likely ancestral. P. pacificus therefore provides an opportunity to further our understanding of the mechanisms regulating these complex behaviours and how they evolved within an ecologically relevant system. Using CRISPR/Cas9 engineering, we show that Ppa-npr-1 has minimal influence on social behaviours, indicating independent evolutionary pathways compared to C. elegans. Furthermore, solitary P. pacificus strains show an unexpected locomotive response to hyperoxic conditions, suggesting a novel regulatory mechanism counteracting social behaviours. By utilising both forward and reverse genetic approaches we identified 10 genes of the intraflagellar transport machinery in ciliated neurons that are essential for this inhibition. Therefore, a novel cilia-mediated environmental input adds an additional level of complexity to the regulation of hyperoxia-induced social behaviours in P. pacificus, a mechanism unknown in C. elegans.

Duplications and Positive Selection Drive the Evolution of Parasitism-Associated Gene Families in the Nematode Strongyloides papillosus

Gene duplication is a major mechanism playing a role in the evolution of phenotypic complexity and in the generation of novel traits. By comparing parasitic and nonparasitic nematodes, a recent study found that the evolution of parasitism in Strongyloididae is associated with a large expansion in the Astacin and CAP gene families.To gain novel insights into the developmental processes in the sheep parasite Strongyloides papillosus, we sequenced transcriptomes of different developmental stages and sexes. Overall, we found that the majority of genes are developmentally regulated and have one-to-one orthologs in the diverged S. ratti genome. Together with the finding of similar expression profiles between S. papillosus and S. ratti, these results indicate a strong evolutionary constraint acting against change at sequence and expression levels. However, the comparison between parasitic and free-living females demonstrates a quite divergent pattern that is mostly due to the previously mentioned expansion in the Astacin and CAP gene families. More detailed phylogenetic analysis of both gene families shows that most members date back to single expansion events early in the Strongyloides lineage and have undergone subfunctionalization resulting in clusters that are highly expressed either in infective larvae or in parasitic females. Finally, we found increased evidence for positive selection in both gene families relative to the genome-wide expectation.In summary, our study reveals first insights into the developmental transcriptomes of S. papillosus and provides a detailed analysis of sequence and expression evolution in parasitism-associated gene families.

A species-specific nematocide that results in terminal embryogenesis

Nematode-insect interactions are ubiquitous, complex and constantly changing as the host and nematode coevolve. The entomophilic nematode Pristionchus pacificus is found on a myriad beetle species worldwide, although the molecular dynamics of this relationship are largely unknown. To better understand how host cues affect P. pacificus embryogenesis, we characterized the threshold of sensitivity to the pheromone (Z)-7-tetradecen-2-one (ZTDO) by determining the minimum exposure duration and developmental window that results in P. pacificus embryonic lethality. We found early-stage embryos exposed to volatile ZTDO for as few as 4 h all display terminal embryogenesis, characterized by punctuated development up to 48 h later, with abnormal morphology and limited cavity formation. To determine if the pheromone arrests pre-hatching development by suffocating or permeabilizing the eggshells, we raised embryos under anoxic conditions and also examined eggshell permeability using the lipophilic dye FM4-64. We found that asphyxiating the embryos arrested embryogenesis in a reversible manner but did not phenocopy the effects of ZTDO exposure, whereas the ZTDO-induced disruption of embryogenesis did correlate with increased eggshell permeability. The effects of ZTDO are also highly specific, as other lipid insect compounds do not produce any detectable embryocidal effect. The high specificity and unusual teratogenic effect of ZTDO may be important in mediating the host-nematode relationship by regulating P. pacificus development.

The Role of DAF-21/Hsp90 in Mouth-Form Plasticity in Pristionchus pacificus

Phenotypic plasticity is increasingly recognized to facilitate adaptive change in plants and animals, including insects, nematodes, and vertebrates. Plasticity can occur as continuous or discrete (polyphenisms) variation. In social insects, for example, in ants, some species have workers of distinct size classes while in other closely related species variation in size may be continuous. Despite the abundance of examples in nature, how discrete morphs are specified remains currently unknown. In theory, polyphenisms might require robustness, whereby the distribution of morphologies would be limited by the same mechanisms that execute buffering from stochastic perturbations, a function attributed to heat-shock proteins of the Hsp90 family. However, this possibility has never been directly tested because plasticity and robustness are considered to represent opposite evolutionary principles. Here, we used a polyphenism of feeding structures in the nematode Pristionchus pacificus to test the relationship between robustness and plasticity using geometric morphometrics of 20 mouth-form landmarks. We show that reducing heat-shock protein activity, which reduces developmental robustness, increases the range of mouth-form morphologies. Specifically, elevated temperature led to a shift within morphospace, pharmacological inhibition of all Hsp90 genes using radicicol treatment increased shape variability in both mouth-forms, and CRISPR/Cas9-induced Ppa-daf-21/Hsp90 knockout had a combined effect. Thus, Hsp90 canalizes the morphologies of plastic traits resulting in discrete polyphenism of mouth-forms.

Unexpected sex-specific post-reproductive lifespan in the free-living nematode Pristionchus exspectatus

Patterns of senescence (or aging) can vary among life history traits and between the sexes, providing an opportunity to study variation in the aging process within a single species. We previously found that females of the nematode Pristionchus exspectatus outlive males by a substantial margin under laboratory conditions. Here, we show that sex-specific reproductive senescence unfolds in the opposite direction in this species, resulting in a prolonged period of female-specific post-reproductive survival: females lost the ability to reproduce at approximately 4.7 weeks despite a median lifespan of about 12.3 weeks under lab conditions, whereas males lost the ability to reproduce at approximately 6.6 weeks, roughly in line with their median lifespan of around 7.6 weeks. Interestingly, somatic senescence (declining crawling speed) only explained reproductive senescence in males, whereas females lost the ability to reproduce regardless of condition. However, we found that housing females with males significantly increased their mortality rate, indicating that female-specific post-reproductive survival is unlikely to occur in the wild. We discuss our results in light of evolutionary theories of post-reproductive survival and previous studies of nematode behavioral ecology, arguing that premature reproductive senescence may stem from sex-specific condition-dependent survival during the reproductive period. Given the proven lab tractability of Prisitonchus nematodes, our findings provide a foundation for integrative research that combines evolutionary ecology and molecular genetics in the study of sex-specific senescence and post-reproductive survival.

Variation in rates of spontaneous male production within the nematode species Pristionchus pacificus supports an adaptive role for males and outcrossing

Background

The nematode species Pristionchus pacificus has an androdioecious mating system in which populations consist of self-fertilizing hermaphrodites and relatively few males. The prevalence of males in such a system is likely to depend on the relative pros and cons of outcrossing. While outcrossing generates novel allelic combinations and can therefore increase adaptive potential, it may also disrupt the potentially beneficial consequences of repeated generations of selfing. These include purging of deleterious alleles, inheritance of co-adapted allele complexes, improved hermaphrodite fitness and increased population growth. Here we use experimental and population genetic approaches to test hypotheses relating to male production and outcrossing in laboratory and natural populations of P. pacificus sampled from the volcanic island of La Réunion.

Results

We find a significant interaction between sampling locality and temperature treatment influencing rates of spontaneous male production in the laboratory. While strains isolated at higher altitude, cooler localities produce a higher proportion of male offspring at 25 °C relative to 20 or 15 °C, the reverse pattern is seen in strains isolated from warmer, low altitude localities. Linkage disequilibrium extends across long physical distances, but fails to approach levels reported for the partially selfing nematode species Caenorhabditis elegans. Finally, we find evidence for admixture between divergent genetic lineages.

Conclusions

Elevated rates of laboratory male generation appear to occur under environmental conditions which differ from those experienced by populations in nature. Such elevated male generation may result in higher outcrossing rates, hence driving increased effective recombination and the creation of potentially adaptive novel allelic combinations. Patterns of linkage disequilibrium decay support selfing as the predominant reproductive strategy in P. pacificus. Finally, despite the potential for outcrossing depression, our results suggest admixture has occurred between distinct genetic lineages since their independent colonization of the island, suggesting outcrossing depression may not be uniform in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0873-7) contains supplementary material, which is available to authorized users.

Genome-Wide Analyses of Individual Strongyloides stercoralis (Nematoda: Rhabditoidea) Provide Insights into Population Structure and Reproductive Life Cycles

The helminth Strongyloides stercoralis, which is transmitted through soil, infects 30–100 million people worldwide. S. stercoralis reproduces sexually outside the host as well as asexually within the host, which causes a life-long infection. To understand the population structure and transmission patterns of this parasite, we re-sequenced the genomes of 33 individual S. stercoralis nematodes collected in Myanmar (prevalent region) and Japan (non-prevalent region). We utilised a method combining whole genome amplification and next-generation sequencing techniques to detect 298,202 variant positions (0.6% of the genome) compared with the reference genome. Phylogenetic analyses of SNP data revealed an unambiguous geographical separation and sub-populations that correlated with the host geographical origin, particularly for the Myanmar samples. The relatively higher heterozygosity in the genomes of the Japanese samples can possibly be explained by the independent evolution of two haplotypes of diploid genomes through asexual reproduction during the auto-infection cycle, suggesting that analysing heterozygosity is useful and necessary to infer infection history and geographical prevalence.

Strongyloides stercoralis, one of the most neglected helminths causes strongyloidiasis mainly in tropical and subtropical regions worldwide. The parasite’s complex lifecycle includes sexual and asexual reproduction outside and inside the host, respectively. The parasite can also asexually complete a life cycle within the host's body, which is called autoinfection causing life-long infection. In order to investigate the population structure and transmission patterns of this parasite we sequenced individual nematodes isolated from human faeces in Japan and Myanmar, where the parasite is present at low and high frequencies, respectively. Whole genome sequencing of small parasites is generally difficult because the amount of DNA is limiting. However, we overcame this problem by combining whole genome amplification with next-generation sequencing. Sequence comparisons revealed 0.6% of the genome is variable among samples, and the variants showed clear separation by the location of their origin. We found that heterozygosity within the genomes was higher in Japan, which is likely explained by the predominance of asexual reproduction through auto-infection, suggesting that analyses of heterozygosity are required to better understand the history of a population.

Hybrid crosses and the genetic basis of interspecific divergence in lifespan in Pristionchus nematodes

Characterizing the genetic basis of among-species variation in lifespan is a major goal of evolutionary gerontology research, but the very feature that defines separate species - the inability to interbreed - makes achieving this goal impractical, if not impossible, for most taxa. Pristionchus nematodes provide an intriguing system for tackling this problem, as female lifespan varies among species that can be crossed to form viable (although infertile) hybrids. By conducting reciprocal crosses among three species - two dioecious (long-lived Pristionchus exspectatus and short-lived Pristionchus arcanus) and one androdioecious (short-lived Pristionchus pacificus) - we found that female lifespan was long for all hybrids, consistent with the hypothesis that the relatively short lifespans seen for P. pacificus hermaphrodites and P. arcanus females are caused by independent, recessive alleles that are masked in hybrid genomes. Cross-direction had a small effect on survivorship for crosses involving P. exspectatus, indicating that nuclear-mitochondrial interactions may also influence Pristionchus longevity. Our findings suggest that long lifespan in P. exspectatus reflects the realization of an ancestral potential for extended longevity in the P. pacificus species complex. This work demonstrates the utility of interspecific hybrids for ageing research and provides a foundation for future work on the genetic architecture of interspecific lifespan variation.

Two New Pristionchus Species (Nematoda: Diplogastridae) from Taiwan are Part of a Species-cluster Representing the Closest Known Relatives of the Model Organism P. pacificus

Matthias Herrmann, Christian Weiler, Christian Rödelsperger, Natsumi Kanzaki, and Ralf J. Sommer (2016) Pristionchus pacificus is an important model organism in evolutionary biology and has been used to establish integrative studies that link developmental biology with ecology and population genetics. This species is part of the P. pacificus species-complex of the genus, many members of which occur in East Asia. While P. pacificus is hermaphroditic, the two most closely related known species are P. exspectatus and P. arcanus from Japan, both of which are gonochorists. P. exspectatus is so far the closest known relative of P. pacificus and thus, considered to represent the sister species. Here, we describe two new species of Pristionchus, P. taiwanensis and P. occultus from Taiwan using morphology, morphometrics, mating experiments and genome- wide sequence analysis. Both species are gonochorists and they are morphologically indistinguishable from P. exspectatus, P. arcanus and P. pacificus. However, reproductive isolation, namely the inability to produce interfertile hybrids, separates all species pairs in the species-complex. Phylogeny inferred from more than 700,000 genome-wide variable sites that were genotyped in all species suggest that P. taiwanensis and P. occultus are the sister species of P. arcanus and P. exspectatus, respectively. P. taiwanensis and P. occultus together with P. exspectatus and P. arcanus form a species-cluster with P. pacificus. The identification of these two novel gonochoristic species is invaluable for studies of population genetics, speciation, and macroevolution in the genus. We discuss the biogeography of Pristionchus in East Asia and the origin of hermaphroditism in the P. pacificus species-complex.

Life History Responses and Gene Expression Profiles of the Nematode Pristionchus pacificus Cultured on Cryptococcus Yeasts

Nematodes, the earth's most abundant metazoa are found in all ecosystems. In order to survive in diverse environments, they have evolved distinct feeding strategies and they can use different food sources. While some nematodes are specialists, including parasites of plants and animals, others such as Pristionchus pacificus are omnivorous feeders, which can live on a diet of bacteria, protozoans, fungi or yeast. In the wild, P. pacificus is often found in a necromenic association with beetles and is known to be able to feed on a variety of microbes as well as on nematode prey. However, in laboratory studies Escherichia coli OP50 has been used as standard food source, similar to investigations in Caenorhabditis elegans and it is unclear to what extent this biases the obtained results and how relevant findings are in real nature. To gain first insight into the variation in traits induced by a non-bacterial food source, we study Pristionchus-fungi interactions under laboratory conditions. After screening different yeast strains, we were able to maintain P. pacificus for at least 50-60 generations on Cryptococcus albidus and Cryptococcus curvatus. We describe life history traits of P. pacificus on both yeast strains, including developmental timing, survival and brood size. Despite a slight developmental delay and problems to digest yeast cells, which are both reflected at a transcriptomic level, all analyses support the potential of Cryptococcus strains as food source for P. pacificus. In summary, our work establishes two Cryptococcus strains as alternative food source for P. pacificus and shows change in various developmental, physiological and morphological traits, including the transcriptomic profiles.

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.

Specialization of a polyphenism switch gene following serial duplications in Pristionchus nematodes

Polyphenism is an extreme manifestation of developmental plasticity, requiring distinct developmental programs and the addition of a switch mechanism. Because the genetic basis of polyphenism switches has only begun to be understood, how their mechanisms arise is unclear. In the nematode Pristionchus pacificus, which has a mouthpart polyphenism specialized for alternative diets, a gene (eud-1) executing the polyphenism switch was recently identified as the product of lineage-specific duplications. Here, we infer the role of gene duplications in producing a switch gene. Using reverse genetics and population genetic analyses, we examine evidence for competing scenarios of degeneration and complementation, neutral evolution, and functional specialization. Of the daughter genes, eud-1 alone has assumed switch-like regulation of the mouth polyphenism. Measurements of life-history traits in single, double, and triple sulfatase mutants did not, given a benign environment, identify alternative or complementary roles for eud-1 paralogs. Although possible roles are still unknown, selection analyses of the sister species and 104 natural isolates of P. pacificus detected purifying selection on the genes, suggesting their functionality by their fixation and evolutionary maintenance. Our approach shows the tractability of reverse genetics in a nontraditional model system to study evolution by gene duplication.

A locus in Pristionchus pacificus that is responsible for the ability to give rise to fertile offspring at higher temperatures

Temperature is a stress factor that varies temporally and spatially, and can affect the fitness of cold-blooded organisms, leading to a loss of reproductive output; however, little is understood about the genetics behind the long-term response of organisms to temperature. Here, we approach this problem in the model nematode Pristionchus pacificus by utilising a large collection of natural isolates with diverse phenotypes. From this collection we identify two strains, one from California that can give rise to fertile offspring up to 28°C and one from Japan that is fertile up to 30°C. We show that the optimum temperature and the upper temperature limit for fertility is shifted higher in the Japanese strain suggesting that there is a mechanism that controls the temperature response of fertility across a range of temperatures. By crossing the two strains, and using genetic mapping, we identify a region on chromosome V that is responsible for maintaining fertility at higher temperatures. Thus, we conclude that fitness of P. pacificus at high temperature is under genetic control, suggesting that it could be subject to natural selection.

Draft Genome of the Scarab Beetle Oryctes borbonicus on La Réunion Island

Beetles represent the largest insect order and they display extreme morphological, ecological and behavioral diversity, which makes them ideal models for evolutionary studies. Here, we present the draft genome of the scarab beetle Oryctes borbonicus, which has a more basal phylogenetic position than the two previously sequenced pest species Tribolium castaneum and Dendroctonus ponderosae providing the potential for sequence polarization. Oryctes borbonicus is endemic to La Réunion, an island located in the Indian Ocean, and is the host of the nematode Pristionchus pacificus, a well-established model organism for integrative evolutionary biology. At 518 Mb, the O. borbonicus genome is substantially larger and encodes more genes than T. castaneum and D. ponderosae. We found that only 25% of the predicted genes of O. borbonicus are conserved as single copy genes across the nine investigated insect genomes, suggesting substantial gene turnover within insects. Even within beetles, up to 21% of genes are restricted to only one species, whereas most other genes have undergone lineage-specific duplications and losses. We illustrate lineage-specific duplications using detailed phylogenetic analysis of two gene families. This study serves as a reference point for insect/coleopteran genomics, although its original motivation was to find evidence for potential horizontal gene transfer (HGT) between O. borbonicus and P. pacificus. The latter was previously shown to be the recipient of multiple horizontally transferred genes including some genes from insect donors. However, our study failed to provide any clear evidence for additional HGTs between the two species.

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.

Functional Conservation and Divergence of daf-22 Paralogs in Pristionchus pacificus Dauer Development

Small-molecule signaling in nematode dauer formation has emerged as a major model to study chemical communication in development and evolution. Developmental arrest as nonfeeding and stress-resistant dauer larvae represents the major survival and dispersal strategy. Detailed studies in Caenorhabditis elegans and Pristionchus pacificus revealed that small-molecule communication changes rapidly in evolution resulting in extreme structural diversity of small-molecule compounds. In C. elegans, a blend of ascarosides constitutes the dauer pheromone, whereas the P. pacificus dauer pheromone includes additional paratosides and integrates building blocks from diverse primary metabolic pathways. Despite this complexity of small-molecule structures and functions, little is known about the biosynthesis of small molecules in nematodes outside C. elegans Here, we show that the genes encoding enzymes of the peroxisomal β-oxidation pathway involved in small-molecule biosynthesis evolve rapidly, including gene duplications and domain switching. The thiolase daf-22, the most downstream factor in C. elegans peroxisomal β-oxidation, has duplicated in P. pacificus, resulting in Ppa-daf-22.1, which still contains the sterol-carrier-protein (SCP) domain that was lost in C. elegans daf-22, and Ppa-daf-22.2. Using the CRISPR/Cas9 system, we induced mutations in both P. pacificus daf-22 genes and identified an unexpected complexity of functional conservation and divergence. Under well-fed conditions, ascaroside biosynthesis proceeds exclusively via Ppa-daf-22.1 In contrast, starvation conditions induce Ppa-daf-22.2 activity, resulting in the production of a specific subset of ascarosides. Gene expression studies indicate a reciprocal up-regulation of both Ppa-daf-22 genes, which is, however, independent of starvation. Thus, our study reveals an unexpected functional complexity of dauer development and evolution.