Search

Background The subfamily Phascolostrongylinae (Superfamily Strongyloidea) comprises nematodes that are parasitic in the gastrointestinal tracts of macropodid (Family Macropodidae) and vombatid (Family Vombatidae) marsupials. Currently, nine genera and 20 species have been attributed to the subfamily Phascolostrongylinae. Previous studies using sequence data sets for the internal transcribed spacers (ITS) of nuclear ribosomal DNA showed conflicting topologies between the Phascolostrongylinae and related subfamilies. Therefore, the aim of this study was to validate the phylogenetic relationships within the Phascolostrongylinae and its relationship with the families Chabertiidae and Strongylidae using mitochondrial amino acid sequences. Methods The sequences of all 12 mitochondrial protein-coding genes were obtained by next-generation sequencing of individual adult nematodes (n = 8) representing members of the Phascolostrongylinae. These sequences were conceptually translated and the phylogenetic relationships within the Phascolostrongylinae and its relationship with the families Chabertiidae and Strongylidae were inferred from aligned, concatenated amino acid sequence data sets. Results Within the Phascolostrongylinae, the wombat-specific genera grouped separately from the genera occurring in macropods. Two of the phascolostrongyline tribes were monophyletic, including Phascolostrongylinea and Hypodontinea, whereas the tribe Macropostrongyloidinea was paraphyletic. The tribe Phascolostrongylinea occurring in wombats was closely related to Oesophagostomum spp., also from the family Chabertiidae, which formed a sister relationship with the Phascolostrongylinae. Conclusion The current phylogenetic relationship within the subfamily Phascolostrongylinae supports findings from a previous study based on ITS sequence data. This study contributes also to the understanding of the phylogenetic position of the subfamily Phascolostrongylinae within the Chabertiidae. Future studies investigating the relationships between the Phascolostrongylinae and Cloacininae from macropodid marsupials may advance our knowledge of the phylogeny of strongyloid nematodes in marsupials. Graphical Abstract

Increasing access to next‐generation sequencing (NGS) technologies is revolutionizing the life sciences. In disease ecology, NGS‐based methods have the potential to provide higher‐resolution data on communities of parasites found in individual hosts as well as host populations. Here, we demonstrate how a novel analytical method, utilizing high‐throughput sequencing of PCR amplicons, can be used to explore variation in blood‐borne parasite (Theileria—Apicomplexa: Piroplasmida) communities of African buffalo at higher resolutions than has been obtained with conventional molecular tools. Results reveal temporal patterns of synchronized and opposite fluctuations of prevalence and relative abundance of Theileria spp. within the host population, suggesting heterogeneous transmission across taxa. Furthermore, we show that the community composition of Theileria spp. and their subtypes varies considerably between buffalo, with differences in composition reflected in mean and variance of overall parasitemia, thereby showing potential to elucidate previously unexplained contrasts in infection outcomes for host individuals. Importantly, our methods are generalizable as they can be utilized to describe blood‐borne parasite communities in any host species. Furthermore, our methodological framework can be adapted to any parasite system given the appropriate genetic marker. The findings of this study demonstrate how a novel NGS‐based analytical approach can provide fine‐scale, quantitative data, unlocking opportunities for discovery in disease ecology.

Background The parasitic flatworm Clonorchis sinensis inhabits the biliary tree of humans and other piscivorous mammals. This parasite can survive and thrive in the bile duct, despite exposure to bile constituents and host immune attack. Although the precise biological mechanisms underlying this adaptation are unknown, previous work indicated that Niemann-pick type C2 (NPC2)-like sterol-binding proteins might be integral in the host-parasite interplay. Expansions of this family in some invertebrates, such as arthropods, have shown functional diversification, including novel forms of chemoreception. Thus, here we curated the NPC2-like protein gene complement in C. sinensis, and predicted their conserved and/or divergent functional roles. Methods We used an established comparative genomic-bioinformatic approach to curate NPC2-like proteins encoded in published genomes of Korean and Chinese isolates of C. sinensis. Protein sequence and structural homology, presence of conserved domains and phylogeny were used to group and functionally classify NPC2-like proteins. Furthermore, transcription levels of NPC2-like protein-encoding genes were explored in different developmental stages and tissues. Results Totals of 35 and 32 C. sinensis NPC2-like proteins were predicted to be encoded in the genomes of the Korean and Chinese isolates, respectively. Overall, these proteins had low sequence homology and high variability of sequence alignment coverage when compared with curated NPC2s. Most C. sinensis proteins were predicted to retain a conserved ML domain and a conserved fold conformation, with a large cavity within the protein. Only one protein sequence retained the conserved amino acid residues required in bovine NPC2 to bind cholesterol. Non-canonical C. sinensis NPC2-like protein-coding domains clustered into four distinct phylogenetic groups with members of a group frequently encoded on the same genome scaffolds. Interestingly, NPC2-like protein-encoding genes were predicted to be variably transcribed in different developmental stages and adult tissues, with most being transcribed in the metacercarial stage. Conclusions The results of the present investigation confirms an expansion of NPC2-like proteins in C. sinensis, suggesting a diverse array of functions beyond sterol binding and transport. Functional explorations of this protein family should elucidate the mechanisms enabling the establishment and survival of C. sinensis and related flukes in the biliary systems of mammalian hosts.

Nematodes of the genus Macropostrongyloides inhabit the large intestines or stomachs of macropodid (kangaroos and wallabies) and vombatid (wombats) marsupials. This study established the relationships of seven species of Macropostrongyloides using mitochondrial (mt) protein amino acid sequence data sets. Phylogenetic analyses revealed that species of Macropostrongyloides (M. lasiorhini, M. baylisi, M. yamagutii, M. spearei, M. mawsonae and M. woodi) from the large intestines of their hosts formed a monophyletic assemblage with strong nodal support to the exclusion of M. dissimilis from the stomach of the swamp wallaby. Furthermore, the mitochondrial protein-coding genes provided greater insights into the diversity and phylogeny of the genus Macropostrongyloides, such data sets could potentially be used to elucidate the relationships among other parasitic nematodes of Australian marsupials.

Trichobilharzia species are parasitic flatworms (called schistosomes or flukes) that cause important diseases in birds and humans, but very little is known about their molecular biology. Here, using a transcriptomics-bioinformatics-based approach, we explored molecular aspects pertaining to the nutritional requirements of Trichobilharzia szidati (‘visceral fluke’) and T. regenti (‘neurotropic fluke’) in their avian host. We studied the larvae of each species before they enter (cercariae) and as they migrate (schistosomules) through distinct tissues in their avian (duck) host. Cercariae of both species were enriched for pathways or molecules associated predominantly with carbohydrate metabolism, oxidative phosphorylation and translation of proteins linked to ribosome biogenesis, exosome production and/or lipid biogenesis. Schistosomules of both species were enriched for pathways or molecules associated with processes including signal transduction, cell turnover and motility, DNA replication and repair, molecular transport and/or catabolism. Comparative informatic analyses identified molecular repertoires (within, e.g., peptidases and secretory proteins) in schistosomules that can broadly degrade macromolecules in both T. szidati and T. regenti, and others that are tailored to each species to selectively acquire nutrients from particular tissues through which it migrates. Thus, this study provides molecular evidence for distinct modes of nutrient acquisition between the visceral and neurotropic flukes of birds.

In this study, we explored the molecular alterations in the developmental switch from the L3 to the exsheathed L3 (xL3) and to the L4 stage of Haemonchus contortus in vitro using an integrated transcriptomic, proteomic and bioinformatic approach. Totals of 9,754 mRNAs, 88 microRNAs (miRNAs) and 1,591 proteins were identified, and 6,686 miRNA-mRNA pairs inferred in all larval stages studied. Approximately 16% of transcripts in the combined transcriptome (representing all three larval stages) were expressed as proteins, and there were positive correlations (r = 0.39-0.44) between mRNA transcription and protein expression in the three distinct developmental stages of the parasite. Of the predicted targets, 1,019 (27.0%) mRNA transcripts were expressed as proteins, and there was a negative correlation (r = -0.60 to -0.50) in the differential mRNA transcription and protein expression between developmental stages upon pairwise comparison. The changes in transcription (mRNA and miRNA) and protein expression from the free-living to the parasitic life cycle phase of H. contortus related to enrichments in biological pathways associated with metabolism (e.g., carbohydrate and lipid degradation, and amino acid metabolism), environmental information processing (e.g., signal transduction, signalling molecules and interactions) and/or genetic information processing (e.g., transcription and translation). Specifically, fatty acid degradation, steroid hormone biosynthesis and the Rap1 signalling pathway were suppressed, whereas transcription, translation and protein processing in the endoplasmic reticulum were upregulated during the transition from the free-living L3 to the parasitic xL3 and L4 stages of the nematode in vitro. Dominant post-transcriptional regulation was inferred to elicit these changes, and particular miRNAs (e.g., hco-miR-34 and hco-miR-252) appear to play roles in stress responses and/or environmental adaptations during developmental transitions of H. contortus. Taken together, these integrated results provide a comprehensive insight into the developmental biology of this important parasite at the molecular level in vitro. The approach applied here to H. contortus can be readily applied to other parasitic nematodes.

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans., Author summary Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

Increasing access to next‐generation sequencing (NGS) technologies is revolutionizing the life sciences. In disease ecology, NGS‐based methods have the potential to provide higher‐resolution data on communities of parasites found in individual hosts as well as host populations.Here, we demonstrate how a novel analytical method, utilizing high‐throughput sequencing of PCR amplicons, can be used to explore variation in blood‐borne parasite (Theileria—Apicomplexa: Piroplasmida) communities of African buffalo at higher resolutions than has been obtained with conventional molecular tools.Results reveal temporal patterns of synchronized and opposite fluctuations of prevalence and relative abundance of Theileria spp. within the host population, suggesting heterogeneous transmission across taxa. Furthermore, we show that the community composition of Theileria spp. and their subtypes varies considerably between buffalo, with differences in composition reflected in mean and variance of overall parasitemia, thereby showing potential to elucidate previously unexplained contrasts in infection outcomes for host individuals.Importantly, our methods are generalizable as they can be utilized to describe blood‐borne parasite communities in any host species. Furthermore, our methodological framework can be adapted to any parasite system given the appropriate genetic marker.The findings of this study demonstrate how a novel NGS‐based analytical approach can provide fine‐scale, quantitative data, unlocking opportunities for discovery in disease ecology., We used high‐throughput amplicon sequencing to uncover cryptic variation in blood‐borne parasite communities of African buffalo. Our methods are generalizable and will unlock opportunity for discovery in disease ecology.

Background Here, we created an automated pipeline for the de novoassembly of genomes from Pacific Biosciences long-read and Illumina short-read data using common workflow language (CWL). To evaluate the performance of this pipeline, we assembled the nuclear genomes of the eukaryotes Caenorhabditis elegans (∼100 Mb), Drosophila melanogaster (∼138 Mb), and Plasmodium falciparum (∼23 Mb) directly from publicly accessible nucleotide sequence datasets and assessed the quality of the assemblies against curated reference genomes. Findings We showed a dependency of the accuracy of assembly on sequencing technology and GC content and repeatedly achieved assemblies that meet the high standards set by the National Human Genome Research Institute, being applicable to gene prediction and subsequent genomic analyses. Conclusions This CWL pipeline overcomes current challenges of achieving repeatability and reproducibility of assembly results and offers a platform for the re-use of the workflow and the integration of diverse datasets. This workflow is publicly available via GitHub (https://github.com/vetscience/Assemblosis) and is currently applicable to the assembly of haploid and diploid genomes of eukaryotes.

Trichobilharzia species are parasitic flatworms (called schistosomes or flukes) that cause important diseases in birds and humans, but very little is known about their molecular biology. Here, using a transcriptomics-bioinformatics-based approach, we explored molecular aspects pertaining to the nutritional requirements of Trichobilharzia szidati (‘visceral fluke’) and T. regenti (‘neurotropic fluke’) in their avian host. We studied the larvae of each species before they enter (cercariae) and as they migrate (schistosomules) through distinct tissues in their avian (duck) host. Cercariae of both species were enriched for pathways or molecules associated predominantly with carbohydrate metabolism, oxidative phosphorylation and translation of proteins linked to ribosome biogenesis, exosome production and/or lipid biogenesis. Schistosomules of both species were enriched for pathways or molecules associated with processes including signal transduction, cell turnover and motility, DNA replication and repair, molecular transport and/or catabolism. Comparative informatic analyses identified molecular repertoires (within, e.g., peptidases and secretory proteins) in schistosomules that can broadly degrade macromolecules in both T. szidati and T. regenti, and others that are tailored to each species to selectively acquire nutrients from particular tissues through which it migrates. Thus, this study provides molecular evidence for distinct modes of nutrient acquisition between the visceral and neurotropic flukes of birds.

The prevalent porcine helminth, Ascaris suum, compromises pig health and reduces farm productivity worldwide. The closely related human parasite, A. lumbricoides, infects more than 800 million people representing a disease burden of 1.31 million disability-adjusted life years. The infections are often chronic in nature, and the parasites have a profound ability to modulate their hosts’ immune responses. This study provides the first in-depth characterisation of extracellular vesicles (EVs) from different developmental stages and body parts of A. suum and proposes the role of these vesicles in the host–parasite interplay. The release of EVs from the third- (L3) and fourth-stage (L4) larvae and adults was demonstrated by transmission electron microscopy (TEM), and sequencing of EV-derived RNA identified a number of microRNAs (miRNAs) and transcripts of potential host immune targets, such as IL-13, IL-25 and IL-33, were identified. Furthermore, proteomics of EVs identified several proteins with immunomodulatory properties and other proteins previously shown to be associated with parasite EVs. Taken together, these results suggest that A. suum EVs and their cargo may play a role in host–parasite interactions. This knowledge may pave the way to novel strategies for helminth infection control and knowledge of their immune modulatory potential.

Background Blood flukes of the genus Schistosoma cause schistosomiasis—a neglected tropical disease (NTD) that affects more than 200 million people worldwide. Studies of schistosome genomes have improved our understanding of the molecular biology of flatworms, but most of them have focused largely on protein-coding genes. Small non-coding RNAs (sncRNAs) have been explored in selected schistosome species and are suggested to play essential roles in the post-transcriptional regulation of genes, and in modulating flatworm-host interactions. However, genome-wide small RNA data are currently lacking for key schistosomes including Schistosoma haematobium—the causative agent of urogenital schistosomiasis of humans. Methodology MicroRNAs (miRNAs) and other sncRNAs of male and female adults of S. haematobium and small RNA transcription levels were explored by deep sequencing, genome mapping and detailed bioinformatic analyses. Principal findings In total, 89 transcribed miRNAs were identified in S. haematobium—a similar complement to those reported for the congeners S. mansoni and S. japonicum. Of these miRNAs, 34 were novel, with no homologs in other schistosomes. Most miRNAs (n = 64) exhibited sex-biased transcription, suggestive of roles in sexual differentiation, pairing of adult worms and reproductive processes. Of the sncRNAs that were not miRNAs, some related to the spliceosome (n = 21), biogenesis of other RNAs (n = 3) or ribozyme functions (n = 16), whereas most others (n = 3798) were novel (‘orphans’) with unknown functions. Conclusions This study provides the first genome-wide sncRNA resource for S. haematobium, extending earlier studies of schistosomes. The present work should facilitate the future curation and experimental validation of sncRNA functions in schistosomes to enhance our understanding of post-transcriptional gene regulation and of the roles that sncRNAs play in schistosome reproduction, development and parasite-host cross-talk., Author summary Human schistosomiasis is a chronic, neglected tropical disease (NTD) that is predominantly caused by the blood flukes Schistosoma haematobium, S. mansoni and S. japonicum. Infections by S. haematobium and/or S. mansoni are highly prevalent in Africa, affecting ~ 200 million people. The decoding of schistosome draft genomes has, to some extent, improved our understanding of the molecular biology of these parasites and now allows for non-protein-coding regions in these genomes to be characterised. Here, we explored small RNAs in adult S. haematobium by deep sequencing, reference genome mapping and detailed bioinformatic analyses. This study provides the first genome-wide miRNA and sncRNA resource for S. haematobium, extending earlier work on schistosomes and facilitating future curation efforts and functional investigations of schistosome sncRNAs. These efforts should enable a better understanding of post-transcriptional RNA modifications, gene regulation and novel aspects of parasite development, parasite-host cross-talk and disease at the molecular level.

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. Electronic supplementary material The online version of this article (10.1186/s13071-017-2462-4) contains supplementary material, which is available to authorized users.

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.

Compounded by a massive global food shortage, many parasitic diseases have a devastating, long-term impact on animal and human health and welfare worldwide. Parasitic helminths (worms) affect the health of billions of animals. Unlocking the systems biology of these neglected pathogens will underpin the design of new and improved interventions against them. Currently, the functional annotation of genomic and transcriptomic sequence data for socio-economically important parasitic worms relies almost exclusively on comparative bioinformatic analyses using model organism- and other databases. However, many genes and gene products of parasitic helminths (often > 50%) cannot be annotated using this approach, because they are specific to parasites and/or do not have identifiable homologs in other organisms for which sequence data are available. This inability to fully annotate transcriptomes and predicted proteomes is a major challenge and constrains our understanding of the biology of parasites, interactions with their hosts and of parasitism and the pathogenesis of disease on a molecular level. In the present article, we compiled transcriptomic data sets of key, socioeconomically important parasitic helminths, and constructed and validated a curated database, called HelmDB (www.helmdb.org). We demonstrate how this database can be used effectively for the improvement of functional annotation by employing data integration and clustering. Importantly, HelmDB provides a practical and user-friendly toolkit for sequence browsing and comparative analyses among divergent helminth groups (including nematodes and trematodes), and should be readily adaptable and applicable to a wide range of other organisms. This web-based, integrative database should assist ‘systems biology’ studies of parasitic helminths, and the discovery and prioritization of novel drug and vaccine targets. This focus provides a pathway toward developing new and improved approaches for the treatment and control of parasitic diseases, with the potential for important biotechnological outcomes.

Angiostrongylus vasorum is a metastrongyloid nematode of dogs and other canids of major clinical importance in many countries. In order to gain first insights into the molecular biology of this worm, we conducted the first large-scale exploration of its transcriptome, and predicted essential molecules linked to metabolic and biological processes as well as host immune responses. We also predicted and prioritized drug targets and drug candidates. Following Illumina sequencing (RNA-seq), 52.3 million sequence reads representing adult A. vasorum were assembled and annotated. The assembly yielded 20,033 contigs, which encoded proteins with 11,505 homologues in Caenorhabditis elegans, and additional 2252 homologues in various other parasitic helminths for which curated data sets were publicly available. Functional annotation was achieved for 11,752 (58.6%) proteins predicted for A. vasorum, including peptidases (4.5%) and peptidase inhibitors (1.6%), protein kinases (1.7%), G protein-coupled receptors (GPCRs) (1.5%) and phosphatases (1.2%). Contigs encoding excretory/secretory and immuno-modulatory proteins represented some of the most highly transcribed molecules, and encoded enzymes that digest haemoglobin were conserved between A. vasorum and other blood-feeding nematodes. Using an essentiality-based approach, drug targets, including neurotransmitter receptors, an important chemosensory ion channel and cysteine proteinase-3 were predicted in A. vasorum, as were associated small molecular inhibitors/activators. Future transcriptomic analyses of all developmental stages of A. vasorum should facilitate deep explorations of the molecular biology of this important parasitic nematode and support the sequencing of its genome. These advances will provide a foundation for exploring immuno-molecular aspects of angiostrongylosis and have the potential to underpin the discovery of new methods of intervention.

The apicoplast (ap) is a unique, non-photosynthetic organelle found in most apicomplexan parasites. Due to the essential roles that this organelle has, it has been widely considered as target for drugs against diseases caused by apicomplexans. Exploring the ap genomes of such parasites would provide a better understanding of their systematics and their basic molecular biology for therapeutics. However, there is limited information available on the ap genomes of apicomplexan parasites. In the present study, the ap genomes of two operational taxonomic units of Babesia (known as Babesia sp. Lintan [Bl] and Babesia sp. Xinjiang [Bx]) from sheep were sequenced, assembled and annotated using a massive parallel sequencing-based approach. Then, the gene content and gene order in these ap genomes (∼30.7kb in size) were defined and compared, and the genetic differences were assessed. In addition, a phylogenetic analysis of ap genomic data sets was carried out to assess the relationships of these taxonomic units with other apicomplexan parasites for which complete ap genomic data sets were publicly available. The results showed that the ap genomes of Bl and Bx encode 59 and 57 genes, respectively, including 2 ribosomal RNA genes, 25 transfer RNA genes and 30-32 protein-encoding genes, being similar in content to those of Babesia bovis and B. orientalis. Ap gene regions that might serve as markers for future epidemiological and population genetic studies of Babesia species were identified. Using sequence data for a subset of six protein-encoding genes, a close relationship of Bl and Bx with Babesia bovis from cattle and B. orientalis from water buffalo was inferred. Although the focus of the present study was on Babesia, we propose that the present sequencing-bioinformatic approach should be applicable to organellar genomes of a wide range of apicomplexans of veterinary importance.

Here, we sequenced, assembled and annotated the mitochondrial (mt) genomes of two operational taxonomic units of Babesia from sheep from China using a deep sequencing-coupled approach. Then, we defined and compared the gene order of these mt genomes (~5.8 to 6.2kb in size), assessed sequence differences in mt genes among Babesia taxa and evaluated genetic relationships among these taxa and related apicomplexans (Theileria) for which mt genomic data sets were available. We also identified mt genetic regions that might be useful as markers for future population genetic and molecular epidemiological studies of Babesia from small ruminants. We propose that the sequencing-bioinformatic approach used here should be applicable to a wide range of protists of veterinary importance.

The advent and integration of high-throughput '-omics' technologies (e.g. genomics, transcriptomics, proteomics, metabolomics, glycomics and lipidomics) are revolutionizing the way biology is done, allowing the systems biology of organisms to be explored. These technologies are now providing unique opportunities for global, molecular investigations of parasites. For example, studies of a transcriptome (all transcripts in an organism, tissue or cell) have become instrumental in providing insights into aspects of gene expression, regulation and function in a parasite, which is a major step to understanding its biology. The purpose of this article was to review recent applications of next-generation sequencing technologies and bioinformatic tools to large-scale investigations of the transcriptomes of parasitic nematodes of socio-economic significance (particularly key species of the order Strongylida) and to indicate the prospects and implications of these explorations for developing novel methods of parasite intervention.

Liver flukes of animals are parasitic flatworms (Platyhelminthes: Digenea) of major socioeconomic importance in many countries. Key representatives, such as Fasciola hepatica and F. gigantica, cause "liver fluke disease" (= fascioliasis), which is of major animal health significance worldwide. In particular, F. hepatica is a leading cause of production losses to the livestock (mainly sheep and cattle) and meat industries due to clinical disease, reduced weight gain and milk production, and deaths. This parasite is also a major food-borne pathogen of humans throughout parts of the Middle East, Asia and South America. Currently, there is a significant focus on the development of new approaches for the prevention and control of fascioliasis in livestock. Recent technological advances in genomics and bioinformatics provide unique opportunities for the identification and prevalidation of drug targets and vaccines through a better understanding of the biology of F. hepatica and related species as well as their relationship with their hosts at the molecular level. Surprisingly, despite the widespread socioeconomic impact of fascioliasis, genomic datasets for F. hepatica are scant, limiting the molecular biological research of this parasite. The present article explores specifically the transcriptome of the adult stage of F. hepatica using an integrated genomic-bioinformatic platform. The analysis of the current data reveals numerous molecules of biological relevance, some of which are inferred to be involved in key biological processes or pathways that could serve as targets for new trematocidal drugs or vaccines. Improved insights into the transcriptome of F. hepatica should pave the way for future, comparative analysis of the transcriptomes of other developmental stages of this and related parasites, such as F. gigantica, cancer-causing flatworms (Clonorchis sinensis and Opisthorchis viverrini) and blood flukes (Schistosoma mansoni and S. japonicum). Prediction of the essentiality of genes and their products, molecular network connectivity of trematode genes as well as experimental exploration of function should also add value to the genomic discovery efforts in the future, focused on biotechnological outcomes.

Mitochondrial (mt) genomics represents an understudied but important field of molecular biology. Increasingly, mt dysfunction is being linked to a range of human diseases, including neurodegenerative disorders, diabetes and impairment of childhood development. In addition, mt genomes provide important markers for systematic, evolutionary and population genetic studies. Some technological limitations have prevented the expanded generation and utilization of mt genomic data for some groups of organisms. These obstacles most acutely impede, but are not limited to, studies requiring the determination of complete mt genomic data from minute amounts of material (e.g. biopsy samples or microscopic organisms). Furthermore, post-sequencing bioinformatic annotation and analyses of mt genomes are time consuming and inefficient. Herein, we describe a high-throughput sequencing and bioinformatic pipeline for mt genomics, which will have implications for the annotation and analysis of other organellar (e.g. plastid or apicoplast genomes) and virus genomes as well as long, contiguous regions in nuclear genomes. We utilize this pipeline to sequence and annotate the complete mt genomes of 12 species of parasitic nematode (order Strongylida) simultaneously, each from an individual organism. These mt genomic data provide a rich source of markers for studies of the systematics and population genetics of a group of socioeconomically important pathogens of humans and other animals.

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. Electronic supplementary material The online version of this article (doi:10.1186/s13071-016-1428-2) contains supplementary material, which is available to authorized users.

Parasitic worm proteins that belong to the cysteine-rich secretory proteins, antigen 5 and pathogenesis-related 1 (CAP) superfamily are proposed to play key roles in the infection process and the modulation of immune responses in host animals. However, there is limited information on these proteins for most socio-economically important worms. Here, we review the CAP protein superfamily of Haemonchus contortus (barber's pole worm), a highly significant parasitic roundworm (order Strongylida) of small ruminants. To do this, we mined genome and transcriptomic datasets, predicted and curated full-length amino acid sequences (n = 45), undertook systematic phylogenetic analyses of these data and investigated transcription throughout the life cycle of H. contortus. We inferred functions for selected Caenorhabditis elegans orthologs (including vap-1, vap-2, scl-5 and lon-1) based on genetic networking and by integrating data and published information, and were able to infer that a subset of orthologs and their interaction partners play pivotal roles in growth and development via the insulin-like and/or the TGF-beta signalling pathways. The identification of the important and conserved growth regulator LON-1 led us to appraise the three-dimensional structure of this CAP protein by comparative modelling. This model revealed the presence of different topological moieties on the canonical fold of the CAP domain, which coincide with an overall charge separation as indicated by the electrostatic surface potential map. These observations suggest the existence of separate sites for effector binding and receptor interactions, and thus support the proposal that these worm molecules act in similar ways as venoms act as ligands for chemokine receptors or G protein-coupled receptor effectors. In conclusion, this review should guide future molecular studies of these molecules, and could support the development of novel interventions against haemonchosis.

The blow fly, Lucilia cuprina (Wiedemann, 1830) is a parasitic insect of major global economic importance. Maggots of this fly parasitize the skin of animal hosts, feed on excretions and tissues, and cause severe disease (flystrike or myiasis). Although there has been considerable research on L. cuprina over the years, little is understood about the molecular biology, biochemistry and genetics of this parasitic fly, as well as its relationship with its hosts and the disease that it causes. This situation might change with the recent report of the draft genome and transcriptome of this blow fly, which has given new and global insights into its biology, interactions with the host animal and aspects of insecticide resistance at the molecular level. This genomic resource will likely enable many fundamental and applied research areas in the future. The present article gives a background on L. cuprina and myiasis, a brief account of past and current treatment, prevention and control approaches, and provides a perspective on the impact that the L. cuprina genome should have on future research of this and related parasitic flies, and the design of new and improved interventions for myiasis.

Although substantial research has been focused on the 'hidden antigen' H11 of Haemonchus contortus as a vaccine against haemonchosis in small ruminants, little is know about this and related aminopeptidases. In the present article, we reviewed genomic and transcriptomic data sets to define, for the first time, the complement of aminopeptidases (designated Hc-AP-1 to Hc-AP-13) of the family M1 with homologues in Caenorhabditis elegans, characterised by zinc-binding (HEXXH) and exo-peptidase (GAMEN) motifs. The three previously published H11 isoforms (accession nos. X94187, FJ481146 and AJ249941) had most sequence similarity to Hc-AP-2 and Hc-AP-8, whereas unpublished isoforms (accession nos. AJ249942 and AJ311316) were both most similar to Hc-AP-3. The aminopeptidases characterised here had homologues in C. elegans. Hc-AP-1 to Hc-AP-8 were most similar in amino acid sequence (28-41%) to C. elegans T07F10.1; Hc-AP-9 and Hc-AP-10 to C. elegans PAM-1 (isoform b) (53-54% similar); Hc-AP-11 and Hc-AP-12 to C. elegans AC3.5 and Y67D8C.9 (26% and 50% similar, respectively); and Hc-AP-13 to C. elegans C42C1.11 and ZC416.6 (50-58% similar). Comparative analysis suggested that Hc-AP-1 to Hc-AP-8 play roles in digestion, metabolite excretion, neuropeptide processing and/or osmotic regulation, with Hc-AP-4 and Hc-AP-7 having male-specific functional roles. The analysis also indicated that Hc-AP-9 and Hc-AP-10 might be involved in the degradation of cyclin (B3) and required to complete meiosis. Hc-AP-11 represents a leucyl/cystinyl aminopeptidase, predicted to have metallopeptidase and zinc ion binding activity, whereas Hc-AP-12 likely encodes an aminopeptidase Q homologue also with these activities and a possible role in gonad function. Finally, Hc-AP-13 is predicted to encode an aminopeptidase AP-1 homologue of C. elegans with hydrolase activity, suggested to operate, possibly synergistically with a PEPT-1 ortholog, as an oligopeptide transporter in the gut for protein uptake and normal development and/or reproduction of the worm. An appraisal of structure-based amino acid sequence alignments revealed that all conceptually translated Hc-AP proteins, with the exception of Hc-AP-12, adopt a topology similar to those observed for the two subgroups of mammalian M1 aminopeptidases, which possess either three (I, II and IV) or four (I-IV) domains. In contrast, Hc-AP-12 lacks the N-terminal domain (I), but possesses a substantially expanded domain III. Although further work needs to be done to assess amino acid sequence conservation of the different aminopeptidases among individual worms within and among H. contortus populations, we hope that these insights will support future localisation, structural and functional studies of these molecules in H. contortus as well as facilitate future assessments of a recombinant subunit or cocktail vaccine against haemonchosis.

Lucilia cuprina is a parasitic fly of major economic importance worldwide. Larvae of this fly invade their animal host, feed on tissues and excretions and progressively cause severe skin disease (myiasis). Here we report the sequence and annotation of the 458-megabase draft genome of Lucilia cuprina. Analyses of this genome and the 14,544 predicted protein-encoding genes provide unique insights into the fly's molecular biology, interactions with the host animal and insecticide resistance. These insights have broad implications for designing new methods for the prevention and control of myiasis., Lucilia cuprina is a parasitic blowfly of major economic importance worldwide that feeds on the tissues of animals such as sheep. Here, the authors sequence the genome of L. cuprina and provide insights into the fly's molecular biology, interactions with the host animal and insecticide resistance.

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.

Giardia duodenalis is the most common gastrointestinal protozoan parasite of humans and a significant contributor to the global burden of both diarrheal disease and post-infectious chronic disorders. Although G. duodenalis can be cultured axenically, significant gaps exist in our understanding of the molecular biology and metabolism of this pathogen. The present study employed RNA sequencing to characterize the mRNA transcriptome of G. duodenalis trophozoites in axenic culture, at log (48 h of growth), stationary (60 h), and declining (96 h) growth phases. Using ~400-times coverage of the transcriptome, we identified 754 differentially transcribed genes (DTGs), mainly representing two large DTG groups: 438 that were down-regulated in the declining phase relative to log and stationary phases, and 281 that were up-regulated. Differential transcription of prominent antioxidant and glycolytic enzymes implicated oxygen tension as a key factor influencing the transcriptional program of axenic trophozoites. Systematic bioinformatic characterization of numerous DTGs encoding hypothetical proteins of unknown function was achieved using structural homology searching. This powerful approach greatly informed the differential transcription analysis and revealed putative novel antioxidant-coding genes, and the presence of a near-complete two-component-like signaling system that may link cytosolic redox or metabolite sensing to the observed transcriptional changes. Motif searching applied to promoter regions of the two large DTG groups identified different putative transcription factor-binding motifs that may underpin global transcriptional regulation. This study provides new insights into the drivers and potential mediators of transcriptional variation in axenic G. duodenalis and provides context for static transcriptional studies., Author Summary Giardia is the most common gastrointestinal protozoan parasite of humans. This parasite causes diarrheal disease and is correlated with post-infectious conditions such as irritable bowel syndrome. In the absence of a vaccine, treatment is limited to drugs such as metronidazole, against which clinical resistance is reported. Effective control of Giardia requires a detailed understanding of its biology, and in turn, complete characterization of the standard in vitro culture system. Using RNA sequencing assisted by informatics to functionally annotate hypothetical proteins, we investigated transcriptional changes in axenic Giardia trophozoites at three growth phases over 96 hours. We found two large groups of differentially transcribed genes that indicate changes in the antioxidant system and central carbon metabolism over time. A putative novel signaling pathway may act together with putative transcription factor-binding motifs to regulate these transcriptional changes. Our results suggest that dissolved oxygen in Giardia culture medium may cause oxidative stress early during in vitro growth and that oxygen depletion may limit the efficiency of glycolysis in the declining phase. This work enhances our understanding of the transcriptional flexibility and metabolism of Giardia in vitro.

Opisthorchiasis is a neglected, tropical disease caused by the carcinogenic Asian liver fluke, Opisthorchis viverrini. This hepatobiliary disease is linked to malignant cancer (cholangiocarcinoma, CCA) and affects millions of people in Asia. No vaccine is available, and only one drug (praziquantel) is used against the parasite. Little is known about O. viverrini biology and the diseases that it causes. Here we characterize the draft genome (634.5 Mb) and transcriptomes of O. viverrini, elucidate how this fluke survives in the hostile environment within the bile duct and show that metabolic pathways in the parasite are highly adapted to a lipid-rich diet from bile and/or cholangiocytes. We also provide additional evidence that O. viverrini and other flukes secrete proteins that directly modulate host cell proliferation. Our molecular resources now underpin profound explorations of opisthorchiasis/CCA and the design of new interventions., The Asian liver fluke is a parasitic worm that is linked to an increased risk of malignant cancer. Here, the authors sequence the draft genome and transcriptome of this fluke and provide insight into how the species has adapted to be able to survive in the bile duct.

In the present study we sequenced or re-sequenced, assembled and annotated 15 mitochondrial genomes representing the 12 currently recognised taxa of Trichinella using a deep sequencing-coupled approach. We then defined and compared the gene order in individual mitochondrial genomes (∼14 to 17.7 kb), evaluated genetic differences among species/genotypes and re-assessed the relationships among these taxa using the mitochondrial nucleic acid or amino acid sequence data sets. In addition, a rich source of mitochondrial genetic markers was defined that could be used in future systematic, epidemiological and population genetic studies of Trichinella. The sequencing-bioinformatic approach employed herein should be applicable to a wide range of eukaryotic parasites.

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.

Parasitic nematodes cause substantial morbidity and mortality in animals and people globally and major losses to food production annually. Ascaris is among the commonest geohelminths of swine and people worldwide, and causes major disease and socioeconomic losses, particularly in developing countries. The control of ascariasis has become a global health and welfare priority, but current treatment programs carry a significant risk of inducing anthelmintic resistance. Therefore, there is a need to work toward the sustainable control of Ascaris/ascariasis, built on a solid understanding of its molecular biology and genetics. Recently, we reported the 273 megabase (Mb) draft genome of Ascaris suum (sequenced from the reproductive tract of a single adult female worm) and explored transcription in different organs, stages, and both sexes of this nematode using advanced sequencing and computer technologies. We characterized key genes and biological pathways linked to the parasite's migration in the host, and its immunobiology, reproduction, and development. We also predicted and prioritized drug targets in A. suum, providing a basis for discovering new groups of nematocides. The present chapter provides an account of these recent advances, describes new methodologies established, and emphasizes prospects for profound investigations into the comparative genomics, genetics, evolution, immunobiology, epidemiology, and ecology of Ascaris from both pig and human hosts as well as for the development of new interventions against ascariasis and other helminthiases. © 2013 Copyright

Fasciola hepatica is a parasitic trematode that infects a wide range of mammalian hosts, including livestock and humans, in temperate and tropical regions globally. This trematode causes the disease fascioliasis, which consists of an acute phase (≤ 12 weeks) during which juvenile parasites migrate through the host liver tissues, and a chronic phase (> 12 weeks) following the establishment of adult parasites in the liver bile ducts. Few studies have explored the progression of the host response over the course of Fasciola infection in the same animals. In this study, we characterized transcriptomic changes in peripheral blood mononuclear cells (PBMCs) collected from sheep at three time points over the first eight weeks of infection relative to uninfected controls. In total, 183 and 76 genes were found to be differentially transcribed at two and eight weeks post-infection respectively. Functional and pathway analysis of differentially transcribed genes revealed changes related to T-cell activation that may underpin a Th2-biased immune response against this parasite. This first insight into the dynamics of host responses during the early stages of infection improves the understanding of the pathogenesis of acute fascioliasis, informs vaccine development and presents a set of PBMC markers with diagnostic potential.

Aelurostrongylus abstrusus (Railliet, 1898) is a metastrongylid nematode of major clinical relevance in felids, causing aelurostrongylosis. In spite of its clinical importance in cats, the genetics, epidemiology and biology of this parasite are not entirely understood. mt DNA can provide markers for studies of these areas, but genetic data are scant for A. abstrusus and related lungworms. Here, the mt genome was amplified by long-range polymerase chain reaction (long-PCR) from a single male adult of A. abstrusus, sequenced using 454 technology and annotated using an established bioinformatic pipeline. This circular mt genome is 13,913 bp and contains two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with most other chromadorean nematodes. This genome should provide a source of markers for future investigations of the epidemiology and ecology of A. abstrusus. Molecular tools, employing such mt markers, are likely to find utility for explorations into the epidemiology, biology and systematics of this parasite, and the diagnosis of feline aelurostrongylosis.

Background The Southeast Asian liver fluke (Opisthorchis viverrini) chronically infects and affects tens of millions of people in regions of Asia, leading to chronic illness and, importantly, inducing malignant cancer ( = cholangiocarcinoma). In spite of this, little is known, at the molecular level, about the parasite itself, its interplay with its hosts or the mechanisms of disease and/or carcinogenesis. Methodology/Principal Findings Here, we generated extensive RNA-Seq data (Illumina) representing adult and juvenile stages of O. viverrini, and combined these sequences with previously published transcriptomic data (454 technology) for this species, yielding a combined assembly of significantly increased quality and allowing quantitative assessment of transcription in the juvenile and adult stage. Conclusions This enhanced assembly reveals that, despite the substantial biological similarities between the human liver flukes, O. viverinni and Clonorchis sinensis, there are previously unrecognized differences in major aspects of their molecular biology. Most notable are differences among the C13 and cathepsin L-like cysteine peptidases, which play key roles in tissue migration, immune evasion and feeding, and, thus, represent potential drug and/or vaccine targets. Furthermore, these data indicate that major lineages of cysteine peptidases of socioeconomically important trematodes have evolved through a process of gene loss rather than independent radiation, contrasting previous proposals., Author Summary Opistorchis viverrini is an important and neglected parasite affecting ∼9 million people in South-east Asia. The parasite has a complex life-cycle which involves an intermediate phase in cyprinoid fishes. Consumption of raw or under-cooked fish infected with the metacercarial (larval) stage of O. viverreni results in infection, with adult worms living primarily in the intra-hepatic bile duct. In addition to the affects of the infection itself, O. viverrini is directly carcinogenic, with up to 70% of infected individuals in endemic regions developing malignant cholangiocarcinomas. Control of the parasite relies exclusively on the use of praziquantel and little is known about the mechanisms through which O. viverrini stimulates carcinogenesis. An improved understanding of the molecular biology of O. viverrini is urgently needed. In our study, we employed RNAseq technology to assess changes in gene transcription during the development of O. viverrini within the definitive host, and significantly improved the characterization of the transcriptome of this parasite. In so doing, we shed new light on the evolution of a major group proteins (i.e., the cysteine peptidases) which, given their important function roles as excreted/secreted molecules, have been proposed as attractive drug/vaccine targets for a wide-range of neglected flukes, including species of Opistorchis, Clonorchis, Schistosoma and Fasciola.

Parasitic diseases have a devastating, long-term impact on human health, welfare and food production worldwide. More than two billion people are infected with geohelminths, including the roundworms Ascaris (common roundworm), Necator and Ancylostoma (hookworms), and Trichuris (whipworm), mainly in developing or impoverished nations of Asia, Africa and Latin America. In humans, the diseases caused by these parasites result in about 135,000 deaths annually, with a global burden comparable with that of malaria or tuberculosis in disability-adjusted life years. Ascaris alone infects around 1.2 billion people and, in children, causes nutritional deficiency, impaired physical and cognitive development and, in severe cases, death. Ascaris also causes major production losses in pigs owing to reduced growth, failure to thrive and mortality. The Ascaris-swine model makes it possible to study the parasite, its relationship with the host, and ascariasis at the molecular level. To enable such molecular studies, we report the 273 megabase draft genome of Ascaris suum and compare it with other nematode genomes. This genome has low repeat content (4.4%) and encodes about 18,500 protein-coding genes. Notably, the A. suum secretome (about 750 molecules) is rich in peptidases linked to the penetration and degradation of host tissues, and an assemblage of molecules likely to modulate or evade host immune responses. This genome provides a comprehensive resource to the scientific community and underpins the development of new and urgently needed interventions (drugs, vaccines and diagnostic tests) against ascariasis and other nematodiases.

Fasciola gigantica (Digenea) is an important foodborne trematode that causes liver fluke disease (fascioliasis) in mammals, including ungulates and humans, mainly in tropical climatic zones of the world. Despite its socioeconomic impact, almost nothing is known about the molecular biology of this parasite, its interplay with its hosts, and the pathogenesis of fascioliasis. Modern genomic technologies now provide unique opportunities to rapidly tackle these exciting areas. The present study reports the first transcriptome representing the adult stage of F. gigantica (of bovid origin), defined using a massively parallel sequencing-coupled bioinformatic approach. From >20 million raw sequence reads, >30,000 contiguous sequences were assembled, of which most were novel. Relative levels of transcription were determined for individual molecules, which were also characterized (at the inferred amino acid level) based on homology, gene ontology, and/or pathway mapping. Comparisons of the transcriptome of F. gigantica with those of other trematodes, including F. hepatica, revealed similarities in transcription for molecules inferred to have key roles in parasite-host interactions. Overall, the present dataset should provide a solid foundation for future fundamental genomic, proteomic, and metabolomic explorations of F. gigantica, as well as a basis for applied outcomes such as the development of novel methods of intervention against this neglected parasite., Author Summary Fasciola gigantica (Digenea) is a socioeconomically important liver fluke of humans and other mammals. It is the predominant cause of fascioliasis in the tropics and has a serious impact on the lives of tens of millions of people and other animals; yet, very little is known about this parasite and its relationship with its hosts at the molecular level. Here, advanced sequencing and bioinformatic technologies were employed to explore the genes transcribed in the adult stage of F. gigantica. From >20 million raw reads, >30,000 contiguous sequences were assembled. Relative levels of transcription were estimated; and molecules were characterized based on homology, gene ontology, and/or pathway mapping. Comparisons of the transcriptome of F. gigantica with those of other trematodes, including F. hepatica, showed similarities in transcription for molecules predicted to play roles in parasite-host interactions. The findings of the present study provide a foundation for a wide range of fundamental molecular studies of this neglected parasite, as well as research focused on developing new methods for the treatment, diagnosis, and control of fascioliasis.