Identification of differentially expressed genes in Oncomelania hupensis chronically infected with Schistosoma japonicum

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.

Identification of Toll-like receptor family members in Oncomelania hupensis and their role in defense against Schistosoma japonicum

The amphibious snail, Oncomelania hupensis, primarily distributed in the Far East, is the only intermediate host of Schistosoma japonicum, which causes the most virulent form of schistosomiasis. Obligatory parasitism of snails is the main vehicle for human and livestock infection and depends primarily on parasite infectivity, snail defense capacity and specificity, and parasite-snail compatibility. Therefore, the schistosome-snail interaction is biomedically significant, particularly the molecular mechanisms involved in the innate immune response against S. japonicum. Several immune effectors and signaling pathways have been successfully identified in mollusks, especially in Biomphalaria glabrata, the intermediate snail host of S. mansoni; however, limited information is available for O. hupensis. Here, we identified 16 Toll-like receptors (TLRs) in O. hupensis. These O. hupensis TLRs (OhTLRs) are highly expressed in haemocytes, the primary immune cell of mollusks. Most of the OhTLRs were more highly expressed in female gonads than in other tissues, which may suggest maternal immune transfer in O. hupensis. After S. japonicum challenge, the expression levels of all of the OhTLRs were significantly up-regulated at 6 h post-challenge; many of the OhTLR expression levels were inhibited at later time points in haemocytes, while they were inhibited and fluctuated to varying degrees in other tissues. Additionally, we further determined the tissue-specific expression and dynamic response against S. japonicum of one of the TLR signaling adaptors, myeloid differentiation factor 88 (MyD88), from O. hupensis. Three OhMyD88 genes were highly expressed in haemocytes, and were up-regulated in haemocytes and inhibited in the head-foot muscle at the early time-point after S. japonicum challenge; however, these had slower changes and longer durations compared to OhTLRs. These results provide evidence suggesting that immune effectors are involved in innate immune responses of O. hupensis against S. japonicum and may play a role in the activation of different haemocytes, and not limited for the early response to S. japonicum invasion. Further investigation into the varied expression of OhTLRs in other tissues after S. japonicum challenge will improve our understanding of TLR function in innate immunity of O. hupensis.

Transcriptomic responses of Biomphalaria pfeifferi to Schistosoma mansoni: Investigation of a neglected African snail that supports more S. mansoni transmission than any other snail species

Background

Biomphalaria pfeifferi is highly compatible with the widespread human-infecting blood fluke Schistosoma mansoni and transmits more cases of this parasite to people than any other snail species. For these reasons, B. pfeifferi is the world’s most important vector snail for S. mansoni, yet we know relatively little at the molecular level regarding the interactions between B. pfeifferi and S. mansoni from early-stage sporocyst transformation to the development of cercariae.

Methodology/Principal findings

We sought to capture a portrait of the response of B. pfeifferi to S. mansoni as it occurs in nature by undertaking Illumina dual RNA-Seq on uninfected control B. pfeifferi and three intramolluscan developmental stages (1- and 3-days post infection and patent, cercariae-producing infections) using field-derived west Kenyan specimens. A high-quality, well-annotated de novo B. pfeifferi transcriptome was assembled from over a half billion non-S. mansoni paired-end reads. Reads associated with potential symbionts were noted. Some infected snails yielded fewer normalized S. mansoni reads and showed different patterns of transcriptional response than others, an indication that the ability of field-derived snails to support and respond to infection is variable. Alterations in transcripts associated with reproduction were noted, including for the oviposition-related hormone ovipostatin and enzymes involved in metabolism of bioactive amines like dopamine or serotonin. Shedding snails exhibited responses consistent with the need for tissue repair. Both generalized stress and immune factors immune factors (VIgLs, PGRPs, BGBPs, complement C1q-like, chitinases) exhibited complex transcriptional responses in this compatible host-parasite system.

Significance

This study provides for the first time a large sequence data set to help in interpreting the important vector role of the neglected snail B. pfeifferi in transmission of S. mansoni, including with an emphasis on more natural, field-derived specimens. We have identified B. pfeifferi targets particularly responsive during infection that enable further dissection of the functional role of these candidate molecules.

Biomphalaria pfeifferi is the world’s most important snail vector for the widespread human-infecting blood fluke Schistosoma mansoni. Despite this, we know relatively little about the biology of this highly compatible African snail host of S. mansoni, especially for specimens from the field. Using an Illumina-based dual-seq approach, we captured a portrait of the transcriptional responses of Kenyan snails that were either uninfected with S. mansoni, or that harbored 1-day, 3-day, or cercariae-producing infections. Responses to infection were influenced both by the extent of schistosome gene expression and infection duration. We note and discuss several alterations in transcriptional activity in immune, stress and reproduction related genes in infected snails and the B. pfeifferi symbionts detected. Several host genes were highly up-regulated following infection and these might comprise excellent candidates for disruption to diminish compatibility. This study provides for the first time a large sequence dataset to help in interpreting the important vector role of B. pfeifferi in transmission of S. mansoni, including with an emphasis on more natural, field-derived specimens.

De Novo transcriptome analysis of Oncomelania hupensis after molluscicide treatment by next-generation sequencing: implications for biology and future snail interventions

The freshwater snail Oncomelania hupensis is the only intermediate host of Schistosoma japonicum, which causes schistosomiasis. This disease is endemic in the Far East, especially in mainland China. Because niclosamide is the only molluscicide recommended by the World Health Organization, 50% wettable powder of niclosamide ethanolamine salt (WPN), the only chemical molluscicide available in China, has been widely used as the main snail control method for over two decades. Recently, a novel molluscicide derived from niclosamide, the salt of quinoid-2',5-dichloro-4'-nitro-salicylanilide (Liu Dai Shui Yang An, LDS), has been developed and proven to have the same molluscicidal effect as WPN, with lower cost and significantly lower toxicity to fish than WPN. The mechanism by which these molluscicides cause snail death is not known. Here, we report the next-generation transcriptome sequencing of O. hupensis; 145,008,667 clean reads were generated and assembled into 254,286 unigenes. Using GO and KEGG databases, 14,860 unigenes were assigned GO annotations and 4,686 unigenes were mapped to 250 KEGG pathways. Many sequences involved in key processes associated with biological regulation and innate immunity have been identified. After the snails were exposed to LDS and WPN, 254 unigenes showed significant differential expression. These genes were shown to be involved in cell structure defects and the inhibition of neurohumoral transmission and energy metabolism, which may cause snail death. Gene expression patterns differed after exposure to LDS and WPN, and these differences must be elucidated by the identification and annotation of these unknown unigenes. We believe that this first large-scale transcriptome dataset for O. hupensis will provide an opportunity for the in-depth analysis of this biomedically important freshwater snail at the molecular level and accelerate studies of the O. hupensis genome. The data elucidating the molluscicidal mechanism will be of great benefit in future snail control efforts.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction

Background

The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis.

Results

A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed.

Conclusions

This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.