De novo assembly and characterization of the transcriptome of the pancreatic fluke Eurytrema pancreaticum (trematoda: Dicrocoeliidae) using Illumina paired-end sequencing

Label-free quantitative proteomics and immunoblotting identifies immunoreactive and other excretory-secretory (E/S) proteins of Anoplocephala perfoliata

Anoplocephala perfoliata is a common tapeworm in horses causing colic and even mortalities. Current diagnostic tests to detect A. perfoliata infections have their limitations and an improved method is needed. Immunoreactive excretory/secretory proteins (E/S proteome) of this parasite can provide promising candidates for diagnostic tests. We compared E/S proteins produced by small (length < 20 mm, width < 5 mm) and large (length 20 to 40 mm, width 5 to 10 mm) A. perfoliata worms in vitro by label-free quantitative proteomics using a database composed of related Hymenolepis diminuta, Echinococcus multilocularis/granulosus and Taenia aseatica proteins for protein identifications. Altogether, 509 E/S proteins were identified after incubating the worms in vitro for three and eight hours. The greatest E/S proteome changes suggested both worm size- and time-dependent changes in cytoskeleton remodeling, apoptosis, and production of antigens/immunogens. The E/S proteins collected at the three-hour time point represented the natural conditions better than those collected at the eight-hour time point, and thereby contained the most relevant diagnostic targets. Immunoblotting using antibodies from horses tested positive/negative for A. perfoliata indicated strongest antigenicity/immunogenicity with 13-, 30- and 100-kDa proteins, involving a thioredoxin, heat-shock chaperone 90 (Hsp90), dynein light chain component (DYNLL), tubulin-specific chaperone A (TBCA) and signaling pathway modulators (14-3-3 and Sj-Ts4). This is among the first studies identifying new diagnostic targets and A. perfoliata antigens eliciting a IgG-response in horses.

Pathological findings of bovine pancreatic lesions induced by Eurytrema pancreticum in Aceh cattle, Sumatra

Great losses of cattle and other ruminants due to Eurytrema pancreaticum have been recorded in Indonesia and other countries in Southeast Asia. The objective of this study was to examine histological alterations after E. pancreaticum infection in cattle. This study analysed the pathological changes of pancreas in 153 cattle sampled randomly at slaughterhouses in Banda Aceh, Indonesia. Samples were obtained during one year as part of routine meat inspection. The samples of cattle pancreas were obtained from slaughterhouse in Banda Aceh and examined for eurytrematosis. The number of infected pancreases was 64 out of 153 (41.8%). The pancreases indicated some alterations including macroscopic colour changes of the pancreatic capsular surface, followed by the production of mucus on the surface. Some adult flukes were found in the capsular surface of the pancreas together with fluke eggs. The dissected pancreases were prepared for histopathological study for each segment, and then observed under microscope. The implications of findings from histopathologic analyses of the pancreas are discussed.

Integrative Transcriptomics and Proteomics Analyses to Reveal the Developmental Regulation of Metorchis orientalis: A Neglected Trematode With Potential Carcinogenic Implications

Metorchis orientalis is a neglected zoonotic parasite of the gallbladder and bile duct of poultry, mammals, and humans. It has been widely reported in Asian, including China, Japanese, and Korea, where it is a potential threat to public health. Despite its significance as an animal and human pathogen, there are few published transcriptomic and proteomics data available. Transcriptome Illumina RNA sequencing and label-free protein quantification were performed to compare the gene and protein expression of adult and metacercariae-stage M. orientalis, resulting in 100,234 unigenes and 3,530 proteins. Of these, 13,823 differentially expressed genes and 1,445 differentially expressed proteins were identified in adult versus metacercariae. In total, 570 genes were differentially expressed consistent with the mRNA and protein level in the adult versus metacercariae stage. Differential gene transcription analyses revealed 34,228 genes to be expressed in both stages, whereas 66,006 genes showed stage-specific expression. Compared with adults, the metacercariae stage was highly transcriptional. GO and KEGG analyses based on transcriptome and proteome revealed numerous up-regulated genes in adult M. orientalis related to microtubule-based processes, microtubule motor activity, and nucleocytoplasmic transport. The up-regulated genes in metacercariae M. orientalis were mainly related to transmembrane receptor protein serine/threonine kinase activity, transmembrane receptor protein serine/threonine kinase signaling pathway. Transcriptome and proteome comparative analyses showed numerous up-regulated genes in adult stage were mainly enriched in actin filament capping, spectrin, and glucose metabolic process, while up-regulated genes in metacercariae stage were mainly related to cilium assembly, cilium movement, and motile cilium. These results highlight changes in protein and gene functions during the development of metacercariae into adults, and provided evidence for the mechanisms involved in morphological and metabolic changes at both the protein and gene levels. Interestingly, many genes had been proved associated with liver fibrosis and carcinogenic factors were identified highly expressed in adult M. orientalis, which suggests that M. orientalis is a neglected trematode with potential carcinogenic implications. These data provide attractive targets for the development of therapeutic or diagnostic interventions for controlling M. orientalis.

Evidence of Immune Modulators in the Secretome of the Equine Tapeworm Anoplocephala perfoliata

Anoplocephala perfoliata is a neglected gastro-intestinal tapeworm, commonly infecting horses worldwide. Molecular investigation of A. perfoliata is hampered by a lack of tools to better understand the host-parasite interface. This interface is likely influenced by parasite derived immune modulators released in the secretome as free proteins or components of extracellular vesicles (EVs). Therefore, adult RNA was sequenced and de novo assembled to generate the first A. perfoliata transcriptome. In addition, excretory secretory products (ESP) from adult A. perfoliata were collected and EVs isolated using size exclusion chromatography, prior to proteomic analysis of the EVs, the EV surface and EV depleted ESP. Transcriptome analysis revealed 454 sequences homologous to known helminth immune modulators including two novel Sigma class GSTs, five α-HSP90s, and three α-enolases with isoforms of all three observed within the proteomic analysis of the secretome. Furthermore, secretome proteomics identified common helminth proteins across each sample with known EV markers, such as annexins and tetraspanins, observed in EV fractions. Importantly, 49 of the 454 putative immune modulators were identified across the secretome proteomics contained within and on the surface of EVs in addition to those identified in free ESP. This work provides the molecular tools for A. perfoliata to reveal key players in the host-parasite interaction within the horse host.

Molecular signatures of the rediae, cercariae and adult stages in the complex life cycles of parasitic flatworms (Digenea: Psilostomatidae)

Background

Parasitic flatworms (Trematoda: Digenea) represent one of the most remarkable examples of drastic morphological diversity among the stages within a life cycle. Which genes are responsible for extreme differences in anatomy, physiology, behavior, and ecology among the stages? Here we report a comparative transcriptomic analysis of parthenogenetic and amphimictic generations in two evolutionary informative species of Digenea belonging to the family Psilostomatidae.

Methods

In this study the transcriptomes of rediae, cercariae and adult worm stages of Psilotrema simillimum and Sphaeridiotrema pseudoglobulus, were sequenced and analyzed. High-quality transcriptomes were generated, and the reference sets of protein-coding genes were used for differential expression analysis in order to identify stage-specific genes. Comparative analysis of gene sets, their expression dynamics and Gene Ontology enrichment analysis were performed for three life stages within each species and between the two species.

Results

Reference transcriptomes for P. simillimum and S. pseudoglobulus include 21,433 and 46,424 sequences, respectively. Among 14,051 orthologous groups (OGs), 1354 are common and specific for two analyzed psilostomatid species, whereas 13 and 43 OGs were unique for P. simillimum and S. pseudoglobulus, respectively. In contrast to P. simillimum, where more than 60% of analyzed genes were active in the redia, cercaria and adult worm stages, in S. pseudoglobulus less than 40% of genes had such a ubiquitous expression pattern. In general, 7805 (36.41%) and 30,622 (65.96%) of genes were preferentially expressed in one of the analyzed stages of P. simillimum and S. pseudoglobulus, respectively. In both species 12 clusters of co-expressed genes were identified, and more than a half of the genes belonging to the reference sets were included into these clusters. Functional specialization of the life cycle stages was clearly supported by Gene Ontology enrichment analysis.

Conclusions

During the life cycles of the two species studied, most of the genes change their expression levels considerably, consequently the molecular signature of a stage is not only a unique set of expressed genes, but also the specific levels of their expression. Our results indicate unexpectedly high level of plasticity in gene regulation between closely related species. Transcriptomes of P. simillimum and S. pseudoglobulus provide high quality reference resource for future evolutionary studies and comparative analyses.

Complex and dynamic transcriptional changes allow the helminth Fasciola gigantica to adjust to its intermediate snail and definitive mammalian hosts

BACKGROUND

The tropical liver fluke, Fasciola gigantica causes fasciolosis, an important disease of humans and livestock. We characterized dynamic transcriptional changes associated with the development of the parasite in its two hosts, the snail intermediate host and the mammalian definitive host.

RESULTS

Differential gene transcription analysis revealed 7445 unigenes transcribed by all F. gigantica lifecycle stages, while the majority (n = 50,977) exhibited stage-specific expression. Miracidia that hatch from eggs are highly transcriptionally active, expressing a myriad of genes involved in pheromone activity and metallopeptidase activity, consistent with snail host finding and invasion. Clonal expansion of rediae within the snail correlates with increased expression of genes associated with transcription, translation and repair. All intra-snail stages (miracidia, rediae and cercariae) require abundant cathepsin L peptidases for migration and feeding and, as indicated by their annotation, express genes putatively involved in the manipulation of snail innate immune responses. Cercariae emerge from the snail, settle on vegetation and become encysted metacercariae that are infectious to mammals; these remain metabolically active, transcribing genes involved in regulation of metabolism, synthesis of nucleotides, pH and endopeptidase activity to assure their longevity and survival on pasture. Dramatic growth and development following infection of the mammalian host are associated with high gene transcription of cell motility pathways, and transport and catabolism pathways. The intra-mammalian stages temporally regulate key families of genes including the cathepsin L and B proteases and their trans-activating peptidases, the legumains, during intense feeding and migration through the intestine, liver and bile ducts. While 70% of the F. gigantica transcripts share homology with genes expressed by the temperate liver fluke Fasciola hepatica, gene expression profiles of the most abundantly expressed transcripts within the comparable lifecycle stages implies significant species-specific gene regulation.

CONCLUSIONS

Transcriptional profiling of the F. gigantica lifecycle identified key metabolic, growth and developmental processes the parasite undergoes as it encounters vastly different environments within two very different hosts. Comparative analysis with F. hepatica provides insight into the similarities and differences of these parasites that diverged > 20 million years ago, crucial for the future development of novel control strategies against both species.

Dicrocoeliidae Family: Major Species Causing Veterinary Diseases

This chapter analyses the taxonomic position of Dicrocoeliidae family and several of its genus and species. The biology of the major species causing veterinary diseases such Dicrocoelium dendriticum, Dicrocoelium hospes, Dicrocoelium chinensis, Eurytrema pancreaticum and Platynosomum fastosum, has been reviewed. All these species have an indirect life cycle, involving two intermediate hosts (molluscs as first and ants, grasshoppers and lizards as second). Dicrocoelium dendriticum is a very widespread hepatic trematode in the ruminants of many countries in Europe, Asia, North Africa and North America, even affecting humans. Dicrocoelium hospes is widely distributed in the savanna areas of Africa south of the Sahara, whilst D. chinensis has mainly been found in ruminants in East Asia and some European countries (probably imported from Asia). Eurytrema pancreaticum is a common parasite whose adults live in ruminant bile ducts, gall bladder, pancreatic ducts and intestines in Europe, Madagascar, Asia and South America. Adult P. fastosum live in the liver, gall bladder and pancreas of birds and mammals in Europe, Africa, Asia, North, Central and South America. Information about the epidemiology, pathology, clinical aspect, diagnosis, treatment, control, prevention and economic impact mainly of Dicrocoeliosis produced by D. dendriticum, as well as of Eurytrematodosis and Platynosomiosis, has been included.

De novo assembly and functional annotations of the transcriptome of Metorchis orientalis (trematoda: Opisthorchiidae)

Metorchis orientalis is a neglected zoonotic parasite, living in the gallbladder and bile duct of poultry and some mammals as well as humans. In spite of its economic and medical importance, the information known about the transcriptome and genome of M. orientalis is limited. In this study, we performed de novo sequencing, transcriptome assembly and functional annotations of the adult M. orientalis, obtained about 77.4 million high-quality clean reads, among which the length of the transcript contigs ranged from 100 to 11,249 nt with mean length of 373 nt and N50 length of 919 nt. We then assembled 31,943 unigenes, of which 20,009 (62.6%) were annotated by BLASTn and BLASTx searches against the available database. Among these unigenes, 19,795 (62.0%), 3407 (10.7%), 10,620 (33.2%) of them had significant similarity in the NR, NT and Swiss-Prot databases, respectively; 5744 (18.0%) and 4678 (14.6%) unigenes were assigned to GO and COG, respectively; and 9099 (28.5%) unigenes were identified and mapped onto 256 pathways in the KEGG Pathway database. Furthermore, we found that 98 (1.08%) unigenes were related to bile secretion and 5 (0.05%) to primary bile acid biosynthesis pathways category. The characterization of these transcriptomic data has implications for the better understanding of the biology of M. orientalis, and will facilitate the development of intervention agents for this and other pathogenic flukes of human and animal health significance.

Characterization of the complete nuclear ribosomal DNA sequences of Eurytrema pancreaticum

Eurytrema pancreaticum is one of the most common trematodes of cattle and sheep, and also infects humans occasionally, causing great economic losses and medical costs. In this study, the sequences of the complete nuclear ribosomal DNA (rDNA) repeat units of five E. pancreaticum individuals were determined for the first time. They were 8306–8310 bp in length, including the small subunit (18S) rDNA, internal transcribed spacer 1 (ITS1), 5.8S rDNA, internal transcribed spacer 2 (ITS2), large subunit (28S) rDNA and intergenic spacer (IGS). There were no length variations in any of the investigated 18S (1996 bp), ITS1 (1103 bp), 5.8S (160 bp), ITS2 (231 bp) or 28S (3669 bp) rDNA sequences, whereas the IGS rDNA sequences of E. pancreaticum had a 4-bp length variation, ranging from 1147 to 1151 bp. The intraspecific variations within E. pancreaticum were 0–0.2% for 18S rDNA, 0–0.5% for ITS1, 0% for 5.8S rDNA and ITS2, 0–0.2% for 28S rDNA and 2.9–20.2% for IGS. There were nine types of repeat sequences in ITS1, two types in 28S rDNA, but none in IGS. A phylogenetic analysis based on the 18S rDNA sequences classified E. pancreaticum in the family Dicrocoeliidae of Plagiorchiata, closely related to the suborder Opisthorchiata. These results provide useful information for the further study of Dicrocoeliidae trematodes.

De novo transcriptome sequencing and analysis of the juvenile and adult stages of Fasciola gigantica

Fasciola gigantica is regarded as the major liver fluke causing fasciolosis in livestock in tropical countries. Despite the significant economic and public health impacts of F. gigantica there are few studies on the pathogenesis of this parasite and our understanding is further limited by the lack of genome and transcriptome information. In this study, de novo Illumina RNA sequencing (RNA-seq) was performed to obtain a comprehensive transcriptome profile of the juvenile (42days post infection) and adult stages of F. gigantica. A total of 49,720 unigenes were produced from juvenile and adult stages of F. gigantica, with an average length of 1286 nucleotides (nt) and N50 of 2076nt. A total of 27,862 (56.03%) unigenes were annotated by BLAST similarity searches against the NCBI non-redundant protein database. Because F. gigantica needs to feed and/or digest host tissues, some proteases (including cysteine proteases and aspartic proteases), which play a role in the degradation of host tissues (protein), have been paid more attention in the present study. A total of 6511 distinct genes were found differentially expressed between juveniles and adults, of which 3993 genes were up-regulated and 2518 genes were down-regulated in adults versus juveniles, respectively. Moreover, stage-specific differentially expressed genes were identified in juvenile (17,009) and adult (6517) F. gigantica. The significantly divergent pathways of differentially expressed genes included cAMP signaling pathway (226; 4.12%), proteoglycans in cancer (256; 4.67%) and focal adhesion (199; 3.63%). The transcription pattern also revealed two egg-laying-associated pathways: cGMP-PKG signaling pathway and TGF-β signaling pathway. This study provides the first comparative transcriptomic data concerning juvenile and adult stages of F. gigantica that will be of great value for future research efforts into understanding parasite pathogenesis and developing vaccines against this important parasite.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

De novo transcriptomic analysis of the female and male adults of the blood fluke Schistosoma turkestanicum

Background

Schistosoma turkestanicum is a parasite of considerable veterinary importance as an agent of animal schistosomiasis in many countries, including China. The S. turkestanicum cercariae can also infect humans, causing cercarial dermatitis in many countries and regions of the world. In spite of its significance as a pathogen of animals and humans, there is little transcriptomic and genomic data in the public databases.

Methods

Herein, we performed the transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females of S. turkestanicum and de novo transcriptome assembly.

Results

Approximately 81.1 (female) and 80.5 (male) million high-quality clean reads were obtained and then 29,526 (female) and 41,346 (male) unigenes were assembled. A total of 34,624 unigenes were produced from S. turkestanicum females and males, with an average length of 878 nucleotides (nt) and N50 of 1480 nt. Of these unigenes, 25,158 (72.7 %) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 21,995 (63.5 %), 22,189 (64.1 %) and 13,754 (39.7 %) of the unigenes had significant similarity in the NCBI non-redundant protein (NR), non-redundant nucleotide (NT) and Swiss-Prot databases, respectively. In addition, 3150 unigenes were identified to be expressed specifically in females and 1014 unigenes were identified to be expressed specifically in males. Interestingly, several pathways associated with gonadal development and sex maintenance were found, including the Wnt signaling pathway (103; 2 %) and progesterone-mediated oocyte maturation (77; 1.5 %).

Conclusions

The present study characterized and compared the transcriptomes of adult female and male blood fluke, S. turkestanicum. These results will not only serve as valuable resources for future functional genomics studies to understand the molecular aspects of S. turkestanicum, but also will provide essential information for ongoing whole genome sequencing efforts on this pathogenic blood fluke.