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.

A novel Tetrahymena thermophila sterol C-22 desaturase belongs to the fatty acid hydroxylase/desaturase superfamily.

Sterols in eukaryotic cells play important roles in modulating membrane fluidity and in cell signaling and trafficking. During evolution, a combination of gene losses and acquisitions gave rise to an extraordinary diversity of sterols in different organisms. The sterol C-22 desaturase identified in plants and fungi as a cytochrome P-450 monooxygenase evolved from the first eukaryotic cytochrome P450 and was lost in many lineages. Although the ciliate Tetrahymena thermophila desaturates sterols at the C-22 position, no cytochrome P-450 orthologs are present in the genome. Here, we aim to identify the genes responsible for the desaturation as well as their probable origin. We used gene knockout and yeast heterologous expression approaches to identify two putative genes, retrieved from a previous transcriptomic analysis, as sterol C-22 desaturases. Furthermore, we demonstrate using bioinformatics and evolutionary analyses that both genes encode a novel type of sterol C-22 desaturase that belongs to the large fatty acid hydroxylase/desaturase superfamily and the genes originated by genetic duplication prior to functional diversification. These results stress the widespread existence of nonhomologous isofunctional enzymes among different lineages of the tree of life as well as the suitability for the use of T. thermophila as a valuable model to investigate the evolutionary process of large enzyme families.

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.

Molecular features similarities between SARS-CoV-2, SARS, MERS and key human genes could favour the viral infections and trigger collateral effects.

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.

Fcγ Receptor IIa (FCGR2A) Polymorphism Is Associated With Severe Respiratory Syncytial Virus Disease in Argentinian Infants.

Background: Most patients with respiratory syncytial virus (RSV) infection requiring hospitalization have no risk factors for severe disease. Genetic variation in the receptor for the Fc portion of IgG (FcγR) determines their affinity for IgG subclasses driving innate and adaptive antiviral immunity. We investigated the relationship between FcγRIIa-H131R polymorphism and RSV disease. Methods: Blood samples were collected from 182 infants ≤24-month-old (50 uninfected, 114 RSV-infected with moderate course and 18 suffering severe disease). FcγRIIa-H131R SNP genotypic frequencies (HH, HR, RR) and anti-RSV IgG1, IgG2 and IgG3 levels were studied. Results: Genotypic frequencies for FcγRIIa-H131R SNP were comparable between uninfected and RSV-infected infants. In contrast, we found a significant higher frequency of HH genotype in severe RSV-infected children compared to moderate patients. Among severe group, HH infants presented more factors associated to severity than HR or RR patients did. Furthermore, compared to moderate RSV-infected infants, severe patients showed higher levels of anti-RSV IgG1 and IgG3. Conclusions: We found an association between an FcγRIIa (H131) polymorphism and severe RSV disease, which points towards a critical role for interactions between FcγRs and immune complexes in RSV pathogenesis. This genetic factor could also predict the worse outcome and identify those infants at risk during hospitalization.

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.

Whole genome analysis of codon usage in Echinococcus.

The species of the genus Echinococcus are parasitic platyhelminths that cause echinococcosis and exert a global burden on public and animal health. Here we performed codon usage bias and comparative genomic analyses using whole genome and expression data of three Echinococcus species. The study of 4,710,883 codons, two orders of magnitude more than in previous research works, showed that the codon usage in Echinococcus genes is biased towards the pyrimidines T and C ending codons, with an average effective number of codons equal to 57 revealing a low codon usage bias. The gene annotations and the expression profile of 7613 genes allowed to accurately determine 27 optimal codons for the Echinococcus species, most of them ending in G/C. Approximately the 30% of Echinococcus genes analysed exhibits higher codon usage bias as well as a higher expression profile. Neutrality-plots demonstrated that the selection pressure is the main evolutionary force shaping the codon usage with a contribution of 80%. Comparative genome analyses among several tapeworm species revealed that codon usage patterns are a conserved trait in cestodes parasites. Since cestodes parasites take advantage of the host protein synthesis pathways, this study could provide valuable information associated with the parasite-host relationship that would be useful to determine which host's factors are relevant for shaping the codon usage.

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.

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

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.