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.

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.