A transcriptomic analysis of the phylum Nematoda

Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channels (GluCls), but it is unknown whether GluCl genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of-function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

A Minimal Kynurenine Pathway Was Preserved for Rhodoquinone but Not for De Novo NAD+ Biosynthesis in Parasitic Worms: The Essential Role of NAD+ Rescue Pathways

Aims: To determine the role of the kynurenine (KYN) pathway in rhodoquinone (RQ) and de novo NAD+ biosynthesis and whether NAD+ rescue pathways are essential in parasitic worms (helminths). Results: We demonstrate that RQ, the key electron transporter used by helminths under hypoxia, derives from the tryptophan (Trp) catabolism even in the presence of a minimal KYN pathway. We show that of the KYN pathway genes only the kynureninase and tryptophan/indoleamine dioxygenases are essential for RQ biosynthesis. Metabolic labeling with Trp revealed that the lack of the formamidase and kynurenine monooxygenase genes did not preclude RQ biosynthesis in the flatworm Mesocestoides corti. In contrast, a minimal KYN pathway prevented de novo NAD+ biosynthesis, as revealed by metabolic labeling in M. corti, which also lacks the 3-hydroxyanthranilate 3,4-dioxygenase gene. Our results indicate that most helminths depend solely on NAD+ rescue pathways, and some lineages rely exclusively on the nicotinamide salvage pathway. Importantly, the inhibition of the NAD+ recycling enzyme nicotinamide phosphoribosyltransferase with FK866 led cultured M. corti to death. Innovation: We use comparative genomics of more than 100 hundred helminth genomes, metabolic labeling, HPLC-mass spectrometry targeted metabolomics, and enzyme inhibitors to define pathways that lead to RQ and NAD+ biosynthesis in helminths. We identified the essential enzymes of these pathways in helminth lineages, revealing new potential pharmacological targets for helminthiasis. Conclusion: Our results demonstrate that a minimal KYN pathway was evolutionary maintained for RQ and not for de novo NAD+ biosynthesis in helminths and shed light on the essentiality of NAD+ rescue pathways in helminths.

rec-1 loss of function increases recombination in the central gene clusters at the expense of autosomal pairing centers

Meiotic control of crossover (CO) number and position is critical for homologous chromosome segregation and organismal fertility, recombination of parental genotypes, and the generation of novel genetic combinations. We here characterize the recombination rate landscape of a rec-1 loss of function modifier of CO position in Caenorhabditis elegans, one of the first ever modifiers discovered. By averaging CO position across hermaphrodite and male meioses and by genotyping 203 single-nucleotide variants covering about 95% of the genome, we find that the characteristic chromosomal arm-center recombination rate domain structure is lost in the loss of function rec-1 mutant. The rec-1 loss of function mutant smooths the recombination rate landscape but is insufficient to eliminate the nonuniform position of CO. Lower recombination rates in the rec-1 mutant are particularly found in the autosomal arm domains containing the pairing centers. We further find that the rec-1 mutant is of little consequence for organismal fertility and egg viability and thus for rates of autosomal nondisjunction. It nonetheless increases X chromosome nondisjunction rates and thus male appearance. Our findings question the maintenance of recombination rate heritability and genetic diversity among C. elegans natural populations, and they further suggest that manipulating genetic modifiers of CO position will help find quantitative trait loci located in low-recombining genomic regions normally refractory to discovery.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole and ivermectin are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channel (GluCl) genes, but it is unknown whether these genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required loss of two GluCl genes (avr-14 and avr-15) or loss of three GluCl genes (avr-14, avr-15, and glc-1). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Associate plant parasitic nematodes to weed species in some newly reclaimed lands

Plant-parasitic nematodes (PPNs) are vital soil organisms well-known to damage and reduce crop yield worldwide. Surveys were attempts to determine the impact of weed species on the communities and composition of nematodes in barley, wheat, quinoa, eggplant, and tomato crops in Alexandria and Ismailia regions of Egypt. During the surveys, eight occurring genera of nematodes were found namely; Meloidogyne spp, Pratylenchus spp, Helicotylenchus spp, Rotylenchulus spp, Xiphinema spp, Criconemoides spp, Ditylenchus spp, and Longidorus spp associated with the soil's rhizosphere of 28 weed species belonging to 12 families. Among these weeds, Hordeum marinum and Sonchus oleraceus were good hosts to nematode species. Both wheat and barley had higher nematode diversity than quinoa in the winter season. Pratylenchus spp, Meloidogyne spp and Rotylenchulus spp can be considered vital potential PPNs with economic importance. Nematode abundances and structural indices varied greatly based on the host weed species, crop types and soil characteristics. A positive correlation was monitored among weeds, nematode frequencies and relative abundances as well as their crops. Finally, weed species are critical components in nematode communities that may increase the incidence and severity of nematode risks based on crop type and soil characteristics. Therefore weeds should be managed properly to diminish reservoir sites when developing nematode management options.

daf-42 is an evolutionarily young gene essential for dauer development in Caenorhabditis elegans

Under adverse environmental conditions, nematodes arrest into dauer, an alternative developmental stage for diapause. Dauer endures unfavorable environments and interacts with host animals to access favorable environments, thus playing a critical role in survival. Here, we report that in Caenorhabditis elegans, daf-42 is essential for development into the dauer stage, as the null mutant of daf-42 exhibited a "no viable dauer" phenotype in which no viable dauers were obtained in any dauer-inducing conditions. Long-term time lapse microscopy of synchronized larvae revealed that daf-42 is involved in developmental changes from the pre-dauer L2d stage to the dauer stage. daf-42 encodes large, disordered proteins of various sizes that are expressed in and secreted from the seam cells within a narrow time window shortly before the molt into dauer stage. Transcriptome analysis showed that the transcription of genes involved in larval physiology and dauer metabolism is highly affected by the daf-42 mutation. Contrary to the notion that essential genes that control the life and death of an organism may be well conserved across diverse species, daf-42 is an evolutionarily young gene conserved only in the Caenorhabditis genus. Our study shows that dauer formation is a vital process that is controlled not only by conserved genes but also by newly emerged genes, providing important insights into evolutionary mechanisms.

An efficient method for viable cryopreservation and recovery of hookworms and other gastrointestinal nematodes in the laboratory

Hookworms (genera Ancylostoma and Necator) are amongst the most prevalent and important parasites of humans globally. These intestinal parasites ingest blood, resulting in anemia, growth stunting, malnutrition, and adverse pregnancy outcomes. They are also critical parasites of dogs and other animals. In addition, hookworms and hookworm products are being explored for their use in treatment of autoimmune and inflammatory diseases. There is thus a significant and growing interest in these mammalian host-obligate parasites. Laboratory research is hampered by the lack of good means of cryopreservation and recovery of parasites. Here, we describe a robust method for long-term (≥3 year) cryopreservation and recovery of both Ancylostoma and Necator hookworms that is also applicable to two other intestinal parasites that passage through the infective L3 stage, Strongyloides ratti and Heligmosomoides polygyrus bakeri. The key is a revised recovery method, in which cryopreserved L1s are thawed and raised to the infective L3 stage using activated charcoal mixed with uninfected feces from a permissive host. This technique will greatly facilitate research on and availability of gastrointestinal parasitic nematodes with great importance to global health, companion animal health, and autoimmune/inflammatory disease therapies.

An efficient method for viable cryopreservation of hookworms and other gastrointestinal nematodes in the laboratory

Hookworms (genera Ancylostoma and Necator ) are amongst of the most prevalent and important parasites of humans globally. These intestinal parasites ingest blood, resulting in anemia, growth stunting, malnutrition, and adverse pregnancy outcomes. They are also critical parasites of dogs and other animals. In addition, hookworms and hookworm products are being explored for their use in treatment of autoimmune and inflammatory diseases. There is thus a significant and growing interest in these mammalian host-obligate parasites. Laboratory research is hampered by the lack of good means of cryopreservation. Here, we describe a robust method for long-term (≥3 year) cryoprotection and recovery of both Ancylostoma and Necator hookworms that is also applicable to two other intestinal parasites that passages through the infective third larval stage, Strongyloides ratti and H eligmosomoides polygyrus bakeri . The key is the use cryo-preserved first-staged larvae raised to the infective third larval stage using activated charcoal mixed with uninfected feces from a permissive host. This technique will greatly facilitate research on and availability of gastrointestinal parasitic nematodes with great importance to global health, companion animal health, and autoimmune/inflammatory disease therapies.

Transcriptomic and Proteomic Analysis of Marine Nematode Litoditis marina Acclimated to Different Salinities

Salinity is a critical abiotic factor for all living organisms. The ability to adapt to different salinity environments determines an organism’s survival and ecological niches. Litoditis marina is a euryhaline marine nematode widely distributed in coastal ecosystems all over the world, although numerous genes involved in its salinity response have been reported, the adaptive mechanisms underlying its euryhalinity remain unexplored. Here, we utilized worms which have been acclimated to either low-salinity or high-salinity conditions and evaluated their basal gene expression at both transcriptomic and proteomic levels. We found that several conserved regulators, including osmolytes biosynthesis genes, transthyretin-like family genes, V-type H⁺-transporting ATPase and potassium channel genes, were involved in both short-term salinity stress response and long-term acclimation processes. In addition, we identified genes related to cell volume regulation, such as actin regulatory genes, Rho family small GTPases and diverse ion transporters, which might contribute to hyposaline acclimation, while the glycerol biosynthesis genes gpdh-1 and gpdh-2 accompanied hypersaline acclimation in L. marina. This study paves the way for further in-depth exploration of the adaptive mechanisms underlying euryhalinity and may also contribute to the study of healthy ecosystems in the context of global climate change.

Advances in our understanding of nematode ion channels as potential anthelmintic targets

Ion channels are specialized multimeric proteins that underlie cell excitability. These channels integrate with a variety of neuromuscular and biological functions. In nematodes, the physiological behaviors including locomotion, navigation, feeding and reproduction, are regulated by these protein entities. Majority of the antinematodal chemotherapeutics target the ion channels to disrupt essential biological functions. Here, we have summarized current advances in our understanding of nematode ion channel pharmacology. We review cys-loop ligand gated ion channels (LGICs), including nicotinic acetylcholine receptors (nAChRs), acetylcholine-chloride gated ion channels (ACCs), glutamate-gated chloride channels (GluCls), and GABA (γ-aminobutyric acid) receptors, and other ionotropic receptors (transient receptor potential (TRP) channels and potassium ion channels). We have provided an update on the pharmacological properties of these channels from various nematodes. This article catalogs the differences in ion channel composition and resulting pharmacology in the phylum Nematoda. This diversity in ion channel subunit repertoire and pharmacology emphasizes the importance of pursuing species-specific drug target research. In this review, we have provided an overview of recent advances in techniques and functional assays available for screening ion channel properties and their application.

Phylogenomic Analysis of the Phylum Nematoda: Conflicts and Congruences With Morphology, 18S rRNA, and Mitogenomes

Phylogenetic relationships within many lineages of the phylum Nematoda remain unresolved, despite numerous morphology-based and molecular analyses. We performed several phylogenomic analyses using 286 published genomes and transcriptomes and 19 new transcriptomes by focusing on Trichinellida, Spirurina, Rhabditina, and Tylenchina separately, and by analyzing a selection of species from the whole phylum Nematoda. The phylogeny of Trichinellida supported the division of Trichinella into encapsulated and non-encapsulated species and placed them as sister to Trichuris. The Spirurina subtree supported the clades formed by species from Ascaridomorpha and Spiruromorpha respectively, but did not support Dracunculoidea. The analysis of Tylenchina supported a clade that included all sampled species from Tylenchomorpha and placed it as sister to clades that included sampled species from Cephalobomorpha and Panagrolaimomorpha, supporting the hypothesis that postulates the single origin of the stomatostylet. The Rhabditina subtree placed a clade composed of all sampled species from Diplogastridae as sister to a lineage consisting of paraphyletic Rhabditidae, a single representative of Heterorhabditidae and a clade composed of sampled species belonging to Strongylida. It also strongly supported all suborders within Strongylida. In the phylum-wide analysis, a clade composed of all sampled species belonging to Enoplia were consistently placed as sister to Dorylaimia + Chromadoria. The topology of the Nematoda backbone was consistent with previous studies, including polyphyletic placement of sampled representatives of Monhysterida and Araeolaimida.

Haemonchus contortus Transthyretin-Like Protein TTR-31 Plays Roles in Post-Embryonic Larval Development and Potentially Apoptosis of Germ Cells

Transthyretin (TTR)-like proteins play multi-function roles in nematode and are important component of excretory/secretory product in Haemonchus contortus. In this study, we functionally characterised a secretory transthyretin-like protein in the barber's pole worm H. contortus. A full-length of transthyretin-like protein-coding gene (Hc-ttr-31) was identified in this parasitic nematode, representing a counterpart of Ce-ttr-31 in Caenorhabditis elegans. High transcriptional levels of Hc-ttr-31 were detected in the egg and early larval stages of H. contortus, with the lowest level measured in the adult stage, indicating a decreased transcriptional pattern of this gene during nematode development. Localisation analysis indicated a secretion of TTR-31 from the intestine to the gonad, suggesting additional roles of Hc-ttr-31 in nematode reproduction. Expression of Hc-ttr-31 and Ce-ttr-31 in C. elegans did not show marked influence on the nematode development and reproduction, whereas Hc-ttr-31 RNA interference-mediated gene knockdown of Ce-ttr-31 shortened the lifespan, decreased the brood size, slowed the pumping rate and inhibited the growth of treated worms. Particularly, gene knockdown of Hc-ttr-31 in C. elegans was linked to activated apoptosis signalling pathway, increased general reactive oxygen species (ROS) level, apoptotic germ cells and facultative vivipary phenotype, as well as suppressed germ cell removal signalling pathways. Taken together, Hc-ttr-31 appears to play roles in regulating post-embryonic larval development, and potentially in protecting gonad from oxidative stress and mediating engulfment of apoptotic germ cells. A better knowledge of these aspects should contribute to a better understanding of the developmental biology of H. contortus and a discovery of potential targets against this and related parasitic worms.

Co-evolution based machine-learning for predicting functional interactions between human genes

Over the next decade, more than a million eukaryotic species are expected to be fully sequenced. This has the potential to improve our understanding of genotype and phenotype crosstalk, gene function and interactions, and answer evolutionary questions. Here, we develop a machine-learning approach for utilizing phylogenetic profiles across 1154 eukaryotic species. This method integrates co-evolution across eukaryotic clades to predict functional interactions between human genes and the context for these interactions. We benchmark our approach showing a 14% performance increase (auROC) compared to previous methods. Using this approach, we predict functional annotations for less studied genes. We focus on DNA repair and verify that 9 of the top 50 predicted genes have been identified elsewhere, with others previously prioritized by high-throughput screens. Overall, our approach enables better annotation of function and functional interactions and facilitates the understanding of evolutionary processes underlying co-evolution. The manuscript is accompanied by a webserver available at: https://mlpp.cs.huji.ac.il.

With the rise in number of eukaryotic species being fully sequenced, large scale phylogenetic profiling can give insights on gene function, Here, the authors describe a machine-learning approach that integrates co-evolution across eukaryotic clades to predict gene function and functional interactions among human genes.

Comparative evaluation of different molecular methods for DNA extraction from individual Teladorsagia circumcincta nematodes

Background

The purpose of this study was to develop a reliable DNA extraction protocol to use on individual Teladorsagia circumcincta nematode specimens to produce high quality DNA for genome sequencing and phylogenetic analysis. Pooled samples have been critical in providing the groundwork for T. circumcincta genome construction, but there is currently no standard method for extracting high-quality DNA from individual nematodes. 11 extraction kits were compared based on DNA quality, yield, and processing time.

Results

11 extraction protocols were compared, and the concentration and purity of the extracted DNA was quantified. Median DNA concentration among all methods measured on NanoDrop 2000™ ranged between 0.45–11.5 ng/μL, and on Qubit™ ranged between undetectable – 0.962 ng/μL. Median A260/280 ranged between 0.505–3.925, and median A260/230 ranged − 0.005 – 1.545. Larval exsheathment to remove the nematode cuticle negatively impacted DNA concentration and purity.

Conclusions

A Schistosoma sp. DNA extraction method was determined as most suitable for individual T. circumcincta nematode specimens due to its resulting DNA concentration, purity, and relatively fast processing time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00695-6.

Genome-Wide Transcriptional Responses of Marine Nematode Litoditis marina to Hyposaline and Hypersaline Stresses

Maintenance of osmotic homeostasis is essential for all organisms, especially for marine animals in the ocean with 3% salinity or higher. However, the underlying molecular mechanisms that how marine animals adapt to high salinity environment compared to their terrestrial relatives, remain elusive. Here, we investigated marine animal’s genome-wide transcriptional responses to salinity stresses using an emerging marine nematode model Litoditis marina. We found that the transthyretin-like family genes were significantly increased in both hyposaline and hypersaline conditions, while multiple neurotransmitter receptor and ion transporter genes were down-regulated in both conditions, suggesting the existence of conserved strategies for response to stressful salinity environments in L. marina. Unsaturated fatty acids biosynthesis related genes, neuronal related tubulins and intraflagellar transport genes were specifically up-regulated in hyposaline treated worms. By contrast, cuticle related collagen genes were enriched and up-regulated for hypersaline response. Given a wide range of salinity tolerance of the marine nematodes, this study and further genetic analysis of key gene(s) of osmoregulation in L. marina will likely provide important insights into biological evolution and environmental adaptation mechanisms in nematodes and other invertebrate animals in general.

The production of excretory-secretory molecules from Heligmosomoides polygyrus bakeri fourth stage larvae varies between mixed and single sex cultures

Background

Excretory-secretory (ES) products are crucial in maintaining helminths in the host. Consequently, the proteins of ES are potential vaccine molecules and potential therapeutic agents for autoimmune diseases. Heligmosomoides polygyrus bakeri, a gastrointestinal parasite of mice, is a model of hookworm infection in humans. ES produced by both sexes of H. polygyrus bakeri L4 stage cultured separately shows different immunomodulatory properties than ES obtained when both sexes are cultured together. Accordingly, the objective of this study was to identify and compare the excretory-secretory molecules from single-sex and mixed cultures.

Methods

The composition of ES of male and female L4 stage nematodes in the presence (cultured together) or absence (cultured alone) of the opposite sex was examined. Proteins were identified using mass spectrometry. The functions of identified proteins were explored with Blast2GO.

Results

A total of 258 proteins derived from mixed larval culture in the presence of sex pheromones were identified, 160 proteins from pure female cultures and 172 from pure male cultures. Exposure of nematodes to the sex pheromones results in abundant production of proteins with immunomodulatory properties such as Val proteins, acetylcholinesterases, TGF-β mimic 9 and HpARI. Proteins found only in ES from mixed larval cultures were TGF-β mimics 6 and 7 as well as galectin.

Conclusions

The presence of the opposite sex strongly influences the composition of ES products, probably by chemical (pheromone) communication between individuals. However, examination of the composition of ES from various conditions gives an opportunity for searching for new potentially therapeutic compounds and anthelminthics as well as components of vaccines. Manipulation of the nematode environment might be important for the studies on the immunomodulatory potential of nematodes.

Evolution of Developmental GATA Factors in Nematodes

GATA transcription factors are found in animals, plants, and fungi. In animals, they have important developmental roles in controlling specification of cell identities and executing tissue-specific differentiation. The Phylum Nematoda is a diverse group of vermiform animals that inhabit ecological niches all over the world. Both free-living and parasitic species are known, including those that cause human infectious disease. To date, GATA factors in nematodes have been studied almost exclusively in the model system C. elegans and its close relatives. In this study, we use newly available sequences to identify GATA factors across the nematode phylum. We find that most species have fewer than six GATA factors, but some species have 10 or more. Comparisons of gene and protein structure suggest that there were at most two GATA factors at the base of the phylum, which expanded by duplication and modification to result in a core set of four factors. The high degree of structural similarity with the corresponding orthologues in C. elegans suggests that the nematode GATA factors share similar functions in development.

Cathepsin F of Teladorsagia circumcincta is a recently evolved cysteine protease

Parasitic cysteine proteases are involved in parasite stage transition, invasion of host tissues, nutrient uptake, and immune evasion. The cysteine protease cathepsin F is the most abundant protein produced by fourth-stage larvae (L4) of the nematode Teladorsagia circumcincta, while its transcript is only detectable in L4 and adults. T. circumcincta cathepsin F is a recently evolved cysteine protease that does not fall clearly into either of the cathepsin L or F subfamilies. This protein exhibits characteristics of both cathepsins F and L, and its phylogenetic relationship to its closest homologs is distant, including proteins of closely related nematodes of the same subfamily.

How do Caenorhabditis elegans worms survive in highly viscous habitats?

The nematode Caenorhabditis elegans is a filter feeder that lives in various viscous habitats such as soil, the intestines of slugs, and rotting materials such as fruits and stems. Caenorhabditis elegans draws in suspensions of bacteria and separates bacteria from water using the pharyngeal pump. Although these worms often live in highly viscous habitats, it is still unclear how they survive in these environments by eating bacteria. In this study, we investigated the effects of suspension viscosity on the survival rate of malnourished worms by combining live imaging and scaling analyses. We found that survival rate decreased with increases in viscosity because the high viscosity suppressed the amount of food ingested. The same tendency was found in two feeding-defective mutants, eat-6(ad467) and eat-6(ad997). We also found that the high viscosity weakened pump function, but the velocities in the pharynx were not zero, even in the most viscous suspensions. Finally, we estimated the amount of ingested food using scaling analyses, which provided a master curve of the experimental survival rates. These results illustrate that the survival rate of C. elegans worms is strongly dependent on the ingested bacteria per unit time associated with physical environments, such as the viscosity of food suspensions and the cell density of bacteria. The pump function of the C. elegans pharynx is not completely lost even in fluids that have 105 times higher viscosity than water, which may contribute to their ability to survive around the world in highly viscous environments.

Characterization of the Biogenic Volatile Organic Compounds (BVOCs) and Analysis of the PR1 Molecular Marker in Vitis vinifera L. Inoculated with the Nematode Xiphinema index

Upon pathogen attack, plants very quickly undergo rather complex physico-chemical changes, such as the production of new chemicals or alterations in membrane and cell wall properties, to reduce disease damages. An underestimated threat is represented by root parasitic nematodes. In Vitis vinifera L., the nematode Xiphinema index is the unique vector of Grapevine fanleaf virus, responsible for fanleaf degeneration, one of the most widespread and economically damaging diseases worldwide. The aim of this study was to investigate changes in the emission of biogenic volatile organic compounds (BVOCs) in grapevines attacked by X. index. BVOCs play a role in plant defensive mechanisms and are synthetized in response to biotic damages. In our study, the BVOC profile was altered by the nematode feeding process. We found a decrease in β-ocimene and limonene monoterpene emissions, as well as an increase in α-farnesene and α-bergamotene sesquiterpene emissions in nematode-treated plants. Moreover, we evaluated the PR1 gene expression. The transcript level of PR1 gene was higher in the nematode-wounded roots, while in the leaf tissues it showed a lower expression compared to control grapevines.

Anthelmintic and metabolomic analyses of chicory (Cichorium intybus) identify an industrial by-product with potent in vitro antinematodal activity

Chicory (Cichorium intybus) is a bioactive forage rich in sesquiterpene lactones (SLs) with reported in vitro and in vivo anthelmintic activity in livestock. However, the on-farm adoption of chicory as an anthelmintic crop is limited and may be facilitated by using standardised industrial chicory material. Chicory root pulp is a by-product obtained from industrial chicory roots after inulin extraction and can potentially retain SLs. However, SL content and associated anthelmintic activity of chicory root pulp have not been investigated. Here, we evaluated the anthelmintic activity of SL-enriched extracts from chicory root pulp and forage chicory, and used untargeted metabolomics and molecular networking to identify potential anthelmintic molecules. Six different sources of chicory material were used: fresh chicory root pulp (from industrial chicory roots C. intybus var. sativum; "Root Pulp"), fresh leaves from chicory cv. Spadona (sampled on four occasions) and fresh leaves from chicory cv. Choice. The resulting extracts were tested for anthelmintic activity against the free-living nematode Caenorhabditis elegans and the pig nematode Ascaris suum. The cytotoxicity of the chicory extracts was evaluated on mammalian (Vero) cells. In the C. elegans assays, the Root Pulp was the most potent extract and induced paralysis in >95% of worms exposed to >250  μg extract/mL (EC50 = 64.2 μg/mL). In the A. suum assays, the Root Pulp was also the most potent chicory extract to inhibit worm motility (EC50 = 87.6  μg/mL), followed closely by two of the Spadona leaf extracts (EC50 = 89.8  μg/mL and 112.2  μg/mL) The Root Pulp extract had the lowest cytotoxicity of all tested extracts towards mammalian cells, with a selectivity index of 5.37. Untargeted metabolomics revealed that chicory Root Pulp had a markedly different chemical profile in comparison with forage chicory extracts. Molecular networking confirmed several SLs and SL-derivatives mainly present in chicory root pulp, that may be responsible of its potent anti-parasitic activity. Bioactivity-based molecular networking of chicory root pulp and the most potent forage chicory extracts revealed a high predicted anthelmintic score for the guaianolide SL 11,13-dihydro-lactucopicrin. In conclusion, chicory root pulp showed potent and selective in vitro anthelmintic activity against C. elegans and A. suum, with low cytotoxicity in mammalian cells. The promising anthelmintic activity of chicory root pulp should be confirmed in vivo to further explore the potential of this agro-industrial by-product as a nutraceutical anthelmintic for livestock and as novel source of anti-parasitic compounds.

Drug discovery technologies: Caenorhabditis elegans as a model for anthelmintic therapeutics

Helminthiasis is one of the gravest problems worldwide. There is a growing concern on less available anthelmintics and the emergence of resistance creating a major threat to human and livestock health resources. Novel and broad-spectrum anthelmintics are urgently needed. The free-living nematode Caenorhabditis elegans could address this issue through automated high-throughput technologies for the screening of large chemical libraries. This review discusses the strong advantages and limitations for using C elegans as a screening method for anthelmintic drug discovery. C elegans is the best model available for the validation of novel effective drugs in treating most, if not all, helminth infections, and for the elucidation the mode of action of anthelmintic candidates. This review also focuses on available technologies in the discovery of anthelmintics published over the last 15 years with particular attention to high-throughput technologies over conventional screens. On the other hand, this review highlights how combinatorial and nanomedicine strategies could prolong the use of anthelmintics and control resistance problems.

Vital roles of PCNA K165 modification during C. elegans gametogenesis and embryogenesis

In eukaryotes, the DNA damage bypass pathway is promoted by ubiquitylation of PCNA at the conserved lysine 164. Using CRISPR-Cas9 system, we introduced amino acid substitution at K165 of C. elegans PCNA that corresponds to K164 in other characterised organisms and examined the contribution of this residue at a variety of stages during development. In the presence of UV-induced DNA lesions, PCNA-K165 is crucial for not only the early embryonic stages but also during larval development, implicating its functions for a broad time period during animal development. We also show that, without induction of DNA damage, concomitant inhibition of PCNA ubiquitylation and checkpoint activation causes abnormal gametogenesis events and severely impairs reproduction of worms. Our findings suggest a conserved function of PCNA ubiquitylation in tolerance of UV-induced damage and also propose that PCNA ubiquitylation contributes to gametogenesis during unperturbed C. elegans development.

Improved phylogenomic sampling of free-living nematodes enhances resolution of higher-level nematode phylogeny

Background

Nematodes are among the most diverse and abundant metazoans on Earth, but research on them has been biased toward parasitic taxa and model organisms. Free-living nematodes, particularly from the clades Enoplia and Dorylaimia, have been underrepresented in genome-scale phylogenetic analyses to date, leading to poor resolution of deep relationships within the phylum.

Results

We supplemented publicly available data by sequencing transcriptomes of nine free-living nematodes and two important outgroups and conducted a phylum-wide phylogenomic analysis including a total of 108 nematodes. Analysis of a dataset generated using a conservative orthology inference strategy resulted in a matrix with a high proportion of missing data and moderate to weak support for branching within and placement of Enoplia. A less conservative orthology inference approach recovered more genes and resulted in higher support for the deepest splits within Nematoda, recovering Enoplia as the sister taxon to the rest of Nematoda. Relationships within major clades were similar to those found in previously published studies based on 18S rDNA.

Conclusions

Expanded transcriptome sequencing of free-living nematodes has contributed to better resolution among deep nematode lineages, though the dataset is still strongly biased toward parasites. Inclusion of more free-living nematodes in future phylogenomic analyses will allow a clearer understanding of many interesting aspects of nematode evolution, such as morphological and molecular adaptations to parasitism and whether nematodes originated in a marine or terrestrial environment.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1444-x) contains supplementary material, which is available to authorized users.

Haem Biology in Metazoan Parasites - 'The Bright Side of Haem'

Traditionally, host haem has been recognized as a cytotoxic molecule that parasites need to eliminate or detoxify in order to survive. However, recent evidence indicates that some lineages of parasites have lost genes that encode enzymes involved specifically in endogenous haem biosynthesis. Such lineages thus need to acquire and utilize haem originating from their host animal, making it an indispensable molecule for their survival and reproduction. In multicellular parasites, host haem needs to be systemically distributed throughout their bodies to meet the haem demands in all cell and tissue types. Host haem also gets deposited in parasite eggs, enabling embryogenesis and reproduction. Clearly, a better understanding of haem biology in multicellular parasites should elucidate organismal adaptations to obligatory blood-feeding.

Comprehensive Evolutionary Analysis of the Major RNA-Induced Silencing Complex Members

RNA-induced silencing complex (RISC) plays a critical role in small interfering RNA (siRNA) and microRNAs (miRNA) pathways. Accumulating evidence has demonstrated that the major RISC members (AGO, DICER, TRBP, PACT and GW182) represent expression discrepancies or multiple orthologues/paralogues in different species. To elucidate their evolutionary characteristics, an integrated evolutionary analysis was performed. Here, animal and plant AGOs were divided into three classes (multifunctional AGOs, siRNA-associated AGOs and piRNA-associated AGOs for animal AGOs and multifunctional AGOs, siRNA-associated AGOs and complementary functioning AGOs for plant AGOs). Animal and plant DICERs were grouped into one class (multifunctional DICERs) and two classes (multifunctional DICERs and siRNA-associated DICERs), respectively. Protista/fungi AGOs or DICERs were specifically associated with the siRNA pathway. Additionally, TRBP/PACT/GW182 were identified only in animals, and all of them functioned in the miRNA pathway. Mammalian AGOs, animal DICERs and chordate TRBP/PACT were found to be monophyletic. A large number of gene duplications were identified in AGO and DICER groups. Taken together, we provide a comprehensive evolutionary analysis, describe a phylogenetic tree-based classification of the major RISC members and quantify their gene duplication events. These findings are potentially useful for classifying RISCs, optimizing species-specific RISCs and developing research model organisms.

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

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

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

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

Cysteine proteases as digestive enzymes in parasitic helminths

We briefly review cysteine proteases (orthologs of mammalian cathepsins B, L, F, and C) that are expressed in flatworm and nematode parasites. Emphasis is placed on enzyme activities that have been functionally characterized, are associated with the parasite gut, and putatively contribute to degrading host proteins to absorbable nutrients [1–4]. Often, gut proteases are expressed as multigene families, as is the case with Fasciola [5] and Haemonchus [6], presumably expanding the range of substrates that can be degraded, not least during parasite migration through host tissues [5]. The application of the free-living planarian and Caenorhabditis elegans as investigative models for parasite cysteine proteases is discussed. Finally, because of their central nutritive contribution, targeting the component gut proteases with small-molecule chemical inhibitors and understanding their utility as vaccine candidates are active areas of research [7].

Genomic Identification and Functional Characterization of Essential Genes in Caenorhabditis elegans

Using combined genetic mapping, Illumina sequencing, bioinformatics analyses, and experimental validation, we identified 60 essential genes from 104 lethal mutations in two genomic regions of Caenorhabditis elegans totaling ∼14 Mb on chromosome III(mid) and chromosome V(left). Five of the 60 genes had not previously been shown to have lethal phenotypes by RNA interference depletion. By analyzing the regions around the lethal missense mutations, we identified four putative new protein functional domains. Furthermore, functional characterization of the identified essential genes shows that most are enzymes, including helicases, tRNA synthetases, and kinases in addition to ribosomal proteins. Gene Ontology analysis indicated that essential genes often encode for enzymes that conduct nucleic acid binding activities during fundamental processes, such as intracellular DNA replication, transcription, and translation. Analysis of essential gene shows that they have fewer paralogs, encode proteins that are in protein interaction hubs, and are highly expressed relative to nonessential genes. All these essential gene traits in C. elegans are consistent with those of human disease genes. Most human orthologs (90%) of the essential genes in this study are related to human diseases. Therefore, functional characterization of essential genes underlines their importance as proxies for understanding the biological functions of human disease genes.

The genomic basis of nematode parasitism

Nematodes are highly abundant animals, and many species have a parasitic lifestyle. Nematode parasites are important pathogens of humans and other animals, and there is considerable interest in understanding their molecular and genomic adaptations to nematode parasitism. This has been approached in three main ways: comparing the genomes of closely related parasitic and free-living taxa, comparing the gene expression of parasitic and free-living life cycle stages of parasitic nematode species, and analysing the molecules that parasitic nematodes excrete and secrete. To date, these studies show that many species of parasitic nematodes have genomes that have large gene families coding for proteases/peptidases, protease inhibitors, SCP/TAPS proteins and acetylcholinesterases, and in many cases there is evidence that these appear to be used by parasitic stages inside hosts, and are often secreted. Many parasitic nematodes have taxa-restricted gene families that also appear to be involved in parasitism, emphasizing that there is still much to be discovered about what it takes to be a parasitic nematode.

Caenorhabditis elegans: nature and nurture gift to nematode parasitologists

The free-living nematode Caenorhabditis elegans is the simplest animal model organism to work with. Substantial knowledge and tools have accumulated over 50 years of C. elegans research. The use of C. elegans relating to parasitic nematodes from a basic biology standpoint or an applied perspective has increased in recent years. The wealth of information gained on the model organism, the use of the powerful approaches and technologies that have advanced C. elegans research to parasitic nematodes and the enormous success of the omics fields have contributed to bridge the divide between C. elegans and parasite nematode researchers. We review key fields, such as genomics, drug discovery and genetics, where C. elegans and nematode parasite research have convened. We advocate the use of C. elegans as a model to study helminth metabolism, a neglected area ready to advance. How emerging technologies being used in C. elegans can pave the way for parasitic nematode research is discussed.

Molecular expression and characterization of GCP7 gene of Haemonchus contortus

Haemonchus contortus is a highly pathogenic and most economically important parasite of sheep and goats worldwide. The cysteine proteases from H. contortus are prime targets for vaccine development. In the present communication we report the molecular expression and characterization of cathepsin B-like cysteine protease, GCP7 gene of H. contortus to study its efficiency as target protein for immunoprophylaxis against haemonchosis in sheep and goats. The complete ORF of GCP7 gene, devoid of the signal sequence, was amplified by RT-PCR from mRNA isolated from H. contortus and was cloned initially into pTZ57R/T cloning vector and then sub-cloned in the pET32a(+) expression vector to produce GCP7 antigen. The nucleotide and deduced amino acid sequence of the GCP7 was aligned against the related sequences of H. contortus available in public domain for in silico analysis by DNA STAR and MEGA version 4.0 softwares. The nucleotide sequence revealed that the GCP7 gene of H.contortus (Indian isolate) encodes 324 amino acids (devoid of signal sequence) and its nucleotide sequence had 95.9% to 99.4% sequence homology with that of U.S.A. and previously published Indian isolates. A high level expression of recombinant (r) GCP7 protein was observed in the molecular range (Mr) of 55 kDa. The rGCP7 protein was confirmed by its specific immunoreactivity against known reference positive sheep sera.

Signatures of co-evolutionary host-pathogen interactions in the genome of the entomopathogenic nematode Steinernema carpocapsae

Background

The entomopathogenic nematode Steinernema carpocapsae has been used worldwide as a biocontrol agent for insect pests, making it an interesting model for understanding parasite-host interactions. Two models propose that these interactions are co-evolutionary processes in such a way that equilibrium is never reached. In one model, known as “arms race”, new alleles in relevant genes are fixed in both host and pathogens by directional positive selection, producing recurrent and alternating selective sweeps. In the other model, known as“trench warfare”, persistent dynamic fluctuations in allele frequencies are sustained by balancing selection. There are some examples of genes evolving according to both models, however, it is not clear to what extent these interactions might alter genome-level evolutionary patterns and intraspecific diversity. Here we investigate some of these aspects by studying genomic variation in S. carpocapsae and other pathogenic and free-living nematodes from phylogenetic clades IV and V.

Results

To look for signatures of an arms-race dynamic, we conducted massive scans to detect directional positive selection in interspecific data. In free-living nematodes, we detected a significantly higher proportion of genes with sites under positive selection than in parasitic nematodes. However, in these genes, we found more enriched Gene Ontology terms in parasites. To detect possible effects of dynamic polymorphisms interactions we looked for signatures of balancing selection in intraspecific genomic data. The observed distribution of Tajima’s D values in S. carpocapsae was more skewed to positive values and significantly different from the observed distribution in the free-living Caenorhabditis briggsae. Also, the proportion of significant positive values of Tajima’s D was elevated in genes that were differentially expressed after induction with insect tissues as compared to both non-differentially expressed genes and the global scan.

Conclusions

Our study provides a first portrait of the effects that lifestyle might have in shaping the patterns of selection at the genomic level. An arms-race between hosts and pathogens seems to be affecting specific genetic functions but not necessarily increasing the number of positively selected genes. Trench warfare dynamics seem to be acting more generally in the genome, likely focusing on genes responding to the interaction, rather than targeting specific genetic functions.

Electronic supplementary material

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

Entomopathogenic nematodes in agricultural areas in Brazil

Entomopathogenic nematodes (EPNs) (Steinernematidae and Heterorhabditidae) can control pests due to the mutualistic association with bacteria that kill the host by septicemia and make the environment favorable for EPNs development and reproduction. The diversity of EPNs in Brazilian soils requires further study. The identification of EPNs, adapted to environmental and climatic conditions of cultivated areas is important for sustainable pest suppression in integrated management programs in agricultural areas of Brazil. The objective was to identify EPNs isolated from agricultural soils with annual, fruit and forest crops in Brazil. Soil samples were collected and stored in 250 ml glass vials. The nematodes were isolated from these samples with live bait traps ([Galleria mellonella L. (Lepidoptera: Pyralidae) larvae]. Infective juveniles were collected with White traps and identified by DNA barcoding procedures by sequencing the D2/D3 expansion of the 28S rDNA region by PCR. EPNs identified in agricultural areas in Brazil were Heterorhabditis amazonensis, Metarhabditis rainai, Oscheios tipulae and Steinernema rarum. These species should be considered pest biocontrol agents in Brazilian agricultural areas.

Recent Duplication and Functional Divergence in Parasitic Nematode Levamisole-Sensitive Acetylcholine Receptors

Helminth parasites rely on fast-synaptic transmission in their neuromusculature to experience the outside world and respond to it. Acetylcholine plays a pivotal role in this and its receptors are targeted by a wide variety of both natural and synthetic compounds used in human health and for the control of parasitic disease. The model, Caenorhabditis elegans is characterized by a large number of acetylcholine receptor subunit genes, a feature shared across the nematodes. This dynamic family is characterized by both gene duplication and loss between species. The pentameric levamisole-sensitive acetylcholine receptor has been characterized from C. elegans, comprised of five different subunits. More recently, cognate receptors have been reconstituted from multiple parasitic nematodes that are found to vary in subunit composition. In order to understand the implications of receptor composition change and the origins of potentially novel drug targets, we investigated a specific example of subunit duplication based on analysis of genome data for 25 species from the 50 helminth genome initiative. We found multiple independent duplications of the unc-29, acetylcholine receptor subunit, where codon substitution rate analysis identified positive, directional selection acting on amino acid positions associated with subunit assembly. Characterization of four gene copies from a model parasitic nematode, Haemonchus contortus, demonstrated that each copy has acquired unique functional characteristics based on phenotype rescue of transgenic C. elegans and electrophysiology of receptors reconstituted in Xenopus oocytes. We found evidence that a specific incompatibility has evolved for two subunits co-expressed in muscle. We demonstrated that functional divergence of acetylcholine receptors, driven by directional selection, can occur more rapidly than previously thought and may be mediated by alteration of receptor assembly. This phenomenon is common among the clade V parasitic nematodes and this work provides a foundation for understanding the broader context of changing anthelmintic drug targets across the parasitic nematodes.

Genomes of parasitic nematodes (Meloidogyne hapla, Meloidogyne incognita, Ascaris suum and Brugia malayi) have a reduced complement of small RNA interference pathway genes: knockdown can reduce host infectivity of M. incognita

The discovery of RNA interference (RNAi) as an endogenous mechanism of gene regulation in a range of eukaryotes has resulted in its extensive use as a tool for functional genomic studies. It is important to study the mechanisms which underlie this phenomenon in different organisms, and in particular to understand details of the effectors that modulate its effectiveness. The aim of this study was to identify and compare genomic sequences encoding genes involved in the RNAi pathway of four parasitic nematodes: the plant parasites Meloidogyne hapla and Meloidogyne incognita and the animal parasites Ascaris suum and Brugia malayi because full genomic sequences were available-in relation to those of the model nematode Caenorhabditis elegans. The data generated was then used to identify some potential targets for control of the root knot nematode, M. incognita. Of the 84 RNAi pathway genes of C. elegans used as model in this study, there was a 42-53 % reduction in the number of effectors in the parasitic nematodes indicating substantial differences in the pathway between species. A gene each from six functional groups of the RNAi pathway of M. incognita was downregulated using in vitro RNAi, and depending on the gene (drh-3, tsn-1, rrf-1, xrn-2, mut-2 and alg-1), subsequent plant infection was reduced by up to 44 % and knockdown of some genes (i.e. drh-3, mut-2) also resulted in abnormal nematode development. The information generated here will contribute to defining targets for more robust nematode control using the RNAi technology.

Understanding Haemonchus contortus Better Through Genomics and Transcriptomics

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.

Early origin and adaptive evolution of the GW182 protein family, the key component of RNA silencing in animals

The GW182 proteins are a key component of the miRNA-dependent post-transcriptional silencing pathway in animals. They function as scaffold proteins to mediate the interaction of Argonaute (AGO)-containing complexes with cytoplasmic poly(A)-binding proteins (PABP) and PAN2-PAN3 and CCR4-NOT deadenylases. The AGO-GW182 complexes mediate silencing of the target mRNA through induction of translational repression and/or mRNA degradation. Although the GW182 proteins are a subject of extensive experimental research in the recent years, very little is known about their origin and evolution. Here, based on complex functional annotation and phylogenetic analyses, we reveal 448 members of the GW182 protein family from the earliest animals to humans. Our results indicate that a single-copy GW182/TNRC6C progenitor gene arose with the emergence of multicellularity and it multiplied in the last common ancestor of vertebrates in 2 rounds of whole genome duplication (WGD) resulting in 3 genes. Before the divergence of vertebrates, both the AGO- and CCR4-NOT-binding regions of GW182s showed significant acceleration in the accumulation of amino acid changes, suggesting functional adaptation toward higher specificity to the molecules of the silencing complex. We conclude that the silencing ability of the GW182 proteins improves with higher position in the taxonomic classification and increasing complexity of the organism. The first reconstruction of the molecular journey of GW182 proteins from the ancestral metazoan protein to the current mammalian configuration provides new insight into development of the miRNA-dependent post-transcriptional silencing pathway in animals.

A preliminary proteomic characterisation of extracellular vesicles released by the ovine parasitic nematode, Teladorsagia circumcincta

Teladorsagia circumcincta is a major cause of ovine parasitic gastroenteritis in temperate climatic regions. The development of high levels of anthelmintic resistance in this nematode species challenges its future control. Recent research indicates that many parasite species release extracellular vesicles into their environment, many of which have been classified as endocytic in origin, termed exosomes. These vesicles are considered to play important roles in the intercellular communication between parasites and their hosts, and thus represent potentially useful targets for novel control strategies. Here, we demonstrate that exosome-like extracellular vesicles can be isolated from excretory-secretory (ES) products released by T. circumcincta fourth stage larvae (Tci-L4ES). Furthermore, we perform a comparative proteomic analysis of vesicle-enriched and vesicle-free Tci-L4ES. Approximately 73% of the proteins identified in the vesicle-enriched fraction were unique to this fraction, whilst the remaining 27% were present in both vesicle-enriched and vesicle-free fraction. These unique proteins included structural proteins, nuclear proteins, metabolic proteins, proteolytic enzymes and activation-associated secreted proteins. Finally, we demonstrate that molecules present within the vesicles-enriched material are targets of the IgA and IgG response in T. circumcincta infected sheep, and could potentially represent useful targets for future vaccine intervention studies.

The genomic basis of parasitism in the Strongyloides clade of nematodes

Soil-transmitted nematodes, including the Strongyloides genus, cause one of the most prevalent neglected tropical diseases. Here we compare the genomes of four Strongyloides species, including the human pathogen Strongyloides stercoralis, and their close relatives that are facultatively parasitic (Parastrongyloides trichosuri) and free-living (Rhabditophanes sp. KR3021). A significant paralogous expansion of key gene families--families encoding astacin-like and SCP/TAPS proteins--is associated with the evolution of parasitism in this clade. Exploiting the unique Strongyloides life cycle, we compare the transcriptomes of the parasitic and free-living stages and find that these same gene families are upregulated in the parasitic stages, underscoring their role in nematode parasitism.

Recent Advances in Elucidating Nematode Moulting - Prospects of Using Oesophagostomum dentatum as a Model

There are major gaps in our knowledge of many molecular biological processes that take place during the development of parasitic nematodes, in spite of the fact that understanding such processes could lead to new ways of treating and controlling parasitic diseases via the disruption of one or more biological pathways in the parasites. Progress in genomics, transcriptomics, proteomics and bioinformatics now provides unique opportunities to investigate the molecular basis of key developmental processes in parasitic nematodes. The porcine nodule worm, Oesophagostomum dentatum, represents a large order (Strongylida) of socioeconomically important nematodes, and provides a useful platform for exploring molecular developmental processes, particularly given that this nematode can be grown and maintained in culture in vitro for periods longer than most other nematodes of this order. In this article, we focus on the moulting process (ecdysis) in nematodes; review recent advances in our understanding of molecular aspects of moulting in O. dentatum achieved by using integrated proteomic-bioinformatic tools and discuss key implications and future prospects for research in this area, also with respect to developing new anti-nematode interventions and biotechnological outcomes.

The Validation of Nematode-Specific Acetylcholine-Gated Chloride Channels as Potential Anthelmintic Drug Targets

New compounds are needed to treat parasitic nematode infections in humans, livestock and plants. Small molecule anthelmintics are the primary means of nematode parasite control in animals; however, widespread resistance to the currently available drug classes means control will be impossible without the introduction of new compounds. Adverse environmental effects associated with nematocides used to control plant parasitic species are also motivating the search for safer, more effective compounds. Discovery of new anthelmintic drugs in particular has been a serious challenge due to the difficulty of obtaining and culturing target parasites for high-throughput screens and the lack of functional genomic techniques to validate potential drug targets in these pathogens. We present here a novel strategy for target validation that employs the free-living nematode Caenorhabditis elegans to demonstrate the value of new ligand-gated ion channels as targets for anthelmintic discovery. Many successful anthelmintics, including ivermectin, levamisole and monepantel, are agonists of pentameric ligand-gated ion channels, suggesting that the unexploited pentameric ion channels encoded in parasite genomes may be suitable drug targets. We validated five members of the nematode-specific family of acetylcholine-gated chloride channels as targets of agonists with anthelmintic properties by ectopically expressing an ivermectin-gated chloride channel, AVR-15, in tissues that endogenously express the acetylcholine-gated chloride channels and using the effects of ivermectin to predict the effects of an acetylcholine-gated chloride channel agonist. In principle, our strategy can be applied to validate any ion channel as a putative anti-parasitic drug target.

Molecular cloning and characterization of arginine kinase gene of Toxocara canis

Toxocara canis is an important gastrointestinal nematode of dogs and also a causative agent of visceral larva migrans in humans. Arginine kinase (AK) gene is one of the important biomolecule of phosphagen kinase of T. canis which is emerging as an exciting novel diagnostic target in toxocarosis. The present study was carried out to clone and characterize AK gene of T. canis for future utilization as a diagnostic molecule. Total RNA was extracted from intact adult worms and reverse transcription was done with oligo dT primers to obtain complementary DNA (cDNA). Polymerase chain reaction (PCR) was carried out using cDNA as template with specific primers which amplified a product of 1,202 bp. The amplicon was cloned into pDrive cloning vector and clone was confirmed by colony PCR and restriction endonuclease analysis. Sequence analysis of the gene showed 99.8 and 77.9 % homology with the published AK gene of T. canis (EF015466.1) and Ascaris suum respectively. Structural analysis shown that the mature AK protein consist of 400 amino acids with a molecular wt of 45360.73 Da. Further expression studies are required for producing the recombinant protein for its evaluation in the diagnosis of T. canis infection in humans as well as in adult dogs.

Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes

The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.

Cracking the nodule worm code advances knowledge of parasite biology and biotechnology to tackle major diseases of livestock

Many infectious diseases caused by eukaryotic pathogens have a devastating, long-term impact on animal health and welfare. Hundreds of millions of animals are affected by parasitic nematodes of the order Strongylida. Unlocking the molecular biology of representatives of this order, and understanding nematode-host interactions, drug resistance and disease using advanced technologies could lead to entirely new ways of controlling the diseases that they cause. Oesophagostomum dentatum (nodule worm; superfamily Strongyloidea) is an economically important strongylid nematode parasite of swine worldwide. The present article reports recent advances made in biology and animal biotechnology through the draft genome and developmental transcriptome of O. dentatum, in order to support biological research of this and related parasitic nematodes as well as the search for new and improved interventions. This first genome of any member of the Strongyloidea is 443 Mb in size and predicted to encode 25,291 protein-coding genes. Here, we review the dynamics of transcription throughout the life cycle of O. dentatum, describe double-stranded RNA interference (RNAi) machinery and infer molecules involved in development and reproduction, and in inducing or modulating immune responses or disease. The secretome predicted for O. dentatum is particularly rich in peptidases linked to interactions with host tissues and/or feeding activity, and a diverse array of molecules likely involved in immune responses. This research progress provides an important resource for future comparative genomic and molecular biological investigations as well as for biotechnological research toward new anthelmintics, vaccines and diagnostic tests.

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.