Draft genome of the filarial nematode parasite Brugia malayi

Draft genome of neurotropic nematode parasite Angiostrongylus cantonensis, causative agent of human eosinophilic meningitis

Angiostrongylus cantonensis is a bursate nematode parasite that causes eosinophilic meningitis (or meningoencephalitis) in humans in many parts of the world. The genomic data from A. cantonensis will form a useful resource for comparative genomic and chemogenomic studies to aid the development of diagnostics and therapeutics. We have sequenced, assembled and annotated the genome of A. cantonensis. The genome size is estimated to be ∼260 Mb, with 17,280 genomic scaffolds, 91X coverage, 81.45% for complete and 93.95% for partial score based on CEGMA analysis of genome completeness. The number of predicted genes of ≥300 bp was 17,482. A total of 7737 predicted protein-coding genes of ≥50 amino acids were identified in the assembled genome. Among the proteins of known function, kinases are the most abundant followed by transferases. The draft genome contains 34 excretory-secretory proteins (ES), a minimum of 44 Nematode Astacin (NAS) metalloproteases, 12 Homeobox (HOX) genes, and 30 neurotransmitters. The assembled genome size (260 Mb) is larger than those of Pristionchus pacificus, Caenorhabditis elegans, Necator americanus, Caenorhabditis briggsae, Trichinella spiralis, Brugia malayi and Loa loa, but smaller than Haemonchus contortus and Ascaris suum. The repeat content (25%) is similar to H. contortus. The GC content (41.17%) is lower compared to P. pacificus (42.7%) and H. contortus (43.1%) but higher compared to C. briggsae (37.69%), A. suum (37.9%) and N. americanus (40.2%) while the scaffold N50 is 42,191. This draft genome will facilitate the understanding of many unresolved issues on the parasite and the disorder it causes.

Genes conserved in bilaterians but jointly lost with Myc during nematode evolution are enriched in cell proliferation and cell migration functions

Animals use a stereotypical set of developmental genes to build body architectures of varying sizes and organizational complexity. Some genes are critical to developmental patterning, while other genes are important to physiological control of growth. However, growth regulator genes may not be as important in small-bodied “micro-metazoans” such as nematodes. Nematodes use a simplified developmental strategy of lineage-based cell fate specifications to produce an adult bilaterian body composed of a few hundreds of cells. Nematodes also lost the MYC proto-oncogenic regulator of cell proliferation. To identify additional regulators of cell proliferation that were lost with MYC, we computationally screened and determined 839 high-confidence genes that are conserved in bilaterians/lost in nematodes (CIBLIN genes). We find that 30 % of all CIBLIN genes encode transcriptional regulators of cell proliferation, epithelial-to-mesenchyme transitions, and other processes. Over 50 % of CIBLIN genes are unnamed genes in Drosophila, suggesting that there are many understudied genes. Interestingly, CIBLIN genes include many Myc synthetic lethal (MycSL) hits from recent screens. CIBLIN genes include key regulators of heparan sulfate proteoglycan (HSPG) sulfation patterns, and lysyl oxidases involved in cross-linking and modification of the extracellular matrix (ECM). These genes and others suggest the CIBLIN repertoire services critical functions in ECM remodeling and cell migration in large-bodied bilaterians. Correspondingly, CIBLIN genes are co-expressed with Myc in cancer transcriptomes, and include a preponderance of known determinants of cancer progression and tumor aggression. We propose that CIBLIN gene research can improve our understanding of regulatory control of cellular growth in metazoans.

A histochemical study of the Nras/let-60 activity in filarial nematodes

BACKGROUND

Control and elimination of filarial pathogens is a central focus of major global health efforts directed at parasitic diseases of developing countries. Accomplishment of these goals would be markedly enhanced by the enhanced destruction of the adult stage of filariae. The identification of new, more quantitative biomarkers that correlate with mortality or chemotherapeutic damage to adult filariae, would greatly facilitate, for example, the development of new macrofilaricides.

METHODS

An immunocytochemical approach using an antibody against human Nras was used to identify and detect changes in the nematode homolog let-60 that is associated with cell growth and maintenance. Single Onchocerca volvulus nodules were removed from each of 13 patients treated with ivermectin (as part of a community-wide mass drug administration programme), and from each of 13 untreated individuals; these 26 nodules were stained with the anti-Nras antibody. The localization and degree of positivity of Nras/let-60 staining were assessed subjectively and compared between the two groups; the positivity of staining was also quantified, using image analysis, in a subgroup of these nodules. In addition, the specific morphological association between Nras/let-60 and the Wolbachia endosymbiont present in these parasites was also observed in 4 additional filarial species using an anti-Wolbachia surface protein (WSP) antibody under light and confocal microscopy.

RESULTS

Nras/let-60 is present in many structures within the adult female worms. A statistically significant decrease in the general staining intensity of Nras/let-60 was observed in adult female O. volvulus treated with ivermectin when compared with parasites from untreated patients. Nras/let-60 staining was frequently observed to be co-localized with WSP in O.volvulus, Brugia malayi, Litomosoides sigmodontis and Dirofilaria immitis. Nras/let60 is also present in Onchocerca ochengi.

CONCLUSION

Nras/let-60, as detected by immunocytochemical staining, is decreased in ivermectin-treated adult female O. volvulus relative to untreated control specimens, suggesting a suppressive effect of ivermectin on the overall biochemical activity of these parasites. Co-localization of Nras/let-60 and WSP suggests the possibility that the endosymbiont utilizes this nematode protein as part of a mutualistic relationship. Nras/let60 appears to be a useful biomarker for assessing the health of filariae.

The helminth parasite proteome at the host-parasite interface - Informing diagnosis and control

Helminth parasites are a significant health burden for humans in the developing world and also cause substantial economic losses in livestock production across the world. The combined lack of vaccines for the major human and veterinary helminth parasites in addition to the development of drug resistance to anthelmintics in sheep and cattle mean that controlling helminth infection and pathology remains a challenge. However, recent high throughput technological advances mean that screening for potential drug and vaccine candidates is now easier than in previous decades. A better understanding of the host-parasite interactions occurring during infection and pathology and identifying pathways that can be therapeutically targeted for more effective and 'evolution proof' interventions is now required. This review highlights some of the advances that have been made in understanding the host-parasite interface in helminth infections using studies of the temporal expression of parasite proteins, i.e. the parasite proteome, and discuss areas for potential future research and translation.

Expression of five acetylcholine receptor subunit genes in Brugia malayi adult worms

Acetylcholine receptors (AChRs) are required for body movement in parasitic nematodes and are targets of “classical” anthelmintic drugs such as levamisole and pyrantel and of newer drugs such as tribendimidine and derquantel. While neurotransmission explains the effects of these drugs on nematode movement, their effects on parasite reproduction are unexplained. The levamisole AChR type (L-AChRs) in Caenorhabditis elegans is comprised of five subunits: Cel-UNC-29, Cel-UNC-38, Cel-UNC-63, Cel-LEV-1 and Cel-LEV-8. The genome of the filarial parasite Brugia malayi contains nine AChRs subunits including orthologues of Cel-unc-29, Cel-unc-38, and Cel-unc-63. We performed in situ hybridization with RNA probes to localize the expression of five AChR genes (Bm1_35890-Bma-unc-29, Bm1_20330-Bma-unc-38, Bm1_38195-Bma-unc-63, Bm1_48815-Bma-acr-26 and Bm1_40515-Bma-acr-12) in B. malayi adult worms. Four of these genes had similar expression patterns with signals in body muscle, developing embryos, spermatogonia, uterine wall adjacent to stretched microfilariae, wall of Vas deferens, and lateral cord. Three L-AChR subunit genes (Bma-unc-29, Bma-unc-38 and Bma-unc-63) were expressed in body muscle, which is a known target of levamisole. Bma-acr-12 was co-expressed with these levamisole subunit genes in muscle, and this suggests that its protein product may form receptors with other alpha subunits. Bma-acr-26 was expressed in male muscle but not in female muscle. Strong expression signals of these genes in early embryos and gametes in uterus and testis suggest that AChRs may have a role in nervous system development of embryogenesis and spermatogenesis. This would be consistent with embryotoxic effects of drugs that target these receptors in filarial worms. Our data show that the expression of these receptor genes is tightly regulated with regard to localization in adult worms and developmental stage in embryos and gametes. These results may help to explain the broad effects of drugs that target AChRs in filarial worms.

•

Expression patterns of Brugia malayi AChR subunit genes studied by in situ hybridization.

•

All genes highly expressed in developing embryos and sperm precursors.

•

Highly expressed in the walls of uterus and Vas deferens with mature offspring.

•

Four of five genes expressed in body muscle of adult worms.

•

Expression patterns shed new light on the action of anthelmintics in filarial parasites.

Application of small RNA technology for improved control of parasitic helminths

•

MicroRNAs and siRNAs in helminth post-transcriptional gene regulation are reviewed.

•

Many parasitic helminth miRNAs are unique and developmentally expressed.

•

miRNAs released by parasites have diagnostic potential, particularly for filarial and schistosome spp.

•

Parasite and host miRNAs may regulate immune responses.

•

Improvements to siRNA-mediated gene silencing are important for functional genomics.

Over the last decade microRNAs (miRNAs) and small interfering RNAs (siRNAs) have emerged as important regulators of post-transcriptional gene expression. miRNAs are short, non-coding RNAs that regulate a variety of processes including cancer, organ development and immune function. This class of small RNAs bind with partial complementarity to their target mRNA sequences, most often in the 3′UTR, to negatively regulate gene expression. In parasitic helminths, miRNAs are being increasingly studied for their potential roles in development and host-parasite interactions. The availability of genome data, combined with small RNA sequencing, has paved the way to profile miRNAs expressed at particular developmental stages for many parasitic helminths. While some miRNAs are conserved across species, others appear to be unique to specific parasites, suggesting important roles in adaptation and survival in the host environment. Some miRNAs are released from parasites, in exosomes or in protein complexes, and the potential effects of these on host immune function are being increasingly studied. In addition, release of miRNAs from schistosome and filarial parasites into host plasma can be exploited for the development of specific and sensitive diagnostic biomarkers of infection. Interfering with miRNA function, as well as silencing key components of the pathways they regulate, will progress our understanding of parasite development and provide a novel approach to therapeutic control. RNA interference (RNAi) by siRNAs has proven to be inconsistent in parasitic nematodes. However, the recent successes reported for schistosome and liver fluke RNAi, encourage further efforts to enhance delivery of RNA and improve in vitro culture systems and assays to monitor phenotypic effects in nematodes. These improvements are important for the establishment of reliable functional genomic platforms for novel drug and vaccine development. In this review we focus on the important roles of miRNAs and siRNAs in post-transcriptional gene regulation in veterinary parasitic helminths and the potential value of these in parasite diagnosis and control.

Syntenin is involved in the bacteria clearance response of kuruma shrimp (Marsupenaeus japonicus)

Syntenin is a multifunctional cytosolic adaptor protein that contributes to cell migration, proliferation, attachment, and apoptosis, as well as immune response to virus, in vertebrates. However, the functions of syntenin in the antibacterial response of invertebrates remain unclear. In this study, we identified a syntenin-like gene (MjSyn) from the kuruma shrimp (Marsupenaeus japonicus) and detected its function in the antibacterial immunity of shrimp. The full-length MjSyn was 1223 bp with a 963 bp open reading frame that encodes 320 amino acids. The deduced MjSyn proteins contained two atypical PDZ domains (sequence repeat that was first reported in the postsynaptic density protein or PSD-95, DlgA, and ZO-1 protein), an N-terminal domain, and a C-terminal domain. Reverse transcription (RT)-PCR results showed that MjSyn was expressed in all tested tissues. Quantitative real-time PCR analysis revealed that MjSyn transcripts in the hemocyte, gill, and intestine were significantly induced at various time points after infection with Staphylococcus aureus and Vibrio anguillarum. The knockdown of the expression of MjSyn by RNA interference resulted in a significant decrease in the phagocytic ability and increased bacteria number in vivo of shrimp. Moreover, the expression of MjCnx, a cytoplasma and membrane location lectin chaperone protein, was inhibited in the MjSyn-knocked down shrimp, which indicated a possible calnexin-related way. Thus, the MjSyn participates in the bacterial clearance response of kuruma shrimp, thereby providing new insight into the function of this kind of important adaptor protein.

Phylum-Level Conservation of Regulatory Information in Nematodes despite Extensive Non-coding Sequence Divergence

Gene regulatory information guides development and shapes the course of evolution. To test conservation of gene regulation within the phylum Nematoda, we compared the functions of putative cis-regulatory sequences of four sets of orthologs (unc-47, unc-25, mec-3 and elt-2) from distantly-related nematode species. These species, Caenorhabditis elegans, its congeneric C. briggsae, and three parasitic species Meloidogyne hapla, Brugia malayi, and Trichinella spiralis, represent four of the five major clades in the phylum Nematoda. Despite the great phylogenetic distances sampled and the extensive sequence divergence of nematode genomes, all but one of the regulatory elements we tested are able to drive at least a subset of the expected gene expression patterns. We show that functionally conserved cis-regulatory elements have no more extended sequence similarity to their C. elegans orthologs than would be expected by chance, but they do harbor motifs that are important for proper expression of the C. elegans genes. These motifs are too short to be distinguished from the background level of sequence similarity, and while identical in sequence they are not conserved in orientation or position. Functional tests reveal that some of these motifs contribute to proper expression. Our results suggest that conserved regulatory circuitry can persist despite considerable turnover within cis elements.

To explore the phylogenetic limits of conservation of cis-regulatory elements, we used transgenesis to test the functions of enhancers of four genes from several species spanning the phylum Nematoda. While we found a striking degree of functional conservation among the examined cis elements, their DNA sequences lacked apparent conservation with the C. elegans orthologs. In fact, sequence similarity between C. elegans and the distantly related nematodes was no greater than would be expected by chance. Short motifs, similar to known regulatory sequences in C. elegans, can be detected in most of the cis elements. When tested, some of these sites appear to mediate regulatory function. However, they seem to have originated through motif turnover, rather than to have been preserved from a common ancestor. Our results suggest that gene regulatory networks are broadly conserved in the phylum Nematoda, but this conservation persists despite substantial reorganization of regulatory elements and could not be detected using naïve comparisons of sequence similarity.

A novel member of the let-7 microRNA family is associated with developmental transitions in filarial nematode parasites

Background

Filarial nematodes are important pathogens in the tropics transmitted to humans via the bite of blood sucking arthropod vectors. The molecular mechanisms underpinning survival and differentiation of these parasites following transmission are poorly understood. microRNAs are small non-coding RNA molecules that regulate target mRNAs and we set out to investigate whether they play a role in the infection event.

Results

microRNAs differentially expressed during the early post-infective stages of Brugia pahangi L3 were identified by microarray analysis. One of these, bpa-miR-5364, was selected for further study as it is upregulated ~12-fold at 24 hours post-infection, is specific to clade III nematodes, and is a novel member of the let-7 family, which are known to have key developmental functions in the free-living nematode Caenorhabditis elegans. Predicted mRNA targets of bpa-miR-5364 were identified using bioinformatics and comparative genomics approaches that relied on the conservation of miR-5364 binding sites in the orthologous mRNAs of other filarial nematodes. Finally, we confirmed the interaction between bpa-miR-5364 and three of its predicted targets using a dual luciferase assay.

Conclusions

These data provide new insight into the molecular mechanisms underpinning the transmission of third stage larvae of filarial nematodes from vector to mammal. This study is the first to identify parasitic nematode mRNAs that are verified targets of specific microRNAs and demonstrates that post-transcriptional control of gene expression via stage-specific expression of microRNAs may be important in the success of filarial infection.

Electronic supplementary material

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

A Drosophila-centric view of protein tyrosine phosphatases

Most of our knowledge on protein tyrosine phosphatases (PTPs) is derived from human pathologies and mouse knockout models. These models largely correlate well with human disease phenotypes, but can be ambiguous due to compensatory mechanisms introduced by paralogous genes. Here we present the analysis of the PTP complement of the fruit fly and the complementary view that PTP studies in Drosophila will accelerate our understanding of PTPs in physiological and pathological conditions. With only 44 PTP genes, Drosophila represents a streamlined version of the human complement. Our integrated analysis places the Drosophila PTPs into evolutionary and functional contexts, thereby providing a platform for the exploitation of the fly for PTP research and the transfer of knowledge onto other model systems.

Exome and transcriptome sequencing of Aedes aegypti identifies a locus that confers resistance to Brugia malayi and alters the immune response

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

Parasitic success without sex – the nematode experience

Asexual reproduction is usually considered as an evolutionary dead end, and difficulties for asexual lineages to adapt to a fluctuating environment are anticipated due to the lack of sufficient genetic plasticity. Yet, unlike their sexual congeners, mitotic parthenogenetic root-knot nematode species, Meloidogyne spp., are remarkably widespread and polyphagous, with the ability to parasitize most flowering plants. Although this may reflect in part the short-term stability of agricultural environments, the extreme parasitic success of these clonal species points them as an outstanding evolutionary paradox regarding current theories on the benefits of sex. The discovery that most of the genome of the clonal species M. incognita is composed of pairs of homologous but divergent segments that have presumably been evolving independently in the absence of sexual recombination has shed new light on this evolutionary paradox. Together with recent studies on other biological systems, including the closely related sexual species M. hapla and the ancient asexual bdelloid rotifers, this observation suggests that functional innovation could emerge from such a peculiar genome architecture, which may in turn account for the extreme adaptive capacities of these asexual parasites. Additionally, the higher proportion of transposable elements in M. incognita compared to M. hapla and other nematodes may also be responsible in part for genome plasticity in the absence of sexual reproduction. We foresee that ongoing sequencing efforts should lead soon to a genomic framework involving genetically diverse Meloidogyne species with various different reproductive modes. This will undoubtedly promote the entire genus as a unique and valuable model system to help deciphering the evolution of asexual reproduction in eukaryotes.

Genetic blueprint of the zoonotic pathogen Toxocara canis

Toxocara canis is a zoonotic parasite of major socioeconomic importance worldwide. In humans, this nematode causes disease (toxocariasis) mainly in the under-privileged communities in developed and developing countries. Although relatively well studied from clinical and epidemiological perspectives, to date, there has been no global investigation of the molecular biology of this parasite. Here we use next-generation sequencing to produce a draft genome and transcriptome of T. canis to support future biological and biotechnological investigations. This genome is 317 Mb in size, has a repeat content of 13.5% and encodes at least 18,596 protein-coding genes. We study transcription in a larval, as well as adult female and male stages, characterize the parasite’s gene-silencing machinery, explore molecules involved in development or host–parasite interactions and predict intervention targets. The draft genome of T. canis should provide a useful resource for future molecular studies of this and other, related parasites.

Toxocara canis is a zoonotic parasite of major worldwide socioeconomic importance. Here, the authors sequence the genome and transcriptome of T. canis, and highlight potential mechanisms involved in development and host–parasite interactions that could support the pursuit of new drug interventions.

Identification and ultrastructural characterization of the Wolbachia symbiont in Litomosoides chagasfilhoi

BACKGROUND

Filarial nematodes are arthropod-transmitted parasites of vertebrates that affect more than 150 million people around the world and remain a major public health problem throughout tropical and subtropical regions. Despite the importance of these nematodes, the current treatment strategies are not efficient in eliminating the parasite. The main strategy of control is based on chemotherapy with diethylcarbamazine, albendazole and ivermectin. In the 1970s, it was found that some filarids possess endosymbiotic bacteria that are important for the development, survival and infectivity of the nematodes. These bacteria belong to the genus Wolbachia, which is a widespread and abundant intracellular symbiont in worms. Knowledge about the structure of the bacteria and their relationship with their nematode hosts may allow new perspectives for the control of filarial nematodes.

METHODS

In this study, we used transmission electron microscopy combined with three-dimensional approaches to observe the structure of the endosymbiont of the filarial nematode Litomosoides chagasfilhoi, an experimental model for the study of lymphatic filariasis. In addition, the bacterium was classified based on PCR analyses.

RESULTS

The bacterium was mainly found in the hypodermis and in the female reproductive system in close association with host cell structures, such as the nucleus and endoplasmic reticulum. Our ultrastructural data also showed that the symbiont envelope is composed of two membrane units and is enclosed in a cytoplasmic vacuole, the symbiosome. Molecular data revealed that the bacterium of L. chagasfilhoi shares 100% identity with the Wolbachia endosymbiont of Litomosoides galizai.

CONCLUSIONS

Here we described ultrastructural aspects of the relationship of the Wolbachia with the filarial nematode Litomosoides chagasfilhoi and the findings lead us to consider this relationship as a mutualistic symbiosis.

Mitochondrial genome plasticity among species of the nematode genus Meloidogyne (Nematoda: Tylenchina)

The mitochondrial (mt) genomes of the plant-parasitic root-knot nematodes Meloidogyne arenaria, Meloidogyne enterolobii and Meloidogyne javanica were sequenced and compared with those of three other root-knot nematode species in order to explore the mt genome plasticity within Meloidogyne. The mt genomes of M. arenaria, M. enterolobii and M. javanica are circular, with an estimated size of 18.8, 18.9 and 19.6 kb, respectively. Compared to other nematodes these mt genomes are larger, due to the presence of large non-coding regions. The mt genome architecture within the genus Meloidogyne varied in the position of trn genes and in the position, length and nucleotide composition of non-coding regions. These variations were observed independent of the species' natural environments or reproductive modes. M. enterolobii showed three main non-coding regions whereas Meloidogyne chitwoodi, Meloidogyne incognita, M. javanica and M. arenaria had two non-coding regions, and Meloidogyne graminicola had a unique large non-coding region interrupted by two trn genes. trn genes were positioned in different regions of the mt genomes in M. chitwoodi, M. enterolobii and M. graminicola, whereas the trn gene order was identical between M. arenaria, M. incognita and M. javanica. Importantly, M. graminicola had extra copies of trnV and trnS2. High divergence levels between the two copies of each trn might indicate duplication events followed by random loss and mutations in the anticodon. Tree-based methods based on amino acid sequences of 12 mt protein-coding genes support the monophyly for the tropical and mitotic parthenogenetic species, M. arenaria, M. enterolobii, M. incognita and M. javanica and for a clade that includes the meiotic parthenogenetic species, M. chitwoodi and M. graminicola. A comparison of the mt genome architecture in plant-parasitic nematodes and phylogenetic analyses support that Pratylenchus is the most recent ancestor of root-knot nematodes.

Genome mining offers a new starting point for parasitology research

Parasites including helminthes, protozoa, and medical arthropod vectors are a major cause of global infectious diseases, affecting one-sixth of the world's population, which are responsible for enormous levels of morbidity and mortality important and remain impediments to economic development especially in tropical countries. Prevalent drug resistance, lack of highly effective and practical vaccines, as well as specific and sensitive diagnostic markers are proving to be challenging problems in parasitic disease control in most parts of the world. The impressive progress recently made in genome-wide analysis of parasites of medical importance, including trematodes of Clonorchis sinensis, Opisthorchis viverrini, Schistosoma haematobium, S. japonicum, and S. mansoni; nematodes of Brugia malayi, Loa loa, Necator americanus, Trichinella spiralis, and Trichuris suis; cestodes of Echinococcus granulosus, E. multilocularis, and Taenia solium; protozoa of Babesia bovis, B. microti, Cryptosporidium hominis, Eimeria falciformis, E. histolytica, Giardia intestinalis, Leishmania braziliensis, L. donovani, L. major, Plasmodium falciparum, P. vivax, Trichomonas vaginalis, Trypanosoma brucei and T. cruzi; and medical arthropod vectors of Aedes aegypti, Anopheles darlingi, A. sinensis, and Culex quinquefasciatus, have been systematically covered in this review for a comprehensive understanding of the genetic information contained in nuclear, mitochondrial, kinetoplast, plastid, or endosymbiotic bacterial genomes of parasites, further valuable insight into parasite-host interactions and development of promising novel drug and vaccine candidates and preferable diagnostic tools, thereby underpinning the prevention and control of parasitic diseases.

Characterization of the Ca2+-gated and voltage-dependent K+-channel Slo-1 of nematodes and its interaction with emodepside

The cyclooctadepsipeptide emodepside and its parent compound PF1022A are broad-spectrum nematicidal drugs which are able to eliminate nematodes resistant to other anthelmintics. The mode of action of cyclooctadepsipeptides is only partially understood, but involves the latrophilin Lat-1 receptor and the voltage- and calcium-activated potassium channel Slo-1. Genetic evidence suggests that emodepside exerts its anthelmintic activity predominantly through Slo-1. Indeed, slo-1 deficient Caenorhabditis elegans strains are completely emodepside resistant. However, direct effects of emodepside on Slo-1 have not been reported and these channels have only been characterized for C. elegans and related Strongylida. Molecular and bioinformatic analyses identified full-length Slo-1 cDNAs of Ascaris suum, Parascaris equorum, Toxocara canis, Dirofilaria immitis, Brugia malayi, Onchocerca gutturosa and Strongyloides ratti. Two paralogs were identified in the trichocephalids Trichuris muris, Trichuris suis and Trichinella spiralis. Several splice variants encoding truncated channels were identified in Trichuris spp. Slo-1 channels of trichocephalids form a monophyletic group, showing that duplication occurred after the divergence of Enoplea and Chromadorea. To explore the function of a representative protein, C. elegans Slo-1a was expressed in Xenopus laevis oocytes and studied in electrophysiological (voltage-clamp) experiments. Incubation of oocytes with 1-10 µM emodepside caused significantly increased currents over a wide range of step potentials in the absence of experimentally increased intracellular Ca2+, suggesting that emodepside directly opens C. elegans Slo-1a. Emodepside wash-out did not reverse the effect and the Slo-1 inhibitor verruculogen was only effective when applied before, but not after, emodepside. The identification of several splice variants and paralogs in some parasitic nematodes suggests that there are substantial differences in channel properties among species. Most importantly, this study showed for the first time that emodepside directly opens a Slo-1 channel, significantly improving the understanding of the mode of action of this drug class.

What helminth genomes have taught us about parasite evolution

The genomes of more than 20 helminths have now been sequenced. Here we perform a meta-analysis of all sequenced genomes of nematodes and Platyhelminthes, and attempt to address the question of what are the defining characteristics of helminth genomes. We find that parasitic worms lack systems for surface antigenic variation, instead maintaining infections using their surfaces as the first line of defence against the host immune system, with several expanded gene families of genes associated with the surface and tegument. Parasite excretory/secretory products evolve rapidly, and proteases even more so, with each parasite exhibiting unique modifications of its protease repertoire. Endoparasitic flatworms show striking losses of metabolic capabilities, not matched by nematodes. All helminths do however exhibit an overall reduction in auxiliary metabolism (biogenesis of co-factors and vitamins). Overall, the prevailing pattern is that there are few commonalities between the genomes of independently evolved parasitic worms, with each parasite having undergone specific adaptations for their particular niche.

Recent advances in candidate-gene and whole-genome approaches to the discovery of anthelmintic resistance markers and the description of drug/receptor interactions

Anthelmintic resistance has a great impact on livestock production systems worldwide, is an emerging concern in companion animal medicine, and represents a threat to our ongoing ability to control human soil-transmitted helminths. The Consortium for Anthelmintic Resistance and Susceptibility (CARS) provides a forum for scientists to meet and discuss the latest developments in the search for molecular markers of anthelmintic resistance. Such markers are important for detecting drug resistant worm populations, and indicating the likely impact of the resistance on drug efficacy. The molecular basis of resistance is also important for understanding how anthelmintics work, and how drug resistant populations arise. Changes to target receptors, drug efflux and other biological processes can be involved. This paper reports on the CARS group meeting held in August 2013 in Perth, Australia. The latest knowledge on the development of molecular markers for resistance to each of the principal classes of anthelmintics is reviewed. The molecular basis of resistance is best understood for the benzimidazole group of compounds, and we examine recent work to translate this knowledge into useful diagnostics for field use. We examine recent candidate-gene and whole-genome approaches to understanding anthelmintic resistance and identify markers. We also look at drug transporters in terms of providing both useful markers for resistance, as well as opportunities to overcome resistance through the targeting of the transporters themselves with inhibitors. Finally, we describe the tools available for the application of the newest high-throughput sequencing technologies to the study of anthelmintic resistance.

The evolution of pentameric ligand-gated ion-channels and the changing family of anthelmintic drug targets

SUMMARY Pentameric ligand-gated ion-channels rapidly transduce synaptic neurotransmitter signals to an electrical response in post-synaptic neuronal or muscle cells and control the neuromusculature of a majority of multicellular animals. A wide range of pharmaceuticals target these receptors including ethanol, nicotine, anti-depressants and other mood regulating drugs, compounds that control pain and mobility and are targeted by a majority of anthelmintic drugs used to control parasitic infection of humans and livestock. Major advances have been made in recent years to our understanding of the structure, function, activity and the profile of compounds that can activate specific receptors. It is becoming clear that these anthelmintic drug targets are not fixed, but differ in significant details from one nematode species to another. Here we review what is known about the evolution of the pentameric ligand-gated ion-channels, paying particular attention to the nematodes, how we can infer the origins of such receptors and understand the factors that determine how they change both over time and from one species to another. Using this knowledge provides a biological framework in which to understand these important drug targets and avenues to identify new receptors and aid the search for new anthelmintic drugs.

Iron necessity: the secret of Wolbachia's success?

The bacterium Wolbachia (order Rickettsiales) is probably the world's most successful vertically-transmitted symbiont, distributed among a staggering 40% of terrestrial arthropod species. Wolbachia has great potential in vector control due to its ability to manipulate its hosts' reproduction and to impede the replication and dissemination of arboviruses and other pathogens within haematophagous arthropods. In addition, the unexpected presence of Wolbachia in filarial nematodes of medical and veterinary importance has provided an opportunity to target the adult worms of Wuchereria bancrofti, Onchocerca volvulus, and Dirofilaria immitis with safe drugs such as doxycycline. A striking feature of Wolbachia is its phenotypic plasticity between (and sometimes within) hosts, which may be underpinned by its ability to integrate itself into several key processes within eukaryotic cells: oxidative stress, autophagy, and apoptosis. Importantly, despite significant differences in the genomes of arthropod and filarial Wolbachia strains, these nexuses appear to lie on a continuum in different hosts. Here, we consider how iron metabolism may represent a fundamental aspect of host homeostasis that is impacted by Wolbachia infection, connecting disparate pathways ranging from the provision of haem and ATP to programmed cell death, aging, and the recycling of intracellular resources. Depending on how Wolbachia and host cells interact across networks that depend on iron, the gradient between parasitism and mutualism may shift dynamically in some systems, or alternatively, stabilise on one or the other end of the spectrum.

Effects of doxycycline on heartworm embryogenesis, transmission, circulating microfilaria, and adult worms in microfilaremic dogs

Tetracycline treatment of animals or humans infected with filariae that harbor Wolbachia endosymbionts blocks further embryogenesis, and existing microfilariae gradually die. This treatment also kills developing larvae and has a slow-kill effect on adult filariae, all presumably due to elimination of the Wolbachia. Also, Dirofilaria immitis microfilariae in blood collected from dogs up to 25 days after the last dose of doxycycline developed to infective L3 that were normal in appearance and motility in mosquitoes but did not continue to develop or migrate normally after subcutaneous (SC) injection into dogs. The present study was designed to determine whether heartworm microfilariae collected at later times after treatment would regain the ability to continue normal development in a dog. The study also was expected to yield valuable data on the effects of treatment on microfilariae and antigen levels and adult worms. The study was conducted in 16 dogs as two separate replicates at different times. A total of five dogs (two in Replicate A and three in Replicate B) infected either by SC injection of L3 or intravenous transplantation of adult heartworms were given doxycycline orally at 10mg/kg twice daily for 30 days, with three untreated controls. Microfilarial counts in the five treated dogs gradually declined during the 12-13 months after treatment initiation. Two dogs were amicrofilaremic before necropsy and three had 13 or fewer microfilariae/ml. Only one treated dog was negative for heartworm antigen before necropsy. Overall, treated dogs generally had fewer live adult heartworms than controls, and most of their live worms were moribund. All three control dogs remained positive for microfilariae and antigen and had many live worms. L3 from mosquitoes fed on blood collected 73-77 or 161-164 days after initiation of doxycycline treatments were injected SC into five dogs. None of the dogs injected with L3 from mosquitoes fed on blood from doxycycline-treated dogs were ever positive for microfilariae or antigen, and none had worms at necropsy; three control dogs were positive for microfilariae and antigen and had many live worms. These data indicate that doxycycline treatment of microfilaremic dogs gradually reduces numbers of microfilariae and blocks further transmission of heartworms. This latter effect should be highly effective in reducing the rate of selection of heartworms with genes that confer resistance to macrocyclic lactone preventives and microfilaricides. The data also suggest that doxycycline has a slow-kill effect on adult heartworms.

Single molecule sequencing and genome assembly of a clinical specimen of Loa loa, the causative agent of loiasis

Background

More than 20% of the world’s population is at risk for infection by filarial nematodes and >180 million people worldwide are already infected. Along with infection comes significant morbidity that has a socioeconomic impact. The eight filarial nematodes that infect humans are Wuchereria bancrofti, Brugia malayi, Brugia timori, Onchocerca volvulus, Loa loa, Mansonella perstans, Mansonella streptocerca, and Mansonella ozzardi, of which three have published draft genome sequences. Since all have humans as the definitive host, standard avenues of research that rely on culturing and genetics have often not been possible. Therefore, genome sequencing provides an important window into understanding the biology of these parasites. The need for large amounts of high quality genomic DNA from homozygous, inbred lines; the availability of only short sequence reads from next-generation sequencing platforms at a reasonable expense; and the lack of random large insert libraries has limited our ability to generate high quality genome sequences for these parasites. However, the Pacific Biosciences single molecule, real-time sequencing platform holds great promise in reducing input amounts and generating sufficiently long sequences that bypass the need for large insert paired libraries.

Results

Here, we report on efforts to generate a more complete genome assembly for L. loa using genetically heterogeneous DNA isolated from a single clinical sample and sequenced on the Pacific Biosciences platform. To obtain the best assembly, numerous assemblers and sequencing datasets were analyzed, combined, and compared. Quiver-informed trimming of an assembly of only Pacific Biosciences reads by HGAP2 was selected as the final assembly of 96.4 Mbp in 2,250 contigs. This results in ~9% more of the genome in ~85% fewer contigs from ~80% less starting material at a fraction of the cost of previous Roche 454-based sequencing efforts.

Conclusions

The result is the most complete filarial nematode assembly produced thus far and demonstrates the utility of single molecule sequencing on the Pacific Biosciences platform for genetically heterogeneous metazoan genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-788) contains supplementary material, which is available to authorized users.

Potential involvement of Brugia malayi cysteine proteases in the maintenance of the endosymbiotic relationship with Wolbachia

•

Bm-cpl-3 and Bm-cpl-6 are involved during development and endosymbiosis.

•

In tetracycline-treated worms Bm-cpl-3 and -6 are regulated in a bimodal pattern.

•

Reduction in Bm-cpl-3 and -6 levels resulted in hindered microfilarial development.

•

Reduction in Bm-cpl-3 and -6 levels resulted in reduced Wolbachia DNA levels.

•

These enzymes might be strong drug target candidates.

Brugia malayi, a parasitic nematode that causes lymphatic filariasis, harbors endosymbiotic intracellular bacteria, Wolbachia, that are required for the development and reproduction of the worm. The essential nature of this endosymbiosis led to the development of anti-Wolbachia chemotherapeutic approaches for the treatment of human filarial infections. Our study is aimed at identifying specific proteins that play a critical role in this endosymbiotic relationship leading to the identification of potential targets in the adult worms. Filarial cysteine proteases are known to be involved in molting and embryogenesis, processes shown to also be Wolbachia dependent. Based on the observation that cysteine protease transcripts are differentially regulated in response to tetracycline treatment, we focused on defining their role in symbiosis. We observe a bimodal regulation pattern of transcripts encoding cysteine proteases when in vitro tetracycline treated worms were examined. Using tetracycline-treated infertile female worms and purified embryos we established that the first peak of the bimodal pattern corresponds to embryonic transcripts while the second takes place within the hypodermis of the adult worms. Localization studies of the native proteins corresponding to Bm-cpl-3 and Bm-cpl-6 indicate that they are present in the area surrounding Wolbachia, and, in some cases, the proteins appear localized within the bacteria. Both proteins were also found in the inner bodies of microfilariae. The possible role of these cysteine proteases during development and endosymbiosis was further characterized using RNAi. Reduction in Bm-cpl-3 and Bm-cpl-6 transcript levels was accompanied by hindered microfilarial development and release, and reduced Wolbachia DNA levels, making these enzymes strong drug target candidates.

The revolution of whole genome sequencing to study parasites

Genome sequencing has revolutionized the way in which we approach biological research from fundamental molecular biology to ecology and epidemiology. In the last 10 years the field of genomics has changed enormously as technology has improved and the tools for genomic sequencing have moved out of a few dedicated centers and now can be performed on bench-top instruments. In this review we will cover some of the key discoveries that were catalyzed by some of the first genome projects and discuss how this field is developing, what the new challenges are and how this may impact on research in the near future.

Operons are a conserved feature of nematode genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and "spliced leader" (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Bioinformatic exploration of RIO protein kinases of parasitic and free-living nematodes

Despite right open reading frame kinases (RIOKs) being essential for life, their functions, substrates and cellular pathways remain enigmatic. In the present study, gene structures were characterised for 26 RIOKs from draft genomes of parasitic and free-living nematodes. RNA-seq transcription profiles of riok genes were investigated for selected parasitic nematodes and showed that these kinases are transcribed in developmental stages that infect their mammalian host. Three-dimensional structural models of Caenorhabditis elegans RIOKs were predicted, and elucidated functional domains and conserved regions in nematode homologs. These findings provide prospects for functional studies of riok genes in C. elegans, and an opportunity for the design and validation of nematode-specific inhibitors of these enzymes in socioeconomic parasitic worms.

Structure of the trehalose-6-phosphate phosphatase from Brugia malayi reveals key design principles for anthelmintic drugs

Parasitic nematodes are responsible for devastating illnesses that plague many of the world's poorest populations indigenous to the tropical areas of developing nations. Among these diseases is lymphatic filariasis, a major cause of permanent and long-term disability. Proteins essential to nematodes that do not have mammalian counterparts represent targets for therapeutic inhibitor discovery. One promising target is trehalose-6-phosphate phosphatase (T6PP) from Brugia malayi. In the model nematode Caenorhabditis elegans, T6PP is essential for survival due to the toxic effect(s) of the accumulation of trehalose 6-phosphate. T6PP has also been shown to be essential in Mycobacterium tuberculosis. We determined the X-ray crystal structure of T6PP from B. malayi. The protein structure revealed a stabilizing N-terminal MIT-like domain and a catalytic C-terminal C2B-type HAD phosphatase fold. Structure-guided mutagenesis, combined with kinetic analyses using a designed competitive inhibitor, trehalose 6-sulfate, identified five residues important for binding and catalysis. This structure-function analysis along with computational mapping provided the basis for the proposed model of the T6PP-trehalose 6-phosphate complex. The model indicates a substrate-binding mode wherein shape complementarity and van der Waals interactions drive recognition. The mode of binding is in sharp contrast to the homolog sucrose-6-phosphate phosphatase where extensive hydrogen-bond interactions are made to the substrate. Together these results suggest that high-affinity inhibitors will be bi-dentate, taking advantage of substrate-like binding to the phosphoryl-binding pocket while simultaneously utilizing non-native binding to the trehalose pocket. The conservation of the key residues that enforce the shape of the substrate pocket in T6PP enzymes suggest that development of broad-range anthelmintic and antibacterial therapeutics employing this platform may be possible.

Here, we describe the structure of trehalose-6-phosphate phosphatase (T6PP) from Brugia malayi. This enzyme is essential to the organism; deletion of the gene encoding T6PP results in toxic accumulation of trehalose 6-phosphate. Structure-guided mutagenesis coupled with kinetic analyses revealed residues important for binding and catalysis. The model for substrate binding suggests a binding mode in which shape complementarity plays a major role. Conservation of binding residues among T6PP orthologs present in pathogenic nematodes and bacteria favors T6PP as a suitable target for broad-range anthelmintic and antibacterial drug design.

Cofactor independent phosphoglycerate mutase of Brugia malayi induces a mixed Th1/Th2 type immune response and inhibits larval development in the host

Lymphatic filariasis is a major debilitating disease, endemic in 72 countries putting more than 1.39 billion people at risk and 120 million are already infected. Despite the significant progress in chemotherapeutic advancements, there is still need for other measures like development of an effective vaccine or discovery of novel drug targets. In this study, structural and immunological characterization of independent phosphoglycerate mutase of filarial parasite Brugia malayi was carried out. Protein was found to be expressed in all major parasite life stages and as an excretory secretory product of adult parasites. Bm-iPGM also reacted to all the categories of human bancroftian patient's sera including endemic normals. In vivo immunological behaviour of protein was determined in immunized BALB/c mice followed by prophylactic analysis in BALB/c mice and Mastomys coucha. Immunization with Bm-iPGM led to generation of a mixed Th1/Th2 type immune response offering 58.2% protection against larval challenge in BALB/c and 65–68% protection in M. coucha. In vitro studies confirmed participation of anti-Bm-iPGM antibodies in killing of B. malayi infective larvae and microfilariae through ADCC mechanism. The present findings reveal potential immunoprotective nature of Bm-iPGM advocating its worth as an antifilarial vaccine candidate.

The evolution of parasitism in Nematoda

Nematodes are abundant and diverse, and include many parasitic species. Molecular phylogenetic analyses have shown that parasitism of plants and animals has arisen at least 15 times independently. Extant nematode species also display lifestyles that are proposed to be on the evolutionary trajectory to parasitism. Recent advances have permitted the determination of the genomes and transcriptomes of many nematode species. These new data can be used to further resolve the phylogeny of Nematoda, and identify possible genetic patterns associated with parasitism. Plant-parasitic nematode genomes show evidence of horizontal gene transfer from other members of the rhizosphere, and these genes play important roles in the parasite-host interface. Similar horizontal transfer is not evident in animal parasitic groups. Many nematodes have bacterial symbionts that can be essential for survival. Horizontal transfer from symbionts to the nematode is also common, but its biological importance is unclear. Over 100 nematode species are currently targeted for sequencing, and these data will yield important insights into the biology and evolutionary history of parasitism. It is important that these new technologies are also applied to free-living taxa, so that the pre-parasitic ground state can be inferred, and the novelties associated with parasitism isolated.

Comparative analysis of the secretome from a model filarial nematode (Litomosoides sigmodontis) reveals maximal diversity in gravid female parasites

Filarial nematodes (superfamily Filarioidea) are responsible for an annual global health burden of ∼6.3 million disability-adjusted life-years, which represents the greatest single component of morbidity attributable to helminths affecting humans. No vaccine exists for the major filarial diseases, lymphatic filariasis and onchocerciasis; in part because research on protective immunity against filariae has been constrained by the inability of the human-parasitic species to complete their lifecycles in laboratory mice. However, the rodent filaria Litomosoides sigmodontis has become a popular experimental model, as BALB/c mice are fully permissive for its development and reproduction. Here, we provide a comprehensive analysis of excretory-secretory products from L. sigmodontis across five lifecycle stages and identifications of host proteins associated with first-stage larvae (microfilariae) in the blood. Applying intensity-based quantification, we determined the abundance of 302 unique excretory-secretory proteins, of which 64.6% were present in quantifiable amounts only from gravid adult female nematodes. This lifecycle stage, together with immature microfilariae, released four proteins that have not previously been evaluated as vaccine candidates: a predicted 28.5 kDa filaria-specific protein, a zonadhesin and SCO-spondin-like protein, a vitellogenin, and a protein containing six metridin-like ShK toxin domains. Female nematodes also released two proteins derived from the obligate Wolbachia symbiont. Notably, excretory-secretory products from all parasite stages contained several uncharacterized members of the transthyretin-like protein family. Furthermore, biotin labeling revealed that redox proteins and enzymes involved in purinergic signaling were enriched on the adult nematode cuticle. Comparison of the L. sigmodontis adult secretome with that of the human-infective filarial nematode Brugia malayi (reported previously in three independent published studies) identified differences that suggest a considerable underlying diversity of potential immunomodulators. The molecules identified in L. sigmodontis excretory-secretory products show promise not only for vaccination against filarial infections, but for the amelioration of allergy and autoimmune diseases.

Genome and transcriptome of the porcine whipworm Trichuris suis

Trichuris (whipworm) infects 1 billion people worldwide and causes a disease (trichuriasis) that results in major socioeconomic losses in both humans and pigs. Trichuriasis relates to an inflammation of the large intestine manifested in bloody diarrhea, and chronic disease can cause malnourishment and stunting in children. Paradoxically, Trichuris of pigs has shown substantial promise as a treatment for human autoimmune disorders, including inflammatory bowel disease (IBD) and multiple sclerosis. Here we report whole-genome sequencing at ∼140-fold coverage of adult male and female T. suis and ∼80-Mb draft assemblies. We explore stage-, sex- and tissue-specific transcription of mRNAs and small noncoding RNAs.

Whipworm genome and dual-species transcriptome analyses provide molecular insights into an intimate host-parasite interaction

Matthew Berriman and colleagues report the whole-genome sequences of the human-infective whipworm Trichuris trichiura and the mouse-infective laboratory model Trichuris muris. Their transcriptome analyses and examination of T. muris infection in mice provide insights into host response to infection and potential drug targets for this major soil-transmitted helminth.

Whipworms are common soil-transmitted helminths that cause debilitating chronic infections in man. These nematodes are only distantly related to Caenorhabditis elegans and have evolved to occupy an unusual niche, tunneling through epithelial cells of the large intestine. We report here the whole-genome sequences of the human-infective Trichuris trichiura and the mouse laboratory model Trichuris muris. On the basis of whole-transcriptome analyses, we identify many genes that are expressed in a sex- or life stage–specific manner and characterize the transcriptional landscape of a morphological region with unique biological adaptations, namely, bacillary band and stichosome, found only in whipworms and related parasites. Using RNA sequencing data from whipworm-infected mice, we describe the regulated T helper 1 (T_H1)-like immune response of the chronically infected cecum in unprecedented detail. In silico screening identified numerous new potential drug targets against trichuriasis. Together, these genomes and associated functional data elucidate key aspects of the molecular host-parasite interactions that define chronic whipworm infection.

Characterization of the reniform nematode genome by shotgun sequencing

The reniform nematode (RN), a major agricultural pest particularly on cotton in the United States, is among the major plant-parasitic nematodes for which limited genomic information exists. In this study, over 380 Mb of sequence data were generated from pooled DNA of four adult female RNs and assembled into 67,317 contigs, including 25,904 (38.5%) predicted coding contigs and 41,413 (61.5%) noncoding contigs. Most of the characterized repeats were of low complexity (88.9%), and 0.9% of the contigs matched with 53.2% of GenBank ESTs. The most frequent Gene Ontology (GO) terms for molecular function and biological process were protein binding (32%) and embryonic development (20%). Further analysis showed that 741 (1.1%), 94 (0.1%), and 169 (0.25%) RN genomic contigs matched with 1328 (13.9%), 1480 (5.4%), and 1330 (7.4%) supercontigs of Meloidogyne incognita, Brugia malayi, and Pristionchus pacificus, respectively. Chromosome 5 of Caenorhabditis elegans had the highest number of hits to the RN contigs. Seven putative detoxification genes and three carbohydrate-active enzymes (CAZymes) involved in cell wall degradation were studied in more detail. Additionally, kinases, G protein-coupled receptors, and neuropeptides functioning in physiological, developmental, and regulatory processes were identified in the RN genome.

Enzymology of the nematode cuticle: A potential drug target?

All nematodes possess an external structure known as the cuticle, which is crucial for their development and survival. This structure is composed primarily of collagen, which is secreted from the underlying hypodermal cells. Extensive studies using the free-living nematode Caenorhabditis elegans demonstrate that formation of the cuticle requires the activity of an extensive range of enzymes. Enzymes are required both pre-secretion, for synthesis of component proteins such as collagen, and post-secretion, for removal of the previous developmental stage cuticle, in a process known as moulting or exsheathment. The excretion/secretion products of numerous parasitic nematodes contain metallo-, serine and cysteine proteases, and these proteases are conserved across the nematode phylum and many are involved in the moulting/exsheathment process. This review highlights the enzymes required for cuticle formation, with a focus on the post-secretion moulting events. Where orthologues of the C. elegans enzymes have been identified in parasitic nematodes these may represent novel candidate targets for future drug/vaccine development.

Hc-daf-2 encodes an insulin-like receptor kinase in the barber's pole worm, Haemonchus contortus, and restores partial dauer regulation

Infective L3s (iL3s) of parasitic nematodes share common behavioural, morphological and developmental characteristics with the developmentally arrested (dauer) larvae of the free-living nematode Caenorhabditis elegans. It is proposed that similar molecular mechanisms regulate entry into or exit from the dauer stage in C. elegans, and the transition from free-living to parasitic forms of parasitic nematodes. In C. elegans, one of the key factors regulating the dauer transition is the insulin-like receptor (designated Ce-DAF-2) encoded by the gene Ce-daf-2. However, nothing is known about DAF-2 homologues in most parasitic nematodes. Here, using a PCR-based approach, we identified and characterised a gene (Hc-daf-2) and its inferred product (Hc-DAF-2) in Haemonchus contortus (a socioeconomically important parasitic nematode of ruminants). The sequence of Hc-DAF-2 displays significant sequence homology to insulin receptors (IR) in both vertebrates and invertebrates, and contains conserved structural domains. A sequence encoding an important proteolytic motif (RKRR) identified in the predicted peptide sequence of Hc-DAF-2 is consistent with that of the human IR, suggesting that it is involved in the formation of the IR complex. The Hc-daf-2 gene was transcribed in all life stages of H. contortus, with a significant up-regulation in the iL3 compared with other stages. To compare patterns of expression between Hc-daf-2 and Ce-daf-2, reporter constructs fusing the Ce-daf-2 or Hc-daf-2 promoter to sequence encoding GFP were microinjected into the N2 strain of C. elegans, and transgenic lines were established and examined. Both genes showed similar patterns of expression in amphidial (head) neurons, which relate to sensation and signal transduction. Further study by heterologous genetic complementation in a daf-2-deficient strain of C. elegans (CB1370) showed partial rescue of function by Hc-daf-2. Taken together, these findings provide a first insight into the roles of Hc-daf-2/Hc-DAF-2 in the biology and development of H. contortus, particularly in the transition to parasitism.

The Conqueror Worm: recent advances with cholinergic anthelmintics and techniques excite research for better therapeutic drugs

The following account is based on a review lecture given recently at the British Society of Parasitology. We point out that nematode parasites cause very widespread infections of humans, particularly in economically underdeveloped areas where sanitation and hygiene are not adequate. In the absence of adequate clean water and effective vaccines, control and prophylaxis relies on anthelmintic drugs. Widespread use of anthelmintics to control nematode parasites of animals has given rise to the development of resistance and so there is a concern that similar problems will occur in humans if mass drug administration is continued. Recent research on the cholinergic anthelmintic drugs has renewed enthusiasm for the further development of cholinergic anthelmintics. Here we illustrate the use of three parasite nematode models, Ascaris suum, Oesophagostomum dentatum and Brugia malayi, microfluidic techniques and the Xenopus oocyte expression system for testing and examining the effects of cholinergic anthelmintics. We also show how the combination of derquantel, the selective nematode cholinergic antagonist and abamectin produce increased inhibition of the nicotinic acetylcholine receptors on the nematode body muscle. We are optimistic that new compounds and combinations of compounds can limit the effects of drug resistance, allowing anthelmintics to be continued to be used for effective treatment of human and animal helminth parasites.

Diversity and expression of microRNAs in the filarial parasite, Brugia malayi

Human filarial parasites infect an estimated 120 million people in 80 countries worldwide causing blindness and the gross disfigurement of limbs and genitals. An understanding of RNA-mediated regulatory pathways in these parasites may open new avenues for treatment. Toward this goal, small RNAs from Brugia malayi adult females, males and microfilariae were cloned for deep-sequencing. From ∼30 million sequencing reads, 145 miRNAs were identified in the B. malayi genome. Some microRNAs were validated using the p19 RNA binding protein and qPCR. B. malayi miRNAs segregate into 99 families each defined by a unique seed sequence. Sixty-one of the miRNA families are highly conserved with homologues in arthropods, vertebrates and helminths. Of those miRNAs not highly conserved, homologues of 20 B. malayi miRNA families were found in vertebrates. Nine B. malayi miRNA families appear to be filarial-specific as orthologues were not found in other organisms. The miR-2 family is the largest in B. malayi with 11 members. Analysis of the sequences shows that six members result from a recent expansion of the family. Library comparisons found that 1/3 of the B. malayi miRNAs are differentially expressed. For example, miR-71 is 5–7X more highly expressed in microfilariae than adults. Studies suggest that in C.elegans, miR-71 may enhance longevity by targeting the DAF-2 pathway. Characterization of B. malayi miRNAs and their targets will enhance our understanding of their regulatory pathways in filariads and aid in the search for novel therapeutics.

De novo transcriptome sequencing and analysis of the cereal cyst nematode, Heterodera avenae

The cereal cyst nematode (CCN, Heterodera avenae) is a major pest of wheat (Triticum spp) that reduces crop yields in many countries. Cyst nematodes are obligate sedentary endoparasites that reproduce by amphimixis. Here, we report the first transcriptome analysis of two stages of H. avenae. After sequencing extracted RNA from pre parasitic infective juvenile and adult stages of the life cycle, 131 million Illumina high quality paired end reads were obtained which generated 27,765 contigs with N50 of 1,028 base pairs, of which 10,452 were annotated. Comparative analyses were undertaken to evaluate H. avenae sequences with those of other plant, animal and free living nematodes to identify differences in expressed genes. There were 4,431 transcripts common to H. avenae and the free living nematode Caenorhabditis elegans, and 9,462 in common with more closely related potato cyst nematode, Globodera pallida. Annotation of H. avenae carbohydrate active enzymes (CAZy) revealed fewer glycoside hydrolases (GHs) but more glycosyl transferases (GTs) and carbohydrate esterases (CEs) when compared to M. incognita. 1,280 transcripts were found to have secretory signature, presence of signal peptide and absence of transmembrane. In a comparison of genes expressed in the pre-parasitic juvenile and feeding female stages, expression levels of 30 genes with high RPKM (reads per base per kilo million) value, were analysed by qRT-PCR which confirmed the observed differences in their levels of expression levels. In addition, we have also developed a user-friendly resource, Heterodera transcriptome database (HATdb) for public access of the data generated in this study. The new data provided on the transcriptome of H. avenae adds to the genetic resources available to study plant parasitic nematodes and provides an opportunity to seek new effectors that are specifically involved in the H. avenae-cereal host interaction.

Dual RNA-seq of parasite and host reveals gene expression dynamics during filarial worm-mosquito interactions

BACKGROUND

Parasite biology, by its very nature, cannot be understood without integrating it with that of the host, nor can the host response be adequately explained without considering the activity of the parasite. However, due to experimental limitations, molecular studies of parasite-host systems have been predominantly one-sided investigations focusing on either of the partners involved. Here, we conducted a dual RNA-seq time course analysis of filarial worm parasite and host mosquito to better understand the parasite processes underlying development in and interaction with the host tissue, from the establishment of infection to the development of infective-stage larva.

METHODOLOGY/PRINCIPAL FINDINGS

Using the Brugia malayi-Aedes aegypti system, we report parasite gene transcription dynamics, which exhibited a highly ordered developmental program consisting of a series of cyclical and state-transitioning temporal patterns. In addition, we contextualized these parasite data in relation to the concurrent dynamics of the host transcriptome. Comparative analyses using uninfected tissues and different host strains revealed the influence of parasite development on host gene transcription as well as the influence of the host environment on parasite gene transcription. We also critically evaluated the life-cycle transcriptome of B. malayi by comparing developmental stages in the mosquito relative to those in the mammalian host, providing insight into gene expression changes underpinning the mosquito-borne parasitic lifestyle of this heteroxenous parasite.

CONCLUSIONS/SIGNIFICANCE

The data presented herein provide the research community with information to design wet lab experiments and select candidates for future study to more fully dissect the whole set of molecular interactions of both organisms in this mosquito-filarial worm symbiotic relationship. Furthermore, characterization of the transcriptional program over the complete life cycle of the parasite, including stages within the mosquito, could help devise novel targets for control strategies.

Mitochondrial genomes of Meloidogyne chitwoodi and M. incognita (Nematoda: Tylenchina): comparative analysis, gene order and phylogenetic relationships with other nematodes

Root-knot nematodes (Meloidogyne spp.) are among the most important plant pathogens. In this study, the mitochondrial (mt) genomes of the root-knot nematodes, M. chitwoodi and M. incognita were sequenced. PCR analyses suggest that both mt genomes are circular, with an estimated size of 19.7 and 18.6-19.1kb, respectively. The mt genomes each contain a large non-coding region with tandem repeats and the control region. The mt gene arrangement of M. chitwoodi and M. incognita is unlike that of other nematodes. Sequence alignments of the two Meloidogyne mt genomes showed three translocations; two in transfer RNAs and one in cox2. Compared with other nematode mt genomes, the gene arrangement of M. chitwoodi and M. incognita was most similar to Pratylenchus vulnus. Phylogenetic analyses (Maximum Likelihood and Bayesian inference) were conducted using 78 complete mt genomes of diverse nematode species. Analyses based on nucleotides and amino acids of the 12 protein-coding mt genes showed strong support for the monophyly of class Chromadorea, but only amino acid-based analyses supported the monophyly of class Enoplea. The suborder Spirurina was not monophyletic in any of the phylogenetic analyses, contradicting the Clade III model, which groups Ascaridomorpha, Spiruromorpha and Oxyuridomorpha based on the small subunit ribosomal RNA gene. Importantly, comparisons of mt gene arrangement and tree-based methods placed Meloidogyne as sister taxa of Pratylenchus, a migratory plant endoparasitic nematode, and not with the sedentary endoparasitic Heterodera. Thus, comparative analyses of mt genomes suggest that sedentary endoparasitism in Meloidogyne and Heterodera is based on convergent evolution.

Nematode Hsp90: highly conserved but functionally diverse

Nematodes are amongst the most successful and abundant organisms on the planet with approximately 30 000 species described, although the actual number of species is estimated to be one million or more. Despite sharing a relatively simple and invariant body plan, there is considerable diversity within the phylum. Nematodes have evolved to colonize most ecological niches, and can be free-living or can parasitize plants or animals to the detriment of the host organism. In this review we consider the role of heat shock protein 90 (Hsp90) in the nematode life cycle. We describe studies on Hsp90 in the free-living nematode Caenorhabditis elegans and comparative work on the parasitic species Brugia pahangi, and consider whether a dependence upon Hsp90 can be exploited for the control of parasitic species.

Identification and characterization of alternative splicing in parasitic nematode transcriptomes

Background

Alternative splicing (AS) of mRNA is a vital mechanism for enhancing genomic complexity in eukaryotes. Spliced isoforms of the same gene can have diverse molecular and biological functions and are often differentially expressed across various tissues, times, and conditions. Thus, AS has important implications in the study of parasitic nematodes with complex life cycles. Transcriptomic datasets are available from many species, but data must be revisited with splice-aware assembly protocols to facilitate the study of AS in helminthes.

Methods

We sequenced cDNA from the model worm Caenorhabditis elegans using 454/Roche technology for use as an experimental dataset. Reads were assembled with Newbler software, invoking the cDNA option. Several combinations of parameters were tested and assembled transcripts were verified by comparison with previously reported C. elegans genes and transcript isoforms and with Illumina RNAseq data.

Results

Thoughtful adjustment of program parameters increased the percentage of assembled transcripts that matched known C. elegans sequences, decreased mis-assembly rates (i.e., cis- and trans-chimeras), and improved the coverage of the geneset. The optimized protocol was used to update de novo transcriptome assemblies from nine parasitic nematode species, including important pathogens of humans and domestic animals. Our assemblies indicated AS rates in the range of 20-30%, typically with 2-3 transcripts per AS locus, depending on the species. Transcript isoforms from the nine species were translated and searched for similarity to known proteins and functional domains. Some 21 InterPro domains, including several involved in nucleotide and chromatin binding, were statistically correlated with AS genetic loci. In most cases, the Roche/454 data explored in this study are the only sequences available from the species in question; however, the recently published genome of the human hookworm Necator americanus provided an additional opportunity to validate our results.

Conclusions

Our optimized assembly parameters facilitated the first survey of AS among parasitic nematodes. The nine transcriptome assemblies, their protein translations, and basic annotations are available from Nematode.net as a resource for the research community. These should be useful for studies of specific genes and gene families of interest as well as for curating draft genome assemblies as they become available.

microRNAs of parasitic helminths - Identification, characterization and potential as drug targets

•

Importance of microRNAs in helminth post-transcriptional gene regulation is reviewed.

•

Increasing helminth miRNA data are available from deep sequencing.

•

Some miRNAs are helminth-specific, many are novel to each species.

•

miRNAs may regulate parasite and host gene expression.

•

Uptake of miRNA inhibitors and mimics is feasible for functional analysis.

microRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene regulation. They were first identified in the free-living nematode Caenorhabditis elegans, where the miRNAs lin-4 and let-7 were shown to be essential for regulating correct developmental progression. The sequence of let-7 was subsequently found to be conserved in higher organisms and changes in expression of let-7, as well as other miRNAs, are associated with certain cancers, indicating important regulatory roles. Some miRNAs have been shown to have essential functions, but the roles of many are currently unknown. With the increasing availability of genome sequence data, miRNAs have now been identified from a number of parasitic helminths, by deep sequencing of small RNA libraries and bioinformatic approaches. While some miRNAs are widely conserved in a range of organisms, others are helminth-specific and many are novel to each species. Here we review the potential roles of miRNAs in regulating helminth development, in interacting with the host environment and in development of drug resistance. Use of fluorescently-labeled small RNAs demonstrates uptake by parasites, at least in vitro. Therefore delivery of miRNA inhibitors or mimics has potential to alter miRNA activity, providing a useful tool for probing the roles of miRNAs and suggesting novel routes to therapeutics for parasite control.

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B. trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764bp, 13,771bp and 13,803bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A+T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.

Bad guys turned nice? A critical assessment of Wolbachia mutualisms in arthropod hosts

Wolbachia are the most abundant bacterial endosymbionts among arthropods. Although maternally inherited, they do not conform to the widespread view that vertical transmission inevitably selects for beneficial symbionts. Instead, Wolbachia are notorious for their reproductive parasitism which, although lowering host fitness, ensures their spread. However, even for reproductive parasites it can pay to enhance host fitness. Indeed, there is a recent upsurge of reports on Wolbachia-associated fitness benefits. Therefore, the question arises how such instances of mutualism are related to the phenotypes of reproductive parasitism. Here, we review the evidence of Wolbachia mutualisms in arthropods, including both facultative and obligate relationships, and critically assess their biological relevance. Although many studies report anti-pathogenic effects of Wolbachia, few actually prove these effects to be relevant to field conditions. We further show that Wolbachia frequently have beneficial and detrimental effects at the same time, and that reproductive manipulations and obligate mutualisms may share common mechanisms. These findings undermine the idea of a clear-cut distinction between Wolbachia mutualism and parasitism. In general, both facultative and obligate mutualisms can have a strong, and sometimes unforeseen, impact on the ecology and evolution of Wolbachia and their arthropod hosts. Acknowledging this mutualistic potential might be the key to a better understanding of some unresolved issues in the study of Wolbachia-host interactions.