Comparative genomics of gene expression in the parasitic and free-living nematodes Strongyloides stercoralis and Caenorhabditis elegans

Model nematodes as a practical innovation to promote high throughput screening of natural products for anthelmintics discovery in South Asia: Current challenges, proposed practical and conceptual solutions

With the increasing prevalence of anthelmintic resistance in animals recorded globally, and the threat of resistance in human helminths, the need for novel anthelmintic drugs is greater than ever. Most research aimed at discovering novel anthelmintic leads relies on high throughput screening (HTS) of large libraries of synthetic small molecules in industrial and academic settings in developed countries, even though it is the tropical countries that are most plagued by helminth infections. Tropical countries, however, have the advantage of possessing a rich flora that may yield natural products (NP) with promising anthelmintic activity. Focusing on South Asia, which produces one of the world's highest research outputs in NP and NP-based anthelmintic discovery, we find that limited basic research and funding, a lack of awareness of the utility of model organisms, poor industry-academia partnerships and lack of technological innovations greatly limit anthelmintics research in the region. Here we propose that utilizing model organisms including the free-living nematode Caenorhabditis elegans, that can potentially allow rapid target identification of novel anthelmintics, and Oscheius tipulae, a closely related, free-living nematode which is found abundantly in soil in hotter temperatures, could be a much-needed innovation that can enable cost-effective and efficient HTS of NPs for discovering compounds with anthelmintic/antiparasitic potential in South Asia and other tropical regions that historically have devoted limited funding for such research. Additionally, increased collaborations at the national, regional and international level between parasitologists and pharmacologists/ethnobotanists, setting up government-industry-academia partnerships to fund academic research, creating a centralized, regional collection of plant extracts or purified NPs as a dereplication strategy and HTS library, and holding regional C. elegans/O. tipulae-based anthelmintics workshops and conferences to share knowledge and resources regarding model organisms may collectively promote and foster a NP-based anthelmintics landscape in South Asia and beyond.

Anthelmintic potency of Rumex crispus L. extracts against Caenorhabditis elegans and non-targeted identification of the bioactive compounds

Traditional healers and ethnoveterinary therapists use several medicinal plants, such as Rumex crispus L., to treat endoparasite infections. R. crispus has been established by researchers to be effective agasint a few parasitic worms. In this study, we evaluated the potency of R. crispus extracts on the model organism, Caenorhabditis elegans and the bioactive compounds of the extracts were also identified. The solvent extracts of R. crispus were tested against C. elegans for up to 72 h. The effect of the extracts on C. elegans was examined using light microscopy (LM) and scanning electron microscopy (SEM). LM and SEM analysis showed damage on the body wall, reduced body and slight modifications of the nematode organs. The lethality test reveals a significant reduction in the viability of the nematode with the water extract of leaf (LF-WAE), among others, having the strongest potency against the nematode, with 83% lethality. Anlysis done with Fourier-transform infrared spectroscopy (FTIR) spectra reveals various characteristic vibration bands and fingerprint bands at 3400–600 cm⁻¹, identifying phenols, organic acids, aromatics, amines, among others in the plant. The compounds were identified with liquid chromatography-mass spectrometry (LC-MS), under the categories of flavonoids, steroidal alkaloids and proanthocyanidin. In conclusion, this study confirmed that R. crispus has anthelmintic potential, using standardised C. elegans models as a tool and suggests that there could be novel compounds yet to be explored in the studied plant that could be of great benefit to livestock and humans.

Combating Parasitic Nematode Infections, Newly Discovered Antinematode Compounds from Marine Epiphytic Bacteria

Parasitic nematode infections cause debilitating diseases and impede economic productivity. Antinematode chemotherapies are fundamental to modern medicine and are also important for industries including agriculture, aquaculture and animal health. However, the lack of suitable treatments for some diseases and the rise of nematode resistance to many available therapies necessitates the discovery and development of new drugs. Here, marine epiphytic bacteria represent a promising repository of newly discovered antinematode compounds. Epiphytic bacteria are ubiquitous on marine surfaces where they are under constant pressure of grazing by bacterivorous predators (e.g., protozoans and nematodes). Studies have shown that these bacteria have developed defense strategies to prevent grazers by producing toxic bioactive compounds. Although several active metabolites against nematodes have been identified from marine bacteria, drug discovery from marine microorganisms remains underexplored. In this review, we aim to provide further insight into the need and potential for marine epiphytic bacteria to become a new source of antinematode drugs. We discuss current and emerging strategies, including culture-independent high throughput screening and the utilization of Caenorhabditis elegans as a model target organism, which will be required to advance antinematode drug discovery and development from marine microbial sources.

Comparative characterization of two galectins excreted-secreted from intestine-dwelling parasitic versus free-living females of the soil-transmitted nematode Strongyloides

Helminths are complex pathogens that ensure their long-term survival by influencing the immune responses of their host. Excretory/secretory products (ESP) can exert immunoregulatory effects which foster parasite survival. Galectins represent a widespread group of β-galactoside-binding proteins which are involved in a multitude of biological processes operative in parasite-host interaction. We had earlier identified seven galectins in Strongyloides ratti, four of them detected in the ESP of distinct developmental stages of the parasite. In the present report, we focused on the characterization of two of them, Sr-galectin-1 (Sr-Gal-1) and Sr-galectin-3 (Sr-Gal-3). While Sr-Gal-3 expression was strongest in parasitic females, Sr-Gal-1 was predominantly expressed in free-living females. Both proteins were cloned and recombinantly expressed in an E. coli expression system. Their glycan-binding activity was verified by haemagglutination and glycan array analysis. Furthermore, primary immunological activities of the Sr-galectins were initially investigated by the application of an in vitro mucosal 3D-culture model, comprising of mucosa-associated epithelial and dendritic cells. The Sr-galectins stimulated preferentially the release of the type 2 cytokines thymic stromal lymphopoietin and IL-22, a first indication for immunoregulatory activity. In addition, the Sr-galectins dose-dependently fostered cell migration. Our results confirm the importance of these carbohydrate-binding proteins in host-parasite-interaction by indicating possible interaction with the host mucosa-associated cells.

The genome ofStrongyloidesspp. gives insights into protein families with a putative role in nematode parasitism

Parasitic nematodes are important and abundant parasites adapted to live a parasitic lifestyle, with these adaptations all aimed at facilitating their survival and reproduction in their hosts. The recently sequenced genomes of fourStrongyloidesspecies, gastrointestinal parasites of humans and other animals, alongside transcriptomic and proteomic analysis of free-living and parasitic stages of their life cycles have revealed a number of protein families with a putative role in their parasitism. Many of these protein families have also been associated with parasitism in other parasitic nematode species, suggesting that these proteins may play a fundamental role in nematode parasitism more generally. Here, we review key protein families that have a putative role inStrongyloides’ parasitism – acetylcholinesterases, astacins, aspartic proteases, prolyl oligopeptidases, proteinase inhibitors (trypsin inhibitors and cystatins), SCP/TAPS and transthyretin-like proteins – and the evidence for their key, yet diverse, roles in the parasitic lifestyle.

Transcriptional and morphological changes in the transition from mycetophagous to phytophagous phase in the plant-parasitic nematode Bursaphelenchus xylophilus

Drastic physiological and morphological changes in parasites are crucial for the establishment of a successful infection. The nematode Bursaphelenchus xylophilus is the pathogenic agent of pine wilt disease, and little is known about the physiology and morphology in this nematode at the initial stage of infection. In this study, we devised an infection system using pine stem cuttings that allowed us to observe transcriptional and morphological changes in the host-infecting phytophagous phase. We found that 60 genes enriched in xenobiotic detoxification were up-regulated in two independent post-inoculation events, whereas down-regulation was observed in multiple members of collagen gene families. After 48 h of inoculation, the tails in some of the adult females exposed to the host changed in morphology. These results suggest that B. xylophilus may change its physiology and morphology to protect itself and to adapt to the host pine wood environment.

RNA-Seq Based Identification of Candidate Parasitism Genes of Cereal Cyst Nematode (Heterodera avenae) during Incompatible Infection to Aegilops variabilis

One of the reasons for the progressive yield decline observed in cereals production is the rapid build-up of populations of the cereal cyst nematode (CCN, Heterodera avenae). These nematodes secrete so-call effectors into their host plant to suppress the plant defense responses, alter plant signaling pathways and then induce the formation of syncytium after infection. However, little is known about its molecular mechanism and parasitism during incompatible infection. To gain insight into its repertoire of parasitism genes, we investigated the transcriptome of the early parasitic second-stage (30 hours, 3 days and 9 days post infection) juveniles of the CCN as well as the CCN infected tissue of the host Aegilops variabilis by Illumina sequencing. Among all assembled unigenes, 681 putative genes of parasitic nematode were found, in which 56 putative effectors were identified, including novel pioneer genes and genes corresponding to previously reported effectors. All the 681 CCN unigenes were mapped to 229 GO terms and 200 KEGG pathways, including growth, development and several stimulus-related signaling pathways. Sixteen clusters were involved in the CCN unigene expression atlas at the early stages during infection process, and three of which were significantly gene-enriched. Besides, the protein-protein interaction network analysis revealed 35 node unigenes which may play an important role in the plant-CCN interaction. Moreover, in a comparison of differentially expressed genes between the pre-parasitic juveniles and the early parasitic juveniles, we found that hydrolase activity was up-regulated in pre J2s whereas binding activity was upregulated in infective J2s. RT-qPCR analysis on some selected genes showed detectable expression, indicating possible secretion of the proteins and putative role in infection. This study provided better insights into the incompatible interaction between H. avenae and the host plant Ae. varabilis. Moreover, RNAi targets with potential lethality were screened out and primarily validated, which provide candidates for engineering-based control of cereal cyst nematode in crops breeding.

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.

Strongyloidiasis with emphasis on human infections and its different clinical forms

Strongyloidiasis (caused by Strongyloides stercoralis, and to a lesser extent by Strongyloides fuelleborni) is one of the most neglected tropical diseases with endemic areas and affecting more than 100 million people worldwide. Chronic infections in endemic areas can be maintained for decades through the autoinfective cycle with the L3 filariform larvae. In these endemic areas, misdiagnosis, inadequate treatment and the facilitation of the hyperinfection syndrome by immunosuppression are frequent and contribute to a high mortality rate. Despite the serious health impact of strongyloidiasis, it is a neglected disease and very little is known about this parasite and the disease when compared to other helminth infections. Control of the disease is difficult because of the many gaps in our knowledge of strongyloidiasis. We examine the recent literature on different aspects of strongyloidiasis with emphasis in those aspects that need further research.

Systems biology studies of adult paragonimus lung flukes facilitate the identification of immunodominant parasite antigens

Background

Paragonimiasis is a food-borne trematode infection acquired by eating raw or undercooked crustaceans. It is a major public health problem in the far East, but it also occurs in South Asia, Africa, and in the Americas. Paragonimus worms cause chronic lung disease with cough, fever and hemoptysis that can be confused with tuberculosis or other non-parasitic diseases. Treatment is straightforward, but diagnosis is often delayed due to a lack of reliable parasitological or serodiagnostic tests. Hence, the purpose of this study was to use a systems biology approach to identify key parasite proteins that may be useful for development of improved diagnostic tests.

Methodology/Principal Findings

The transcriptome of adult Paragonimus kellicotti was sequenced with Illumina technology. Raw reads were pre-processed and assembled into 78,674 unique transcripts derived from 54,622 genetic loci, and 77,123 unique protein translations were predicted. A total of 2,555 predicted proteins (from 1,863 genetic loci) were verified by mass spectrometric analysis of total worm homogenate, including 63 proteins lacking homology to previously characterized sequences. Parasite proteins encoded by 321 transcripts (227 genetic loci) were reactive with antibodies from infected patients, as demonstrated by immunoaffinity purification and high-resolution liquid chromatography-mass spectrometry. Serodiagnostic candidates were prioritized based on several criteria, especially low conservation with proteins in other trematodes. Cysteine proteases, MFP6 proteins and myoglobins were abundant among the immunoreactive proteins, and these warrant further study as diagnostic candidates.

Conclusions

The transcriptome, proteome and immunolome of adult P. kellicotti represent a major advance in the study of Paragonimus species. These data provide a powerful foundation for translational research to develop improved diagnostic tests. Similar integrated approaches may be useful for identifying novel targets for drugs and vaccines in the future.

Paragonimiasis is a food-borne trematode infection that people acquire when they eat raw or undercooked crustaceans. Disease symptoms (including cough, fever, blood in sputum, etc.) can be similar to those observed in patients with tuberculosis or bacterial pneumonia, frequently resulting in misdiagnosis. Although the infection is relatively easy to treat, diagnosis is complicated. Available diagnostic assays rely on total parasite homogenate to facilitate the detection of Paragonimus-specific antibodies in patients. Though these blot-based assays have shown high sensitivity and specificity, they are inconvenient because total parasite homogenate is not readily available. This study used next generation genomic and proteomic methods to identify transcripts and proteins expressed in adult Paragonimus flukes. We then used sera from patients infected with P. kellicotti to isolate immunoreactive proteins, and these were analyzed by mass spectrometry. The annotated transcriptome and the associated proteome of the antibody immune response represent a significant advance in research on Paragonimus. This information will be a valuable resource for further research on Paragonimus and paragonimiasis. Thus this project illustrates the potential power of employing systems biology for translational research in parasitology.

Genome-wide tissue-specific gene expression, co-expression and regulation of co-expressed genes in adult nematode Ascaris suum

BACKGROUND

Caenorhabditis elegans has traditionally been used as a model for studying nematode biology, but its small size limits the ability for researchers to perform some experiments such as high-throughput tissue-specific gene expression studies. However, the dissection of individual tissues is possible in the parasitic nematode Ascaris suum due to its relatively large size. Here, we take advantage of the recent genome sequencing of Ascaris suum and the ability to physically dissect its separate tissues to produce a wide-scale tissue-specific nematode RNA-seq datasets, including data on three non-reproductive tissues (head, pharynx, and intestine) in both male and female worms, as well as four reproductive tissues (testis, seminal vesicle, ovary, and uterus). We obtained fundamental information about the biology of diverse cell types and potential interactions among tissues within this multicellular organism.

METHODOLOGY/PRINCIPAL FINDINGS

Overexpression and functional enrichment analyses identified many putative biological functions enriched in each tissue studied, including functions which have not been previously studied in detail in nematodes. Putative tissue-specific transcriptional factors and corresponding binding motifs that regulate expression in each tissue were identified, including the intestine-enriched ELT-2 motif/transcription factor previously described in nematode intestines. Constitutively expressed and novel genes were also characterized, with the largest number of novel genes found to be overexpressed in the testis. Finally, a putative acetylcholine-mediated transcriptional network connecting biological activity in the head to the male reproductive system is described using co-expression networks, along with a similar ecdysone-mediated system in the female.

CONCLUSIONS/SIGNIFICANCE

The expression profiles, co-expression networks and co-expression regulation of the 10 tissues studied and the tissue-specific analysis presented here are a valuable resource for studying tissue-specific biological functions in nematodes.

The dauer hypothesis and the evolution of parasitism: 20 years on and still going strong

How any complex trait has evolved is a fascinating question, yet the evolution of parasitism among the nematodes is arguably one of the most arresting. How did free-living nematodes cross that seemingly insurmountable evolutionary chasm between soil dwelling and survival inside another organism? Which of the many finely honed responses to the varied and harsh environments of free-living nematodes provided the material upon which natural selection could act? Although several complementary theories explain this phenomenon, I will focus on the dauer hypothesis. The dauer hypothesis posits that the arrested third-stage dauer larvae of free-living nematodes such as Caenorhabditis elegans are, due to their many physiological similarities with infective third-stage larvae of parasitic nematodes, a pre-adaptation to parasitism. If so, then a logical extension of this hypothesis is that the molecular pathways which control entry into and recovery from dauer formation by free-living nematodes in response to environmental cues have been co-opted to control the processes of infective larval arrest and activation in parasitic nematodes. The molecular machinery that controls dauer entry and exit is present in a wide range of parasitic nematodes. However, the developmental outputs of the different pathways are both conserved and divergent, not only between populations of C. elegans or between C. elegans and parasitic nematodes but also between different species of parasitic nematodes. Thus the picture that emerges is more nuanced than originally predicted and may provide insights into the evolution of such an interesting and complex trait.

Transcriptome analyses reveal protein and domain families that delineate stage-related development in the economically important parasitic nematodes, Ostertagia ostertagi and Cooperia oncophora

Background

Cooperia oncophora and Ostertagia ostertagi are among the most important gastrointestinal nematodes of cattle worldwide. The economic losses caused by these parasites are on the order of hundreds of millions of dollars per year. Conventional treatment of these parasites is through anthelmintic drugs; however, as resistance to anthelmintics increases, overall effectiveness has begun decreasing. New methods of control and alternative drug targets are necessary. In-depth analysis of transcriptomic data can help provide these targets.

Results

The assembly of 8.7 million and 11 million sequences from C. oncophora and O. ostertagi, respectively, resulted in 29,900 and 34,792 transcripts. Among these, 69% and 73% of the predicted peptides encoded by C. oncophora and O. ostertagi had homologues in other nematodes. Approximately 21% and 24% were constitutively expressed in both species, respectively; however, the numbers of transcripts that were stage specific were much smaller (~1% of the transcripts expressed in a stage). Approximately 21% of the transcripts in C. oncophora and 22% in O. ostertagi were up-regulated in a particular stage. Functional molecular signatures were detected for 46% and 35% of the transcripts in C. oncophora and O. ostertagi, respectively. More in-depth examinations of the most prevalent domains led to knowledge of gene expression changes between the free-living (egg, L1, L2 and L3 sheathed) and parasitic (L3 exsheathed, L4, and adult) stages. Domains previously implicated in growth and development such as chromo domains and the MADF domain tended to dominate in the free-living stages. In contrast, domains potentially involved in feeding such as the zinc finger and CAP domains dominated in the parasitic stages. Pathway analyses showed significant associations between life-cycle stages and peptides involved in energy metabolism in O. ostertagi whereas metabolism of cofactors and vitamins were specifically up-regulated in the parasitic stages of C. oncophora. Substantial differences were observed also between Gene Ontology terms associated with free-living and parasitic stages.

Conclusions

This study characterized transcriptomes from multiple life stages from both C. oncophora and O. ostertagi. These data represent an important resource for studying these parasites. The results of this study show distinct differences in the genes involved in the free-living and parasitic life cycle stages. The data produced will enable better annotation of the upcoming genome sequences and will allow future comparative analyses of the biology, evolution and adaptation to parasitism in nematodes.

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

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

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

Transcriptomic analysis of four developmental stages of Strongyloides venezuelensis

Strongyloides venezuelensis is one of some 50 species of genus Strongyloides, obligate gastrointestinal parasites of vertebrates, responsible for strongyloidiasis in humans and other domestic/companion animals. Although S. venezuelensis has been widely used as a model species for studying human/animal strongyloidiasis, the sequence information for this species has been quite limited. To create a more comprehensive catalogue of expressed genes for identification of genes potentially involved in animal parasitism, we conducted a de novo sequencing analysis of the transcriptomes from four developmental stages of S. venezuelensis, using a Roche 454 GS FLX Titanium pyrosequencing platform. A total of 14,573 contigs were produced after de novo assemblies of over 2 million sequencing reads and formed a dataset "Vene454". BLAST homology search of Vene454 against proteome and transcriptome data from other animal-parasitic and non-animal-parasitic nematode species revealed several interesting genes, which may be potentially related to animal parasitism, including nicotinamide phosphoribosyltransferase and ferrochelatase. The Vene454 dataset analysis also enabled us to identify transcripts that are specifically enriched in each developmental stage. This work represents the first large-scale transcriptome analysis of S. venezuelensis and the first study to examine the transcriptome of the lung L3 developmental stage of any Strongyloides species. The results not only will serve as valuable resources for future functional genomics analyses to understand the molecular aspects of animal parasitism, but also will provide essential information for ongoing whole genome sequencing efforts in this species.

Comparison of the transcriptional profiles of head and body lice

Head and body lice are both blood-feeding parasites of humans although only the body louse is a potent disease vector. In spite of numerous morphological and life history differences, head and body lice have recently been hypothesized to be ecotypes of the same species. We took a comparative genomics approach to measure nucleotide diversity by comparing expressed sequence tag data sets from head and body lice. A total of 10 771 body louse and 10 770 head louse transcripts were predicted from a combined assembly of Roche 454 and Illumina sequenced cDNAs from whole body tissues collected at all life stages and during pesticide exposure and bacterial infection treatments. Illumina reads mapped to the 10 775 draft body louse gene models from the whole genome assembly predicted nine presence/absence differences, but PCR confirmation resulted in a single gene difference. Read per million base pair estimates indicated that 14 genes showed significant differential expression between head and body lice under our treatment conditions. One novel microRNA was predicted in both lice species and 99% of the 544 transcripts from Candidatus riesia indicate that they share the same endosymbiont. Overall, few differences exist, which supports the hypothesis that these two organisms are ecotypes of the same species.

Functional reconstitution of Haemonchus contortus acetylcholine receptors in Xenopus oocytes provides mechanistic insights into levamisole resistance

BACKGROUND AND PURPOSE

The cholinergic agonist levamisole is widely used to treat parasitic nematode infestations. This anthelmintic drug paralyses worms by activating a class of levamisole-sensitive acetylcholine receptors (L-AChRs) expressed in nematode muscle cells. However, levamisole efficacy has been compromised by the emergence of drug-resistant parasites, especially in gastrointestinal nematodes such as Haemonchus contortus. We report here the first functional reconstitution and pharmacological characterization of H. contortus L-AChRs in a heterologous expression system.

EXPERIMENTAL APPROACH

In the free-living nematode Caenorhabditis elegans, five AChR subunit and three ancillary protein genes are necessary in vivo and in vitro to synthesize L-AChRs. We have cloned the H. contortus orthologues of these genes and expressed them in Xenopus oocytes. We reconstituted two types of H. contortus L-AChRs with distinct pharmacologies by combining different receptor subunits.

KEY RESULTS

The Hco-ACR-8 subunit plays a pivotal role in selective sensitivity to levamisole. As observed with C. elegans L-AChRs, expression of H. contortus receptors requires the ancillary proteins Hco-RIC-3, Hco-UNC-50 and Hco-UNC-74. Using this experimental system, we demonstrated that a truncated Hco-UNC-63 L-AChR subunit, which was specifically detected in a levamisole-resistant H. contortus isolate, but not in levamisole-sensitive strains, hampers the normal function of L-AChRs, when co-expressed with its full-length counterpart.

CONCLUSIONS AND IMPLICATIONS

We provide the first functional evidence for a putative molecular mechanism involved in levamisole resistance in any parasitic nematode. This expression system will provide a means to analyse molecular polymorphisms associated with drug resistance at the electrophysiological level.

Developmental plasticity and the evolution of parasitism in an unusual nematode, Parastrongyloides trichosuri

BACKGROUND

Parasitism is an important life history strategy in many metazoan taxa. This is particularly true of the Phylum Nematoda, in which parasitism has evolved independently at least nine times. The apparent ease with which parasitism has evolved amongst nematodes may, in part, be due to a feature of nematode development acting as a pre-adaptation for the transition from a free-living to a parasitic life history. One candidate pre-adaptive feature for evolution in terrestrial nematodes is the dauer larva, a developmentally arrested morph formed in response to environmental signals.

RESULTS

We investigated the role of dauer development in the nematode, Parastrongyloides trichosuri, which has retained a complete free-living life cycle in addition to a life cycle as a mammalian gastrointestinal parasite. We show that the developmental switch between these life histories is sensitive to the same environmental cues as dauer arrest in free-living nematodes, including sensitivity to a chemical cue produced by the free-living stages. Furthermore, we show that genetic variation for the sensitivity of the cue(s) exists in natural populations of P. trichosuri, such that we derived inbred lines that were largely insensitive to the cue and other lines that were supersensitive to the cue.

CONCLUSIONS

For this parasitic clade, and perhaps more widely in the phylum, the evolution of parasitism co-opted the dauer switch of a free-living ancestor. This lends direct support to the hypothesis that the switch to developmental arrest in the dauer larva acted as a pre-adaptation for the evolution of parasitism, and suggests that the sensory transduction machinery downstream of the cue may have been similarly co-opted and modified.

In silico secretome analysis approach for next generation sequencing transcriptomic data

BACKGROUND

Excretory/secretory proteins (ESPs) play a major role in parasitic infection as they are present at the host-parasite interface and regulate host immune system. In case of parasitic helminths, transcriptomics has been used extensively to understand the molecular basis of parasitism and for developing novel therapeutic strategies against parasitic infections. However, none of transcriptomic studies have extensively covered ES protein prediction for identifying novel therapeutic targets, especially as parasites adopt non-classical secretion pathways.

RESULTS

We developed a semi-automated computational approach for prediction and annotation of ES proteins using transcriptomic data from next generation sequencing platforms. For the prediction of non-classically secreted proteins, we have used an improved computational strategy, together with homology matching to a dataset of experimentally determined parasitic helminth ES proteins. We applied this protocol to analyse 454 short reads of parasitic nematode, Strongyloides ratti. From 296231 reads, we derived 28901 contigs, which were translated into 20877 proteins. Based on our improved ES protein prediction pipeline, we identified 2572 ES proteins, of which 407 (1.9%) proteins have classical N-terminal signal peptides, 923 (4.4%) were computationally identified as non-classically secreted while 1516 (7.26%) were identified by homology to experimentally identified parasitic helminth ES proteins. Out of 2572 ES proteins, 2310 (89.8%) ES proteins had homologues in the free-living nematode Caenorhabditis elegans and 2220 (86.3%) in parasitic nematodes. We could functionally annotate 1591 (61.8%) ES proteins with protein families and domains and establish pathway associations for 691 (26.8%) proteins. In addition, we have identified 19 representative ES proteins, which have no homologues in the host organism but homologous to lethal RNAi phenotypes in C. elegans, as potential therapeutic targets.

CONCLUSION

We report a comprehensive approach using freely available computational tools for the secretome analysis of NGS data. This approach has been applied to S. ratti 454 transcriptomic data for in silico excretory/secretory proteins prediction and analysis, providing a foundation for developing new therapeutic solutions for parasitic infections.

Life cycle stage-resolved proteomic analysis of the excretome/secretome from Strongyloides ratti--identification of stage-specific proteases

A wide range of biomolecules, including proteins, are excreted and secreted from helminths and contribute to the parasite's successful establishment, survival, and reproduction in an adverse habitat. Excretory and secretory proteins (ESP) are active at the interface between parasite and host and comprise potential targets for intervention. The intestinal nematode Strongyloides spp. exhibits an exceptional developmental plasticity in its life cycle characterized by parasitic and free-living generations. We investigated ESP from infective larvae, parasitic females, and free-living stages of the rat parasite Strongyloides ratti, which is genetically very similar to the human pathogen, Strongyloides stercoralis. Proteomic analysis of ESP revealed 586 proteins, with the largest number of stage-specific ESP found in infective larvae (196), followed by parasitic females (79) and free-living stages (35). One hundred and forty proteins were identified in all studied stages, including anti-oxidative enzymes, heat shock proteins, and carbohydrate-binding proteins. The stage-selective ESP of (1) infective larvae included an astacin metalloproteinase, the L3 Nie antigen, and a fatty acid retinoid-binding protein; (2) parasitic females included a prolyl oligopeptidase (prolyl serine carboxypeptidase), small heat shock proteins, and a secreted acidic protein; (3) free-living stages included a lysozyme family member, a carbohydrate-hydrolyzing enzyme, and saponin-like protein. We verified the differential expression of selected genes encoding ESP by qRT-PCR. ELISA analysis revealed the recognition of ESP by antibodies of S. ratti-infected rats. A prolyl oligopeptidase was identified as abundant parasitic female-specific ESP, and the effect of pyrrolidine-based prolyl oligopeptidase inhibitors showed concentration- and time-dependent inhibitory effects on female motility. The characterization of stage-related ESP from Strongyloides will help to further understand the interaction of this unique intestinal nematode with its host.

HelmCoP: an online resource for helminth functional genomics and drug and vaccine targets prioritization

A vast majority of the burden from neglected tropical diseases result from helminth infections (nematodes and platyhelminthes). Parasitic helminthes infect over 2 billion, exerting a high collective burden that rivals high-mortality conditions such as AIDS or malaria, and cause devastation to crops and livestock. The challenges to improve control of parasitic helminth infections are multi-fold and no single category of approaches will meet them all. New information such as helminth genomics, functional genomics and proteomics coupled with innovative bioinformatic approaches provide fundamental molecular information about these parasites, accelerating both basic research as well as development of effective diagnostics, vaccines and new drugs. To facilitate such studies we have developed an online resource, HelmCoP (Helminth Control and Prevention), built by integrating functional, structural and comparative genomic data from plant, animal and human helminthes, to enable researchers to develop strategies for drug, vaccine and pesticide prioritization, while also providing a useful comparative genomics platform. HelmCoP encompasses genomic data from several hosts, including model organisms, along with a comprehensive suite of structural and functional annotations, to assist in comparative analyses and to study host-parasite interactions. The HelmCoP interface, with a sophisticated query engine as a backbone, allows users to search for multi-factorial combinations of properties and serves readily accessible information that will assist in the identification of various genes of interest. HelmCoP is publicly available at: http://www.nematode.net/helmcop.html.

Microarray-based analysis of differential gene expression between infective and noninfective larvae of Strongyloides stercoralis

BACKGROUND

Differences between noninfective first-stage (L1) and infective third-stage (L3i) larvae of parasitic nematode Strongyloides stercoralis at the molecular level are relatively uncharacterized. DNA microarrays were developed and utilized for this purpose.

METHODS AND FINDINGS

Oligonucleotide hybridization probes for the array were designed to bind 3,571 putative mRNA transcripts predicted by analysis of 11,335 expressed sequence tags (ESTs) obtained as part of the Nematode EST project. RNA obtained from S. stercoralis L3i and L1 was co-hybridized to each array after labeling the individual samples with different fluorescent tags. Bioinformatic predictions of gene function were developed using a novel cDNA Annotation System software. We identified 935 differentially expressed genes (469 L3i-biased; 466 L1-biased) having two-fold expression differences or greater and microarray signals with a p value<0.01. Based on a functional analysis, L1 larvae have a larger number of genes putatively involved in transcription (p = 0.004), and L3i larvae have biased expression of putative heat shock proteins (such as hsp-90). Genes with products known to be immunoreactive in S. stercoralis-infected humans (such as SsIR and NIE) had L3i biased expression. Abundantly expressed L3i contigs of interest included S. stercoralis orthologs of cytochrome oxidase ucr 2.1 and hsp-90, which may be potential chemotherapeutic targets. The S. stercoralis ortholog of fatty acid and retinol binding protein-1, successfully used in a vaccine against Ancylostoma ceylanicum, was identified among the 25 most highly expressed L3i genes. The sperm-containing glycoprotein domain, utilized in a vaccine against the nematode Cooperia punctata, was exclusively found in L3i biased genes and may be a valuable S. stercoralis target of interest.

CONCLUSIONS

A new DNA microarray tool for the examination of S. stercoralis biology has been developed and provides new and valuable insights regarding differences between infective and noninfective S. stercoralis larvae. Potential therapeutic and vaccine targets were identified for further study.

NEMBASE4: the nematode transcriptome resource

Nematode parasites are of major importance in human health and agriculture, and free-living species deliver essential ecosystem services. The genomics revolution has resulted in the production of many datasets of expressed sequence tags (ESTs) from a phylogenetically wide range of nematode species, but these are not easily compared. NEMBASE4 presents a single portal into extensively functionally annotated, EST-derived transcriptomes from over 60 species of nematodes, including plant and animal parasites and free-living taxa. Using the PartiGene suite of tools, we have assembled the publicly available ESTs for each species into a high-quality set of putative transcripts. These transcripts have been translated to produce a protein sequence resource and each is annotated with functional information derived from comparison with well-studied nematode species such as Caenorhabditis elegans and other non-nematode resources. By cross-comparing the sequences within NEMBASE4, we have also generated a protein family assignment for each translation. The data are presented in an openly accessible, interactive database. To demonstrate the utility of NEMBASE4, we have used the database to examine the uniqueness of the transcriptomes of major clades of parasitic nematodes, identifying lineage-restricted genes that may underpin particular parasitic phenotypes, possible viral pathogens of nematodes, and nematode-unique protein families that may be developed as drug targets.

Transcripts analysis of infective larvae of an intestinal nematode, Strongyloides venezuelensis

Free-living infective larvae of Strongyloides nematodes fulfill a number of requirements for the successful infection. They need to endure a long wait in harsh environmental conditions, like temperature, salinity, and pH, which might change drastically from time to time. Infective larvae also have to deal with pathogens and potentially hazardous free-living microbes in the environment. In addition, infective larvae must recognize the adequate host properly, and start skin penetration as quickly as possible. All these tasks are essentially important for the survival of Strongyloides nematodes, however, our knowledge is extremely limited in any one of these aspects. In order to understand how Strongyloides infective larvae meet these requirements, we examined transcripts of infective larvae by randomly sequencing cDNA clones constructed from S. venezuelensis infective larvae. After assembling successfully sequenced clones, we obtained 162 unique singletons and contigs, of which 84 had been significantly annotated. Annotated genes included those for respiratory enzymes, heat-shock proteins, neuromuscular proteins, proteases, and immunodominant antigens. Genes for lipase, small heat-shock protein, globin-like protein and cytochrome c oxidase were most abundantly transcribed, though genes of unknown functions were also abundantly transcribed. There were no hits found against NCBI or NEMABASE4 for 37 (22.3%) EST out of the total 162 EST. Although most of the transcripts were not infective larva-specific, the expression of respiration related proteins was most actively transcribed in the infective larva stage. The expression of astacin-like metalloprotease, small heat-shock protein, S. stercoralis L3Nie antigen homologue, and one unannotated and 2 novel genes was highly specific for the infective larva stage.

Transcript analysis of parasitic females of the sedentary semi-endoparasitic nematode Rotylenchulus reniformis

Rotylenchulus reniformis, the reniform nematode, is a sedentary semi-endoparasitic nematode capable of infecting >300 plant species, including a large number of crops such as cotton, soybean, and pineapple. In contrast to other economically important plant-parasitic nematodes, molecular genetic data regarding the R. reniformis transcriptome is virtually nonexistant. Herein, we present a survey of R. reniformis ESTs that were sequenced from a sedentary parasitic female cDNA library. Cluster analysis of 2004 high quality ESTs produced 123 contigs and 508 singletons for a total of 631 R. reniformis unigenes. BLASTX analyses revealed that 39% of all unigenes showed similarity to known proteins (E<or=1.0e-04). R. reniformis genes homologous to known parasitism genes were identified and included beta-1,4-endoglucanase, fatty acid- and retinol-binding proteins, and an esophageal gland cell-specific gene from Heterodera glycines. Furthermore, a putative ortholog of an enzyme involved in thiamin biosynthesis, thought to exist solely in prokaryotes, fungi, and plants, was identified. Lastly, 114 R. reniformis unigenes orthologous to RNAi-lethal Caenorhabditis elegans genes were discovered. The work described here offers a glimpse into the transcriptome of a sedentary semi-endoparasitic nematode which (i) provides the transcript sequence data necessary for investigating engineered resistance against R. reniformis and (ii) hints at the existance of a thiamin biosynthesis pathway in an animal.

Analysis of expressed sequence tags and identification of genes encoding cell-wall-degrading enzymes from the fungivorous nematode Aphelenchus avenae

Background

The fungivorus nematode, Aphelenchus avenae is widespread in soil and is found in association with decaying plant material. This nematode is also found in association with plants but its ability to cause plant disease remains largely undetermined. The taxonomic position and intermediate lifestyle of A. avenae make it an important model for studying the evolution of plant parasitism within the Nematoda. In addition, the exceptional capacity of this nematode to survive desiccation makes it an important system for study of anhydrobiosis. Expressed sequence tag (EST) analysis may therefore be useful in providing an initial insight into the poorly understood genetic background of A. avenae.

Results

We present the generation, analysis and annotation of over 5,000 ESTs from a mixed-stage A. avenae cDNA library. Clustering of 5,076 high-quality ESTs resulted in a set of 2,700 non-redundant sequences comprising 695 contigs and 2,005 singletons. Comparative analyses indicated that 1,567 (58.0%) of the cluster sequences had homologues in Caenorhabditis elegans, 1,750 (64.8%) in other nematodes, 1,321(48.9%) in organisms other than nematodes, and 862 (31.9%) had no significant match to any sequence in current protein or nucleotide databases. In addition, 1,100 (40.7%) of the sequences were functionally classified using Gene Ontology (GO) hierarchy. Similarity searches of the cluster sequences identified a set of genes with significant homology to genes encoding enzymes that degrade plant or fungal cell walls. The full length sequences of two genes encoding glycosyl hydrolase family 5 (GHF5) cellulases and two pectate lyase genes encoding polysaccharide lyase family 3 (PL3) proteins were identified and characterized.

Conclusion

We have described at least 2,214 putative genes from A. avenae and identified a set of genes encoding a range of cell-wall-degrading enzymes. This EST dataset represents a starting point for studies in a number of different fundamental and applied areas. The presence of genes encoding a battery of cell-wall-degrading enzymes in A. avenae and their similarities with genes from other plant parasitic nematodes suggest that this nematode can act not only as a fungal feeder but also a plant parasite. Further studies on genes encoding cell-wall-degrading enzymes in A. avenae will accelerate our understanding of the complex evolutionary histories of plant parasitism and the use of genes obtained by horizontal gene transfer from prokaryotes.

Transcripts analysis of the entomopathogenic nematode Steinernema carpocapsae induced in vitro with insect haemolymph

Steinernema carpocapsae is an insect parasitic nematode widely used in pest control programs. The efficacy of this nematode in controlling insects has been found to be related to the pathogenicity of the infective stage. In order to study the parasitic mechanisms exhibited by this parasite, a cDNA library of the induced S. carpocapsae parasitic phase was generated. A total of 2500 clones were sequenced and 2180 high-quality ESTs were obtained from this library. Cluster analysis generated a total of 1592 unique sequences including 1393 singletons. About 63% of the unique sequences had significant hits (e</=1e-05) to the non-redundant protein database. The remaining sequences most likely represent putative novel protein coding genes. Comparative analysis identified 377 homologs in C. elegans, 431 in C. briggsae and 75 in other nematodes. Classification of the predicted proteins revealed involvement in diverse cellular, metabolic and extracellular functions. One hundred and nineteen clusters were predicted to encode putative secreted proteins such as proteases, proteases inhibitors, lectins, saposin-like proteins, acetylcholinesterase, anti-oxidants, and heat-shock proteins, which can possibly have host interactions. This dataset provides a basis for genomic studies towards a better understanding of the events that occur in the parasitic process of this entomopathogenic nematode, including invasion of the insect haemocoelium, adaptations to insect innate immunity and stress responses, and production of virulence factors. The identification of key genes in the parasitic process provides useful tools for the improvement of S. carpocapsae as a biological agent.

Helminth genomics: The implications for human health

More than two billion people (one-third of humanity) are infected with parasitic roundworms or flatworms, collectively known as helminth parasites. These infections cause diseases that are responsible for enormous levels of morbidity and mortality, delays in the physical development of children, loss of productivity among the workforce, and maintenance of poverty. Genomes of the major helminth species that affect humans, and many others of agricultural and veterinary significance, are now the subject of intensive genome sequencing and annotation. Draft genome sequences of the filarial worm Brugia malayi and two of the human schistosomes, Schistosoma japonicum and S. mansoni, are now available, among others. These genome data will provide the basis for a comprehensive understanding of the molecular mechanisms involved in helminth nutrition and metabolism, host-dependent development and maturation, immune evasion, and evolution. They are likely also to predict new potential vaccine candidates and drug targets. In this review, we present an overview of these efforts and emphasize the potential impact and importance of these new findings.

From parasite genomes to one healthy world: Are we having fun yet?

In 1990, the Human Genome Sequencing Project was established. This laid the ground work for an explosion of sequence data that has since followed. As a result of this effort, the first complete genome of an animal, Caenorhabditis elegans was published in 1998. The sequence of Drosophila melanogaster was made available in March, 2000 and in the following year, working drafts of the human genome were generated with the completed sequence (92%) being released in 2003. Recent advancements and next-generation technologies have made sequencing common place and have infiltrated every aspect of biological research, including parasitology. To date, sequencing of 32 apicomplexa and 24 nematode genomes are either in progress or near completion, and over 600k nematode EST and 200k apicomplexa EST submissions fill the databases. However, the winds have shifted and efforts are now refocusing on how best to store, mine and apply these data to problem solving. Herein we tend not to summarize existing X-omics datasets or present new technological advances that promise future benefits. Rather, the information to follow condenses up-to-date-applications of existing technologies to problem solving as it relates to parasite research. Advancements in non-parasite systems are also presented with the proviso that applications to parasite research are in the making.

Transcriptomes and pathways associated with infectivity, survival and immunogenicity in Brugia malayi L3

BACKGROUND

Filarial nematode parasites cause serious diseases such as elephantiasis and river blindness in humans, and heartworm infections in dogs. Third stage filarial larvae (L3) are a critical stage in the life cycle of filarial parasites, because this is the stage that is transmitted by arthropod vectors to initiate infections in mammals. Improved understanding of molecular mechanisms associated with this transition may provide important leads for development of new therapies and vaccines to prevent filarial infections. This study explores changes in gene expression associated with the transition of Brugia malayi third stage larvae (BmL3) from mosquitoes into mammalian hosts and how these changes are affected by radiation. Radiation effects are especially interesting because irradiated L3 induce partial immunity to filarial infections. The underlying molecular mechanisms responsible for the efficacy of such vaccines are unkown.

RESULTS

Expression profiles were obtained using a new filarial microarray with 18, 104 64-mer elements. 771 genes were identified as differentially expressed in two-way comparative analyses of the three L3 types. 353 genes were up-regulated in mosquito L3 (L3i) relative to cultured L3 (L3c). These genes are important for establishment of filarial infections in mammalian hosts. Other genes were up-regulated in L3c relative to L3i (234) or irradiated L3 (L3ir) (22). These culture-induced transcripts include key molecules required for growth and development. 165 genes were up-regulated in L3ir relative to L3c; these genes encode highly immunogenic proteins and proteins involved in radiation repair. L3ir and L3i have similar transcription profiles for genes that encode highly immunogenic proteins, antioxidants and cuticle components.

CONCLUSION

Changes in gene expression that normally occur during culture under conditions that support L3 development and molting are prevented or delayed by radiation. This may explain the enhanced immunogenicity of L3ir. Gene Ontology and KEGG analyses revealed altered pathways between L3 types. Energy and "immune pathways" are up-regulated and may be needed for L3i invasion and survival, while growth and development are priorities for L3c. This study has improved our understanding of molecules involved in parasite invasion and immune evasion, potential targets of protective immunity, and molecules required for parasite growth and development.

Advanced in silico analysis of expressed sequence tag (EST) data for parasitic nematodes of major socio-economic importance--fundamental insights toward biotechnological outcomes

Parasitic nematodes infect humans, other animals and plants, and impose a significant public health and economic burden worldwide due to the diseases that they cause. A better understanding of parasite genomes, host-parasite relationships and the molecular biology of parasites themselves will enable the rational development of diagnostic tests and/or safe anti-parasitic compounds, following the functional annotation of parasite genomic sequences. With only a few completely sequenced nematode genomes, expressed sequence tag (EST) datasets provide a low-cost alternative ("poor man's genome") to whole genome sequences and a glimpse of the transcriptome of an organism. EST data require a number of computational methods for their pre-processing, clustering, assembly and annotation to yield biologically relevant information. In this article, we review the steps involved in EST data analysis, the development of new semi-automated bioinformatic pipelines and their application to parasitic nematodes of major socio-economic significance, focused on identifying molecules involved in key biological processes or pathways that might serve as targets for new drugs or vaccines.

Molecular determinants archetypical to the phylum Nematoda

Background

Nematoda diverged from other animals between 600–1,200 million years ago and has become one of the most diverse animal phyla on earth. Most nematodes are free-living animals, but many are parasites of plants and animals including humans, posing major ecological and economical challenges around the world.

Results

We investigated phylum-specific molecular characteristics in Nematoda by exploring over 214,000 polypeptides from 32 nematode species including 27 parasites. Over 50,000 nematode protein families were identified based on primary sequence, including ~10% with members from at least three different species. Nearly 1,600 of the multi-species families did not share homology to Pfam domains, including a total of 758 restricted to Nematoda. Majority of the 462 families that were conserved among both free-living and parasitic species contained members from multiple nematode clades, yet ~90% of the 296 parasite-specific families originated only from a single clade. Features of these protein families were revealed through extrapolation of essential functions from observed RNAi phenotypes in C. elegans, bioinformatics-based functional annotations, identification of distant homology based on protein folds, and prediction of expression at accessible nematode surfaces. In addition, we identified a group of nematode-restricted sequence features in energy-generating electron transfer complexes as potential targets for new chemicals with minimal or no toxicity to the host.

Conclusion

This study identified and characterized the molecular determinants that help in defining the phylum Nematoda, and therefore improved our understanding of nematode protein evolution and provided novel insights for the development of next generation parasite control strategies.

The transcriptomes of the cattle parasitic nematode Ostertagia ostartagi

Ostertagia ostertagi is a gastrointestinal parasitic nematode that affects cattle and leads to a loss of production. In this study, we present the first large-scale genomic survey of O. ostertagi by the analysis of expressed transcripts from three stages of the parasite: third-stage larvae, fourth-stage larvae and adult worms. Using an in silico approach, 2284 genes were identified from over 7000 expressed sequence tags and abundant transcripts were analyzed and characterized by their functional profile. Of the 2284 genes, 66% had similarity to other known or predicted genes while the rest were novel and potentially represent genes specific to the species and/or stages. Furthermore, a subset of the novel proteins were structurally annotated and assigned putative function based on orthologs in Caenorhabditis elegans and corresponding RNA interference phenotypes. Hence, over 70% of the genes were annotated using protein sequences, domains and pathway databases. Differentially expressed transcripts from the two larval stages and their functional profiles were also studied leading to a more detailed understanding of the parasite's life-cycle. The identified transcripts are a valuable resource for genomic studies of O. ostertagi and can facilitate the design of control strategies and vaccine programs.

Sequence mining and transcript profiling to explore cyst nematode parasitism

BACKGROUND

Cyst nematodes are devastating plant parasites that become sedentary within plant roots and induce the transformation of normal plant cells into elaborate feeding cells with the help of secreted effectors, the parasitism proteins. These proteins are the translation products of parasitism genes and are secreted molecular tools that allow cyst nematodes to infect plants.

RESULTS

We present here the expression patterns of all previously described parasitism genes of the soybean cyst nematode, Heterodera glycines, in all major life stages except the adult male. These insights were gained by analyzing our gene expression dataset from experiments using the Affymetrix Soybean Genome Array GeneChip, which contains probeset sequences for 6,860 genes derived from preparasitic and parasitic H. glycines life stages. Targeting the identification of additional H. glycines parasitism-associated genes, we isolated 633 genes encoding secretory proteins using algorithms to predict secretory signal peptides. Furthermore, because some of the known H. glycines parasitism proteins have strongest similarity to proteins of plants and microbes, we searched for predicted protein sequences that showed their highest similarities to plant or microbial proteins and identified 156 H. glycines genes, some of which also contained a signal peptide. Analyses of the expression profiles of these genes allowed the formulation of hypotheses about potential roles in parasitism. This is the first study combining sequence analyses of a substantial EST dataset with microarray expression data of all major life stages (except adult males) for the identification and characterization of putative parasitism-associated proteins in any parasitic nematode.

CONCLUSION

We have established an expression atlas for all known H. glycines parasitism genes. Furthermore, in an effort to identify additional H. glycines genes with putative functions in parasitism, we have reduced the currently known 6,860 H. glycines genes to a pool of 788 most promising candidate genes (including known parasitism genes) and documented their expression profiles. Using our approach to pre-select genes likely involved in parasitism now allows detailed functional analyses in a manner not feasible for larger numbers of genes. The generation of the candidate pool described here is an important enabling advance because it will significantly facilitate the unraveling of fascinating plant-animal interactions and deliver knowledge that can be transferred to other pathogen-host systems. Ultimately, the exploration of true parasitism genes verified from the gene pool delineated here will identify weaknesses in the nematode life cycle that can be exploited by novel anti-nematode efforts.

Strategies for undertaking expressed sequence tag (EST) projects

Complementary DNA (cDNA) sequencing can be used to sample an organism's transcriptome, and the generated EST sequences can be used for a variety of purposes. They are especially important for enhancing the utility of a genome sequence or for providing a gene catalog for a genome that has not or will not be sequenced. In planning and executing a cDNA project, several criteria must be considered. One should clearly define the project purpose, including organism tissue(s) choice, whether those tissues should be pooled, ability to acquire adequate amounts of clean and well-preserved tissue, choice of type(s) of library, and construction of a library (or libraries) that is compatible with project goals. In addition, one must possess the skills to construct the library (or libraries), keeping in mind the number of clones that will be necessary to meet the project requirements. If one is inexperienced in cDNA library construction, it might be wise to outsource the library production and/or sequence and analysis to a sequencing center or to a company that specializes in those activities. One should also be aware that new sequencing platforms are being marketed that may offer simpler protocols that can produce cDNA data in a more rapid and economical manner. Of course, the bioinformatics tools will have to be in place to de-convolute and aid in data analysis for these newer technologies. Possible funding sources for these projects include well-justified grant proposals, private funding, and/or collaborators with available funds.

Role of horizontal gene transfer in the evolution of plant parasitism among nematodes

Horizontal gene transfer (HGT) implies the non-sexual exchange of genetic material between species - in some cases even across kingdoms. Although common among Bacteria and Archaea, HGTs from pro- to eukaryotes and between eukaryotes were thought to be extremely rare. Recent studies on intracellular bacteria and their hosts seriously question this view. Recipient organisms could benefit from HGT as new gene packages could allow them to broaden or change their diet, colonize new habitats, or survive conditions that previously would have been lethal.About a decade ago, plant parasitic nematodes were shown to produce and secrete cellulases. Prior to this, animals were thought to fully depend on microbial symbionts for the breakdown of plant cell walls. This discovery prompted Keen and Roberts (1) to hypothesize that the ability of nematodes to parasitize plants was acquired by HGT from soil bacteria to (ancestral) bacterivorous nematodes. Since the identification of the first nematode cellulases, many more plant cell wall-degrading enzymes (CWDE) have been identified in a range of plant parasitic nematode species.Here we discuss a number of criteria that can be used to underpin an HGT claim. HGT requires close physical contact between donor and recipient, and this could be achieved in, for example, a symbiont-host, or a trophic relationship. The former type of relationship was indeed shown to potentially result in the transfer of genetic material (e.g., Brugia malayi and Wolbachia). However, currently known endosymbionts of nematodes may not be the source of CWDEs. Remarkably, all cellulases discovered so far within the order Tylenchida belong to a single glycoside hydrolase family (GHF5). A range of soil bacteria harbours GHF5 cellulases, but of course nothing can be said about the gene content of soil bacteria at the time HGT took place (if at all). We suggest that characterisation of cellulases (and other CWDEs) and their genomic organisation in more basal (facultative) plant parasitic Tylenchida is needed to find out if CWDEs were indeed acquired via HGT from bacteria. A more complete picture about the evolution of CWDEs among plant parasitic Tylenchida will require a detailed characterisation of two - so far - fully unexplored basal suborders, Tylenchina and Criconematina. Finally, we performed a computational high-throughput identification of potential HGT candidates (including ones unrelated to CWDEs) in plant parasitic nematodes using a genomics approach.

NemaPath: online exploration of KEGG-based metabolic pathways for nematodes

Background

Nematode.net is a web-accessible resource for investigating gene sequences from parasitic and free-living nematode genomes. Beyond the well-characterized model nematode C. elegans, over 500,000 expressed sequence tags (ESTs) and nearly 600,000 genome survey sequences (GSSs) have been generated from 36 nematode species as part of the Parasitic Nematode Genomics Program undertaken by the Genome Center at Washington University School of Medicine. However, these sequencing data are not present in most publicly available protein databases, which only include sequences in Swiss-Prot. Swiss-Prot, in turn, relies on GenBank/Embl/DDJP for predicted proteins from complete genomes or full-length proteins.

Description

Here we present the NemaPath pathway server, a web-based pathway-level visualization tool for navigating putative metabolic pathways for over 30 nematode species, including 27 parasites. The NemaPath approach consists of two parts: 1) a backend tool to align and evaluate nematode genomic sequences (curated EST contigs) against the annotated Kyoto Encyclopedia of Genes and Genomes (KEGG) protein database; 2) a web viewing application that displays annotated KEGG pathway maps based on desired confidence levels of primary sequence similarity as defined by a user. NemaPath also provides cross-referenced access to nematode genome information provided by other tools available on Nematode.net, including: detailed NemaGene EST cluster information; putative translations; GBrowse EST cluster views; links from nematode data to external databases for corresponding synonymous C. elegans counterparts, subject matches in KEGG's gene database, and also KEGG Ontology (KO) identification.

Conclusion

The NemaPath server hosts metabolic pathway mappings for 30 nematode species and is available on the World Wide Web at . The nematode source sequences used for the metabolic pathway mappings are available via FTP , as provided by the Genome Center at Washington University School of Medicine.

Nematode.net update 2008: improvements enabling more efficient data mining and comparative nematode genomics

Nematode.net (http://nematode.net) is a publicly available resource dedicated to the study of parasitic nematodes. In 2000, the Genome Center at Washington University (GC) joined a consortium including the Nematode Genomics group in Edinburgh, and the Pathogen Sequencing Unit of the Sanger Institute to generate expressed sequence tags (ESTs) as an inexpensive and efficient solution for gene discovery in parasitic nematodes. As of 2008 the GC, sampling key parasites of humans, animals and plants, has generated over 500 000 ESTs and 1.2 million genome survey sequences from more than 30 non-Caenorhabditis elegans nematodes. Nematode.net was implemented to offer user-friendly access to data produced by this project. In addition to sequence data, the site hosts: assembled NemaGene clusters in GBrowse views characterizing composition and protein homology, functional Gene Ontology annotations presented via the AmiGO browser, KEGG-based graphical display of NemaGene clusters mapped to metabolic pathways, codon usage tables, NemFam protein families which represent conserved nematode-restricted coding sequences not found in public protein databases, a web-based WU-BLAST search tool that allows complex querying and other assorted resources. The primary aim of Nematode.net is the dissemination of this diverse collection of information to the broader scientific community in a way that is useful, consistent, centralized and enduring.

Exploring transcriptional conservation between Ancylostoma caninum and Haemonchus contortus by oligonucleotide microarray and bioinformatic analyses

In this study, we identified, using an established oligonucleotide microarray platform for the parasitic nematode Haemonchus contortus, transcripts that are 'conserved' between serum-activated and non-activated L3s of Ancylostoma caninum (aL3 and L3, respectively) and H. contortus by cross-species hybridization (CSH) at high stringency and conducted extensive bioinformatic analyses of the cross-hybridizing expressed sequence tags (ESTs). The microarray analysis revealed significant differential hybridization between aL3 and L3 for 32 molecules from A. caninum, of which 29 were shown to have homologues/orthologues in the free-living nematode Caenorhabditis elegans and/or A. caninum and the other three molecules had no homologues in current gene databases. 'Non-wildtype' RNAi phenotypes were recorded for 13 of the C. elegans homologues. A subset of 16 C. elegans homologues/orthologues (i.e. genes abce-1, act-2, C08H9.2, C55F2.1, calu-1, col-181, cpr-6, elo-2, asp-1, K07E3.4, rpn-2, sel-9, T28C12.4, hsb-1, Y57G11C.15 and ZK593.1) were predicted to interact genetically with a total of 156 (range 1-88) other genes. Gene ontology (GO) analysis of the interacting genes revealed that the most common subcategories were signal transduction (7%), intracellular protein transport and glycolysis (6.2%) within 'biological process'; nuclear (25.7%) and intracellular (19.8%) within 'cellular component'; and ATP-binding (14.4%) and protein-binding (8.4%) within 'molecular function'. The potential roles of key molecules in the two blood-feeding parasitic nematodes are discussed in relation to the known roles of their homologues/orthologues in C. elegans. The CSH approach used may provide a tool for the screening of genes conserved across a range of different taxa of parasites for which DNA microarray platforms are not available.

Intestinal transcriptomes of nematodes: comparison of the parasites Ascaris suum and Haemonchus contortus with the free-living Caenorhabditis elegans

Background

The nematode intestine is a major organ responsible for nutrient digestion and absorption; it is also involved in many other processes, such as reproduction, innate immunity, stress responses, and aging. The importance of the intestine as a target for the control of parasitic nematodes has been demonstrated. However, the lack of detailed knowledge on the molecular and cellular functions of the intestine and the level of its conservation across nematodes has impeded breakthroughs in this application.

Methods and Findings

As part of an extensive effort to investigate various transcribed genomes from Ascaris suum and Haemonchus contortus, we generated a large collection of intestinal sequences from parasitic nematodes by identifying 3,121 A. suum and 1,755 H. contortus genes expressed in the adult intestine through the generation of expressed sequence tags. Cross-species comparisons to the intestine of the free-living C. elegans revealed substantial diversification in the adult intestinal transcriptomes among these species, suggesting lineage- or species-specific adaptations during nematode evolution. In contrast, significant conservation of the intestinal gene repertories was also evident, despite the evolutionary distance of ∼350 million years separating them. A group of 241 intestinal protein families (IntFam-241), each containing members from all three species, was identified based on sequence similarities. These conserved proteins accounted for ∼20% of the sampled intestinal transcriptomes from the three nematodes and are proposed to represent conserved core functions in the nematode intestine. Functional characterizations of the IntFam-241 suggested important roles in molecular functions such as protein kinases and proteases, and biological pathways of carbohydrate metabolism, energy metabolism, and translation. Conservation in the core protein families was further explored by extrapolating observable RNA interference phenotypes in C. elegans to their parasitic counterparts.

Conclusions

Our study has provided novel insights into the nematode intestine and lays foundations for further comparative studies on biology, parasitism, and evolution within the phylum Nematoda.

Biological properties of the nematode intestine warrant in-depth investigation, the results of which can be utilized in the control of parasitic nematodes that infect humans, livestock, and plants. Both the importance of intestinal antigens from Haemonchus contortus in immunity and the damage to H. contortus intestine by anthelmintic fenbendazole have highlighted the versatility of the intestine as an emerging target. However, biological information regarding fundamental intestinal cell functions and mechanisms is currently limited. Conserved intestinal genes across nematode pathogens could offer molecular targets for broad parasite control. Furthermore, qualitative and quantitative comparisons on intestinal gene expression among species and lineages can identify basic adaptations relative to a critical selective force, the nutrient acquisition. This study begins to identify intestinal cell characteristics that are conserved across representatives of two clades of nematodes (V and III) and further clarifies diversities that likely reflect species- or lineage-specific adaptations. Results consistent with functional data on digestive enzymes from H. contortus and RNAi in Caenorhabditis elegans, as examples, support the potential for the comparative genomics approach to produce practical applications. This study provides a platform on which extensive investigation of intestinal genes and a more comprehensive understanding of the Nematoda can be gained.

Four abundant novel transcript genes from Toxocara canis with unrelated coding sequences share untranslated region tracts implicated in the control of gene expression

The Toxocara canis "abundant novel transcripts" (ant) are four highly expressed products, constituting >18% of ESTs from the infective stage of this widely prevalent nematode parasite. Using 5' RACE, we determined full-length sequences for each ant gene, between 1.8 and 2.8kb. The four genes (termed ant-3, -5, -30 and -34), share no coding sequence similarity, although their 3'UTRs (untranslated regions) are homologous. Predicted ANT-5 and ANT-30 proteins show distant similarity to RNA regulatory proteins, RNA-dependent RNA polymerase and DEAH-box helicase, respectively. Surprisingly, ant-3 appears to be bi-cistronic, encoding two ORFs (ANT-3.1 and -3.2), each with a predicted N-terminal signal sequence. Antibodies raised to recombinant proteins did not react with native parasite products, indicating that protein expression did not accord with transcript abundance. However, antibody reactivity to two gene products (ANT-3.1 and ANT-34) was present in patient sera, suggesting that these proteins are synthesized later in infection. To test whether 3'UTRs may regulate expression, the ant-34 3'UTR sequence was inserted adjacent to enhanced green fluorescent protein (EGFP) for transformation of Caenorhabditis elegans. The ant-34 3'UTR greatly reduced EGFP expression, inhibiting both transcription and translation. We identified a tract in this UTR with significant sequence complementarity to the C. elegans micro-RNA lin-4. While infective stage parasites stockpile high levels of the ant transcripts, we suggest that translation is repressed, possibly by a mechanism involving 3' UTR motifs shared by the four genes.

On the extent and origins of genic novelty in the phylum Nematoda

BACKGROUND

The phylum Nematoda is biologically diverse, including parasites of plants and animals as well as free-living taxa. Underpinning this diversity will be commensurate diversity in expressed genes, including gene sets associated specifically with evolution of parasitism.

METHODS AND FINDINGS

Here we have analyzed the extensive expressed sequence tag data (available for 37 nematode species, most of which are parasites) and define over 120,000 distinct putative genes from which we have derived robust protein translations. Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae, these proteins have been grouped into 65,000 protein families that in turn contain 40,000 distinct protein domains. We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla. Gene loss is common, and in particular we identify nearly 5,000 genes that may have been lost from the lineage leading to the model nematode C. elegans. We find a preponderance of novelty, including 56,000 nematode-restricted protein families and 26,000 nematode-restricted domains. Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty. A number of genes from parasitic species had signatures of horizontal transfer from their host organisms, and parasitic species had a greater proportion of novel, secreted proteins than did free-living ones.

CONCLUSIONS

These classes of genes may underpin parasitic phenotypes, and thus may be targets for development of effective control measures.

Transgenesis and neuronal ablation in parasitic nematodes: revolutionary new tools to dissect host-parasite interactions

Ease of experimental gene transfer into viral and prokaryotic pathogens has made transgenesis a powerful tool for investigating the interactions of these pathogens with the host immune system. Recent advances have made this approach feasible for more complex protozoan parasites. By contrast, the lack of a system for heritable transgenesis in parasitic nematodes has hampered progress toward understanding the development of nematode-specific cellular responses. Recently, however, significant strides towards such a system have been made in several parasitic nematodes, and the possible applications of these in immunological research should now be contemplated. In addition, methods for targeted cell ablation have been successfully adapted from Caenorhabditis elegans methodology and applied to studies of neurobiology and behaviour in Strongyloides stercoralis. Together, these new technical developments offer exciting new tools to interrogate multiple aspects of the host-parasite interaction following nematode infection.

Mining the secretome of the root-knot nematode Meloidogyne chitwoodi for candidate parasitism genes

Parasite proteins secreted at the interface of nematode and host are believed to play an essential role in parasitism. Here, we present an efficient pipeline of bio-informatic algorithms and laboratory experiments to identify candidate parasitism genes within nematode secretomes, i.e. the repertoire of secreted proteins in an organism. We performed our approach on 12 218 expressed sequence tags (ESTs) originating from three life stages of the plant parasitic nematode Meloidogyne chitwoodi--a molecularly unexplored root-knot nematode species. The ESTs from M. chitwoodi were assembled into 5880 contigs and open reading frames translated from the consensus sequences were searched for features of putative signal peptides for protein secretion and trans-membrane regions, resulting in the identification of 398 secretome members. The products of parasitism genes are secreted by a range of organs, including the oesophageal, amphidial and rectal glands, the intestine, and the hypodermis. To localize the site of expression in M. chitwoodi, we subjected the most abundant secretome members to in situ hybridization microscopy. We found hybridization of one tag in the dorsal oesophageal gland, seven in the two subventral oesophageal glands, two in the intestine and one tag hybridized to the tail tip in the proximity of the phasmids. Four sequences showed similarity to putative parasitism genes from other nematode species, whereas seven represented pioneering sequences. Our approach presents an efficient method to identify candidate parasitism genes, which does not require sophisticated cDNA isolation and selection protocols, and can therefore be used as a powerful starting point for the molecular investigation of parasites.

Stage-specific gene expression in Teladorsagia circumcincta (Nematoda: Strongylida) infective larvae and early parasitic stages

Suppression subtractive hybridisation was used to enrich genes expressed in a stage-specific manner in infective, exsheathed L3s (xL3) versus early L4s of the ovine nematode, Teladorsagia circumcincta prior to gene expression profiling by microarray. The 769 cDNA sequences obtained from the xL3-enriched library contained 361 unique sequences, with 292 expressed sequence tags (ESTs) being represented once ("singletons") and 69 sequences which were represented more than once (overlapping and non-overlapping "contigs"). The L4-enriched EST dataset contained 472 unique sequences, with 314 singletons and 158 contigs. Of these 833 sequences, 85% of the xL3 sequences and 86% of the L4 sequences exhibited homology to known genes or ESTs derived from other species of nematode. Quantitative differential expression (P<0.05) was demonstrated for 563 (68%) of the ESTs by microarray. Within the L3-specific dataset, more than 30% of the transcripts represented the enzyme, guanosine-5'-triphosphate (GTP)-cyclohydrolase, which is the first and rate-limiting enzyme of the tetrahydrobiopterin synthesis pathway and may be involved in critical elements of larval development. In L4s, proteolytic enzymes were highly up-regulated, as were collagens and a number of previously characterised secretory proteins, reflecting the rapid growth of these larvae in abomasal glands. Nucleotide sequence data reported in this paper are available in the EMBL, GenBank and DDJB databases under accession numbers AM 743198-AM 744942.