Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators.

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

SOD1G93A Astrocyte-Derived Extracellular Vesicles Induce Motor Neuron Death by a miRNA-155-5p-Mediated Mechanism.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by upper and lower motor neuron (MN) degeneration. Astrocytes surrounding MNs are known to modulate ALS progression. When cocultured with astrocytes overexpressing the ALS-linked mutant Cu/Zn superoxide dismutase (SOD1G93A) or when cultured with conditioned medium from SOD1G93A astrocytes, MN survival is reduced. The exact mechanism of this neurotoxic effect is unknown. Astrocytes secrete extracellular vesicles (EVs) that transport protein, mRNA, and microRNA species from one cell to another. The size and protein markers characteristic of exosomes were observed in the EVs obtained from cultured astrocytes, indicating their abundance in exosomes. Here, we analyzed the microRNA content of the exosomes derived from SOD1G93A astrocytes and evaluated their role in MN survival. Purified MNs exposed to SOD1G93A astrocyte-derived exosomes showed reduced survival and neurite length compared to those exposed to exosomes derived from non-transgenic (non-Tg) astrocytes. Analysis of the miRNA content of the exosomes revealed that miR-155-5p and miR-582-3p are differentially expressed in SOD1G93A exosomes compared with exosomes from non-Tg astrocytes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicates that miR-155-5p and miR-582-3p predicted targets are enriched in the neurotrophin signaling pathway. Importantly, when levels of miR-155-5p were reduced by incubation with a specific antagomir, SOD1G93A exosomes did not affect MN survival or neurite length. These results demonstrate that SOD1G93A-derived exosomes are sufficient to induce MN death, and miRNA-155-5p contributes to this effect. miRNA-155-5p may offer a new therapeutic target to modulate disease progression in ALS.

Special considerations for studies of extracellular vesicles from parasitic helminths: A community-led roadmap to increase rigour and reproducibility.

Over the last decade, research interest in defining how extracellular vesicles (EVs) shape cross-species communication has grown rapidly. Parasitic helminths, worm species found in the phyla Nematoda and Platyhelminthes, are well-recognised manipulators of host immune function and physiology. Emerging evidence supports a role for helminth-derived EVs in these processes and highlights EVs as an important participant in cross-phylum communication. While the mammalian EV field is guided by a community-agreed framework for studying EVs derived from model organisms or cell systems [e.g., Minimal Information for Studies of Extracellular Vesicles (MISEV)], the helminth community requires a supplementary set of principles due to the additional challenges that accompany working with such divergent organisms. These challenges include, but are not limited to, generating sufficient quantities of EVs for descriptive or functional studies, defining pan-helminth EV markers, genetically modifying these organisms, and identifying rigorous methodologies for in vitro and in vivo studies. Here, we outline best practices for those investigating the biology of helminth-derived EVs to complement the MISEV guidelines. We summarise community-agreed standards for studying EVs derived from this broad set of non-model organisms, raise awareness of issues associated with helminth EVs and provide future perspectives for how progress in the field will be achieved.

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.

Role of Fasciola hepatica Small RNAs in the Interaction With the Mammalian Host.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression being involved in many different biological processes and play a key role in developmental timing. Additionally, recent studies have shown that miRNAs released from parasites are capable of regulating the expression of host genes. In the present work, we studied the expression patterns of ncRNAs of various intra-mammalian life-cycle stages of the liver fluke, Fasciola hepatica, as well as those packaged into extracellular vesicles and shed by the adult fluke. The miRNA expression profile of the intra-mammalian stages shows important variations, despite a set of predominant miRNAs that are highly expressed across all stages. No substantial variations in miRNA expression between dormant and activated metacercariae were detected, suggesting that they might not be central players in regulating fluke gene expression during this crucial step in the invasion of the definitive host. We generated a curated pipeline for the prediction of putative target genes that reports only sites conserved between three different prediction approaches. This pipeline was tested against an iso-seq curated database of the 3' UTR regions of F. hepatica genes to detect miRNA regulation networks within liver fluke. Several functions related to the host immune response or modulation were enriched among the targets of the most highly expressed parasite miRNAs, stressing that they might be key players during the establishment and maintenance of infection. Additionally, we detected fragments derived from the processing of tRNAs, in all developmental stages analyzed, and documented the presence of novel long tRNA fragments enriched in vesicles. We confirmed the presence of at least 5 putative vault RNAs (vtRNAs), that are expressed across different stages and enriched in vesicles. The presence of tRNA fragments and vtRNAs in vesicles raise the possibility that they could be involved in the host-parasite interaction.

Adaptive Radiation of the Flukes of the Family Fasciolidae Inferred from Genome-Wide Comparisons of Key Species.

Liver and intestinal flukes of the family Fasciolidae cause zoonotic food-borne infections that impact both agriculture and human health throughout the world. Their evolutionary history and the genetic basis underlying their phenotypic and ecological diversity are not well understood. To close that knowledge gap, we compared the whole genomes of Fasciola hepatica, Fasciola gigantica, and Fasciolopsis buski and determined that the split between Fasciolopsis and Fasciola took place ∼90 Ma in the late Cretaceous period, and that between 65 and 50 Ma an intermediate host switch and a shift from intestinal to hepatic habitats occurred in the Fasciola lineage. The rapid climatic and ecological changes occurring during this period may have contributed to the adaptive radiation of these flukes. Expansion of cathepsins, fatty-acid-binding proteins, protein disulfide-isomerases, and molecular chaperones in the genus Fasciola highlights the significance of excretory-secretory proteins in these liver-dwelling flukes. Fasciola hepatica and Fasciola gigantica diverged ∼5 Ma near the Miocene-Pliocene boundary that coincides with reduced faunal exchange between Africa and Eurasia. Severe decrease in the effective population size ∼10 ka in Fasciola is consistent with a founder effect associated with its recent global spread through ruminant domestication. G-protein-coupled receptors may have key roles in adaptation of physiology and behavior to new ecological niches. This study has provided novel insights about the genome evolution of these important pathogens, has generated genomic resources to enable development of improved interventions and diagnosis, and has laid a solid foundation for genomic epidemiology to trace drug resistance and to aid surveillance.

Compositional Analysis of Flatworm Genomes Shows Strong Codon Usage Biases Across All Classes.

In the present work, we performed a comparative genome-wide analysis of 22 species representative of the main clades and lifestyles of the phylum Platyhelminthes. We selected a set of 700 orthologous genes conserved in all species, measuring changes in GC content, codon, and amino acid usage in orthologous positions. Values of 3rd codon position GC spanned over a wide range, allowing to discriminate two distinctive clusters within freshwater turbellarians, Cestodes and Trematodes respectively. Furthermore, a hierarchical clustering of codon usage data differs remarkably from the phylogenetic tree. Additionally, we detected a synonymous codon usage bias that was more dramatic in extreme GC-poor or GC-rich genomes, i.e., GC-poor Schistosomes preferred to use AT-rich terminated synonymous codons, while GC-rich M. lignano showed the opposite behavior. Interestingly, these biases impacted the amino acidic usage, with preferred amino acids encoded by codons following the GC content trend. These are associated with non-synonymous substitutions at orthologous positions. The detailed analysis of the synonymous and non-synonymous changes provides evidence for a two-hit mechanism where both mutation and selection forces drive the diverse coding strategies of flatworms.

Pleiotropic alterations in gene expression in Latin American Fasciola hepatica isolates with different susceptibility to drugs.

BACKGROUND: Fasciola hepatica is the main agent of fasciolosis, a zoonotic disease affecting livestock worldwide, and an emerging food-borne disease in humans. Even when effective treatments are available, drugs are costly and can result in tolerance, liver damage and normally they do not prevent reinfection. Drug-resistant strains in livestock have been reported in various countries and, more worryingly, drug resistance in human cases has emerged in South America. The present study aims to characterize the transcriptome of two South American resistant isolates, the Cajamarca isolate from Peru, resistant to both triclabendazole and albendazole (TCBZR/ABZR) and the Rubino isolate from Uruguay, resistant to ABZ (TCBZS/ABZR), and compare them to a sensitive strain (Cenapa, Mexico, TCBZS/ABZS) to reveal putative molecular mechanisms leading to drug resistance. RESULTS: We observed a major reduction in transcription in the Cajamarca TCBZR/ABZR isolate in comparison to the other isolates. While most of the differentially expressed genes are still unannotated, several trends could be detected. Specific reduction in the expression levels of cytoskeleton proteins was consistent with a role of tubulins as putative targets of triclabendazole (TCBZ). A marked reduction of adenylate cyclase might be underlying pleiotropic effects on diverse metabolic pathways of the parasite. Upregulation of GST mu isoforms suggests this detoxifying mechanism as one of the strategies associated with resistance. CONCLUSIONS: Our results stress the value of transcriptomic approaches as a means of providing novel insights to advance the understanding of drug mode of action and drug resistance. The results provide evidence for pleiotropic variations in drug-resistant isolates consistent with early observations of TCBZ and ABZ effects and recent proteomic findings.

Conservation and diversification of small RNA pathways within flatworms.

BACKGROUND: Small non-coding RNAs, including miRNAs, and gene silencing mediated by RNA interference have been described in free-living and parasitic lineages of flatworms, but only few key factors of the small RNA pathways have been exhaustively investigated in a limited number of species. The availability of flatworm draft genomes and predicted proteomes allowed us to perform an extended survey of the genes involved in small non-coding RNA pathways in this phylum. RESULTS: Overall, findings show that the small non-coding RNA pathways are conserved in all the analyzed flatworm linages; however notable peculiarities were identified. While Piwi genes are amplified in free-living worms they are completely absent in all parasitic species. Remarkably all flatworms share a specific Argonaute family (FL-Ago) that has been independently amplified in different lineages. Other key factors such as Dicer are also duplicated, with Dicer-2 showing structural differences between trematodes, cestodes and free-living flatworms. Similarly, a very divergent GW182 Argonaute interacting protein was identified in all flatworm linages. Contrasting to this, genes involved in the amplification of the RNAi interfering signal were detected only in the ancestral free living species Macrostomum lignano. We here described all the putative small RNA pathways present in both free living and parasitic flatworm lineages. CONCLUSION: These findings highlight innovations specifically evolved in platyhelminths presumably associated with novel mechanisms of gene expression regulation mediated by small RNA pathways that differ to what has been classically described in model organisms. Understanding these phylum-specific innovations and the differences between free living and parasitic species might provide clues to adaptations to parasitism, and would be relevant for gene-silencing technology development for parasitic flatworms that infect hundreds of million people worldwide.

Genomes of Fasciola hepatica from the Americas Reveal Colonization with Neorickettsia Endobacteria Related to the Agents of Potomac Horse and Human Sennetsu Fevers.

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis' gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans.

The miRnome of Fasciola hepatica juveniles endorses the existence of a reduced set of highly divergent micro RNAs in parasitic flatworms.

The liver fluke Fasciola hepatica is a foodborne zoonotic parasite affecting livestock worldwide, with increasing relevance in human health. The first developmental stage that the host meets after ingestion of the parasite is the newly excysted juvenile, that actively transverses the gut wall and migrates to its final location in the liver. The regulation of the early developmental events in newly excysted juveniles is still poorly understood and a relevant target for control strategies. Here we investigated the putative involvement of small regulatory RNAs in the invasion process. The small RNA population of the newly excysted juvenile fall into two classes, one represented by micro (mi)RNAs and a secondary group of larger (32-33 nucleotides) tRNA-derived sequences. We identified 40 different miRNAs, most of those belonging to ancient miRNAs conserved in protostomes and metazoans, notably with a highly predominant miR-125b variant. Remarkably, several protostomian and metazoan conserved families were not detected in consonance with previous reports of drastic miRnome reduction in parasitic flatworms. Additionally, a set of five novel miRNAs was identified, probably associated with specific gene regulation expression needs in F. hepatica. While sequence conservation in mature miRNA is high across the metazoan tree, we observed that flatworm miRNAs are more divergent, suggesting that mutation rates in parasitic flatworms could be high. Finally, the distinctive presence of tRNA-derived sequences, mostly 5' tRNA halves of selected tRNAs in the small RNA population of newly excysted juveniles, raises the possibility that both miRNA and tRNA fragments participate in the regulation of gene expression in this parasite.