MicroRNAs (miRNAs), a class of small non-coding RNAs, are key regulators of gene expression at post-transcriptional level and play essential roles in fundamental biological processes such as development and metabolism. The particular developmental and metabolic characteristics of cestode parasites highlight the importance of studying miRNA gene regulation in these organisms. Here, we perform a comprehensive analysis of miRNAs in the parasitic cestode Echinococcus canadensis G7, one of the causative agents of the neglected zoonotic disease cystic echinococcosis. Small RNA libraries from protoscoleces and cyst walls of E. canadensis G7 and protoscoleces of E. granulosus sensu stricto G1 were sequenced using Illumina technology. For miRNA prediction, miRDeep2 core algorithm was used. The output list of candidate precursors was manually curated to generate a high confidence set of miRNAs. Differential expression analysis of miRNAs between stages or species was estimated with DESeq. Expression levels of selected miRNAs were validated using poly-A RT-qPCR. In this study we used a high-throughput approach and found transcriptional evidence of 37 miRNAs thus expanding the miRNA repertoire of E. canadensis G7. Differential expression analysis showed highly regulated miRNAs between life cycle stages, suggesting a role in maintaining the features of each developmental stage or in the regulation of developmental timing. In this work we characterize conserved and novel Echinococcus miRNAs which represent 30 unique miRNA families. Here we confirmed the remarkable loss of conserved miRNA families in E. canadensis, reflecting their low morphological complexity and high adaptation to parasitism. We performed the first in-depth study profiling of small RNAs in the zoonotic parasite E. canadensis G7. We found that miRNAs are the preponderant small RNA silencing molecules, suggesting that these small RNAs could be an essential mechanism of gene regulation in this species. We also identified both parasite specific and divergent miRNAs which are potential biomarkers of infection. This study will provide valuable information for better understanding of the complex biology of this parasite and could help to find new potential targets for therapy and/or diagnosis.

Echinococcosis represents a major public health problem worldwide and is considered a neglected disease by the World Health Organization. The etiological agents are Echinococcus tapeworms, which display elaborate developmental traits that imply a complex control of gene expression. MicroRNAs (miRNAs), a class of small regulatory RNAs, are involved in the regulation of many biological processes such as development and metabolism. They act through the repression of messenger RNAs (mRNAs) usually by binding to the 3’ untranslated region (3’UTR). Previously, we described the miRNome of several Echinococcus species and found that miRNAs are highly expressed in all life cycle stages, suggesting an important role in gene expression regulation. However, studying the role of miRNAs in helminth biology remains a challenge. To develop methodology for functional analysis of miRNAs in tapeworms, we performed miRNA knockdown experiments in primary cell cultures of Echinococcus multilocularis, which mimic the development of metacestode vesicles from parasite stem cells in vitro. First, we analysed the miRNA repertoire of E. multilocularis primary cells by small RNA-seq and found that miR-71, a bilaterian miRNA absent in vertebrate hosts, is one of the top five most expressed miRNAs. Using genomic information and bioinformatic algorithms for miRNA binding prediction, we found a high number of potential miR-71 targets in E. multilocularis. Inhibition of miRNAs can be achieved by transfection of antisense oligonucleotides (anti-miRs) that block miRNA function. To this end, we evaluated a variety of chemically modified anti-miRs for miR-71 knockdown. Electroporation of primary cells with 2’-O-methyl modified anti-miR-71 led to significantly reduced miR-71 levels. Transcriptomic analyses showed that several predicted miR-71 targets were up-regulated in anti-miR-treated primary cells, including genes potentially involved in parasite development, host parasite interaction, and several genes of as yet unknown function. Notably, miR-71-silenced primary cell cultures showed a strikingly different phenotype as from control cells and did not develop into fully mature metacestodes. These findings indicate an important function of miR-71 in Echinococcus development and provide, for the first time, methodology to functionally study miRNAs in a tapeworm. Echinococcosis, caused by the larval stages of tapeworms of the genus Echinococcus, is a neglected disease that affects millions of people world-wide. These parasites show elaborate developmental features that rely on a complex control of gene expression. microRNAs are small molecules which have been discovered in the last decades and control gene expression in animals, plants and viruses. microRNAs are highly expressed in several tapeworms but their biological function in these parasites is unknown. Assuming that microRNAs will be important for parasite development, we analysed the function of these molecules in Echinococcus multilocularis, employing an in vitro model that mimics the first developmental transitions which occur in the human host. By artificially decreasing the concentration of the highest expressed microRNA, we observed phenotypic alterations and inhibition of development. In addition, we identified possible mRNA molecules targeted by microRNAs and found that some of these are known for being involved in developmental processes in other organisms. This work provides novel methodology to study microRNA function in tapeworms. Furthermore, highly expressed parasite microRNAs that are absent in the host but fulfil an important role in parasite developmental processes can serve as selective drug targets against the underlying diseases.

The parasite Echinococcus canadensis (G7) (phylum Platyhelminthes, class Cestoda) is one of the causative agents of echinococcosis. Echinococcosis is a worldwide chronic zoonosis affecting humans as well as domestic and wild mammals, which has been reported as a prioritized neglected disease by the World Health Organisation. No genomic data, comparative genomic analyses or efficient therapeutic and diagnostic tools are available for this severe disease. The information presented in this study will help to understand the peculiar biological characters and to design species-specific control tools. We sequenced, assembled and annotated the 115-Mb genome of E. canadensis (G7). Comparative genomic analyses using whole genome data of three Echinococcus species not only confirmed the status of E. canadensis (G7) as a separate species but also demonstrated a high nucleotide sequences divergence in relation to E. granulosus (G1). The E. canadensis (G7) genome contains 11,449 genes with a core set of 881 orthologs shared among five cestode species. Comparative genomics revealed that there are more single nucleotide polymorphisms (SNPs) between E. canadensis (G7) and E. granulosus (G1) than between E. canadensis (G7) and E. multilocularis. This result was unexpected since E. canadensis (G7) and E. granulosus (G1) were considered to belong to the species complex E. granulosus sensu lato. We described SNPs in known drug targets and metabolism genes in the E. canadensis (G7) genome. Regarding gene regulation, we analysed three particular features: CpG island distribution along the three Echinococcus genomes, DNA methylation system and small RNA pathway. The results suggest the occurrence of yet unknown gene regulation mechanisms in Echinococcus. This is the first work that addresses Echinococcus comparative genomics. The resources presented here will promote the study of mechanisms of parasite development as well as new tools for drug discovery. The availability of a high-quality genome assembly is critical for fully exploring the biology of a pathogenic organism. The E. canadensis (G7) genome presented in this study provides a unique opportunity to address the genetic diversity among the genus Echinococcus and its particular developmental features. At present, there is no unequivocal taxonomic classification of Echinococcus species; however, the genome-wide SNPs analysis performed here revealed the phylogenetic distance among these three Echinococcus species. Additional cestode genomes need to be sequenced to be able to resolve their phylogeny.

Serotonin (5-hydroxytryptamine, 5-HT) is an important neuroactive and morphogenetic molecule in several metazoan phyla, including flatworms. Serotoninergic nervous system studies are incomplete and 5-HT function/s are unknown in Echinococcus spp., the flatworm parasites that cause hydatid disease. In the present work, we searched for genes of the serotoninergic pathway and performed immunocytochemical and functional analyses of 5-HT in Echinococcus spp. Bioinformatic analysis using the recently available Echinococcus multilocularis and Echinococcus granulosus genomes suggests the presence of genes encoding enzymes, receptors and transporters participating in 5-HT synthesis, sensing and transport in these parasites. However, some components of the pathway could not be identified, suggesting loss or divergence of parasite homologous genes. The serotoninergic neuroanatomy study performed by confocal scanning laser microscopy on different E. granulosus stages showed an increasing level of complexity when the protoscolex develops towards the adult stage and a progressive diminution when the parasite develops towards the metacestode stage. The role of 5-HT as a neurotransmitter in E. granulosus was evaluated by determining the effect of this substance on protoscolex motility. The addition of 5-HT to protoscoleces induced a significant increase in motility for short time periods. Preincubation with 100 μM citalopram, a known 5-HT transporter inhibitor, abolished the 5-HT-induced increase in motility, indicating that the effect could be mediated by a 5-HT transporter. Incubation of protoscoleces with 5-HT for time periods of several days induced a progressive differentiation towards the metacestode stage. The results indicate that 5-HT could have nervous and prenervous roles during Echinococcus spp. development. Taking into account the important roles of 5-HT in parasite biology and the divergence of 5-HT pathway genes with respect to human counterparts, the serotoninergic system could be considered as an amenable drug target against hydatid disease.