Search

Background: Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method., Methods: We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow., Results: We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441-445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat., Conclusions: We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such 'new' repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

Background: Signalling pathways have been extensively investigated in the free-living nematode Caenorhabditis elegans, but very little is known about these pathways in parasitic nematodes. Here, we constructed a model for the dauer-associated signalling pathways in an economically highly significant parasitic worm, Haemonchus contortus., Methods: Guided by data and information available for C. elegans, we used extensive genomic and transcriptomic datasets to infer gene homologues in the dauer-associated pathways, explore developmental transcriptomic, proteomic and phosphoproteomic profiles in H. contortus and study selected molecular structures., Results: The canonical cyclic guanosine monophosphate (cGMP), transforming growth factor-β (TGF-β), insulin-like growth factor 1 (IGF-1) and steroid hormone signalling pathways of H. contortus were inferred to represent a total of 61 gene homologues. Compared with C. elegans, H. contortus has a reduced set of genes encoding insulin-like peptides, implying evolutionary and biological divergences between the parasitic and free-living nematodes. Similar transcription profiles were found for all gene homologues between the infective stage of H. contortus and dauer stage of C. elegans. High transcriptional levels for genes encoding G protein-coupled receptors (GPCRs), TGF-β, insulin-like ligands (e.g. ins-1, ins-17 and ins-18) and transcriptional factors (e.g. daf-16) in the infective L3 stage of H. contortus were suggestive of critical functional roles in this stage. Conspicuous protein expression patterns and extensive phosphorylation of some components of these pathways suggested marked post-translational modifications also in the L3 stage. The high structural similarity in the DAF-12 ligand binding domain among nematodes indicated functional conservation in steroid (i.e. dafachronic acid) signalling linked to worm development., Conclusions: Taken together, this pathway model provides a basis to explore hypotheses regarding biological processes and regulatory mechanisms (via particular microRNAs, phosphorylation events and/or lipids) associated with the development of H. contortus and related nematodes as well as parasite-host cross talk, which could aid the discovery of new therapeutic targets.

Background: Toxocara canis is quite closely related to Ascaris suum but its biology is more complex, involving a phase of arrested development (diapause or hypobiosis) in tissues as well as transplacental and transmammary transmission routes. In the present study, we explored and compared dauer-like signalling pathways of T. canis and A. suum to infer which components in these pathways might associate with, or regulate, this added complexity in T. canis., Methods: Guided by information for Caenorhabditis elegans, we bioinformatically inferred and compared components of dauer-like signalling pathways in T. canis and A. suum using genomic and transcriptomic data sets. In these two ascaridoids, we also explored endogenous dafachronic acids (DAs), which are known to be critical in regulating larval developmental processes in C. elegans and other nematodes, by liquid chromatography-mass spectrometry (LC-MS)., Results: Orthologues of C. elegans dauer signalling genes were identified in T. canis (n = 55) and A. suum (n = 51), inferring the presence of a dauer-like signalling pathway in both species. Comparisons showed clear differences between C. elegans and these ascaridoids as well as between T. canis and A. suum, particularly in the transforming growth factor-β (TGF-β) and insulin-like signalling pathways. Specifically, in both A. suum and T. canis, there was a paucity of genes encoding SMAD transcription factor-related protein (daf-3, daf-5, daf-8 and daf-14) and insulin/insulin-like peptide (daf-28, ins-4, ins-6 and ins-7) homologues, suggesting an evolution and adaptation of the signalling pathway in these parasites. In T. canis, there were more orthologues coding for homologues of antagonist insulin-like peptides (Tc-ins-1 and Tc-ins-18), an insulin receptor substrate (Tc-ist-1) and a serine/threonine kinase (Tc-akt-1) than in A. suum, suggesting potentiated functional roles for these molecules in regulating larval diapause and reactivation. A relatively conserved machinery was proposed for DA synthesis in the two ascaridoids, and endogenous Δ4- and Δ7-DAs were detected in them by LC-MS analysis. Differential transcription analysis between T. canis and A. suum suggests that ins-17 and ins-18 homologues are specifically involved in regulating development and migration in T. canis larvae in host tissues., Conclusion: The findings of this study provide a basis for functional explorations of insulin-like peptides, signalling hormones (i.e. DAs) and related nuclear receptors, proposed to link to development and/or parasite-host interactions in T. canis. Elucidating the functional roles of these molecules might contribute to the discovery of novel anthelmintic targets in ascaridoids.

Background: The accurate tracking of Cryptosporidium in faecal, water and/or soil samples in water catchment areas is central to developing strategies to manage the potential risk of cryptosporidiosis transmission to humans. Various PCR assays are used for this purpose. Although some assays achieve specific amplification from Cryptosporidium DNA in animal faecal samples, some do not. Indeed, we have observed non-specificity of some oligonucleotide primers in the small subunit of nuclear ribosomal RNA gene (SSU), which has presented an obstacle to the identification and classification of Cryptosporidium species and genotypes (taxa) from faecal samples., Results: Using a novel bioinformatic approach, we explored all available Cryptosporidium genome sequences for new and diagnostically-informative, multi-copy regions to specifically design oligonucleotide primers in the large subunit of nuclear ribosomal RNA gene (LSU) as a basis for an effective nested PCR-based sequencing method for the identification and/or classification of Cryptosporidium taxa., Conclusion: This newly established PCR, which has high analytical specificity and sensitivity, is now in routine use in our laboratory, together with other assays developed by various colleagues. Although the present bioinformatic workflow used here was for the specific design of primers in nuclear DNA of Cryptosporidium, this approach should be broadly applicable to many other microorganisms.

Background: Babesiosis is a socioeconomically important tick-borne disease of animals (including humans) caused by haemoprotozoan parasites. The severity of babesiosis relates to host and parasite factors, particularly virulence/pathogenicity. Although Babesia bovis is a particularly pathogenic species of cattle, there are species of Babesia of ruminants that have limited pathogenicity. For instance, the operational taxonomic unit Babesia sp. Xinjiang (abbreviated here as Bx) of sheep from China is substantially less virulent/pathogenic than B. bovis is in cattle. Although the reason for this distinctiveness is presently unknown, it is possible that Bx has a reduced ability to adhere to cells or evade/suppress immune responses, which might relate to particular proteins, such as the variant erythrocyte surface antigens (VESAs)., Results: We sequenced and annotated the 8.4 Mb nuclear draft genome of Bx and compared it with those of B. bovis and B. bigemina by synteny analysis; we also investigated the genetic relationship of Bx with selected Babesia species and related apicomplexans for which genomic datasets are available, and explored the VESA complement in Bx., Conclusions: The availability of the Bx genome now provides unique opportunities to elucidate aspects of the molecular biology, biochemistry and physiology of Bx, and to explore the reason(s) for its limited virulence and/or apparent ability to evade immune attack by the host animal. Moreover, the present genomic resource and an in vitro culture system for Bx raises the prospect of establishing a functional genomic platform to explore essential genes as new intervention targets against babesiosis.

Background: Proteins of the cysteine-rich secretory proteins, antigen 5 and pathogenesis-related 1 (CAP) superfamily are recognized or proposed to play roles in parasite development and reproduction, and in modulating host immune attack and infection processes. However, little is known about these proteins for most parasites., Results: In the present study, we explored CAP proteins of Toxocara canis, a socioeconomically important zoonotic roundworm. To do this, we mined and curated transcriptomic and genomic data, predicted and curated full-length protein sequences (n = 28), conducted analyses of these data and studied the transcription of respective genes in different developmental stages of T. canis. In addition, based on information available for Caenorhabditis elegans, we inferred that selected genes (including lon-1, vap-1, vap-2, scl-1, scl-8 and scl-11 orthologs) of T. canis and their interaction partners likely play central roles in this parasite's development and/or reproduction via TGF-beta and/or insulin-like signaling pathways, or via host interactions., Conclusion: In conclusion, this study could provide a foundation to guide future studies of CAP proteins of T. canis and related parasites, and might assist in finding new interventions against diseases caused by these parasites.

Background: Toxocara canis is the causative agent of toxocariasis of humans and other animals. This parasitic nematode (roundworm) has a complex life cycle, in which substantial developmental changes and switches occur. As small non-coding RNAs (sRNAs) are key regulators of gene expression in a wide range of organisms, we explored these RNAs in T. canis to provide a basis for future studies of its developmental biology as well as host interactions and disease at the molecular level., Methods: We conducted high-throughput RNA sequencing and bioinformatic analyses to define sRNAs in individual male and female adults of T. canis., Results: Apart from snRNA and snoRNA, 560 and 619 microRNAs (miRNAs), including 5 and 2 novel miRNAs, were identified in male and female worms, respectively, without piRNAs being detected in either sex. An analysis of transcriptional profiles showed that, of 564 miRNAs predicted as being differentially transcribed between male and female individuals of T. canis, 218 miRNAs were transcribed exclusively in male and 277 in female worms. Functional enrichment analysis predicted that both male and female miRNAs were mainly involved in regulating embryonic morphogenesis, hemidesmosome assembly and genetic information processing. The miRNAs differentially transcribed between the sexes were predicted to be associated with sex determination, embryonic morphogenesis and nematode larval development. The roles of miRNAs were predicted based on gene ontology (GO) and KEGG pathway annotations. The miRNAs Tc-miR-2305 and Tc-miR-6090 are proposed to have roles in reproduction, embryo development and larval development, and Tc-let-7-5p, Tc-miR-34 and Tc-miR-100 appear to be involved in host-parasite interactions. Together with published information from previous studies, some miRNAs (such as Tc-miR-2861, Tc-miR-2881 and Tc-miR-5126) are predicted to represent drug targets and/or associated with drug resistance., Conclusions: This is the first exploration of miRNAs in T. canis, which could provide a basis for fundamental investigations of the developmental biology of the parasite, parasite-host interactions and toxocariasis as well as applied areas, such as the diagnosis of infection/disease, drug target discovery and drug resistance detection.

Background: Ion channels are well characterised in model organisms, principally because of the availability of functional genomic tools and datasets for these species. This contrasts the situation, for example, for parasites of humans and animals, whose genomic and biological uniqueness means that many genes and their products cannot be annotated. As ion channels are recognised as important drug targets in mammals, the accurate identification and classification of parasite channels could provide major prospects for defining unique targets for designing novel and specific anti-parasite therapies. Here, we established a reliable bioinformatic pipeline for the identification and classification of ion channels encoded in the genome of the cancer-causing liver fluke Opisthorchis viverrini, and extended its application to related flatworms affecting humans., Methods: We built an ion channel identification + classification pipeline (called MuSICC), employing an optimised support vector machine (SVM) model and using the Kyoto Encyclopaedia of Genes and Genomes (KEGG) classification system. Ion channel proteins were first identified and grouped according to amino acid sequence similarity to classified ion channels and the presence and number of ion channel-like conserved and transmembrane domains. Predicted ion channels were then classified to sub-family using a SVM model, trained using ion channel features., Results: Following an evaluation of this pipeline (MuSICC), which demonstrated a classification sensitivity of 95.2 % and accuracy of 70.5 % for known ion channels, we applied it to effectively identify and classify ion channels in selected parasitic flatworms., Conclusions: MuSICC provides a practical and effective tool for the identification and classification of ion channels of parasitic flatworms, and should be applicable to a broad range of organisms that are evolutionarily distant from taxa whose ion channels are functionally characterised.

Background: Phosphoinositide-dependent protein kinase-1 (PDK-1), which functions downstream of phosphoinositide 3-kinase (AGE-1) and activates protein kinases of the AGC family, plays critical roles in regulating biology processes, such as metabolism, growth, development and survival. In the free-living nematode Caenorhabditis elegans, PDK-1 is a key component of the insulin-like signalling pathway, regulating the entry into and exit from dauer (arrested development). Although it is proposed that similar molecular mechanisms control the transition from the free-living to the parasitic stages of nematodes, nothing is known about PDK-1 in Haemonchus contortus, a socioeconomically important gastric nematode of ruminants., Methods: Here, we isolated and characterized the pdk-1 gene (Hc-pdk-1) and its inferred product (Hc-PDK-1) from H. contortus. Using in vitro and in vivo methods, we then studied the transcriptional profiles of Hc-pdk-1 and anatomical gene expression patterns of Hc-PDK-1 in different developmental stages of C. elegans., Results: In silico analysis of Hc-PDK-1 displayed conserved functional domains, such as protein kinase and pleckstrin homology (PH) domains and two predicted phosphorylation sites (Thr226/Tyr229), which are crucial for the phosphorylation of downstream signalling. The Hc-pdk-1 gene is transcribed in all of the main developmental stages of H. contortus, with its highest transcription in the infective third-stage larvae (iL3) compared with other stages. Transgene constructs, in which respective promoters were fused to the coding sequence for green fluorescent protein (GFP), were used to transform C. elegans, and to localize and compare the expression of Hc-pdk-1 and Ce-pdk-1. The expression of GFP under the control of the Hc-pdk-1 promoter was localized to the intestine, and head and tail neurons, contrasting somewhat the profile for the C. elegans ortholog, which is expressed in pharynx, intestine and head and tail neurons., Conclusions: This is the first characterization of pdk-1/PDK-1 from a trichostrongyloid nematode. Taken together, the findings from this study provide a first glimpse of the involvement of Hc-pdk-1 in the insulin-like signalling pathway in H. contortus.

Background: Protein kinases regulate a plethora of essential signalling and other biological pathways in all eukaryotic organisms, but very little is known about them in most parasitic nematodes., Methods: Here, we defined, for the first time, the entire complement of protein kinases (kinome) encoded in the barber's pole worm (Haemonchus contortus) through an integrated analysis of transcriptomic and genomic datasets using an advanced bioinformatic workflow., Results: We identified, curated and classified 432 kinases representing ten groups, 103 distinct families and 98 subfamilies. A comparison of the kinomes of H. contortus and Caenorhabditis elegans (a related, free-living nematode) revealed considerable variation in the numbers of casein kinases, tyrosine kinases and Ca(2+)/calmodulin-dependent protein kinases, which likely relate to differences in biology, habitat and life cycle between these worms. Moreover, a suite of kinase genes was selectively transcribed in particular developmental stages of H. contortus, indicating central roles in developmental and reproductive processes. In addition, using a ranking system, drug targets (n = 13) and associated small-molecule effectors (n = 1517) were inferred., Conclusions: The H. contortus kinome will provide a useful resource for fundamental investigations of kinases and signalling pathways in this nematode, and should assist future anthelmintic discovery efforts; this is particularly important, given current drug resistance problems in parasitic nematodes.

Background: Phosphoinositide 3-kinases (PI3Ks) are relatively conserved and important intracellular lipid kinases involved in signalling and other biological pathways. In the free-living nematode Caenorhabditis elegans, the heterodimeric form of PI3K consists of catalytic (AGE-1) and regulatory (AAP-1) subunits. These subunits are key components of the insulin-like signalling pathway and play roles in the regulation of the entry into and exit from dauer. Although, in parasitic nematodes, similar components are proposed to regulate the transition from free-living or arrested stages to parasitic larvae, nothing is known about PI3Ks in relation to the transition of third-stage larvae (L3s) to parasitism in Haemonchus contortus., Methods: An integrated molecular approach was used to investigate age-1 and aap-1 of H. contortus (Hc-age-1 and Hc-aap-1) in C. elegans., Results: The two genes Hc-age-1 and Hc-aap-1 were transcribed in all life stages, with the highest levels in the egg, infective L3 and adult female of H. contortus. The expression of these genes was localized to the intestine, contrasting the pattern of their orthologues in C. elegans (where they are expressed in both head neurons and the intestine). The yeast two-hybrid analysis demonstrated that the adaptor-binding domain of Hc-AGE-1 interacted strongly with the Hc-AAP-1; however, this complex did not rescue the function of its orthologue in age-1-deficient C. elegans., Conclusions: This is the first time that the PI3K-encoding genes have been characterized from a strongylid parasitic nematode. The findings provide insights into the role of the PI3K heterodimer represented by Hc-age-1 and Hc-aap-1 in the developmental biology of H. contortus.

Background: Schistosomiasis is a parasitic disease affecting ~200 million people worldwide. Schistosoma haematobium and S. mansoni are two relatively closely related schistosomes (blood flukes), and the causative agents of urogenital and hepatointestinal schistosomiasis, respectively. The availability of genomic, transcriptomic and proteomic data sets for these two schistosomes now provides unprecedented opportunities to explore their biology, host interactions and schistosomiasis at the molecular level. A particularly important group of molecules involved in a range of biological and developmental processes in schistosomes and other parasites are the G protein-coupled receptors (GPCRs). Although GPCRs have been studied in schistosomes, there has been no detailed comparison of these receptors between closely related species. Here, using a genomic-bioinformatic approach, we identified and characterised key GPCRs in S. haematobium and S. mansoni (two closely related species of schistosome)., Methods: Using a Hidden Markov Model (HMM) and Support Vector Machine (SVM)-based pipeline, we classified and sub-classified GPCRs of S. haematobium and S. mansoni, combined with phylogenetic and transcription analyses., Results: We identified and classified classes A, B, C and F as well as an unclassified group of GPCRs encoded in the genomes of S. haematobium and S. mansoni. In addition, we characterised ligand-specific subclasses (i.e. amine, peptide, opsin and orphan) within class A (rhodopsin-like)., Conclusions: Most GPCRs shared a high degree of similarity and conservation, except for members of a particular clade (designated SmGPR), which appear to have diverged between S. haematobium and S. mansoni and might explain, to some extent, some of the underlying biological differences between these two schistosomes. The present set of annotated GPCRs provides a basis for future functional genomic studies of cellular GPCR-mediated signal transduction and a resource for future drug discovery efforts in schistosomes.