Transcriptomic Features of Echinococcus granulosus Protoscolex during the Encystation Process

Genome-wide transcriptome analysis of Echinococcus multilocularis larvae and germinative cell cultures reveals genes involved in parasite stem cell function

The lethal zoonosis alveolar echinococcosis is caused by tumour-like growth of the metacestode stage of the tapeworm Echinococcus multilocularis within host organs. We previously demonstrated that metacestode proliferation is exclusively driven by somatic stem cells (germinative cells), which are the only mitotically active parasite cells that give rise to all differentiated cell types. The Echinococcus gene repertoire required for germinative cell maintenance and differentiation has not been characterised so far. We herein carried out Illumina sequencing on cDNA from Echinococcus metacestode vesicles, from metacestode tissue depleted of germinative cells, and from Echinococcus primary cell cultures. We identified a set of ~1,180 genes associated with germinative cells, which contained numerous known stem cell markers alongside genes involved in replication, cell cycle regulation, mitosis, meiosis, epigenetic modification, and nucleotide metabolism. Interestingly, we also identified 44 stem cell associated transcription factors that are likely involved in regulating germinative cell differentiation and/or pluripotency. By in situ hybridization and pulse-chase experiments, we also found a new general Echinococcus stem cell marker, EmCIP2Ah, and we provide evidence implying the presence of a slow cycling stem cell sub-population expressing the extracellular matrix factor Emkal1. RNA-Seq analyses on primary cell cultures revealed that metacestode-derived Echinococcus stem cells display an expanded differentiation capability and do not only form differentiated cell types of the metacestode, but also cells expressing genes specific for protoscoleces, adult worms, and oncospheres, including an ortholog of the schistosome praziquantel target, EmTRPMPZQ. Finally, we show that primary cell cultures contain a cell population expressing an ortholog of the tumour necrosis factor α receptor family and that mammalian TNFα accelerates the development of metacestode vesicles from germinative cells. Taken together, our analyses provide a robust and comprehensive characterization of the Echinococcus germinative cell transcriptome, demonstrate expanded differentiation capability of metacestode derived stem cells, and underscore the potential of primary germinative cell cultures to investigate developmental processes of the parasite. These data are relevant for studies into the role of Echinococcus stem cells in parasite development and will facilitate the design of anti-parasitic drugs that specifically act on the parasite germinative cell compartment.

In vitro and in silico scolicidal effect of sanguinarine on the hydatid cyst protoscoleces

We aimed to investigate the scolicidal effects of sanguinarine on hydatid cyst protoscoleces (PSCs) in vitro and in silico. Different targets were docked into the active sites of sanguinarine. Molecular docking processes and visualization of interactions were performed using AutoDock Vina and Discovery Studio Visualizer. Binding energy was calculated and compared (kcal/mol). PSCs were aspirated from the hydatid cysts and washed. The sediments of PSCs were then exposed to various concentrations (50, 25, 12, 6, 3, and 1 μg/mL) of sanguinarine. The viability test was finally evaluated by the Trypan blue solution 4%. Levels of malondialdehyde (MDA), superoxide dismutase (SOD), glutathione (GSH), glutathione peroxidase (GPX), and catalase were analyzed to assess the level of oxidative stress-treated PSCs. Caspase-3 activity rate was determined to evaluate cell apoptosis in treated PSCs. Among the receptors, acetylcholinesterase was identified as the excellent target, with Vina score of -11.8. Sanguinarine showed high scolicidal effects after 12, 24, and 48 h. Also, in the first hour of exposure to the drug, caspase-3 activity and MDA level significantly increased, but the levels of GSH and GPx had a significant reduction after 12, 24, and 48 h (P < 0.05). The findings of this study revealed that sanguinarine have potent scolicidal effects in vitro and in silico and could be considered an opportunity for the introduction of a novel and safe therapeutic agent for the treatment of cystic echinococcosis. However, supplementary studies will be desired to prove the current findings by examining sanguinarine in a clinical setting.

Bioinformatics analysis and experimental verification of Notch signalling pathway-related miRNA-mRNA subnetwork in extracellular vesicles during Echinococcus granulosus encystation

Background

Encystation of the protoscoleces (PSCs) of Echinococcus granulosus is the main cause of secondary hydatid dissemination in the intermediate host. Extracellular vesicles (EVs) can transfer miRNAs into parasite cells to regulate mRNA expression. However, loading of developmental pathway-related miRNAs, such as those related to the Notch signalling pathway in EVs is unclear. Thus, we screened the miRNA-mRNA subnetwork involved in the Notch pathway during E. granulosus encystation in vitro and assessed changes in expression in the parasite and EVs.

Methods

mRNAs and miRNAs differentially expressed (DE) between PSCs and microcysts (MCs) were screened using high-throughput sequencing. DE mRNAs obtained from transcriptome analysis were intersected with mRNAs predicted to be targets of the conserved DE miRNAs of a small RNA library. DE miRNA functions were analysed using public databases, and a miRNA–mRNA subnetwork related to the Notch pathway was established. Notch pathway-related mRNA and miRNA expression of worms and EVs at different times was verified.

Results

In total, 1445 DE mRNAs between MCs and PSCs were screened after the intersection between 1586 DE mRNAs from the transcriptome and 9439 target mRNAs predicted using 39 DE miRNAs from the small RNA library. The DE mRNAs were clustered into 94 metabolic pathways, including the Notch pathway. Five DE miRNAs, including the most significantly expressed new DE miRNA, egr-new-mir0694-3p, corresponding to four target mRNAs (EgrG_000892700, EgrG_001029400, EgrG_001081400 and EgrG_000465800) were all enriched in the Notch pathway. The expression of the above mRNAs and miRNAs was consistent with the results of high-throughput sequencing, and the expression of each miRNA in EVs was verified. Annotated as ADAM17/TACE in the Notch pathway, EgrG_000892700 was down-regulated during PSC encystation. egr-miR-4989-3p and egr-miR-277a-3p expression in EVs after encystation was nearly five times that in EVs before encystation, which might regulate the expression of EgrG_000892700.

Conclusions

Five miRNAs corresponding to four target mRNAs may be involved in regulating the Notch pathway during the PSC encystation. EVs may regulate the expression of EgrG_000892700 in PSCs because of continuous targeting of egr-miR-4989-3p and egr-miR-277a-3p and participate in the regulation the Notch pathway. The study might expand new ideas for blocking the secondary infection of E. granulosus PSCs via EVs miRNAs.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05391-8.

Evolutionary Adaptations of Parasitic Flatworms to Different Oxygen Tensions

During the evolution of the Earth, the increase in the atmospheric concentration of oxygen gave rise to the development of organisms with aerobic metabolism, which utilized this molecule as the ultimate electron acceptor, whereas other organisms maintained an anaerobic metabolism. Platyhelminthes exhibit both aerobic and anaerobic metabolism depending on the availability of oxygen in their environment and/or due to differential oxygen tensions during certain stages of their life cycle. As these organisms do not have a circulatory system, gas exchange occurs by the passive diffusion through their body wall. Consequently, the flatworms developed several adaptations related to the oxygen gradient that is established between the aerobic tegument and the cellular parenchyma that is mostly anaerobic. Because of the aerobic metabolism, hydrogen peroxide (H₂O₂) is produced in abundance. Catalase usually scavenges H₂O₂ in mammals; however, this enzyme is absent in parasitic platyhelminths. Thus, the architecture of the antioxidant systems is different, depending primarily on the superoxide dismutase, glutathione peroxidase, and peroxiredoxin enzymes represented mainly in the tegument. Here, we discuss the adaptations that parasitic flatworms have developed to be able to transit from the different metabolic conditions to those they are exposed to during their life cycle.

The Genomes of Two Strains of Taenia crassiceps the Animal Model for the Study of Human Cysticercosis

Human cysticercosis by Taenia solium is the major cause of neurological illness in countries of Africa, Southeast Asia, and the Americas. Publication of four cestode genomes (T. solium, Echinococcus multilocularis, E. granulosus and Hymenolepis microstoma) in the last decade, marked the advent of novel approaches on the study of the host-parasite molecular crosstalk for cestode parasites of importance for human and animal health. Taenia crassiceps is another cestode parasite, closely related to T. solium, which has been used in numerous studies as an animal model for human cysticercosis. Therefore, characterization of the T. crassiceps genome will also contribute to the understanding of the human infection. Here, we report the genome of T. crassiceps WFU strain, reconstructed to a noncontiguous finished resolution and performed a genomic and differential expression comparison analysis against ORF strain. Both strain genomes were sequenced using Oxford Nanopore (MinION) and Illumina technologies, achieving high quality assemblies of about 107 Mb for both strains. Dotplot comparison between WFU and ORF demonstrated that both genomes were extremely similar. Additionally, karyotyping results for both strains failed to demonstrate a difference in chromosome composition. Therefore, our results strongly support the concept that the absence of scolex in the ORF strain of T. crassiceps was not the result of a chromosomal loss as proposed elsewhere. Instead, it appears to be the result of subtle and extensive differences in the regulation of gene expression. Analysis of variants between the two strains identified 2,487 sites with changes distributed in 31 of 65 scaffolds. The differential expression analysis revealed that genes related to development and morphogenesis in the ORF strain might be involved in the lack of scolex formation.

Chromosome-scale Echinococcus granulosus (genotype G1) genome reveals the Eg95 gene family and conservation of the EG95-vaccine molecule

Cystic echinococcosis is a socioeconomically important parasitic disease caused by the larval stage of the canid tapeworm Echinococcus granulosus, afflicting millions of humans and animals worldwide. The development of a vaccine (called EG95) has been the most notable translational advance in the fight against this disease in animals. However, almost nothing is known about the genomic organisation/location of the family of genes encoding EG95 and related molecules, the extent of their conservation or their functions. The lack of a complete reference genome for E. granulosus genotype G1 has been a major obstacle to addressing these areas. Here, we assembled a chromosomal-scale genome for this genotype by scaffolding to a high quality genome for the congener E. multilocularis, localised Eg95 gene family members in this genome, and evaluated the conservation of the EG95 vaccine molecule. These results have marked implications for future explorations of aspects such as developmentally-regulated gene transcription/expression (using replicate samples) for all E. granulosus stages; structural and functional roles of non-coding genome regions; molecular ‘cross-talk’ between oncosphere and the immune system; and defining the precise function(s) of EG95. Applied aspects should include developing improved tools for the diagnosis and chemotherapy of cystic echinococcosis of humans.

A high-quality genome for the parasitic tapeworm, Echinococcus granulosus, provides further insight into the EG95 vaccine target for cystic echinococcosis.

Transcriptome analysis of Echinococcus granulosus sensu stricto protoscoleces reveals differences in immune modulation gene expression between cysts found in cattle and sheep

Cystic Echinococcosis (CE), a zoonotic parasitic disease, is caused by the cestode Echinococcus granulosus sensu lato. CE inflicts severe damage in cattle, sheep, and human hosts worldwide. Fertile CE cysts are characterized by the presence of viable protoscoleces. These parasite forms are studied with minimal contamination with host molecules. Hosts, cattle and sheep, show differences in their CE cyst fertility. The effect of the host in protoscolex transcriptome is not known. We genotyped and performed transcriptomic analysis on sheep protoscoleces obtained from liver and lung CE cysts. The transcriptomic data of Echinococcus granulosus sensu stricto protoscoleces from 6 lung CE cysts and 6 liver CE cysts were Collected. For host comparison analysis, 4 raw data files belonging to Echinococcus granulosus sensu stricto protoscoleces from cattle liver CE cysts were obtained from the NCBI SRA database. Principal component and differential expression analysis did not reveal any statistical differences between protoscoleces obtained from liver or lung cysts, either within the same sheep or different sheep hosts. Conversely, there are significant differences between cattle and sheep protoscolex samples. We found differential expression of immune-related genes. In cattle, 7 genes were upregulated in protoscoleces from liver cysts. In sheep, 3 genes were upregulated in protoscoleces from liver and lung CE cysts. Noteworthy, are the differential expression of antigen B, tegument antigen, and arginase-2 in samples obtained from sheep CE cysts, and basigin in samples from cattle CE cysts. These findings suggest that the host species is an important factor involved in the differential expression of immune related genes, which in turn is possibly related to the fertility of Echinococcus granulosus sensu stricto cysts.

Synergism therapeutic and immunoregulatory effects of Albendazole + rAd-mIL-28B against Echinococcosis in experiment-infected mice with protoscoleces

The metacestode stage of Echinococcus granulosus can cause cystic echinococcosis (CE), which still widely occurs around the world. Since the early 1970s, benzimidazoles have been shown to inhibit the growth of cysts and used to treat CE. However, benzimidazoles are still ineffective in 20%-40% of cases. In order to explore the new agents against CE, we have investigated the therapeutic effect of the recombinant adenoviral vector expressing mouse IL-28B (rAd-mIL-28B) on protoscoleces-infected mice. In our study, we successfully established the model mice which infected with protoscoleces intraperitoneally. At 18 weeks post-infection, the mice received rAd-mIL-28B (1×10⁷ PFU) weekly by intramuscular injection for 6 weeks. Compared with the untreated control (13.1 ± 2.2 g), there was a significant reduction in cysts wet weight in rAd-mIL-28B group (8.3 ± 3.5 g) (P < 0.05), especially in Albendazole (ABZ) + rAd-mIL-28B group (5.8 ± 1.4 g) (P < 0.01). We also observed the severe damage of the germinal layer and the laminated layer of cysts after treatment. rAd-mIL-28B group showed a prominent increase in the level of Th1 type cytokines (such as IFN-γ, IL-2 and TNF-α). Meanwhile, the frequency of Foxp3⁺ T cells was decreased in the rAd-mIL-28B group (4.83 ± 0.81%) and ABZ + rAd-mIL-28B group (4.60 ± 0.51%), comparing with the untreated group (8.13 ± 2.60%) (P < 0.05). In addition, compared with the untreated control (122.14 ± 81.09 pg/ml), the level of IFN-γ significantly increased in peritoneal fluid in the rAd-mIL-28B group (628.87 ± 467.16 pg/ml) (P < 0.05) and ABZ + rAd-mIL-28B group (999.76 ± 587.60 pg/ml) (P < 0.001). Taken together, it suggested that ABZ + IL-28B may be a potential therapeutic agent against CE.

Echinococcosis is a chronic zoonotic parasitic disease, which is caused by the larval stage of Echinococcus granulosus (E. granulosus) and Echinococcus multiocularis (E. multiocularis). This disease is still widely prevalent in the world and seriously endangers human health and life, causing heavy burdens and economic losses to agriculture and animal husbandry. China is also one of the high incidence areas of the disease. At present, the preferred treatment is surgical excision of the parasitic mass, but patients with multiple cysts appear in multiple organs have to receive drug therapy. The approved chemotherapeutic drugs in clinic, such as albendazole (ABZ) and mebendazole, often do not work expectedly with a high rate of recurrence. Therefore, it is urgent to develop the new anti-echinococcal drug. IL-28B is an important member of type-III IFNs, which is equipped with the capacity of anti-viral activity and anti-tumor. In our previous study, we found that IL-28B could inhibit the proliferation of cervical cancer cells via down-regulating Treg cells in mice. This study mainly studied the therapeutic effect of rAd-mIL-28B on E. granulosus-infected mice, and the results showed that rAd-mIL-28B could relieve the parasitic burden and inhibit Treg cells meanwhile improve the Th1 and Th17 immune responses. It may contribute to another choice for the anti-echinococcal treatment.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BACKGROUND

The complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.

METHODS

Whole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.

RESULTS

At least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.

CONCLUSIONS

The high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.

iTRAQ-based comparative proteomic analysis in different developmental stages of Echinococcus granulosus

Cystic echinococcosis, caused by infection with the larval stage of the cestode Echinococcus granulosus, is a chronic zoonosis. The lifecycle of the E. granulosus parasite includes three consecutive stages that require specific gene regulation or protein expression to survive environmental shifts between definitive hosts and intermediate hosts. The aim of the present study is to screen and analyze the stage differential antigens to be considered for vaccine development against E. granulosus. By using the iTRAQ (isobaric tags for relative and absolute quantification) method, the differentially expressed proteins were selected from the three consecutive developmental stages of E. granulosus: oncosphere, adult tapeworms, and protoscolex. Through a bioinformatics analysis including Clusters of Orthologous Groups (COG), Gene Ontology (GO), and pathway metabolic annotation, we identified some proteins of interest from each stage. The results showed that a large number of differentially expressed proteins (375: oncosphere vs. adult, 346: oncosphere vs. protoscolex, and 391: adult vs. protoscolex) were identified from the three main lifecycle stages. Analysis of the differential protein pathways showed that these differential proteins are mainly enriched in metabolic pathways, Huntington’s diseases, Alzheimer’s diseases, and ribosome metabolic pathways. Interestingly, among these differential proteins, expression levels of paramyosin, HSP60, HSP70, HSP90, cathepsin L1, cathepsin D, casein kinase, and calmodulin were significantly higher in the oncosphere than in the adult or protoscolex (p < 0.05). We hope our findings will help to identify potential targets for diagnosis or for therapeutic and prophylactic intervention.