Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

A Novel Toll-Like Receptor (TLR) Influences Compatibility between the Gastropod Biomphalaria glabrata, and the Digenean Trematode Schistosoma mansoni

Abstract

Schistosomiasis, a devastating disease caused by parasitic flatworms of the genus Schistosoma, affects over 260 million people worldwide especially in tropical and sub-tropical regions. Schistosomes must undergo their larval development within specific species of snail intermediate hosts, a trait that is shared among almost all digenean trematodes. This unique and long-standing host-parasite relationship presents an opportunity to study both the importance of conserved immunological features in novel immunological roles, as well as new immunological adaptations that have arisen to combat a very specific type of immunological challenge. While it is well supported that the snail immune response is important for protecting against schistosome infection, very few specific snail immune factors have been identified and even fewer have been functionally characterized. Here, we provide the first functional report of a snail Toll-like receptor, which we demonstrate as playing an important role in the cellular immune response of the snail Biomphalaria glabrata following challenge with Schistosoma mansoni. This TLR (BgTLR) was identified as part of a peptide screen of snail immune cell surface proteins that differed in abundance between B. glabrata snails that differ in their compatibility phenotype to challenge by S. mansoni. The S. mansoni-resistant strain of B. glabrata (BS-90) displayed higher levels of BgTLR compared to the susceptible (M-line) strain. Transcript expression of BgTLR was found to be very responsive in BS-90 snails when challenged with S. mansoni, increasing 27 fold relative to β-actin (non-immune control gene); whereas expression in susceptible M-line snails was not significantly increased. Knockdown of BgTLR in BS-90 snails via targeted siRNA oligonucleotides was confirmed using a specific anti-BgTLR antibody and resulted in a significant alteration of the resistant phenotype, yielding patent infections in 43% of the normally resistant snails, which shed S. mansoni cercariae 1-week before the susceptible controls. Our results represent the first functional characterization of a gastropod TLR, and demonstrate that BgTLR is an important snail immune receptor that is capable of influencing infection outcome following S. mansoni challenge.

Author Summary

The freshwater snail Biomphalaria glabrata is the subject of intensive research, primarily due to its biomedical importance as an intermediate host for the parasitic flatworm Schistosoma mansoni, which is a causative agent of the disease schistosomiasis–a chronic, debilitating condition that affects over 260 million people worldwide. Studies of this snail have led to the identification of many factors, some with high homology to known immune molecules in other organisms and many more that are unique with no known homology. However, research into the functional and mechanistic roles of these factors has only recently progressed. In this study, we have demonstrated the functional relevance of a Toll-like receptor (TLR) that we identified in B. glabrata. Transcriptional expression of this TLR was rapidly induced in a snail strain that is resistant to S. mansoni, in a pattern that is consistent with their phenotype. Furthermore, knockdown of the TLR resulted in 43% of resistant snails becoming infected with S. mansoni to patency. This advances our understanding of the mechanistic basis of snail-schistosome compatibility significantly, and may one day facilitate the development of tools for improving the control of schistosomiasis. It also contributes novel functionality to TLRs, one of the most evolutionarily conserved pattern-recognition receptors and cognate signalling pathways in immunity.

Citation: Pila EA, Tarrabain M, Kabore AL, Hanington PC (2016) A Novel Toll-Like Receptor (TLR) Influences Compatibility between the Gastropod Biomphalaria glabrata, and the Digenean Trematode Schistosoma mansoni. PLoS Pathog 12(3): e1005513. https://doi.org/10.1371/journal.ppat.1005513

Editor: Matty Knight, George Washington University School of Medicine and Health Sciences, UNITED STATES

Received: October 8, 2015; Accepted: February 29, 2016; Published: March 25, 2016

Copyright: © 2016 Pila et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was supported by the Natural Sciences and Engineering Council of Canada, Discovery Grant #418540. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Schistosomiasis is a devastating disease caused by parasitic flatworms of the genus Schistosoma. It affects over 260 million people worldwide especially in tropical and sub-tropical regions [1,2]. The freshwater snail Biomphalaria glabrata acts an obligate intermediate host for Schistosoma mansoni, which is the causative agent of intestinal schistosomiasis in Africa and South America. This association to human schistosomiasis has made B. glabrata one of the most extensively studied gastropods in terms of immunobiology and host-parasite interactions. This snail continues to play an important role as a model for studying the intra-molluscan aspects of the parasite lifecycle, which has gained popularity as a possible target for disease control purposes [3]. Strains of B. glabrata have been bred that display differing compatibility phenotypes to S. mansoni infection. Strains displaying resistance such as BS-90, 13-16-R1, 10-R2 [4,5] or susceptibility such as the M-line and NMRI [6,7], serve as a means to evaluate and better understand the driving mechanisms underpinning snail resistance and susceptibility. They provide valuable tools for elucidating the specifics of these naturally occurring processes [8].

An improved understanding of the molecular basis for susceptibility/resistance is considered important due to the potential for the development of novel control strategies for schistosomiasis. Moreover, such knowledge would contribute significantly to the field of evolutionary and invertebrate immunology, particularly of the Lophotrochozoa, a superphylum of Metazoa to which molluscs belong. Compared to the other two superphyla [Ecdysozoa (represented by fruit flies and nematodes) and Deuterostoma (represented by vertebrates)], with respect to our understanding of immunological capability and function, this group remains under-represented in the literature [9].

Studies utilizing the B. glabrata/S. mansoni and other Biomphalaria/trematode models have made significant progress in the past two decades, with the focus shifting from classic comparative immunology to one with molecular insight [8,10]. These advances have been made possible by the availability and integration of techniques such as proteomics [11–14], functional genomics [15–19] and population genetics [20–24]. Gene discovery, comparative proteomics and differential gene expression studies with these snails, their haemocytes, or various tissues, have led to the identification of hundreds of immune factors with putative influence on snail-schistosome compatibility [12,14,16, 25–30].

Few studies have dealt with the functional and mechanistic roles of these factors using defined bioassays. Of those that have been characterized, the fibrinogen-related proteins (FREPs) are perhaps the best understood. FREPs comprise a diverse family of secreted lectins that have been shown to be important in the anti-trematode response [31,32]. They are known to act as opsonins, recognizing schistosome surface proteins [31] and excretory/secretory products [33], forming complexes with another diverse family of proteins on larval schistosomes–the Schistosoma mansoni polymorphic mucins [34]. Biomphalaria glabrata migration inhibitory factor (BgMIF) is another factor for which functional involvement in the anti-schistosome response has been demonstrated. BgMIF was shown to stimulate cellular proliferation by activating extracellular signal-regulated kinase and to suppress nitric oxide-induced apoptosis in Biomphalaria glabrata embryonic cells. Knockdown of BgMIF reduced encapsulation of sporocysts by these cells and increased parasite burden in the snails [35]. Other molecules that have been studied functionally include the copper/zinc superoxide dismutase and its products [36–38], and the beta pore-forming toxin biomphalysin [39]. To date, a mechanism or receptor demonstrating how these factors might engage with and elicit the cellular encapsulation response that appears to be necessary for clearance of an invading trematode has not been identified.

Transition into studies of this nature is required in order to piece together the big picture of how the snail immune response works, and to fully understand the molecular basis of snail-trematode interactions. The evidence thus far indicates that a number of conserved pathways involved in the vertebrate immune response or their interacting components are also present in B. glabrata. These include pathways of the reactive oxygen and nitrogen intermediates mediated by superoxide dismutase, and macrophage migration inhibitory factor, mentioned above. The mitogen-activated protein kinase pathway is also conserved in the snail and shown to be relevant in haemocyte motility, spreading, phagocytosis, and encapsulation [40–43].

One of the major evolutionarily conserved pathways in innate immunity is the Toll/Toll-like receptor (TLR) pathway. TLRs are trans-membrane proteins composed of an extracellular leucine-rich repeat domain, responsible for pathogen recognition and a conserved cytoplasmic Toll/IL-1 (TIR) domain, which is responsible for signal transduction and activation of effector functions. In 2007, Zhang and colleagues [44] published homologues of Gram-negative bacteria binding protein and peptidoglycan recognition protein in B. glabrata. Both are extracellular components that can activate TLR and the related immune-deficiency (IMD) pathways respectively. Later, they also identified the nuclear factor kappa B (NF-kB), a major downstream transcription factor in the TLR pathway. NF-kB expression patterns in resistant snails infected with schistosome parasites was found to be consistent with the early immune response patterns critical for parasite killing [9]. However, it was not known whether, and how, TLR was involved in regulating parasite infection in this snail. While comparing haemocyte surface proteins between S. mansoni-susceptible and resistant B. glabrata snails, we discovered peptides of a TLR-like protein, some of which displayed a high amino acid sequence identity with the TIR domain of known vertebrate TLRs. Peptides associated with this protein were at a significantly higher abundance in the resistant snail strain, prompting us to further investigate this putative TLR.

In this study, we have conducted the first functional characterization of a TLR in the snail B. glabrata. TLRs play a key role in the innate immune response by directly recognizing a variety of pathogens (typically bacteria, viruses and fungi) or factors evidencing their presence and transducing signals to the immune cells [45]. The leucine-rich repeat (LRR) region of TLRs recognizes conserved motifs on pathogens known as the pathogen-associated molecular patterns (or PAMPs) such as bacterial lipopolysaccharide, flagellin protein and genomic DNA containing unmethylated CpG motifs [46–48], viral DNA [49], fungal zymosan [50] and the plasmodium pigment hemozoin [51]. The Toll/TLR pathway can also be activated by host endogenous ligands or molecules that signal tissue damage–the so-called damage-associated molecular patterns (DAMPs). These include molecules such as uric acid, fibrinogen, mitochondrial DNA, heat-shock proteins 60 and 70 and fibronectin [52]. Drosophila Toll, the founding member of the TLR family is not a direct recognition receptor. Instead, pathogen components trigger the activation of protease cascades that then lead to the cleavage of the Toll ligand–Spaetzle [53]. The cytoplasmic region of TLRs which shares homology with the interleukin-1 receptor (known as a Toll/interleukin-1 receptor (TIR) domain) [45] is responsible for transducing signals from pathogen recognition to the immune cells leading to the activation of the effector functions. This pathway is evolutionarily conserved from nematodes to mammals [54].

The B. glabrata TLR we have characterized in this study (GenBank accession number: JX014259.1, herein after referred to as BgTLR) possesses complete LRR and TIR domains, and we demonstrate here the involvement of BgTLR in the immune response of B. glabrata against S. mansoni. This implies that BgTLR has a novel functionality with respect to helminth parasites, making it unique among known TLRs, one of the most evolutionarily conserved pattern-recognition receptors and cognate signalling pathways in immunity.

Results

BgTLR was identified as part of a haemocyte protein expression screen comparing resistant and susceptible snails

Membrane-associated proteins displaying differential expression in haemocytes of BS-90 B. glabrata when compared to the M-line strain were identified using iTRAQ and LC/MS/MS. Including BgTLR, 16 proteins were identified with confidence, as being present at higher abundance in BS-90 snail haemocytes (S1 Table). Four proteins could not be matched with any factors that have a known function, however, of the 12 identified proteins superoxide dismutase Cu/Zn [36, 55, 56], dermatopontin 2, matrilin, an elastase-like protein [30, 31], and a protein with significant similarity to protein tyrosine phosphatase domain-containing proteins [57], have all been associated with resistant snail phenotypes in other studies.

BgTLR displayed increased transcript expression in BS-90 snails following challenge with S. mansoni

To determine whether BgTLR transcript expression was responsive to S. mansoni challenge, we used quantitative polymerase chain reaction (qPCR) to measure its expression patterns in schistosome-challenged snails, which were compared to control non-challenged snails at each experimental time point, and the endogenous control β-actin. BgTLR expression was rapidly induced in the S. mansoni-resistant (BS-90) snails responding to S. mansoni challenge, with transcript abundance significantly increasing 14 fold as early as 12 hours post challenge compared to the 0-hour time point. The highest value recorded was a 27 fold increase at 1 day post challenge, which was significantly different from all other time points (P < 0.05). This was followed by reduced BgTLR transcript abundance that returned to pre-exposure levels by day 8 post challenge (dpc) (Fig 1A). In contrast, S. mansoni challenge did not induce significant changes in BgTLR transcript expression in the susceptible (M-line) snails (Fig 1A). Non-exposed BS-90 and M-line snails expressed very similar levels of BgTLR transcript as determined by comparison of absolute expression (Ct values) which ranged approximately 7.2 to 7.5.

thumbnail
Download:
Fig 1. BgTLR displayed increased transcript expression in BS-90 snails following challenge with S. mansoni.

(A) BgTLR transcript expression with and without S. mansoni challenge. Snails (BS-90 and M-line strains) were individually exposed to ~5 miracidia or left unexposed (control). Five snails were collected at indicated time points over the incubation period of the parasite. RNA was extracted from whole snails, converted to cDNA and BgTLR expression was measured by quantitative PCR. Expression was quantified in fold changes normalized to time 0-hour controls. Bars represent standard error (n = 5). Asterisk (*) indicates significant difference (P < 0.05) between experimental and control samples, while hash (#) indicates significant difference between BS-90 and M-line snails at the respective time points. (B) Knockdown of BgTLR or GFP (control) in BS-90 snails responding to S. mansoni parasite challenge. (C) Knockdown of BgTLR in unexposed M-line and BS-90 snails. Five snails were collected for each time point in B and C for RNA extraction and cDNA synthesis. Bars represent standard error (n = 5). Asterisks indicate significant difference (P < 0.05) from 0-hour time points.

https://doi.org/10.1371/journal.ppat.1005513.g001

Predicted structure of BgTLR

The full-length coding sequence of BgTLR is composed of 3546 nucleotides. The open reading frame predicts a protein of 1181 amino acid residues. Its ectodomain contains 23 LRR motifs spanning residues 3–725 as predicted by LRRfinder [58] and ScanProsite tools [59]. The TIR domain is encoded by residues 899–1037, preceded by a transmembrane region on residues 849–870 (S1 Fig). There is no identifiable signal peptide in BgTLR based on current predictive software. BlastP analyses [60] indicate that the TIR region of BgTLR shares the highest identity to other molluscan TLR TIR regions (65–88%), followed by those of arthropods (51–58%). Among the mammals, TLRs having closest TIR region identities to BgTLR are TLR 4 (40%), TLR 3 (38%) and TLR 13 (37%). The LRR region of BgTLR is much more variable, with identity range of 27–43% only to other molluscan and arthropod TLRs.

Amplification specificity was confirmed using qPCR melt curve and customized BLAST sequence analyses (S2 Fig and S3 Fig) while siRNA knockdown specificity was assessed by measuring the abundance of the three sequences in the B. glabrata genome that shared the highest nucleotide identity with BgTLR using BgTLR knockdown cDNA as a template (S5 Fig). The peptide fragment used for the synthesis of the antibody used in this study for Western blot and immunocytochemistry analyses was also confirmed to be specific by testing the binding of peptide pre-incubated antibody to BgTLR (S7 Fig).

siRNA-mediated knockdown of BgTLR

Injection of an siRNA oligonucleotide cocktail targeting the extracellular LRR region of BgTLR (Fig 2A) into M-line and BS-90 snails induced a measurable knockdown of BgTLR expression in both snail strains. Maximum knockdown was observed at 3 days post injection with transcript levels of 0.09 and 0.04 fold compared to the time 0-hour control for M-line and BS-90 strains respectively. Both were significantly different from time 0-hour control (P < 0.05) (Fig 1C). Knockdown effect of BgTLR siRNA was sustained in BS-90 snails following S. mansoni challenge with BgTLR transcript levels reducing down to 0.4 fold at 3-days post injection. In comparison, snails injected with siRNA oligonucleotides targeted to GFP as a control displayed a 3 fold BgTLR expression increase at 3 days post injection (Fig 1B). BgTLR knockdown kinetics at the protein level in BS-90 snails lagged behind transcript knockdown with an observable reduction in BgTLR protein abundance at 4 days post injection as determined by Western blots using both haemocytes isolated from BS-90 snails and using B. glabrata embryonic (Bge) cell line (Fig 2B).

thumbnail
Download:
Fig 2. siRNA-mediated knockdown of BgTLR transcripts.

(A) Graphic view of BgTLR (not drawn to scale) showing the leucine-repeat motifs and Toll/IL-1 domain. Red bars indicate the approximate positions targeted by siRNA oligonucleotides used to inject the snails, while the blue bar indicates the targeted position of the antibody used. (B) Confirmation of BgTLR (~135 kDa) knockdown in BS-90 haemocytes and Bge cells respectively. Protein was detected via Western blot with primary antibody developed against the extracellular region of BgTLR. BgActin (~42 kDa) served as the loading control.

https://doi.org/10.1371/journal.ppat.1005513.g002

BgTLR immunolocalization in haemocytes and Bge cells

Immunocytochemical studies using an antibody developed against the extracellular region of BgTLR (Fig 2A) indicate that it is expressed in the majority of haemocytes found in the circulation of B. glabrata, although some haemocytes appeared not to express BgTLR (S6 Fig). In the Bge cell line, all cells appear to be BgTLR-positive with differences in staining intensity among cells (Fig 3). BgTLR positive signals were detected primarily on the extracellular surface of the plasma membrane, and faint immunolocalization was observed in the cytoplasm of both haemocytes and Bge cells (Fig 3 and S1 Video). Intracellular localization may be due to BgTLR trafficking between the cellular membrane and intracellular compartments similar to TLR 4 expression in human monocytes [61] or it may indicate that BgTLR has a role in the intracellular environment. We currently have studies underway focusing on identifying the BgTLR ligand(s) which might clarify its subcellular localization and possible intracellular role.

thumbnail
Download:
Fig 3. BgTLR immunolocalization in haemocytes and Bge cells.

Haemocytes (A-C) and Bge cells (D-E) were labelled with DAPI and anti-BgTLR primary antibody. Control haemocyte (C) and Bge (E) samples were stained with DAPI but the primary antibody step omitted. Scale bars represent 20 μM.

https://doi.org/10.1371/journal.ppat.1005513.g003

Knockdown of BgTLR in BS-90 snails decreases haemocyte phagocytic response

Knockdown of BgTLR significantly (P<0.05) decreased the mean number of S. mansoni sporocyst excretory/secretory/lysate-coated beads phagocytosed by individual haemocytes of BS-90 snails. Snails injected with siRNA oligonucleotides targeting BgTLR and 96 hours later assessed for phagocytic activity averaged 13.1 beads per haemocyte, whereas, GFP knockdown control snails averaged 17 beads per haemocyte (Fig 4A). A breakdown of the number of beads observed in each haemocyte demonstrated that haemocytes from BgTLR KD snails were much less likely to take in more than 25 beads per haemocyte compared to the GFP controls (Fig 4B).

thumbnail
Download:
Fig 4. Knockdown of BgTLR in BS-90 snails decreases haemocyte phagocytic response.

(A) Mean number of phagocytosed beads per haemocyte. Five BS-90 snails each were injected with siRNA targeting BgTLR or GFP (control) and 96 hours later, haemolymph was extracted from the snails and immediately mixed with ~1 x 106 1μM FITC-labelled streptavidin-coated beads pre-incubated with biotinylated S. mansoni excretory/secretory products and sporocyst. After 3 hours, haemocytes from each snail were counted from a random field of view on the slide, and 30 haemocytes for each snail were assessed for the number of beads within each cell from which mean number of beads per haemocyte was calculated. Asterisk (*) indicates significant reduction (P < 0.05) in the mean number of beads per haemocyte in BgTLR knockdown snails. (B) Frequency of number of beads observed in 30 haemocytes from one field of view (FOV). Asterisk (*) indicates significant difference (P < 0.05) in the average number of beads phagocytosed between BgTLR and GFP knockdown.

https://doi.org/10.1371/journal.ppat.1005513.g004

Knockdown of BgTLR in resistant BS-90 snails yields cercariae-producing S. mansoni infections

Knockdown of BgTLR abrogated the resistant phenotype of BS-90 snails to S. mansoni infection. Snails injected with siRNA oligonucleotides and challenged with S. mansoni 2-days later developed patent infections that led to cercariae shedding in ~ 43% of BS-90 snails over two trials, whereas none of the control GFP-siRNA injected BS-90 snails shed cercariae (Fig 5). Knockdown of BgTLR resulted in cercariae shedding in the BS-90 snails one week prior to the susceptible M-line snails, which served as a benchmark for infection success. At week 4 post exposure, 2% of the BS-90 BgTLR knockdown snails shed cercariae. This increased to 43% at week 5 and then decreased to 22% at week 7 when the experiment was terminated. The proportions of BS-90 BgTLR knockdown snails shedding cercariae at weeks 5, 6 and 7 were statistically significant (P<0.05) compared to BS-90 GFP knockdown snails. In comparison, M-line snails started shedding cercariae only at week 5 (46%), increasing to 56% at week 7 (Fig 5). There was no statistically significant difference at week 5 in the proportions of snails shedding cercariae between BS-90 BgTLR knockdown and the M-line susceptible control. Injection of siRNA oligonucleotides resulted in an overall mortality of 27%, 17% and 19% for BS-90 BgTLR knockdown, BS-90 GFP knockdown and the M-line susceptible controls respectively. It is worth noting that these mortalities might have impacted the percentages of shedding snails which were calculated out of the total starting number. Since snails were not tracked individually over the course of the experiment, these percentages likely underestimate shedding snails that were lost in the mortalities.

thumbnail
Download:
Fig 5. Knockdown of BgTLR in resistant BS-90 snails yields cercariae-producing S. mansoni infections.

BS-90 snails were injected with BgTLR or GFP (knockdown control) siRNA oligonucleotide mix and after 48 hours exposed individually to ~5 miracidia. Non-injected M-line snails served as the susceptible controls for parasite viability. The number of cercariae-producing snails was expressed as a percentage of the starting number (n) of snails. Figure represents pooled data from two independent experiments. Asterisk (*) indicates significant difference from BS-90 GFP knockdown control for each week while hash (#) indicates significant difference between BS-90 BgTLR knockdown and M-line susceptible control for each week.

https://doi.org/10.1371/journal.ppat.1005513.g005

Discussion

Biomphalaria glabrata has received tremendous research attention due to its biomedical importance as an intermediate host of the causative agent of human schistosomiasis [1,2]. In the wake of the World Health Organization’s call for global elimination of schistosomiasis by the year 2025 [62], there is intensified interest to understand the molecular basis of the interactions that take place between schistosomes, and their snail hosts, with a view to identifying targets or mechanisms that could be exploited in disease prevention and control efforts. As a primary determinant of compatibility between this host and parasite, the snail internal defense response constitutes one area worthy of in-depth investigation in this context.

In this study, we have characterized a TLR in B. glabrata that we demonstrate is immunologically relevant during the anti-schistosome immune response of the snail. Our study provides significant evidence that BgTLR is involved in regulating parasite infection in this snail. BgTLR expression was increased as early as 3 hours following S. mansoni challenge in the resistant (BS-90) snails, lasting up to 3 days before returning to pre-exposure levels (Fig 1A). This pattern of expression is consistent with the involvement of BgTLR in the immune response because the duration coincides with the critical period for the clearance of parasite infection in incompatible strains [63–65]. Moreover, the means by which BgTLR was initially identified suggested that it was present in higher abundance in BS-90 snails when compared to M-line. This is further supported by siRNA-mediated knockdown of BgTLR, which also had maximal effect within this duration (Fig 1B and 1C). Forty-three percent of resistant snails that received BgTLR-specific siRNA oligonucleotides before S. mansoni challenge went on to shed cercariae after 4 weeks incubation period. On the other hand, resistant snails that did not receive BgTLR-specific siRNA but challenged with S. mansoni cleared the infection within 4 days as determined by a qPCR assay targeting the parasite glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene which failed to detect parasite product beyond day 4 (S4 Fig). Conversely, in susceptible (M-line) snails, GAPDH expression was still detectable after 4 days post-infection and increased throughout the incubation period.

The siRNA oligonucleotides used to knockdown BgTLR transcript are specific as shown with qPCR amplification (S5 Fig) of the three predicted transcripts sharing the highest nucleotide identity with BgTLR available on the Vectorbase database for B. glabrata (www.vectorbase.org). However, their cross-reaction with BGLB008602 leaves the possibility that the phenotypic knockdown effects observed may be due to either or both of these proteins. BGLB008602 is 100% identical to BgTLR but lacks a 1195 bp region of BgTLR and cannot be functionally distinguished from BgTLR using our siRNA knockdown approach. We have reason to believe that both may be splice variants of the same gene based on our finding that the ‘missing’ 1195 bp region is currently annotated as introns 8 and 9 in the B. glabrata genome (www.vectorbase.org). As shown by Western blot (Fig 2 and S7 Fig), no band representative of the size predicted for BGLB008602 (~ 127 kDa) is observed, suggesting that its effect if any is minimal, or in a different context than the anti-S. mansoni immune response. Should these transcripts be true splice variants, the difference in size would most likely impact their interaction with pathogen-associated molecular patterns (PAMPs) or endogenous ligands since the differential region is in the extracellular domain and contains leucine-rich motifs commonly required for ligand recognition.

Currently, we do not know the precise recognition mechanism by which BgTLR is involved in the snail immune response. That is, what PAMPs or perhaps endogenous ligands it recognizes and how this signal is transduced to the haemocyte. TLRs have been characterized in a number of other molluscs such as the disk abalone, Hawaiian bobtail squid, Zhikong scallop, soft-shell clam, mussel and oyster [66–70], however this is the first functional identification of a TLR in a gastropod. To date, all molluscan TLRs that have been functionally assessed are thought to be involved in antimicrobial immune responses. This is typical of most TLRs as they recognize PAMPs ranging from lipids, proteins, lipoproteins and nucleic acids [54]. An interesting feature of BgTLR is that it is induced in response to S. mansoni challenge. This suggests that it might be involved in recognizing PAMPs that are shared between schistosomes and microbes. However, an intriguing possibility is that BgTLR may be exclusive for schistosome (or digenean trematode)-specific epitopes. Should this be the case, it will mark a significant finding in the field of invertebrate immunology, and add a new PAMP to the TLR ligand cannon. To the best of our knowledge, no anti-parasitic recognition/response involving TLRs have been characterized in an invertebrate.

There are two possible sources of exogenous ligand(s) that might form the basis of TLR-based recognition in B. glabrata with respect to schistosomes. The first includes the complex carbohydrates which are major components of larval teguments [71–73]. Secondly, during larval transformation, complex glycoproteins are also released with the other larval transformation products [74] which might be recognized by BgTLR. Larval carbohydrates have been shown to elicit functional immune responses such as the induction of reactive oxygen species production in haemocytes [75–78] and phagocytosis [79]. Both parasite larval carbohydrates and glycoproteins can bind to snail plasma and haemocyte proteins in vitro [80–83].

It also remains possible that the BgTLR ligand(s) might be endogenous. This is also the case for the founding member of the TLR family (Drosophila Toll) and for some mammalian TLRs. In Drosophila, the Toll does not function as a direct recognition receptor of a microbial pattern. Instead, its ligand (spaetzle) is activated first through proteolytic cleavage by pathogen-triggered protease cascades [53]. Examples of molecules that serve as endogenous ligands abound in mammals. These include uric acid, mitochondrial DNA, heat-shock proteins 60 and 70, fibronectin and fibrinogen [52]. In mice, evidence suggests that TLR4 recognizes fibrinogen to induce chemokine secretion in macrophages [84]. Based on this evidence, one candidate group of molecules that may serve as an endogenous ligand for BgTLR is the fibrinogen-related proteins, which are highly diversified family of lectins functionally linked to the anti-parasite response in B. glabrata [31,32]. We are currently undertaking studies with the objective to screen for and identify BgTLR ligand(s).

In terms of functional immune responses driven by BgTLR, our experimental results support the involvement of phagocytosis as a downstream functional response elicited following engagement, and ultimately resulting in the prevention of parasite infection. Knockdown of BgTLR in BS-90 snails significantly reduced the mean number of beads phagocytosed per haemocyte from 17 to 13 (Fig 4A), although the overall percentage of phagocytic haemocytes did not differ from GFP knockdown control snails. This, and the fact that parasites are too large to be phagocytosed by a single haemocyte indicate that other effector responses might also be activated by BgTLR. Phagocytosis together with cytotoxicity and encapsulation are the major haemocyte responses that are induced as a result of receptor recognition of PAMPS in the snail [85]. It is likely that the role of BgTLR in the snail immune response against S. mansoni is through the combined effects of these processes. However, the involvement of encapsulation and cytotoxicity remains to be investigated in the context of BgTLR.

The major pathways that are known to be activated downstream of TLR are present in B. glabrata, and for some, their active involvement in the snail immune response have been demonstrated. These major pathways include the NF-kB, MAPK/ERK and phosphatidylinositol-3 kinase (PI3K). MAPK, ERK and PI3K pathways have been shown to be involved in cellular adhesion, motility and spreading required for phagocytosis and encapsulation [86] as well as in regulating the release of hydrogen peroxide and nitric oxide molecules used in cytotoxicity [77,78]. No functional role has been demonstrated for NF-kB although its expression is consistent with involvement in an immune response [9]. BgTLR may be associated with any of these pathways or their combination in order to activate the effector responses in haemocytes, and this is an area actively under investigation.

BgTLR appears to be selectively expressed in certain haemocytes only (S6 Fig), and differentially on cells of the Bge cell line. The relationship between BgTLR-negative haemocytes with the other haemocyte subsets has not yet been explored. However, expression of pattern recognition receptors on haemocyte subsets is not unique to BgTLR. A member of the fibrinogen-related proteins (FREP3) has been shown to be expressed only by subsets of haemocytes which increase in number during an immune response in the snail [31]. The fact that these immunologically relevant factors are not equally present on all haemocytes of known morphological subset alludes to the likelihood that there is far greater functional diversity within a haemocyte morphotype than is currently understood.

In conclusion, our studies presented here have demonstrated a connection between BgTLR and the resistance phenotype of B. glabrata with respect to S. mansoni challenge. By characterizing the receptor which is a major component of the TLR pathway, our findings complement others that have characterized both upstream and downstream components of the pathway [9,77,78,86] and provides conclusive evidence that this pathway is fully conserved and functionally relevant in the immune response of B. glabrata. It also paves the way for further studies not only to identify the ligands, specific immune responses induced and the haemocyte subsets involved, but also to determine the role of the other B. glabrata TLRs and leucine rich repeat containing molecules.

Materials and Methods

Snails and parasite

Two strains of Biomphalaria glabrata snails were used in this study. The BS-90 strain is resistant to Schistosoma mansoni infection [4,5] while the M-line strain is susceptible [6,7]. Snails were maintained in aerated artificial spring water at 23–25°C, 12-hour day/night cycle and fed red-leaf lettuce as needed. All snail exposures were performed with the NMRI strain [6] of S. mansoni, which was maintained at the University of Alberta, cycling between mice and M-line B. glabrata.

Identification of BgTLR

BgTLR was identified via a peptide screen of B. glabrata haemocyte surface proteins using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS analysis) conducted at the University of Victoria, Canada. Haemocytes were isolated from 10 M-line and 10 BS-90 strain B. glabrata snails by the head-foot retraction method and immediately mixed 1:1 with chilled 2x Sterile Snail Saline (SSS) and placed on ice. The cells were then centrifuged for 10 minutes at 10,000 x g. The supernatants were aspirated and the haemocytes were resuspended in 0.5 mL 1x SSS. The haemocytes were then sonicated 5 times using 30-second pulses and the resulting homogenate was centrifuged for 30 minutes at 4°C and 10,000 x g. The supernatant was collected and ultracentrifuged at 100,000 x g for 1 hour at 4°C to purify cell membrane. The pellets were resuspended in a buffer containing 8 M urea, 30 mM HEPES, 0.5% SDS, 1 mM PMSF, 2 mM EDTA and 10 mM DTT. This solution was sonicated for 3 minutes using 5-second pulses, interspersed with 3-second breaks. Protein concentration was assessed using a Qubit (Life Technologies) following the manufacturers protocols.

Equal amounts of isolated haemocyte membrane proteins were digested and labelled following the iTRAQ manufacturer’s protocols (Applied Biosystems). 100 μg solutions of each haemocyte population were precipitated using chilled acetone and stored at -20°C for 3 hours. The samples were then centrifuged at 20,000 x g for 45 minutes before resuspension in 20 μL of supplied dissolution buffer. Samples were then processed using the supplied reducing/denaturing kits and cysteines were blocked prior to digestion of the membrane proteins using 1μg/μL sequencing grade trypsin (Promega) overnight at 37°C. Digested protein samples were then labelled with iTRAQ reagents 114 (M-line haemocytes) and 117 (BS-90 haemocytes) following the manufacturers protocols. Upon completion of the labeling procedure, the samples were combined and sent for LC-MS/MS analysis. Peptides identified during analysis were initially compared to a custom Mascot database comprised of known proteins and predicted transcript translations from B. glabrata, Lottia gigantica, and Aplysia californica. If no matches were found using this approach, peptide fragments were used in a BlastP search of the GenBank database in order to find the match with the highest amino acid identity. Differential abundance of the membrane proteins between haemocytes of M-line or BS-90 snails was measured as a ratio of the 114:117 labels associated with each analyzed spectra peak.

Peptides relevant to BgTLR, particularly the TIR domain, were used to design primers for rapid amplification of cDNA ends (RACE) PCR (Table 1). RACE successfully amplified the complete BgTLR transcript, which was sequenced and confirmed to be expressed by B. glabrata using RT-PCR.

thumbnail
Download:
Table 1. Primers used in this study.

https://doi.org/10.1371/journal.ppat.1005513.t001

Measurement of BgTLR expression in B. glabrata snails during parasite challenge

Quantitative polymerase chain reaction (qPCR) was used to measure BgTLR expression patterns in schistosome-challenged and control snails. Resistant (BS-90) and susceptible (M-line) snails were exposed individually to ~5 miracidia in 12-well plates. Five snails were collected at selected time points (0, 1, 3, 6, 12 hours and 1, 2, 3, 4, 8, 16, and 35 days post challenge) that represent important milestones in the life cycle of S. mansoni, and encompass the time from miracidia infection to cercarial shedding in the snail. RNA was extracted from whole snails using a commercial spin column-based kit according to the manufacturer’s instructions (Life Technologies). RNA concentration was determined using a UV/visible spectrophotometer and 1 μg was used for first strand cDNA synthesis. The cDNA was diluted five-fold and 5 μL was used as template in qPCR using primers specific for BgTLR and B. glabrata β-actin (Table 1), and a dye-based detection system (Quanta Biosciences). Primers were used at final concentrations of 0.6 μM in a reaction volume of 25 μL. The B. glabrata β-actin gene was used as endogenous control. Snails were confirmed to be S. mansoni positive or negative using a qPCR assay (S4 Fig) with primers [87] targeting the parasite glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene. All quantitative PCRs were performed on the ABI 7500 Fast Real-Time PCR system (Applied Biosystems) using the following thermo cycling conditions: initial hold at 95°C for 10 minutes, followed by 40 cycles of 95°C for 15 seconds and 60°C for 1 minute, with data collection every cycle. Specificity for the qPCR amplicons was confirmed by continuous melt curve analysis.

The relative expression of TLR was calculated using the delta-delta cycle threshold method (ΔΔCt). In order to generate values for relative expression of TLR transcript normalized to time 0-hour, the cycle threshold values (Ct) for β-actin were subtracted from TLR Ct for the same sample to generate ΔCt values for all samples. Then mean ΔCt for time 0-hour was subtracted from those of the other time points of the same snail strain and treatment to generate ΔΔCt values. Relative quantification (RQ) values were derived from ΔΔCt values using the formula 2-ΔΔCt.

BgTLR knockdown and phenotypic influence on resistance

Juvenile snails (~ 8 mm shell diameter) were injected with a cocktail of 27-mer siRNA oligonucleotides designed to specifically target 4 different regions of the BgTLR transcript. The oligonucleotide sequences were unique to the BgTLR used in this study (S2 Fig) and their specificity was confirmed through qPCR amplifying the three sequences representing transcripts with highest shared nucleotide identity with BgTLR (S5 Fig). The oligonucleotide mix was suspended in Xfect transfection reagent (Clone Tech) in order to enhance delivery and 10–20 μL was injected directly into the snail haemocoel at an approximate final concentration of 6 nM, which was determined by estimating the volume of haemolymph within the snail. Control snails received siRNA oligonucleotides targeting the green fluorescent protein (GFP). Sequences for siRNA oligonucleotides are shown in Table 2.

thumbnail
Download:
Table 2. siRNA oligonucleotide sequences.

https://doi.org/10.1371/journal.ppat.1005513.t002

SiRNA-mediated knockdown of BgTLR was demonstrated in both BS-90 and M-line snails without parasite challenge, and BS-90 that have been exposed to S. mansoni. Five snails were collected for each time point at 3, 12, 24, 48, 72 96 and 120 hours post-injection. RNA extraction, cDNA synthesis and qPCR were performed as described above. BgTLR knockdown was also confirmed at the protein level at the same time points by Western blot using both haemocytes isolated from BS-90 snails and using the B. glabrata embryonic (Bge) cell line which has haemocyte-like properties [17,88].

Phenotypic influence of BgTLR transcript knockdown was assessed in BS-90 snails. Snails were injected with the BgTLR-specific siRNA oligonucleotide mix and after 48 hours exposed individually to ~5 miracidia. Four weeks later, snails were assessed on a weekly basis, for cercariae shedding. The number of cercariae-producing infections was used as an assessment of the influence of TLR knockdown on the resistance phenotype. BgTLR knockdown snails were compared to BS-90 snails injected with the GFP-specific siRNA oligonucleotides, and to M-line snail controls challenged at the same time.

Phagocytosis assay

Five BS-90 snails each were injected with siRNA targeting BgTLR or GFP (control) and 96 hours later, haemolymph was extracted from the snails by the head-foot retraction method and immediately mixed with ~1 x 106 1μM FITC-labelled streptavidin-coated microspheres that were previously incubated with biotinylated S. mansoni excretory/secretory products and sporocysts [89] following the manufacturers protocols (Bang Laboratories Inc.). Following mixing, the haemolymph was quickly deposited on a microscope slide and allowed to sit in a humidified chamber for 3 hours. The slides were then washed thrice in 1X PBS and examined under a fluorescent microscope. Haemocytes from each snail were counted from a random field of view on the slide, and 30 haemocytes for each snail were assessed for the number of beads within each cell from which percentage phagocytosis (% of cells with one or more beads) and mean number of beads per haemocyte were calculated.

Immunocytochemistry

To visualize BgTLR in B. glabrata haemocytes and Bge cells, immunocytochemistry was done using an anti-BgTLR polyclonal antibody raised in rabbits against the extracellular domain of BgTLR (GenScript), and a goat anti-rabbit IgG secondary antibody conjugated to Alexa555 fluorophore (Life Technologies). Snails were bled by head-foot retraction and 100–200 μL of haemolymph was mixed with equal volume of 4% paraformaldehyde prepared in 1X PBS. After 10 minutes of fixation at room temperature, haemocytes were spun at 700 rpm for 5 minutes onto a coverslip and placed in the well of a 6-well plate. Haemocytes were washed with 1 mL of 1X antibody staining buffer (ASB) [0.05% sodium azide, 1% bovine serum albumin (BSA) in 1X phosphate buffered saline (PBS)], then blocked at room temperature for 1 hour in 1% BSA prepared in 1X ASB before staining using the primary antibody (suspended in blocking buffer at 1:200 dilution) at room temperature for 1 hour. Haemocytes were washed three times with 1 mL 1X ASB. Staining in secondary antibody (1:500 dilution) was performed for 1 hour at room temperature and washed as described for the primary antibody. Haemocytes were then mounted in a solution containing DAPI for 5 minutes for the staining of nuclei (GeneTex). Cells from 200 μL of confluent Bge cultures were fixed and treated using the same protocol described for the haemocytes.

Control slides were treated similarly but with primary antibody staining step omitted, or both primary and secondary antibody staining steps omitted. Observation and imaging was done under the LSM710 confocal microscope (Carl Zeiss Microimaging, Germany) at the Cross-Cancer Institute, University of Alberta. Images were processed with the accompanying ZEN 2011 software, version 7.0.0.285 and Photoshop CS5, version 12.0 x64 (Adobe Systems Incorporated, USA).

Bge cell line

The Biomphalaria glabrata embryonic (Bge) cell line is the only existing molluscan cell line. It shares behavioural and molecular attributes with Biomphalaria haemocytes such as ability to recognize and phagocytose or encapsulate foreign material including trematode targets [17,88]. Cells were cultured as described by Odoemelam et al. [90] and passaged once every 2 weeks by firm tapping or use of cell scrapper to release the cells and reseeding at 1:10 dilution.

Statistical analysis

To determine significant differences in BgTLR transcript levels, one-way analysis of variance (ANOVA) with Tukey’s post-hoc tests were performed using GraphPad Prism version 6.0f for Mac OS X (GraphPad Software, California USA, www.graphpad.com). Significant differences in the proportions of snails shedding cercariae were determined using the z-test. Statistical significance threshold was set at P ≤ 0.05.

Supporting Information

S1 Fig. Annotated BgTLR coding sequence.

The entire open reading frame is shown with alternating lines of nucleotide codons and corresponding amino acid residues. Colour codes represent the start codon (bright green), leucine-rich repeat motifs (pink), transmembrane region (grey), TIR domain (blue) and the stop codon (red).

https://doi.org/10.1371/journal.ppat.1005513.s001

(TIF)

S2 Fig. Nucleotide BLAST of siRNA positions targeted in this study.

Sequences representing transcripts with high shared nucleotide identity with BgTLR were retrieved from VectorBase (www.vectorbase.org) and used to create a custom database against which the siRNA sequences were searched. BgTLR is highlighted in grey for each siRNA target results. Percentages in brackets represent the query coverage. Note that BgTLR and BGLB008602-RA are 100% identical in all aspects except that the latter lacks the region between nucleotides 608–1489 (possibly splice variants). Most of the retrieved transcripts have nucleotide conservations of less than 50% to any of the siRNA targets. Shown in the figure is the alignment of 10 top sequences.

https://doi.org/10.1371/journal.ppat.1005513.s002

(TIF)

S3 Fig. Protein and Nucleotide BLAST of BgTLR antibody and qPCR primers respectively.

Sequences were analyzed similarly as described for S2 Fig. BgTLR is highlighted in grey for each alignment, displaying up to 10 top sequences. Custom database and BLAST analyses were done using Geneious version 6.1.6 (www.geneious.com) [91].

https://doi.org/10.1371/journal.ppat.1005513.s003

(TIF)

S4 Fig. Schistosoma mansoni GAPDH expression.

M-line (A) and BS-90 (B) snails were individually exposed to ~5 miracidia (challenged) or left unexposed (not challenged). Five and three snails respectively were collected at indicated time points over the incubation period of the parasite. RNA was extracted from whole snails, converted to cDNA and S. mansoni GAPDH expression was measured by quantitative PCR. All snails having a cycle threshold (Ct) value above zero were considered infected. RNA extracted from miracidia was used as a positive control template.

https://doi.org/10.1371/journal.ppat.1005513.s004

(TIF)

S5 Fig. siRNA-mediated knockdown of BgTLR does not affect other transcripts with a high shared nucleotide identity with BgTLR.

Quantitative PCR was performed targeting 3 transcripts that appear to be the most similar to BgTLR: BGLB008602 (A), BGLB010031 (B) and BGLB011379 (C) using cDNA generated from BS-90 BgTLR and GFP siRNA knockdown samples as templates. BgTLR siRNA did not have knockdown effect on these TLRs except on the putative splice variant (BGLB008602) which displayed a pattern similar to BgTLR knockdown.

https://doi.org/10.1371/journal.ppat.1005513.s005

(TIF)

S6 Fig. BgTLR is not expressed on all haemocytes.

Bright field view (A) of two adjacent haemocytes labelled with the nuclear stain DAPI (B) and anti-BgTLR primary antibody (C). The merged view (D) shows that BgTLR protein was only expressed on the haemocyte on the right of the panel. Scale bars represent 20 μM.

https://doi.org/10.1371/journal.ppat.1005513.s006

(TIF)

S7 Fig. BgTLR detection antibody is specific to its cognate peptide.

Protein extracts from BS-90 haemocytes were ran on duplicate SDS-PAGE gels, then transferred to nitrocellulose membranes and probed with BgTLR antibody used in this study without pre-incubation with its cognate peptide [–], pre-incubation with the peptide at 2:1 (peptide:antibody) molar ratio [3(2:1)], equal molar ratio [3(1:1)] or with an alternative BgTLR antibody targeting a different peptide [1(2:1)]. B. glabrata actin served as protein loading control. M = molecular marker.

https://doi.org/10.1371/journal.ppat.1005513.s007

(TIF)

S1 Table. Membrane-associated proteins displaying differential expression in haemocytes of BS-90 B. glabrata.

https://doi.org/10.1371/journal.ppat.1005513.s008

(DOCX)

S1 Video. Video of images acquired from multiple planes (z-stacks) supporting the description of BgTLR as plasma membrane-associated.

Samples were processed as described in materials and methods and 0.5 μM image slices were acquired under the LSM710 confocal microscope and processed into a video with the accompanying ZEN 2011 software, version 7.0.0.285.

https://doi.org/10.1371/journal.ppat.1005513.s009

(AVI)

Acknowledgments

We thank the Cross-Cancer Institute Microscopy Core Facility, University of Alberta and Valerie Phillips for technical assistance, and Michelle Gordy for helpful discussion on data presentation.

Author Contributions

Conceived and designed the experiments: EAP ALK PCH. Performed the experiments: EAP MT ALK PCH. Analyzed the data: EAP PCH MT ALK. Contributed reagents/materials/analysis tools: PCH. Wrote the paper: EAP PCH.

References

  1. 1. World Health Organization. Schistosomiasis: number of people treated worldwide in 2014. Weekly Epidemiological Record 2016 Feb; 91(5):53–60. pmid:26852429
  2. 2. Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. The Lancet Infectious Diseases 2006 7;6(7):411–425. pmid:16790382
  3. 3. Adema CM, Bayne CJ, Bridger JM, Knight M, Loker ES, Yoshino TP, et al. Will all scientists working on snails and the diseases they transmit please stand up? PLoS Neglected Tropical Diseases 2012(12).
  4. 4. Paraense WL, Correa LR. Variation in susceptibility of populations of Australorbis glabratus to a strain of Schistosoma mansoni. Rev Inst Med Trop Sao Paulo 1963 Jan-Feb;5:15–22. pmid:13941355
  5. 5. Richards CS, Merritt JW Jr. Genetic factors in the susceptibility of juvenile Biomphalaria glabrata to Schistosoma mansoni infection. Am J Trop Med Hyg 1972 Jul;21(4):425–434. pmid:5050093
  6. 6. Newton WL. The establishment of a strain of Australorbis glabratus which combines albinism and high susceptibility to infection with Schistosoma mansoni. J Parasitol 1955 Oct;41(5):526–528. pmid:13264025
  7. 7. Cooper LA, Richards CS, Lewis FA, Minchella DJ. Schistosoma mansoni: relationship between low fecundity and reduced susceptibility to parasite infection in the snail Biomphalaria glabrata. Exp Parasitol 1994 Aug;79(1):21–28. pmid:8050522
  8. 8. Bayne CJ. Successful parasitism of vector snail Biomphalaria glabrata by the human blood fluke (trematode) Schistosoma mansoni: a 2009 assessment. Mol Biochem Parasitol 2009 May;165(1):8–18. pmid:19393158
  9. 9. Zhang SM, Coultas KA. Identification and characterization of five transcription factors that are associated with evolutionarily conserved immune signaling pathways in the schistosome-transmitting snail Biomphalaria glabrata. Mol Immunol 2011 Sep;48(15–16):1868–1881. pmid:21696828
  10. 10. Knight M, Arican-Goktas HD, Ittiprasert W, Odoemelam EC, Miller AN, Bridger JM. Schistosomes and snails: a molecular encounter. Front Genet 2014 Jul 21;5:230. pmid:25101114
  11. 11. Guillou F, Roger E, Mone Y, Rognon A, Grunau C, Theron A, et al. Excretory-secretory proteome of larval Schistosoma mansoni and Echinostoma caproni, two parasites of Biomphalaria glabrata. Mol Biochem Parasitol 2007 Sep;155(1):45–56. pmid:17606306
  12. 12. Bouchut A, Sautiere PE, Coustau C, Mitta G. Compatibility in the Biomphalaria glabrata/Echinostoma caproni model: Potential involvement of proteins from hemocytes revealed by a proteomic approach. Acta Trop 2006 Jul;98(3):234–246. pmid:16792992
  13. 13. Vergote D, Bouchut A, Sautiere PE, Roger E, Galinier R, Rognon A, et al. Characterisation of proteins differentially present in the plasma of Biomphalaria glabrata susceptible or resistant to Echinostoma caproni. Int J Parasitol 2005 Feb;35(2):215–224. pmid:15710442
  14. 14. Roger E, Mitta G, Mone Y, Bouchut A, Rognon A, Grunau C, et al. Molecular determinants of compatibility polymorphism in the Biomphalaria glabrata/Schistosoma mansoni model: new candidates identified by a global comparative proteomics approach. Mol Biochem Parasitol 2008 Feb;157(2):205–216. pmid:18083248
  15. 15. Raghavan N, Miller AN, Gardner M, FitzGerald PC, Kerlavage AR, Johnston DA, et al. Comparative gene analysis of Biomphalaria glabrata hemocytes pre- and post-exposure to miracidia of Schistosoma mansoni. Mol Biochem Parasitol 2003 Feb;126(2):181–191. pmid:12615317
  16. 16. Mitta G, Galinier R, Tisseyre P, Allienne JF, Girerd-Chambaz Y, Guillou F, et al. Gene discovery and expression analysis of immune-relevant genes from Biomphalaria glabrata hemocytes. Dev Comp Immunol 2005;29(5):393–407. pmid:15707661
  17. 17. Coustau C, Mitta G, Dissous C, Guillou F, Galinier R, Allienne JF, et al. Schistosoma mansoni and Echinostoma caproni excretory-secretory products differentially affect gene expression in Biomphalaria glabrata embryonic cells. Parasitology 2003 Dec;127(Pt 6):533–542. pmid:14700189
  18. 18. Guillou F, Mitta G, Galinier R, Coustau C. Identification and expression of gene transcripts generated during an anti-parasitic response in Biomphalaria glabrata. Dev Comp Immunol 2007;31(7):657–671. pmid:17166585
  19. 19. Guillou F, Mitta G, Dissous C, Pierce R, Coustau C. Use of individual polymorphism to validate potential functional markers: case of a candidate lectin (BgSel) differentially expressed in susceptible and resistant strains of Biomphalaria glabrata. Comp Biochem Physiol B Biochem Mol Biol 2004 Jun;138(2):175–181. pmid:15193273
  20. 20. Oliveira AL, Da Silva D, Zanotti-Magalhaes EM, Abdel-Hamid AZ, Ribeiro-Paes JT. Schistosome/mollusk: genetic compatibility. Genet Mol Res 2008;7(2):518–526. pmid:18561386
  21. 21. Webster JP, Davies CM, Ndamba J, Noble LR, Jones CS, Woolhouse ME. Spatio-temporal genetic variability in the schistosome intermediate host Biomphalaria pfeifferi. Ann Trop Med Parasitol 2001 Jul;95(5):515–527. pmid:11487373
  22. 22. Zavodna M, Sandland GJ, Minchella DJ. Effects of intermediate host genetic background on parasite transmission dynamics: a case study using Schistosoma mansoni. Exp Parasitol 2008 Sep;120(1):57–61. pmid:18538767
  23. 23. Thiele EA, Sorensen RE, Gazzinelli A, Minchella DJ. Genetic diversity and population structuring of Schistosoma mansoni in a Brazilian village. Int J Parasitol 2008 Mar;38(3–4):389–399. pmid:17825305
  24. 24. Theron A, Sire C, Rognon A, Prugnolle F, Durand P. Molecular ecology of Schistosoma mansoni transmission inferred from the genetic composition of larval and adult infrapopulations within intermediate and definitive hosts. Parasitology 2004 Nov;129(Pt 5):571–585. pmid:15552402
  25. 25. Zahoor Z, Lockyer AE, Davies AJ, Kirk RS, Emery AM, Rollinson D, et al. Differences in the gene expression profiles of haemocytes from schistosome-susceptible and -resistant Biomphalaria glabrata exposed to Schistosoma mansoni excretory-secretory products. PLoS One 2014 Mar 24;9(3):e93215. pmid:24663063
  26. 26. Lockyer AE, Spinks JN, Walker AJ, Kane RA, Noble LR, Rollinson D, et al. Biomphalaria glabrata transcriptome: identification of cell-signalling, transcriptional control and immune-related genes from open reading frame expressed sequence tags (ORESTES). Dev Comp Immunol 2007;31(8):763–782. pmid:17208299
  27. 27. Lockyer AE, Emery AM, Kane RA, Walker AJ, Mayer CD, Mitta G, et al. Early Differential Gene Expression in Haemocytes from Resistant and Susceptible Biomphalaria glabrata Strains in Response to Schistosoma mansoni. PLoS One 2012 Dec 26;7(12):.
  28. 28. Lockyer AE, Noble LR, Rollinson D, Jones CS. Schistosoma mansoni: resistant specific infection-induced gene expression in Biomphalaria glabrata identified by fluorescent-based differential display. Exp Parasitol 2004 May-Jun;107(1–2):97–104. pmid:15208043
  29. 29. Lockyer AE, Spinks J, Noble LR, Rollinson D, Jones CS. Identification of genes involved in interactions between Biomphalaria glabrata and Schistosoma mansoni by suppression subtractive hybridization. Mol Biochem Parasitol 2007 Jan;151(1):18–27. pmid:17081633
  30. 30. Hanington PC, Lun CM, Adema CM, Loker ES. Time series analysis of the transcriptional responses of Biomphalaria glabrata throughout the course of intramolluscan development of Schistosoma mansoni and Echinostoma paraensei. Int J Parasitol 2010 Jun;40(7):819–831. pmid:20083115
  31. 31. Hanington PC, Forys MA, Dragoo JW, Zhang S, Adema CM, Loker ES, et al. Role for a somatically diversified lectin in resistance of an invertebrate to parasite infection. Proc Natl Acad Sci U S A 2010(49):21087. pmid:21084634
  32. 32. Hanington PC, Forys MA, Loker ES. A Somatically Diversified Defense Factor, FREP3, is a Determinant of Snail Resistance to Schistosome Infection. PLoS Negl Trop Dis 2012 Mar 27;6(3):.
  33. 33. Adema CM, Hertel LA, Miller RD, Loker ES. A family of fibrinogen-related proteins that precipitates parasite-derived molecules is produced by an invertebrate after infection. Proc Natl Acad Sci U S A 1997 Aug 5;94(16):8691–8696. pmid:9238039
  34. 34. Moné Y, Gourbal B, Duval D, Du Pasquier L, Kieffer-Jaquinod S, Mitta G. A Large Repertoire of Parasite Epitopes Matched by a Large Repertoire of Host Immune Receptors in an Invertebrate Host/Parasite Model. PLoS Negl Trop Dis 2010 Sep 7;4(9):.
  35. 35. Baeza Garcia A, Pierce RJ, Gourbal B, Werkmeister E, Colinet D, Reichhart JM, et al. Involvement of the cytokine MIF in the snail host immune response to the parasite Schistosoma mansoni. PLoS Pathog 2010 Sep 23;6(9):e1001115. pmid:20886098
  36. 36. Goodall CP, Bender RC, Broderick EJ, Bayne CJ. Constitutive differences in Cu/Zn superoxide dismutase mRNA levels and activity in hemocytes of Biomphalaria glabrata (Mollusca) that are either susceptible or resistant to Schistosoma mansoni (Trematoda). Molecular & Biochemical Parasitology 2004;137:321–328.
  37. 37. Bender RC, Broderick EJ, Goodall CP, Bayne CJ. Respiratory burst of Biomphalaria glabrata hemocytes: Schistosoma mansoni-resistant snails produce more extracellular H2O2 than susceptible snails. J Parasitol 2005 Apr;91(2):275–279. pmid:15986600
  38. 38. Goodall CP, Bender RC, Brooks JK, Bayne CJ. Biomphalaria glabrata cytosolic copper/zinc superoxide dismutase (SOD1) gene: association of SOD1 alleles with resistance/susceptibility to Schistosoma mansoni. Mol Biochem Parasitol 2006 Jun;147(2):207–210. pmid:16564582
  39. 39. Galinier R, Portela J, Moné Y, Allienne JF, Henri H, Delbecq S, et al. Biomphalysin, a New β Pore-forming Toxin Involved in Biomphalaria glabrata Immune Defense against Schistosoma mansoni. PLoS Pathog 2013 Mar;9(3):. Epub 2013 Mar 21
  40. 40. Plows LD, Cook RT, Davies AJ, Walker AJ. Integrin engagement modulates the phosphorylation of focal adhesion kinase, phagocytosis, and cell spreading in molluscan defence cells. Biochim Biophys Acta 2006 Aug;1763(8):779–786. pmid:16766054
  41. 41. Humphries JE, Yoshino TP. Schistosoma mansoni excretory-secretory products stimulate a p38 signalling pathway in Biomphalaria glabrata embryonic cells. Int J Parasitol 2006 Jan;36(1):37–46. pmid:16194541
  42. 42. Walker AJ. Do trematode parasites disrupt defence-cell signalling in their snail hosts? Trends Parasitol 2006 Apr;22(4):154–159. pmid:16488192
  43. 43. Salamat Z, Sullivan JT. Involvement of protein kinase C signalling and mitogen-activated protein kinase in the amebocyte-producing organ of Biomphalaria glabrata (Mollusca). Dev Comp Immunol 2009 Jun;33(6):725–727. pmid:19183562
  44. 44. Zhang SM, Zeng Y, Loker ES. Characterization of immune genes from the schistosome host snail Biomphalaria glabrata that encode peptidoglycan recognition proteins and gram-negative bacteria binding protein. Immunogenetics 2007 Nov;59(11):883–898. pmid:17805526
  45. 45. Bowie A, O'Neill LA. The interleukin-1 receptor/Toll-like receptor superfamily: signal generators for pro-inflammatory interleukins and microbial products. J Leukoc Biol 2000 Apr;67(4):508–514. pmid:10770283
  46. 46. Hoshino K, Takeuchi O, Kawai T, Sanjo H, Ogawa T, Takeda Y, et al. Cutting edge: Toll-like receptor 4 (TLR4)-deficient mice are hyporesponsive to lipopolysaccharide: evidence for TLR4 as the Lps gene product. J Immunol 1999 Apr 1;162(7):3749–3752. pmid:10201887
  47. 47. Hayashi F, Smith KD, Ozinsky A, Hawn TR, Yi EC, Goodlett DR, et al. The innate immune response to bacterial flagellin is mediated by Toll-like receptor 5. Nature 2001 Apr 26;410(6832):1099–1103. pmid:11323673
  48. 48. Hemmi H, Takeuchi O, Kawai T, Kaisho T, Sato S, Sanjo H, et al. A Toll-like receptor recognizes bacterial DNA. Nature 2000 Dec 7;408(6813):740–745. pmid:11130078
  49. 49. Hochrein H, Schlatter B, O'Keeffe M, Wagner C, Schmitz F, Schiemann M, et al. Herpes simplex virus type-1 induces IFN-alpha production via Toll-like receptor 9-dependent and -independent pathways. Proc Natl Acad Sci U S A 2004 Aug 3;101(31):11416–11421. pmid:15272082
  50. 50. Takahara K, Yashima Y, Omatsu Y, Yoshida H, Kimura Y, Kang YS, et al. Functional comparison of the mouse DC-SIGN, SIGNR1, SIGNR3 and Langerin, C-type lectins. Int Immunol 2004 Jun;16(6):819–829. pmid:15096474
  51. 51. Coban C, Ishii KJ, Kawai T, Hemmi H, Sato S, Uematsu S, et al. Toll-like receptor 9 mediates innate immune activation by the malaria pigment hemozoin. J Exp Med 2005 Jan 3;201(1):19–25. pmid:15630134
  52. 52. Bryant CE, Gay NJ, Heymans S, Sacre S, Schaefer L, Midwood KS. Advances in Toll-like receptor biology: Modes of activation by diverse stimuli. Crit Rev Biochem Mol Biol 2015 Apr 10:1–21.
  53. 53. Michel T, Reichhart JM, Hoffmann JA, Royet J. Drosophila Toll is activated by Gram-positive bacteria through a circulating peptidoglycan recognition protein. Nature 2001 Dec 13;414(6865):756–759. pmid:11742401
  54. 54. Akira S, Uematsu S, Takeuchi O. Pathogen recognition and innate immunity. Cell 2006 Feb 24;124(4):783–801. pmid:16497588
  55. 55. Bender RC, Goodall CP, Blouin MS, Bayne CJ. Variation in expression of Biomphalaria glabrata SOD1: a potential controlling factor in susceptibility/resistance to Schistosoma mansoni. Dev Comp Immunol 2007; 31(9):874–878. pmid:17292470
  56. 56. Bonner KM, Bayne CJ, Larson MK, Blouin MS. Effects of Cu/Zn superoxide dismutase (sod1) genotype and genetic background on growth, reproduction and defense in Biomphalaria glabrata. PLoS Negl Trop Dis 2012; 6(6): e1701. pmid:22724037
  57. 57. Tennessen JA, Théron A, Marine M, Yeh JY, Rognon A, Blouin MS. Hyperdiverse gene cluster in snail host conveys resistance to human schistosome parasites. PLoS Genet 2015; 11(3): e1005067. pmid:25775214
  58. 58. Offord V, Werling D. LRRfinder2.0: a webserver for the prediction of leucine-rich repeats. Innate Immun 2013;19(4):398–402. pmid:23178228
  59. 59. de Castro E, Sigrist CJ, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, et al. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res 2006 Jul 1;34(Web Server issue):W362–5. pmid:16845026
  60. 60. Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schaffer AA, et al. Protein database searches using compositionally adjusted substitution matrices. FEBS J 2005 Oct;272(20):5101–5109. pmid:16218944
  61. 61. Latz E, Visintin A, Lien E, Fitzgerald KA, Monks BG, Kurt-Jones E, et al. Lipopolysaccharide rapidly traffics to and from the golgi apparatus with the toll-like receptor 4-MD-2-CD14 complex in a process that is distinct from the initiation of signal transduction. J Biol Chem 2002;277(49):47834–47843. pmid:12324469
  62. 62. World Health Organization. Accelerating work to overcome the global impact of neglected tropical diseases—a roadmap for implementation. Geneva: Department of Control of Neglected Tropical Diseases. 2012; Available at: http://www.who.int/entity/neglected_diseases/NTD_RoadMap_2012_Fullversion.pdf. Accessed June 24, 2015.
    • 63. Sullivan JT, Spence JV, Nunez JK. Killing of Schistosoma mansoni sporocysts in Biomphalaria glabrata implanted with amoebocyte-producing organ allografts from resistant snails. J Parasitol 1995 Oct;81(5):829–833. pmid:7472893
    • 64. Hahn UK, Bender RC, Bayne CJ. Killing of Schistosoma mansoni sporocysts by hemocytes from resistant Biomphalaria glabrata: role of reactive oxygen species. J Parasitol 2001 Apr;87(2):292–299. pmid:11318558
    • 65. Lie KJ, Jeong KH, Heyneman D. Tissue reactions induced by Schistosoma mansoni in Biomphalaria glabrata. Ann Trop Med Parasitol 1980 Apr;74(2):157–166. pmid:7436601
    • 66. Elvitigala DA, Premachandra HK, Whang I, Nam BH, Lee J. Molecular insights of the first gastropod TLR counterpart from disk abalone (Haliotis discus discus), revealing its transcriptional modulation under pathogenic stress. Fish Shellfish Immunol 2013 Aug;35(2):334–342. pmid:23669649
    • 67. Qiu L, Song L, Xu W, Ni D, Yu Y. Molecular cloning and expression of a Toll receptor gene homologue from Zhikong Scallop, Chlamys farreri. Fish Shellfish Immunol 2007 May;22(5):451–466. pmid:17158063
    • 68. Mateo DR, Greenwood SJ, Araya MT, Berthe FC, Johnson GR, Siah A. Differential gene expression of gamma-actin, Toll-like receptor 2 (TLR-2) and interleukin-1 receptor-associated kinase 4 (IRAK-4) in Mya arenaria haemocytes induced by in vivo infections with two Vibrio splendidus strains. Dev Comp Immunol 2010 Jul;34(7):710–714. pmid:20156478
    • 69. Toubiana M, Gerdol M, Rosani U, Pallavicini A, Venier P, Roch P. Toll-like receptors and MyD88 adaptors in Mytilus: complete cds and gene expression levels. Dev Comp Immunol 2013 Jun;40(2):158–166. pmid:23485525
    • 70. Zhang Y, He X, Yu F, Xiang Z, Li J, Thorpe KL, et al. Characteristic and Functional Analysis of Toll-like Receptors (TLRs) in the lophotrocozoan, Crassostrea gigas, Reveals Ancient Origin of TLR-Mediated Innate Immunity. PLoS ONE 2013;8(10):e76464. pmid:24098508
    • 71. Lehr T, Geyer H, Maass K, Doenhoff MJ, Geyer R. Structural characterization of N-glycans from the freshwater snail Biomphalaria glabrata cross-reacting with Schistosoma mansoni glycoconjugates. Glycobiology 2007 Jan;17(1):82–103. pmid:16971380
    • 72. Lehr T, Beuerlein K, Doenhoff MJ, Grevelding CG, Geyer R. Localization of carbohydrate determinants common to Biomphalaria glabrata as well as to sporocysts and miracidia of Schistosoma mansoni. Parasitology 2008 Jul;135(8):931–942. pmid:18507884
    • 73. Peterson NA, Hokke CH, Deelder AM, Yoshino TP. Glycotope analysis in miracidia and primary sporocysts of Schistosoma mansoni: differential expression during the miracidium-to-sporocyst transformation. Int J Parasitol 2009 Oct;39(12):1331–1344. pmid:19545571
    • 74. Wu XJ, Sabat G, Brown JF, Zhang M, Taft A, Peterson N, et al. Proteomic analysis of Schistosoma mansoni proteins released during in vitro miracidium-to-sporocyst transformation. Mol Biochem Parasitol 2009 Mar;164(1):32–44. pmid:19095013
    • 75. Hahn UK, Bender RC, Bayne CJ. Production of reactive oxygen species by hemocytes of Biomphalaria glabrata: carbohydrate-specific stimulation. Dev Comp Immunol 2000 Sep-Oct;24(6–7):531–541. pmid:10831788
    • 76. Lacchini AH, Davies AJ, Mackintosh D, Walker AJ. Beta-1, 3-glucan modulates PKC signalling in Lymnaea stagnalis defence cells: a role for PKC in H2O2 production and downstream ERK activation. J Exp Biol 2006 Dec;209(Pt 24):4829–4840. pmid:17142672
    • 77. Zelck UE, Gege BE, Schmid S. Specific inhibitors of mitogen-activated protein kinase and PI3-K pathways impair immune responses by hemocytes of trematode intermediate host snails. Dev Comp Immunol 2007;31(4):321–331. pmid:16926049
    • 78. Humphries JE, Yoshino TP. Regulation of hydrogen peroxide release in circulating hemocytes of the planorbid snail Biomphalaria glabrata. Dev Comp Immunol 2008;32(5):554–562. pmid:17981329
    • 79. Fryer SE, Dykes-Hoberg M, Bayne CJ. Changes in plasma opsonization of yeast after isolation of Biomphalaria glabrata in small volumes of water. J Invertebr Pathol 1989;54(2):275–276. pmid:2768831
    • 80. Johnston LA, Yoshino TP. Analysis of lectin- and snail plasma-binding glycopeptides associated with the tegumental surface of the primary sporocysts of Schistosoma mansoni. Parasitology 1996 May;112 (Pt 5)(Pt 5):469–479. pmid:8677136
    • 81. Johnston LA, Yoshino TP. Larval Schistosoma mansoni excretory-secretory glycoproteins (ESPs) bind to hemocytes of Biomphalaria glabrata (Gastropoda) via surface carbohydrate binding receptors. J Parasitol 2001 Aug;87(4):786–793. pmid:11534642
    • 82. Castillo MG, Yoshino TP. Carbohydrate inhibition of Biomphalaria glabrata embryonic (Bge) cell adhesion to primary sporocysts of Schistosoma mansoni. Parasitology 2002 Dec;125(Pt 6):513–525. pmid:12553570
    • 83. Castillo MG, Wu XJ, Dinguirard N, Nyame AK, Cummings RD, Yoshino TP. Surface membrane proteins of Biomphalaria glabrata embryonic cells bind fucosyl determinants on the tegumental surface of Schistosoma mansoni primary sporocysts. J Parasitol 2007 Aug;93(4):832–840. pmid:17918362
    • 84. Smiley ST, King JA, Hancock WW. Fibrinogen stimulates macrophage chemokine secretion through toll-like receptor 4. J Immunol 2001 Sep 1;167(5):2887–2894. pmid:11509636
    • 85. Yoshino TP, Coustau C. Immunobiology of Biomphalaria–Trematode Interactions. In: Toledo R, Fried B, editors. Biomphalaria Snails and Larval Trematodes. New York: Springer; 2011:159–180.
      • 86. Humphries JE, Elizondo L, Yoshino TP. Protein kinase C regulation of cell spreading in the molluscan Biomphalaria glabrata embryonic (Bge) cell line. Biochim Biophys Acta 2001 Sep 26;1540(3):243–252. pmid:11583819
      • 87. Boyle JP, Wu XJ, Shoemaker CB, Yoshino TP. Using RNA interference to manipulate endogenous gene expression in Schistosoma mansoni sporocysts. Mol Biochem Parasitol 2003 May;128(2):205–215. pmid:12742587
      • 88. Yoshino TP, Coustau C, Modat S, Castillo MG. The Biomphalaria glabrata Embryonic (BGE) Molluscan Cell Line: Establishment of an In Vitro Cellular Model for the Study of Snail Host-Parasite Interactions. Malacologia-Philadelphia 1999;41:331–344.
      • 89. Zahoor Z, Davies AJ, Kirk RS, David Rollinson D,Walker AJ. Nitric oxide production by Biomphalaria glabrata haemocytes: effects of Schistosoma mansoni ESPs and regulation through the extracellular signal-regulated kinase pathway. Parasit Vectors 2009; 2(1): 18. pmid:19386102
      • 90. Odoemelam E, Raghavan N, Miller A, Bridger JM, Knight M. Revised karyotyping and gene mapping of the Biomphalaria glabrata embryonic (Bge) cell line. Int J Parasitol 2009;39:675–681. pmid:19133265
      • 91. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012 06/15;28(12):1647–1649. pmid:22543367