# A shift from cellular to humoral responses contributes to innate immune memory in the vector snail

## *Biomphalaria glabrata*

Silvain Pinaud, Julien Portela, David Duval, Fanny Nowacki, Marie-Aude Olive, Jean-François Allienne, Richard Galinier, Nolwenn M. Dheilly, Sylvie Kieffer-Jaquinod, Guillaume Mitta, et al.

<figure><img src="image_1.png" /></figure>

### ► To cite this version:

Silvain Pinaud, Julien Portela, David Duval, Fanny Nowacki, Marie-Aude Olive, et al. A shift from cellular to humoral responses contributes to innate immune memory in the vector snail *Biomphalaria glabrata*. PLoS Pathogens, 2016, 12 (1), pp.e1005361. 10.1371/journal.ppat.1005361. hal-03778962

HAL Id: hal-03778962

https://hal.science/hal-03778962

Submitted on 20 Sep 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

CrossMark

click for updates

RESEARCH ARTICLE

# A Shift from Cellular to Humoral Responses Contributes to Innate Immune Memory in the Vector Snail *Biomphalaria glabrata*

Silvain Pinaud¹, Julien Portela¹, David Duval¹, Fanny C. Nowacki¹,², Marie-Aude Olive¹, Jean-François Allienne¹, Richard Galinier¹, Nolwenn M. Dheilly¹,³, Sylvie Kieffer-Jaquinod⁴, Guillaume Mitta¹, André Théron¹, Benjamin Gourbal¹*

1 University of Perpignan Via Domitia, IHPE UMR 5244, CNRS, IFREMER, University of Montpellier, Perpignan, France, 2 Institute of Biological, Environmental & Rural Sciences, Aberystwyth University, Aberystwyth Ceredigion, United Kingdom, 3 School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, New York, United States of America, 4 Plate-forme d'analyses protéomiques EDyP-Service, Laboratoire de Biologie à Grande Echelle UMR-S 1038 Inserm/CEA/UJF CEA, Grenoble, France

* benjamin.gourbal@univ-perp.fr

## Abstract

Discoveries made over the past ten years have provided evidence that invertebrate anti-parasitic responses may be primed in a sustainable manner, leading to the failure of a secondary encounter with the same pathogen. This phenomenon called "immune priming" or "innate immune memory" was mainly phenomenological. The demonstration of this process remains to be obtained and the underlying mechanisms remain to be discovered and exhaustively tested with rigorous functional and molecular methods, to eliminate all alternative explanations. In order to achieve this ambitious aim, the present study focuses on the Lophotrochozoan snail, *Biomphalaria glabrata*, in which innate immune memory was recently reported. We provide herein the first evidence that a shift from a cellular immune response (encapsulation) to a humoral immune response (biomphalysin) occurs during the development of innate memory. The molecular characterisation of this process in *Biomphalaria/Schistosoma* system was undertaken to reconcile mechanisms with phenomena, opening the way to a better comprehension of innate immune memory in invertebrates. This prompted us to revisit the artificial dichotomy between innate and memory immunity in invertebrate systems.

## Author Summary

Schistosomiasis is the second most widespread tropical parasitic disease after malaria. It is caused by flatworms of the genus *Schistosoma*. Its life cycle is complex and requires certain freshwater snail species as the intermediate host. Given the limited options for treating S. mansoni infections, much research has focused on a better understanding of the immunobiological interactions between the invertebrate host *Biomphalaria glabrata* and its parasite S. mansoni. Recently, we demonstrated the existence of a time-dependent and

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Pinaud S, Portela J, Duval D, Nowacki FC, Olive M-A, Allienne J-F, et al. (2016) A Shift from Cellular to Humoral Responses Contributes to Innate Immune Memory in the Vector Snail *Biomphalaria glabrata*. PLoS Pathog 12(1): e1005361. doi:10.1371/journal.ppat.1005361

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

**Received:** September 15, 2015

**Accepted:** December 4, 2015

**Published:** January 6, 2016

**Copyright:** © 2016 Pinaud 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 Statement:** All relevant data are within the paper and its Supporting Information files. Transcriptome is available on IHPE website bioinformatics resources (http://2ei.univ-perp.fr/?page_id=89).

**Funding:** BG acknowledges support from ANR JCJC INVIMORY (number ANR-13-JSV7-0009). 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.

genotype-dependent acquired innate immune memory in *B. glabrata* snails. A primo-infection of the Lophotrochozoan vector snail, *Biomphalaria glabrata*, with *Schistosoma mansoni* totally protected the snail against a secondary challenge. Learning more about the immunobiological interactions between *B. glabrata* and *S. mansoni* could have important socioeconomic and public health impacts by changing the way we attempt to eradicate parasitic diseases and prevent or control Schistosomiasis in the field.

## Introduction

The environment of an invertebrate is filled with complex and changing populations of micro-organisms, including potential pathogens. This engenders selective pressures comparable to those experienced by Gnathostomes [1], and means that invertebrates should possess sophisticated immune systems capable of dealing with such pathogens. Indeed, recent studies have shown that the immune defenses of invertebrates are more complex and specific than previously thought, and the existence of innate immune memory or priming has been suggested [2–6].

To date, the observations of invertebrate innate immune memory have been mainly phenomenological and based on ecological or phenotypic studies, and little work has addressed the potential molecular and cellular mechanisms underlying these processes. The innate immune system of invertebrates includes the barrier functions of the epithelium, the cellular immune response and the humoral defense response [7,8]: the former refers to mucosal immunity and hemocyte responses (phagocytosis, encapsulation, melanization), while the latter includes pathogen recognition receptors (PRRs), antimicrobial peptides, coagulation, and the production of cytolytic molecules or reactive intermediates of oxygen and nitrogen [7,8]. During the primary infection of an invertebrate, the pathogen is recognized at the first encounter and the cellular and humoral defense responses are coordinated to neutralize the intruder. To the best of our knowledge, the existing studies investigating the molecular mechanisms of innate immune memory in invertebrates have all suggested that the cellular immune response and/or hemocyte phagocytosis is/are improved upon a subsequent encounter with the same pathogen [9–11]. For example, in *Porcellio scaber*, enhanced phagocytosis was demonstrated after a first encounter with different bacteria species (*Bacillus thuringiensis/Escherichia coli*) [12]. For *Anopheles gambiae*, hemocytic differentiation evidenced by increased mRNA levels of hemocyte-specific genes, such as thioester-containing protein 1 (TEP1) and leucine rich repeat immune protein 1 (LRIM1) was reported following initial exposures to *Plasmodium* and bacteria [5]. In shrimp, bacterial challenge was followed by an enhanced cellular immunity characterized by a significant increase in the percentage of phagocytic cells [13]. The innate immune memory response in invertebrates has been previously described as being involved in two mechanisms, namely:

i. a process of acquired resistance or sustained response that consists of a long-lasting protection against a later challenge that persists even if the pathogen is neutralized; or

ii. a recalled response that consists of the ability to store information of previously met pathogens and recall it later to generate a faster and more powerful response against a subsequent exposure to the same pathogen.

However, the international community working on invertebrate innate immunity believes that the existing innate immune memory observations cannot be used in isolation. Instead, they should be used to construct hypotheses that need now to be exhaustively tested with
rigorous functional, cellular, biochemical and molecular methods, until all alternative explanations are eliminated [10,14]. To fully understand the capabilities of invertebrate immune systems, we must use global molecular approaches at the whole-organism level to investigate the mechanisms that form the basis for innate immune memory, while being guided by the existing observations and phenotypic descriptions. Without this demonstration, the existence of immune memory process in invertebrates will remain controversial and doubted by many immunologists [for the polemic, see [14,15]].

Recently, we demonstrated that a primo-infection of the Lophotrochozoan vector snail, *Biomphalaria glabrata*, with *Schistosoma mansoni* protected completely the snail against a homologous secondary challenge [4]. Total protection was reached 10 days after primo-infection, and was maintained for the rest of the snail's lifespan. Our findings provided evidence for the existence of a time-dependent acquired innate immune memory in *B. glabrata* snails [4,16]. Moreover, we used homologous and heterologous challenges to demonstrate that this innate immune memory was genotype-dependent, in that the protection decreased with increasing neutral genetic distance between the parasites used for the primo-infection and the secondary challenge [4].

Here, we sought to investigate the molecular mechanisms of innate immune memory in *B. glabrata* in response to *S. mansoni* infections. We precisely describe the snail's immune response phenotypes, then characterize the molecular determinants involved in *B. glabrata* innate immune memory via a global transcriptomic approach and a targeted proteomic analysis of plasmatic factors performed in combination with RNA interference. This work provides the first global investigation of the molecular processes supporting innate immune memory in an invertebrate model. Moreover, *Biomphalaria* snails have an important role in the transmission of Schistosomiasis the second most widespread human parasitic disease after malaria causing 200,000 deaths annually. More emphasis on snail-related research, on the role of snails and parasite intramolluscan larval stages in transmission is essential for changing the way we attempt to eradicate parasitic diseases [17,18]. Learning more about the immunobiological interactions between *B. glabrata* and *S. mansoni* could have important socioeconomic and public health impacts by contributing to the discovery of new ways to prevent and/or control Schistosomiasis diseases by limiting the parasite in the field.

## Materials and Methods

Further details are provided in S1 Appendix and Fig 1.

### Snail and parasite strains

The present work utilized a strain of *Biomphalaria glabrata* originating from Brazil (BgBRE), along with its 100% compatible sympatric strain of *Schistosoma mansoni* (SmBRE). Both were as previously described [19].

### Innate immune memory protocol

One hundred and forty BgBRE snails were primo-infected with 10 miracidia of SmBRE each; 60 were secondarily challenged 25 days later with 10 miracidia of SmBRE each. For the RNA-seq approach, pools of 20 BgBRE snails were recovered at 1, 4, 15, and 25 days post primo infection (DPPI); these were designated 1DPPI, 4DPPI, 15DPPI and 25DPPI, respectively. Twenty more snails were recovered at each of 1, 4 and 15 days after the secondary challenge; equimolar amounts of each of these experimental groups were joined together into a single sample designated as days post-secondary challenge (DPC) sample. Two pools of 20 uninfected snails (designated naive1 and naive2) were sampled and used as control replicates. For the
<figure><img src="image_5.png" /><figcaption>**Fig 1. Overview of experimental procedures.** Innate immune memory experiments were carried out. For primo-infection, Brazilian *Biomphalaria glabrata* (BgBRE) snails were individually exposed to 10 miracidia of their sympatric Brazilian *Schistosoma mansoni* trematode parasite (SmBRE). Following infection depending on the compatibility status of the snail/parasite couples, some of the miracidia were encapsulated by the hemocytes (snail immune cells) or developed into primary sporocysts (intra-molluscan stage of the parasite). Intramolluscan parasite stages include two generations of sporocysts (primary sporocyst (SPI) and secondary sporocyst (SPII)) and the production of cercariae. SPII developed inside SPI and migrated to reach the snail hepato-pancrea. Cercariae developed inside SPI and migrated back into the snail to reach the aquatic environment. Twenty-five days after primo-infection, the snails were challenged for a second time with again 10 SmBRE miracidia. In this case all miracidia degenerated in snail tissues, demonstrating the activation of a humoral immune response. Immune phenotypes observed during innate immune memory process were analyzed using a histological approach (see Fig 2A, 2B & 2C). In order to explore the molecular mechanisms of innate immune memory several experimental procedures were designed. A RNAseq experiment was realized with samples recovered from uninfected snails (Naive 1, Naive 2), samples recovered at 1, 4, 15 and 25 days post primo-infection (DPPI) and at 1, 4 and 15 days post-secondary challenge (DPC) (see Fig 3). Based on RNaseq results, functional validation of the FREP immune recognition receptor was undergone. First, individual quantification were made for all FREPs annotated on transcriptomic analysis (see Fig 4A). FREP knockdown was then carried out by siRNA injection, normalized by siGFP and monitored by Q-RT-PCR (see Fig 4B & 4C). Finally, to confirm the involvement of plasmatic factors in innate immune memory, snail hemolymph was recovered (Naive, 15, 25 DPPI and 15 DPC) and plasmatic fraction was characterized by 2D-gel electrophoresis (see Fig 5A & 5B). Plasma samples were also injected to naïve snails to demonstrate that immune protection could be acquired following primed snail plasma transfer (see Fig 5C).</figcaption></figure>

doi:10.1371/journal.ppat.1005361.g001

proteomic approach, 200 BgBRE snails were primo-infected and secondarily challenged with 10 miracidia (per round) of SmBRE. Plasma sampling was performed on the same schedule described for the RNAseq experiments. Five biological replicates (5 independent pooled plasmas from 10 snails) were recovered from naïve snails (controls), and 15DPPI, 25DPPI and 15DPC snails.
## Histological procedures

Forty-eight hours after primo-infection or secondary challenge with 10 miracidia of SmBRE, snails were fixed in Halmi's fixative (mercuric chloride 4.5%, sodium chloride 0.5%, trichloroacetic acid 2%, formol 20%, acetic acid 4% and 10% of picric acid-saturated aqueous solution). Fixed mollusks were dehydrated and embedded in paraffin as previously described [20,21]. Transverse histological sections (10-µm thick) were cut and stained using azocarmine G and Heidenhain's azan. The following serial steps were used: (i) re-hydration (toluene, 95, 70, 30% ethanol and distilled water); (ii) coloration (azocarmine G, 70% ethanol / 1% aniline, 1% acetic alcohol, distilled water, 5% phosphotungstic acid, distilled water, Heidenhain's azan) and (iii) dehydration (95% ethanol, absolute ethanol, toluene). The preparations were then mounted with Entellan (Sigma Life Science, St. Louis Missouri, USA) and subjected to microscopic examination.

## Whole-body RNA sequencing

Total RNA was extracted from naive1, naive2, 1DPPI, 4DPPI, 15DPPI, 25DPPI and DPC samples using nitrogen and the TRIzol reagent (Sigma Life Science, St. Louis Missouri, USA). mRNAs were sequenced in paired-end, 72-bp read length, with three samples multiplexed per lane, using an Illumina Genome Analyzer 2 (Montpellier Genomix (MGX), Montpellier, France).

## Transcriptome assembly

*De novo* transcriptomes were assembled using high-quality reads (phred > 29) from all seven sequenced samples and an in-house pipeline created using the Velvet (v1.2.01), Oases (v0.2.04) and CDhit EST (v4.5.4) programs [22]. A consensus reference transcriptome was optimized using various parameters, including k-mer length, insert length and expected coverage, as previously described [22]. To suppress non-molluscan transcripts, a BLAST-based transcriptome subtraction strategy was used. Transcripts were compared against the *B. glabrata* draft genome, and transcripts with identities and coverages of less than 80% were deleted. This led to a final subtraction of 45.1% of the transcripts, most of which corresponded to *S. mansoni*. The final size of the transcriptome was 159,711 transcripts.

## Differential expression analysis of transcripts

Quality reads (Phred score >29) were aligned to the assembled transcriptome using the C++ script, Bowtie2 (v2.0.2) (mapping quality score 255), which was run locally using the Galaxy Project server [23]. The DESeq2 software [24] (v2.12; http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html) was run under the default settings to compare duplicate samples from uninfected snails (naive1 and naive2) versus infected samples to quantify differential gene expression (P value < 0.1). A heatmap was constructed to analyze transcript expression patterns (log2 fold change) using Hierarchical Ascending Clustering (HAC) with Pearson correlation, as applied by the Cluster 3.0 [25] and JavaTreeView software packages. We used quantitative real time PCR (Q-RT-PCR) to ascertain RNAseq results. For that purpose the correlation between the rpkm (reads per kilobase per million mapped reads) of RNAseq and the ct (cycle threshold) of Q-RT-PCR was tested and confirmed (see S1 Fig). Primer sequences used in Q-RT-PCR were available in S3 Table.

## FREP RNA interference

Small interfering RNAs (siRNAs) were used to knock down FREP2 (GenBank accession number gi|16303186), FREP3 (gi|18389116) and FREP4 (gi|16303188). Three siRNA duplexes (Eurogentec) were used in conjunction with Invivofectamine transfecting agent (Invitrogen, CA). Two hundred ng of pooled FREP duplex siRNA or GFP siRNA (used as control) were
injected into the cardiac sinus of *B. glabrata* snails as previously described [26]. One, Two and four days following siRNA injection, the knock-down efficiency was confirmed by Q-RT-PCR analysis of FREP2, -3 and -4 mRNA expression. For phenotypic analysis, 90 BgBRE snails were infected with 10 miracidia of SmBRE. Twenty-one days after primo-infection, these snails were individually injected with or without siFREPs or siGFP. Four days later, each snail was subjected to a secondary challenge with 10 miracidia of SmBRE. Fifteen days post-challenge, snails were fixed and the parasite prevalence was determined by identifying schistosome primary sporocysts (SpI), in snails tissue as previously described [27,28].

## Proteomic screening of primed snail plasma

Hemolymph was collected from the head-foot region of 50 BgBRE snails (5 pools of 10 snails, 5 biological replicates) at each time point of the infection process (Naïve, 15DPPI, 25DPPI, and 15DPC). Hemocytes were pelleted by centrifugation at 2500 rpm for 15 min at 4°C, and discarded. Hemoglobin was separated from the plasma by ultracentrifugation at 55,000 rpm for 2.5 h at 4°C and the supernatant plasma samples were recovered. Protein concentration of the samples were estimated using a 2D Quant kit (GE Healthcare life sciences). The plasma (5 biological replicates for each of the 4 time points) was subjected to 2D gel electrophoresis on 12% SDS-PAGE gels using 100 µg of plasma, denaturing buffer, and 17 cm IPG Strips, pH 3-10 non-linear gradient (BioRad). Gels were stained with mass spectrometry (MS)-compatible silver staining, and comparative analysis of digitized proteome maps was performed using the PDQuest 7.4.0 image analysis software (Bio-Rad). Spots showing obvious qualitative and quantitative (at least 2-fold) differences were excised from the gel and characterized by nanoscale capillary liquid chromatography on an Ultimate 3000 coupled to a LTQ-Orbitrap tandem mass spectrometer (nanoLC-MS/MS) mapped to Swiss prot-trembl and the *Biomphalaria glabrata* Brazil transcriptome (http://ihpe.univ-perp.fr/enseignement/). A protein was considered to be correctly identified if at least two peptides were confidently matched to database sequences with an overall Mascot score greater than 50 [29].

## Plasma transfer

Plasma samples were recovered from 15DPPI or naïve snails as described above. PBS-snail solution (8.41 mM Na₂HPO₄, 1.65 mM NaH₂PO₄, 45.34 mM NaCl) was used as the control injection. Twenty µl of each sample (naïve plasma, 15DPPI plasma, or TBS Tween solution) were injected into 25 BgBRE snails per group. One group of 48 BgBRE snails was left untreated and used as a control for infection. Fifteen days after injections, the snails in the four experimental groups were infected with 10 miracidia/snail of SmBRE. This period of 15 days was used to ensure that an observed phenotype was not a direct effect of the injected molecules, but rather reflected their ability to activate the snail's immune response pathways. Fifteen days after infections, snails were dissected and the parasite prevalence was quantified as described above.

## Statistical analysis

To test for significant differences in prevalence, Fisher's exact test was used, with $P \le 0.05$. For the proteomic approach, statistically significant quantitative differences between spots were tested using a Mann-Whitney U test, as applied with the PDQuest 7.4.0 software (BioRad).

## Ethical statement

The laboratory and experimenters possessed an official certificate of the French Ministry of National Education, Research, and Technology, CNRS and DRAAF Languedoc Roussillon for
<figure><img src="image_6.png" /><figcaption>**Fig 2.** The immune response of *B. glabrata* to *S. mansoni* infection. The Brazilian strain of albino *B. glabrata* (BgBRE) is 100% susceptible (for 10 miracidia and upwards) to its corresponding strain of *S. mansoni* (SmBRE). When a snail is infected with 10 miracidia of *S. mansoni* within the same individual compatible and incompatible interactions occur, 3 to 4 miracidia develop normally in the snail's tissues while the others are recognized and encapsulated by the snail's cellular immune response. A. Six-day-old sporocyst in a compatible interaction. B. Encapsulated sporocyst 48 h after primary infection in an incompatible interaction. C. Six-day-old sporocyst in a primed snail. Primed BgBRE are 100% protected against a secondary challenge with SmBRE. Sporocysts from secondary challenge were neutralized by immune humoral factors.</figcaption></figure>

doi:10.1371/journal.ppat.1005361.g002

experiments on animals, animal housing, and animal breeding (# A66040; decree # 87-848, October 19, 1987; and authorization # 007083).

## Results

### Histological inspection of *B. glabrata* immune response phenotypes

Three types of immuno-biological interactions were observed following infection of *B. glabrata* by *S. mansoni*. After primo-infection, a compatible interaction was characterized by the ability of miracidia to penetrate, transform into sporocysts and develop normally in snail tissues (30 to 40% of the entering miracidia produce SPI) (Fig 2A). In an incompatible interaction, miracidia were immediately recognized, encapsulated and killed by hemocytes; in this case, a multi-cellular capsule could be observed surrounding the parasite (60 to 70% of the entering miracidia were killed) (Fig 2B). Finally, in primed snails subjected to a secondary challenge, the sporocysts degenerated in the snail tissues. Encapsulation was never observed nor was the accumulation of hemocytes observed near the parasite (100% of the entering miracidia were killed by humoral factors) (Fig 2C). These findings indicate that a cellular immune response was the main outcome following a primo-infection, whereas an exclusively humoral immune response was found in primed snails subjected to a secondary challenge (Fig 2).

### A whole-snail transcriptomic approach for investigating the molecular basis of innate immune memory

The RNAseq values were validated by a correlation analysis performed with our quantitative Q-RT-PCR data (S1 Fig, S3 Table). We obtained a good correlation (R² = 0.768), indicating that the representations of transcripts in our RNAseq results appeared to be proportional to those obtained through Q-RT-PCR amplifications.

This RNAseq experiment allowed us to identify 1887 differentially expressed transcripts, most of which were found in the 1DPPI, 25DPPI and DPC samples (Fig 3, S1 Table). At 1DPPI, most of the differentially expressed transcripts were under-represented (Fig 3). This
<figure><img src="image_7.png" /><figcaption>**Fig 3.** RNAseq analysis of the innate immune memory response of *B. glabrata* to *S. mansoni*. Heatmap showing differentially represented transcripts compared to naïve snails, as identified by DESeq2 analysis (p < 0.1). Color scale indicates the Log2FC ratio from under-represented (blue) to over-represented (red) transcripts. Transcripts were grouped into six clusters based on their expression patterns during the process of innate immune memory. Samples were recovered at 1DPPI, 4DPPI, 15DPPI and 25DPPI following primo-infection. Following secondary challenge samples were recovered at 1 day, 4 days and 15 days and pooled into DPC sample. Six clusters are identified: Cluster 1: transcripts over represented more than once all along infection and challenge. Cluster 2: transcripts exclusively over represented in single one condition. Cluster 3: transcripts exclusively over represented after immune challenge (DPC). Cluster 4: transcripts exclusively under represented after immune challenge (DPC). Cluster 5: transcripts exclusively under-represented in single one condition. Cluster 6: transcripts under represented more than once all along infection and challenge. FREP: Fibrinogen-related protein, HSP: Heat-shock protein, PGRP 1-like: Pathogenesis-related protein 1-like, LBP/BPI: lipopolysaccharide-binding protein/bactericidal/permeability-increasing protein, BgLBP/BPI: *Biomphalaria glabrata* LBP/BPI, TEP: thioester-containing protein.</figcaption></figure>

doi:10.1371/journal.ppat.1005361.g003

could reflect the immuno-suppression/modulation induced by *S. mansoni* parasites within the first hours of an infection, when the parasites develop their molecular mimicry strategy [30,31]. At 25DPPI and DPC, the proportions of over- and under-represented transcripts were more comparable (Fig 3). Hierarchical clustering was used to sort the differentially represented transcripts into six clusters, three each corresponding to over-represented transcripts (clusters 1-3) and under-represented transcripts (clusters 4-6) (Fig 3, S1 Table). Cluster 2 included transcripts that were specific to the primo-infection and not activated following the secondary challenge. Clusters 1 and 3 appeared to be the most promising possible sources of innate
immune memory candidates. Cluster 1 comprised transcripts that were over-represented fol-
lowing primo-infection and remained highly expressed throughout infection (sustained
response); it included immune molecules known to be involved in the cellular immune
response (macrophage mannose receptors, thrombospondin), as well as certain pathogen rec-
ognition receptors (PRRs) (selectins and C-type lectins). In contrast, cluster 3 comprised tran-
scripts that were over-represented exclusively following the secondary challenge (Fig 3),
including certain PRRs (FREPs, C-type lectins) and circulating immune effectors (anti-micro-
bial peptides such as mytimacin and LBP/BPI, and biomphalysin). Interestingly, clusters 4, 5,
and 6, which contained the under-represented transcripts, included immune molecules belong-
ing to the same immune functional groups represented in clusters 1 to 3 (Fig 3). This may indi-
cate that there is a trade-off among the immune variants involved in the specific response to S.
mansoni parasites, with some members of the same family (i.e., variants or isoforms) being
down-regulated to permit the over-expression of other variants. In cluster 4, notably, most of
the under-represented transcripts were involved in the cellular and epithelial immune response
(e.g., macrophage mannose receptors, extracellular matrix compounds, mucins). This indicates
a down regulation of epithelial (mucosal immunity) and cellular immunity two major compo-
nents of innate immune response following secondary challenge. The complete list of differen-
tially represented transcripts is presented in S1 Table.

Loss-of-function, FREP-mediated RNA interference

To estimate total FREP expression (Fig 4A) over the entire innate immune memory process based on our RNAseq data, we summed the fold changes for all the differentially represented FREP transcripts (identified in Fig 3). A huge induction (5.096-fold) of FREPs was observed following the secondary challenge (DPC, Fig 4A). SiRNA-mediated knockdown of FREPs were used in combination to ascertain the role played by these FREP candidates in innate immune memory. Knockdown of FREPs was confirmed by Q-RT-PCR; the results were normalized with respect to the mRNA expression levels in naïve snails, and compared to the corresponding levels in siGFP-control-injected snails (Fig 4B). The best knockdown level was reached at 96 h post-siRNA injection, when we observed a decrease of 7.98-, 1.56- and 3.60-fold change for FREP 2, 3 and 4, respectively (Fig 4B). The effect of FREP knockdown on the innate immune memory phenotype was assessed in a typical homologous innate immune memory assay. In untreated snails, the parasite prevalence after secondary challenge was 0% (Fig 4C), indicating that innate immune memory had effectively protected them against a subsequent exposure to S. mansoni. Injection of siGFP did not significantly change the parasite prevalence after sec- ondary challenge (4%; Fig 4C). Following injection of the FREP siRNA, however, 15% of the snails became infected with S. mansoni following secondary challenge (Fisher's exact test P < 0.05) (Fig 4C). Thus, FREP knockdown partially altered the innate immune memory phe- notype, rendering primed snails more susceptible to S. mansoni infection.

Global proteomic comparison of primed snail plasma

To validate our RNAseq results, we focused our attention on the plasmatic compartment of *B. glabrata*, which should contain the actors of the humoral innate immune memory response.
We performed global comparative 2D gel electrophoresis on plasma samples obtained from naive, 15DPPI, 25DPPI, and 15DPC snails. Our quantitative and qualitative bioinformatic analyses identified 29 differentially expressed spots corresponding to 62 different proteins (Fig 5A, S2 Table).

A hit map of the quantitated expression ratios was generated for the entire infection process (Fig 5B, S2 Table). Four clusters were identified (Fig 5B): the proteins in cluster 1 were up-
<figure><img src="image_8.png" /><figcaption>**Fig 4.** FREPs knock-down mediated by RNA interference. A. Cumulative expression [Log2FC (fold change) from the DESeq2 analysis] of FREP transcripts showed that FREPs were over-represented after the secondary challenge (DPC; 5.096 log2 fold change enrichment of FREPs transcripts). Green points corresponded to the differentially represented FREPs transcripts in each samples. Purple bars represent the cumulated Log2FC of FREP transcripts. At 4DPPI no FREP transcripts were differentially expressed, thus no value appeared in the graph. B. siRNA injection against FREP2, FREP3 & FREP4 was carried out and mRNA abundance was monitored during 4 days by Q-RT-PCR. Snails were injected with siRNAs against FREP 2, 3, and 4 or GFP (control), the relevant mRNA levels were assessed following normalization with respect to the S19 gene in siGFP injected snails versus siFREPs injected snails. Knock-down of the three FREPs tested was confirmed at 96h. C. Naïve *B. glabrata* and siFREP-injected snails were subjected to a typical priming experiment: Snails were infected with 10 miracidia of *S. mansoni* as a primo-infection, 21 days later they were injected with siGFP, or SiFREP or not treated and 4 days later they were infected with another 10 miracidia as a secondary challenge. FREP siRNA-injected snails show a significant proportion of non-primed snails (15%; *, binomial test, P < 0.05).</figcaption></figure>

doi:10.1371/journal.ppat.1005361.g004

regulated following the secondary challenge; those in cluster 2 were up-regulated after the primo-infection and remained highly expressed throughout infection (sustained response); those in cluster 3 were up-regulated after the primo-infection and down-regulated following secondary challenge those of cluster 4 were down-regulated. The differentially expressed proteins were also subjected to functional classification based on their putative functions and/or structural features. Five functional groups were identified: innate immune response proteins (receptors, effectors, and regulators), scavengers of reactive oxygen species, cell signaling proteins, gametogenesis-related proteins, and glycolysis-related proteins. Thirty seven percent of the identified molecules belonged to the innate immune response group (23 of 62 proteins) (Fig 5B). Consistent with our RNAseq results, we observed an important enrichment of potential humoral immune candidates in the plasma of primed snails. The molecules involved in immune recognition and opsonization included numerous isoforms belonging to different families of circulating immune receptors that were also identified in our transcriptomic analysis (e.g., C-type lectins, β-1,3-glucan binding protein, and hemagglutinin/amebocyte
A

<figure><img src="image_9.png" /></figure>

B

<figure><img src="image_10.png" /></figure>

C

<figure><img src="image_11.png" /><figcaption>**Fig 5. Role of B. glabrata plasmatic factors in innate immune memory response.** A. 2D gel electrophoresis of plasma proteins. One gel of each plasma sample analysed was shown. Spot numbers of qualitative and quantitative differences were indicated. Four plasma samples were analysed from naïve (uninfected snails), 15DPPI and 25DPPI (recovered at 15 and 25 days after primo-infection) and 15DPC (recovered at 15 days after secondary challenge). B. Heat-Map of the qualitative and quantitative ratio versus naïve sample. Ratios were calculated using PDQuest software between all differentially regulated spots. Blue to red scale indicate ratio values from lower to higher represented spots. Four clusters are identified: Cluster 1: higher-represented proteins exclusively following secondary challenge (15 DPC). Cluster 2: sustained response: higher-represented proteins after the primo-infection and secondary infection. Cluster 3: higher-represented proteins at 15DPPI and thereafter down regulated at 25DPPI and 15DPC. Cluster 4: lower-represented proteins. C. Plasma transfer and effect on prevalence of S. mansoni infection. Four conditions were tested: untreated snails (Control group, n = 48); saline injected snails (control of injection, n = 25); naïve-plasma injected snails (n = 22); and primed-plasma injected snails (n = 25). For all the experimental groups, 15 days following injection, snails were infected with 10 miracidia of SmBRE. * indicated significant differences (P< 0.05).</figcaption></figure>

doi:10.1371/journal.ppat.1005361.g005

aggregation factors). Most of these isoforms were found in clusters 1 and 2, suggesting that they may actively contribute to the innate immune memory response (Fig 5B). Among the immune effector molecules, chitinase (spot 3002, Fig 5B), fatty acid binding protein (spots 4004 and 1001, Fig 5B), and biomphalysin (spot 7315, Fig 5B) were found in clusters 1, 2 and 4, respectively. Biomphalysin was over-represented following primo-infection and secondary challenge in the transcriptomic data, but down regulated at 15DPPI and 15DPC in our proteomic data (cluster 4, spot 7315; Fig 5B). Thus, it seems that biomphalysin is highly transcribed during the innate immune memory response against S. mansoni parasites, but is thereafter consumed at the protein level. Some known regulators of the immune response (e.g., thymosin, spot 3002; serpin, spot 1314; and β chain acetylcholine binding protein, spots 2206 and 3008) were associated with protein clusters 2 and 4 (Fig 5B). In addition, other proteins known to be
involved in host/parasite interactions were identified as being differentially expressed, includ-
ing some anti-oxidant molecules, which were identified in protein clusters 1 and 2 (e.g., gluta-
thione peroxidase 3 precursor, spot 1803; glutathione-S-transferase mu 3-like, spots 7009 and
8006; and thioredoxin peroxidase, spot 4004; Fig 5B).

Proteins involved in gametogenesis were also identified by our global proteomic analysis. They included schistosomin (spot 1004), which is involved in the regulation of gametogenesis [32]; pathogen-related protein (spots 2002 and 1009), which is a seminal fluid protein [33,34]; and ovipostatin (spot 1106), which is involved in oviposition and hatching [35]. These proteins were found in clusters 3 and 4, and were thus down-regulated following primo-infection and/or secondary challenge (Fig 5B) (see S2 Appendix for details).

Plasma transfer

A plasma experiment was carried out to confirm the role played by humoral factors in innate
immune memory. Plasmas were recovered from naïve or 15DPPI snails and injected into naïve
snails. The snails were subjected to infection 15 days post-plasma-injection, and analyzed for
phenotypes at 15 days post-infection (Fig 5C). Injection of saline solution has no significant
effect on the prevalence of parasitic infection compared to the untreated control (88% and
90%, respectively) (Fig 5C). Transfer of naive snail plasma also failed to have any significant
effect (prevalence of 77%) (Fig 5C). However, a statistically significant decrease in the preva-
lence of infection was observed for snails injected with primed snail plasma (from 90% to 68%;
P = 0.04). This protection remained efficient at 15 days post-transfer, confirming that plas-
matic humoral factors function to activate the immune system and prepare the snail to answer
a subsequent encounter with the parasite.

Discussion

Innate immune memory, which is a process through which an organism acquires a (more or less) specific and long-lasting protection against later challenges that persists even after the pathogen is neutralized, has been described in diverse invertebrate phyla [2,3,5,6,14,36-39]. This suggests that innate immune memory exists in a broad range of invertebrates.

A form of innate immune memory was previously demonstrated in *B. glabrata* in response to *S. mansoni* infection, and shown to be dependent on the genotype of the pathogen [4]. Here, we provide evidence that innate immune memory in *B. glabrata* is not associated with an enhancement of the cellular immune response, as has been suggested in other invertebrate species [5,12,13]. Instead, our histological analysis (Fig 2) revealed that in primed snails, the parasites of a secondary challenge fail to develop into sporocysts and are killed by the host, without any observable cellular immune response/hemocytic reaction. Together, these findings show that, after a primo-infection, each successive encounter with a similar parasite initiates an exclusive humoral immune defense response in *B. glabrata*.

Accordingly, we set out to identify the humoral factors through which primed snails are
able to neutralize the sporocysts of a secondary challenge. A whole-snail RNA sequencing
approach was conducted to identify molecular candidates that might be involved in the innate
immune memory response. For RNAseq approach the cluster 1 (sustained-response tran-
scripts) and cluster 3 (secondary challenge-specific transcripts) contained the most promising
candidates (Fig 3). Numerous circulating or cell-surface PRRs were found. These immune
receptors, which recognize terminal sugar residues on the glycans that are attached to the sur-
face proteins of some microorganisms, act in the pathogen recognition and clearance processes
of innate immunity. The identified PRRs belong to various families, including the macrophage
mannose receptors, selectins, C-type lectins and fibrinogen-related proteins (FREPs).
Cytotoxic and cytolytic molecules were also identified, as were anti-oxidant molecules. The latter finding suggests that the reactive oxygen species pathway may be activated as part of the humoral innate immune memory response of *B. glabrata*. Moreover, RNAseq cluster 4 (molecules down-regulated upon secondary challenge) included transcripts whose products are involved in the epithelial and cellular immune response, such as macrophage mannose receptors, extracellular matrix compounds and mucins. This strengthens our contention that there is a shift from cellular or mucosal immunity to humoral immune response during the innate immune memory response in *B. glabrata* snails.

An unexpected result of this study was the identification of RNAseq cluster 3 and proteomic cluster 1. Until now, innate immune memory was believed to be the result of two potential regulatory processes: (i) a sustained immune response consisting of the long-lasting up-regulation of immune molecules after a primo-infection (as observed in RNAseq cluster 1 and proteomic cluster 2); or (ii) a recalled response that consists of the ability to store information and recall it later for a faster and more powerful response against subsequent pathogenic exposure (never been observed in any of the RNAseq or proteomic clusters) [36,39]. Here, we provide the first report of a novel unexpected regulatory process for invertebrate innate immune memory, in which transcripts are specifically up regulated following the secondary challenge (RNAseq cluster 3, proteomic cluster 1) without having been previously activated by the primo-infection. This new observation warrants further study, such as through the identification of its potential molecular hallmarks, including the genetic (activation of transcription factors) and/or epigenetic (DNA methylation, chromatin markers, lncRNA, miRNA, etc.) factors that form the basis for this exclusive secondary challenge innate immune memory response.

To validate the potential roles played by the identified molecular candidates in the observed innate immune memory phenotype, we performed three sets of experiments: we used RNAi to knockdown certain FREPs and examined the priming response; we characterized the snail plasma proteome; and we examined the effect of injecting naïve snails with plasma from primed snails.

As FREP lectins showed extensive differential representation following the secondary challenge (see Fig 4A), and they are known to play key roles in the immunobiological interactions between *B. glabrata* and *S. mansoni*, they were chosen as candidates for RNAi invalidation. The FREPs comprise one or two amino-terminal immunoglobulin domains (IgSF) and a carboxyl-terminal fibrinogen domain (FBG). They belong to a multigene family of at least 14 members, and undergo somatic rearrangements leading to remarkable diversification [40,41]. FREPs can precipitate soluble antigens derived from trematodes. Their expression levels increase in response to infection with *S. mansoni* [42,43], and they form immune complexes with similarly highly polymorphic and individually variable mucins (the SmPoMucs) that act as antigens of *S. mansoni* [44]. The high level of diversification within FREP family members results in a huge proportion of partial FREP transcript sequences generated by the RNAseq *de novo* assembly [22]. We were unable to precisely identify the FREP variants/isoforms present in the various RNAseq clusters. Thus, to confirm the role of FREPs in innate immune memory, we invalidated FREP 2, 3, and 4, as they had previously been demonstrated to be involved in the response of *B. glabrata* to trematodes. More specifically, FREP 2 is involved in immune complexes [44]; FREP 3 knock-down reverts the snail resistance status to trematodes [45]; and microarray analysis showed that FREP 4 is over-expressed following infection by *S. mansoni* [43]. SiRNA-mediated knock-down of FREP 2, 3 and 4 (Fig 4B) was found to reduce the innate immune memory phenotype by 15% (Fig 4C). These results indicate that these FREPs are involved in *B. glabrata* innate immune memory, but demonstrate also that additional molecular partners (and/or other FREP variants) also play roles in the pathogen recognition and innate immune memory of this snail.
Our characterization of the plasmatic proteome confirmed that most of the molecules identified in the RNAseq data were present and differentially expressed at the protein level. The presence of these molecules in snail plasma also confirms that they act as circulating humoral factors. The most promising candidates for participation in innate immune memory were found in protein cluster 1 (secondary challenge-specific proteins) and protein cluster 2 (sustained-response proteins) (Fig 5). Numerous PRRs were identified (e.g., macrophage mannose receptors, C-type lectin, hemagglutinin, and β-1, 3-glucan binding proteins). Interestingly, we did not recover any FREP in this proteomic analysis of primed snail plasma. We speculate that the FREPs were lost from our analysis due to their precipitation with parasitic antigens [44]. Alternatively, and as suggested previously for other families of highly variable molecules [46], the different variants may be expressed at such low levels that they are not visible on 2D gel electrophoresis. Indeed, we recently demonstrated that FREPs are expressed at a low level compared to other PRRs in naive *B. glabrata* snails [22]. Plasmatic proteome was also composed of cytotoxic and cytolytic molecules, including antimicrobial peptides [AMPs; e.g., hydramicin, mytimacin and lipopolysaccharide binding protein/bactericidal permeability-increasing protein (LBP/BPI)], and the biomphalysin. We recently reported the molecular cloning and functional characterization of the *B. glabrata* biomphalysin and its involvement in the killing of S. mansoni [47]. Biomphalysin, which shares a common architecture with proteins belonging to the aerolysin superfamily, is strictly expressed in immune-competent cells. Recombinant biomphalysin was shown to bind to parasitic membranes and exhibit cytotoxic activity toward S. mansoni sporocysts. Our RNAseq and proteomic analyses primarily recovered immune recognition factors, many of which had multiple isoforms. This suggests the involvement of both translational and post-translational regulation. Additional work is now needed to clarify the functions of these molecules in the humoral innate immune memory response of *B. glabrata*. For example, C-type lectins are known to be involved in the innate immune memory responses of mollusks; their expression levels increased in the scallop, *Chlamys farreri*, following vaccination with heat-killed *Vibrio anguillarum*, and successive challenges with *V. anguillarum* or *Micrococcus luteus* enhanced this protection [11]. In addition to their roles in pathogen recognition and opsonization, lectins also possess direct cytotoxic activities. Examples of cytotoxic lectins include plant-derived ricin from *Ricinus communis* beans [48], the fungus-derived N-acetyl-D-galactosamine-specific lectin from *Schizophyllum commune* [49], and the invertebrate-derived hemolytic lectin, CELIII, from *Cucumaria echinata* [50]. Thus, it seems reasonable to hypothesize that some of the lectins identified in the present study may participate in the humoral innate immune memory response both as recognition receptors and as cytotoxic/cytolytic molecules involved in killing the S. mansoni SpIs of the secondary challenge.

Lastly, we explored whether the transfer of plasma from primo-infected snails could provide recipient naïve snails with enhanced immunity. In a previous study, plasma transfer from S. mansoni-resistant *Biomphalaria tenagophila* to susceptible snails was associated with a transfer of resistance [51]. Here, we found that transfer of primed snail plasma to naive snails significantly reduced the prevalence of S. mansoni infection by more than 20% compared to controls (Fig 5). A comparable experiment was previously performed using the mosquito, *Anopheles gambiae*, and its bacterial pathogens [5]. The transfer of cell-free hemolymph from infected mosquitoes into healthy mosquitoes triggered increases in the hemocyte populations of transf erred mosquitoes, indicating that humoral factors could promote a cellular immune response that protected mosquitoes against subsequent bacterial challenges [5]. In *Biomphalaria*, plasma transfer experiments demonstrated that soluble humoral factors are released into the hemo lymph of primed snails, and that the transfer of such factors activates/regulates the humoral immune response in recipients and confers enhanced antischistosomal immunity against sub sequent encounters with the pathogen.
To conclude, we evidence for the first time the molecular basis of innate immune memory in an invertebrate model, and we demonstrate the role of humoral factors in such phenomenon. The previous studies addressing the molecular mechanisms of innate immune memory in invertebrates have mainly suggested the involvement of elevated hemocyte phagocytosis [9,10]. In *B. glabrata*, innate immune memory seems to be supported by humoral factors that trigger the degeneration/death of the parasite. Moreover the innate immune memory protection was previously shown to decrease with increasing neutral genetic distance between the parasite used for the primo-infection and the one used for the secondary challenge [4]. We hypothesized that better protection against a homologous (vs. heterologous) secondary infection involves the mobilization of specific repertoires of *B. glabrata* immune receptors to target certain subsets of *S. mansoni* genotypes. Our present findings prompt us to speculate that the genotype-dependent innate immune memory of *B. glabrata* may be supported by a diverse repertoire of FREPs and other PRRs (Figs 4 and 5). As previously suggested [22], PRRs might serve as collaborative recognition factors that can be processed as homologous or heterologous multimers, which then act as immune recognition complexes to increase the host's PRR repertoire and mediating anti-pathogen responses. Once recognized, the pathogen is neutralized by the release of cytotoxic/cytolytic circulating factors. Here, we demonstrate that biomphalysin is highly transcribed following secondary challenge and consumed at the protein level during innate immune memory response, suggesting that biomphalysin may play a major role in neutralizing the pathogen following immune recognition. Clearly, future studies are needed to fully understand how these molecules mediate and regulate the specificity of the innate immune memory defense seen in *B. glabrata*.

To reconcile mechanisms with phenomena a molecular characterisation of innate immune memory in *Biomphalaria/Schistosoma* model was undertook. Transcriptomic, proteomic, and functional validation lead us to a new molecular comprehension of innate immune memory processes and prompted us to revisit the artificial dichotomy between innate and memory immunity in invertebrate systems.

**Supporting Information**

**S1 Fig. Correlation between RNAseq and Q-RT-PCR.** The correlations between the RNAseq reads per kilobase per million (log2 RPKM) values and the corresponding Q-RT-PCR cycle threshold (Ct) values were calculated and shown for 36 randomly selected transcripts. RPKM values were calculated using a local alignment with Bowtie2 on *de novo* assembled *B. glabrata* transcriptome. Cts were calculated by Q-RT-PCR using the Naive1/Naive2 biological samples previously used for RNAseq. Log2 RPKM values were an average of two biological replicates of two independent snail pools and represented in log2 scale. Horizontal lines indicate the standard error between these two points. Non-normalized Ct values were an average of three technical replicates of Q-RT-PCR. The vertical lines indicate the standard error between these three points. A dependency between variables, increasing of RPKM and decreasing of Ct is shown by least-square regression analysis and its correlation coefficient R² = 0,76804. (TIF)

**S1 Table. List of differentially represented transcripts in RNAseq clusters.** Quality reads (Phred score >29) were aligned on the transcriptome assembly using the C++ script Bowtie2 (v2.0.2) (255 score) running thanks local engine using Galaxy Project server (Giardine, Riemer et al. 2005). The DESeq2 software (Love, Huber et al. 2014) (v2.12; http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html) (defaults settings) allows for quantifying the differential gene expression with comparing two biological duplicates from uninfected snails sample (Bre1 and Bre2) against infected samples (Pvalue<0.1). For each cluster transcript ID,
Blast2GO annotation and Log2FC results were indicated.
(XLSX)

**S2 Table. All cluster identified in Fig 5.** Ratio of expression and mapping annotation were indicated for all the sequences identified as differentially regulated in the plasma comparative proteomic approach.
(XLSX)

**S3 Table. Primers used for Q-RT-PCR.** Identification numbers (iD), Blast2GO annotation and sequences of forward and reverse primers were indicated for 36 randomly selected transcripts of Naive1/Naive2 RNAseq biological samples. These sequences were used to analyse the correlation between the RNAseq reads per kilobase per million (log2 RPKM) values and the corresponding Q-RT-PCR cycle threshold (Ct) see S1 Fig.
(XLSX)

**S1 Appendix. Expanded Materials and Methods section.**
(DOCX)

**S2 Appendix. Trade-off between reproduction and immunity.**
(DOCX)

## Acknowledgments

We thank Raphaël Palisse for participating in this project as a Master student at Perpignan University. We thank Nathalie Arancibia and Bernard Dejean for their technical assistance. This paper is dedicated to our colleague, Bernard Dejean.

## Author Contributions

Conceived and designed the experiments: SP JP DD BG. Performed the experiments: SP JP DD FCN MAO JFA RG NMD SKJ BG. Analyzed the data: SP JP DD RG GM AT BG. Wrote the paper: SP JP DD RG NMD GM AT BG.

## References

1. Loker ES, Adema CM, Zhang SM, Kepler TB (2004) Invertebrate immune systems—not homogeneous, not simple, not well understood. *Immunol Rev* 198: 10–24. PMID: 15199951

2. Kurtz J (2004) Memory in the innate and adaptive immune systems. *Microbes Infect* 6: 1410–1417. PMID: 15596128

3. Kurtz J, Franz K (2003) Innate defence: evidence for memory in invertebrate immunity. *Nature* 425: 37–38. PMID: 12955131

4. Portela J, Duval D, Rognon A, Galinier R, Boissier J, et al. (2013) Evidence for specific genotype-dependent immune priming in the lophotrochozoan Biomphalaria glabrata snail. *J Innate Immun* 5: 261–276. doi: 10.1159/000345909 PMID: 23343530

5. Rodrigues J, Brayner FA, Alves LC, Dixit R, Barillas-Mury C (2010) Hemocyte differentiation mediates innate immune memory in Anopheles gambiae mosquitoes. *Science* 329: 1353–1355. doi: 10.1126/science.1190689 PMID: 20829487

6. Sadd BM, Schmid-Hempel P (2006) Insect immunity shows specificity in protection upon secondary pathogen exposure. *Curr Biol* 16: 1206–1210. PMID: 16782011

7. Lavine MD, Strand MR (2002) Insect hemocytes and their role in immunity. *Insect Biochem Mol Biol* 32: 1295–1309. PMID: 12225920

8. Zanker KS (2010) Immunology of Invertebrates: Humoral. *Encyclopedia of Life Sciences* John Wiley & Sons, Ltd: Chichester.

9. Netea MG, Quintin J, van der Meer JW (2011) Trained immunity: a memory for innate host defense. *Cell Host Microbe* 9: 355–361. doi: 10.1016/j.chom.2011.04.006 PMID: 21575907
10. Rowley AF, Powell A (2007) Invertebrate immune systems specific, quasi-specific, or nonspecific? *J Immunol* 179: 7209-7214. PMID: 18025161

11. Wang J, Wang L, Yang C, Jiang Q, Zhang H, et al. (2013) The response of mRNA expression upon secondary challenge with Vibrio anguillarum suggests the involvement of C-lectins in the immune priming of scallop Chlamys farreri. *Dev Comp Immunol* 40: 142-147. doi: 10.1016/j.dci.2013.02.003 PMID: 23434464

12. Roth O, Kurtz J (2009) Phagocytosis mediates specificity in the immune defence of an invertebrate, the woodlouse Porcellio scaber (Crustacea: Isopoda). *Dev Comp Immunol* 33: 1151-1155. doi: 10.1016/j.dci.2009.04.005 PMID: 19416736

13. Pope EC, Powell A, Roberts EC, Shields RJ, Wardle R, et al. (2011) Enhanced cellular immunity in shrimp (Litopenaeus vannamei) after 'vaccination'. *PLoS One* 6: e20960. doi: 10.1371/journal.pone.0020960 PMID: 21698190

14. Hauton C, Smith VJ (2007) Adaptive immunity in invertebrates: a straw house without a mechanistic foundation. *Bioessays* 29: 1138-1146. PMID: 17935208

15. Little TJ, Colegrave N, Sadd BM, Schmid-Hempel P (2008) Studying immunity at the whole organism level. *Bioessays* 30: 404-405; author reply 406. doi: 10.1002/bies.20738 PMID: 18348252

16. Sire C, Rognon A, Théron A (1998) Failure of *Schistosoma mansoni* to reinfect *Biomphalaria glabrata* snails: acquired humoral resistance or intra-specific larval antagonism? *Parasitology* 117: 117-122. PMID: 9778633

17. Adema CM, Bayne CJ, Bridger JM, Knight M, Loker ES, et al. (2012) Will all scientists working on snails and the diseases they transmit please stand up? *PLoS Negl Trop Dis* 6: e1835. doi: 10.1371/journal.pntd.0001835 PMID: 23301104

18. Sanchez Alvarado A (2013) On the trail of a tropical disease. *Elife* 2: e01115. doi: 10.7554/eLife.01115 PMID: 23908770

19. Theron A, Pages JR, Rognon A (1997) Schistosoma mansoni: Distribution patterns of miracidia among Biomphalaria glabrata snail as related to host susceptibility and sporocyst regulatory processes. *Experimental Parasitology* 85: 1-9. PMID: 9024196

20. Roger E, Mitta G, Mone Y, Bouchut A, Rognon A, et al. (2008) Molecular determinants of compatibility polymorphism in the Biomphalaria glabrata/Schistosoma mansoni model: New candidates identified by a global comparative proteomics approach. *Molecular and Biochemical Parasitology* 157: 205-216. PMID: 18083248

21. Gourbal BEF, Gabrion C (2006) Histomorphological study of the preputial and clitoral glands in BALB/c mice with experimental Taenia crassiceps infections. *Journal of Parasitology* 92: 189-192. PMID: 16629336

22. Dheilly NM, Duval D, Mouahid G, Emans R, Allienne JF, et al. (2015) A family of variable immunoglobulin and lectin domain containing molecules in the snail Biomphalaria glabrata. *Dev Comp Immunol* 48: 234-243. doi: 10.1016/j.dci.2014.10.009 PMID: 25451302

23. Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, et al. (2005) Galaxy: a platform for interactive large-scale genome analysis. *Genome Res* 15: 1451-1455. PMID: 16169926

24. Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol* 15: 550. PMID: 25516281

25. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. *Proc Natl Acad Sci U S A* 95: 14863-14868. PMID: 9843981

26. Baeza Garcia A, Pierce RJ, Gourbal B, Werkmeister E, Colinet D, et al. (2010) Involvement of the cytokine MIF in the snail host immune response to the parasite Schistosoma mansoni. *PLoS Pathog* 6.

27. Allienne JF, Theron A, Gourbal B (2011) Recovery of primary sporocysts in vivo in the Schistosoma mansoni/Biomphalaria glabrata model using a simple fixation method suitable for extraction of genomic DNA and RNA. *Exp Parasitol* 129: 11-16. doi: 10.1016/j.exppara.2011.06.003 PMID: 21726555

28. Moné Y, Mitta G, Duval D, Gourbal BE (2010) Effect of amphotericin B on the infection success of *Schistosoma mansoni* in *Biomphalaria glabrata*. *Exp Parasitol* 125: 70-75. doi: 10.1016/j.exppara.2009.12.024 PMID: 20067790

29. Beltran S, Gourbal B, Boissier J, Duval D, Kieffer-Jaquinod S, et al. (2010) Vertebrate host protective immunity drives genetic diversity and antigenic polymorphism in Schistosoma mansoni. *Journal of Evolutionary Biology* 24: 554-572. doi: 10.1111/j.1420-9101.2010.02190.x PMID: 21159003

30. Damian RT (1997) Parasite immune evasion and exploitation: reflections and projections. *Parasitology* 115 Suppl: S169-175. PMID: 9571701

31. Guillou F, Roger E, Mone Y, Rognon A, Grunau C, et al. (2007) Excretory-secretory proteome of larval Schistosoma mansoni and Echinostoma caproni, two parasites of Biomphalaria glabrata. *Mol Biochem Parasitol* 155: 45-56. PMID: 17606306
32. Zhang SM, Nian H, Wang B, Loker ES, Adema CM (2009) Schistosomin from the snail Biomphalaria glabrata: expression studies suggest no involvement in trematode-mediated castration. Mol Biochem Parasitol 165: 79-86. doi: 10.1016/j.molbiopara.2009.01.006 PMID: 19393164

33. Davies SJ, Chapman T (2006) Identification of genes expressed in the accessory glands of male Mediterranean Fruit Flies (Ceratitis capitata). Insect Biochem Mol Biol 36: 846-856. PMID: 17046598

34. Ram KR, Wolfner MF (2009) A network of interactions among seminal proteins underlies the long-term postmating response in Drosophila. Proc Natl Acad Sci U S A 106: 15384-15389. doi: 10.1073/pnas.0902923106 PMID: 19706411

35. Koene JM, Sloot W, Montagne-Wajer K, Cummins SF, Degnan BM, et al. (2010) Male accessory gland protein reduces egg laying in a simultaneous hermaphrodite. PLoS One 5: e10117. doi: 10.1371/journal.pone.0010117 PMID: 20404934

36. Kurtz J (2005) Specific memory within innate immune systems. Trends Immunol 26: 186-192. PMID: 15797508

37. Vazquez L, Alpuche J, Maldonado G, Agundis C, Pereyra-Morales A, et al. (2009) Review: Immunity mechanisms in crustaceans. Innate Immun 15: 179-188. doi: 10.1177/1753425909102876 PMID: 19474211

38. Cooper EL, Rinkevich B, Uhlenbruck G, Valembois P (1992) Invertebrate immunity: another viewpoint. Scand J Immunol 35: 247-266. PMID: 1535984

39. Little TJ, Kraaijeveld AR (2004) Ecological and evolutionary implications of immunological priming in invertebrates. Trends Ecol Evol 19: 58-60. PMID: 16701227

40. Zhang SM, Adema CM, Kepler TB, Loker ES (2004) Diversification of Ig superfamily genes in an invertebrate. Science 305: 251-254. PMID: 15247481

41. Zhang SM, Loker ES (2003) The FREP gene family in the snail Biomphalaria glabrata: additional members, and evidence consistent with alternative splicing and FREP retrosequences. Fibrinogen-related proteins. Developmental and Comparative Immunology 27: 175-187. PMID: 12590969

42. Adema CM, Hertel LA, Miller RD, Loker ES (1997) 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 94: 8691-8696. PMID: 9238039

43. Hanington PC, Lun CM, Adema CM, Loker ES (2010) 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 40: 819-831. doi: 10.1016/j.ijpara.2009.12.005 PMID: 20083115

44. Moné Y, Gourbal B, Duval D, Du Pasquier L, Kieffer-Jaquinod S, et al. (2010) 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 4: e813. doi: 10.1371/journal.pntd.0000813 PMID: 20838648

45. Hanington PC, Forsy MA, Dragoo JW, Zhang SM, Adema CM, et al. (2010) Role for a somatically diversified lectin in resistance of an invertebrate to parasite infection. Proc Natl Acad Sci U S A 107: 21087-21092. doi: 10.1073/pnas.1011242107 PMID: 21084634

46. Dheilly NM, Nair SV, Smith LC, Raftos DA (2009) Highly variable immune-response proteins (185/333) from the sea urchin, Strongylocentrotus purpuratus: proteomic analysis identifies diversity within and between individuals. J Immunol 182: 2203-2212. doi: 10.4049/jimmunol.07012766 PMID: 19201874

47. Galinier R, Portela J, Mone Y, Allienne JF, Henri H, et al. (2013) Biomphalysin, a New beta Pore-forming Toxin Involved in Biomphalaria glabrata Immune Defense against Schistosoma mansoni. PLoS Pathog 9: e1003216. doi: 10.1371/journal.ppat.1003216 PMID: 23555242

48. Endo Y, Mitsui K, Motizuki M, Tsurugi K (1987) The mechanism of action of ricin and related toxic lectins on eukaryotic ribosomes. The site and the characteristics of the modification in 28 S ribosomal RNA caused by the toxins. J Biol Chem 262: 5908-5912. PMID: 3571242

49. Chumkhunthod P, Rodtong S, Lambert SJ, Fordham-Skelton AP, Rizkallah PJ, et al. (2006) Purification and characterization of an N-acetyl-D-galactosamine-specific lectin from the edible mushroom Schizophyllum commune. Biochim Biophys Acta 1760: 326-332. PMID: 16507335

50. Oda T, Tsuru M, Hatakeyama T, Nagatomo H, Muramatsu T, et al. (1997) Temperature- and pH-dependent cytotoxic effect of the hemolytic lectin CEL-III from the marine invertebrate Cucumaria echinata on various cell lines. J Biochem 121: 560-567. PMID: 9133626

51. Pereira CA, Martins-Souza RL, Correa A Jr., Coelho PM, Negrao-Correa D (2008) Participation of cell-free haemolymph of Biomphalaria tenagophila in the defence mechanism against Schistosoma mansoni sporocysts. Parasite Immunol 30: 610-619. doi: 10.1111/j.1365-3024.2008.01062.x PMID: 19067842