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.

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, serotonergic G-protein coupled receptors (5-HT GPCRs) play major roles in key functions like movement, development and reproduction in parasites. 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). 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.

In December 2019, rising pneumonia cases caused by a novel β-coronavirus (SARS-CoV-2) occurred in Wuhan, China, which has rapidly spread worldwide, causing thousands of deaths. The WHO declared the SARS-CoV-2 outbreak as a public health emergency of international concern, since then several scientists are dedicated to its study. It has been observed that many human viruses have codon usage biases that match highly expressed proteins in the tissues they infect and depend on the host cell machinery for the replication and co-evolution. In this work, we analysed 91 molecular features and codon usage patterns for 339 viral genes and 463 human genes that consisted of 677,873 codon positions. Hereby, we selected the highly expressed genes from human lung tissue to perform computational studies that permit to compare their molecular features with those of SARS, SARS-CoV-2 and MERS genes. The integrated analysis of all the features revealed that certain viral genes and overexpressed human genes have similar codon usage patterns. The main pattern was the A/T bias that together with other features could propitiate the viral infection, enhanced by a host dependant specialization of the translation machinery of only some of the overexpressed genes. The envelope protein E, the membrane glycoprotein M and ORF7 could be further benefited. This could be the key for a facilitated translation and viral replication conducting to different comorbidities depending on the genetic variability of population due to the host translation machinery. This is the first codon usage approach that reveals which human genes could be potentially deregulated due to the codon usage similarities between the host and the viral genes when the virus is already inside the human cells of the lung tissues. Our work leaded to the identification of additional highly expressed human genes which are not the usual suspects but might play a role in the viral infection and settle the basis for further research in the field of human genetics associated with new viral infections. To identify the genes that could be deregulated under a viral infection is important to predict the collateral effects and determine which individuals would be more susceptible based on their genetic features and comorbidities associated.

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.

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.