Biomphalaria glabrata immunity: Post-genome advances

SlugAtlas, a histological and 3D online resource of the land slugs Deroceras laeve and Ambigolimax valentianus

Due to their distinctive anatomical characteristics, land slugs are appealing research subjects from a variety of angles, including stem cell biology, regeneration, embryonic development, allometry, and neurophysiology. Here we present SlugAtlas, an anatomical and histological online resource of the land slugs Deroceras laeve and Ambigolimax valentianus. The atlas is composed of series of histological sections on the horizontal, sagittal, and transversal planes for both species with 3D viewing tools of their major organs. The atlas was used in this work for a comparative analysis of the major organs and tissues of these two species. We provide a comprehensive histological description of them and an explanation of novel findings of unique features of their anatomy. For D. laeve, we additionally studied its ability for degrowth and regrowth, a feature that characterizes animals with high regeneration potential and obtained initial results on the study of the regeneration of its tail. SlugAtlas is a resource that is also built to accommodate future growth and, along with the experimental techniques that we have developed, will provide the foundation for research projects in a variety of biological domains.

Lipids from a snail host regulate the multicellular behavior of a predator of parasitic schistosomes

Transmission of vector-borne diseases can be slowed by symbionts within the secondary hosts that spread disease. Snails spread schistosomiasis, and the snail symbiont Capsaspora owczarzaki kills schistosome larvae. In studying how Capsaspora colonizes its host snail, we discovered that Capsaspora responded to its host by forming multicellular aggregates. We elucidated the chemical cue for aggregation: hemolymph phosphatidylcholines (PCs). Furthermore, we uncovered that Capsaspora cells aggregate to different degrees in sera from different host snails-and these responses correlate with serum concentrations of PCs. Therefore, Capsaspora senses a host factor that can indicate the identity and physiological state of its host. Since cellular aggregation controls microbial motility, feeding, and immune evasion, this response within host tissue may be important for colonization. If so, snail serum PC and Capsaspora aggregation will be molecular and cellular markers to discern which conditions will favor the colonization of snails (and potential exclusion of schistosomes) by Capsaspora.

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.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Host lipids regulate multicellular behavior of a predator of a human pathogen

As symbionts of animals, microbial eukaryotes benefit and harm their hosts in myriad ways. A model microeukaryote (Capsaspora owczarzaki) is a symbiont of Biomphalaria glabrata snails and may prevent transmission of parasitic schistosomes from snails to humans. However, it is unclear which host factors determine Capsaspora's ability to colonize snails. Here, we discovered that Capsaspora forms multicellular aggregates when exposed to snail hemolymph. We identified a molecular cue for aggregation: a hemolymph-derived phosphatidylcholine, which becomes elevated in schistosome-infected snails. Therefore, Capsaspora aggregation may be a response to the physiological state of its host, and it may determine its ability to colonize snails and exclude parasitic schistosomes. Furthermore, Capsaspora is an evolutionary model organism whose aggregation may be ancestral to animals. This discovery, that a prevalent lipid induces Capsaspora multicellularity, suggests that this aggregation phenotype may be ancient. Additionally, the specific lipid will be a useful tool for further aggregation studies.

Assessing the microbiota of the snail intermediate host of trematodes, Galba truncatula

BACKGROUND

The microbiome is known to play key roles in health and disease, including host susceptibility to parasite infections. The freshwater snail Galba truncatula is the intermediate host for many trematode species, including the liver and rumen flukes Fasciola hepatica and Calicophoron daubneyi, respectively. The snail-parasite system has previously been investigated. However, the specific interaction between the snail-associated microbiota and intra-snail developmental stages of trematodes has yet to be explored.

METHODS

Galba truncatula snails were collected from farms in Northern Ireland and trematode infection was diagnosed using PCR. High-throughput sequencing analysis of the bacterial 16S ribosomal DNA V3-V4 hypervariable regions was subsequently applied to characterise the microbiota of both uninfected and infected snails.

RESULTS

We first showed that the snail harboured microbiota that was distinct for its environment. The microbiota of infected snails was found to differ significantly from that of uninfected snails. In particular, the bacterial genera Mycoplasma and Methylotenera were significantly more abundant in infected snails, while genera Sphingomonas and Nocardioides were predominantly associated with uninfected snails.

CONCLUSION

These findings pave the way to future studies on the functional roles of bacteria in host-parasite relationships.

Immuno-molecular profile for Biomphalaria glabrata/Schistosoma mansoni interaction

The complex innate immune defense of Biomphalaria glabrata, the intermediate host of Schistosoma mansoni, governs the successful development of the intramolluscan stages of the parasite. The interaction between the snail and the parasite involves a complex immune molecular crosstalk between several parasite antigens and the snail immune recognition receptors, evoking different signals and effector molecules. This work seeks to discuss the immune-related molecules that influence compatibility in Biomphalaria glabrata/Schistosoma mansoni interaction and the differential expression of these molecules between resistant and susceptible snails. It also includes the current understanding of the immune molecular determinants that govern the compatibility in sympatric and allopatric interactions, and the expression of these molecules after immune priming and the secondary immune response. Herein, the differences in the immune-related molecules in the interaction of other Biomphalaria species with Schistosoma mansoni compared to the Biomphalaria glabrata model snail are highlighted. Understanding the diverse immune molecular determinants in the snail/schistosome interaction can lead to alternative control strategies for schistosomiasis.

Exploring the immune interactions between Oncomelania hupensis and Schistosoma japonicum, with a cross-comparison of immunological research progress in other intermediate host snails

Schistosomiasis, the second largest parasitic disease in the world after malaria, poses a significant threat to human health and causes public health issues. The disease primarily affects populations in economically underdeveloped tropical regions, earning it the title of "neglected tropical disease". Schistosomiasis is difficult to eradicate globally if medication alone is used. One of the essential elements of thorough schistosomiasis prevention and control is the management and disruption of the life cycle of intermediate host snails. The key approach to controlling the transmission of schistosomiasis is to control the intermediate hosts of the schistosome to disrupt its life cycle. We believe that approaching it from the perspective of the intermediate host's immunity could be an environmentally friendly and potentially effective method. Currently, globally significant intermediate host snails for schistosomes include Oncomelania hupensis, Biomphalaria glabrata, and Bulinus truncatus. The immune interaction research between B. glabrata and Schistosoma mansoni has a history of several decades, and the complete genome sequencing of both B. glabrata and B. truncatus has been accomplished. We have summarized the immune-related factors and research progress primarily studied in B. glabrata and B. truncatus and compared them with several humoral immune factors that O. hupensis research focuses on: macrophage migration inhibitory factor (MIF), Toll-like receptors (TLRs), and thioredoxin (Trx). We believe that continued exploration of the immune interactions between O. hupensis and Schistosoma japonicum is valuable. This comparative analysis can provide some direction and clues for further in-depth research. Comparative immunological studies between them not only expand our understanding of the immune defense responses of snails that act as intermediaries for schistosomes but also facilitate the development of more comprehensive and integrated strategies for schistosomiasis prevention and control. Furthermore, it offers an excellent opportunity to study the immune system of gastropods and their co-evolution with pathogenic organisms.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Status Quo and Future Perspectives of Molecular and Genomic Studies on the Genus Biomphalaria-The Intermediate Snail Host of Schistosoma mansoni

Schistosomiasis, or also generally known as bilharzia or snail fever, is a parasitic disease that is caused by trematode flatworms of the genus Schistosoma. It is considered by the World Health Organisation as the second most prevalent parasitic disease after malaria and affects more than 230 million people in over 70 countries. People are infected via a variety of activities ranging from agricultural, domestic, occupational to recreational activities, where the freshwater snails Biomphalaria release Schistosoma cercariae larvae that penetrate the skin of humans when exposed in water. Understanding the biology of the intermediate host snail Biomphalaria is thus important to reveal the potential spread of schistosomiasis. In this article, we present an overview of the latest molecular studies focused on the snail Biomphalaria, including its ecology, evolution, and immune response; and propose using genomics as a foundation to further understand and control this disease vector and thus the transmission of schistosomiasis.

Using freshwater snail Biomphalaria glabrata (Say, 1818) as a biological model for ecotoxicology studies: a systematic review

Over time, a growing increase in human pollutants in the aquatic environment has been observed. The global presence of residues in water bodies reinforces the need to develop improved methods to detect them and evaluate their ecotoxicological effects in aquatic environments. Thus, this study aimed to present the main assays using Biomphalaria glabrata as a biological model for ecotoxicological studies. We performed a systematic literature review with data published up to June 2022 on the Web of Science, SCOPUS, Science Direct, PubMed, and SciELO databases. Thirty studies were selected for this review after screening. Biomphalaria glabrata has been studied as an ecotoxicological model for different substances through toxicity, embryotoxicity, cytotoxicity, genotoxicity, and bioaccumulation assays. Studies evaluating the impact of B. glabrata exposure to several substances have reported effects on their offspring, as well as toxicity and behavioral and reproductive effects. This review presents various assays using B. glabrata as a biological model for ecotoxicological studies. The use of a representative species of ecosystems from tropical regions is a necessary tool for tropical environmental monitoring. It was observed that the freshwater snail B. glabrata was effective for the evaluation of the ecotoxicity of several types of chemical substances, but further studies are needed to standardize the model.

Schiff bases complexed with iron and their relation with the life cycle and infection by Schistosoma mansoni

Introduction

The trematode Schistosoma mansoni causes schistosomiasis, and this parasite's life cycle depends on the mollusk Biomphalaria glabrata. The most effective treatment for infected people is administering a single dose of Praziquantel. However, there are naturally resistant to treatment. This work has developed, considering this parasite's complex life cycle.

Methods

The synthetics compound were evaluated: i) during the infection of B. glabrata, ii) during the infection of BALB/c mice, and iii) during the treatment of mice infected with S. mansoni.

Results and Discussion

For the first objective, snails infected with miracidia treated with compounds C1 and C3 at concentrations of 25% IC50 and 50% IC50, after 80 days of infection, released fewer cercariae than the infected group without treatment. For the second objective, compounds C1 and C3 did not show significant results in the infected group without treatment. For the third objective, the mice treated with C3 and C1 reduced the global and differential cell count. The results suggest that although the evaluated compounds do not present schistosomicidal properties when placed in cercariae suspension, they can stimulate an immune reaction in snails and decrease mice's inflammatory response. In general, we can conclude that compound C1 and C3 has an anti-schistosomicidal effect both in the larval phase (miracidia) and in the adult form of the parasite.

Single cell RNA sequencing reveals hemocyte heterogeneity in Biomphalaria glabrata: Plasticity over diversity

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.

Molecular characterization of thioester-containing proteins in Biomphalaria glabrata and their differential gene expression upon Schistosoma mansoni exposure

Schistosomiasis is a disease caused by trematode parasites of the genus Schistosoma that affects approximately 200 million people worldwide. Schistosomiasis has been a persistent problem in endemic areas as there is no vaccine available, currently used anti-helmintic medications do not prevent reinfection, and most concerning, drug resistance has been documented in laboratory and field isolates. Thus, alternative approaches to curtail this human disease are warranted. Understanding the immunobiology of the obligate intermediate host of these parasites, which include the freshwater snail Biomphalaria glabrata, may facilitate the development of novel methods to stop or reduce transmission to humans. Molecules from the thioester-containing protein (TEP) superfamily have been shown to be involved in immunological functions in many animals including corals and humans. In this study we identified, characterized, and compared TEP transcripts and their expression upon S. mansoni exposure in resistant and susceptible strains of B. glabrata snails. Results showed the expression of 11 unique TEPs in B. glabrata snails. These transcripts present high sequence identity at the nucleotide and putative amino acid levels between susceptible and resistant strains. Further analysis revealed differences in several TEPs’ constitutive expression levels between resistant and susceptible snail strains, with C3-1, C3-3, and CD109 having higher constitutive expression levels in the resistant (BS90) strain, whereas C3-2 and TEP-1 showed higher constitutive expression levels in the susceptible (NMRI) strain. Furthermore, TEP-specific response to S. mansoni miracidia exposure reiterated their differential expression, with resistant snails upregulating the expression of both TEP-4 and TEP-3 at 2 h and 48 h post-exposure, respectively. Further understanding the diverse TEP genes and their functions in invertebrate animal vectors will not only expand our knowledge in regard to this ancient family of immune proteins, but also offer the opportunity to identify novel molecular targets that could aid in the efforts to develop control methods to reduce schistosomiasis transmission.

Transcriptome profiling based on larvae at different time points after hatching provides a core set of gene resource for understanding the immune response mechanisms of the egg-protecting behavior against Vibrio anguillarum infection in Amphioctopus fangsiao

Mollusks have recently received increasing attention because of their unique immune systems. Mollusks such as Amphioctopus fangsiao are economically important cephalopods, and the effects of their egg-protecting behavior on the larval immune response are unclear. Meanwhile, little research has been done on the resistance response of cephalopod larvae infected with pathogenic bacteria such as Vibrio anguillarum. In this study, V. anguillarum was used to infect the primary hatching A. fangsiao larvae under different egg-protecting behaviors for 24 h, and a total of 7156 differentially expressed genes (DEGs) were identified at four time points after hatching based on transcriptome analysis. GO and KEGG enrichment analyses showed that multiple immune-related GO terms and KEGG signaling pathways were enriched. Protein-protein interaction networks (PPI networks) were used to search functional relationships between immune-related DEGs. Finally, 20 hub genes related to multiple gene functions or involved in multiple signaling pathways were identified, and their accuracy was verified using quantitative RT-PCR. PPI networks were first used to study the effects A. fangsiao larvae after infection with V. anguillarum under different egg-protecting behaviors. The results provide significant genetic resources for exploring invertebrate larval immune processes. The data lays a foundation for further study the immune response mechanisms for invertebrates after infection.

An Overview of Transcriptional Responses of Schistosome-Susceptible (M line) or -Resistant (BS-90) Biomphalaria glabrata Exposed or Not to Schistosoma mansoni Infection

Background

We seek to provide a comprehensive overview of transcriptomics responses of immune-related features of the gastropod Biomphalaria glabrata (Bg) following exposure to Schistosoma mansoni (Sm), a trematode causing human schistosomiasis. Responses of schistosome-susceptible (M line, or SUS) and -resistant (BS-90, or RES) Bg strains were characterized following exposure to Sm for 0.5, 2, 8 or 40 days post-exposure (dpe).

Methods

RNA-Seq and differential expression analysis were undertaken on 56 snails from 14 groups. We considered 7 response categories: 1) constitutive resistance factors; 2) constitutive susceptibility factors; 3) generalized stress responses; 4) induced resistance factors; 5) resistance factors suppressed in SUS snails; 6) suppressed/manipulated factors in SUS snails; and 7) tolerance responses in SUS snails. We also undertook a gene co-expression network analysis. Results from prior studies identifying schistosome resistance/susceptibility factors were examined relative to our findings.

Results

A total of 792 million paired-end reads representing 91.2% of the estimated 31,985 genes in the Bg genome were detected and results for the 7 categories compiled and highlighted. For both RES and SUS snails, a single most supported network of genes with highly correlated expression was found.

Conclusions

1) Several constitutive differences in gene expression between SUS and RES snails were noted, the majority over-represented in RES; 2) There was little indication of a generalized stress response shared by SUS and RES snails at 0.5 or 2 dpe; 3) RES snails mounted a strong, multi-faceted response by 0.5 dpe that carried over to 2 dpe; 4) The most notable SUS responses were at 40 dpe, in snails shedding cercariae, when numerous features were either strongly down-regulated indicative of physiological distress or parasite manipulation, or up-regulated, suggestive of tolerance or survival-promoting effects; 5) Of 55 genes previously identified in genome wide mapping studies, 29 (52.7%) were responsive to Sm, as were many familiar resistance-associated genes (41.0%) identified by other means; 6) Both network analysis and remarkably specific patterns of expression of lectins and G protein-coupled receptors in categories 4, 6 and 7 were indicative of orchestrated responses of different suites of genes in SUS or RES snails following exposure to Sm.

Nuclear genome of Bulinus truncatus, an intermediate host of the carcinogenic human blood fluke Schistosoma haematobium

Some snails act as intermediate hosts (vectors) for parasitic flatworms (flukes) that cause neglected tropical diseases, such as schistosomiases. Schistosoma haematobium is a blood fluke that causes urogenital schistosomiasis and induces bladder cancer and increased risk of HIV infection. Understanding the molecular biology of the snail and its relationship with the parasite could guide development of an intervention approach that interrupts transmission. Here, we define the genome for a key intermediate host of S. haematobium-called Bulinus truncatus-and explore protein groups inferred to play an integral role in the snail's biology and its relationship with the schistosome parasite. Bu. truncatus shared many orthologous protein groups with Biomphalaria glabrata-the key snail vector for S. mansoni which causes hepatointestinal schistosomiasis in people. Conspicuous were expansions in signalling and membrane trafficking proteins, peptidases and their inhibitors as well as gene families linked to immune response regulation, such as a large repertoire of lectin-like molecules. This work provides a sound basis for further studies of snail-parasite interactions in the search for targets to block schistosomiasis transmission.

Comparative proteomics and transcriptomics illustrate the allograft-induced stress response in the pearl oyster (Pinctada fucata martensii)

Implantation of a spherical nucleus into a recipient oyster is a critical step in artificial pearl production. However, the molecular mechanisms underlying the response of the pearl oyster to this operation are poorly understood. In this research, we used transcriptomic and proteomic analyses to examine allograft-induced changes in gene/protein expression patterns in Pinctada fucata martensii 12 h after nucleus implantation. Transcriptome analysis identified 688 differential expression genes (DEGs) (FDR<0.01 and |fold change) ＞ 2). Using a 1.2-fold increase or decrease in protein expression as a benchmark for differentially expressed proteins (DEPs), 108 DEPs were reliably quantified, including 71 up-regulated proteins (DUPs) and 37 down-regulated proteins (DDPs). Further analysis revealed that the GO terms, including "cellular process", "biological regulation" and "metabolic process" were considerably enriched. In addition, the transcriptomics analysis showed that "Neuroactive ligand-receptor interaction", "NF-kappa B signaling pathway", "MAPK signaling pathway", "PI3K-Akt signaling pathway', "Toll-like receptor signaling pathway", and "Notch signaling pathway" were significantly enriched in DEGs. The proteomics analysis showed that "ECM-receptor interaction", "Human papillomavirus infection", and "PI3K-Akt signaling pathway" were significantly enriched in DEPs. The results indicate that these functions could play an important role in response to pear oyster stress at nucleus implantation. To assess the potential relevance of quantitative information between mRNA and proteins, using Ward's hierarchical clustering analysis clustered the protein/gene expression patterns across the experimental and control samples into six groups. To investigate the biological processes associated with the protein in each cluster, we identified the significantly enriched GO terms and KEGG pathways in the proteins in each cluster. Gene set enrichment analysis (GSEA) was used to reveal the potential protein or transcription pathways associated with the response to nuclear implantation. Thus, the study of P. f. martensii is essential to enhance our understanding of the molecular mechanisms involved in pearl biosynthesis and the biology of bivalve molluscs.

Single-cell RNA-seq profiling of individual Biomphalaria glabrata immune cells with a focus on immunologically relevant transcripts

The immune cells of the snail Biomphalaria glabrata are classified into hyalinocyte and granulocyte subtypes. Both subtypes are essential for the proper functioning of the snail immune response, which we understand best within the context of how it responds to challenge with the human parasite Schistosoma mansoni. Granulocytes are adherent phagocytic cells that possess conspicuous granules within the cell cytoplasm. Hyalinocytes, on the other hand, are predominantly non-adherent and are known to produce a handful of anti-S. mansoni immune effectors. While our understanding of these cells has progressed, an in-depth comparison of the functional capabilities of each type of immune cell has yet to be undertaken. Here, we present the results of a single-cell RNA-seq study in which single granulocytes and hyalinocytes from S. mansoni-susceptible M-line B. glabrata and S. mansoni-resistant BS-90 B. glabrata are compared without immune stimulation. This transcriptomic analysis supports a role for the hyalinocytes as producers of immune effectors such as biomphalysin and thioester-containing proteins. It suggests that granulocytes are primarily responsible for producing fibrinogen-related proteins and are armed with various pattern-recognition receptors such as toll-like receptors with a confirmed role in the anti-S. mansoni immune response. This analysis also confirms that the granulocytes and hyalinocytes of BS-90 snails are generally more immunologically prepared than their M-line counterparts. As the first single-cell analysis of the transcriptional profiles of B. glabrata immune cells, this study provides crucial context for understanding the B. glabrata immune response. It sets the stage for future investigations into how each immune cell subtype differs in its response to various immunological threats.

A proteomics evaluation of the primary and secondary immune response of Biomphalaria straminea challenged by Schistosoma mansoni

Biomphalaria spp. snails are intermediary hosts of Schistosoma mansoni, etiologic agent of intestinal schistosomiasis, one of the most important neglected tropical diseases. Biomphalaria straminea is an important intermediary host that possess a different phenotype to parasite infection but shows a large geographic distribution and high capacity of new ecologic niche invasion. Our purpose was to characterize for the first time the differentially expressed proteome in B. straminea during two times intervals after primary and secondary exposure to S. mansoni. The hemolymph was collected at 1 and 15 days after primary and secondary exposure of snails to the parasite. Total proteins were extracted and digested with trypsin. LC-MS/MS label-free quantification was performed and analyzed using Maxquant and Perseus software. Proteins were identified and annotated using Blast2GO tools. After 1 day of exposure, most of upregulated proteins are hemoglobin type 2, C and H type lectins, molecules related to cell adhesion, and response to oxidative stress. After 15 days, we found a similar pattern of upregulated proteins but some fibrinogen-related proteins (FREPs) and TEPs homologs were downregulated. Regarding the differentially expressed proteins during secondary response, the principal immune-related proteins upregulated were C and H type lectins, cellular adhesion molecules, biomphalysin, and FREP3. We noted a several upregulated biological processes during both responses that could be the one of the key points of efficacy in the immune response to parasite. Our data suggests different immune mechanisms used by B. straminea snails challenged with S. mansoni.

PIWI silencing mechanism involving the retrotransposon nimbus orchestrates resistance to infection with Schistosoma mansoni in the snail vector, Biomphalaria glabrata

Background

Schistosomiasis remains widespread in many regions despite efforts at its elimination. By examining changes in the transcriptome at the host-pathogen interface in the snail Biomphalaria glabrata and the blood fluke Schistosoma mansoni, we previously demonstrated that an early stress response in juvenile snails, manifested by induction of heat shock protein 70 (Hsp 70) and Hsp 90 and of the reverse transcriptase (RT) domain of the B. glabrata non-LTR- retrotransposon, nimbus, were critical for B. glabrata susceptibility to S. mansoni. Subsequently, juvenile B. glabrata BS-90 snails, resistant to S. mansoni at 25°C become susceptible by the F2 generation when maintained at 32°C, indicating an epigenetic response.

Methodology/Principal findings

To better understand this plasticity in susceptibility of the BS-90 snail, mRNA sequences were examined from S. mansoni exposed juvenile BS-90 snails cultured either at 25°C (non-permissive temperature) or 32°C (permissive). Comparative analysis of transcriptomes from snails cultured at the non-permissive and permissive temperatures revealed that whereas stress related transcripts dominated the transcriptome of susceptible BS-90 juvenile snails at 32°C, transcripts encoding proteins with a role in epigenetics, such as PIWI (BgPiwi), chromobox protein homolog 1 (BgCBx1), histone acetyltransferase (BgHAT), histone deacetylase (BgHDAC) and metallotransferase (BgMT) were highly expressed in those cultured at 25°C. To identify robust candidate transcripts that will underscore the anti-schistosome phenotype in B. glabrata, further validation of the differential expression of the above transcripts was performed by using the resistant BS-90 (25°C) and the BBO2 susceptible snail stock whose genome has now been sequenced and represents an invaluable resource for molecular studies in B. glabrata. A role for BgPiwi in B. glabrata susceptibility to S. mansoni, was further examined by using siRNA corresponding to the BgPiwi encoding transcript to suppress expression of BgPiwi, rendering the resistant BS-90 juvenile snail susceptible to infection at 25°C. Given transposon silencing activity of PIWI as a facet of its role as guardian of the integrity of the genome, we examined the expression of the nimbus RT encoding transcript at 120 min after infection of resistant BS90 piwi-siRNA treated snails. We observed that nimbus RT was upregulated, indicating that modulation of the transcription of the nimbus RT was associated with susceptibility to S. mansoni in BgPiwi-siRNA treated BS-90 snails. Furthermore, treatment of susceptible BBO2 snails with the RT inhibitor lamivudine, before exposure to S. mansoni, blocked S. mansoni infection concurrent with downregulation of the nimbus RT transcript and upregulation of the BgPiwi encoding transcript in the lamivudine-treated, schistosome-exposed susceptible snails.

Conclusions and significance

These findings support a role for the interplay of BgPiwi and nimbus in the epigenetic modulation of plasticity of resistance/susceptibility in the snail-schistosome relationship.

Progress is being made to eliminate schistosomiasis, a tropical disease that remains endemic in the tropics and neotropics. In 2020, WHO proposed controlling the snail population as part of a strategy toward reducing schistosomiasis, a vector borne disease, by 2025. The life cycle of the causative parasite is, however, complex and in the absence of vaccines, new drugs, and access to clean water and sanitation, reduction of schistosomiasis will remain elusive. To break the parasite’s life cycle during the snail stage of its development, a better understanding of the molecular basis of how schistosomes survive, or not, in the snail is required. By examining changes in the transcriptome at the host-pathogen interface in the snail Biomphalaria glabrata and Schistosoma mansoni, we showed that early stress response, manifested by the induction of Heat Shock Proteins (Hsps) and the RT domain of the non-LTR retrotransposon, nimbus, were critical for snail susceptibility. Subsequently, juvenile B. glabrata BS-90 snails, resistant to S. mansoni at 25°C were observed to become susceptible by the F2 generation when maintained at 32°C, indicating an epigenetic response. This study confirms these earlier results and shows an interplay between PIWI and nimbus in the anti-schistosome response in the snail host.

Current understanding of immune priming phenomena in insects

It may seem that the most important issues related to insect immunity have already been described. However, novel phenomena observed in recent years shed new light on the understanding of the immune response in insects.The adaptive abilities of insects helped them to populate all ecological land niches.One important adaptive ability of insects that facilitates their success is the plasticity of their immune system. Although they only have innate immune mechanisms, insects can increase their resistance after the first encounter with the pathogen. In recent years, this phenomenon,namedimmunepriming, has become a "hot topic" in immunobiology.Priming can occur within or across generations. In the first case, the resistance of a given individual can increase after surviving a previous infection. Transstadial immune priming occurs when infection takes place at one of the initial developmental stages and increased resistance is observed at the pupal or imago stages. Priming across generations (transgenerationalimmune priming, TGIP) relies on the increased resistance of the offspring when one or both parents are infected during their lifetime.Despite the attention that immune priming has received, basic questions remain to be answered, such as regulation of immune priming at the molecular level. Research indicates that pathogen recognition receptors (PRRs) can be involved in the priming phenomenon. Recent studies have highlighted the special role of microRNAs and epigenetics, which can influence expression of genes that can be transmitted through generations although they are not encoded in the nucleotide sequence. Considerable amounts of research are required to fully understand the mechanisms that regulate priming phenomena. The aim of our work is to analyse thoroughly the most important information on immune priming in insects and help raise pertinent questions such that a greater understanding of this phenomenon can be obtained in the future.

Schistosomiasis Drug Discovery in the Era of Automation and Artificial Intelligence

Schistosomiasis is a parasitic disease caused by trematode worms of the genus Schistosoma and affects over 200 million people worldwide. The control and treatment of this neglected tropical disease is based on a single drug, praziquantel, which raises concerns about the development of drug resistance. This, and the lack of efficacy of praziquantel against juvenile worms, highlights the urgency for new antischistosomal therapies. In this review we focus on innovative approaches to the identification of antischistosomal drug candidates, including the use of automated assays, fragment-based screening, computer-aided and artificial intelligence-based computational methods. We highlight the current developments that may contribute to optimizing research outputs and lead to more effective drugs for this highly prevalent disease, in a more cost-effective drug discovery endeavor.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Mobilizing molluscan models and genomes in biology

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

New Insights Into Biomphalysin Gene Family Diversification in the Vector Snail Biomphalaria glabrata

Aerolysins initially characterized as virulence factors in bacteria are increasingly found in massive genome and transcriptome sequencing data from metazoans. Horizontal gene transfer has been demonstrated as the main way of aerolysin-related toxins acquisition in metazoans. However, only few studies have focused on their potential biological functions in such organisms. Herein, we present an extensive characterization of a multigene family encoding aerolysins - named biomphalysin - in Biomphalaria glabrata snail, the intermediate host of the trematode Schistosoma mansoni. Our results highlight that duplication and domestication of an acquired bacterial toxin gene in the snail genome result in the acquisition of a novel and diversified toxin family. Twenty-three biomphalysin genes were identified. All are expressed and exhibited a tissue-specific expression pattern. An in silico structural analysis was performed to highlight the central role played by two distinct domains i) a large lobe involved in the lytic function of these snail toxins which constrained their evolution and ii) a small lobe which is structurally variable between biomphalysin toxins and that matched to various functional domains involved in moiety recognition of targets cells. A functional approach suggests that the repertoire of biomphalysins that bind to pathogens, depends on the type of pathogen encountered. These results underline a neo-and sub-functionalization of the biomphalysin toxins, which have the potential to increase the range of effectors in the snail’s immune arsenal.

Hemocyte siRNA uptake is increased by 5' cholesterol-TEG addition in Biomphalaria glabrata, snail vector of schistosome

Biomphalaria glabrata is one of the snail intermediate hosts of Schistosoma mansoni, the causative agent of intestinal schistosomiasis disease. Numerous molecular studies using comparative approaches between susceptible and resistant snails to S. mansoni infection have helped identify numerous snail key candidates supporting such susceptible/resistant status. The functional approach using RNA interference (RNAi) remains crucial to validate the function of such candidates. CRISPR-Cas systems are still under development in many laboratories, and RNA interference remains the best tool to study B. glabrata snail genetics. Herein, we describe the use of modified small interfering RNA (siRNA) molecules to enhance cell delivery, especially into hemocytes, the snail immune cells. Modification of siRNA with 5′ Cholesteryl TriEthylene Glycol (Chol-TEG) promotes cellular uptake by hemocytes, nearly eightfold over that of unmodified siRNA. FACS analysis reveals that more than 50% of hemocytes have internalized Chol-TEG siRNA conjugated to Cy3 fluorophores, 2 hours only after in vivo injection into snails. Chol-TEG siRNA targeting BgTEP1 (ThioEster-containing Protein), a parasite binding protein, reduced BgTEP1 transcript expression by 70–80% compared to control. The level of BgTEP1 protein secreted in the hemolymph was also decreased. However, despite the BgTEP1 knock-down at both RNA and protein levels, snail compatibility with its sympatric parasite is not affected suggesting functional redundancy among the BgTEP genes family in snail-schistosoma interaction.

RNA-seq: The early response of the snail Physella acuta to the digenetic trematode Echinostoma paraensei

To analyze the response of the snail Physella acuta to Echinostoma paraensei, a compatible digenetic trematode, Illumina RNA-seq data were collected from snails with early infection (5 snails at 2 days post-exposure [DPE]) and established infection (4 snails, 8 DPE), and 7 control (unexposed) snails. A reference transcriptome (325,563 transcripts, including 98% of eukaryotic universal single-copy orthologs; BUSCO) and a draft P. acuta genome (employing available genomic Illumina reads; 799,945 scaffolds, includes 88% BUSCO genes) were assembled to guide RNA-seq analyses. Parasite exposure of P. acuta led to 10,195 differentially expressed (DE) genes at 2 DPE and 8,876 DE genes at 8 DPE with only 18% of up-regulated and 22% of down-regulated sequences shared between these time points. Gene ontology (GO) analysis yielded functional annotation of only 1.2% of DE genes but did not indicate major changes in biological activities of P. acuta between 2 and 8 DPE. Increased insights were achieved by analysis of expression profiles of 460 immune-relevant DE transcripts, identified by BLAST and InterProScan. Physella acuta has expanded gene families that encode immune-relevant domains, including CD109/TEP, GTPase IMAP, Limulus agglutination factor (dermatopontin), FReD (≥82 sequences with fibrinogen-related domains), and transcripts that combine C-type lectin (C-LECT) and C1q domains, novel among metazoa. Notably, P. acuta expressed sequences from these immune gene families at all time points, but the assemblages of unique transcripts from particular immune gene families differed between 2 and 8 DPE. The shift in profiles of DE immune genes, from early exposure to parasite establishment, suggests that compatible P. acuta initially respond to infection but switch to express immune genes that likely are less effective against E. paraensei but counter other types of (opportunistic) pathogens and parasites. We propose that the latter expression profile is part of an extended phenotype of E. paraensei, imposed upon P. acuta through parasite manipulation of the host, following successful parasite establishment in the snail after 2 DPE.

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.

The pond snail Lymnaea stagnalis

The freshwater snail Lymnaea stagnalis has a long research history, but only relatively recently has it emerged as an attractive model organism to study molecular mechanisms in the areas of developmental biology and translational medicine such as learning/memory and neurodegenerative diseases. The species has the advantage of being a hermaphrodite and can both cross- and self-mate, which greatly facilitates genetic approaches. The establishment of body-handedness, or chiromorphogenesis, is a major topic of study, since chirality is evident in the shell coiling. Chirality is maternally inherited, and only recently a gene-editing approach identified the actin-related gene Lsdia1 as the key handedness determinant. This short article reviews the natural habitat, life cycle, major research questions and interests, and experimental approaches.

The hemolymph of Biomphalaria snail vectors of schistosomiasis supports a diverse microbiome

The microbiome - the microorganism community that is found on or within an organism's body - is increasingly recognized to shape many aspects of its host biology and is a key determinant of health and disease. Microbiomes modulate the capacity of insect disease vectors (mosquitoes, tsetse flies, sandflies) to transmit parasites and disease. We investigate the diversity and abundance of microorganisms within the hemolymph (i.e. blood) of Biomphalaria snails, the intermediate host for Schistosoma mansoni, using Illumina MiSeq sequencing of the bacterial 16S V4 rDNA. We sampled hemolymph from five snails from six different laboratory populations of B. glabrata and one population of B. alexandrina. We observed 279.84 ± 0.79 amplicon sequence variants per snail. There were significant differences in microbiome composition at the level of individual snails, snail populations and species. Snail microbiomes were dominated by Proteobacteria and Bacteroidetes while water microbiomes from snail tank were dominated by Actinobacteria. We investigated the absolute bacterial load using qPCR: hemolymph samples contained 2784 ± 339 bacteria/μl. We speculate that the microbiome may represent a critical, but unexplored intermediary in the snail-schistosome interaction as hemolymph is in very close contact with the parasite at each step of its development.

Toll-like signaling pathway in the transcriptome of Littorina littorea

The transcriptome of the caenogastropod mollusk Littorina littorea was scanned for the presence of sequences encoding Toll-like receptors (TLRs) and corresponding proteins involved in downstream TLR signaling pathway. In the transcriptomic snapshots of hemocytes and kidney tissues, 45 complete TLRs encoded by 35 genes were identified. Out of the 59 non-TLR molecules involved in a canonical TLR signaling pathway, 35 genes were classified as homologous and could be placed within the TLR-mediated MyD88-and MAPK-dependent circuitries. No reference vertebrate adapters TIRAP, TRIF and TRAM were identified in the transcriptome. The results of RNA-seq experiments with an immune challenge (rediae of the digenean Himasthla elongata) indicate that four TLRs (LlTLR1, 3, 5 and 8) and a set of upregulated genes involved in signal transduction (LlMyd88, LlTNFα, LlCASP8, LlFADD, LlNFKBIA (IkBα), LlIRAK1, LlSTAT1, LlMAPK14 (P38), LlMAP2K1 (MEK1/2), LlIRF3 and LlIRF5) may participate in the anti-digenean immune response of L. littorea.

NEONATAL SUSCEPTIBILITY TO INFECTION WITH SCHISTOSOMA MANSONI IN RESISTANT BIOMPHALARIA GLABRATA

The laboratory BS-90 strain of the freshwater pulmonate snail Biomphalaria glabrata, progeny of snails collected from Salvador, Brazil, is resistant to infection with Schistosoma mansoni as juveniles or adults, which rapidly kill primary sporocysts with an attack by the internal defense system (IDS). However, neonatal snails are susceptible to infection. Although neonatal susceptibility of Salvador B. glabrata was reported in 1953 and confirmed subsequently, this phenomenon has been largely ignored. In this study, susceptibility was examined in discrete sizes (shell diameters) of BS-90 snails. We found that 1 mm snails are highly susceptible and develop patent infections. Unexpectedly, most infected 1 mm snails contain primary sporocysts in the digestive gland. Snails measuring 2 and 3 mm show reduced prevalence of infection, and 4 mm and larger snails are refractory. In snails larger than 1 mm, sporocysts fail to develop normally, as shown by reduced numbers of germinal cells at 48 hr post exposure. Moreover, in larger snails an increasingly stronger response of the IDS is mounted in the form of increased numbers of sporocysts undergoing encapsulation and destruction by hemocytes, increased layers of encapsulating hemocytes, as well as increased mitotic activity of the hematopoietic amebocyte-producing organ. These results indicate a relatively narrow size range over which resistance develops and suggest that the IDS of 1 mm snails is developmentally immature. The occurrence of infections in neonatal snails may help to explain transmission of schistosomiasis in regions of low snail susceptibility and may complicate future efforts in biological control.