Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 05 July 2021
Sec. Parasitology

Comparative Transcriptome Analyses of the Developmental Stages of Taenia multiceps

\nWen-Hui Li
Wen-Hui Li1*Yang YangYang Yang2Nian-Zhang ZhangNian-Zhang Zhang1Jian-Kui WangJian-Kui Wang3Yin-Ju LiuYin-Ju Liu1Li LiLi Li1Hong-Bin YanHong-Bin Yan1Wan-Zhong JiaWan-Zhong Jia1Baoquan Fu,
Baoquan Fu1,4*
  • 1State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China
  • 2Wuhan Animal Disease Prevention and Control Center, Wuhan, China
  • 3Center for Animal Disease Control and Prevention of Qilihe District, Lanzhou, China
  • 4Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou, China

Cerebral coenurosis, caused by the larvae of Taenia multiceps (Coenurus cerebralis), is a fatal central nervous system disease in sheep and other herbivores and occasionally humans. Comparative transcriptomic profiles of the developmental stages of the parasite remain unknown. In this study, RNA sequencing was used to determine the transcriptome profiles of different stages of the life cycle of T. multiceps, including Oncosphere, Coenurus cerebralis (Pro with Cyst), and Adult (Adu), as well as scolex-neck proglottids (Snp), immature–mature proglottids (Imp), and gravid proglottids (Grp) of the adult stage. A total of 42.6 Gb (average 6.1 Gb) Illumina pair-end reads with a 125-bp read length were generated for seven samples. The total number of differentially expressed genes (DEGs) in the various life stages ranged from 2,577 to 3,879; however, for the tissues of the adult worm, the range was from 1,229 to 1,939. Kyoto Encyclopedia of Genes and Genomes analysis showed that the DEGs mainly participated in cellular and metabolic processes, binding and catalytic activity, genetic information processing, and environmental information processing. In addition, a large number of genes related to development and parasite–host interaction were identified. Quantitative reverse transcription-polymerase chain reaction confirmed that the levels of 28 selected DEGs were consistent with those determined using RNA sequencing. The present study provides insights into the mechanisms of the development and parasitic life of T. multiceps.

Keypoints

- A total of 42.6 Gb (average 6.1 Gb) Illumina pair-end reads with a 125-bp read length were generated for seven samples from T. multiceps.

- The total number of DEGs in the various life stages of T. multiceps ranged from 2,577 to 3,879; however, for the tissues of the adult worm, the range was from 1,229 to 1,939.

- Some differentially expressed genes are involved in the development of T. multiceps and the parasite–host interaction.

Introduction

Taenia multiceps is a tapeworm that inhabits the small intestine of dogs and other canids. The larval stage (known as Coenurus cerebralis) with a fluid-filled cyst is usually found in the nervous system of the intermediate host and often causes encephalitis and neurological symptoms in infected hosts, and sometimes even death (1, 2). Definitive hosts excrete egg-containing feces into the environment; when intermediate herbivorous hosts ingest these eggs during grazing, oncospheres spread from the eggs. Furthermore, through blood circulation, they reach various locations, including the central nervous system, where they could develop into mature Coenurus cysts of the larval stage within 2–3 months (3). The life cycle of T. multiceps is completed after the scoleces contained within the coenurus are ingested by the definitive host. The disease may last several months and is fatal (4, 5). T. multiceps infection in sheep results in significant economic losses in many parts of the world, predominantly in the developing countries of Africa, Europe, and Southeast Asia (1, 6, 7). Cases of human infections due to accidental ingestion of T. multiceps eggs have also been reported frequently and cause serious pathological conditions (812). Cestode development is well-known to involve complex morphological and physiological changes associated with changes in the host species (13, 14). Strobilation development (generation of serially repeated reproductive organs) is an interesting but difficult study subject, as it is difficult to obtain the parasites and maintain them in experimental hosts or in vitro (15). Recently, the genome of T. multiceps was sequenced to understand its biology and evolution (16). Extensive transcriptome analysis of T. multiceps larvae (17) and adults (18) provides substantial and useful information. Additionally, comparative transcriptomic analysis can identify key genes associated with the development, structure, and storage of carbohydrates, thereby contributing to the understanding of the biology of organisms (1921). However, the mechanisms underlying the development of T. multiceps from larva to adult remain unclear. In the present study, comparative transcriptome analyses of the developmental stages of T. multiceps were performed to obtain important insights into the molecular basis of biological processes involved in worm development and reproduction.

Materials and Methods

Sample Collection

(i) The larvae of T. multiceps (named Coenurus cerebralis) were obtained from the brains of naturally infected sheep at a local slaughterhouse in Gansu province, China. Freshly separated Coenurus cerebralis were washed repeatedly in sterile phosphate-buffered saline (PBS, pH 7.4) to remove the host components. Protoscolices (Pro) were then separated from the cyst wall (cyst) of Coenurus cerebralis. (ii) Adult T. multiceps (Adu) were collected from experimentally infected dogs at our laboratory animal center. The dog infection experiment was performed according to a previously reported protocol (18, 22). In brief, two 4-month-old parasite-free beagle dogs were infected orally with 40 T. multiceps protoscolices. When gravid proglottids (Grp) were found in the feces, the infected dogs were exsanguinated, and the worms were collected and washed. Thereafter, scolex-neck proglottids (Snp), immature–mature proglottids (Imp), and Grp were separated from each adult based on the proglottid characteristics. (iii) For the isolation of oncospheres, gravid tapeworm segments were obtained from the above collected adult worms, which were then homogenized, filtered, washed, and purified by Percoll density-gradient separation. Finally, oncospheres (Onc) were hatched according to the procedure described previously (23, 24). To prepare tissue samples from adult T. multiceps, intact adults were segmented into three parts: Snp, Imp, and Grp based on their proglottid characteristics. All prepared samples were washed with sterile PBS, preserved in five volumes of RNA-later (TaKaRa, Japan), and then immediately stored in liquid nitrogen at −80°C for RNA sequencing (RNA-seq) analysis.

Total RNA Isolation, Qualification, and Transcriptomic Library Construction

Total RNA was extracted from each sample using TRIzol Reagent (Invitrogen, USA), and contaminating genomic DNA was removed by treatment with RNase-free DNase I (Takara, Japan). RNA degradation and contamination were monitored on 1% agarose gels, and RNA integrity was assessed using an RNA Nano 6000 Assay Kit and Bioanalyzer 2100 system (Agilent Technologies, USA). RNA with an RNA integrity number (RIN) > 7.0 was considered to be of high enough quality for transcriptomic library construction and RNA sequencing, according to the manufacturer's instructions. Subsequently, mRNA was isolated from total RNA, cut into fragments, and reverse-transcribed into first-strand cDNA using random hexamers. Briefly, double-stranded cDNA was synthesized and purified, and the ends were repaired. Finally, RNA sequencing libraries were constructed using the Illumina TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA), according to the manufacturer's instructions. Library quality was assessed using an Agilent Bioanalyzer 2100 system (Agilent Technologies).

Transcriptome Sequencing and Data Analysis

Seven transcriptomic libraries were pooled and sequenced using the Illumina HiSeq 2000 platform (Illumina, San Diego, CA, USA). Raw reads in fastq format were first processed through in-house Perl scripts to obtain clean reads by discarding reads containing adapters, reads containing poly-N, and low-quality reads. Finally, high-quality reads were generated after RNA-seq and filtering using the FastQC software (25). The filtered reads were then aligned to the draft of the T. multiceps genome sequence, the indexes of Burrows–Wheeler were constructed using Bowtie2 in default parameters (26), and the clean data from seven RNA-seq libraries were mapped to the assembled genome of T. multiceps using the BWT-MEM software.

To assess the differentially expressed genes (DEGs) in each life-cycle stage and tissue, transcript FPKM values from a replicate were plotted against transcript FPKM values from another replicate using cutflinks (27). The false discovery rate (FDR) was used to determine the threshold P-values in multiple tests for comparing the differential expression. We used |log2 (FPKMsampe1/FPKMsampe2)| ≥ 1 and FDR ≤ 0.001 as the threshold to judge the significance of DEGs.

Functional Classification of DEGs

Functional annotation of DEGs was performed using the UniProt database and BLASTP v2.3.0+ (28) and InterProScan (29) analyses. Gene ontology (GO) enrichment analysis of DEGs was performed using the GOseq R package (v1.10.0).

Data Access

The RNA-seq data from this study have been deposited in NCBI Sequence Read Archive (SRA accession number PRJNA307624) and are accessible through accession numbers SRX2357663, SRX2357664, SRX2357665, SRX2357666, SRX2357667, SRX2357668, and SRX2357669.

Quantitative Amplification to Verify DEGs Identified Using RNA-Seq

qPCR assay was performed to validate the DEGs identified using RNA-seq. Twenty-eight DEGs were selected and investigated using qRT-PCR at the transcriptional level. The β-tubulin gene (Tmu007343) was used as an internal control. All primers were designed using Primer Premier 5.0 and synthesized by Sangon (Shanghai, China) (Table 1). Total RNA was isolated as described above for the samples from protoscolices (Pro), cyst wall (Cyst), adult (Adu), scolex and neck proglottids (Snp), Imp, and Grp. cDNA was synthesized using the PrimeScript RT Reagent Kit with gDNA Eraser (TaKaRa). The qRT-PCR was performed on an ABI Stepone plus real-time PCR Detection System following the manufacturer's protocol for SG Fast qPCR Master Mix (High Rox) (2X) (BBI) (ABI). The PCR conditions were as follows: an initial denaturation step at 95°C for 3 min, followed by 45 cycles of denaturation at 95°C for 7 s, annealing at 57°C for 10 s, and extension at 72°C for 15 s. Relative quantification of all selected genes was performed using the 2−ΔΔCT method.

TABLE 1
www.frontiersin.org

Table 1. qRT-PCR primers of differentially expressed genes.

Results

Summary of RNA-Seq Data and Transcript Assembly

The seven RNA samples from the developmental stages of T. multiceps had prominent 18S and 28S ribosomal peaks on agarose gels (data not shown). All RNA samples possessed high integrity and purity with RINs>7.0, according to Bioanalyzer 2100 analysis, and were considered of sufficient quality for further experiments. Raw read sequences were filtered to discard adaptors and low-quality reads (error rate > 0.001) from all seven RNA-seq libraries. The clean data from seven RNA-seq libraries with a matching ratio >97% were mapped to the genome of T. multiceps (Table 2). Overall, a total of 42.6 Gb (average 6.1 Gb) Illumina pair-end reads with a 125-bp read length were generated for seven samples, including oncospheres, protoscolices and cyst wall, Snp, Imp, and Grp. A total of 8,208 to 10,498 predicted genes in all samples with FPKM value > 1, 4,396–6,652 genes with FPKM value >10, 752 to 1, and 290 genes with FPKM value >10 were detected (Table 3).

TABLE 2
www.frontiersin.org

Table 2. RNA-seq reads for transcriptomes aligned to the assembled genome of T. multiceps.

TABLE 3
www.frontiersin.org

Table 3. Statistics of expression genes in seven RNA-seq samples with FPKM > 1, FPKM > 10, and FPKM > 100.

Identification and Annotation of DEGs in the Developmental Stages and Different Tissues of T. multiceps

To comprehensively understand the dynamics of gene expression in T. multiceps, we compared the transcriptome data between different life cycle stages and different tissues of adult worms. The total number of DEGs in the various life stages ranged from 2,577 to 3,879; however, for the tissues of the adult worm, the range was from 1,229 to 1,939 (Supplementary Table 1; Figures 1A,B).

FIGURE 1
www.frontiersin.org

Figure 1. Differentially expressed genes (DEGs) among the developmental stages and tissues of adult T. multiceps. (A) Numbers of DEGs in the different stages. (B) Numbers of DEGs in the different adult tissues.

A total of 2,437 DEGs were selected at different developmental stages (Figure 2A). Of these, 664 were oncosphere stage-specific and upregulated, including oncosphere protein Tm18 (Tmu005436), Tm16 (Tmu007899), 45m (Tmu003810), and fibronectin type III domain-containing protein (Tmu011598). Four hundred eighteen of the DEGs were protoscolices stage-specific and upregulated, including immunogenic proteins (Tmu003750), neuronal nitric oxide synthase protein inhibitor (Tmu000145), and inhibitor of apoptosis protein (Tmu004095). In addition, 237 of the DEGs were cyst wall-specific and upregulated, including 8 kDa glycoproteins (Tmu005741, Tmu005823), diagnostic antigen gp50 (Tmu00752), and Hsp90 (Tmu007808). In addition, 1,118 of the DEGs were adult-specific and upregulated, including the tegumental antigen (Tmu000235), ornithine aminotransferase (Tmu000506), heat shock protein 70 (Tmu01113), and tetraspanin (TSP) (Tmu000409). These phase-specific upregulated genes were annotated using GO. Furthermore, there were 456 specific upregulated genes during the oncosphere stage that could be annotated to 31 GO terms. Two hundred thirty-four protoscolices stage-specific upregulated genes were annotated to 32 GO terms, 60 cyst wall-specific upregulated genes to 23 GO terms, and 446 adult-specific upregulated genes to 31 GO terms.

FIGURE 2
www.frontiersin.org

Figure 2. Specific upregulated genes in each development stage and tissue of T. multiceps. (A) Specific upregulated genes in each developmental stage of T. multiceps. (B) Specific upregulated genes in Snp, Imp, and Grp of adult T. multiceps.

A total of 1,228 DEGs were selected from different developmental tissues (Figure 2B). Ninety-seven of the DEGs were Snp tissue-specific and upregulated, including diagnostic antigens (Tmu006941, Tmu012211, and Tmu012074), TSP (Tmu004996), and profoldin (Tmu010836). Eight hundred forty-seven of the DEGs were Imp tissue-specific and upregulated, including Tmu002740 (similar to ubiquitin/40S ribosomal protein S27a fusion protein) and putative thioredoxin-2 (Tmu000756), whereas 284 of the DEGs were Grp tissue-specific and upregulated, such as the fatty acid-binding protein (FABP) 3 (Tmu003178), insulin-like protein (Tmu004894), epidermal growth factor receptor kinase (Tmu008864), and oncosphere protein Tm18. These tissue-specific upregulated genes were annotated using GO. The 45 Snp tissue-specific upregulated genes were annotated to 21 GO terms, 328 Imp tissue-specific upregulated genes were annotated to 32 GO terms, and 124 Grp tissue-specific upregulated genes were annotated to 27 GO units.

To gain more insight into the specifically upregulated genes in each stage and tissues of the adult T. multiceps, the GO database was used to annotate the functions of these specific and upregulated genes. In summary, GO enrichment of the upregulated DEGs in stage-specific and tissue-specific conditions was analyzed. A total of 456 (68.7%), 234 (56.0%), and 446 (39.9%) specifically upregulated genes were assigned to one or more GO categories for Onc, Pro, and Adu, respectively. A total of 45 (46.4%), 328 (38.7%), and 124 (43.7%) specifically upregulated genes were assigned to one or more GO categories for Snp, Imp, and Grp, respectively, as shown in Figures 2A,B. Notably, the number of DEGs among the different stages was higher than that among the different tissues and was consistent with the extent of the phenotypic differences.

GO classifications enriched in genes for biological process were involved in cellular processes, metabolic processes, pigmentation, and biological regulation; GO classifications enriched in genes for molecular function process were involved in binding, catalytic, and transcription regulation; and GO classifications enriched in genes for the cellular component were involved in the regulation of macromolecular complexes and organelle development (Figure 3). The above observations suggest that the process of development from oncospheres to protoscolices and then to adults involves numerous molecular functions and requires extensive biological processes.

FIGURE 3
www.frontiersin.org

Figure 3. GO analysis of the specific genes upregulated in different developmental stages and tissues of T. multiceps. (A) Functional classification of specific upregulated genes in three stages (Onc, Pro, and Adu) of T. multiceps. (B) Functional classification of specific upregulated genes in three tissues of adult T. multiceps.

In addition, genes coding for an 8-kDa glycoprotein and the oncosphere proteins Tm16 and Tm18 showed obvious stage-specific variation in expression at different developmental stages of T. multiceps. In total, 15 of the 16 genes coding for the 8-kDa glycoprotein were expressed with FPKM > 1 (Supplementary Tables 2, 3).

DEGs Related to Protease Inhibitor

A total of 29 protease inhibitor-associated genes were found in the transcriptomes of T. multiceps and consisted mostly of a trypsin inhibitor, followed by serine and cysteine protease inhibitors, with varying expression levels in the different developmental stages. A comparative analysis of the gene expression levels between Adu, Pro, and Onc showed that Tmu006090, Tmu006091, and Tmu008198 were upregulated in Onc; Tmu010089, Tmu010102, Tmu004311, Tmu010010, Tmu008761, Tmu005998, Tmu002313, Tmu010090, and Tmu007272 were upregulated in Pro; whereas Tmu007271, Tmu012421, Tmu010103, and Tmu004242 were upregulated in Adu. However, Tmu002208, Tmu007199, Tmu005489, Tmu010101, Tmu003186, Tmu002646, Tmu008502, and Tmu004708 were highly expressed in Adu and Pro, in contrast to Onc (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. Heatmap of the expression of inhibitor-associated genes in the transcriptomes of T. multiceps.

KEGG pathway analysis indicated that most of the upregulated DEGs were significantly involved in metabolism, genetic information processing, cellular processes, and environmental information processing, such as glutathione metabolism, nitrogen metabolism, endoplasmic reticulum protein processing, endocytosis, and ubiquitin-mediated proteolysis. Furthermore, upregulated DEGs related to the Wnt signaling pathway, the mitogen-activated protein kinase (MAPK) signaling pathway, and the TGF-β signaling pathway were significantly abundant in the adult stage. Those related to inositol phosphate metabolism, citric acid cycle, and adhesion-related genes were significantly abundant in oncospheres, while progesterone-mediated oocyte maturation and oocyte meiosis-related DEGs were significantly abundant in Imp.

DEG Confirmation Using Real-Time PCR

Twenty-eight genes from the DEGs designated Tmu003811 (oncosphere protein Tso31b), Tmu005434 and Tmu012228 (oncosphere protein Tm18), Tmu007899 (oncosphere protein Tm16), Tmu005154 (8 kDa glycoprotein), Tmu003810 (45 m oncosphere protein), Tmu003173 (cytosolic FABP), Tmu000409 (TSP family), Tmu011357 and Tmu012976 (T-Box domain signature), Tmu004894 (insulin-like growth factor), Tmu007712 and Tmu010856 (Unknown), and heat shock protein (Tmu007808) were selected for qRT-PCR analysis to quantify their transcription levels (Figure 5). Of these genes, Tmu003811, Tmu005434, Tmu012228, Tmu007899, Tmu003810, Tmu012421, and Tmu006091 were downregulated in protoscoleces, cyst wall, and the adult stage compared with those in the oncosphere stage, whereas Tmu005042, Tmu003173, and Tmu003193 were upregulated. Meanwhile, Tmu005434, Tmu012228, Tmu003810, and Tmu004894 were upregulated in Grp compared with those in Snp and Imp. The qRT-PCT results of the 28 genes corroborated their RNA-seq results.

FIGURE 5
www.frontiersin.org

Figure 5. Relative expression levels of the 28 selected DEGs determined using qRT-PCR. (A,B) Confirmation of the DEGs in the developmental stages of T. multiceps. (C) Confirmation of the DEGs in different tissues in the adult stage of T. multiceps.

Discussion

The life cycle of T. multiceps requires two hosts (dog and sheep). When the final host ingests the mature Coenurus cysts, T. multiceps reaches maturity in its intestine around 40–42 days. Thereafter, the dog starts to disseminate 3–4 proglottids daily, which could contain almost 37,000 eggs each (30). T. multiceps eggs are usually discharged/released/disseminated to the environment through feces from the proglottids (31). When taken up by the intermediate host with contaminated food or water upon reaching the intestine, eggs with hexacanth larvae release oncospheres that penetrate the intestinal wall and migrate toward target organs (usually the brain) through the bloodstream (32). Why oncospheres have a predilection for certain sites, such as muscle, brain, and subcutaneous tissue, is not clear. The larval stage, known as Coenurus cerebralis, causes central nervous system disease in sheep, commonly known as coenurosis. The cyst in the brain of sheep is a large, white, translucent structure containing clear fluid with numerous protoscoleces attached to the wall (32). The biology of cestodes is unique in many aspects, especially regarding their unique development, which is difficult to compare to that of other organisms. In the present study, a comprehensive comparative transcriptome analysis of the different developmental stages of T. multiceps was carried out for the first time, and 2,437 DEGs were obtained, including oncosphere-specific proteins, immunogenic proteins, diagnostic antigens, heat shock proteins, and tegumental antigens. Podesta and Siddiqui (33) suggested that protein synthesis in flatworm parasites is controlled post-transcriptionally and that these intracellular regulatory mechanisms are activated or suppressed by effectors of the host's immune response. It is evident that this phase-specific or upregulated gene is important for adapting to parasitic hosts and certain sites. This is similar to our previous study, in which the lifecycle process of the oncosphere, larvae, and adults to adapt to the new host environment is a complex differentiation process that is regulated through various protein expressions (16).

Cestodes generally show impressive adaptation to parasitism, including the complete lack of a digestive system, complex life cycle that typically involves two or more hosts, specialized life stages for invasion of each host, specialized structures for attachment (such as an anterior, which usually contains attachment organs, such as grooves or suckers), and serial repetition of reproductive structures (14). A well-adapted scolex, an elongated neck, and a strobila are characteristic of adult T. multiceps. The scolex is a holdfast organ armed with four suckers and two circles of rostellar hooks numbered between 18 and 34 (34). The scolex appears to act as an organ of attachment and may also have nutritive and/or sensory functions. In the present study, three of the seven Hsp90 genes (Tmu009603, Tmu007811, and Tmu007806) were upregulated in the adult stage; some reports have suggested that the Hsp90 gene controls stage differentiation in Leishmania donovani (35). Moreover, the Hsp90 gene has been found to correlate with strobilization in Echinococcus granulosus (36). Therefore, we hypothesized that the three Hsp90 genes in T. multiceps correlate with strobilization in adult worm proglottids. Our study also revealed that FABP was upregulated in protoscoleces and the cyst wall. Similar to what we observed, Basika et al. (20) showed that FABP transcripts were abundant in Mesocestoides corti. FABP is known to participate in fatty acid transport to obtain essential lipids from the host and to modify them according to the parasite's needs, as cestodes lack the ability to synthesize fatty acids and cholesterol de novo (37). Nie et al. (38) confirmed that FABPs were mainly distributed in the muscle layer in the protoscolices and cells between the body wall and the parenchyma layer of the cestode in various life stages of T. multiceps, suggesting that the FABP is closely related to the formation or development of protoscolices and their adaptation to a parasitic environment. Of note, the TSP (Tmu000409) gene was significantly upregulated in adults. TSPs have been used as vaccine candidates against schistosomiasis and echinococcosis and as diagnostic antigens for cysticercosis (3942). Tran et al. (43) demonstrated that TSPs in the tegument of schistosomula and adult worms can act as receptors for host ligands, including MHC molecules, allowing parasites to mask their non-self-status and escape host immune responses. Furthermore, Huang et al. (44) identified 11 amino acid sequences in the TSP superfamily of E. multilocularis, and the Em-TSP1 homologs were found to be highly expressed in early-stage metacestodes compared with those in non-activated and activated oncospheres; however, the expression level of Em-TSP3 varied within non-activated oncosphere samples. In summary, the stage-specific expression of the TSP gene in T. multiceps requires further verification.

In the present study, three, nine, and four protease inhibitor-associated genes were upregulated in the oncosphere, protoscolices, and adult stage, respectively. Because tapeworms lack a gut, all nutritive materials must pass through the body surface in the same way as waste materials. The complex cestode tegument contains specific systems for molecular and ion transport, especially amino acids, hexose sugars, vitamins, purines, pyrimidines, nucleotides, and lipids. It might also serve as a major site for catalytic activities and contains enzymes of parasitic and possibly of host origin (45, 46). It is well-known that in these developmental stages, the tapeworm must resist the digestive action of gastric acid in the internal environment of the respective hosts to survive; therefore, the protease inhibitor-associated genes mentioned above may play a critical role in protecting them from the host digestion system.

We surprisingly found a highly expressed 8 kDa glycoprotein in Pro; some information is available concerning 8 kDa as a serodiagnosis antigen in T. solium cysticercosis (47). In addition, this 8 kDa glycoprotein from Echinococcus multilocularis has been identified as a potential diagnostic marker in patients with alveolar echinonococcosis (48). We predict that the 8 kDa glycoprotein may be outside for the serodiagnosis of T. multiceps coenurosis in the future.

Although KEGG pathways involved in metabolism, genetic information processing, cellular processes, and environmental information processing were identified in the present study, unique pathways were observed in the different stages of T. multiceps. For instance, genes associated with the Wnt, MAPK, and TGF-β signaling pathways were significantly abundant in the adult stage, which suggested that the growth and development of the adult tapeworm rely on the above signaling pathways as they play key roles in the regulation of cell growth, differentiation, and tissue and organ formation (49, 50).

Conclusions

In summary, we employed an Illumina RNA-Seq platform to identify DEGs in different developmental stages and tissues of T. multiceps. We found that the differential expression of some important genes may contribute to the development of T. multiceps and parasite–host interaction.

Data Availability Statement

The RNAseq data from this study have been deposited in NCBI's Sequence Read Archive and are accessible through accession NumberSRX2357663, SRX2357664, SRX2357665, SRX2357666, SRX2357667, SRX2357668, SRX2357669. https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR5032913.

Ethics Statement

The animal study was reviewed and approved by the Animal Administration and Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

BF and WHL conceived and designed the study. WHL wrote the manuscript. YY performed the experiments. NZZ analyzed the data. YJL and JKW assisted to collect samples. LL, HBY, and WZJ helped in the implementation of the study. All authors contributed to the article and approved the submitted version.

Funding

This study was financially supported by grants from the National Natural Science Foundation of China (Grant No. 31702228).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.677045/full#supplementary-material

Supplementary Table 1. Differentially expressed genes between the developmental stages or tissues.

Supplementary Table 2. Expression of the 8 kDa glycoprotein in the different developmental stages of T. multiceps.

Supplementary Table 3. Expression of the oncosphere proteins, Tm16 and Tm18, in T. multiceps.

References

1. Scala A, Cancedda GM, Varcasia A, Ligios C, Garippa G, Genchi C. A survey of Taenia multiceps coenurosis in Sardinian sheep. Vet Parasitol. (2007) 143:294–8. doi: 10.1016/j.vetpar.2006.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Afonso SMS, Mukaratirwa S, Hajovska K, Capece BPS, Cristofol C, Arboic M, et al. Prevalence and morphological characteristics of Taenia multiceps cysts (Coenurus cerebralis) fromabattoir-slaughtered and experimentally infected goats. J Neuroparasitol. (2011) 2:1–5. doi: 10.4303/jnp/235532

CrossRef Full Text | Google Scholar

3. Scott PR. Diagnosis and treatment of coenurosis in sheep. Vet Parasitol. (2012) 189:75–8. doi: 10.1016/j.vetpar.2012.03.034

CrossRef Full Text | Google Scholar

4. Edwards GT, Herbert IV. Observations on the course of Taenia multiceps infections in sheep: clinical signs and post-mortem findings. Br Vet J. (1982) 138:489–500. doi: 10.1016/S0007-1935(17)30934-X

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Abo-Shehada MN, Jebreen E, Arab B, Mukbel R, Torgerson PR. Prevalence of Taenia multiceps in sheep in northern Jordan. Prev Vet Med. (2002) 55:201–7. doi: 10.1016/S0167-5877(02)00056-9

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Achenef M, Markos T, Feseha G, Hibret A, Tembely S. Coenurus cerebralis infection in Ethiopian highland sheep: incidence and observations on pathogenesis and clinical signs. Trop Anim Health Prod. (1999) 31:15–24. doi: 10.1023/A:1005125316275

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sharma D, Chauhan P. Coenurosis status in Afro-Asian region: a review. Small Rumin Res. (2006) 64:197–202 doi: 10.1016/j.smallrumres.2005.05.021

CrossRef Full Text | Google Scholar

8. Ing MB, Schantz PM, Turner JA. Human coenurosis in North America: case reports and review. Clin Infect Dis. (1998) 27:519–23. doi: 10.1086/514716

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Collomb J, Machouart M, Biava MF, Brizion M, Montagne K, Plénat F, et al. Contribution of NADH dehydrogenase subunit I and cytochrome C oxidase subunit I sequences toward identifying a case of human coenuriasis in France. J Parasitol. (2007) 93:934–7. doi: 10.1645/GE-1160R.1

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Benifla M, Barrelly R, Shelef I, El-On J, Cohen A, Cagnano E. Huge hemispheric intraparenchymal cyst caused by Taenia multiceps in a child. Case report. J Neurosurg. (2007) 107:511–4. doi: 10.3171/PED-07/12/511

PubMed Abstract | CrossRef Full Text | Google Scholar

11. El-On J, Shelef I, Cagnano E, Benifla M. Taenia multiceps: a rare human cestode infection in Israel. Vet Ital. (2008) 44:621–31.

PubMed Abstract | Google Scholar

12. Ambekar S, Prasad C, Dwarakanath S, Mahadevan A. MRS findings in cerebral coenurosis due to Taenia multiceps. J Neuroimaging. (2013) 23:149–51. doi: 10.1111/j.1552-6569.2011.00616.x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Brehm K, Spiliotis M, Zavala-Góngora R, Konrad C, Frosch M. The molecular mechanisms of larval cestode development: first steps into an unknown world. Parasitol Int. (2006) 55(Suppl):S15–S21. doi: 10.1016/j.parint.2005.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Koziol U. Evolutionary developmental biology (evo-devo) of cestodes. Exp Parasitol. (2017) 180:84–100. doi: 10.1016/j.exppara.2016.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Siles-Lucas M, Hemphill A. Cestode parasites: application of in vivo and in vitro models for studies on the host-parasite relationship. Adv Parasitol. (2002) 51:133–230. doi: 10.1016/S0065-308X(02)51005-8

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li WH, Liu B, Yang Y, Ren YW, Wang S, Liu CH, et al. The genome of tapeworm Taenia multiceps sheds light on understanding parasitic mechanism and control of coenurosis disease. DNA Res. (2018) 25:499–510. doi: 10.1093/dnares/dsy020

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Li WH, Zhang NZ, Yue L, Yang Y, Li L, Yan HB, et al. Transcriptomic analysis of the larva Taenia multiceps. Res Vet Sci. (2017) 115:407–11. doi: 10.1016/j.rvsc.2017.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wu XH, Fu Y, Yang DY, Zhang RH, Zheng WP, Nie HM, et al. Detailed transcriptome description of the neglected cestode Taenia multiceps. PLoS One. (2012) 7:e45830. doi: 10.1371/journal.pone.0045830

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhang S. Comparative transcriptomic analysis of the larval and adult stages of Taenia pisiformis. Genes (Basel). (2019) 10:507. doi: 10.3390/genes10070507

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Basika T, Paludo GP, Araujo FM, Salim AC, Pais F, Maldonado L, et al. Transcriptomic profile of two developmental stages of the cestode parasite Mesocestoides corti. Mol Biochem Parasitol. (2019) 229:35–46. doi: 10.1016/j.molbiopara.2019.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhang L, Cao Z, Liang G, Li X, Wu H, Yang G. Comparative transcriptome analysis reveals candidate genes related to structural and storage carbohydrate biosynthesis in kelp Saccharina japonica (Laminariales, Phaeophyceae). J Phycol. (2020) 56:1168–83. doi: 10.1111/jpy.13016

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gauci C, Vural G, Oncel T, Varcasia A, Damian V, Kyngdon CT, et al. Vaccination with recombinant oncosphere antigens reduces the susceptibility of sheep to infection with Taenia multiceps. Int J Parasitol. (2008) 38:1041–50. doi: 10.1016/j.ijpara.2007.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wang IC, Ma YX, Kuo CH, Fan PC. A comparative study on egg hatching methods and oncosphere viability determination for Taenia solium eggs. Int J Parasitol. (1997) 27:1311–4. doi: 10.1016/S0020-7519(97)00087-8

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Rajasekariah GR, Rickard MD, Mitchell GF. Density-gradient separation of Taenia pisiformis oncospheres. J Parasitol. (1980) 66:355–6. doi: 10.2307/3280839

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Andrews DF, Brasher PM, Manchester KE, Percy ME, Rusk AC, Soltan HC, et al. DUCHEN: an interactive computer program for calculating heterozygosity (carrier) risks in X-linked recessive lethal diseases, and its application in Duchenne muscular dystrophy. Am J Med Genet. (1986) 25:211–8. doi: 10.1002/ajmg.1320250203

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. (2012) 9:357–9. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. (2012) 7:562–78. doi: 10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, Boeckmann B, et al. The Universal Protein Resource (UniProt): an expanding universe of protein information. Nucleic Acids Res. (2006) 34(Database issue):D187–D91. doi: 10.1093/nar/gkj161

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, et al. InterProScan: protein domains identifier. Nucleic Acids Res. (2005) 33(Web Server issue):W116–20. doi: 10.1093/nar/gki442

CrossRef Full Text | Google Scholar

30. Willis JM, Herbert IV. Some factors affecting the eggs of Taenia multiceps: their transmission onto pasture and their viability. Ann Trop Med Parasitol. (1984) 78:236–42. doi: 10.1080/00034983.1984.11811808

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Herbert IV, Edwards GT, Willis JM. Some host factors which influence the epidemiology of Taenia multiceps infections in sheep. Ann Trop Med Parasitol. (1984) 78:243–28. doi: 10.1080/00034983.1984.11811809

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Varcasia A, Pipia AP, Arru D, Pes AM, Tamponi C, Dore F, et al. Morphological and molecular characterization of bovine coenurosis in Sardinia, Italy. Parasitol Res. (2013) 112:2079–82. doi: 10.1007/s00436-012-3257-3

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Podesta RB, Siddiqui AA. Regulation of gene expression in parasitic flatworms. Parasitol Today. (1991) 7:285–7. doi: 10.1016/0169-4758(91)90098-9

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Loos-Frank B. An up-date of Verster's (1969) ‘Taxonomic revision of the genus Taenia Linnaeus' (Cestoda) in table format. Syst Parasitol. (2000) 45:155–83. doi: 10.1023/A:1006219625792

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wiesgigl M, Clos J. Heat shock protein 90 homeostasis controls stage differentiation in Leishmania donovani. Mol Biol Cell. (2001) 12:3307–16. doi: 10.1091/mbc.12.11.3307

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zheng H, Zhang W, Zhang L, Zhang Z, Li J, Lu G, et al. The genome of the hydatid tapeworm Echinococcus granulosus. Nat Genet. (2013) 45:1168–75. doi: 10.1038/ng.2757

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Tsai IJ, Zarowiecki M, Holroyd N, Garciarrubio A, Sánchez-Flores A, Brooks KL, et al. The genomes of four tapeworm species reveal adaptations to parasitism. Nature. (2013) 496:57–63. doi: 10.1038/nature12031

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Nie HM, Xie Y, Fu Y, Yang YD, Gu XB, Wang SX, et al. Cloning and characterization of the fatty acid-binding protein gene from the protoscolex of Taenia multiceps. Parasitol Res. (2013) 112:1833–39. doi: 10.1007/s00436-013-3328-0

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Dang Z, Yagi K, Oku Y, et al. Evaluation of Echinococcus multilocularis tetraspanins as vaccine candidates against primary alveolar echinococcosis. Vaccine. (2009) 27:7339–45. doi: 10.1016/j.vaccine.2009.09.045

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Dang Z, Yagi K, Oku Y, Kouguchi H, Kajino K, Matsumoto J, et al. A pilot study on developing mucosal vaccine against alveolar echinococcosis (AE) using recombinant tetraspanin 3: vaccine efficacy and immunology. PLoS Negl Trop Dis. (2012) 6:e1570. doi: 10.1371/journal.pntd.0001570

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Zhu Y, Ren J, Da'dara A, Harn D, Xu M, Si J, et al. The protective effect of a Schistosoma japonicum Chinese strain 23kDa plasmid DNA vaccine in pigs is enhanced with IL-12. Vaccine. (2004) 23:78–83. doi: 10.1016/j.vaccine.2004.04.031

CrossRef Full Text | Google Scholar

42. Hancock K, Pattabhi S, Whitfield FW, Yushak ML, Lane WS, Garcia HH, et al. Characterization and cloning of T24, a Taenia solium antigen diagnostic for cysticercosis. Mol Biochem Parasitol. (2006) 147:109–17. doi: 10.1016/j.molbiopara.2006.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Tran MH, Pearson MS, Bethony JM, Smyth DJ, Jones MK, Duke M, et al. Tetraspanins on the surface of Schistosoma mansoni are protective antigens against schistosomiasis. Nat Med. (2006) 12:835–40. doi: 10.1038/nm1430

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Huang F, Dang Z, Suzuki Y, Horiuchi T, Yagi K, Kouguchi H, et al. Analysis on gene expression profile in oncospheres and early stage metacestodes from Echinococcus multilocularis. PLoS Negl Trop Dis. (2016) 10:e0004634. doi: 10.1371/journal.pntd.0004634

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Schroeder LL, Pappas PW, Means GE. Trypsin inactivation by intact Hymenolepis diminuta (Cestoda): some characteristics of the inactivated enzyme. J Parasitol. (1981) 67:378–85. doi: 10.2307/3280560

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Cumino AC, Nicolao MC, Loos JA, Denegri G, Elissondo MC. Echinococcus granulosus tegumental enzymes as in vitro markers of pharmacological damage: a biochemical and molecular approach. Parasitol Int. (2012) 61:579–85. doi: 10.1016/j.parint.2012.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Hancock K, Khan A, Williams FB, Yushak ML, Pattabhi S, Noh J, et al. Characterization of the 8-kilodalton antigens of Taenia solium metacestodes and evaluation of their use in an enzyme-linked immunosorbent assay for serodiagnosis. J Clin Microbiol. (2003) 41:2577–86. doi: 10.1128/JCM.41.6.2577-2586.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Mamuti W, Yamasaki H, Sako Y, Nakao M, Xiao N, Nakaya K, et al. Molecular cloning, expression, and serological evaluation of an 8-kilodalton subunit of antigen B from Echinococcus multilocularis. J Clin Microbiol. (2004) 42:1082–8. doi: 10.1128/JCM.42.3.1082-1088.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Massagué J. TGF-β signal transduction. Annu Rev Biochem. (2003) 67:753–91. doi: 10.1146/annurev.biochem.67.1.753

CrossRef Full Text | Google Scholar

50. Willert K, Nusse R. Wnt proteins. Cold Spring Harb Perspect Biol. (2012) 4:a007864. doi: 10.1101/cshperspect.a007864

CrossRef Full Text | Google Scholar

Keywords: comparative transcriptome, developmental stages, Coenurus cerebralis, differentially expressed gene, Taenia multiceps

Citation: Li WH, Yang Y, Zhang NZ, Wang JK, Liu YJ, Li L, Yan HB, Jia WZ and Fu B (2021) Comparative Transcriptome Analyses of the Developmental Stages of Taenia multiceps. Front. Vet. Sci. 8:677045. doi: 10.3389/fvets.2021.677045

Received: 07 March 2021; Accepted: 31 May 2021;
Published: 05 July 2021.

Edited by:

Vito Colella, The University of Melbourne, Australia

Reviewed by:

Andreas Lazaros Chryssafidis, University of the State of Santa Catarina, Brazil
Xiangye Liu, Xuzhou Medical University, China

Copyright © 2021 Li, Yang, Zhang, Wang, Liu, Li, Yan, Jia and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wen-Hui Li, lwh0317@163.com; Baoquan Fu, fubaoquan@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.