Expression of translationally controlled tumor protein (TCTP) gene of Dirofilaria immitis guided by transcriptomic screening.

Dirofilaria immitis (heartworm) infections affect domestic dogs, cats, and various wild mammals with increasing incidence in temperate and tropical areas. More sensitive antibody detection methodologies are required to diagnose asymptomatic dirofilariasis with low worm burdens. Applying current transcriptomic technologies would be useful to discover potential diagnostic markers for D. immitis infection. A filarial homologue of the mammalian translationally controlled tumor protein (TCTP) was initially identified by screening the assembled transcriptome of D. immitis (DiTCTP). A BLAST analysis suggested that the DiTCTP gene shared the highest similarity with TCTP from Loa loa at protein level (97%). A histidine-tagged recombinant DiTCTP protein (rDiTCTP) of 40 kDa expressed in Escherichia coli BL21 (DE3) showed immunoreactivity with serum from a dog experimentally infected with heartworms. Localization studies illustrated the ubiquitous presence of rDiTCTP protein in the lateral hypodermal chords, dorsal hypodermal chord, muscle, intestine, and uterus in female adult worms. Further studies on D. immitis-derived TCTP are warranted to assess whether this filarial protein could be used for a diagnostic purpose.

Identification of Dirofilaria immitis miRNA using illumina deep sequencing.

The heartworm Dirofilaria immitis is the causative agent of cardiopulmonary dirofilariosis in dogs and cats, which also infects a wide range of wild mammals and humans. The complex life cycle of D. immitis with several developmental stages in its invertebrate mosquito vectors and its vertebrate hosts indicates the importance of miRNA in growth and development, and their ability to regulate infection of mammalian hosts. This study identified the miRNA profiles of D. immitis of zoonotic significance by deep sequencing. A total of 1063 conserved miRNA candidates, including 68 anti-sense miRNA (miRNA*) sequences, were predicted by computational methods and could be grouped into 808 miRNA families. A significant bias towards family members, family abundance and sequence nucleotides was observed. Thirteen novel miRNA candidates were predicted by alignment with the Brugia malayi genome. Eleven out of 13 predicted miRNA candidates were verified by using a PCR-based method. Target genes of the novel miRNA candidates were predicted by using the heartworm transcriptome dataset. To our knowledge, this is the first report of miRNA profiles in D. immitis, which will contribute to a better understanding of the complex biology of this zoonotic filarial nematode and the molecular regulation roles of miRNA involved. Our findings may also become a useful resource for small RNA studies in other filarial parasitic nematodes.

Identification of neglected cestode Taenia multiceps microRNAs by illumina sequencing and bioinformatic analysis.

BACKGROUND: Worldwide, but especially in developing countries, coenurosis of sheep and other livestock is caused by Taenia multiceps larvae, and zoonotic infections occur in humans. Infections frequently lead to host death, resulting in huge socioeconomic losses. MicroRNAs (miRNAs) have important roles in the post-transcriptional regulation of a large number of animal genes by imperfectly binding target mRNAs. To date, there have been no reports of miRNAs in T. multiceps. RESULTS: In this study, we obtained 12.8 million high quality raw reads from adult T. multiceps small RNA library using Illumina sequencing technology. A total of 796 conserved miRNA families (containing 1,006 miRNAs) from 170,888 unique miRNAs were characterized using miRBase (Release 17.0). Here, we selected three conserved miRNA/miRNA* (antisense strand) duplexes at random and amplified their corresponding precursors using a PCR-based method. Furthermore, 20 candidate novel miRNA precursors were verified by genomic PCR. Among these, six corresponding T. multiceps miRNAs are considered specific for Taeniidae because no homologs were found in other species annotated in miRBase. In addition, 181,077 target sites within T. multiceps transcriptome were predicted for 20 candidate newly miRNAs. CONCLUSIONS: Our large-scale investigation of miRNAs in adult T. multiceps provides a substantial platform for improving our understanding of the molecular regulation of T. multiceps and other cestodes development.

Cloning and characterization of the fatty acid-binding protein gene from the protoscolex of Taenia multiceps.

Taenia multiceps (Cestoda: Taeniidae), a worldwide cestode parasite, is emerging as an important helminthic zoonosis due to serious or fatal central nervous system disease commonly known as coenurosis in domestic and wild ruminants including humans. Herein, a fatty acid-binding protein (FABP) gene was identified from transcriptomic data in T. multiceps. This gene, which contains a complete coding sequence, was amplified by reverse-transcriptase polymerase chain reaction. The corresponding protein, which was named TmFABP, had a molecular weight of 14 kDa, and subsequently was recombinantly expressed in Escherichia coli. The fusion protein was purified on Ni-NTA beads (Bio-Rad). Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and Western blot analyses showed that the purified recombinant protein caused immunogenicity. Immunohistochemical studies showed that TmFABP was expressed at the tegumental level in the protoscolices and in the cells between the body wall and parenchyma layer of the cestode. In sections from gravid proglottids, intense staining was detected in the uterus and eggs. Based on this, TmFABP could be switched on during differentiation of germinative layers to protoscoleces and from metacestodes to adult worms. Taken together, our results already reported for T. multiceps suggest the possibility of TmFABP developing a vaccine to control and prevent coenurosis.

Genetic variation of Taenia pisiformis collected from Sichuan, China, based on the mitochondrial cytochrome B gene.

Taenia pisiformis is one of the most important parasites of canines and rabbits. T. pisiformis cysticercus (the larval stage) causes severe damage to rabbit breeding, which results in huge economic losses. In this study, the genetic variation of T. pisiformis was determined in Sichuan Province, China. Fragments of the mitochondrial cytochrome b (cytb) (922 bp) gene were amplified in 53 isolates from 8 regions of T. pisiformis. Overall, 12 haplotypes were found in these 53 cytb sequences. Molecular genetic variations showed 98.4% genetic variation derived from intra-region. FST and Nm values suggested that 53 isolates were not genetically differentiated and had low levels of genetic diversity. Neutrality indices of the cytb sequences showed the evolution of T. pisiformis followed a neutral mode. Phylogenetic analysis revealed no correlation between phylogeny and geographic distribution. These findings indicate that 53 isolates of T. pisiformis keep a low genetic variation, which provide useful knowledge for monitoring changes in parasite populations for future control strategies.

Comparative efficacy of ivermectin and levamisole for reduction of migrating and encapsulated larvae of Baylisascaris transfuga in mice.

The comparative efficacy of 2 anthelmintics (ivermectin and levamisole) against Baylisascaris transfuga migrating and encapsulated larvae was studied in mice. A total of 60 BALB/c mice inoculated each with about 1,000 embryonated B. transfuga eggs were equally divided into 6 groups (A-F) randomly. Mice of groups A and B were treated with ivermectin and levamisole, respectively, on day 3 post-infection (PI). Mice of groups A-C were killed on day 13 PI. Similarly, groups D and E were treated with ivermectin and levamisole, respectively, on day 14 PI, and all mice of groups D-F were treated on day 24 PI. The groups C and F were controls. Microexamination was conducted to count the larvae recovering from each mouse. The percentages of reduction in the number of migrating larvae recovered from group A (ivermectin) and B (levamisole) were 88.3% and 81.1%, respectively. In addition, the reduction in encapsulated larvae counts achieved by ivermectin (group D) and levamisole (group E) was 75.0% and 49.2%, respectively. The results suggested that, to a certain extent, both anthelmintics appeared to be more effective against migrating larvae than encapsulated larvae. However, in the incipient stage of infection, ivermectin may be more competent than levamisole as a larvicidal drug for B. transfuga.

An ELISA using recombinant TmHSP70 for the diagnosis of Taenia multiceps infections in goats.

Infections with the tapeworm Taenia multiceps are problematic for ruminant farming worldwide. Here we develop a novel and rapid method for serodiagnosis of T. multiceps infections via an indirect ELISA (iELISA) that uses a heat shock protein, namely, TmHSP70. We extracted the total RNA of T. multiceps from the protoscoleces of cysts dissected from the brains of infected goats. Subsequently, we successfully amplified, cloned and expressed the TmHSP70 gene in Escherichia coli BL21 (DE3). Western blot analysis showed that the recombinant protein (∼34 kDa molecular weight) was recognized by the coenurosis positive serum. Given these initial, robust immunogenic properties for recombinant TmHSP protein, we assessed the ELISA-based serodiagnostic potential of this gene. The indirect ELISA was then optimized to 2.70 µg/well dilution for antigen and 1:80 dilution for serum,while the cut-off value is 0.446. We report that our novel TmHSP ELISA detected T. multiceps sera with a sensitivity of 1:10240 and a specificity of 83.3% (5/6). In a preliminary application, this assay correctly confirmed T. multiceps infection in 30 infected goats, consistent with the clinical examination. This study has revealed that our novel iELISA, which uses the rTmHSP protein, provides a rapid test for diagnosing coenurosis.

Genetic variability of Echinococcus granulosus from the Tibetan plateau inferred by mitochondrial DNA sequences.

To analyse genetic variability and population structure, 84 isolates of Echinococcus granulosus (Cestoda: Taeniidae) collected from various host species at different sites of the Tibetan plateau in China were sequenced for the whole mitochondrial nad1 (894 bp) and atp6 (513 bp) genes. The vast majority were classified as G1 genotype (n=82), and two samples from human patients in Sichuan province were identified as G3 genotype. Based on the concatenated sequences of nad1+atp6, 28 different haplotypes (NA1-NA28) were identified. A parsimonious network of the concatenated sequence haplotypes showed star-like features in the overall population, with NA1 as the major haplotype in the population networks. By AMOVA it was shown that variation of E. granulosus within the overall population was the main pattern of the total genetic variability. Neutrality indexes of the concatenated sequence (nad1+atp6) were computed by Tajima's D and Fu's Fs tests and showed high negative values for E. granulosus, indicating significant deviations from neutrality. FST and Nm values suggested that the populations were not genetically differentiated.

Detailed transcriptome description of the neglected cestode Taenia multiceps.

BACKGROUND: The larval stage of Taenia multiceps, a global cestode, encysts in the central nervous system (CNS) of sheep and other livestock. This frequently leads to their death and huge socioeconomic losses, especially in developing countries. This parasite can also cause zoonotic infections in humans, but has been largely neglected due to a lack of diagnostic techniques and studies. Recent developments in next-generation sequencing provide an opportunity to explore the transcriptome of T. multiceps. METHODOLOGY/PRINCIPAL FINDINGS: We obtained a total of 31,282 unigenes (mean length 920 bp) using Illumina paired-end sequencing technology and a new Trinity de novo assembler without a referenced genome. Individual transcription molecules were determined by sequence-based annotations and/or domain-based annotations against public databases (Nr, UniprotKB/Swiss-Prot, COG, KEGG, UniProtKB/TrEMBL, InterPro and Pfam). We identified 26,110 (83.47%) unigenes and inferred 20,896 (66.8%) coding sequences (CDS). Further comparative transcripts analysis with other cestodes (Taenia pisiformis, Taenia solium, Echincoccus granulosus and Echincoccus multilocularis) and intestinal parasites (Trichinella spiralis, Ancylostoma caninum and Ascaris suum) showed that 5,100 common genes were shared among three Taenia tapeworms, 261 conserved genes were detected among five Taeniidae cestodes, and 109 common genes were found in four zoonotic intestinal parasites. Some of the common genes were genes required for parasite survival, involved in parasite-host interactions. In addition, we amplified two full-length CDS of unigenes from the common genes using RT-PCR. CONCLUSIONS/SIGNIFICANCE: This study provides an extensive transcriptome of the adult stage of T. multiceps, and demonstrates that comparative transcriptomic investigations deserve to be further studied. This transcriptome dataset forms a substantial public information platform to achieve a fundamental understanding of the biology of T. multiceps, and helps in the identification of drug targets and parasite-host interaction studies.

Annotation of the transcriptome from Taenia pisiformis and its comparative analysis with three Taeniidae species.

BACKGROUND: Taenia pisiformis is one of the most common intestinal tapeworms and can cause infections in canines. Adult T. pisiformis (canines as definitive hosts) and Cysticercus pisiformis (rabbits as intermediate hosts) cause significant health problems to the host and considerable socio-economic losses as a consequence. No complete genomic data regarding T. pisiformis are currently available in public databases. RNA-seq provides an effective approach to analyze the eukaryotic transcriptome to generate large functional gene datasets that can be used for further studies. METHODOLOGY/PRINCIPAL FINDINGS: In this study, 2.67 million sequencing clean reads and 72,957 unigenes were generated using the RNA-seq technique. Based on a sequence similarity search with known proteins, a total of 26,012 unigenes (no redundancy) were identified after quality control procedures via the alignment of four databases. Overall, 15,920 unigenes were mapped to 203 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Through analyzing the glycolysis/gluconeogenesis and axonal guidance pathways, we achieved an in-depth understanding of the biochemistry of T. pisiformis. Here, we selected four unigenes at random and obtained their full-length cDNA clones using RACE PCR. Functional distribution characteristics were gained through comparing four cestode species (72,957 unigenes of T. pisiformis, 30,700 ESTs of T. solium, 1,058 ESTs of Eg+Em [conserved ESTs between Echinococcus granulosus and Echinococcus multilocularis]), with the cluster of orthologous groups (COG) and gene ontology (GO) functional classification systems. Furthermore, the conserved common genes in these four cestode species were obtained and aligned by the KEGG database. CONCLUSION: This study provides an extensive transcriptome dataset obtained from the deep sequencing of T. pisiformis in a non-model whole genome. The identification of conserved genes may provide novel approaches for potential drug targets and vaccinations against cestode infections. Research can now accelerate into the functional genomics, immunity and gene expression profiles of cestode species.

Novel insights into the transcriptome of Dirofilaria immitis.

BACKGROUND: The heartworm Dirofilaria immitis is the causal agent of cardiopulmonary dirofilariosis in dogs and cats, and also infects a wide range of wild mammals as well as humans. One bottleneck for the design of fundamentally new intervention and management strategies against D. immitis may be the currently limited knowledge of fundamental molecular aspects of D. immitis. METHODOLOGY/PRINCIPAL FINDINGS: A next-generation sequencing platform combining computational approaches was employed to assess a global view of the heartworm transcriptome. A total of 20,810 unigenes (mean length  =  1,270 bp) were assembled from 22.3 million clean reads. From these, 15,698 coding sequences (CDS) were inferred, and about 85% of the unigenes had orthologs/homologs in public databases. Comparative transcriptomic study uncovered 4,157 filarial-specific genes as well as 3,795 genes potentially involved in filarial-Wolbachia symbiosis. In addition, the potential intestine transcriptome of D. immitis (1,101 genes) was mined for the first time, which might help to discover 'hidden antigens'. CONCLUSIONS/SIGNIFICANCE: This study provides novel insights into the transcriptome of D. immitis and sheds light on its molecular processes and survival mechanisms. Furthermore, it provides a platform to discover new vaccine candidates and potential targets for new drugs against dirofilariosis.

Complete mitochondrial genomes of Baylisascaris schroederi, Baylisascaris ailuri and Baylisascaris transfuga from giant panda, red panda and polar bear.

Roundworms of the genus Baylisascaris are the most common parasitic nematodes of the intestinal tracts of wild mammals, and most of them have significant impacts in veterinary and public health. Mitochondrial (mt) genomes provide a foundation for studying epidemiology and ecology of these parasites and therefore may be used to assist in the control of Baylisascariasis. Here, we determined the complete sequences of mtDNAs for Baylisascaris schroederi, Baylisascaris ailuri and Baylisascaris transfuga, with 14,778 bp, 14,657 bp and 14,898 bp in size, respectively. Each mtDNA encodes 12 protein-coding genes, 22 transfer RNAs and 2 ribosomal RNAs, typical for other chromadorean nematodes. The gene arrangements for the three Baylisascaris species are the same as those of the Ascaridata species, but radically different from those of the Spirurida species. Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes from nine nematode species indicated that the three Baylisascaris species are more closely related to Ascaris suum than to the three Toxocara species (Toxocara canis, Toxocara cati and Toxocara malaysiensis) and Anisakis simplex, and that B. ailuri is more closely related to B. transfuga than to B. schroeder. The determination of the complete mt genome sequences for these three Baylisascaris species (the first members of the genus Baylisascaris ever sequenced) is of importance in refining the phylogenetic relationships within the order Ascaridida, and provides new molecular data for population genetic, systematic, epidemiological and ecological studies of parasitic nematodes of socio-economic importance in wildlife.

The mitochondrial genome of Baylisascaris procyonis.

BACKGROUND: Baylisascaris procyonis (Nematoda: Ascaridida), an intestinal nematode of raccoons, is emerging as an important helminthic zoonosis due to serious or fatal larval migrans in animals and humans. Despite its significant veterinary and public health impact, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. Mitochondrial (mt) genomes can provide a foundation for investigations in these areas and assist in the diagnosis and control of B. procyonis. In this study, the first complete mt genome sequence of B. procyonis was determined using a polymerase chain reaction (PCR)-based primer-walking strategy. METHODOLOGY/PRINCIPAL FINDINGS: The circular mt genome (14781 bp) of B. procyonis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes congruent with other chromadorean nematodes. Interestingly, the B. procyonis mtDNA featured an extremely long AT-rich region (1375 bp) and a high number of intergenic spacers (17), making it unique compared with other secernentean nematodes characterized to date. Additionally, the entire genome displayed notable levels of AT skew and GC skew. Based on pairwise comparisons and sliding window analysis of mt genes among the available 11 Ascaridida mtDNAs, new primer pairs were designed to amplify specific short fragments of the genes cytb (548 bp fragment) and rrnL (200 bp fragment) in the B. procyonis mtDNA, and tested as possible alternatives to existing mt molecular beacons for Ascaridida. Finally, phylogenetic analysis of mtDNAs provided novel estimates of the interrelationships of Baylisasaris and Ascaridida. CONCLUSIONS/SIGNIFICANCE: The complete mt genome sequence of B. procyonis sequenced here should contribute to molecular diagnostic methods, epidemiological investigations and ecological studies of B. procyonis and other related ascaridoids. The information will be important in refining the phylogenetic relationships within the order Ascaridida and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of parasitic nematodes of socio-economic importance.