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

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.

Accelerating snail vector genomics

The three most important genera of snails for the transmission of schistosomes are Bulinus, Biomphalaria and Oncomelania. Each of these genera, found in two distantly related families, includes species that act as the intermediate host for one of the three most widespread schistosome species infecting humans, Schistosoma haematobium, S. mansoni and S. japonicum, respectively. An important step in the fight against schistosomiasis in Asia has been taken with the publication of the article "Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum", which means that genomes for all three major genera, including species across three continents, are now available in the public domain. This includes the first genomes of African snail vectors, namely Biomphalaria sudanica, Bi. pfeifferi and Bulinus truncatus, as well as high-quality chromosome level assemblies for South American Bi. glabrata. Most importantly, the wealth of new genomic and transcriptomic data is helping to establish the specific molecular mechanisms that underly compatibility between snails and their schistosomes, which although diverse and complex, may help to identify potential targets dictating host parasite interactions that can be utilised in future transmission control strategies. This new work on Oncomelania hupensis and indeed studies on other snail vectors, which provide deep insights into the genome, will stimulate research that may well lead to new and much needed control interventions.

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.

A genome sequence for Biomphalaria pfeifferi, the major vector snail for the human-infecting parasite Schistosoma mansoni

BACKGROUND

Biomphalaria pfeifferi is the world's most widely distributed and commonly implicated vector snail species for the causative agent of human intestinal schistosomiasis, Schistosoma mansoni. In efforts to control S. mansoni transmission, chemotherapy alone has proven insufficient. New approaches to snail control offer a way forward, and possible genetic manipulations of snail vectors will require new tools. Towards this end, we here offer a diverse set of genomic resources for the important African schistosome vector, B. pfeifferi.

METHODOLOGY/PRINCIPAL FINDINGS

Based largely on PacBio High-Fidelity long reads, we report a genome assembly size of 772 Mb for B. pfeifferi (Kenya), smaller in size than known genomes of other planorbid schistosome vectors. In a total of 505 scaffolds (N50 = 3.2Mb), 430 were assigned to 18 large linkage groups inferred to represent the 18 known chromosomes, based on whole genome comparisons with Biomphalaria glabrata. The annotated B. pfeifferi genome reveals a divergence time of 3.01 million years with B. glabrata, a South American species believed to be similar to the progenitors of B. pfeifferi which undertook a trans-Atlantic colonization < five million years ago.

CONCLUSIONS/SIGNIFICANCE

The genome for this preferentially self-crossing species is less heterozygous than related species known to be preferential out-crossers; its smaller genome relative to congeners may similarly reflect its preference for selfing. Expansions of gene families with immune relevance are noted, including the FReD gene family which is far more similar in its composition to B. glabrata than to Bulinus truncatus, a vector for Schistosoma haematobium. Provision of this annotated genome will help better understand the dependencies of trematodes on snails, enable broader comparative insights regarding factors contributing to susceptibility/ resistance of snails to schistosome infections, and provide an invaluable resource with respect to identifying and manipulating snail genes as potential targets for more specific snail control programs.

The transcriptomic analysis of Planorbarius corneus hemocytes (Gastropoda) naturally infected with Bilharziella polonica (Schistosomatidae)

The study of molluscan innate immunity is essential for understanding the evolution of the immune system. An advance in the knowledge of their immune system can be achieved by increasing the number of model species. Our study focuses on the immunity of Planorbarius corneus, a pulmonate snail widely distributed in Eurasia. These snails are intermediate hosts of many trematodes, including Bilharziella polonica (Schistosomatidae). In this paper we obtained and analyzed transcriptomes of hemocytes of uninfected snails Planorbarius corneus and snails naturally infected with Bilharziella polonica. The transcriptomes were found to contain transcripts encoding all major groups of immune factors previously described for other gastropods. Pathogen-recognition molecules were the most diverse group of immune factors. Comparison of the transcriptomes of the infected and the uninfected molluscs showed that the expression of some genes changed during infection. Our results extend the knowledge of immune responses of pulmonate snails to trematode invasion and promote P. corneus as a new model for the study of molluscan defence reactions.

Compatibility between snails and schistosomes: insights from new genetic resources, comparative genomics, and genetic mapping

The freshwater snail Biomphalaria glabrata is an important intermediate host of the parasite Schistosoma mansoni that causes human intestinal schistosomiasis. To better understand vector snail biology and help advance innovative snail control strategies, we have developed a new snail model consisting of two homozygous B. glabrata lines (iM line and iBS90) with sharply contrasting schistosome-resistance phenotypes. We produced and compared high-quality genome sequences for iM line and iBS90 which were assembled from 255 (N50 = 22.7 Mb) and 346 (N50 = 19.4 Mb) scaffolds, respectively. Using F2 offspring bred from the two lines and the newly generated iM line genome, we constructed 18 linkage groups (representing the 18 haploid chromosomes) covering 96% of the genome and identified three new QTLs (quantitative trait loci), two involved in snail resistance/susceptibility and one relating to body pigmentation. This study provides excellent genomic resources for unveiling complex vector snail biology, reveals genomic difference between resistant and susceptible lines, and offers novel insights into genetic mechanism of the compatibility between snail and schistosome.

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.

Structural insights into the redox regulation of Oncomelania hupensis TRP14 and its potential role in the snail host response to parasite invasion

The freshwater amphibious snail Oncomelania hupensis is the unique intermediate host of Schistosoma japonicum, but little attention has been paid to the interaction between the two. In snails, the production of reactive oxygen species (ROS) by hemocytes has been shown to be vital for snail immune defense against schistosome infection. However, excessive ROS accumulation could lead to oxidative damage, requiring the antioxidant system for maintaining the cellular redox homeostasis. Previously we identified a thioredoxin-related protein of 14 kDa from O. hupensis (OhTRP14), and showed that it was involved in the scavenging of ROS in circulating hemocytes. Here, we confirmed that OhTRP14 plays a potential role in the snail host response to parasite challenge and determined the crystal structures of OhTRP14 in two different states (oxidized and transition state). The overall structure revealed a typical Trx fold and is similar to that of human TRP14 (hTRP14), but there were significant structural differences between the two states. Noticeably, there was a different pair of thiol groups from Cys30 and Cys44 in the transition state of OhTRP14, were with the similar separation of 2.9 Å as that (2.6 Å) between Cys41 and Cys44, but in a different orientation, suggesting that the Cys30 is likely to function as an important molecular switch involved in the oxidoreductase activity of OhTRP14. Comparative studies between OhTRP14 and hTRP14 by analyzing the surface characteristics, charge distribution and oxidoreductase activity toward insulin demonstrated they might have similar substrates. The results are expected to provide structural insights into the redox regulation of OhTRP14 and contribute to better understanding of TRP14 family. DATA DEPOSITION: The atomic coordinates of the structure and the structure factors were deposited in Protein Data Bank with PDB ID codes 7XQ3 and 7XPW.