Differential expression of small RNA pathway genes associated with the Biomphalaria glabrata/Schistosoma mansoni interaction

Application of microphysiological systems to unravel the mechanisms of schistosomiasis egg extravasation

Despite decades of control efforts, the prevalence of schistosomiasis remains high in many endemic regions, posing significant challenges to global health. One of the key factors contributing to the persistence of the disease is the complex life cycle of the Schistosoma parasite, the causative agent, which involves multiple stages of development and intricate interactions with its mammalian hosts and snails. Among the various stages of the parasite lifecycle, the deposition of eggs and their migration through host tissues is significant, as they initiate the onset of the disease pathology by inducing inflammatory reactions and tissue damage. However, our understanding of the mechanisms underlying Schistosoma egg extravasation remains limited, hindering efforts to develop effective interventions. Microphysiological systems, particularly organ-on-a-chip systems, offer a promising approach to study this phenomenon in a controlled experimental setting because they allow the replication of physiological microenvironments in vitro. This review provides an overview of schistosomiasis, introduces the concept of organ-on-a-chip technology, and discusses its potential applications in the field of schistosomiasis research.

ARGONAUTE2 Localizes to Sites of Sporocysts in the Schistosome-Infected Snail, Biomphalaria glabrata

MicroRNAs (miRNAs) are a class of small regulatory RNA that are generated via core protein machinery. The miRNAs direct gene-silencing mechanisms to mediate an essential role in gene expression regulation. In mollusks, miRNAs have been demonstrated to be required to regulate gene expression in various biological processes, including normal development, immune responses, reproduction, and stress adaptation. In this study, we aimed to establishment the requirement of the miRNA pathway as part of the molecular response of exposure of Biomphalaria glabrata (snail host) to Schistosoma mansoni (trematode parasite). Initially, the core pieces of miRNA pathway protein machinery, i.e., Drosha, DGCR8, Exportin-5, Ran, and Dicer, together with the central RNA-induced silencing complex (RISC) effector protein Argonaute2 (Ago2) were elucidated from the B. glabrata genome. Following exposure of B. glabrata to S. mansoni miracidia, we identified significant expression up-regulation of all identified pieces of miRNA pathway protein machinery, except for Exportin-5, at 16 h post exposure. For Ago2, we went on to show that the Bgl-Ago2 protein was localized to regions surrounding the sporocysts in the digestive gland of infected snails 20 days post parasite exposure. In addition to documenting elevated miRNA pathway protein machinery expression at the early post-exposure time point, a total of 13 known B. glabrata miRNAs were significantly differentially expressed. Of these thirteen B. glabrata miRNAs responsive to S. mansoni miracidia exposure, five were significantly reduced in their abundance, and correspondingly, these five miRNAs were determined to putatively target six genes with significantly elevated expression and that have been previously associated with immune responses in other animal species, including humans. In conclusion, this study demonstrates the central importance of a functional miRNA pathway in snails, which potentially forms a critical component of the immune response of snails to parasite exposure. Further, the data reported in this study provide additional evidence of the complexity of the molecular response of B. glabrata to S. mansoni infection: a molecular response that could be targeted in the future to overcome parasite infection and, in turn, human schistosomiasis.

Different metazoan parasites, different transcriptomic responses, with new insights on parasitic castration by digenetic trematodes in the schistosome vector snail Biomphalaria glabrata

BACKGROUND

Gastropods of the genus Biomphalaria (Family Planorbidae) are exploited as vectors by Schistosoma mansoni, the most common causative agent of human intestinal schistosomiasis. Using improved genomic resources, overviews of how Biomphalaria responds to S. mansoni and other metazoan parasites can provide unique insights into the reproductive, immune, and other systems of invertebrate hosts, and their responses to parasite challenges.

RESULTS

Using Illumina-based RNA-Seq, we compared the responses of iM line B. glabrata at 2, 8, and 40 days post-infection (dpi) to single infections with S. mansoni, Echinostoma paraensei (both digenetic trematodes) or Daubaylia potomaca (a nematode parasite of planorbid snails). Responses were compared to unexposed time-matched control snails. We observed: (1) each parasite provoked a distinctive response with a predominance of down-regulated snail genes at all time points following exposure to either trematode, and of up-regulated genes at 8 and especially 40dpi following nematode exposure; (2) At 2 and 8dpi with either trematode, several snail genes associated with gametogenesis (particularly spermatogenesis) were down-regulated. Regarding the phenomenon of trematode-mediated parasitic castration in molluscs, we define for the first time a complement of host genes that are targeted, as early as 2dpi when trematode larvae are still small; (3) Differential gene expression of snails with trematode infection at 40dpi, when snails were shedding cercariae, was unexpectedly modest and revealed down-regulation of genes involved in the production of egg mass proteins and peptide processing; and (4) surprisingly, D. potomaca provoked up-regulation at 40dpi of many of the reproduction-related snail genes noted to be down-regulated at 2 and 8dpi following trematode infection. Happening at a time when B. glabrata began to succumb to D. potomaca, we hypothesize this response represents an unexpected form of fecundity compensation. We also document expression patterns for other Biomphalaria gene families, including fibrinogen domain-containing proteins (FReDs), C-type lectins, G-protein coupled receptors, biomphalysins, and protease and protease inhibitors.

CONCLUSIONS

Our study is relevant in identifying several genes involved in reproduction that are targeted by parasites in the vector snail B. glabrata and that might be amenable to manipulation to minimize their ability to serve as vectors of schistosomes.

Multifunctional Roles of MicroRNAs in Schistosomiasis

Schistosomiasis is a parasitic disease that is caused by helminths of the genus Schistosoma. The dioecious schistosomes mate and lay eggs after undergoing a complex life cycle. Schistosome eggs are mostly responsible for the transmission of schistosomiasis and chronic fibrotic disease induced by egg antigens is the main cause of the high mortality rate. Currently, chemotherapy with praziquantel (PZQ) is the only effective treatment against schistosomiasis, although the potential of drug resistance remains a concern. Hence, there is an urgent demand for new and effective strategies to combat schistosomiasis, which is the second most prevalent parasitic disease after malaria. MicroRNAs (miRNAs) are small non-coding RNAs that play pivotal regulatory roles in many organisms, including the development and sexual maturation of schistosomes. Thus, miRNAs are potential targets for treatment of schistosomiasis. Moreover, miRNAs can serve as multifunctional “nano-tools” for cross-species delivery in order to regulate host-parasite interactions. In this review, the multifunctional roles of miRNAs in the growth and development of schistosomes are discussed. The various regulatory functions of host-derived and worm-derived miRNAs on the progression of schistosomiasis are also thoroughly addressed, especially the promotional and inhibitory effects on schistosome-induced liver fibrosis. Additionally, the potential of miRNAs as biomarkers for the diagnosis and treatment of schistosomiasis is considered.

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasis

BACKGROUND

Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea.

FINDINGS

Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage.

CONCLUSION

This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

Bulinus truncatus transcriptome – a resource to enable molecular studies of snail and schistosome biology

Despite advances in high-throughput sequencing and bioinformatics, molecular investigations of snail intermediate hosts that transmit parasitic trematodes are scant. Here, we report the first transcriptome for Bulinus truncatus – a key intermediate host of Schistosoma haematobium – a blood fluke that causes urogenital schistosomiasis in humans. We assembled this transcriptome from short- and long-read RNA-sequence data. From this transcriptome, we predicted 12,998 proteins, 58% of which had orthologs in Biomphalaria glabrata – an intermediate host of Schistosoma mansoni – a blood fluke that causes hepato-intestinal schistosomiasis. We predicted that select protein groups are involved in signal transduction, cell growth and death, the immune system, environmental adaptation and/or the excretory/secretory system, suggesting roles in immune responses, pathogen defence and/or parasite-host interactions. The transcriptome of Bu. truncatus provides a useful resource to underpin future molecular investigations of this and related snail species, and its interactions with pathogens including S. haematobium. The present resource should enable comparative investigations of other molluscan hosts of socioeconomically important parasites in the future.

•

First transcriptome to represent Bulinus truncatus – a snail intermediate host of Schistosoma haematobium.

•

Select protein groups of Bu. truncatus are inferred to associate with innate immune responses against pathogens.

•

Transcriptome provides a resource for future studies of parasite-host interactions and snail-host resistance to pathogens.

Deep sequencing of small RNAs reveals the repertoire of miRNAs and piRNAs in Biomphalaria glabrata

BACKGROUND

Biomphalaria glabrata snails are widely distributed in schistosomiasis endemic areas like America and Caribe, displaying high susceptibility to infection by Schistosoma mansoni. After the availability of B. glabrata genome and transcriptome data, studies focusing on genetic markers and small non-coding RNAs have become more relevant. The small RNAs have been considered important through their ability to finely regulate the gene expression in several organisms, thus controlling the functions like cell growth, metabolism, and susceptibility/resistance to infection.

OBJECTIVE

The present study aims on identification and characterisation of the repertoire of small non-coding RNAs in B. glabrata (Bgl-small RNAs).

METHODS

By using small RNA sequencing, bioinformatics tools and quantitative reverse transcription polymerase chain reaction (RT-qPCR), we identified, characterised, and validated the presence of small RNAs in B. glabrata.

FINDINGS

89 mature miRNAs were identified and five of them were classified as Mollusk-specific. When compared to model organisms, sequences of B. glabrata miRNAs showed a high degree of conservation. In addition, several target genes were predicted for all the mature miRNAs identified. Furthermore, piRNAs were identified in the genome of B. glabrata for the first time. The B. glabrata piRNAs showed strong conservation of uridine as first nucleotide at 5’ end, besides adenine at 10th position. Our results showed that B. glabrata has diverse repertoire of circulating ncRNAs, several which might be involved in mollusk susceptibility to infection, due to their potential roles in the regulation of S. mansoni development.

MAIN CONCLUSIONS

Further studies are necessary in order to confirm the role of the Bgl-small RNAs in the parasite/host relationship thus opening new perspectives on interference of small RNAs in the organism development and susceptibility to infection.

Temporal transcriptome change of Oncomelania hupensis revealed by Schistosoma japonicum invasion

Background

The freshwater snail Oncomelania hupensis is the obligate intermediate host for Schistosoma japonicum in China. Transcriptomic examination of snail–schistosome interactions can provide valuable information of host response at physiological and immune levels.

Methods

To investigate S. japonicum-induced changes in O. hupensis gene expression, we utilized high-throughput sequencing to identify transcripts that were differentially expressed between infected snails and their uninfected controls at two key time-point, Day 7 and Day 30 after challenge. Time-series transcriptomic profiles were analyzed using R package DESeq 2, followed by GO, KEGG and (weighted gene correlation network analysis) WGCNA analysis to elucidate and identify important molecular mechanism, and subsequently understand host–parasite relationship. The identified unigenes was verified by bioinformatics and real-time PCR. Possible adaptation molecular mechanisms of O. hupensis to S. japonicum challenge were proposed.

Results

Transcriptomic analyses of O. hupensis by S. japonicum invasion yielded billion reads including 92,144 annotated transcripts. Over 5000 differentially expressed genes (DEGs) were identified by pairwise comparisons of infected libraries from two time points to uninfected libraries in O. hupensis. In total, 6032 gene ontology terms and 149 KEGG pathways were enriched. After the snails were infected with S. japonicum on Day 7 and Day 30, DEGs were shown to be involved in many key processes associated with biological regulation and innate immunity pathways. Gene expression patterns differed after exposure to S. japonicum. Using WGCNA, 16 modules were identified. Module-trait analysis identified that a module involved in RNA binding, ribosome, translation, mRNA processing, and structural constituent of ribosome were strongly associated with S. japonicum invasion. Many of the genes from enriched KEGG pathways were involved in lysosome, spliceosome and ribosome, indicating that S. japonicum invasion may activate the regulation of ribosomes and immune response to infection in O. hupensis.

Conclusions

Our analysis provided a temporally dynamic gene expression pattern of O. hupensis by S. japonicum invasion. The identification of gene candidates serves as a foundation for future investigations of S. japonicum infection. Additionally, major DEGs expression patterns and putative key regulatory pathways would provide useful information to construct gene regulatory networks between host-parasite crosstalk.

Biomphalaria glabrata immunity: Post-genome advances

The freshwater snail, Biomphalaria glabrata, is an important intermediate host in the life cycle for the human parasite Schistosoma mansoni, the causative agent of schistosomiasis. Current treatment and prevention strategies have not led to a significant decrease in disease transmission. However, the genome of B. glabrata was recently sequenced to provide additional resources to further our understanding of snail biology. This review presents an overview of recently published, post-genome studies related to the topic of snail immunity. Many of these reports expand on findings originated from the genome characterization. These novel studies include a complementary gene linkage map, analysis of the genome of the B. glabrata embryonic (Bge) cell line, as well as transcriptomic and proteomic studies looking at snail-parasite interactions and innate immune memory responses towards schistosomes. Also included are biochemical investigations on snail pheromones, neuropeptides, and attractants, as well as studies investigating the frontiers of molluscan epigenetics and cell signaling were also included. Findings support the current hypotheses on snail-parasite strain compatibility, and that snail host resistance to schistosome infection is dependent not only on genetics and expression, but on the ability to form multimeric molecular complexes in a timely and tissue-specific manner. The relevance of cell immunity is reinforced, while the importance of humoral factors, especially for secondary infections, is supported. Overall, these studies reflect an improved understanding on the diversity, specificity, and complexity of molluscan immune systems.

The role of microRNAs in the pathogenesis, grading and treatment of hepatic fibrosis in schistosomiasis

Schistosomiasis is a prevalent parasitic disease worldwide. The main pathological changes of hepatosplenic schistosomiasis are hepatic granuloma and fibrosis due to worm eggs. Portal hypertension and ascites induced by hepatic fibrosis are usually the main causes of death in patients with chronic hepatosplenic schistosomiasis. Currently, no effective vaccine exists for preventing schistosome infections. For quite a long time, praziquantel (PZQ) was widely used for the treatment of schistosomiasis and has shown benefit in treating liver fibrosis. However, drug resistance and chemical toxicity from PZQ are being increasingly reported in recent years; therefore, new and effective strategies for treating schistosomiasis-induced hepatic fibrosis are urgently needed. MicroRNA (miRNA), a non-coding RNA, has been proved to be associated with the development of many human diseases, including schistosomiasis. In this review, we present a balanced and comprehensive view of the role of miRNAs in the pathogenesis, grading, and treatment of schistosomiasis-associated hepatic fibrosis. The multiple regulatory roles of miRNAs, such as promoting or inhibiting the development of liver pathology in murine schistosomiasis are also discussed in depth. Additionally, miRNAs may serve as candidate biomarkers for diagnosing liver pathology of schistosomiasis and as novel therapeutic targets for treating schistosomiasis-associated hepatic fibrosis.

Genome-wide identification, characterisation and expression profiling of the ubiquitin-proteasome genes in Biomphalaria glabrata

BACKGROUND

Biomphalaria glabrata is the major species used for the study of schistosomiasis-related parasite-host relationships, and understanding its gene regulation may aid in this endeavor. The ubiquitin-proteasome system (UPS) performs post-translational regulation in order to maintain cellular protein homeostasis and is related to several mechanisms, including immune responses.

OBJECTIVE

The aims of this work were to identify and characterise the putative genes and proteins involved in UPS using bioinformatic tools and also their expression on different tissues of B. glabrata.

METHODS

The putative genes and proteins of UPS in B. glabrata were predicted using BLASTp and as queries reference proteins from model organism. We characterised these putative proteins using PFAM and CDD software describing the conserved domains and active sites. The phylogenetic analysis was performed using ClustalX2 and MEGA5.2. Expression evaluation was performed from 12 snail tissues using RPKM.

FINDINGS

119 sequences involved in the UPS in B. glabrata were identified, which 86 have been related to the ubiquitination pathway and 33 to proteasome. In addition, the conserved domains found were associated with the ubiquitin family, UQ_con, HECT, U-box and proteasome. The main active sites were lysine and cysteine residues. Lysines are responsible and the starting point for the formation of polyubiquitin chains, while the cysteine residues of the enzymes are responsible for binding to ubiquitin. The phylogenetic analysis showed an organised distribution between the organisms and the clades of the sequences, corresponding to the tree of life of the animals, for all groups of sequences analysed. The ubiquitin sequence was the only one with a high expression profile found in all libraries, inferring its wide range of performance.

MAIN CONCLUSIONS

Our results show the presence, conservation and expression profile of the UPS in this mollusk, providing a basis and new knowledge for other studies involving this system. Due to the importance of the UPS and B. glabrata, this work may influence the search for new methodologies for the control of schistosomiasis.

Proteomic Analysis of Biomphalaria glabrata Hemocytes During in vitro Encapsulation of Schistosoma mansoni Sporocysts

Circulating hemocytes of the snail Biomphalaria glabrata, a major intermediate host for the blood fluke Schistosoma mansoni, represent the primary immune effector cells comprising the host's internal defense system. Within hours of miracidial entry into resistant B. glabrata strains, hemocytes infiltrate around developing sporocysts forming multi-layered cellular capsules that results in larval death, typically within 24–48 h post-infection. Using an in vitro model of hemocyte-sporocyst encapsulation that recapitulates in vivo events, we conducted a comparative proteomic analysis on the responses of hemocytes from inbred B. glabrata strains during the encapsulation of S. mansoni primary sporocysts. This was accomplished by a combination of Laser-capture microdissection (LCM) to isolate sections of hemocyte capsules both in the presence and absence of sporocysts, in conjunction with mass spectrometric analyses to establish protein expression profiles. Comparison of susceptible NMRI snail hemocytes in the presence and absence of sporocysts revealed a dramatic downregulation of proteins in during larval encapsulation, especially those involved in protein/CHO metabolism, immune-related, redox and signaling pathways. One of 4 upregulated proteins was arginase, competitor of nitric oxide synthetase and inhibitor of larval-killing NO production. By contrast, when compared to control capsules, sporocyst-encapsulating hemocytes of resistant BS-90 B. glabrata exhibited a more balanced profile with enhanced expression of shared proteins involved in protein synthesis/processing, immunity, and redox, and unique expression of anti-microbial/anti-parasite proteins. A final comparison of NMRI and BS-90 host hemocyte responses to co-cultured sporocysts demonstrated a decrease or downregulation of 77% of shared proteins by NMRI cells during encapsulation compared to those of the BS-90 strain, including lipopolysaccharide-binding protein, thioredoxin reductase 1 and hemoglobins 1 and 2. Overall, using this in vitro model, results of our proteomic analyses demonstrate striking differences in proteins expressed by susceptible NMRI and resistant BS-90 snail hemocytes to S. mansoni sporocysts during active encapsulation, with NMRI hemocytes exhibiting extensive downregulation of protein expression and a lower level of constitutively expressed immune-relevant proteins (e.g., FREP2) compared to BS-90. Our data suggest that snail strain differences in hemocyte protein expression during the encapsulation process account for observed differences in their cytotoxic capacity to interact with and kill sporocysts.