The Same or Not the Same: Lineage-Specific Gene Expansions and Homology Relationships in Multigene Families in Nematodes

Sex-specificity of the C. elegans metabolome

Recent studies of animal metabolism have revealed large numbers of novel metabolites that are involved in all aspects of organismal biology, but it is unclear to what extent metabolomes differ between sexes. Here, using untargeted comparative metabolomics for the analysis of wildtype animals and sex determination mutants, we show that C. elegans hermaphrodites and males exhibit pervasive metabolomic differences. Several hundred small molecules are produced exclusively or in much larger amounts in one sex, including a host of previously unreported metabolites that incorporate building blocks from nucleoside, carbohydrate, lipid, and amino acid metabolism. A subset of male-enriched metabolites is specifically associated with the presence of a male germline, whereas enrichment of other compounds requires a male soma. Further, we show that one of the male germline-dependent metabolites, an unusual dipeptide incorporating N,N-dimethyltryptophan, increases food consumption, reduces lifespan, and accelerates the last stage of larval development in hermaphrodites. Our results serve as a foundation for mechanistic studies of how the genetic sex of soma and germline shape the C. elegans metabolome and provide a blueprint for the discovery of sex-dependent metabolites in other animals.

Comparative genomics and community curation further improve gene annotations in the nematode Pristionchus pacificus

Background

Nematode model organisms such as Caenorhabditis elegans and Pristionchus pacificus are powerful systems for studying the evolution of gene function at a mechanistic level. However, the identification of P. pacificus orthologs of candidate genes known from C. elegans is complicated by the discrepancy in the quality of gene annotations, a common problem in nematode and invertebrate genomics.

Results

Here, we combine comparative genomic screens for suspicious gene models with community-based curation to further improve the quality of gene annotations in P. pacificus. We extend previous curations of one-to-one orthologs to larger gene families and also orphan genes. Cross-species comparisons of protein lengths, screens for atypical domain combinations and species-specific orphan genes resulted in 4311 candidate genes that were subject to community-based curation. Corrections for 2946 gene models were implemented in a new version of the P. pacificus gene annotations. The new set of gene annotations contains 28,896 genes and has a single copy ortholog completeness level of 97.6%.

Conclusions

Our work demonstrates the effectiveness of comparative genomic screens to identify suspicious gene models and the scalability of community-based approaches to improve the quality of thousands of gene models. Similar community-based approaches can help to improve the quality of gene annotations in other invertebrate species, including parasitic nematodes.

Characterisation of a Teladorsagia circumcincta glutathione transferase

A 615 bp full length cDNA encoding a Teladorsagia circumcincta glutathione transferase (TcGST) was cloned, expressed in Escherichia coli and the recombinant protein purified and its kinetic properties determined. The predicted protein consisted of 205 amino acids and was present as a single band of about 24 kDa on SDS-PAGE. Multiple alignments of the protein sequence of TcGST with homologues from other helminths showed that the highest identity of 53-68% with haem-binding nematode proteins designated as members of the nu class of GSTs. Substrate binding sites and conserved regions were identified and were generally conserved. The predicted 3-dimensional structures of TcGST and HcGST revealed highly open binding cavities typical of this class of GST, considered to allow greater accessibility to diverse ligands compared with other classes of GST. At 25 °C, the optimum pH for TcGST activity was pH 7, the V_max was 1535 ± 33 nmoles.min^-1. mg^-1 protein and the apparent K_m for the substrate 1-chloro-2,4-dinitrobenzene (CDNB) was 0.22 ± 0.01 mM (mean ± SD, n = 2). Antibodies in both serum and saliva from field-immune, but not nematode-naïve, sheep, recognised recombinant TcGST in enzyme-linked immunosorbent assays. The recognition of the recombinant protein by antibodies generated by exposure of sheep to the native enzyme indicates similar antigenicity of the two proteins. These findings could aid in the design of novel drugs and vaccine antigens for economically important parasites of livestock.

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

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

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.

Presence or absence? Primary structure, regioselectivity and evolution of Δ12/ω3 fatty acid desaturases in nematodes

For vertebrates, the adequate supply of polyunsaturated fatty acids (PUFA) by the diet, in particular ω3 long-chain PUFA, is considered essential for neural development, growth and reproduction. In contrast to aquatic ecosystems, ω3 long-chain PUFA apparently are not widely available in the terrestrial food chain. Their de novo synthesis requires the presence of Δ12 and ω3 fatty acid desaturase enzymes, which are absent in vertebrates but present, for example, in the nematode Caenorhabditis elegans (FAT-2 and FAT-1). This raises the question if soil-dwelling nematodes offer substantial supply of these valuable nutritional compounds in terrestrial food webs. BLAST searches in available nematode genomes revealed the existence of fat-2 like genes in almost all clade III-V species, but failed to identify orthologs in clade I-II nematodes. An additional RT-PCR screen across soil-dwelling nematode species identified six novel fat-2 like genes. Hints for the genetic basis of a ω3 (fat-1) desaturase activity was found only in selected clade IV-V species, but not in clades I to III nematodes. Fatty acid pattern analyses following a PUFA-free cultivation and enzymatic characterization of six selected fat-2 or fat-1 like desaturases in yeast confirmed the findings from the genetic approaches. Thus, in similar soil habitats, taxa exist that can synthesize ω3 long-chain PUFA (as Panagrolaimus, Mesorhabditis and Caenorhabditis) whereas others are unable to do so (Acrobeloides, Cephalobus and Oscheius). While these nematodes do not differ in trophic position or major diet, distinction in reproduction mode may have led to the observed variations in desaturase genes.

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.

Rapid Gene Family Evolution of a Nematode Sperm Protein Despite Sequence Hyper-conservation

Reproductive proteins are often observed to be the most rapidly evolving elements within eukaryotic genomes. The major sperm protein (MSP) is unique to the phylum Nematoda and is required for proper sperm locomotion and fertilization. Here, we annotate the MSP gene family and analyze their molecular evolution in 10 representative species across Nematoda. We show that MSPs are hyper-conserved across the phylum, having maintained an amino acid sequence identity of 83.5–97.7% for over 500 million years. This extremely slow rate of evolution makes MSPs some of the most highly conserved genes yet identified. However, at the gene family level, we show hyper-variability in both gene copy number and genomic position within species, suggesting rapid, lineage-specific gene family evolution. Additionally, we find evidence that extensive gene conversion contributes to the maintenance of sequence identity within chromosome-level clusters of MSP genes. Thus, while not conforming to the standard expectation for the evolution of reproductive proteins, our analysis of the molecular evolution of the MSP gene family is nonetheless consistent with the widely repeatable observation that reproductive proteins evolve rapidly, in this case in terms of the genomic properties of gene structure, copy number, and genomic organization. This unusual evolutionary pattern is likely generated by strong pleiotropic constraints acting on these genes at the sequence level, balanced against expansion at the level of the whole gene family.

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.

Serotonin Drives Predatory Feeding Behavior via Synchronous Feeding Rhythms in the Nematode Pristionchus pacificus

Feeding behaviors in a wide range of animals are regulated by the neurotransmitter serotonin, although the exact neural circuits and associated mechanism are often unknown. The nematode Pristionchus pacificus can kill other nematodes by opening prey cuticles with movable teeth. Previous studies showed that exogenous serotonin treatment induces a predatory-like tooth movement and slower pharyngeal pumping in the absence of prey; however, physiological functions of serotonin during predation and other behaviors in P. pacificus remained completely unknown. Here, we investigate the roles of serotonin by generating mutations in Ppa-tph-1 and Ppa-bas-1, two key serotonin biosynthesis enzymes, and by genetic ablation of pharynx-associated serotonergic neurons. Mutations in Ppa-tph-1 reduced the pharyngeal pumping rate during bacterial feeding compared with wild-type. Moreover, the loss of serotonin or a subset of serotonergic neurons decreased the success of predation, but did not abolish the predatory feeding behavior completely. Detailed analysis using a high-speed camera revealed that the elimination of serotonin or the serotonergic neurons disrupted the timing and coordination of predatory tooth movement and pharyngeal pumping. This loss of synchrony significantly reduced the efficiency of successful predation events. These results suggest that serotonin has a conserved role in bacterial feeding and in addition drives the feeding rhythm of predatory behavior in Pristionchus.

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.

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.

Determination of Ancylostoma caninum ova viability using metabolic profiling

Differentiation between viable and non-viable hookworm ova in environmental samples is necessary in order to implement strategies to mitigate re-infections in endemic regions. In this study, an untargeted metabolic profiling method was developed that utilised gas chromatography-mass spectrometry (GC-MS) in order to investigate hookworm ova viability. Ancylostoma caninum was used to investigate the metabolites within viable and non-viable ova. Univariate and multivariate statistical analyses of the data resulted in the identification of 53 significant metabolites across all hookworm ova samples. The major compounds observed in viable and non-viable hookworm ova were tetradecanoic acid, commonly known as myristic acid [fold change (FC) = 0.4], and dodecanoic acid, commonly known as lauric acid (FC = 0.388). Additionally, the viable ova had self-protecting metabolites such as prostaglandins, a typical feature absent in non-viable ova. The results of this study demonstrate that metabolic profiling using GC-MS methods can be used to determine the viability of canine hookworm ova. Further studies are needed to assess the applicability of metabolic profiling using GC-MS to detect viable hookworm ova in the mixed (viable and non-viable) populations from environmental samples and identify the metabolites specific to human hookworm species.

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.

Haemonchus contortus: Genome Structure, Organization and Comparative Genomics

One of the first genome sequencing projects for a parasitic nematode was that for Haemonchus contortus. The open access data from the Wellcome Trust Sanger Institute provided a valuable early resource for the research community, particularly for the identification of specific genes and genetic markers. Later, a second sequencing project was initiated by the University of Melbourne, and the two draft genome sequences for H. contortus were published back-to-back in 2013. There is a pressing need for long-range genomic information for genetic mapping, population genetics and functional genomic studies, so we are continuing to improve the Wellcome Trust Sanger Institute assembly to provide a finished reference genome for H. contortus. This review describes this process, compares the H. contortus genome assemblies with draft genomes from other members of the strongylid group and discusses future directions for parasite genomics using the H. contortus model.

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.