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

A Multi-Omics Database for Parasitic Nematodes and Trematodes

Helminth.net ( www.helminth.net ) is a web-based resource that was launched in 2000 as simply " Nematode.net " to host and investigate gene sequences from nematode genomes. Over the years it has evolved to become the moniker for a collection of databases: Nematode.net and Trematode.net . These databases host information for 73 nematode (roundworms) and 17 trematode (flatworms) species and serve as backbone for a number of tools that allow users to query slices of the data for multifactorial combinations of species-omics properties. Recent focus has been on inclusion of gene and protein expression data, population genomics and cross-kingdom interactions (metagenomics datasets). This chapter describes the website, the available tools and some of the new features.

Exploiting Helminth-Host Interactomes through Big Data

Helminths facilitate their parasitic existence through the production and secretion of different molecules, including proteins. Some helminth proteins can manipulate the host's immune system, a phenomenon that is now being exploited with a view to developing therapeutics for inflammatory diseases. In recent years, hundreds of helminth genomes have been sequenced, but as a community we are still taking baby steps when it comes to identifying proteins that govern host-helminth interactions. The information generated from genomic, immunomic, and proteomic studies, as well as from cutting-edge approaches such as proteogenomics, is leading to a substantial volume of big data that can be utilised to shed light on fundamental biology and provide solutions for the development of bioactive-molecule-based therapeutics.

Profiling the macrofilaricidal effects of flubendazole on adult female Brugia malayi using RNAseq

The use of microfilaricidal drugs for the control of onchocerciasis and lymphatic filariasis (LF) necessitates prolonged yearly dosing. Prospects for elimination or eradication of these diseases would be enhanced by the availability of a macrofilaricidal drug. Flubendazole (FLBZ), a benzimidazole anthelmintic, is an appealing candidate. FLBZ has demonstrated potent macrofilaricidal effects in a number of experimental rodent models and in one human trial. Unfortunately, FLBZ was deemed unsatisfactory for use in mass drug administration campaigns due to its limited oral bioavailability. A new formulation that enables sufficient bioavailability following oral administration could render FLBZ an effective treatment for onchocerciasis and LF. Identification of drug-derived effects is important in ascertaining a dosage regimen which is predicted to be lethal to the parasite in situ. In previous histological studies, exposure to FLBZ induced damage to tissues required for reproduction and survival at pharmacologically relevant concentrations. However, more precise and quantitative indices of drug effects are needed. This study assessed drug effects using a transcriptomic approach to confirm effects observed histologically and to identify genes which were differentially expressed in treated adult female Brugia malayi. Comparative analysis across different concentrations (1 μM and 5 μM) and durations (48 and 120 h) provided an overview of the processes which are affected by FLBZ exposure. Genes with dysregulated expression were consistent with the reproductive effects observed via histology in our previous studies. This study revealed transcriptional changes in genes involved in embryo development. Additionally, significant downregulation was observed in genes encoding cuticle components, which may reflect changes in developing embryos, the adult worm cuticle or both. These data support the hypothesis that FLBZ acts predominantly on rapidly dividing cells, and provides a basis for selecting molecular markers of drug-induced damage which may be of use in predicting efficacious FLBZ regimens.

The Effects of Ivermectin on Brugia malayi Females In Vitro: A Transcriptomic Approach

BACKGROUND

Lymphatic filariasis and onchocerciasis are disabling and disfiguring neglected tropical diseases of major importance in developing countries. Ivermectin is the drug of choice for mass drug administration programs for the control of onchocerciasis and lymphatic filariasis in areas where the diseases are co-endemic. Although ivermectin paralyzes somatic and pharyngeal muscles in many nematodes, these actions are poorly characterized in adult filariae. We hypothesize that paralysis of pharyngeal pumping by ivermectin in filariae could result in deprivation of essential nutrients, especially iron, inducing a wide range of responses evidenced by altered gene expression, changes in metabolic pathways, and altered developmental states in embryos. Previous studies have shown that ivermectin treatment significantly reduces microfilariae release from females within four days of exposure in vivo, while not markedly affecting adult worms. However, the mechanisms responsible for reduced production of microfilariae are poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed transcriptomic profiles from Brugia malayi adult females, an important model for other filariae, using RNAseq technology after exposure in culture to ivermectin at various concentrations (100 nM, 300 nM and 1 μM) and time points (24, 48, 72 h, and 5 days). Our analysis revealed drug-related changes in expression of genes involved in meiosis, as well as oxidative phosphorylation, which were significantly down-regulated as early as 24 h post-exposure. RNA interference phenotypes of the orthologs of these down-regulated genes in C. elegans include "maternal sterile", "embryonic lethal", "larval arrest", "larval lethal" and "sick".

CONCLUSION/SIGNIFICANCE

These changes provide insight into the mechanisms involved in ivermectin-induced reduction in microfilaria output and impaired fertility, embryogenesis, and larval development.

The Effect of In Vitro Cultivation on the Transcriptome of Adult Brugia malayi

BACKGROUND

Filarial nematodes cause serious and debilitating infections in human populations of tropical countries, contributing to an entrenched cycle of poverty. Only one human filarial parasite, Brugia malayi, can be maintained in rodents in the laboratory setting. It has been a widely used model organism in experiments that employ culture systems, the impact of which on the worms is unknown.

METHODOLOGY/PRINCIPAL FINDINGS

Using Illumina RNA sequencing, we characterized changes in gene expression upon in vitro maintenance of adult B. malayi female worms at four time points: immediately upon removal from the host, immediately after receipt following shipment, and after 48 h and 5 days in liquid culture media. The dramatic environmental change and the 24 h time lapse between removal from the host and establishment in culture caused a globally dysregulated gene expression profile. We found a maximum of 562 differentially expressed genes based on pairwise comparison between time points. After an initial shock upon removal from the host and shipping, a few stress fingerprints remained after 48 h in culture and until the experiment was stopped. This was best illustrated by a strong and persistent up-regulation of several genes encoding cuticle collagens, as well as serpins.

CONCLUSIONS/SIGNIFICANCE

These findings suggest that B. malayi can be maintained in culture as a valid system for pharmacological and biological studies, at least for several days after removal from the host and adaptation to the new environment. However, genes encoding several stress indicators remained dysregulated until the experiment was stopped.

A globin domain in a neuronal transmembrane receptor of Caenorhabditis elegans and Ascaris suum: molecular modeling and functional properties

We report the structural and biochemical characterization of GLB-33, a putative neuropeptide receptor that is exclusively expressed in the nervous system of the nematode Caenorhabditis elegans. This unique chimeric protein is composed of a 7-transmembrane domain (7TM), GLB-33 7TM, typical of a G-protein-coupled receptor, and of a globin domain (GD), GLB-33 GD. Comprehensive sequence similarity searches in the genome of the parasitic nematode, Ascaris suum, revealed a chimeric protein that is similar to a Phe-Met-Arg-Phe-amide neuropeptide receptor. The three-dimensional structures of the separate domains of both species and of the full-length proteins were modeled. The 7TM domains of both proteins appeared very similar, but the globin domain of the A. suum receptor surprisingly seemed to lack several helices, suggesting a novel truncated globin fold. The globin domain of C. elegans GLB-33, however, was very similar to a genuine myoglobin-type molecule. Spectroscopic analysis of the recombinant GLB-33 GD showed that the heme is pentacoordinate when ferrous and in the hydroxide-ligated form when ferric, even at neutral pH. Flash-photolysis experiments showed overall fast biphasic CO rebinding kinetics. In its ferrous deoxy form, GLB-33 GD is capable of reversibly binding O2 with a very high affinity and of reducing nitrite to nitric oxide faster than other globins. Collectively, these properties suggest that the globin domain of GLB-33 may serve as a highly sensitive oxygen sensor and/or as a nitrite reductase. Both properties are potentially able to modulate the neuropeptide sensitivity of the neuronal transmembrane receptor.

The map-1 gene family in root-knot nematodes, Meloidogyne spp.: a set of taxonomically restricted genes specific to clonal species

Taxonomically restricted genes (TRGs), i.e., genes that are restricted to a limited subset of phylogenetically related organisms, may be important in adaptation. In parasitic organisms, TRG-encoded proteins are possible determinants of the specificity of host-parasite interactions. In the root-knot nematode (RKN) Meloidogyne incognita, the map-1 gene family encodes expansin-like proteins that are secreted into plant tissues during parasitism, thought to act as effectors to promote successful root infection. MAP-1 proteins exhibit a modular architecture, with variable number and arrangement of 58 and 13-aa domains in their central part. Here, we address the evolutionary origins of this gene family using a combination of bioinformatics and molecular biology approaches. Map-1 genes were solely identified in one single member of the phylum Nematoda, i.e., the genus Meloidogyne, and not detected in any other nematode, thus indicating that the map-1 gene family is indeed a TRG family. A phylogenetic analysis of the distribution of map-1 genes in RKNs further showed that these genes are specifically present in species that reproduce by mitotic parthenogenesis, with the exception of M. floridensis, and could not be detected in RKNs reproducing by either meiotic parthenogenesis or amphimixis. These results highlight the divergence between mitotic and meiotic RKN species as a critical transition in the evolutionary history of these parasites. Analysis of the sequence conservation and organization of repeated domains in map-1 genes suggests that gene duplication(s) together with domain loss/duplication have contributed to the evolution of the map-1 family, and that some strong selection mechanism may be acting upon these genes to maintain their functional role(s) in the specificity of the plant-RKN interactions.

Genomics and transcriptomics across the diversity of the Nematoda

The diversity of biology in nematodes is reflected in the diversity of their genomes. Parasitic species in particular have evolved mechanisms to invade and outwit their hosts, and these offer opportunities for the development of control measures. Genomic analyses can reveal the molecular underpinnings of phenotypes such as parasitism and thus, initiate and support research programmes that explore the manipulation of host and parasite physiologies to achieve favourable outcomes. Wide sampling across nematode diversity allows phylogenetically informed formulation of research hypotheses, identification of core features shared by all species or important evolutionary novelties present in isolated clades. Many nematode species have been investigated through the use of the expressed sequence tag approach, which samples from the transcribed genome. Gene catalogues generated in this way can be explored to reveal the patterns of expression associated with parasitism and candidates for testing as drug targets or vaccine components. Analysis environments, such as NEMBASE facilitate exploitation of these data. The development of new high-throughput DNA-sequencing technologies has facilitated transcriptomic and genomic approaches to parasite biology. Whole genome sequencing offers more complete catalogues of genes and assists a systems approach to phenotype dissection. These efforts are being coordinated through the 959 Nematode Genomes initiative.

MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects

Background

Second-generation sequencing technologies are precipitating major shifts with regards to what kinds of genomes are being sequenced and how they are annotated. While the first generation of genome projects focused on well-studied model organisms, many of today's projects involve exotic organisms whose genomes are largely terra incognita. This complicates their annotation, because unlike first-generation projects, there are no pre-existing 'gold-standard' gene-models with which to train gene-finders. Improvements in genome assembly and the wide availability of mRNA-seq data are also creating opportunities to update and re-annotate previously published genome annotations. Today's genome projects are thus in need of new genome annotation tools that can meet the challenges and opportunities presented by second-generation sequencing technologies.

Results

We present MAKER2, a genome annotation and data management tool designed for second-generation genome projects. MAKER2 is a multi-threaded, parallelized application that can process second-generation datasets of virtually any size. We show that MAKER2 can produce accurate annotations for novel genomes where training-data are limited, of low quality or even non-existent. MAKER2 also provides an easy means to use mRNA-seq data to improve annotation quality; and it can use these data to update legacy annotations, significantly improving their quality. We also show that MAKER2 can evaluate the quality of genome annotations, and identify and prioritize problematic annotations for manual review.

Conclusions

MAKER2 is the first annotation engine specifically designed for second-generation genome projects. MAKER2 scales to datasets of any size, requires little in the way of training data, and can use mRNA-seq data to improve annotation quality. It can also update and manage legacy genome annotation datasets.

Nematode.net update 2011: addition of data sets and tools featuring next-generation sequencing data

Nematode.net (http://nematode.net) has been a publicly available resource for studying nematodes for over a decade. In the past 3 years, we reorganized Nematode.net to provide more user-friendly navigation through the site, a necessity due to the explosion of data from next-generation sequencing platforms. Organism-centric portals containing dynamically generated data are available for over 56 different nematode species. Next-generation data has been added to the various data-mining portals hosted, including NemaBLAST and NemaBrowse. The NemaPath metabolic pathway viewer builds associations using KOs, rather than ECs to provide more accurate and fine-grained descriptions of proteins. Two new features for data analysis and comparative genomics have been added to the site. NemaSNP enables the user to perform population genetics studies in various nematode populations using next-generation sequencing data. HelmCoP (Helminth Control and Prevention) as an independent component of Nematode.net provides an integrated resource for storage, annotation and comparative genomics of helminth genomes to aid in learning more about nematode genomes, as well as drug, pesticide, vaccine and drug target discovery. With this update, Nematode.net will continue to realize its original goal to disseminate diverse bioinformatic data sets and provide analysis tools to the broad scientific community in a useful and user-friendly manner.

Characterization of a secreted macrophage migration inhibitory factor homologue of the parasitic nematode Strongyloides acting at the parasite-host cell interface

Strongyloidiasis is a tropical parasitosis characterized by an alternation between free-living and parasitic stages, and by long-term infection via autoinfection. Since invasion and evasion processes of helminth parasites are substantially attained by the involvement of excretory-secretory products, we identified and characterized the 13.5 kDa macrophage migration inhibitory factor (MIF)-like protein in Strongyloides ratti. Sra-MIF is mainly secreted from the infective stage larvae (iL3), while the transcript was found at lower levels in parasitic and free-living females. Sequence analysis of the full-length cDNA showed the highest homology to the human pathogen Strongyloides stercoralis, and both are related to the MIF type-2. Unlike other mif genes, the Sra-mif includes no intron. The protein was recombinantly expressed in Escherichia coli and purified. Sra-MIF exhibited no in vitro tautomerase activity. The exposure of Sra-MIF to the host immune system is confirmed by high IgG reactivities found in the hosts' sera following infection or immunization. Flow cytometric analysis indicated the binding of Sra-MIF to the monocytes/macrophage lineage but not to peripheral lymphocytes. After exposure to Sra-MIF, monocytes released IL-10 but not TNF-alpha suggesting the involvement of the secreted parasite MIF in host immune responses.

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.

Stage-specific excretory-secretory small heat shock proteins from the parasitic nematode Strongyloides ratti--putative links to host's intestinal mucosal defense system

In a search for molecules involved in the interaction between intestinal nematodes and mammalian mucosal host cells, we performed MS to identify excretory-secretory proteins from Strongyloides ratti. In the excretory-secretory proteins of the parasitic female stage, we detected, in addition to other peptides, peptides homologous with the Caenorhabditis elegans heat shock protein (HSP)-17, named Sra-HSP-17.1 (∼ 19 kDa) and Sra-HSP-17.2 (∼ 18 kDa), with 49% amino acid identity. The full-length cDNAs (483 bp and 474 bp, respectively) were identified, and the genomic organization was analyzed. To allow further characterization, the proteins were recombinantly expressed and purified. Profiling of transcription by quantitative real-time-PCR and of protein by ELISA in various developmental stages revealed parasitic female-specific expression. Sequence analyses of both the DNA and amino acid sequences showed that the two proteins share a conserved α-crystallin domain and variable N-terminals. The Sra-HSP-17s showed the highest homology with the deduced small HSP sequence of the human pathogen Strongyloides stercoralis. We observed strong immunogenicity of both proteins, leading to strong IgG responses following infection of rats. Flow cytometric analysis indicated the binding of Sra-HSP-17s to the monocyte-macrophage lineage but not to peripheral lymphocytes or neutrophils. A rat intestinal epithelial cell line showed dose-dependent binding to Sra-HSP-17.1, but not to Sra-HSP-17.2. Exposed monocytes released interleukin-10 but not tumor necrosis factor-α in response to Sra-HSP-17s, suggesting the possible involvement of secreted female proteins in host immune responses.

Computational archaeology of the Pristionchus pacificus genome reveals evidence of horizontal gene transfers from insects

Background

The recent sequencing of nematode genomes has laid the basis for comparative genomics approaches to study the impact of horizontal gene transfer (HGT) on the adaptation to new environments and the evolution of parasitism. In the beetle associated nematode Pristionchus pacificus HGT events were found to involve cellulase genes of microbial origin and Diapausin genes that are known from beetles, but not from other nematodes. The insect-to-nematode horizontal transfer is of special interest given that P. pacificus shows a tight association with insects.

Results

In this study we utilized the observation that horizontally transferred genes often exhibit codon usage patterns more similar to that of the donor than that of the acceptor genome. We introduced GC-normalized relative codon frequencies as a measure to detect characteristic features of P. pacificus orphan genes that show no homology to other nematode genes. We found that atypical codon usage is particularly prevalent in P. pacificus orphans. By comparing codon usage profiles of 71 species, we detected the most significant enrichment in insect-like codon usage profiles. In cross-species comparisons, we identified 509 HGT candidates that show a significantly higher similarity to insect-like profiles than genes with nematode homologs. The most abundant gene family among these genes are non-LTR retrotransposons. Speculating that retrotransposons might have served as carriers of foreign genetic material, we found a significant local clustering tendency of orphan genes in the vicinity of retrotransposons.

Conclusions

Our study combined codon usage bias, phylogenetic analysis, and genomic colocalization into a general picture of the computational archaeology of the P. pacificus genome and suggests that a substantial fraction of the gene repertoire is of insect origin. We propose that the Pristionchus-beetle association has facilitated HGT and discuss potential vectors of these events.

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.

RNAi effector diversity in nematodes

While RNA interference (RNAi) has been deployed to facilitate gene function studies in diverse helminths, parasitic nematodes appear variably susceptible. To test if this is due to inter-species differences in RNAi effector complements, we performed a primary sequence similarity survey for orthologs of 77 Caenorhabditis elegans RNAi pathway proteins in 13 nematode species for which genomic or transcriptomic datasets were available, with all outputs subjected to domain-structure verification. Our dataset spanned transcriptomes of Ancylostoma caninum and Oesophagostomum dentatum, and genomes of Trichinella spiralis, Ascaris suum, Brugia malayi, Haemonchus contortus, Meloidogyne hapla, Meloidogyne incognita and Pristionchus pacificus, as well as the Caenorhabditis species C. brenneri, C. briggsae, C. japonica and C. remanei, and revealed that: (i) Most of the C. elegans proteins responsible for uptake and spread of exogenously applied double stranded (ds)RNA are absent from parasitic species, including RNAi-competent plant-nematodes; (ii) The Argonautes (AGOs) responsible for gene expression regulation in C. elegans are broadly conserved, unlike those recruited during the induction of RNAi by exogenous dsRNA; (iii) Secondary Argonautes (SAGOs) are poorly conserved, and the nuclear AGO NRDE-3 was not identified in any parasite; (iv) All five Caenorhabditis spp. possess an expanded RNAi effector repertoire relative to the parasitic nematodes, consistent with the propensity for gene loss in nematode parasites; (v) In spite of the quantitative differences in RNAi effector complements across nematode species, all displayed qualitatively similar coverage of functional protein groups. In summary, we could not identify RNAi effector deficiencies that associate with reduced susceptibility in parasitic nematodes. Indeed, similarities in the RNAi effector complements of RNAi refractory and competent nematode parasites support the broad applicability of this research genetic tool in nematodes.

Many organisms regulate gene expression through an RNA interference (RNAi) pathway, first characterized in the nematode Caenorhabditis elegans. This pathway can be triggered experimentally using double-stranded (ds)RNA to selected gene targets, thereby allowing researchers to ‘silence’ individual genes and so investigate their function. It is hoped that this technology will facilitate gene silencing in important parasitic nematodes that impose a considerable health and economic burden on mankind. Unfortunately, differences in RNAi susceptibility have been observed between species. Here we investigated the possibility that differences in the complement of effector proteins involved in the RNAi pathway are responsible for these differences in susceptibility. Our data revealed that most facets of the RNAi pathway are well represented across parasitic nematodes, although there were fewer pathway proteins in other nematodes compared to C. elegans. In contrast, the proteins responsible for uptake and spread of dsRNA are not well represented in parasitic nematodes. However, the importance of these differences is undermined by our observation that the protein complements in all the parasites were qualitatively similar, regardless of RNAi-susceptibility. Clearly, differences in the RNAi pathway of parasitic nematodes do not explain the variations in susceptibility to experimental RNAi.

Targeting protein-protein interactions for parasite control

Finding new drug targets for pathogenic infections would be of great utility for humanity, as there is a large need to develop new drugs to fight infections due to the developing resistance and side effects of current treatments. Current drug targets for pathogen infections involve only a single protein. However, proteins rarely act in isolation, and the majority of biological processes occur via interactions with other proteins, so protein-protein interactions (PPIs) offer a realm of unexplored potential drug targets and are thought to be the next-generation of drug targets. Parasitic worms were chosen for this study because they have deleterious effects on human health, livestock, and plants, costing society billions of dollars annually and many sequenced genomes are available. In this study, we present a computational approach that utilizes whole genomes of 6 parasitic and 1 free-living worm species and 2 hosts. The species were placed in orthologous groups, then binned in species-specific orthologous groups. Proteins that are essential and conserved among species that span a phyla are of greatest value, as they provide foundations for developing broad-control strategies. Two PPI databases were used to find PPIs within the species specific bins. PPIs with unique helminth proteins and helminth proteins with unique features relative to the host, such as indels, were prioritized as drug targets. The PPIs were scored based on RNAi phenotype and homology to the PDB (Protein DataBank). EST data for the various life stages, GO annotation, and druggability were also taken into consideration. Several PPIs emerged from this study as potential drug targets. A few interactions were supported by co-localization of expression in M. incognita (plant parasite) and B. malayi (H. sapiens parasite), which have extremely different modes of parasitism. As more genomes of pathogens are sequenced and PPI databases expanded, this methodology will become increasingly applicable.

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.

The draft genome of the parasitic nematode Trichinella spiralis

Genome-based studies of metazoan evolution are most informative when phylogenetically diverse species are incorporated in the analysis. As such, evolutionary trends within and outside the phylum Nematoda have been less revealing by focusing only on comparisons involving Caenorhabditis elegans. Herein, we present a draft of the 64 megabase nuclear genome of Trichinella spiralis, containing 15,808 protein coding genes. This parasitic nematode is an extant member of a clade that diverged early in the evolution of the phylum enabling identification of archetypical genes and molecular signatures exclusive to nematodes. Comparative analyses support intrachromosomal rearrangements across the phylum, disproportionate numbers of protein family deaths over births in parasitic vs. a non-parasitic nematode, and a preponderance of gene loss and gain events in nematodes relative to Drosophila melanogaster. This sequence and the panphylum characteristics identified herein will advance evolutionary studies and strategies to combat global parasites of humans, food animals and crops.

Gender-associated genes in filarial nematodes are important for reproduction and potential intervention targets

Background

A better understanding of reproductive processes in parasitic nematodes may lead to development of new anthelmintics and control strategies for combating disabling and disfiguring neglected tropical diseases such as lymphatic filariasis and onchocerciasis. Transcriptomatic analysis has provided important new insights into mechanisms of reproduction and development in other invertebrates. We have performed the first genome-wide analysis of gender-associated (GA) gene expression in a filarial nematode to improve understanding of key reproductive processes in these parasites.

Methodology/Principal Findings

The Version 2 Filarial Microarray with 18,104 elements representing ∼85% of the filarial genome was used to identify GA gene transcripts in adult Brugia malayi worms. Approximately 19% of 14,293 genes were identified as GA genes. Many GA genes have potential Caenorhabditis elegans homologues annotated as germline-, oogenesis-, spermatogenesis-, and early embryogenesis- enriched. The potential C. elegans homologues of the filarial GA genes have a higher frequency of severe RNAi phenotypes (such as lethal and sterility) than other C. elegans genes. Molecular functions and biological processes associated with GA genes were gender-segregated. Peptidase, ligase, transferase, regulator activity for kinase and transcription, and rRNA and lipid binding were associated with female GA genes. In contrast, catalytic activity from kinase, ATP, and carbohydrate binding were associated with male GA genes. Cell cycle, transcription, translation, and biological regulation were increased in females, whereas metabolic processes of phosphate and carbohydrate metabolism, energy generation, and cell communication were increased in males. Significantly enriched pathways in females were associated with cell growth and protein synthesis, whereas metabolic pathways such as pentose phosphate and energy production pathways were enriched in males. There were also striking gender differences in environmental information processing and cell communication pathways. Many proteins encoded by GA genes are secreted by Brugia malayi, and these encode immunomodulatory molecules such as antioxidants and host cytokine mimics. Expression of many GA genes has been recently reported to be suppressed by tetracycline, which blocks reproduction in female Brugia malayi. Our localization of GA transcripts in filarial reproductive organs supports the hypothesis that these genes encode proteins involved in reproduction.

Conclusions/Significance

Genome-wide expression profiling coupled with a robust bioinformatics analysis has greatly expanded our understanding of the molecular biology of reproduction in filarial nematodes. This study has highlighted key molecules and pathways associated with reproductive and other biological processes and identified numerous potential candidates for rational drug design to target reproductive processes.

Lymphatic filariasis is a neglected tropical disease that is caused by thread-like parasitic worms that live and reproduce in lymphatic vessels of the human host. There are no vaccines to prevent filariasis, and available drugs are not effective against all stages of the parasite. In addition, recent reports suggest that the filarial nematodes may be developing resistance to key medications. Therefore, there is an urgent need to identify new drug targets in filarial worms. The purpose of this study was to perform a genome-wide analysis of gender-associated gene transcription to improve understanding of key reproductive processes in filarial nematodes. Our results indicate that thousands of genes are differentially expressed in male and female adult worms. Many of those genes are involved in specific reproductive processes such as embryogenesis and spermatogenesis. In addition, expression of some of those genes is suppressed by tetracycline, a drug that leads to sterilization of adult female worms in many filarial species. Thus, gender-associated genes represent priority targets for design of vaccines and drugs that interfere with reproduction of filarial nematodes. Additional work with this type of integrated systems biology approach should lead to important new tools for controlling filarial diseases.

Identification of hookworm DAF-16/FOXO response elements and direct gene targets

BACKGROUND

The infective stage of the parasitic nematode hookworm is developmentally arrested in the environment and needs to infect a specific host to complete its life cycle. The canine hookworm (Ancylostoma caninum) is an excellent model for investigating human hookworm infections. The transcription factor of A. caninum, Ac-DAF-16, which has a characteristic fork head or "winged helix" DNA binding domain (DBD), has been implicated in the resumption of hookworm development in the host. However, the precise roles of Ac-DAF-16 in hookworm parasitism and its downstream targets are unknown. In the present study, we combined molecular techniques and bioinformatics to identify a group of Ac-DAF-16 binding sites and target genes.

METHODOLOGY/PRINCIPAL FINDINGS

The DNA binding domain of Ac-DAF-16 was used to select genomic fragments by in vitro genomic selection. Twenty four bound genomic fragments were analyzed for the presence of the DAF-16 family binding element (DBE) and possible alternative Ac-DAF-16 bind motifs. The 22 genes linked to these genomic fragments were identified using bioinformatics tools and defined as candidate direct gene targets of Ac-DAF-16. Their developmental stage-specific expression patterns were examined. Also, a new putative DAF-16 binding element was identified.

CONCLUSIONS/SIGNIFICANCE

Our results show that Ac-DAF-16 is involved in diverse biological processes throughout hookworm development. Further investigation of these target genes will provide insights into the molecular basis by which Ac-DAF-16 regulates its downstream gene network in hookworm infection.

Genome-wide analysis reveals novel genes essential for heme homeostasis in Caenorhabditis elegans

Heme is a cofactor in proteins that function in almost all sub-cellular compartments and in many diverse biological processes. Heme is produced by a conserved biosynthetic pathway that is highly regulated to prevent the accumulation of heme—a cytotoxic, hydrophobic tetrapyrrole. Caenorhabditis elegans and related parasitic nematodes do not synthesize heme, but instead require environmental heme to grow and develop. Heme homeostasis in these auxotrophs is, therefore, regulated in accordance with available dietary heme. We have capitalized on this auxotrophy in C. elegans to study gene expression changes associated with precisely controlled dietary heme concentrations. RNA was isolated from cultures containing 4, 20, or 500 µM heme; derived cDNA probes were hybridized to Affymetrix C. elegans expression arrays. We identified 288 heme-responsive genes (hrgs) that were differentially expressed under these conditions. Of these genes, 42% had putative homologs in humans, while genomes of medically relevant heme auxotrophs revealed homologs for 12% in both Trypanosoma and Leishmania and 24% in parasitic nematodes. Depletion of each of the 288 hrgs by RNA–mediated interference (RNAi) in a transgenic heme-sensor worm strain identified six genes that regulated heme homeostasis. In addition, seven membrane-spanning transporters involved in heme uptake were identified by RNAi knockdown studies using a toxic heme analog. Comparison of genes that were positive in both of the RNAi screens resulted in the identification of three genes in common that were vital for organismal heme homeostasis in C. elegans. Collectively, our results provide a catalog of genes that are essential for metazoan heme homeostasis and demonstrate the power of C. elegans as a genetic animal model to dissect the regulatory circuits which mediate heme trafficking in both vertebrate hosts and their parasites, which depend on environmental heme for survival.

Heme is an iron-containing cofactor for proteins involved in many critical cellular processes. However, free heme is toxic to cells, suggesting that heme synthesis, acquisition, and transport is highly regulated. Efforts to understand heme trafficking in multicellular organisms have failed primarily due to the inability to separate the processes of endogenous heme synthesis from heme uptake and transport. Caenorhabditis elegans is unique among model organisms because it cannot synthesize heme but instead eats environmental heme to grow and develop normally. Thus, worms are an ideal genetic animal model to study heme homeostasis. This work identifies a novel list of 288 heme-responsive genes (hrgs) in C. elegans and a number of related genes in humans and medically relevant parasites. Knocking down the function of each of these hrgs reveals roles for several in heme uptake, transport, and detection within the organism. Our study provides insights into metazoan regulation of organismal heme homeostasis. The identification of parasite-specific hrg homologs may permit the selective design and screening of drugs that specifically target heme uptake pathways in parasites without affecting the host. Thus, this work has therapeutic implications for the treatment of human iron deficiency, one of the top ten mortality factors world-wide.

Characterizing Ancylostoma caninum transcriptome and exploring nematode parasitic adaptation

BACKGROUND

Hookworm infection is one of the most important neglected diseases in developing countries, with approximately 1 billion people infected worldwide. To better understand hookworm biology and nematode parasitism, the present study generated a near complete transcriptome of the canine hookworm Ancylostoma caninum to a very high coverage using high throughput technology, and compared it to those of the free-living nematode Caenorhabditis elegans and the parasite Brugia malayi.

RESULTS

The generated transcripts from four developmental stages, infective L3, serum stimulated L3, adult male and adult female, covered 93% of the A. caninum transcriptome. The broad diversity among nematode transcriptomes was confirmed, and an impact of parasitic adaptation on transcriptome diversity was inferred. Intra-population analysis showed that A. caninum has higher coding sequence diversity than humans. Examining the developmental expression profiles of A. caninum revealed major transitions in gene expression from larval stages to adult. Adult males expressed the highest number of selectively expressed genes, but adult female expressed the highest number of selective parasitism-related genes. Genes related to parasitism adaptation and A. caninum specific genes exhibited more expression selectivity while those conserved in nematodes tend to be consistently expressed. Parasitism related genes were expressed more selectively in adult male and female worms. The comprehensive analysis of digital expression profiles along with transcriptome comparisons enabled identification of a set of parasitism genes encoding secretory proteins in animal parasitic nematode.

CONCLUSIONS

This study validated the usage of deep sequencing for gene expression profiling. Parasitic adaptation of the canine hookworm is related to its diversity and developmental dynamics. This comprehensive comparative genomic and expression study substantially improves our understanding of the basic biology and parasitism of hookworms and, is expected, in the long run, to accelerate research toward development of vaccines and novel anthelmintics.

Alternative oxidase in animals: unique characteristics and taxonomic distribution

Alternative oxidase (AOX), a ubiquinol oxidase, introduces a branch point into the respiratory electron transport chain, bypassing complexes III and IV and resulting in cyanide-resistant respiration. Previously, AOX was thought to be limited to plants and some fungi and protists but recent work has demonstrated the presence of AOX in most kingdoms of life, including animals. In the present study we identified AOX in 28 animal species representing nine phyla. This expands the known taxonomic distribution of AOX in animals by 10 species and two phyla. Using bioinformatics we found AOX gene sequences in members of the animal phyla Porifera, Placozoa, Cnidaria, Mollusca, Annelida, Nematoda, Echinodermata, Hemichordata and Chordata. Using reverse-transcriptase polymerase chain reaction (RT-PCR) with degenerate primers designed to recognize conserved regions of animal AOX, we demonstrated that AOX genes are transcribed in several animals from different phyla. An analysis of full-length AOX sequences revealed an amino acid motif in the C-terminal region of the protein that is unique to animal AOXs. Animal AOX also lacks an N-terminal cysteine residue that is known to be important for AOX enzyme regulation in plants. We conclude that the presence of AOX is the ancestral state in animals and hypothesize that its absence in some lineages, including vertebrates, is due to gene loss events.

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.

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.

Systematic analysis of insertions and deletions specific to nematode proteins and their proposed functional and evolutionary relevance

Background

Amino acid insertions and deletions in proteins are considered relatively rare events, and their associations with the evolution and adaptation of organisms are not yet understood. In this study, we undertook a systematic analysis of over 214,000 polypeptides from 32 nematode species and identified insertions and deletions unique to nematode proteins in more than 1000 families and provided indirect evidence that these alterations are linked to the evolution and adaptation of nematodes.

Results

Amino acid alterations in sequences of nematodes were identified by comparison with homologous sequences from a wide range of eukaryotic (metzoan) organisms. This comparison revealed that the proteins inferred from transcriptomic datasets for nematodes contained more deletions than insertions, and that the deletions tended to be larger in length than insertions, indicating a decreased size of the transcriptome of nematodes compared with other organisms. The present findings showed that this reduction is more pronounced in parasitic nematodes compared with the free-living nematodes of the genus Caenorhabditis. Consistent with a requirement for conservation in proteins involved in the processing of genetic information, fewer insertions and deletions were detected in such proteins. On the other hand, more insertions and deletions were recorded for proteins inferred to be involved in the endocrine and immune systems, suggesting a link with adaptation. Similarly, proteins involved in multiple cellular pathways tended to display more deletions and insertions than those involved in a single pathway. The number of insertions and deletions shared by a range of plant parasitic nematodes were higher for proteins involved in lipid metabolism and electron transport compared with other nematodes, suggesting an association between metabolic adaptation and parasitism in plant hosts. We also identified three sizable deletions from proteins found to be specific to and shared by parasitic nematodes, which, given their uniqueness, might serve as target candidates for drug design.

Conclusion

This study illustrates the significance of using comparative genomics approaches to identify molecular elements unique to parasitic nematodes, which have adapted to a particular host organism and mode of existence during evolution. While the focus of this study was on nematodes, the approach has applicability to a wide range of other groups of organisms.