The Echinococcus canadensis (G7) genome: a key knowledge of parasitic platyhelminth human diseases

Molecular characterization of EcCLP1, a new putative cathepsin L protease from Echinococcus canadensis

Echinococcus granulosus sensu lato is a platyhelminth parasite and the etiological cause of cystic echinococcosis (CE), a zoonotic and neglected disease that infects animals and humans worldwide. As a part of the biological arsenal of the parasite, cathepsin L proteases are a group of proteins that are believed to be essential for parasite penetration, immune evasion, and establishment in the tissues of the host. In this work, we have cloned and sequenced a new putative cathepsin L protease from Echinococcus canadensis (EcCLP1). The bioinformatic analysis suggests that EcCLP1 could be synthesized as a zymogen and activated after proteolytic cleavage. The multiple sequence alignment with other cathepsin proteases reveals important functional conserved features like a conserved active site, an N-linked glycosylation residue, a catalytic triad, an oxyanion hole, and three putative disulfide bonds. The phylogenetic analysis suggests that EcCLP1 could indeed be a cathepsin L cysteine protease from clade 1 as it grouped with cathepsins from other species in this clade. Modeling studies suggest that EcCLP1 has two domains forming a cleft where the active site is located and an occluding role for the propeptide. The transcriptomic analysis reveals different levels of cathepsin transcript expression along the different stages of the parasite life cycle. The whole-mount immunohistochemistry shows an interesting superficial punctate pattern of staining which suggests a secretory pattern of expression. The putative cathepsin L protease characterized here may represent an interesting tool for diagnostic purposes, vaccine design, or a new pharmacological target for antiparasitic intervention.

Countering Zoonotic Diseases: Current Scenario and Advances in Diagnostics, Monitoring, Prophylaxis and Therapeutic Strategies

Human life and health have interacted reciprocally with the surrounding environment and animal fauna for ages. This relationship is evident in developing nations, where human life depends more on the animal population for food, transportation, clothing, draft power, and fuel sources, among others. This inseparable link is a potent source of public health issues, especially in outbreaks of zoonotic diseases transmitted from animals to humans. Zoonotic diseases are referred to as diseases that are naturally transmitted between vertebrate animals and humans. Among the globally emerging diseases in the last decade, 75% are of animal origin, most of which are life-threatening. Since most of them are caused by potent new pathogens capable of long-distance transmission, the impact is widespread and has serious public health and economic consequences. Various other factors also contribute to the transmission, spread, and outbreak of zoonotic diseases, among which industrialization-led globalization followed by ecological disruption and climate change play a critical role. In this regard, all the possible strategies, including advances in rapid and confirmatory disease diagnosis and surveillance/monitoring, immunization/vaccination, therapeutic approaches, appropriate prevention and control measures to be adapted, and awareness programs, need to be adopted collaboratively among different health sectors in medical, veterinary, and concerned departments to implement the necessary interventions for the effective restriction, minimization, and timely control of zoonotic threats. The present review focuses on the current scenario of zoonotic diseases and their counteracting approaches to safeguard their health impact on humans.

Human cystic echinococcosis: first molecular identification of Echinococcus canadensis G7 in Brazil

Echinococcus granulosus sensu lato (s.l.) is a species complex with the potential to cause cystic echinococcosis (CE). Contact with the feces of domestic dogs (Canis familiaris) fed with raw viscera of intermediate livestock hosts is a risk factor for this infection in the southern region of Brazil. Although the region has been considered endemic to CE for many years, molecular data regarding the species of the complex causing CE in humans are scarce. This study aimed to perform a molecular analysis of the biological fluid from a human liver cyst to investigate the species responsible for CE. Genetic material obtained from the hydatid fluid of a hepatic cyst from a human with CE was subjected to PCR to amplify mitochondrial and nuclear DNA sequences. The phylogenetic analysis confirmed the human infection by Echinococcus canadensis G7 in the state of Paraná, Brazil. This is the first molecular record of E. canadensis G7 infecting a human in Brazil, and it is important to reiterate the risk of human CE caused by this species in South America, as reported by a previous study in Patagonia, Argentina. From the epidemiological point of view, this finding is of great relevance for the southern region of Brazil, since this parasite has previously only been detected in pigs in the state of Rio Grande do Sul, neighboring Paraná. The finding points to the importance of this identification in the molecular epidemiology of E. granulosus s.l., especially in South America.

A Comparative Analysis of the Protein Cargo of Extracellular Vesicles from Helminth Parasites

Helminth parasites cause debilitating-sometimes fatal-diseases in humans and animals. Despite their impact on global health, mechanisms underlying host-parasite interactions are still poorly understood. One such mechanism involves the exchange of extracellular vesicles (EVs), which are membrane-enclosed subcellular nanoparticles. To date, EV secretion has been studied in helminth parasites, including EV protein content. However, information is highly heterogeneous, since it was generated in multiple species, using varied protocols for EV isolation and data analysis. Here, we compared the protein cargo of helminth EVs to identify common markers for each taxon. For this, we integrated published proteomic data and performed a comparative analysis through an orthology approach. Overall, only three proteins were common in the EVs of the seven analyzed species. Additionally, varied repertoires of proteins with moonlighting activity, vaccine antigens, canonical and non-canonical proteins related to EV biogenesis, taxon-specific proteins of unknown function and RNA-binding proteins were observed in platyhelminth and nematode EVs. Despite the lack of consensus on EV isolation protocols and protein annotation, several proteins were shown to be consistently detected in EV preparations from organisms at different taxa levels, providing a starting point for a selective biochemical characterization.

Circulating Small RNA Profiling of Patients with Alveolar and Cystic Echinococcosis

Alveolar (AE) and cystic (CE) echinococcosis are two parasitic diseases caused by the tapeworms Echinococcus multilocularis and E. granulosus sensu lato (s. l.), respectively. Currently, AE and CE are mainly diagnosed by means of imaging techniques, serology, and clinical and epidemiological data. However, no viability markers that indicate parasite state during infection are available. Extracellular small RNAs (sRNAs) are short non-coding RNAs that can be secreted by cells through association with extracellular vesicles, proteins, or lipoproteins. Circulating sRNAs can show altered expression in pathological states; hence, they are intensively studied as biomarkers for several diseases. Here, we profiled the sRNA transcriptomes of AE and CE patients to identify novel biomarkers to aid in medical decisions when current diagnostic procedures are inconclusive. For this, endogenous and parasitic sRNAs were analyzed by sRNA sequencing in serum from disease negative, positive, and treated patients and patients harboring a non-parasitic lesion. Consequently, 20 differentially expressed sRNAs associated with AE, CE, and/or non-parasitic lesion were identified. Our results represent an in-depth characterization of the effect E. multilocularis and E. granulosus s. l. exert on the extracellular sRNA landscape in human infections and provide a set of novel candidate biomarkers for both AE and CE detection.

The polymorphic analysis of cox1 and cob genes of Echinococcus granulosus in the Ngari region of Tibet in China

Cystic echinococcosis (CE) is an important zoonotic parasitic disease caused by Echinococcus granulosus (E. granulosus). CE seriously threatens human health and the development of animal husbandry. The Ngari region is one of the world's highest endemic regions for CE, while genetic polymorphisms of E. granulosus were unclear. Paraffin slices of liver Cyst were collected from seventy-nine surgical patients with echinococcosis in the Ngari region. DNA was extracted from samples. The cox1 and cob genes of mitochondrial DNA of E. granulosus were simultaneously amplified and sequenced. The sequencing results were compared with the standard sequence (KU925397.1and HF947574.1). Phylogenetic trees and the haplotype network of cob and cox1 genes were constructed and analyzed genotypes of E. granulosus isolated from humans in the Ngari Region of Tibet. Out of 79 hydatid cyst samples collected from surgery patients, 60 isolates were identified as G1/ G3, and two isolates were identified as G6/ G7. Analysis of the cob/ cox1 genes revealed 9/7 mutations resulting in 8/6 haplotypes, respectively. The cob and cox1 neutrality indices computed by Tajima's D and Fu's Fs tests showed high negative values in Echinococcus granulosus sensu stricto (E. granulosus s. s.). The result suggested that E. granulosus in the Ngari region experienced population expansion or a negative selection. We found that G1/ G3 was still the main genotype, and G6/ G7 was found occasionally in humans of the Ngari region. Therefore, we recommend future surveys and control efforts to investigate G1/ G3 and G6/ G7 transmission in the Ngari region.

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.

Comparative Characterization of Mitogenomes From Five Orders of Cestodes (Eucestoda: Tapeworms)

The recognized potential of using mitogenomics in phylogenetics and the more accessible use of high-throughput sequencing (HTS) offer an opportunity to investigate groups of neglected organisms. Here, we leveraged HTS to execute the most comprehensive documentation of mitogenomes for cestodes based on the number of terminals sequenced. We adopted modern approaches to obtain the complete mitogenome sequences of 86 specimens representing five orders of cestodes (three reported for the first time: Phyllobothriidea, “Tetraphyllidea” and Trypanorhyncha). These complete mitogenomes represent an increase of 41% of the mitogenomes available for cestodes (61–147) and an addition of 33% in the representativeness of the cestode orders. The complete mitochondrial genomes are conserved, circular, encoded in the same strand, and transcribed in the same direction, following the pattern observed previously for tapeworms. Their length varies from 13,369 to 13,795 bp, containing 36 genes in total. Except for the Trypanorhyncha specimen, the gene order of the other four cestode orders sequenced here suggests that it could be a synapomorphy for the acetabulate group (with a reversion for taenids). Our results also suggest that no single gene can tell all the evolutionary history contained in the mitogenome. Therefore, cestodes phylogenies based on a single mitochondrial marker may fail to capture their evolutionary history. We predict that such phylogenies would be improved if conducted under a total evidence framework. The characterization of the new mitochondrial genomes is the first step to provide a valuable resource for future studies on the evolutionary relationships of these groups of parasites.

Characterization of a new type of neuronal 5-HT G- protein coupled receptor in the cestode nervous system

Cestodes are platyhelminth parasites with a wide range of hosts that cause neglected diseases. Neurotransmitter signaling is of critical importance for these parasites which lack circulatory, respiratory and digestive systems. For example, serotonin (5-HT) and serotonergic G-protein coupled receptors (5-HT GPCRs) play major roles in cestode motility, development and reproduction. In previous work, we deorphanized a group of 5-HT7 type GPCRs from cestodes. However, little is known about another type of 5-HT GPCR, the 5-HT1 clade, which has been studied in several invertebrate phyla but not in platyhelminthes. Three putative 5-HT GPCRs from Echinococcus canadensis, Mesocestoides vogae (syn. M. corti) and Hymenolepis microstoma were cloned, sequenced and bioinformatically analyzed. Evidence grouped these new sequences within the 5-HT1 clade of GPCRs but differences in highly conserved GPCR motifs were observed. Transcriptomic analysis, heterologous expression and immunolocalization studies were performed to characterize the E. canadensis receptor, called Eca-5-HT_1a. Functional heterologous expression studies showed that Eca-5-HT_1a is highly specific for serotonin. 5-Methoxytryptamine and α-methylserotonin, both known 5-HT GPCR agonists, give stimulatory responses whereas methysergide, a known 5-HT GPCR ligand, give an antagonist response in Eca-5-HT_1a. Mutants obtained by the substitution of key predicted residues resulted in severe impairment of receptor activity, confirming that indeed, these residues have important roles in receptor function. Immunolocalization studies on the protoscolex stage from E. canadensis, showed that Eca-5-HT_1a is localized in branched fibers which correspond to the nervous system of the parasite. The patterns of immunoreactive fibers for Eca-5-HT_1a and for serotonin were intimately intertwined but not identical, suggesting that they are two separate groups of fibers. These data provide the first functional, pharmacological and localization report of a serotonergic receptor that putatively belongs to the 5-HT1 type of GPCRs in cestodes. The serotonergic GPCR characterized here may represent a new target for antiparasitic intervention.

Advances in the treatment, diagnosis, control and scientific understanding of taeniid cestode parasite infections over the past 50 years

In the past 50 years, enormous progress has been made in the diagnosis, treatment and control of taeniid cestode infections/diseases and in the scientific understanding thereof. Most interest in this group of parasites stems from the serious diseases that they cause in humans. It is through this lens that we summarize here the most important breakthroughs that have made a difference to the treatment of human diseases caused by these parasites, reduction in transmission of the taeniid species associated with human disease, or understanding of the parasites' biology likely to impact diagnosis or treatment in the foreseeable future. Key topics discussed are the introduction of anti-cestode drugs, including benzimidazoles and praziquantel, and the development of new imaging modalities that have transformed the diagnosis and post-treatment monitoring of human echinococcoses and neurocysticercosis. The availability of new anti-cestode drugs for use in dogs and a detailed understanding of the transmission dynamics of Echinococcus granulosus sensu lato have underpinned successful programs that have eliminated cystic echinococcosis in some areas of the world and greatly reduced the incidence of infection in others. Despite these successes, cystic and alveolar echinococcosis and neurocysticercosis continue to be prevalent in many parts of the world, requiring new or renewed efforts to prevent the associated taeniid infections. Major advances made in the development of practical vaccines against E. granulosus and Taenia solium will hopefully assist in this endeavour, as might the understanding of the parasites' biology that have come from an elucidation of the nuclear genomes of each of the most important taeniid species causing human diseases.

Cestodes in the genomic era

The first cestode genomes were obtained by an international consortium led by the Wellcome Sanger Institute that included representative institutions from countries where the sequenced parasites have been studied for decades, in part because they are etiological agents of endemic diseases (Argentina, Uruguay, Mexico, Canada, UK, Germany, Switzerland, Ireland, USA, Japan, and China). After this, several complete genomes were obtained reaching 16 species to date. Cestode genomes have smaller relative size compared to other animals including free-living flatworms. Moreover, the features genome size and repeat content seem to differ in the two analyzed orders. Cyclophyllidean species have smaller genomes and with fewer repetitive content than Diphyllobothriidean species. On average, cestode genomes have 13,753 genes with 6 exons per gene and 41% GC content. More than 5,000 shared cestode proteins were accurately annotated by the integration of gene predictions and transcriptome evidence being more than 40% of these proteins of unknown function. Several gene losses and reduction of gene families were found and could be related to the extreme parasitic lifestyle of these species. The application of cutting-edge sequencing technology allowed the characterization of the terminal sequences of chromosomes that possess unique characteristics. Here, we review the current status of knowledge of complete cestode genomes and place it within a comparative genomics perspective. Multidisciplinary work together with the implementation of new technologies will provide valuable information that can certainly improve our chances to finally eradicate or at least control diseases caused by cestodes.

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.

Regulation of DNA methylation on key parasitism genes of Cysticercus cellulosae revealed by integrative epigenomic-transcriptomic analyses

Background

The life cycle of Taenia solium is characterized by different stages of development, requiring various kinds of hosts that can appropriately harbor the eggs (proglottids), the oncospheres, the larvae and the adults. Similar to other metazoan pathogens, T. solium undergoes transcriptional and developmental regulation via epigenetics during its complex lifecycle and host interactions.

Result

In the present study, we integrated whole-genome bisulfite sequencing and RNA-seq technologies to characterize the genome-wide DNA methylation and its effect on transcription of Cysticercus cellulosae of T. solium. We confirm that the T. solium genome in the cysticercus stage is epigenetically modified by DNA methylation in a pattern similar to that of other invertebrate genomes, i.e., sparsely or moderately methylated. We also observed an enrichment of non-CpG methylation in defined genetic elements of the T. solium genome. Furthermore, an integrative analysis of both the transcriptome and the DNA methylome indicated a strong correlation between these two datasets, suggesting that gene expression might be tightly regulated by DNA methylation. Importantly, our data suggested that DNA methylation might play an important role in repressing key parasitism-related genes, including genes encoding excretion-secretion proteins, thereby raising the possibility of targeting DNA methylation processes as a useful strategy in therapeutics of cysticercosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-021-00195-9.

Lost and Found: Piwi and Argonaute Pathways in Flatworms

Platyhelminthes comprise one of the major phyla of invertebrate animals, inhabiting a wide range of ecosystems, and one of the most successful in adapting to parasitic life. Small non-coding RNAs have been implicated in regulating complex developmental transitions in model parasitic species. Notably, parasitic flatworms have lost Piwi RNA pathways but gained a novel Argonaute gene. Herein, we analyzed, contrasted and compared the conservation of small RNA pathways among several free-living species (a paraphyletic group traditionally known as ‘turbellarians’) and parasitic species (organized in the monophyletic clade Neodermata) to disentangle possible adaptations during the transition to parasitism. Our findings showed that complete miRNA and RNAi pathways are present in all analyzed free-living flatworms. Remarkably, whilst all ‘turbellarians’ have Piwi proteins, these were lost in parasitic Neodermantans. Moreover, two clusters of Piwi class Argonaute genes are present in all ‘turbellarians’. Interestingly, we identified a divergent Piwi class Argonaute in free living flatworms exclusively, which we named ‘Fliwi’. In addition, other key proteins of the Piwi pathways were conserved in ‘turbellarians’, while none of them were detected in Neodermatans. Besides Piwi and the canonical Argonaute proteins, a flatworm-specific class of Argonautes (FL-Ago) was identified in the analyzed species confirming its ancestrallity to all Platyhelminthes. Remarkably, this clade was expanded in parasitic Neodermatans, but not in free-living species. These phyla-specific Argonautes showed lower sequence conservation compared to other Argonaute proteins, suggesting that they might have been subjected to high evolutionary rates. However, key residues involved in the interaction with the small RNA and mRNA cleavage in the canonical Argonautes were more conserved in the FL-Agos than in the Piwi Argonautes. Whether this is related to specialized functions and adaptations to parasitism in Neodermatans remains unclear. In conclusion, differences detected in gene conservation, sequence and structure of the Argonaute family suggest tentative biological and evolutionary diversifications that are unique to Platyhelminthes. The remarkable divergencies in the small RNA pathways between free-living and parasitic flatworms indicate that they may have been involved in the adaptation to parasitism of Neodermatans.

Deep amplicon sequencing highlights low intra-host genetic variability of Echinococcus multilocularis and high prevalence of the European-type haplotypes in coyotes and red foxes in Alberta, Canada

Echinococcus multilocularis (Em) is a zoonotic parasite considered a global emergent pathogen. Recent findings indicate that the parasite is expanding its range in North America and that European-type haplotypes are circulating in western Canada. However, genetic analyses are usually conducted only on a few parasites out of thousands of individuals within each definitive host, likely underestimating the prevalence of less common haplotypes. Moreover, mixed infections with several mtDNA haplotypes in the same host have been reported, but their relative abundance within the host was never estimated. We aimed to 1) estimate the frequency of co-infections of different Em haplotypes in coyotes (Canis latrans) and red foxes (Vulpes vulpes) from western Canada and their relative abundance within the definitive hosts, 2) detect less prevalent haplotypes by sampling a larger proportion of the parasite subpopulation per host, and 3) investigate differences in the distribution of Em haplotypes in these main definitive hosts; foxes and coyotes. We extracted DNA from ~10% of the worm subpopulation per host (20 foxes and 47 coyotes) and used deep amplicon sequencing (NGS technology) on four loci, targeting the most polymorphic regions from the mitochondrial genes cox1 (814 bp), nad1 (344 bp), and cob (387 bp). We detected the presence of mixed infections with multiple Em haplotypes and with different Echinococcus species including Em and E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the European-type haplotypes in both hosts. Our results suggest a population expansion of the European over the North American strain in Alberta and a limited distribution of some European-type haplotypes. Our findings indicate that deep amplicon sequencing represents a valuable tool to characterize Em in multiple hosts, to assess the current distribution and possible origins of the European strain in North America. The potential use of next-generation sequencing technologies is particularly important to understand the patterns of geographic expansion of this parasite.

Identification of hydatidosis-related modules and key regulatory genes

Background

Echinococcosis caused by larval of Echinococcus is prevalent all over the world. Although clinical experience showed that the presence of tapeworms could not be found in liver lesions, the repeated infection and aggravation of lesions still occur in the host. Here, this study constructed a multifactor-driven disease-related dysfunction network to explore the potential molecular pathogenesis mechanism in different hosts after E.multilocularis infection.

Method

First, iTRAQ sequencing was performed on human liver infected with E.multilocularis. Second, obtained microRNAs(miRNAs) expression profiles of humans and canine infected with Echinococcus from the GEO database. In addition, we also performed differential expression analysis, protein interaction network analysis, enrichment analysis, and crosstalk analysis to obtain genes and modules related to E.multilocularis infection. Pivot analysis is used to calculate the potential regulatory effects of multiple factors on the module and identify related non-coding RNAs(ncRNAs) and transcription factors(TFs). Finally, we screened the target genes of miRNAs of Echinococcus to further explore its infection mechanism.

Results

A total of 267 differentially expressed proteins from humans and 3,635 differentially expressed genes from canine were obtained. They participated in 16 human-related dysfunction modules and five canine-related dysfunction modules, respectively. Both human and canine dysfunction modules are significantly involved in BMP signaling pathway and TGF-beta signaling pathway. In addition, pivot analysis found that 1,129 ncRNAs and 110 TFs significantly regulated human dysfunction modules, 158 ncRNAs and nine TFs significantly regulated canine dysfunction modules. Surprisingly, the Echinococcus miR-184 plays a role in the pathogenicity regulation by targeting nine TFs and one ncRNA in humans. Similarly, miR-184 can also cause physiological dysfunction by regulating two transcription factors in canine.

Conclusion

The results show that the miRNA-184 of Echinococcus can regulate the pathogenic process through various biological functions and pathways. The results laid a solid theoretical foundation for biologists to further explore the pathogenic mechanism of Echinococcosis.

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.

Evidence for densovirus integrations into tapeworm genomes

BACKGROUND

Tapeworms lack a canonical piRNA-pathway, raising the question of how they can silence existing mobile genetic elements (MGE). Investigation towards the underlying mechanisms requires information on tapeworm transposons which is, however, presently scarce.

METHODS

The presence of densovirus-related sequences in tapeworm genomes was studied by bioinformatic approaches. Available RNA-Seq datasets were mapped against the Echinococcus multilocularis genome to calculate expression levels of densovirus-related genes. Transcription of densovirus loci was further analyzed by sequencing and RT-qPCR.

RESULTS

We herein provide evidence for the presence of densovirus-related elements in a variety of tapeworm genomes. In the high-quality genome of E. multilocularis we identified more than 20 individual densovirus integration loci which contain the information for non-structural and structural virus proteins. The majority of densovirus loci are present as head-to-tail concatemers in isolated repeat containing regions of the genome. In some cases, unique densovirus loci have integrated close to histone gene clusters. We show that some of the densovirus loci of E. multilocularis are actively transcribed, whereas the majority are transcriptionally silent. RT-qPCR data further indicate that densovirus expression mainly occurs in the E. multilocularis stem cell population, which probably forms the germline of this organism. Sequences similar to the non-structural densovirus genes present in E. multilocularis were also identified in the genomes of E. canadensis, E. granulosus, Hydatigera taeniaeformis, Hymenolepis diminuta, Hymenolepis microstoma, Hymenolepis nana, Taenia asiatica, Taenia multiceps, Taenia saginata and Taenia solium.

CONCLUSIONS

Our data indicate that densovirus integration has occurred in many tapeworm species. This is the first report on widespread integration of DNA viruses into cestode genomes. Since only few densovirus integration sites were transcriptionally active in E. multilocularis, our data are relevant for future studies into gene silencing mechanisms in tapeworms. Furthermore, they indicate that densovirus-based vectors might be suitable tools for genetic manipulation of cestodes.

Genome skimming is a low-cost and robust strategy to assemble complete mitochondrial genomes from ethanol preserved specimens in biodiversity studies

Global loss of biodiversity is an ongoing process that concerns both local and global authorities. Studies of biodiversity mainly involve traditional methods using morphological characters and molecular protocols. However, conventional methods are a time consuming and resource demanding task. The development of high-throughput sequencing (HTS) techniques has reshaped the way we explore biodiversity and opened a path to new questions and novel empirical approaches. With the emergence of HTS, sequencing the complete mitochondrial genome became more accessible, and the number of genome sequences published has increased exponentially during the last decades. Despite the current state of knowledge about the potential of mitogenomics in phylogenetics, this is still a relatively under-explored area for a multitude of taxonomic groups, especially for those without commercial relevance, non-models organisms and with preserved DNA. Here we take the first step to assemble and annotate the genomes from HTS data using a new protocol of genome skimming which will offer an opportunity to extend the field of mitogenomics to under-studied organisms. We extracted genomic DNA from specimens preserved in ethanol. We used Nextera XT DNA to prepare indexed paired-end libraries since it is a powerful tool for working with diverse samples, requiring a low amount of input DNA. We sequenced the samples in two different Illumina platform (MiSeq or NextSeq 550). We trimmed raw reads, filtered and had their quality tested accordingly. We performed the assembly using a baiting and iterative mapping strategy, and the annotated the putative mitochondrion through a semi-automatic procedure. We applied the contiguity index to access the completeness of each new mitogenome. Our results reveal the efficiency of the proposed method to recover the whole mitogenomes of preserved DNA from non-model organisms even if there are gene rearrangement in the specimens. Our findings suggest the potential of combining the adequate platform and library to the genome skimming as an innovative approach, which opens a new range of possibilities of its use to obtain molecular data from organisms with different levels of preservation.

Revisiting the Phylogenetic History of Helminths Through Genomics, the Case of the New Echinococcus oligarthrus Genome

The first parasitic helminth genome sequence was published in 2007; since then, only ∼200 genomes have become available, most of them being draft assemblies. Nevertheless, despite the medical and economical global impact of helminthic infections, parasite genomes in public databases are underrepresented. Recently, through an integrative approach involving morphological, genetic, and ecological aspects, we have demonstrated that the complete life cycle of Echinococcus oligarthrus (Cestoda: Taeniidae) is present in South America. The neotropical E. oligarthrus parasite is capable of developing in any felid species and producing human infections. Neotropical echinococcosis is poorly understood yet and requires a complex medical examination to provide the appropriate intervention. Only a few cases of echinococcosis have been unequivocally identified and reported as a consequence of E. oligarthrus infections. Regarding phylogenetics, the analyses of mitogenomes and nuclear datasets have resulted in discordant topologies, and there is no unequivocal taxonomic classification of Echinococcus species so far. In this work, we sequenced and assembled the genome of E. oligarthrus that was isolated from agoutis (Dasyprocta azarae) naturally infected and performed the first comparative genomic study of a neotropical Echinococcus species. The E. oligarthrus genome assembly consisted of 86.22 Mb which showed ∼90% identity and 76.3% coverage with Echinococcus multilocularis and contained the 85.0% of the total expected genes. Genetic variants analysis of whole genome revealed a higher rate of intraspecific genetic variability (23,301 SNPs; 0.22 SNPs/kb) rather than for the genomes of E. multilocularis and Echinococcus canadensis G7 but lower with respect to Echinococcus granulosus G1. Comparative genomics against E. multilocularis, E. granulosus G1, and E. canadensis G7 revealed 38,762, 125,147, and 170,049 homozygous polymorphic sites, respectively, indicating a higher genetic distance between E. oligarthrus and E. granulosus sensu lato species. The SNP distribution in chromosomes revealed a higher SNP density in the longest chromosomes. Phylogenetic analysis using whole-genome SNPs demonstrated that E. oligarthrus is one of the basal species of the genus Echinococcus and is phylogenetically closer to E. multilocularis. This work sheds light on the Echinococcus phylogeny and settles the basis to study sylvatic Echinococcus species and their developmental evolutionary features.

First detection of zoonotic tapeworm Echinococcus granulosus sensu lato genotype G7 in continental Italy

The larval stage of the species complex Echinococcus granulosus sensu lato (s.l.) is the cause of a widespread zoonotic disease known as cystic echinococcosis (CE). The disease is highly prevalent in southern Italy and represents a serious public health issue. The main aim of this study was to characterize E. granulosus s.l. genotypes from wild boar on a continental area of Italy (Campania region), using recently developed mtDNA markers of nad2 and nad5 for reliable identification of different genotypes. Here, nad5 (680 bp) allowed for a clear identification of G1 and G3, whereas a combination of nad2 (714 bp) and nad5 (1394 bp in total) did the same for genotype G7 and its haplogroups G7a and G7b. The results of this study revealed for the first time the presence of genotype G7 in continental Italy. While haplogroup G7b was previously shown to be restricted to the islands of Corsica and Sardinia, here we demonstrate that haplogroup G7b is also present on the mainland of Italy. This work has implications in designing future strategies to reduce CE in Italy.

Histone deacetylase enzymes as potential drug targets of Neglected Tropical Diseases caused by cestodes

Cestode parasites cause neglected diseases, such as echinococcosis and cysticercosis, which represent a significant problem in human and animal health. Benzimidazoles and praziquantel are the only available drugs for chemotherapy and it is therefore important to identify new alternative drugs against cestode parasites. Histone deacetylases (HDACs) are validated drug targets for the treatment of cancer and other diseases, including neglected diseases. However, knowledge of HDACs in cestodes is very scarce. In this work, we investigated cestode HDACs as potential drug targets to develop new therapies against neglected diseases caused by cestodes. Here we showed the full repertoire of HDAC coding genes in several members of the class Cestoda. Between 6 and 7 zinc-dependent HDAC coding genes were identified in the genomes of species from Echinococcus, Taenia, Mesocestoides and Hymenolepis genera. We classified them as Class I and II HDACs and analyzed their transcriptional expression levels throughout developmental stages of Echinococcus spp. We confirmed for the first time the complete HDAC8 nucleotide sequences from Echinococcus canadensis G7 and Mesocestoides corti. Homology models for these proteins showed particular structural features which differentiate them from HDAC8 from Homo sapiens. Furthermore, we showed that Trichostatin A (TSA), a pan-HDAC inhibitor, decreases the viability of M. corti, alters its tegument and morphology and produces an increment of the total amount of acetylated proteins, including acetylated histone H4. These results suggest that HDAC from cestodes are functional and might play important roles on survival and development. The particular structural features observed in cestode HDAC8 proteins suggest that these enzymes could be selectively targeted. This report provides the basis for further studies on cestode HDAC enzymes and for discovery of new HDAC inhibitors for the treatment of neglected diseases caused by cestode parasites.

Transcriptomic profile of two developmental stages of the cestode parasite Mesocestoides corti

Cestode development involves complex morphological and physiological changes. Here, we performed a differential expression analysis of gene transcripts between two developmental stages of the model cestode Mesocestoides corti. A RNA-seq-based approach was used to compare the transcriptomes of the tetrathyridium (larval, TT) and strobilated worm (ST) stages of the parasite. We found 19,053 transcripts, from which ∼45% were complete matches to genes previously annotated in the available M. corti draft genome sequence, ∼24% were considered novel isoforms, and ∼24% were considered potential novel transcripts. Stage-specific transcripts were found for both TTs (66) and STs (136), along with shared transcripts significantly overrepresented in one stage (342 in TTs, and 559 in STs). Differential expression and Gene Ontology term enrichment analyzes provided evidence of upregulation of different sets of transcripts associated with 'cytoskeleton', 'metabolism' and 'oxidation-reduction' processes in each stage, suggesting functional involvement of the corresponding genes with stage-specific features. Transcripts and processes enriched in the TT reflect typical larval processes that occur with the parasite in the intermediate host, such as asexual reproduction and budding, as well as active migration from the peritoneum to the liver and vice versa. In STs, transcripts associated with 'development', 'cell growth', and 'morphogenesis' were enriched, along with processes related to sexual reproduction, represented by the upregulation of numerous transcription factors, protein kinases, and histones. Overall, our results contributed to significantly increase the knowledge on the M. corti gene repertoire and expression profile in two developmental stages. Functional implications for the biology of larval and adult cestode parasites and for host-parasite interactions are discussed.

Echinococcus granulosus genotypes in Iran: a systematic review

Cystic echinococcosis (CE) caused by Echinococcus granulosus sensu lato (s.l.) is a significant zoonosis, especially in developing countries of the Middle East, with many studies focusing on CE genotypes in Iran. We performed a systematic review to determine the exact status of E. granulosus genotypes in the country. We explored English (Pubmed, Scopus, ISI Web of Science and Science Direct) and Persian (Magiran, Iran Medex and Scientific Information Database) databases along with Google Scholar. Our review included 73 studies published prior to the end of 2015. In total, 2952 animal (intermediate and definitive) hosts were examined, and the prevalent genotypes comprised G1 (92.75%) and G6 (4.53%) in sheep, cattle, camels, goats and buffaloes; G3 (2.43%) in five herbivore hosts and dogs; G7 (0.2%) in sheep and goats; and G2 (0.06%) in dogs. G1 was mostly dominant in West Azerbaijan, whereas G3 and G6 were identified most frequently in the provinces of Isfahan and Fars, respectively. Regarding human CE infection, 340 cases were reported from Iran, with the identified genotypes G1 (n = 320), G6 (n = 13) and G3 (n = 7). Most CE-infected humans originated from Isfahan province (168 cases), whereas the lowest number of infected persons was noted in Kerman province (two cases). The information obtained from this systematic review is central to better understanding the biological and epidemiological characteristics of E. granulosus s.l. genotypes in Iran, leading to more comprehensive control strategies.

Unique pharmacological properties of serotoninergic G-protein coupled receptors from cestodes

Background

Cestodes are a diverse group of parasites, some of them being agents of neglected diseases. In cestodes, little is known about the functional properties of G protein coupled receptors (GPCRs) which have proved to be highly druggable targets in other organisms. Notably, serotoninergic G-protein coupled receptors (5-HT GPCRs) play major roles in key functions like movement, development and reproduction in parasites.

Methodology/Principal findings

Three 5-HT GPCRs from Echinococcus granulosus and Mesocestoides corti were cloned, sequenced, bioinformatically analyzed and functionally characterized. Multiple sequence alignment with other GPCRs showed the presence of seven transmembrane segments and conserved motifs but interesting differences were also observed. Phylogenetic analysis grouped these new sequences within the 5-HT7 clade of GPCRs. Molecular modeling showed a striking resemblance in the spatial localization of key residues with their mammalian counterparts. Expression analysis using available RNAseq data showed that both E. granulosus sequences are expressed in larval and adult stages. Localization studies performed in E. granulosus larvae with a fluorescent probe produced a punctiform pattern concentrated in suckers. E. granulosus and M. corti larvae showed an increase in motility in response to serotonin. Heterologous expression revealed elevated levels of cAMP production in response to 5-HT and two of the GPCRs showed extremely high sensitivity to 5-HT (picomolar range). While each of these GPCRs was activated by 5-HT, they exhibit distinct pharmacological properties (5-HT sensitivity, differential responsiveness to ligands).

Conclusions/Significance

These data provide the first functional report of GPCRs in parasitic cestodes. The serotoninergic GPCRs characterized here may represent novel druggable targets for antiparasitic intervention.

Cestode parasites are flatworms with the ability to parasitize almost every vertebrate species. Several of these parasites are etiological agents of neglected diseases prioritized by WHO, such as hydatid disease, or hydatidosis, a zoonosis caused by species of the genus Echinococcus that affects millions of people worldwide. Due to the scarcity of anthelmintic drugs available and the emergence of resistant parasites, the discovery of new anthelmintic drugs is mandatory. Neuromuscular function has been the target of commonly used drugs against parasitic diseases to impact movement, parasite development and reproduction. Here we describe three new proteins, some of them highly expressed in cestodes which could be relevant for motility. Using different approaches, the three proteins were identified as G protein coupled receptors for serotonin, an important neurotransmitter and a known modulator of cestode motility. These new receptors exhibit unique characteristics including a particular sensitivity to serotonin as well as a distinctive pharmacology, which will assist their targeting for chemotherapeutic intervention.

Cestode parasites release extracellular vesicles with microRNAs and immunodiagnostic protein cargo

Intercellular communication is crucial in multiple aspects of cell biology. This interaction can be mediated by several mechanisms including extracellular vesicle (EV) transfer. EV secretion by parasites has been reported in protozoans, trematodes and nematodes. Here we report that this mechanism is present in three different species of cestodes, Taenia crassiceps, Mesocestoides corti and Echinococcus multilocularis. To confirm this we determined, in vitro, the presence of EVs in culture supernatants by transmission electron microscopy. Interestingly, while T. crassiceps and M. corti metacestodes secrete membranous structures into the culture media, similar vesicles were observed in the interface of the germinal and laminated layers of E. multilocularis metacestodes and were hardly detected in culture supernatants. We then determined the protein cargo in the EV-enriched secreted fractions of T. crassiceps and M. corti conditioned media by LC-MS/MS. Among the identified proteins, eukaryotic vesicle-enriched proteins were identified as expected, but also proteins used for cestode disease diagnosis, proteins related to neurotransmission, lipid binding proteins as well as host immunoglobulins and complement factors. Finally, we confirmed by capillary electrophoresis the presence of intravesicular RNA for both parasites and detected microRNAs by reverse transcription-PCR. This is the first report of EV secretion in cestode parasites and of an RNA secretion mechanism. These findings will provide valuable data not only for basic cestode biology but also for the rational search for new diagnostic targets.

Genetic diversity and phylogenetic analysis of EG95 sequences of Echinococcus granulosus: Implications for EG95 vaccine application

OBJECTIVE

To analyse the genetic variability of EG95 sequences and provide guidance for EG95 vaccine application against Echinococcus granulosus (E. granulosus).

METHODS

We analysed EG95 polymorphism by collecting total 97 different E. granulosus isolates from 12 different host species that originated from 10 different countries. Multiple sequence alignments and the homology were performed by Lasergene 1 (DNASTAR Inc., Madison, WI), and the phylogenetic analysis was performed by using MEGA5.1 (CEMI, Tempe, AZ, USA). In addition, linear and conformational epitopes were analysed, including secondary structure, NXT/S glycosylation, fibronectin type III (FnIII) domain and glycosylphosphatidylinositol anchor signal (GPI-anchor). The secondary structure was predicted by PSIPRED method.

RESULTS

Our results indicated that most isolates overall shared 72.6-100% identity in EG95 gene sequence with the published standard EG95 sequence, X90928. However, EG95 gene indeed has polymorphism in different isolates. Phylogenetic analysis showed that different isolates could be divided into three subgroups. Subgroup 1 contained 87 isolates while Subgroup 2 and Subgroup 3 consisted of 3 and 7 isolates, respectively. Four sequences cloned from oncosphere shared a high identity with the parental sequence of the current vaccine, X90928, and they belonged to Subgroup 1. However, in comparison to X90928, several amino acid mutations occurred in most isolates besides oncosphere, which potentially altered the immunodominant linear epitopes, glycosylation sites and secondary structures in EG95 genes. All these variations might change their previous antigenicity and thereby affecting the efficacy of current EG95 vaccine.

CONCLUSIONS

This study reveals the genetic variability of EG95 sequences in different E. granulosus isolates, and proposed that more vaccination trials would be needed to test the effectiveness of current EG95 vaccine against distinct isolates in different countries.

Genome-wide identification of microRNA targets in the neglected disease pathogens of the genus Echinococcus

MicroRNAs (miRNAs), a class of small non-coding RNAs, are key regulators of gene expression at post-transcriptional level and play essential roles in biological processes such as development. MiRNAs silence target mRNAs by binding to complementary sequences in the 3'untranslated regions (3'UTRs). The parasitic helminths of the genus Echinococcus are the causative agents of echinococcosis, a zoonotic neglected disease. In previous work, we performed a comprehensive identification and characterization of Echinococcus miRNAs. However, current knowledge about their targets is limited. Since target prediction algorithms rely on complementarity between 3'UTRs and miRNA sequences, a major limitation is the lack of accurate sequence information of 3'UTR for most species including parasitic helminths. We performed RNA-seq and developed a pipeline that integrates the transcriptomic data with available genomic data of this parasite in order to identify 3'UTRs of Echinococcus canadensis. The high confidence set of 3'UTRs obtained allowed the prediction of miRNA targets in Echinococcus through a bioinformatic approach. We performed for the first time a comparative analysis of miRNA targets in Echinococcus and Taenia. We found that many evolutionarily conserved target sites in Echinococcus and Taenia may be functional and under selective pressure. Signaling pathways such as MAPK and Wnt were among the most represented pathways indicating miRNA roles in parasite growth and development. Genome-wide identification and characterization of miRNA target genes in Echinococcus provide valuable information to guide experimental studies in order to understand miRNA functions in the parasites biology. miRNAs involved in essential functions, especially those being absent in the host or showing sequence divergence with respect to host orthologs, might be considered as novel therapeutic targets for echinococcosis control.