The genome of the hydatid tapeworm Echinococcus granulosus

Novel biomarkers for the early prediction of pediatric cystic echinococcosis post-surgical outcomes

OBJECTIVE

This study aims to search for reliable serological biomarkers allowing the early prediction of cystic echinococcosis (CE) post-operative outcomes.

METHODS

We applied immunoprecipitation (IP) of Echinococcus granulosus protoscolex antigens with pediatric CE patients' plasma collected at 1-month and 1-year post-surgery, followed by Liquid Chromatography with tandem mass spectrometry (LC-MS/MS). We compared IP proteomic content from relapsed patients within the first-year post-surgery (RCE) to cases with no relapses until 3 post-operative years (NRCE). Selected proteins were recombinantly synthesized and assessed for their prognostic performance by Enzyme-linked immunosorbent assay (ELISA).

RESULTS

A total of 305 immunoreactive parasitic proteins were identified, 59 of which were significantly more abundant in RCE than NRCE for both time-points. Four proteins showed the most promising characteristics for predicting CE outcomes: cytoplasmic malate dehydrogenase (Eg-cMDH), citrate synthase (Eg-CS), annexin A6 and severin. ELISA-IgG against the four markers were significantly lower at 1-year post-surgery than 1-month in NRCE, in contrast to RCE that displayed either stable or higher levels. The Eg-cMDH and Eg-CS showed the best prognostic performance, with respective probabilities of being "relapse-free" of 83% and 81%, if a decrease of IgG levels occurred between 1-month and 1-year post-surgery.

CONCLUSION

The Eg-cMDH and Eg-CS are promising biomarkers to predict early CE post-surgical outcomes.

Molecular characterization and immune protection of the 3-hydroxyacyl-CoA dehydrogenase gene in Echinococcus granulosus

BACKGROUND

Cystic echinococcosis (CE) is a serious parasitic zoonosis caused by the larvae of the tapeworm Echinococcus granulosus. The development of an effective vaccine is one of the most promising strategies for controlling CE.

METHODS

The E. granulosus 3-hydroxyacyl-CoA dehydrogenase (EgHCDH) gene was cloned and expressed in Escherichia coli. The distribution of EgHCDH in protoscoleces (PSCs) and adult worms was analyzed using immunofluorescence. The transcript levels of EgHCDH in PSCs and adult worms were analyzed using quantitative real-time reverse transcription PCR (RT-qPCR). The immune protective effects of the rEgHCDH were evaluated.

RESULTS

The 924-bp open reading frame sequence of EgHCDH, which encodes a protein of approximately 34 kDa, was obtained. RT-qPCR analysis revealed that EgHCDH was expressed in both the PSCs and adult worms of E. granulosus. Immunofluorescence analysis showed that EgHCDH was mainly localized in the tegument of PSCs and adult worms. Western blot analysis showed that the recombinant protein was recognized by E. granulosus-infected dog sera. Animal challenge experiments demonstrated that dogs immunized with recombinant (r)EgHCDH had significantly higher serum IgG, interferon gamma and interleukin-4 concentrations than the phosphate-buffered saline (PBS) control group. The rEgHCDH vaccine was able to significantly reduce the number of E. granulosus and inhibit the segmental development of E. granulosus compared to the PBS control group.

CONCLUSIONS

The results suggest that rEgHCDH can induce partial immune protection against infection with E. granulosus and could be an effective candidate for the development of new vaccines.

Metabolomic profiling of active and inactive liver cystic echinococcosis

Cystic Echinococcosis (CE) is one of the life-threatening diseases worldwide. It is a parasitic zoonosis caused by tapeworms of the species Echinococcus granulosus sensu lato (s.l). The treatment options of CE vary from simple "watch and wait" approach to invasive treatment, based on the type and especially the nature of the cyst (active/inactive). Serological tests are inadequate to distinguish between active and inactive CE. A diagnostic reference that can determine whether the cyst is active or inactive can easily guide the treatment strategy. We aimed to test whether gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-quadropole time of flight mass spectrometry (LC-qTOF-MS) based metabolomics can establish a plasma metabolic fingerprint of CE patients and identify a diagnostic reference to discriminate active and inactive CE cysts. Metabolite concentrations were measured in plasma samples of 36 active CE patients, 17 inactive CE patients and 31 healthy controls. Multivariate statistical analysis on 232 identified metabolites obtained from two analytical platforms was performed by using principle component analysis (PCA) and partial least square-discriminant analysis (PLS-DA) methods. The PLS-DA scores plot of the combined data set demonstrated a good separation between the groups. Compared to the healthy control group, decreased levels of squalene and increased levels of glyceric acid, 3-phosphoglycerate, glutamic acid, palmitoleic acid and oleic acid were determined in the CE patients. However, decreased levels of 3-phosphoglycerate and increased levels of 4-hydroxyphenylacetylglutamine, docosahexanoic acid were determined in active CE patients compared to the inactive CE patients. Determination of differences in metabolites may provide detailed understandings of potential metabolic process associated with active and inactive CE patients, and altered specific metabolic changes may provide some clues to obtain diagnostic reference for CE. This study has certain limitations: a. various factors affecting results of metabolomic studies such as lifestyle and dietary habits of the patients could not be fully controlled b. other infectious or malignant diseases of the liver should also be included as a positive control to evaluate the specificity of the diagnostic references.

Expression and distribution of glutathione transferases in protoscoleces of Echinococcus granulosus sensu lato

Glutathione transferases (GSTs) belong to a diverse superfamily of multifunctional proteins involved in metabolic detoxification. In helminth parasite, GSTs are particularly relevant since they are also involved in host immunomodulation. Echinococcus granulosus sensu lato (s.l.) is a cestode parasite known to express at least three phylogenetically distant cytosolic GSTs: EgGST1 and EgGST2 previously grouped within Mu and Sigma classes, respectively; and EgGST3 related to both Omega and Sigma classes. To better characterize E. granulosus s.l. GSTs, herein their expression and distribution were assessed in the pre-adult protoscolex (PSC) parasite stage. Potential transcriptional regulatory mechanisms of the corresponding EgGST genes were also explored. Firstly, the transcription of the three EgGSTs was significantly induced during the early stages of the murine model of infection, suggesting a potential role during parasite establishment. EgGST1 was detected in the parenchyma of PSCs and its expression increased after H₂O₂ exposure, supporting its role in detoxification. EgGST2 was mainly detected on the PSCs tegument, strategically localized for potential immunoregulation functions due to its Sigma-class characteristics. In addition, its expression increased after anthelmintic treatment, suggesting a role in chemotherapy resistance. Finally, the Omega-related EgGST3 was localized throughout the entire PSC body, including suckers and tegument, and since its expression also increased after H₂O₂ treatment, a potential role in oxidative stress response could also be ascribed. On the other hand, known cis-acting regulatory motifs were detected in EgGST genes, suggesting similar transcription processes to other eukaryotes. The results herein reported provide additional data regarding the roles of EgGSTs in E. granulosus s.l. biology, contributing to a better understanding of its host-parasite interaction.

Identification and evolution of nuclear receptors in Platyhelminths

Since the first complete set of Platyhelminth nuclear receptors (NRs) from Schistosoma mansoni were identified a decade ago, more flatworm genome data is available to identify their NR complement and to analyze the evolutionary relationship of Platyhelminth NRs. NRs are important transcriptional modulators that regulate development, differentiation and reproduction of animals. In this study, NRs are identified in genome databases of thirty-three species including in all Platyhelminth classes (Rhabditophora, Monogenea, Cestoda and Trematoda). Phylogenetic analysis shows that NRs in Platyhelminths follow two different evolutionary lineages: 1) NRs in a free-living freshwater flatworm (Schmidtea mediterranea) and all parasitic flatworms share the same evolutionary lineage with extensive gene loss. 2) NRs in a free-living intertidal zone flatworm (Macrostomum lignano) follow a different evolutionary lineage with a feature of multiple gene duplication and gene divergence. The DNA binding domain (DBD) is the most conserved region in NRs which contains two C4-type zinc finger motifs. A novel zinc finger motif is identified in parasitic flatworm NRs: the second zinc finger of parasitic Platyhelminth HR96b possesses a CHC2 motif which is not found in NRs of all other animals studied to date. In this study, novel NRs (members of NR subfamily 3 and 6) are identified in flatworms, this result demonstrates that members of all six classical NR subfamilies are present in the Platyhelminth phylum. NR gene duplication, loss and divergence in Platyhelminths are analyzed along with the evolutionary relationship of Platyhelminth NRs.

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.

Transcriptome analysis uncovers the key pathways and candidate genes related to the treatment of Echinococcus granulosus protoscoleces with the repurposed drug pyronaridine

Background

Cystic echinococcosis (CE) is a life-threatening zoonosis caused by the larval form of Echinococcus granulosus tapeworm. Our previous study showed that an approved drug pyronaridine (PND) is highly effective against CE, both in vitro and in an animal model. To identify possible target genes, transcriptome analysis was performed with E. granulosus sensu stricto protoscoleces treated with PND.

Results

A total of 1,321 genes were differentially expressed in protoscoleces treated with PND, including 541 upregulated and 780 downregulated genes. Gene ontology and KEGG analyses revealed that the spliceosome, mitogen-activated protein kinase (MAPK) pathway and ATP-binding cassette (ABC) transporters were the top three enriched pathways. Western blot analysis showed that PND treatment resulted in a dose-dependent increase in protein expression levels of EgMKK1 (MKK3/6-like) and EgMKK2 (MEK1/2-like), two members of MAPK cascades. Interestingly, several heat shock protein (HSP) genes were greatly downregulated including stress-inducible HSPs and their constitutive cognates, and some of them belong to Echinococcus-specific expansion of HSP70.

Conclusions

PND has a great impact on the spliceosome, MAPK pathway and ABC transporters, which may underline the mechanisms by which PND kills E. granulosus protoscoleces. In addition, PND downregulates HSPs expression, suggesting a close relationship between the drug and HSPs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07875-w.

First Evidence of Function for Schistosoma japonicumriok-1 and RIOK-1

Protein kinases are known as key molecules that regulate many biological processes in animals. The right open reading frame protein kinase (riok) genes are known to be essential regulators in model organisms such as the free-living nematode Caenorhabditis elegans. However, very little is known about their function in parasitic trematodes (flukes). In the present study, we characterized the riok-1 gene (Sj-riok-1) and the inferred protein (Sj-RIOK-1) in the parasitic blood fluke, Schistosoma japonicum. We gained a first insight into function of this gene/protein through double-stranded RNA interference (RNAi) and chemical inhibition. RNAi significantly reduced Sj-riok-1 transcription in both female and male worms compared with untreated control worms, and subtle morphological alterations were detected in the ovaries of female worms. Chemical knockdown of Sj-RIOK-1 with toyocamycin (a specific RIOK-1 inhibitor/probe) caused a substantial reduction in worm viability and a major accumulation of mature oocytes in the seminal receptacle (female worms), and of spermatozoa in the sperm vesicle (male worms). These phenotypic alterations indicate that the function of Sj-riok-1 is linked to developmental and/or reproductive processes in S. japonicum.

Impact on Bile Acid Concentrations by Alveolar Echinococcosis and Treatment with Albendazole in Mice

Alveolar echinococcosis (AE) caused by Echinococcus multilocularis is a chronic, progressive liver disease widely distributed in the Northern Hemisphere. The main treatment options include surgical interventions and chemotherapy with benzimidazole albendazole (ABZ). To improve the current diagnosis and therapy of AE, further investigations into parasite-host interactions are needed. This study used liquid chromatography-tandem mass spectrometry (LC-MS/MS) to assess serum and liver tissue bile acid profiles in the i.p. chronic E. multilocularis-infected mouse model and evaluated the effects of the anthelmintic drug ABZ. Additionally, hepatic mRNA and protein expression of enzymes and transporters regulating bile acid concentrations were analyzed. AE significantly decreased unconjugated bile acids in serum and liver tissue. Taurine-conjugated bile salts were unchanged or increased in the serum and unchanged or decreased in the liver. Ratios of unconjugated to taurine-conjugated metabolites are proposed as useful serum markers of AE. The expression of the bile acid synthesis enzymes cytochrome P450 (CYP) 7A1 and aldo-keto reductase (AKR) 1D1 tended to decrease or were decreased in mice with AE, along with decreased expression of the bile acid transporters Na+/taurocholate cotransporting polypeptide (NTCP) and bile salt efflux pump (BSEP). Importantly, treatment with ABZ partially or completely reversed the effects induced by E. multilocularis infection. ABZ itself had no effect on the bile acid profiles and the expression of relevant enzymes and transporters. Further research is needed to uncover the exact mechanism of the AE-induced changes in bile acid homeostasis and to test whether serum bile acids and ratios thereof can serve as biomarkers of AE and for monitoring therapeutic efficacy.

De Novo Sequencing and High-Contiguity Genome Assembly of Moniezia expansa Reveals Its Specific Fatty Acid Metabolism and Reproductive Stem Cell Regulatory Network

Moniezia expansa (M. expansa) parasitizes the small intestine of sheep and causes inhibited growth and development or even death. Being globally distributed, it causes considerable economic losses to the animal husbandry industry. Here, using Illumina, PacBio and BioNano techniques, we obtain a high-quality genome assembly of M. expansa, which has a total length of 142 Mb, a scaffold N50 length of 7.27 Mb and 8,104 coding genes. M. expansa has a very high body fat content and a specific type of fatty acid metabolism. It cannot synthesize any lipids due to the loss of some key genes involved in fatty acid synthesis, and it may can metabolize most lipids via the relatively complete fatty acid β-oxidation pathway. The M. expansa genome encodes multiple lipid transporters and lipid binding proteins that enable the utilization of lipids in the host intestinal fluid. Although many of its systems are degraded (with the loss of homeobox genes), its reproductive system is well developed. PL10, AGO, Nanos and Pumilio compose a reproductive stem cell regulatory network. The results suggest that the high body lipid content of M. expansa provides an energy source supporting the high fecundity of this parasite. Our study provides insight into host interaction, adaptation, nutrient acquisition, strobilization, and reproduction in this parasite and this is also the first genome published in Anoplocephalidae.

Comparative Transcriptome Analyses of the Developmental Stages of Taenia multiceps

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.

Parasitic Helminth-Derived microRNAs and Extracellular Vesicle Cargos as Biomarkers for Helminthic Infections

As an adaption to their complex lifecycles, helminth parasites garner a unique repertoire of genes at different developmental stages with subtle regulatory mechanisms. These parasitic worms release differential components such as microRNAs (miRNAs) and extracellular vesicles (EVs) as mediators which participate in the host-parasite interaction, immune regulation/evasion, and in governing processes associated with host infection. MiRNAs are small (~ 22-nucleotides) non-coding RNAs that regulate gene expression at the post-transcriptional level, and can exist in stable form in bodily fluids such as serum/plasma, urine, saliva and bile. In addition to reports focusing on the identification of miRNAs or in the probing of differentially expressed miRNA profiles in different development stages/sexes or in specific tissues, a number of studies have focused on the detection of helminth-derived miRNAs in the mammalian host circulatory system as diagnostic biomarkers. Extracellular vesicles (EVs), small membrane-surrounded structures secreted by a wide variety of cell types, contain rich cargos that are important in cell-cell communication. EVs have attracted wide attention due to their unique functional relevance in host-parasite interactions and for their potential value in translational applications such as biomarker discovery. In the current review, we discuss the status and potential of helminth parasite-derived circulating miRNAs and EV cargos as novel diagnostic tools.

Genome of the fatal tapeworm Sparganum proliferum uncovers mechanisms for cryptic life cycle and aberrant larval proliferation

The cryptic parasite Sparganum proliferum proliferates in humans and invades tissues and organs. Only scattered cases have been reported, but S. proliferum infection is always fatal. However, S. proliferum’s phylogeny and life cycle remain enigmatic. To investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the mechanisms underlying pathogenicity, we sequenced the entire genomes of S. proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei. Additionally, we performed larvae transcriptome analyses of S. proliferum plerocercoid to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum has experienced a clearly distinct evolutionary history from S. erinaceieuropaei. Moreover, we found that nonordinal extracellular matrix coordination allows asexual reproduction in the host, and loss of sexual maturity in S. proliferum are responsible for its fatal pathogenicity to humans. Our high-quality reference genome sequences should be valuable for future studies of pseudophyllidean tapeworm biology and parasitism.

Kikuchi et al. sequence the genome of the fatal tapeworm Sparganum proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei, and describe its genomic features suggesting the natural history and molecular mechanisms underlying pathogenicity. Their findings indicate that nonordinal extracellular matrix coordination is important for its asexual reproduction, and suggest that loss of sexual maturity contributes to the fatal pathogenicity of S. proliferum to humans.

Screening, construction, and serological identification of the diagnostic antigen molecule EG-06283 for the diagnosis of Echinococcus granulosus

Several strategies exist to prevent and control echinococcosis, a global parasitic disease. However, most treatments are ineffective and adverse effects are common. Therefore, we aimed to screen protoscolex antigen molecules of Echinococcus granulosus to identify a diagnostic biomarker for hydatid disease. Published E. granulosus transcriptome sequencing data were analyzed to screen for antigen molecules that are highly expressed in protoscoleces but not in oncospheres. The membrane protein EG-06283 (annotated as Frizzled-4) was selected from 16 antigens, and its gene fragment was subjected to codon optimization and synthesis. rEG-06283 expression was induced in the pET-24a/EG-06283/BL21 strain; subsequently, the protein was purified and subcutaneously injected into ICR mice at weeks 0, 2, 4, and 6. Blood sampling occurred periodically to quantify serum immunoglobulin G (IgG) levels via enzyme-linked immunosorbent assays (ELISA). Immunogenicity was determined by western blot assays using sera from normal mice and mice with secondary hydatid infections. The antigen's immune reactivity and diagnostic value were validated using sera of patients with hydatid disease. ELISA results confirmed that the antigen molecule induced specific IgG production in mice, resulting in significantly higher levels than those in the adjuvant and control groups (P < 0.05). The western blot results indicated that the protein was recognized by antibodies in the sera of mice with hydatid infection and the antisera of immunized mice. Quantification of protein levels in the sera of patients with hydatid disease significantly differed from levels in healthy participants (P < 0.05). These results indicate that rEG-06283 is a potential diagnostic antigen for E. granulosus infections.

In vitro effects of tropisetron and granisetron against Echinococcus granulosus (s.s.) protoscoleces by involvement of calcineurin and calmodulin

Background

Cystic echinococcosis (CE) is a disease caused by the larval stage of Echinococcus granulosus sensu lato  (s.l.). The treatment of CE mainly relies on the use of benzimidazoles, which can commonly cause adverse side effects. Therefore, more efficient treatment options are needed. Drug repurposing is a useful approach for advancing drug development. We have evaluated the in vitro protoscolicidal effects of tropisetron and granisetron in E. granulosus sensu stricto (s.s.) and assessed the expression of the calcineurin (CaN) and calmodulin (CaM) genes, both of which have been linked to cellular signaling activities and thus are potentially promising targets for the development of drugs.

Methods

Protoscoleces (PSC) of E. granulosus (s.s.) (genotype G1) obtained from sheep hepatic hydatid cysts were exposed to tropisetron and granisetron at concentrations of 50, 150 and 250 µM for various periods of time up to 10 days. Cyclosporine A (CsA) and albendazole sulfoxide were used for comparison. Changes in the morphology of PSC were investigated by light microscopy and scanning electron microscopy. Gene expression was assessed using real-time PCR at the mRNA level for E. granulosus calcineurin subunit A (Eg-CaN-A), calcineurin subunit B (Eg-CaN-B) and calmodulin (Eg-CaM) after a 24-h exposure at 50 and 250 µM, respectively.

Results

At 150 and 250 µM, tropisetron had the highest protoscolicidal effect, whereas CsA was most effective at 50 µM. Granisetron, however, was less effective than tropisetron at all three concentrations. Examination of morphological alterations revealed that the rate at which PSC were killed increased with increasing rate of PSC evagination, as observed in PSC exposed to tropisetron. Gene expression analysis revealed that tropisetron at 50 μM significantly upregulated Eg-CaN-B and Eg-CaM expression while at 250 μM it significantly downregulated both Eg-CaN-B and Eg-CaM expressions; in comparison, granisetron decreased the expression of all three genes at both concentrations.

Conclusions

Tropisetron exhibited a higher efficacy than granisetron against E. granulosus (s.s.) PSC, which is probably due to the different mechanisms of action of the two drugs. The concentration-dependent effect of tropisetron on calcineurin gene expression might reflect its dual functions, which should stimulate future research into its mechanism of action and evaluation of its potential therapeutical effect in the treatment of CE.

Graphic Abstract

Intestinal transcriptomes in Kazakh sheep with different haplotypes after experimental Echinococcus granulosus infection

Cystic echinococcosis (CE) is a chronic zoonosis caused by infection with the larval stage of the cestode Echinococcus granulosus. As the intermediate host, sheep are highly susceptible to this disease. Our previous studies have shown that sheep with haplotype MHC Mva Ibc-Sac IIab-Hin1I ab were resistant to CE infection, while their counterparts without this haplotype were not. In order to reveal the molecular mechanism of resistance in Kazakh sheep, after selecting the differential miRNA in our previous study, herein, transcriptome analyses were conducted to detect the differential expression genes in the intestinal tissue of Kazakh sheep with resistant and non-resistant MHC haplotypes, after peroral infection with E. granulosus eggs. A total of 3835 differentially expressed genes were identified between the two groups, with 2229 upregulated and 1606 downregulated. Further function analysis showed that the most significant genes were related to both innate immune response and adaptive response participating in the defense against E. granulosus infection and the metabolic changes associated with it. The results suggest that genes related to lectin receptors, NK cells activation, chemokines, and tumor necrosis factor, may play important roles in the response of intestinal tissue to E. granulosus.

Dipylidium caninum draft genome - a new resource for comparative genomic and genetic explorations of flatworms

Here, we present a draft genome of the tapeworm Dipylidium caninum (family Dipylidiidae) and compare it with other cestode genomes. This draft genome of D. caninum is 110 Mb in size, has a repeat content of ~13.4% and is predicted to encode ~10,000 protein-coding genes. We inferred excretory/secretory molecules (representing the secretome), other key groups of proteins (including peptidases, kinases, phosphatases, GTPases, receptors, transporters and ion-channels) and predicted potential intervention targets for future evaluation. Using 144 shared single-copy orthologous sequences, we investigated the genetic relationships of cestodes for which nuclear genomes are available. This study provides first insights into the molecular biology of D. caninum and a new resource for comparative genomic and genetic explorations of this and other flatworms.

BRAKER2: automatic eukaryotic genome annotation with GeneMark-EP+ and AUGUSTUS supported by a protein database

The task of eukaryotic genome annotation remains challenging. Only a few genomes could serve as standards of annotation achieved through a tremendous investment of human curation efforts. Still, the correctness of all alternative isoforms, even in the best-annotated genomes, could be a good subject for further investigation. The new BRAKER2 pipeline generates and integrates external protein support into the iterative process of training and gene prediction by GeneMark-EP+ and AUGUSTUS. BRAKER2 continues the line started by BRAKER1 where self-training GeneMark-ET and AUGUSTUS made gene predictions supported by transcriptomic data. Among the challenges addressed by the new pipeline was a generation of reliable hints to protein-coding exon boundaries from likely homologous but evolutionarily distant proteins. In comparison with other pipelines for eukaryotic genome annotation, BRAKER2 is fully automatic. It is favorably compared under equal conditions with other pipelines, e.g. MAKER2, in terms of accuracy and performance. Development of BRAKER2 should facilitate solving the task of harmonization of annotation of protein-coding genes in genomes of different eukaryotic species. However, we fully understand that several more innovations are needed in transcriptomic and proteomic technologies as well as in algorithmic development to reach the goal of highly accurate annotation of eukaryotic genomes.

Molecular characterization and expression analysis of annexin B3 and B38 as secretory proteins in Echinococcus granulosus

Background

Cystic echinococcosis is a parasitic zoonotic disease, which poses a threat to public health and animal husbandry, and causes significant economic losses. Annexins are a family of phospholipid-binding proteins with calcium ion-binding activity, which have many functions.

Methods

Two annexin protein family genes [Echinococcus granulosus annexin B3 (EgAnxB3) and EgAnxB38] were cloned and molecularly characterized using bioinformatic analysis. The immunoreactivity of recombinant EgAnxB3 (rEgAnxB3) and rEgAnxB38 was investigated using western blotting. The distribution of EgAnxB3 and EgAnxB38 in protoscoleces (PSCs), the germinal layer, 18-day strobilated worms and 45-day adult worms was analyzed by immunofluorescence localization, and their secretory characteristics were analyzed preliminarily; in addition, quantitative real-time reverse transcription polymerase chain reaction was used to analyze their transcript levels in PSCs and 28-day strobilated worms stages. The phospholipid-binding activities of rEgAnxB3 and rEgAnxB38 were also analyzed.

Results

EgAnxB3 and EgAnxB38 are conserved and contain calcium-binding sites. Both rEgAnxB3 and rEgAnxB38 could be specifically recognized by the serum samples from E. granulosus-infected sheep, indicating that they had strong immunoreactivity. EgAnxB3 and EgAnxB38 were distributed in all stages of E. granulosus, and their transcript levels were high in the 28-day strobilated worms. They were found in liver tissues near the cysts. In addition, rEgAnxB3 has Ca²⁺-dependent phospholipid-binding properties.

Conclusions

EgAnxB3 and EgAnxB38 contain calcium-binding sites, and rEgAnxB3 has Ca²⁺-dependent phospholipid-binding properties. EgAnxB3 and EgAnxB38 were transcribed in PSCs and 28-day strobilated worms. They were expressed in all stages of E. granulosus, and distributed in the liver tissues near the hydatid cyst, indicating that they are secreted proteins that play a crucial role in the development of E. granulosus.

Molecular characterization and serodiagnostic potential of Echinococcus granulosus hexokinase

Background

Cystic echinococcosis (CE), caused by the larval stage of Echinococcus granulosus (sensu stricto), is a life-threatening but neglected zoonosis. Glycolytic enzymes are crucial molecules for the survival and development of E. granulosus. The aim of this study was to investigate the molecular characterization, immunogenicity, tissue distribution and serodiagnostic potential of E. granulosus hexokinase (EgHK), the first key enzyme in the glycolytic pathway.

Methods

EgHK was cloned and expressed in Escherichia coli. Specific serum antibodies were evaluated in mice immunized with recombinant EgHK (rEgHK). The location of EgHK in the larval stage of E. granulosus was determined using fluorescence immunohistochemistry, and the potential of rEgHK as a diagnostic antigen was investigated in patients with CE using indirect enzyme-linked immunosorbent assay (ELISA).

Results

Recombinant EgHK could be identified in the sera of patients with CE and in mouse anti-rEgHK sera. High titers of specific immunoglobulin G were induced in mice after immunization with rEgHK. EgHK was mainly located in the tegument, suckers and hooklets of protoscoleces and in the germinal layer and laminated layer of the cyst wall. The sensitivity and specificity of the rEgHK-ELISA reached 91.3% (42/46) and 87.8% (43/49), respectively.

Conclusions

We have characterized the sequence, structure and location of EgHK and investigated the immunoreactivity, immunogenicity and serodiagnostic potential of rEgHK. Our results suggest that EgHK may be a promising candidate for the development of vaccines against E. granulosus and an effective antigen for the diagnosis of human CE.

Stage-specific transcriptomic analysis of the model cestode Hymenolepis microstoma

Most parasitic flatworms go through different life stages with important physiological and morphological changes. In this work, we used a transcriptomic approach to analyze the main life-stages of the model tapeworm Hymenolepis microstoma (eggs, cysticercoids, and adults). Our results showed massive transcriptomic changes in this life cycle, including key gene families that contribute substantially to the expression load in each stage. In particular, different members of the cestode-specific hydrophobic ligand-binding protein (HLBP) family are among the most highly expressed genes in each life stage. We also found the transcriptomic signature of major metabolic changes during the transition from cysticercoids to adult worms. Thus, this work contributes to uncovering the gene expression changes that accompany the development of this important cestode model species, and to the best of our knowledge represents the first transcriptomic study with robust replicates spanning all of the main life stages of a tapeworm.

Inhibition of inflammatory cytokine production and proliferation in macrophages by Kunitz-type inhibitors from Echinococcus granulosus

The genus Echinococcus of cestode parasites includes important pathogens of humans and livestock animals. Transcriptomic and genomic studies on E. granulosus and E. multilocularis uncovered striking expansion of monodomain Kunitz proteins. This expansion is accompanied by the specialization of some family members away from the ancestral protease inhibition function to fulfill cation channel blockade functions. Since cation channels are involved in immune processes, we tested the effects on macrophage physiology of two E. granulosus Kunitz-type inhibitors of voltage-activated cation channels (K_v) that are close paralogs. Both inhibitors, EgKU-1 and EgKU-4, inhibited production of the Th1/Th17 cytokine subunit IL-12/23p40 by macrophages stimulated with the TLR4 agonist LPS. In addition, EgKU-4 but not EgKU-1 inhibited production of the inflammatory cytokine IL-6. These activities were not displayed by EgKU-3, a family member that is a protease inhibitor without known activity on cation channels. EgKU-4 potently inhibited macrophage proliferation in response to M-CSF, whereas EgKU-1 displayed similar activity but with much lower potency, similar to EgKU-3. We discuss structural differences, including a heavily cationic C-terminal extension present in EgKU-4 but not in EgKU-1, that may explain the differential activities of the two close paralogs.

Expression profile of circular RNAs in cystic echinococcosis pericystic tissue

Objective

The aim of the present study was to investigate the expression profiles of circular RNAs (circRNAs) in the pericystic tissue of patients with cystic echinococcosis (CE).

Patients and methods

CircRNA expression profiles were obtained by circRNA microarray of four matched pairs of pericystic tissues affected by CE and adjacent normal liver tissues. qRT‐PCR was used to validate the differential expression of some circRNAs identified by the microarray analysis. The potential functions of the differentially expressed circRNAs in the CE pericystic tissues were predicted by bioinformatic analysis.

Results

Compared with the adjacent normal liver tissues, 177 circRNAs were upregulated and 166 circRNAs were downregulated in CE pericystic tissues based on a ≥2.0‐fold change. The top 10 upregulated circRNAs were hsa_circRNA_001654,hsa_circRNA_103361,hsa_circRNA_001490,hsa_circRNA_104310,hsa_circRNA_100395,hsa_circRNA_102485,hsa_circRNA_001459,hsa_circRNA_104193,hsa_circRNA_400043, and hsa_circRNA_006773; The top 10 downregulated circRNAs were hsa_circRNA_400633,hsa_circRNA_404974,hsa_circRNA_068482 ,hsa_circRNA_100974,hsa_circRNA_049637,hsa_circRNA_404798,hsa_circRNA_400064,hsa_circRNA_004045,hsa_circRNA_101379, and hsa_circRNA_016771;The circRNA‐seq results for 15 selected differentially expressed circRNAs were validated by qRT‐PCR. The qRT‐PCR analysis showed that hsa_circRNA_006773, hsa_circRNA_049637, hsa_circRNA_104349, and hsa_circRNA_406281 were differentially expressed in CE pericystic tissues when compared with their expression in the adjacent normal liver tissues. Interestingly, 319 miRNAs and 52 mRNAs were predicted to be adsorbed by these four differentially expressed circRNAs. Gene Ontology analysis revealed that these circRNAs may be involved in the response to organic cyclic compounds and endogenous stimuli and in cellular organismal processes.

Conclusion

Differential expression of circRNAs may be associated with the development and progression of CE, and these circRNAs might be useful as biomarkers for prognosis prediction and as treatment targets.

Guidelines for the use and interpretation of assays for monitoring autophagy (4th edition)1

In 2008, we published the first set of guidelines for standardizing research in autophagy. Since then, this topic has received increasing attention, and many scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Thus, it is important to formulate on a regular basis updated guidelines for monitoring autophagy in different organisms. Despite numerous reviews, there continues to be confusion regarding acceptable methods to evaluate autophagy, especially in multicellular eukaryotes. Here, we present a set of guidelines for investigators to select and interpret methods to examine autophagy and related processes, and for reviewers to provide realistic and reasonable critiques of reports that are focused on these processes. These guidelines are not meant to be a dogmatic set of rules, because the appropriateness of any assay largely depends on the question being asked and the system being used. Moreover, no individual assay is perfect for every situation, calling for the use of multiple techniques to properly monitor autophagy in each experimental setting. Finally, several core components of the autophagy machinery have been implicated in distinct autophagic processes (canonical and noncanonical autophagy), implying that genetic approaches to block autophagy should rely on targeting two or more autophagy-related genes that ideally participate in distinct steps of the pathway. Along similar lines, because multiple proteins involved in autophagy also regulate other cellular pathways including apoptosis, not all of them can be used as a specific marker for bona fide autophagic responses. Here, we critically discuss current methods of assessing autophagy and the information they can, or cannot, provide. Our ultimate goal is to encourage intellectual and technical innovation in the field.

Prediction and Identification of Epitopes in the Emy162 Antigen of Echinococcus multilocularis

INTRODUCTION

Alveolar echinococcosis (AE) is a zoonotic disease caused by the parasitism of Echinococcus multilocularis larvae in the intermediate host or the final host. This study aims to identify and analyze the B-cell and T-cell (Th1, Th2 and Th17) epitopes of E. multilocularis antigen Emy162.

METHODS

(1) The secondary structural characteristics of the Emy162 protein were predicted by bioinformatics software to further predict the potential T- and B-cell epitopes. (2) The dominant antigen epitopes were detected by ELISA through the reaction of patient serum with small B-cell antigen peptide and assessing the proliferation of splenic lymphocytes of mice immunized with Emy162. (3) The expression of cytokines in splenic lymphocytes of mice stimulated by small T-cell antigen peptides was detected by ELISA, ELISpot and flow cytometry to enable the identification of the T-cell epitopes.

RESULTS

(1) The high-scored T-cell epitopes were located at positions E7-13, E36-41, E80-89, E87-96, E97-106 and E129-139, while B-cell epitopes were located at positions E7-13, E19-27, E28-36, E37-48, E78-83, E101-109, E112-121 and E129-139. (2) The three advanced antigen epitopes of Emy162 were E19-27, E112-121 and E129-139. (3) The four Th1 advanced antigen epitopes of Emy162 were E7-13, E36-41, E80-89 and E129-139. The three Th2 advanced antigen epitopes were E36-41, E87-96 and E97-106. The three Th17 advanced antigen epitopes were E36-41, E87-96 and E97-106.

CONCLUSION

(1) The Emy162 protein has advanced antigenicity and numerous potential epitopes. Six T-cell and eight B-cell dominant epitopes were revealed using bioinformatics methods. (2) There are three dominant B-cell epitopes, four dominant Th1 epitopes, three dominant Th2 epitopes, and three dominant Th17 epitopes in the Emy162 antigen.

Metformin Suppresses Development of the Echinococcus multilocularis Larval Stage by Targeting the TOR Pathway

Alveolar echinococcosis (AE) is a severe disease caused by the larval stage of the tapeworm Echinococcus multilocularis Current chemotherapeutic treatment options based on benzimidazoles are of limited effectiveness, which underlines the need to find new antiechinococcosis drugs. Metformin is an antihyperglycemic and antiproliferative agent that shows activity against the related parasite Echinococcus granulosus Hence, we assessed the in vitro and in vivo effects of the drug on E. multilocularis Metformin exerted significant dose-dependent killing effects on in vitro cultured parasite stem cells and protoscoleces and significantly reduced the dedifferentiation of protoscoleces into metacestodes. Likewise, oral administration of metformin (50 mg/kg of body weight/day for 8 weeks) was effective in achieving a significant reduction of parasite weight in a secondary murine AE model. Our results revealed mitochondrial membrane depolarization, activation of Em-AMPK, suppression of Em-TOR, and overexpression of Em-Atg8 in the germinal layer of metformin-treated metacestode vesicles. The opposite effects on the level of active Em-TOR in response to exogenous insulin and rapamycin suggest that Em-TOR is part of the parasite's insulin signaling pathway. Finally, the presence of the key lysosomal pathway components, through which metformin reportedly acts, was confirmed in the parasite by in silico assays. Taken together, these results introduce metformin as a promising candidate for AE treatment. Although our study highlights the importance of those direct mechanisms by which metformin reduces parasite viability, it does not necessarily preclude any additional systemic effects of the drug that might reduce parasite growth in vivo.

Echinococcus granulosus (sensu stricto) (G1, G3) and E. ortleppi (G5) in Pakistan: phylogeny, genetic diversity and population structural analysis based on mitochondrial DNA

Background

Cystic echinococcosis (CE) is a serious tapeworm infection caused by Echinococcus granulosus (sensu lato) which infects a wide range of animals and humans worldwide. Despite the millions of livestock heads reared in Pakistan, only a few reports on CE prevalence and even fewer on the genetic diversity are available for the country. Meanwhile, the available reports on the genetic diversity are predominantly based on short sequences of the cox1 gene.

Methods

To close this knowledge gap, this study was designed to investigate the genetic diversity and population structure of Echinococcus spp. in Pakistan using the complete mitochondrial cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1) genes.

Results

Based on BLAST searches of the generated cox1 and nad1 gene sequences from a total of 60 hydatid cysts collected from cattle (n = 40) and buffalo (n = 20), 52 isolates were identified as E. granulosus (s.s.) (G1, G3) and 8 as E. ortleppi (G5). The detection of the G5 genotype represents the first in Pakistan. The phylogeny inferred by the Bayesian method using nucleotide sequences of cox1-nad1 further confirmed their identity. The diversity indices indicated a high haplotype diversity and a low nucleotide diversity. The negative values of Tajima’s D and Fu’s Fs test demonstrated deviation from neutrality suggesting a recent population expansion.

Conclusions

To the best of our knowledge, this report described the genetic variation of E. granulosus population for the first time in Pakistan using the complete cox1 and nad1 mitochondrial genes and confirms E. ortleppi as one of the causative agents of CE among livestock in Pakistan. While this report will contribute to baseline information for CE control, more studies considering species diversity and distribution in different hosts across unstudied regions of Pakistan are highly needed.

Transcriptomic Analysis of the Early Strobilar Development of Echinococcus granulosus

Echinococcus granulosus has a complex life cycle involving two mammalian hosts. The transition from one host to another is accompanied by changes in gene expression, and the transcriptional events that underlie this transition have not yet been fully characterized. In this study, RNA-seq was used to compare the transcription profiles of samples from E. granulosus protoscoleces induced in vitro to strobilar development at three time points. We identified 818 differentially expressed genes, which were divided into eight expression clusters formed over the entire 24 h period. An enrichment of gene transcripts with molecular functions of signal transduction, enzymes, and protein modifications was observed upon induction and developmental progression. This transcriptomic study provides insights for understanding the complex life cycle of E. granulosus and contributes for searching for the key genes correlating with the strobilar development, which can be used to identify potential candidates for the development of anthelmintic drugs.

Study on the effect of an ion channel inhibitor "Fluralaner" on Echinococcus granulosus protoscolices and metacestode layers in vitro

Hydatid disease has a great impact on public health, causing high morbidity and mortality. Main lines of treatment include surgery, which mostly requires the installation of a scolicidal agent into hydatid cysts to prevent dissemination. Alternatively, medical treatment involves the use of benzimidazole drugs; however, the results are not satisfactory, and new drug compounds are urgently needed. Fluralaner is a potent inhibitor of GABA-gated chloride channels and L-glutamate-gated chloride channels (GluCls) providing immediate and persistent flea, tick and mite control in dogs after a single oral dose. Researches previously identified different genes encoding ion channels in Echinococcus granulosus, making ion channel inhibitors a promising target for treating hydatid disease. Thus, the present study aimed to evaluate the effect of fluralaner on protoscolices and metacestode layers. Parasite materials (Protoscolices, Metacestodes layers) were exposed to different concentrations of the drug ranging from "12.5-100 ug/ml" and examined for viability after 1, 6 and 24 h. Morphological and ultrastructural alterations were recorded by both light and electron microscopies. Immunohistochemical staining confirmed caspase-3 activation as an indicator of apoptosis- induced therapy. The treated protoscolices and metacestode layers showed loss of the viability, the formation of vacuoles and lipid droplets, separation of the germinal layer, and damage in the laminated layer; apoptosis was prominent after treatment. These findings revealed that fluralaner has a potent scolicidal activity and suggested its therapeutic potential against hydatid disease. Further evaluations for animals and human use in the treatment and prevention of hydatid disease are needed.

International consensus on terminology to be used in the field of echinococcoses

Echinococcoses require the involvement of specialists from nearly all disciplines; standardization of the terminology used in the field is thus crucial. To harmonize echinococcosis terminology on sound scientific and linguistic grounds, the World Association of Echinococcosis launched a Formal Consensus process. Under the coordination of a Steering and Writing Group (SWG), a Consultation and Rating Group (CRG) had the main missions of (1) providing input on the list of terms drafted by the SWG, taking into account the available literature and the participants’ experience; and (2) providing independent rating on all debated terms submitted to vote. The mission of the Reading and Review Group (RRG) was to give an opinion about the recommendation paper in terms of readability, acceptability and applicability. The main achievements of this process were: (1) an update of the current nomenclature of Echinococcus spp.; (2) an agreement on three names of diseases due to Echinococcus spp.: Cystic Echinococcosis (CE), Alveolar Echinococcosis (AE) and Neotropical Echinococcosis (NE), and the exclusion of all other names; (3) an agreement on the restricted use of the adjective “hydatid” to refer to the cyst and fluid due to E. granulosus sensu lato; and (4) an agreement on a standardized description of the surgical operations for CE, according to the “Approach, cyst Opening, Resection, and Completeness” (AORC) framework. In addition, 95 “approved” and 60 “rejected” terms were listed. The recommendations provided in this paper will be applicable to scientific publications in English and communication with professionals. They will be used for translation into other languages spoken in endemic countries.

Transcriptomic Features of Echinococcus granulosus Protoscolex during the Encystation Process

Cystic echinococcosis (CE) is a zoonotic infection caused by Echinococcus granulosus larvae. It seriously affects the development of animal husbandry and endangers human health. Due to a poor understanding of the cystic fluid formation pathway, there is currently a lack of innovative methods for the prevention and treatment of CE. In this study, the protoscoleces (PSCs) in the encystation process were analyzed by high-throughput RNA sequencing. A total of 32,401 transcripts and 14,903 cDNAs revealed numbers of new genes and transcripts, stage-specific genes, and differently expressed genes. Genes encoding proteins involved in signaling pathways, such as putative G-protein coupled receptor, tyrosine kinases, and serine/threonine protein kinase, were predominantly up-regulated during the encystation process. Antioxidant enzymes included cytochrome c oxidase, thioredoxin glutathione, and glutathione peroxidase were a high expression level. Intriguingly, KEGG enrichment suggested that differentially up-regulated genes involved in the vasopressin-regulated water reabsorption metabolic pathway may play important roles in the transport of proteins, carbohydrates, and other substances. These results provide valuable information on the mechanism of cystic fluid production during the encystation process, and provide a basis for further studies on the molecular mechanisms of growth and development of PSCs.

Molecular and Functional Characterization of Inhibitor of Apoptosis Proteins (IAP, BIRP) in Echinococcus granulosus

The larval stage of Echinococcus granulosus sensu lato, resulting in cystic echinococcosis, a parasitic zoonosis, causes huge economic losses to the livestock industry and poses a threat to public health. Inhibitor of apoptosis proteins (IAPs) is a class of endogenous anti-apoptotic family, which plays a significant functional role in the regulation of organism’s development. Herein, to explore potential functions of IAPs in E. granulosus, two members of IAPs from E. granulosus (Eg-IAP and Eg-BIRP) were cloned, expressed, and molecularly characterized. Eg-IAP and Eg-BIRP encoded putative 331 and 168 residue proteins, respectively. Bioinformatic analysis showed that both proteins contained a type II BIR domain-the essential functional domain of IAPs. Fluorescence immunohistochemistry revealed that both proteins were ubiquitously localized in all life-cycle stages of E. granulosus. Our fluorescent quantitative PCR (RT-qPCR) results revealed relatively higher transcription levels of two Eg-IAPs in protoscoleces (PSCs) compared to the 18-day strobilated worms. We further used different concentrations of LCL161, a Smac-mimetic pan-IAPs inhibitor, to induce the apoptosis in PSCs in vitro, and revealed that the survival rate of PSCs and transcription levels of both genes were negatively correlated with the concentration of LCL161. While the results of light microscopy, transmission electron microscopy (TEM), and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay also showed a higher apoptotic rate in PSCs with the increasing concentrations of LCL161. Taken together, our findings provide the reasonable evidence that both Eg-IAP and Eg-BIRP have potential implication in critical anti-apoptotic roles during the development of E. granulosus.

Comprehensive characterization of plasma cell-free Echinococcus spp. DNA in echinococcosis patients using ultra-high-throughput sequencing

Background

Echinococcosis is a life-threatening parasitic disease caused by Echinococcus spp. tapeworms with over one million people affected globally at any time. The Echinococcus spp. tapeworms in the human body release DNA to the circulatory system, which can be a biomarker for echinococcosis. Cell-free DNA (cfDNA) is widely used in medical research and has been applied in various clinical settings. As for echinococcosis, several PCR-based tests had been trialed to detect cell-free Echinococcus spp. DNA in plasma or serum, but the sensitivity was about 20% to 25%. Low sensitivity of PCR-based methods might be related to our limited understanding of the features of cell-free Echinococcus spp. DNA in plasma, including its concentration, fragment pattern and release source. In this study, we applied ultra-high-throughput sequencing to comprehensively investigate the characteristics of cell-free Echinococcus spp. DNA in plasma of echinococcosis patients.

Methodology/Principal findings

We collected plasma samples from 23 echinococcosis patients. Total plasma cfDNA was extracted and sequenced with a high-throughput sequencing platform. An average of 282 million read pairs were obtained for each plasma sample. Sequencing data were analyzed with bioinformatics workflow combined with Echinococcus spp. sequence database. After identification of cell-free Echinococcus spp. reads, we found that the cell-free Echinococcus spp. reads accounted for 1.8e-5 to 4.0e-9 of the total clean reads. Comparing fragment length distribution of cfDNA between Echinococcus spp. and humans showed that cell-free Echinococcus spp. DNA of cystic echinococcosis (CE) had a broad length range, while that of alveolar echinococcosis (AE) had an obvious peak at about 135 bp. We found that most of the cell-free Echinococcus spp. DNA reads were from the nuclear genome with an even distribution, which might indicate a random release pattern of cell-free Echinococcus spp. DNA.

Conclusions/Significance

With ultra-high-throughput sequencing technology, we analyzed the concentration, fragment length, release source, and other characteristics of cell-free Echinococcus spp. DNA in the plasma of echinococcosis patients. A better understanding of the characteristics of cell-free Echinococcus spp. DNA in plasma may facilitate their future application as a biomarker for diagnosis.

Echinococcosis is one of the most neglected tropical diseases caused by the metacestodes of Echinococcus spp. tapeworms, which affect both humans and livestock. Plasma cell-free DNA (cfDNA) consists of nucleic acid fragments found extracellularly and may contain DNA released from the parasites. Research shows that a variety of parasites can be detected from plasma cfDNA. Cell-free Echinococcus spp. DNA in plasma or serum had been tested with PCR-based methods, but these PCR methods had low sensitivity ranged from 20% to 25%. Low sensitivity may be due to our limited understanding of cell-free Echinococcus spp. DNA in plasma. Here, we take advantage of high-throughput sequencing to get a comprehensive characterization of cell-free Echinococcus spp. DNA. Our results showed that with high-throughput sequencing we could detect cell-free Echinococcus spp. DNA in all samples, though at a very low level. Based on the sequencing data, we found that cell-free Echinococcus spp. DNA in plasma had a different fragment length distribution to cell-free human DNA, and fragment length distribution of cell-free Echinococcus spp. DNA is also different between cystic echinococcosis (CE) and alveolar echinococcosis (AE). The sequencing data can also help trace the release source of cell-free Echinococcus spp. DNA from the genome. According to the mapping results of cell-free Echinococcus spp. DNA reads, we found that most of them were from the nuclear genome rather than the mitochondrial genome, and their release position showed an even distribution on the genome. These characteristics of cell-free Echinococcus spp. DNA in echinococcosis patients’ plasma could facilitate their future application in research or clinical settings.

Biological and morphological consequences of dsRNA-induced suppression of tetraspanin mRNA in developmental stages of Echinococcus granulosus

BACKGROUND

Cystic echinococcosis, caused by the cestode Echinococcus granulosus, is a neglected tropical disease with remarkable morbidity in humans and a problem of worldwide economic importance in livestock industry. Understanding the molecular basis of the parasite growth and development is essential for the disease diagnosis, management and control. The tetraspanin (TSP) family of proteins are transmembrane proteins with a role in many physiological processes of eukaryotic organisms. TSPs present in the tegumental surface of platyhelminths play pivotal roles in host-parasite interaction. However, little is known about the role of TSPs in growth and development in the Platyhelminthes. To understand the role of TSP1 in the growth and development of E. granulosus we investigated the effect of EgTSP1-specific long dsRNA in different in vitro stages of the parasite.

METHODS

Different stages of E. granulosus, protoscoleces and strobilated worms, were cultivated In vitro in di-phasic media. Using long dsRNA and two delivery methods, i.e. electroporation and electro-soaking, EgTSP1 silencing was performed with an EgTSP1-specific dsRNA. The TSP1 expression profile was assessed as well as the biological and ultrastructural properties of the parasites.

RESULTS

After three days of dsRNA treatment, EgTSP1 expression was significantly reduced in both stages of E. granulosus as compared to irrelevant/unrelated dsRNA and untreated controls. Silencing expression of EgTSP1 in different stages of E. granulosus resulted in reduced viability and body contractions, inhibition of protoscoleces evagination and distinctive tegumental changes. Ultrastructural morphology of the strobilated worms treated with EgTSP1-specific dsRNA was indicative of the microtriches impairments and vacuolated tegument compared to the control helminths.

CONCLUSIONS

Results of the present study suggest that EgTSP1 plays important structural roles in tegument configuration in E. granulosus. EgTSP1 is proved to be a potential target for the development of vaccines and RNAi-based drugs.

Dynamic Changes in the Global Transcriptome and MicroRNAome Reveal Complex miRNA-mRNA Regulation in Early Stages of the Bi-Directional Development of Echinococcus granulosus Protoscoleces

Background

Cystic echinococcosis is a life-threatening disease caused by the larval stages of the dog tapeworm Echinococcus granulosus. Protoscoleces (PSCs) of this worm have the ability of bi-directional development to either larval cysts or strobilar adult worms. However, the molecular mechanisms underlying this development process are unknown.

Results

RNA and small RNAs sequencing was employed to characterize the gene and miRNA expression at 0–24 h and 7–14 days in the bi-directional development of PSCs. A total of 963 genes and 31 miRNAs were differentially expressed in the early development of PSCs to adult worms whereas 972 genes and 27 miRNAs were differentially expressed in the early development of PSCs to cysts. Pairwise comparison between the two developmental patterns showed that 172 genes and 15 miRNAs were differentially expressed at three time-points. Most of these genes were temporally changed at 24 h or 7 days. GO enrichment analysis revealed that the differentially expressed genes in early adult worm development are associated with nervous system development and carbohydrate metabolic process; whereas, the differentially expressed genes in early cystic development are associated with transmembrane transporter activity and nucleoside triphosphatase activity. In addition, miR-71 and miR-219 regulated genes are likely involved in oxidation reduction in adult worm development.

Conclusion

The early stages of bi-directional development in E. granulosus PSCs are controlled by miRNAs and genes likely associated with nervous system development and carbohydrate metabolic process. ATP-dependent transporter genes are associated with cystic development. These results may be important for exploring the mechanisms underlying early development in E. granulosus providing novel information that can be used to discover new therapeutics for controlling cystic echinococcosis.

Old drug repurposing for neglected disease: Pyronaridine as a promising candidate for the treatment of Echinococcus granulosus infections

BACKGROUND

Cystic echinococcosis (CE), a condition caused by the larval stage of the dog tapeworm Echinococcus granulosus sensu stricto, is a globally distributed zoonotic disease. Current treatment options for CE are limited, and an effective and safe anti-echinococcal drug is urgently required.

METHODS

Drug repurposing strategy was employed to identify new therapeutic agents against echinococcal cysts. An in vitro protoscolicidal assay along with in vivo murine models was applied in the drug screening. A microinjection procedure was employed to mimic the clinical PAIR (puncture, aspiration, injection and reaspiration) technique to evaluate the potential application of the candidate drug in clinical practice.

FINDINGS

We repurposed pyronaridine, an approved antimalarial drug, for the treatment of CE. Following a three-dose intraperitoneal regimen (57 mg/kg, q.d. for 3 days), pyronaridine caused 100% cyst mortality. Oral administration of pyronaridine at 57 mg/kg, q.d. for 30 days significantly reduced the parasitic burden in the pre-infected mice compared with albendazole group (p < 0.001). Using a microinjection of drug into cysts, pyronaridine (200 μM) showed highly effective in term of inhibition of cyst growth (p < 0.05, compared with saline group). Pharmacokinetic analysis revealed that pyronaridine was highly distributed in the liver and lungs, the most affected organs of CE. Function analysis showed that pyronaridine inhibited the activity of topoisomerase I (IC50 = 209.7 ± 1.1 μM). In addition, classical apoptotic hallmarks, including DNA fragmentation and caspase activation, were triggered.

INTERPRETATION

Given its approved clinical safety, the repurposing of pyronaridine offers a rapidly translational option for treating CE including PAIR. FUND: National Natural Science Foundation of China and International Cooperation Project of the Qinghai Science and Technology Department.

Revisited - Failure of tetrodotoxin to protect red-spotted newts, Notophthalmus viridescens, from endoparasites

Red-spotted newts, Notophthalmus viridescens, contain tetrodotoxin (TTX) and its analogue 6-epiTTX in variable concentrations. In a follow-up study, newts were sampled from a pond in Pennsylvania, USA, in 2010, 2014, and 2018. Their toxin levels were assayed by liquid-chromatography-fluorescence detection (LC-FLD), and assessment of their infection with endoparasites such as nematodes and helminths was performed by histological examination of internal organs. In the 2010 and 2014 samples, average prevalence of parasite infection was 53 and 60%, respectively, but reached 100% in the 2018 sample, where metacercaria stages of the digenean trematode genus Australapatemon/Apatemon (family: Strigeidae) were predominant causing severe tissue damage in liver and kidney. Mean values of TTX and 6-epiTTX were not significantly different in parasitized or parasite-free newts over the study period, confirming previous findings that host toxicity and parasite load are not negatively correlated. Whereas the role of TTX in defence against predators is undisputed, its efficacy to prevent parasitic infections is less obvious. Toxin-resistance of various metazoan parasites may promote their widespread occurrence in poisonous newts.

In vitro anti-echinococcal activity of octreotide: Additive effect of metformin linked to autophagy

Cystic echinococcosis (CE) is a worldwide zoonosis caused by the Echinococcus granulosus larval stage. The currently available therapy for this disease is based on benzimidazoles, which are rarely curative and cause several adverse effects. Therefore, new treatment options are needed. Octreotide (Oct) is a somatostatin analogue which exhibits anti-proliferative and anti-secretory effects over several cancer cell lines expressing somatostatin receptors. Here, we assessed the in vitro pharmacological effect of Oct against the E. granulosus larval stage. The drug caused a significant dose-dependent decrease in the viability of both protoscoleces and metacestodes. SEM and TEM analysis showed ultrastructural damage in both larval forms under drug treatment. Based on this, we investigated the possible presence of an Oct binding receptor in the parasite. The putative somatostatin/allatostatin-like receptor (Eg-s/ast) conserves the characteristic topology and signature sequences of the prototype somatostatin receptor common to vertebrates and is expressed in both metacestodes and protoscoleces. Moreover, Oct treated-parasites showed the presence of autophagic structures and a significant increase in transcriptional expression of autophagy key genes such as Eg-atg6, Eg-atg8, Eg-atg12 and Eg-atg16. In addition, by in toto immunolocalization assays, an increase in the punctate pattern and Eg-Atg8 protein expression was detected in Oct-treated metacestodes. Subsequently, the combination of Oct and Met had an additive effect on the viability of both larval forms. Our results provide additional evidence for the participation of PI3K/AKT/TOR/autophagy pathway in the Echinococcus survival and suggest the concomitant use of these drugs as potential therapeutic agents in treating of CE.

Review of Cystic Echinococcosis in Nigeria: A Story of Neglect

Purpose

Cystic echinococcosis (CE) caused by Echinococcus granulosus sensu lato is a widespread zoonotic disease of global concern. In Nigeria, the exact picture/status of CE is unclear, as most of the states are largely uninvestigated. Yet, as with every parasitic zoonosis, the first step towards planning a comprehensive management and control programme involves assessment of available national/regional prevalence data, host range, and risk factors at play in the transmission dynamics.

Methods

Published articles on echinococcosis were searched on PubMed and Africa Journal Online (AJOL) databases. Inclusion criteria were based on studies reporting prevalence of echinococcosis in animals and humans (including case reports) from 1970 to 2018.

Results

In this study, we evaluated and summarized cystic echinococcosis reports in Nigeria and found that post 1970–80s, studies on cystic echinococcosis have remained sparse regardless of the high prevalence recorded in the early years of CE investigation. In addition, information on the genetic population structure and the role of wildlife in CE transmission is still lacking.

Conclusions

This study appraises the prevalence and distribution of CE in Nigeria and identified areas where surveillance and control efforts should be focused and intensified.

Elucidating the molecular and developmental biology of parasitic nematodes: Moving to a multiomics paradigm

In the past two decades, significant progress has been made in the sequencing, assembly, annotation and analyses of genomes and transcriptomes of parasitic worms of socioeconomic importance. This progress has somewhat improved our knowledge and understanding of these pathogens at the molecular level. However, compared with the free-living nematode Caenorhabditis elegans, the areas of functional genomics, transcriptomics, proteomics and metabolomics of parasitic nematodes are still in their infancy, and there are major gaps in our knowledge and understanding of the molecular biology of parasitic nematodes. The information on signalling molecules, molecular pathways and microRNAs (miRNAs) that are known to be involved in developmental processes in C. elegans and the availability of some molecular resources (draft genomes, transcriptomes and some proteomes) for selected parasitic nematodes provide a basis to start exploring the developmental biology of parasitic nematodes. Indeed, some studies have identified molecules and pathways that might associate with developmental processes in related, parasitic nematodes, such as Haemonchus contortus (barber's pole worm). However, detailed information is often scant and 'omics resources are limited, preventing a proper integration of 'omic data sets and comprehensive analyses. Moreover, little is known about the functional roles of pheromones, hormones, signalling pathways and post-transcriptional/post-translational regulations in the development of key parasitic nematodes throughout their entire life cycles. Although C. elegans is an excellent model to assist molecular studies of parasitic nematodes, its use is limited when it comes to explorations of processes that are specific to parasitism within host animals. A deep understanding of parasitic nematodes, such as H. contortus, requires substantially enhanced resources and the use of integrative 'omics approaches for analyses. The improved genome and well-established in vitro larval culture system for H. contortus provide unprecedented opportunities for comprehensive studies of the transcriptomes (mRNA and miRNA), proteomes (somatic, excretory/secretory and phosphorylated proteins) and lipidomes (e.g., polar and neutral lipids) of this nematode. Such resources should enable in-depth explorations of its developmental biology at a level, not previously possible. The main aims of this review are (i) to provide a background on the development of nematodes, with a particular emphasis on the molecular aspects involved in the dauer formation and exit in C. elegans; (ii) to critically appraise the current state of knowledge of the developmental biology of parasitic nematodes and identify key knowledge gaps; (iii) to cover salient aspects of H. contortus, with a focus on the recent advances in genomics, transcriptomics, proteomics and lipidomics as well as in vitro culturing systems; (iv) to review recent advances in our knowledge and understanding of the molecular and developmental biology of H. contortus using an integrative multiomics approach, and discuss the implications of this approach for detailed explorations of signalling molecules, molecular processes and pathways likely associated with nematode development, adaptation and parasitism, and for the identification of novel intervention targets against these pathogens. Clearly, the multiomics approach established recently is readily applicable to exploring a wide range of interesting and socioeconomically significant parasitic worms (including also trematodes and cestodes) at the molecular level, and to elucidate host-parasite interactions and disease processes.

miRNAs and lncRNAs in Echinococcus and Echinococcosis

Echinococcosis are considered to be potentially lethal zoonotic diseases that cause serious damage to hosts. The metacestode of Echinococcus multilocularis and E. granulosus can result in causing the alveolar and cystic echinococcoses, respectively. Recent studies have shown that non-coding RNAs are widely expressed in Echinococcus spp. and hosts. In this review, the two main types of non-coding RNAs—long non-coding RNAs (lncRNAs) and microRNAs (miRNAs)—and the wide-scale involvement of these molecules in these parasites and their hosts were discussed. The expression pattern of miRNAs in Echinococcus spp. is species- and developmental stage-specific. Furthermore, common miRNAs were detected in three Echinococcus spp. and their intermediate hosts. Here, we primarily focus on recent insights from transcriptome studies, the expression patterns of miRNAs and lncRNAs, and miRNA-related databases and techniques that are used to investigate miRNAs in Echinococcus and echinococcosis. This review provides new avenues for screening therapeutic and diagnostic markers.

The local immune response during Echinococcus granulosus growth in a quantitative hepatic experimental model

The local immune mechanisms responsible for the establishment and development of Echinococcus granulosus sensu stricto infection in the liver, have been little explored. We developed a suitable experimental model that mimics naturally infected livers using portal injection of protoscoleces. Opposite to Echinococcus multilocularis infection which is dose-dependent, fully mature hydatid cysts can be established in the liver whatever the injection dose; although most of the infection sites were seen at the establishment phase as inflammatory granulomas associated with fibrosis, they never matured into cysts. At the establishment phase, a strong immune response was composed of T and B cells, with T1-type, T2-type cells and cytokines and IL-10-secreting CD8⁺ T cells in the liver. At the established phase, results suggested a local production of antibodies by B cells, and an involvement of NK and NKT cells. Infection outcome and local immune response in the liver, were different in the mouse models of Echinococcus granulosus sensu stricto and Echinococcus multilocularis respectively; however, only early specificities at the microenvironment level might explain the major differences found between the lesions induced by the two species. Our quantitative experimental model appears fully appropriate to further study this microenvironment and its relationship with each cestode species.

In vitro metabolomic footprint of the Echinococcus multilocularis metacestode

Alveolar echinococcosis (AE) is a zoonotic disease that is deadly if left untreated. AE is caused by the larval metacestode stage of the cestode Echinococcus multilocularis. Better knowledge on the host-parasite interface could yield novel targets for improvement of the treatment against AE. We analyzed culture media incubated with in vitro grown E. multilocularis metacestodes by 1H nuclear magnetic resonance spectroscopy to identify the unknown metabolic footprint of the parasite. Moreover, we quantitatively analyzed all amino acids, acetate, glucose, lactate, and succinate in time-course experiments using liquid chromatography and enzymatic assays. The E. multilocularis metacestodes consumed glucose and, surprisingly, threonine and produced succinate, acetate, and alanine as major fermentation products. The metabolic composition of vesicle fluid (VF) from in vitro grown E. multilocularis metacestodes was different from parasite-incubated culture medium with respect to the abundance, but not the spectrum, of metabolites, and some metabolites, in particular amino acids, accumulated in the VF. Overall, this study presents the first characterization of the in vitro metabolic footprint of E. multilocularis metacestodes and VF composition, and it provides the basis for analyses of potentially targetable pathways for future drug development.

Hybrid de novo whole-genome assembly and annotation of the model tapeworm Hymenolepis diminuta

Despite the use of Hymenolepis diminuta as a model organism in experimental parasitology, a full genome description has not yet been published. Here we present a hybrid de novo genome assembly based on complementary sequencing technologies and methods. The combination of Illumina paired-end, Illumina mate-pair and Oxford Nanopore Technology reads greatly improved the assembly of the H. diminuta genome. Our results indicate that the hybrid sequencing approach is the method of choice for obtaining high-quality data. The final genome assembly is 177 Mbp with contig N50 size of 75 kbp and a scaffold N50 size of 2.3 Mbp. We obtained one of the most complete cestode genome assemblies and annotated 15,169 potential protein-coding genes. The obtained data may help explain cestode gene function and better clarify the evolution of its gene families, and thus the adaptive features evolved during millennia of co-evolution with their hosts.