Divergent gene expression in the conserved dauer stage of the nematodes Pristionchus pacificus and Caenorhabditis elegans

Mass Spectrometry-Driven Discovery of Neuropeptides Mediating Nictation Behavior of Nematodes

Neuropeptides regulate animal physiology and behavior, making them widely studied targets of functional genetics research. While the field often relies on differential -omics approaches to build hypotheses, no such method exists for neuropeptidomics. It would nonetheless be valuable for studying behaviors suspected to be regulated by neuropeptides, especially when little information is otherwise available. This includes nictation, a phoretic strategy of Caenorhabditis elegans dauers that parallels host-finding strategies of infective juveniles of many pathogenic nematodes. We here developed a targeted peptidomics method for the model organism C. elegans and show that 161 quantified neuropeptides are more abundant in its dauer stage compared with L3 juveniles. Many of these have orthologs in the commercially relevant pathogenic nematode Steinernema carpocapsae, in whose infective juveniles, we identified 126 neuropeptides in total. Through further behavioral genetics experiments, we identify flp-7 and flp-11 as novel regulators of nictation. Our work advances knowledge on the genetics of nictation behavior and adds comparative neuropeptidomics as a tool to functional genetics workflows.

Alternative somatic and germline gene-regulatory strategies during starvation-induced developmental arrest

Nutrient availability governs growth and quiescence, and many animals arrest development when starved. Using C. elegans L1 arrest as a model, we show that gene expression changes deep into starvation. Surprisingly, relative expression of germline-enriched genes increases for days. We conditionally degrade the large subunit of RNA polymerase II using the auxin-inducible degron system and analyze absolute expression levels. We find that somatic transcription is required for survival, but the germline maintains transcriptional quiescence. Thousands of genes are continuously transcribed in the soma, though their absolute abundance declines, such that relative expression of germline transcripts increases given extreme transcript stability. Aberrantly activating transcription in starved germ cells compromises reproduction, demonstrating important physiological function of transcriptional quiescence. This work reveals alternative somatic and germline gene-regulatory strategies during starvation, with the soma maintaining a robust transcriptional response to support survival and the germline maintaining transcriptional quiescence to support future reproductive success.

Webster et al. show that the transcriptional response to starvation is mounted early in larval somatic cells supporting survival but that it wanes over time. In contrast, they show that the germline remains transcriptionally quiescent deep into starvation, supporting reproductive potential, while maintaining its transcriptome via transcript stability.

Kinetic characterization and thermostability of C. elegans cytoplasmic and mitochondrial malate dehydrogenases

Malate dehydrogenase (MDH) catalyzes the conversion of NAD⁺ and malate to NADH and oxaloacetate in the citric acid cycle. Eukaryotes have one MDH isozyme that is imported into the mitochondria and one in the cytoplasm. We overexpressed and purified Caenorhabditis elegans cytoplasmic MDH-1 and mitochondrial MDH-2 in E. coli. Our goal was to compare the kinetic and structural properties of these enzymes because C. elegans can survive adverse environmental conditions, such as lack of food and elevated temperatures. In steady-state enzyme kinetics assays, we measured K_M values for oxaloacetate of 54 and 52 μM and K_M values for NADH of 61 and 107 μM for MDH-1 and MDH-2, respectively. We partially purified endogenous MDH-1 and MDH-2 from a mixed population of worms and separated them using anion exchange chromatography. Both endogenous enzymes had a K_M for oxaloacetate similar to that of the corresponding recombinant enzyme. Recombinant MDH-1 and MDH-2 had maximum activity at 40 °C and 35 °C, respectively. In a thermotolerance assay, MDH-1 was much more thermostable than MDH-2. Protein homology modeling predicted that MDH-1 had more intersubunit salt-bridges than mammalian MDH1 enzymes, and these ionic interactions may contribute to its thermostability. In contrast, the MDH-2 homology model predicted fewer intersubunit ionic interactions compared to mammalian MDH2 enzymes. These results suggest that the increased stability of MDH-1 may facilitate its ability to remain active in adverse environmental conditions. In contrast, MDH-2 may use other strategies, such as protein binding partners, to function under similar conditions.

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

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

Rhabditophanes diutinus a parthenogenetic clade IV nematode with dauer larvae

Comparative studies using non-parasitic model species such as Caenorhabditis elegans, have been very helpful in investigating the basic biology and evolution of parasitic nematodes. However, as phylogenetic distance increases, these comparisons become more difficult, particularly when outside of the nematode clade to which C. elegans belongs (V). One of the reasons C. elegans has nevertheless been used for these comparisons, is that closely related well characterized free-living species that can serve as models for parasites of interest are frequently not available. The Clade IV parasitic nematodes Strongyloides are of great research interest due to their life cycle and other unique biological features, as well as their medical and veterinary importance. Rhabditophanes, a closely related free-living genus, forms part of the Strongyloidoidea nematode superfamily. Rhabditophanes diutinus (= R. sp. KR3021) was included in the recent comparative genomic analysis of the Strongyloididae, providing some insight into the genomic nature of parasitism. However, very little is known about this species, limiting its usefulness as a research model. Here we provide a species description, name the species as R. diutinus and investigate its life cycle and subsequently gene expression in multiple life stages. We identified two previously unreported starvation induced life stages: dauer larvae and arrested J2 (J2A) larvae. The dauer larvae are morphologically similar to and are the same developmental stage as dauers in C. elegans and infective larvae in Strongyloides. As in C. elegans and Strongyloides, dauer formation is inhibited by treatment with dafachronic acid, indicating some genetic control mechanisms are conserved. Similarly, the expression patterns of putative dauer/infective larva control genes resemble each other, in particular between R. diutinus and Strongyloides spp. These findings illustrate and increase the usefulness of R. diutinus as a non-parasitic, easy to work with model species for the Strongyloididae for studying the evolution of parasitism as well as many aspects of the biology of Strongyloides spp, in particular the formation of infective larvae.

Parasitic worms are an issue of great medical, veterinary, agricultural and economic importance, yet little is known about how worms become parasites. Comparative studies with non-parasitic model species like C. elegans have been useful, however, this usefulness decreases as the evolutionary distance between the species increases. One way to combat this is by having more well-studied closely related species to parasites of interest. To address this, we provide information about Rhabditophanes diutinus, a free-living nematode that is part of the same superfamily as the medically and veterinary important Strongyloides parasites. We provide analysis on its life cycle, in particular on two starvation induced life stages, along with gene expression data. Overall, this important information illustrates and improves the use of R. diutinus, as a non-parasitic model species for studying parasite evolution and basic biology within Strongyloides.

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.

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.

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.

Key genes involved in desiccation tolerance and dormancy across life forms

Desiccation tolerance (DT, the ability of certain organisms to survive severe dehydration) was a key trait in the evolution of life in terrestrial environments. Likely, the development of desiccation-tolerant life forms was accompanied by the acquisition of dormancy or a dormancy-like stage as a second powerful adaptation to cope with variations in the terrestrial environment. These naturally stress tolerant life forms may be a good source of genetic information to generate stress tolerant crops to face a future with predicted higher occurrence of drought. By mining for key genes and mechanisms related to DT and dormancy conserved across different species and life forms, unique candidate key genes may be identified. Here we identify several of these putative key genes, shared among multiple organisms, encoding for proteins involved in protection, growth and energy metabolism. Mutating a selection of these genes in the model plant Arabidopsis thaliana resulted in clear DT-, dormancy- and other seed-associated phenotypes, showing the efficiency and power of our approach and paves the way for the development of drought-stress tolerant crops. Our analysis supports a co-evolution of DT and dormancy by shared mechanisms that favour survival and adaptation to ever-changing environments with strong seasonal fluctuations.

Analysis of C. elegans muscle transcriptome using trans-splicing-based RNA tagging (SRT)

Current approaches to profiling tissue-specific gene expression in C. elegans require delicate manipulation and are difficult under certain conditions, e.g. from dauer or aging worms. We have developed an easy and robust method for tissue-specific RNA-seq by taking advantage of the endogenous trans-splicing process. In this method, transgenic worms are generated in which a spliced leader (SL) RNA gene is fused with a sequence tag and driven by a tissue-specific promoter. Only in the tissue of interest, the tagged SL RNA gene is transcribed and then trans-spliced onto mRNAs. The tag allows enrichment and sequencing of mRNAs from that tissue only. As a proof of principle, we profiled the muscle transcriptome, which showed high coverage and efficient enrichment of muscle specific genes, with low background noise. To demonstrate the robustness of our method, we profiled muscle gene expression in dauer larvae and aging worms, revealing gene expression changes consistent with the physiology of these stages. The resulting muscle transcriptome also revealed 461 novel RNA transcripts, likely muscle-expressed long non-coding RNAs. In summary, the splicing-based RNA tagging (SRT) method provides a convenient and robust tool to profile trans-spliced genes and identify novel transcripts in a tissue-specific manner, with a low false positive rate.

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

Background

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

Results

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

Conclusions

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

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

A transcriptomic insight into the infective juvenile stage of the insect parasitic nematode, Heterorhabditis indica

Background

Nematodes are the most numerous animals in the soil. Insect parasitic nematodes of the genus Heterorhabditis are capable of selectively seeking, infecting and killing their insect-hosts in the soil. The infective juvenile (IJ) stage of the Heterorhabditis nematodes is analogous to Caenorhabditis elegans dauer juvenile stage, which remains in ‘arrested development’ till it finds and infects a new insect-host in the soil. H. indica is the most prevalent species of Heterorhabditis in India. To understand the genes and molecular processes that govern the biology of the IJ stage, and to create a resource to facilitate functional genomics and genetic exploration, we sequenced the transcriptome of H. indica IJs.

Results

The de-novo sequence assembly using Velvet-Oases pipeline resulted in 13,593 unique transcripts at N50 of 1,371 bp, of which 53 % were annotated by blastx. H. indica transcripts showed higher orthology with parasitic nematodes as compared to free living nematodes. In-silico expression analysis showed 30 % of transcripts expressing with ≥100 FPKM value. All the four canonical dauer formation pathways like cGMP-PKG, insulin, dafachronic acid and TGF-β were active in the IJ stage. Several other signaling pathways were highly represented in the transcriptome. Twenty-four orthologs of C. elegans RNAi pathway effector genes were discovered in H. indica, including nrde-3 that is reported for the first time in any of the parasitic nematodes. An ortholog of C. elegans tol-1 was also identified. Further, 272 kinases belonging to 137 groups, and several previously unidentified members of important gene classes were identified.

Conclusions

We generated high-quality transcriptome sequence data from H. indica IJs for the first time. The transcripts showed high similarity with the parasitic nematodes, M. hapla, and A. suum as opposed to C. elegans, a species to which H. indica is more closely related. The high representation of transcripts from several signaling pathways in the IJs indicates that despite being a developmentally arrested stage; IJs are a hotbed of signaling and are actively interacting with their environment.

Electronic supplementary material

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

Comparative genomics of Steinernema reveals deeply conserved gene regulatory networks

BACKGROUND

Parasitism is a major ecological niche for a variety of nematodes. Multiple nematode lineages have specialized as pathogens, including deadly parasites of insects that are used in biological control. We have sequenced and analyzed the draft genomes and transcriptomes of the entomopathogenic nematode Steinernema carpocapsae and four congeners (S. scapterisci, S. monticolum, S. feltiae, and S. glaseri).

RESULTS

We used these genomes to establish phylogenetic relationships, explore gene conservation across species, and identify genes uniquely expanded in insect parasites. Protein domain analysis in Steinernema revealed a striking expansion of numerous putative parasitism genes, including certain protease and protease inhibitor families, as well as fatty acid- and retinol-binding proteins. Stage-specific gene expression of some of these expanded families further supports the notion that they are involved in insect parasitism by Steinernema. We show that sets of novel conserved non-coding regulatory motifs are associated with orthologous genes in Steinernema and Caenorhabditis.

CONCLUSIONS

We have identified a set of expanded gene families that are likely to be involved in parasitism. We have also identified a set of non-coding motifs associated with groups of orthologous genes in Steinernema and Caenorhabditis involved in neurogenesis and embryonic development that are likely part of conserved protein-DNA relationships shared between these two genera.

Ancient gene duplications have shaped developmental stage-specific expression in Pristionchus pacificus

Background

The development of multicellular organisms is accompanied by gene expression changes in differentiating cells. Profiling stage-specific expression during development may reveal important insights into gene sets that contributed to the morphological diversity across the animal kingdom.

Results

We sequenced RNA-seq libraries throughout a developmental timecourse of the nematode Pristionchus pacificus. The transcriptomes reflect early larval stages, adult worms including late larvae, and growth-arrested dauer larvae and allowed the identification of developmentally regulated gene clusters. Our data reveals similar trends as previous transcriptome profiling of dauer worms and represents the first expression data for early larvae in P. pacificus. Gene expression clusters characterizing early larval stages show most significant enrichments of chaperones, while collagens are most significantly enriched in transcriptomes of late larvae and adult worms. By combining expression data with phylogenetic analysis, we found that developmentally regulated genes are found in paralogous clusters that have arisen through lineage-specific duplications after the split from the Caenorhabditis elegans branch.

Conclusions

We propose that gene duplications of developmentally regulated genes represent a plausible evolutionary mechanism to increase the dosage of stage-specific expression. Consequently, this may contribute to the substantial divergence in expression profiles that has been observed across larger evolutionary time scales.

Electronic supplementary material

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

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

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

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

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

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

Many nematodes form dauer larvae when exposed to unfavorable conditions, representing an example of phenotypic plasticity and a major survival and dispersal strategy. In Caenorhabditis elegans, the regulation of dauer induction is a model for pheromone, insulin, and steroid-hormone signaling. Recent studies in Pristionchus pacificus revealed substantial natural variation in various aspects of dauer development, i.e. pheromone production and sensing and dauer longevity and fitness. One intriguing example is a strain from Ohio, having extremely long-lived dauers associated with very high fitness and often forming the most dauers in response to other strains' pheromones, including the reference strain from California. While such examples have been suggested to represent intraspecific competition among strains, the molecular mechanisms underlying these dauer-associated patterns are currently unknown. We generated recombinant-inbred-lines between the Californian and Ohioan strains and used quantitative-trait-loci analysis to investigate the molecular mechanism determining natural variation in dauer development. Surprisingly, we discovered that the orphan gene dauerless controls dauer formation by copy number variation. The Ohioan strain has one dauerless copy causing high dauer formation, whereas the Californian strain has two copies, resulting in strongly reduced dauer formation. Transgenic animals expressing multiple copies do not form dauers. dauerless is exclusively expressed in CAN neurons, and both CAN ablation and dauerless mutations increase dauer formation. Strikingly, dauerless underwent several duplications and acts in parallel or downstream of steroid-hormone signaling but upstream of the nuclear-hormone-receptor daf-12. We identified the novel or fast-evolving gene dauerless as inhibitor of dauer development. Our findings reveal the importance of gene duplications and copy number variations for orphan gene function and suggest daf-12 as major target for dauer regulation. We discuss the consequences of the novel vs. fast-evolving nature of orphans for the evolution of developmental networks and their role in natural variation and intraspecific competition.

Pan-phylum Comparison of Nematode Metabolic Potential

Nematodes are among the most important causative pathogens of neglected tropical diseases. The increased availability of genomic and transcriptomic data for many understudied nematode species provides a great opportunity to investigate different aspects of their biology. Increasingly, metabolic potential of pathogens is recognized as a critical determinant governing their development, growth and pathogenicity. Comparing metabolic potential among species with distinct trophic ecologies can provide insights on overall biology or molecular adaptations. Furthermore, ascertaining gene expression at pathway level can help in understanding metabolic dynamics over development. Comparison of biochemical pathways (or subpathways, i.e. pathway modules) among related species can also retrospectively indicate potential mistakes in gene-calling and functional annotation. We show with numerous illustrative case studies that comparisons at the level of pathway modules have the potential to uncover biological insights while remaining computationally tractable. Here, we reconstruct and compare metabolic modules found in the deduced proteomes of 13 nematodes and 10 non-nematode species (including hosts of the parasitic nematode species). We observed that the metabolic potential is, in general, concomitant with phylogenetic and/or ecological similarity. Varied metabolic strategies are required among the nematodes, with only 8 out of 51 pathway modules being completely conserved. Enzyme comparison based on topology of metabolic modules uncovered diversification between parasite and host that can potentially guide therapeutic intervention. Gene expression data from 4 nematode species were used to study metabolic dynamics over their life cycles. We report unexpected differential metabolism between immature and mature microfilariae of the human filarial parasite Brugia malayi. A set of genes potentially important for parasitism is also reported, based on an analysis of gene expression in C. elegans and the human hookworm Necator americanus. We illustrate how analyzing and comparing metabolism at the level of pathway modules can improve existing knowledge of nematode metabolic potential and can provide parasitism related insights. Our reconstruction and comparison of nematode metabolic pathways at a pan-phylum and inter-phylum level enabled determination of phylogenetic restrictions and differential expression of pathways. A visualization of our results is available at http://nematode.net and the program for identification of module completeness (modDFS) is freely available at SourceForge. The methods reported will help biologists to predict biochemical potential of any organism with available deduced proteome, to direct experiments and test hypotheses.

We reconstructed metabolic pathways of 23 organisms including 13 nematode species, using their complete deduced protein coding sequences and compared them to 10 non-nematodes. We observed that metabolic potential availability is concomitant with phylogenetic and/or ecological similarity, with the exceptions providing interesting case studies. We also studied changes in metabolic profiles under different developmental stages of 4 nematode species using stage-specific transcriptomic data. A comparison of the variation patterns in these profiles led to recognition of modules that share metabolic profiles at various life-cycle stages or during development. The undertaken analysis improved genome annotation and the obtained results provided insight into parasitism, resulting in identification of taxonomically-restricted pathways and enzymes that may provide new mechanisms for control of nematode infections.

Genome-wide analysis of trans-splicing in the nematode Pristionchus pacificus unravels conserved gene functions for germline and dauer development in divergent operons

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader-specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

A functional genomic screen for evolutionarily conserved genes required for lifespan and immunity in germline-deficient C. elegans

The reproductive system regulates lifespan in insects, nematodes and vertebrates. In Caenorhabditis elegans removal of germline increases lifespan by 60% which is dependent upon insulin signaling, nuclear hormone signaling, autophagy and fat metabolism and their microRNA-regulators. Germline-deficient C. elegans are also more resistant to various bacterial pathogens but the underlying molecular mechanisms are largely unknown. Firstly, we demonstrate that previously identified genes that regulate the extended lifespan of germline-deficient C. elegans (daf-2, daf-16, daf-12, tcer-1, mir-7.1 and nhr-80) are also essential for resistance to the pathogenic bacterium Xenorhabdus nematophila. We then use a novel unbiased approach combining laser cell ablation, whole genome microarrays, RNAi screening and exposure to X. nematophila to generate a comprehensive genome-wide catalog of genes potentially required for increased lifespan and innate immunity in germline-deficient C. elegans. We find 3,440 genes to be upregulated in C. elegans germline-deficient animals in a gonad dependent manner, which are significantly enriched for genes involved in insulin signaling, fatty acid desaturation, translation elongation and proteasome complex function. Using RNAi against a subset of 150 candidate genes selected from the microarray results, we show that the upregulated genes such as transcription factor DAF-16/FOXO, the PTEN homolog lipid phosphatase DAF-18 and several components of the proteasome complex (rpn-6.1, rpn-7, rpn-9, rpn-10, rpt-6, pbs-3 and pbs-6) are essential for both lifespan and immunity of germline deficient animals. We also identify a novel role for genes including par-5 and T12G3.6 in both lifespan-extension and increased survival on X. nematophila. From an evolutionary perspective, most of the genes differentially expressed in germline deficient C. elegans also show a conserved expression pattern in germline deficient Pristionchus pacificus, a nematode species that diverged from C. elegans 250-400 MYA.

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

The hermaphroditic nematode Pristionchus pacificus is an established model system for comparative studies with Caenorhabditis elegans in developmental biology, ecology, and population genetics. In this study, we present whole-genome sequencing data of 104 P. pacificus strains and the draft assembly of the obligate outcrossing sister species P. exspectatus. We characterize genetic diversity within P. pacificus and investigate the population genetic processes shaping this diversity. P. pacificus is 10 times more diverse than C. elegans and exhibits substantial population structure that allows us to probe its evolution on multiple timescales. Consistent with reduced effective recombination in this self-fertilizing species, we find haplotype blocks that span several megabases. Using the P. exspectatus genome as an outgroup, we polarized variation in P. pacificus and found a site frequency spectrum (SFS) that decays more rapidly than expected in neutral models. The SFS at putatively neutral sites is U shaped, which is a characteristic feature of pervasive linked selection. Based on the additional findings (i) that the majority of nonsynonymous variation is eliminated over timescales on the order of the separation between clades, (ii) that diversity is reduced in gene-rich regions, and (iii) that highly differentiated clades show very similar patterns of diversity, we conclude that purifying selection on many mutations with weak effects is a major force shaping genetic diversity in P. pacificus.

New role for DCR-1/dicer in Caenorhabditis elegans innate immunity against the highly virulent bacterium Bacillus thuringiensis DB27

Bacillus thuringiensis produces toxins that target invertebrates, including Caenorhabditis elegans. Virulence of Bacillus strains is often highly specific, such that B. thuringiensis strain DB27 is highly pathogenic to C. elegans but shows no virulence for another model nematode, Pristionchus pacificus. To uncover the underlying mechanisms of the differential responses of the two nematodes to B. thuringiensis DB27 and to reveal the C. elegans defense mechanisms against this pathogen, we conducted a genetic screen for C. elegans mutants resistant to B. thuringiensis DB27. Here, we describe a B. thuringiensis DB27-resistant C. elegans mutant that is identical to nasp-1, which encodes the C. elegans homolog of the nuclear-autoantigenic-sperm protein. Gene expression analysis indicated a substantial overlap between the genes downregulated in the nasp-1 mutant and targets of C. elegans dcr-1/Dicer, suggesting that dcr-1 is repressed in nasp-1 mutants, which was confirmed by quantitative PCR. Consistent with this, the nasp-1 mutant exhibits RNA interference (RNAi) deficiency and reduced longevity similar to those of a dcr-1 mutant. Building on these surprising findings, we further explored a potential role for dcr-1 in C. elegans innate immunity. We show that dcr-1 mutant alleles deficient in microRNA (miRNA) processing, but not those deficient only in RNAi, are resistant to B. thuringiensis DB27. Furthermore, dcr-1 overexpression rescues the nasp-1 mutant's resistance, suggesting that repression of dcr-1 determines the nasp-1 mutant's resistance. Additionally, we identified the collagen-encoding gene col-92 as one of the downstream effectors of nasp-1 that play an important role in resistance to DB27. Taken together, these results uncover a previously unknown role for DCR-1/Dicer in C. elegans antibacterial immunity that is largely associated with miRNA processing.

The nematode Pristionchus pacificus as a model system for integrative studies in evolutionary biology

Comprehensive studies of evolution have historically been hampered by the division among disciplines. Now, as biology moves towards an '-omics' era, it is more important than ever to tackle the evolution of function and form by considering all those research areas involved in the regulation of phenotypes. Here, we review recent attempts to establish the nematode Pristionchus pacificus as a model organism that allows integrative studies of development and evo-devo, with ecology and population genetics. Originally developed for comparative study with the nematode Caenorhabditis elegans, P. pacificus provided insight into developmental pathways including dauer formation, vulva and gonad development, chemosensation, innate immunity and neurobiology. Its subsequent discovery across a wide geographic distribution in association with scarab beetles enabled its evaluation in a biogeographic context. Development of an evolutionary field station on La Réunion Island, where P. pacificus is present in high abundance across a number of widespread habitat types, allows examination of the microfacets of evolution - processes of natural selection, adaptation and drift among populations can now be examined in this island setting. The combination of laboratory-based functional studies with fieldwork in P. pacificus has the long-term prospective to provide both proximate (mechanistic) and ultimate (evolutionary and ecological) causation and might therefore help to overcome the long-term divide between major areas in biology.

RNAseq analysis of the parasitic nematode Strongyloides stercoralis reveals divergent regulation of canonical dauer pathways

The infectious form of many parasitic nematodes, which afflict over one billion people globally, is a developmentally arrested third-stage larva (L3i). The parasitic nematode Strongyloides stercoralis differs from other nematode species that infect humans, in that its life cycle includes both parasitic and free-living forms, which can be leveraged to investigate the mechanisms of L3i arrest and activation. The free-living nematode Caenorhabditis elegans has a similar developmentally arrested larval form, the dauer, whose formation is controlled by four pathways: cyclic GMP (cGMP) signaling, insulin/IGF-1-like signaling (IIS), transforming growth factor β (TGFβ) signaling, and biosynthesis of dafachronic acid (DA) ligands that regulate a nuclear hormone receptor. We hypothesized that homologous pathways are present in S. stercoralis, have similar developmental regulation, and are involved in L3i arrest and activation. To test this, we undertook a deep-sequencing study of the polyadenylated transcriptome, generating over 2.3 billion paired-end reads from seven developmental stages. We constructed developmental expression profiles for S. stercoralis homologs of C. elegans dauer genes identified by BLAST searches of the S. stercoralis genome as well as de novo assembled transcripts. Intriguingly, genes encoding cGMP pathway components were coordinately up-regulated in L3i. In comparison to C. elegans, S. stercoralis has a paucity of genes encoding IIS ligands, several of which have abundance profiles suggesting involvement in L3i development. We also identified seven S. stercoralis genes encoding homologs of the single C. elegans dauer regulatory TGFβ ligand, three of which are only expressed in L3i. Putative DA biosynthetic genes did not appear to be coordinately regulated in L3i development. Our data suggest that while dauer pathway genes are present in S. stercoralis and may play a role in L3i development, there are significant differences between the two species. Understanding the mechanisms governing L3i development may lead to novel treatment and control strategies.

Parasitic nematodes infect over one billion people worldwide and cause many diseases, including strongyloidiasis, filariasis, and hookworm disease. For many of these parasites, including Strongyloides stercoralis, the infectious form is a developmentally arrested and long-lived thirdstage larva (L3i). Upon encountering a host, L3i quickly resume development and mature into parasitic adults. In the free-living nematode Caenorhabditis elegans, a similar developmentally arrested third-stage larva, known as the dauer, is regulated by four key cellular mechanisms. We hypothesized that similar cellular mechanisms control L3i arrest and activation. Therefore, we used deep-sequencing technology to characterize the S. stercoralis transcriptome (RNAseq), which allowed us to identify S. stercoralis homologs of components of these four mechanisms and examine their temporal regulation. We found similar temporal regulation between S. stercoralis and C. elegans for components of two mechanisms, but dissimilar temporal regulation for two others, suggesting conserved as well as novel modes of developmental regulation for L3i. Understanding L3i development may lead to novel control strategies as well as new treatments for strongyloidiasis and other diseases caused by parasitic nematodes.

System wide analysis of the evolution of innate immunity in the nematode model species Caenorhabditis elegans and Pristionchus pacificus

The evolution of genetic mechanisms used to combat bacterial infections is critical for the survival of animals and plants, yet how these genes evolved to produce a robust defense system is poorly understood. Studies of the nematode Caenorhabditis elegans have uncovered a plethora of genetic regulators and effectors responsible for surviving pathogens. However, comparative studies utilizing other free-living nematodes and therefore providing an insight into the evolution of innate immunity have been lacking. Here, we take a systems biology approach and use whole genome microarrays to profile the transcriptional response of C. elegans and the necromenic nematode Pristionchus pacificus after exposure to the four different pathogens Serratia marcescens, Xenorhabdus nematophila, Staphylococcus aureus and Bacillus thuringiensis DB27. C. elegans is susceptible to all four pathogens whilst P. pacificus is only susceptible to S. marcescens and X. nematophila. We show an unexpected level of specificity in host responses to distinct pathogens within and across species, revealing an enormous complexity of effectors of innate immunity. Functional domains enriched in the transcriptomes on different pathogens are similar within a nematode species but different across them, suggesting differences in pathogen sensing and response networks. We find translation inhibition to be a potentially conserved response to gram-negative pathogens in both the nematodes. Further computational analysis indicates that both nematodes when fed on pathogens up-regulate genes known to be involved in other stress responses like heat shock, oxidative and osmotic stress, and genes regulated by DAF-16/FOXO and TGF-beta pathways. This study presents a platform for comparative systems analysis of two nematode model species, and a catalog of genes involved in the evolution of nematode immunity and identifies both pathogen specific and pan-pathogen responses. We discuss the potential effects of ecology on evolution of downstream effectors and upstream regulators on evolution of nematode innate immunity.

Phosphoproteome of Pristionchus pacificus provides insights into architecture of signaling networks in nematode models

Pristionchus pacificus is a nematode that is increasingly used as a model organism in evolutionary biology. The genome of P. pacificus differs markedly from that of C. elegans, with a high number of orphan genes that are restricted to P. pacificus and have no homologs in other species. To gain insight into the architecture of signal transduction networks in model nematodes, we performed a large-scale qualitative phosphoproteome analysis of P. pacificus. Using two-stage enrichment of phosphopeptides from a digest of P. pacificus proteins and their subsequent analysis via high accuracy MS, we detected and localized 6,809 phosphorylation events on 2,508 proteins. We compared the detected P. pacificus phosphoproteome to the recently published phosphoproteome of C. elegans. The overall numbers and functional classes of phosphoproteins were similar between the two organisms. Interestingly, the products of orphan genes were significantly underrepresented among the detected P. pacificus phosphoproteins. We defined the theoretical kinome of P. pacificus and compared it to that of C. elegans. While tyrosine kinases were slightly underrepresented in the kinome of P. pacificus, all major classes of kinases were present in both organisms. Application of our kinome annotation to a recent transcriptomic study of dauer and mixed stage populations showed that Ser/Thr and Tyr kinases show similar expression levels in P. pacificus but not in C. elegans. This study presents the first systematic comparison of phosphoproteomes and kinomes of two model nematodes and, as such, will be a useful resource for comparative studies of their signal transduction networks.

Genome-wide analysis of germline signaling genes regulating longevity and innate immunity in the nematode Pristionchus pacificus

Removal of the reproductive system of many animals including fish, flies, nematodes, mice and humans can increase lifespan through mechanisms largely unknown. The abrogation of the germline in Caenorhabditis elegans increases longevity by 60% due to a signal emitted from the somatic gonad. Apart from increased longevity, germline-less C. elegans is also resistant to other environmental stressors such as feeding on bacterial pathogens. However, the evolutionary conservation of this pathogen resistance, its genetic basis and an understanding of genes involved in producing this extraordinary survival phenotype are currently unknown. To study these evolutionary aspects we used the necromenic nematode Pristionchus pacificus, which is a genetic model system used in comparison to C. elegans. By ablation of germline precursor cells and subsequent feeding on the pathogen Serratia marcescens we discovered that P. pacificus shows remarkable resistance to bacterial pathogens and that this response is evolutionarily conserved across the Genus Pristionchus. To gain a mechanistic understanding of the increased resistance to bacterial pathogens and longevity in germline-ablated P. pacificus we used whole genome microarrays to profile the transcriptional response comparing germline ablated versus un-ablated animals when fed S. marcescens. We show that lipid metabolism, maintenance of the proteasome, insulin signaling and nuclear pore complexes are essential for germline deficient phenotypes with more than 3,300 genes being differentially expressed. In contrast, gene expression of germline-less P. pacificus on E. coli (longevity) and S. marcescens (immunity) is very similar with only 244 genes differentially expressed indicating that longevity is due to abundant gene expression also involved in immunity. By testing existing mutants of Ppa-DAF-16/FOXO and the nuclear hormone receptor Ppa-DAF-12 we show a conserved function of both genes in resistance to bacterial pathogens and longevity. This is the first study to show that the influence of the reproductive system on extending lifespan and innate immunity is conserved in evolution.