Skip to main content
  • Research article
  • Open access
  • Published:

Clam focal and systemic immune responses to QPX infection revealed by RNA-seq technology

Abstract

Background

The hard clam Mercenaria mercenaria is an important seafood species widely exploited along the eastern coasts of the United States and play a crucial role in coastal ecology and economy. Severe hard clam mortalities have been associated with the protistan parasite QPX (Quahog Parasite Unknown). QPX infection establishes in pallial organs with the lesions typically characterized as nodules, which represent inflammatory masses formed by hemocyte infiltration and encapsulation of parasites. QPX infection is known to induce host changes on both the whole-organism level and at specific lesion areas, which imply systemic and focal defense responses, respectively. However, little is known about the molecular mechanisms underlying these alterations.

Results

RNA-seq was performed using Illumina Hiseq 2000 (641 Million 100 bp reads) to characterize M. mercenaria focal and systemic immune responses to QPX. Transcripts were assembled and the expression levels were compared between nodule and healthy tissues from infected clams, and between these and tissues from healthy clams. De novo assembly reconstructed a consensus transcriptome of 62,980 sequences that was functionally-annotated. A total of 3,131 transcripts were identified as differentially expressed in different tissues. Results allowed the identification of host immune factors implicated in the systemic and focal responses against QPX and unraveled the pathways involved in parasite neutralization. Among transcripts significantly modulated upon host-pathogen interactions, those involved in non-self recognition, signal transduction and defense response were over-represented. Alterations in pathways regulating hemocyte focal adhesion, migration and apoptosis were also demonstrated.

Conclusions

Our study is the first attempt to thoroughly characterize M. mercenaria transcriptome and identify molecular features associated with QPX infection. It is also one of the first studies contrasting focal and systemic responses to infections in invertebrates using high-throughput sequencing. Results identified the molecular signatures of clam systemic and focal defense responses, to collectively mediate immune processes such as hemocyte recruitment and local inflammation. These investigations improve our understanding of bivalve immunity and provide molecular targets for probing the biological bases of clam resistance towards QPX.

Background

The hard clam, Mercenaria mercenaria, is an ecologically- and economically-important marine bivalve species that thrives along the northeastern coasts of the United States and Maritime Canada. In the past few decades, the hard clam industry has been severely impacted by a protistan parasite called QPX (Quahog Parasite Unknown), which is responsible for mortality episodes in both wild and cultured clam populations [17]. QPX is believed to be an opportunistic pathogen and has been detected in a wide variety of environmental substrates including sediments, marine aggregates and other organic matrices [810]. Interestingly, previous reports highlight the ability of QPX to sustain very low abundance in clams without causing disease outbreaks until it encounters hosts with reduced immunity or following shifts of environmental conditions that add to the virulence of the parasite, under which conditions QPX can take advantage to establish infection sometimes leading to severe clam mortality events [4, 11].

Lesions caused by QPX, usually associated with the presence of nodules, are commonly found in pallial tissues, such as alongside the inner rim of the mantle or at the base of the siphon [1, 2]. These places are widely considered as the portal of entry for QPX cells acquired from the surrounding environment during suspension-feeding [5, 7]. The QPX nodules represent inflammatory masses containing both parasite cells and abundant clam hemocytes, resulting from a series of comprehensive host immune responses induced by the infection that leads to massive focal hemocyte infiltration, parasite encapsulation and partial necrosis of the affected area [1]. Like other invertebrates, the hard clam lacks the specific immune responses and their defense mechanism mainly relies on the effectors of innate immunity, which is mediated by circulating hemocytes and highly diversified humoral antimicrobial factors. These cellular and humoral immune components work in a synergistic way to initiate the recognition, segregation and ultimately elimination of pathogens and other non-self entities [12, 13]. The launching of innate immune responses involves myriad cellular and humoral events modulated not only at the infection sites (focally) but also at a larger, whole-organism scale (systemically). In general, the focal response represents the alterations driven by direct host-pathogen interactions at the infection sites where direct cell-cell (e.g., molecular patterns) interactions mediate the response; while the systemic response reflects overall modifications within the host as a result of the ongoing infection and is mainly associated with dynamic changes of circulating hemocytes and their secreted immune mediators.

Most of the previous investigations have solely focused on the systemic response of M. mercenaria against QPX during the infection events, where changes in cellular and humoral immune parameters (e.g., anti-QPX activity and lysozyme activity in clam plasma, hemocyte phagocytic activity, reactive oxygen species (ROS) production, etc.) as well as expression of a limited number of immune-related genes in tissues and circulating hemocytes were assessed [11, 1416]. In contrast, no previous studies have focused on the characterization of clam focal response at the infection sites. Given the fact that QPX disease is usually focal with formation of well-delimited lesions, the study of clam immune responses at the infection site in the lesions per se is of specific value as it provides insights to better characterize cellular interactions between the hard clam and QPX upon their contact. In this framework, QPX disease in clams offer a unique opportunity to contrast focal and systemic responses against microbial diseases in invertebrates allowing for a more comprehensive understanding of defense strategies used by these animals to fend microbial attacks.

Our study aimed to characterize the gene regulation features of M. mercenaria during QPX infection by profiling the transcripts at the infection lesion and compare focal clam responses with systemic responses detected in healthy tissues from infected clams (in addition to a parallel comparison with tissues from healthy clams). This study allowed the identification of factors involved in the interactions with the parasite as well as molecular pathways activated by the host to neutralize QPX.

Results and discussion

Illumina sequencing and de novo transcriptome assembly

The main objective of this study was to identify molecular features of M. mercenaria in response to QPX infection and to compare the immune-related pathways involved in the lesion-specific focal response with the whole-organism scale systemic response. The high-throughput Illumina RNA sequencing and de novo assembly employed in this investigation allowed the construction of the transcriptome in the absence of M. mercenaria genome information. A total of 640,596,320 of 100 bp raw reads were generated from the Illumina paired-end sequencing with about 27 to 48 Millions paired-end reads generated from each of the 9 sequenced libraries (Table 1, Fig. 1a). The short read sequences generated from this RNAseq project have been deposited at the NCBI short Read Archive database under the SRA accession number SRP068241. Trimming and filtering procedures yielded 606,021,407 clean reads that were used for the de novo assembly of consensus transcriptome based on all sequenced RNA libraries in order to maximize the diversity of transcripts. This allowed 90.61 to 92.20 % of the reads from the 9 libraries be used for the transcriptome assembly. A total of 62,980 transcripts ranging from 201 to 23,103 bp with average size of 1297.59 bp and median size of 835 bp were produced from the assembly after low FPKM and rare isoforms filtering. The size distribution of all the de novo assembled transcripts is shown in Fig. 1b. Once the transcriptome was constructed, the 9 libraries were individually remapped to the 62,980 transcripts and resulted with 85.27 to 89.05 % of reads remapping. Theses counting data were then used for DE analysis.

Table 1 RNA samples for RNA-seq libraries. Each pool is made with equal amounts of RNA from 3 individual clams. Pools A and B were derived from the same clams
Fig. 1
figure 1

M. mercenaria de novo assembled transcriptome summary. a Transcriptome sequencing, assembly and annotation overview. b Assembled transcripts size distribution. c Distribution of the top 10 species with most homologues to M. mercenaria. Transcripts were searched using BLASTx against NCBI nr database with a cutoff value of E < 10E-03

Transcriptome functional annotation

The transcriptome annotation performed using Blast2GO returned a total of 19,107 transcripts (30.3 %) with significant BlastX homology matches to other sequences in NCBI nr database (E-value < 10E-03) (Fig. 1). Not surprisingly, the top 3 species that had the most similarity to M. mercenaria sequences were mollusks with available genomes and included the Pacific oyster, Crassostrea gigas (7,467), followed by the owl limpet Lottia gigantea (3,539) and the California sea slug Aplysia californica (1,931) (Fig. 1c). KEGG Orthology (KO) terms were assigned to 6,425 sequences and reference pathways were mapped to the KEGG database based on the assigned KO terms (Fig. 1a, Additional file 1). A total of 29,815 sequences were identified to match to at least one conserved protein domain in the InterPro database (Fig. 1a, Additional file 1).

Gene ontology (GO) assignments were used to classify functions of the predicted clam proteins. Based on sequence similarity (E-value of 10E-03), 13,584 sequences were assigned to at least one GO annotation (Fig. 1a, Additional file 1). As summarized in Fig. 2, a total of 8,168, 4,600 and 4,231 sequences were respectively categorized into the three main categories: biological process, cellular component, and molecular function at the second functional annotation level. The most dominant terms presented in the three categories are the “cellular process”, “metabolic process”, “binding”, “catalytic activity”, “cell”, and “organelle”. Very few transcripts were clustered into “rhythmic process”, “cell killing”, “protein tag”, “channel regulator activity”, “nucleoid” or “virion”. It is noticeable that a good fraction of transcripts were clustered into the immune-related categories of response to stimulus (503), immune system process (43) and biological adhesion (38). Those transcripts were of special interest given that they might be involved in the M. mercenaria defense and resistance toward QPX infection.

Fig. 2
figure 2

Gene Ontology (GO) annotations of the M. mercenaria transcriptome. GO terms were identified by Blast2GO and the results were summarized in three main GO categories: biological process (8,168 annotations), cellular component (4,600 annotations), molecular function (4,231 annotations) at level-2

A significant portion (69.7 %) of M. mercenaria transcripts did not match any BlastX hit in NCBI nr database, in agreement with previous transcriptomic studies in mollusks [1720]. Most of the unannotated transcripts may represent transcripts spanning untranslated mRNA regions, or transcripts containing only non-conserved protein domains [21, 22].

Identification of differentially expressed transcripts

The generated transcriptome was used as a reference for downstream investigations of global gene expression in the three different tissues of interest (nodule, non-nodule and healthy) to identify genes associated with M. mercenaria’s focal and systemic immune response against QPX. A gene-isoform relationship was estimated using RSEM over Trinity output isoforms. Results showed that about 43 % (27,021) of all the transcripts had 1 isoform, 19 % (12,307) had 2 isoforms and 38 % (23,652) had 3 isoforms, suggesting extensive isoform diversity in M. mercenaria transcriptome. Transcript isoform variation could affect mRNA stability, localization and translation, as well as the production of protein variants that differ in localization or function [23].

By comparing the number of transcripts expressed in each sample, the contribution of specific samples to the analysis can be estimated. The highest number of expressed transcripts was found in the nodules of infected clams, which were closely followed by that found in the healthy clam samples (Fig. 3a). The lowest number of expressed transcripts came from non-nodule samples of QPX infected clams, with about 1,500 less transcripts expressed than the other two samples. Read coverage, which is critical in accurate determination of fold change, averaged 477, 422 and 509 reads per transcript for nodule, non-nodule and healthy tissue samples, respectively (Fig. 3b).

Fig. 3
figure 3

Number of transcripts expressed and reads coverage in each sample. a The X-axis indicates the sample (A1-A3: nodule; B1-B3: non-nodule; C1-C3: healthy, refer to Table 1 for more details). The Y-axis indicates the number of transcripts expressed in the samples. b Summary statistics of the reads coverage in each sample

Statistical analysis by DEseq identified 3,131 differentially expressed (DE) transcripts from the pair-wise comparisons (|log2 (fold change)| >2, adjusted p-value <0.001) between clam tissue samples (Fig. 4). In nodules, a total of 829 transcripts, including 408 over-expressed and 421 under-expressed transcripts were identified as compared to non-nodule samples of QPX infected clams. Compared to tissues from healthy clams, 1,591 DE transcripts were identified in nodules with 864 over-expressed and 727 under-expressed transcripts. Similarly, 1,681 DE transcripts were obtained from the comparison between healthy and non-nodule clam tissues, of which 513 and 1,168 were over- and under-expressed, respectively (Fig. 4, Additional file 3). A total of 1,694 of these DE transcripts had protein homologs found in NCBI nr database by Blastx searches (e-value <10E-03), which were further examined for their putative functions during M. mercenaria immune response toward QPX. Annotated transcripts were subsequently grouped into curated categories according to their biological functions based on the gene ontology (GO) terms and literature searches highlighting immune functions. DE transcript sets were further examined in reference to the assigned KO terms for the analysis of pathways regulation. Significantly enriched KEGG pathways were identified via the Fisher’s exact test (P < 0.01) [24].

Fig. 4
figure 4

Heatmap (a) and number of differentially expressed transcripts (b and c) across all samples (FDR ≤0.001, and |log2 (fold change)| >2). Over-expressed transcripts are shown in red in A and are enumerated in B. Under-expressed transcripts are shown in green in A and are enumerated in C. Replicate biological samples are displayed in A (see Table 1 for more details). For B and C: the cross-shaded areas inside each bar represent the number of transcripts with higher (B) or lower (C) expression levels in a condition as compared to the other two (e.g., the expression levels of 38 focally over-expressed transcripts in B were higher in nodule than in non-nodule tissues and healthy clams)

Here we specifically focus on DE transcripts drawn from comparisons between “nodule vs. non-nodule” and “non-nodule vs. healthy” tissue samples, which were respectively considered to reflect the transcriptomic changes caused by “focal” and “systemic” immune response of M. mercenaria toward QPX infection, respectively. The overview of DE transcripts drawn from these two responses is presented in Fig. 4 and Additional files 2 and 3.

Differential expression of immune related transcripts

Focal response

Clam focal response reflected the alterations caused by direct clam-QPX interaction at the infection site (Additional file 2). QPX nodules are inflammatory masses resulting from massive hemocytes infiltration and encapsulation of parasite cells [1, 2]. This process largely relies on the motility and adhesion properties of hemocytes, thus allowing these cells to migrate throughout the circulatory system and recruit to the infection site. Hemocytes can sense stimuli in host tissues through an array of cell surface receptors, and use these cues to adjust their behavior accordingly [25]. The activation of hemocytes requires the binding of specific ligands to the cell surface receptors, which subsequently initiate the transduction of extracellular signals into the cytoplasm via a variety of signaling pathways, thus inducing a series of hemocyte-mediated immune response such as phagocytosis, encapsulation (prominent response against QPX in clams), ROS production, as well as secretion of immune effectors and cytokines [2628]. A collection of DE transcripts involved in these defense processes was identified during focal response, suggesting that strong and comprehensive host-pathogen interactions were taking places inside the QPX lesions (Table 2, Additional file 2). A large fraction of the DE transcripts of focal response was annotated as receptors or molecules with receptor activities, which putatively contribute to the host defense against QPX as (1) cell surface receptors expressed on hemocytes that mediate the recognition and phagocytosis or encapsulation of foreign entities through microbe associated molecular patterns (MAMP); (2) signaling receptors activating intracellular signaling cascades or (3) the soluble bridging molecules mediating the linkage between MAMPs and hemocytes [29]. Among those receptors, most are identified as pathogen pattern recognition receptors (PRRs), which include the C-type lectins (CTLs), the scavenger receptors (SRs) and the toll-like receptors (TLRs).

Table 2 Transcripts with annotated functions (GO terms) related to immune recognition, signaling and regulation that were differentially expressed during M. mercenaria focal response against QPX. Additional information on these transcripts is given in Additional file 2. “Inf” designates an infinite fold change calculated for focal response as the expression of that transcript in non-nodule tissue was equal to 0

The C-type mannose receptor-2 (MRC2) identified during the focal response (Table 2) is a member of the C-type lectins (CTLs) superfamily, a large group of Ca2+-dependent carbohydrate-binding proteins that play crucial roles in innate immunity. CTLs recognize pathogens and facilitate their phagocytosis [30, 31] or encapsulation [3234]. MRCs are also key regulators of inflammatory responses and contribute to the removal of harmful inflammatory agents [35, 36]. The 16-fold over-expression of MRC2 during the focal response suggested that active hemocyte encapsulation and local inflammation was induced by QPX at the infection lesions, which is consistent with the results of histopathological observations [1]. Another over-expressed CTL member, the perlucin-like protein, has been previously shown to trigger immune response in Manila clams during microbial infection [37].

Scavenger receptors (SRs) were also among the strongly over-expressed transcripts in infection foci (Table 2). These included somatomedin-b and thrombospondin type-1 domain-containing (RPE spondin), insulin-related peptide receptor, hemicentin-1, lysyl oxidase-like protein 2 and mam domain-containing glycosylphosphatidylinositol anchor protein 1. SRs are structurally diverse PRRs that share the common function of recognizing oxidized or acetylated low-density lipoprotein (LDL) [38]. They contribute to innate immunity by recognizing MAMPs and mediating non-opsonic phagocytosis [29, 39]. They are extensively found on immune cells and are able to interact with both modified-host components and exogenous ligands, which makes SRs a key component in host defense, apoptosis, inflammation and lipoprotein homeostasis [4042]. For example, scallop SRs bind not only to acetylated LDL but also to MAMP including lipopolysaccharides (LPS), peptidoglycans (PGN), mannan and zymosan particles [39]. The sea urchin genome encodes approximately 150 genes consisting of one or more scavenger receptor cysteine-rich (SRCR) domains [43], and the members of this gene family exhibit dynamic shifts in transcription after immune challenge [40, 44].

Our results also show an over-expression of TLR-1 and Toll-8/tollo in nodules (Table 2), which is in agreement with previous investigations showing up-regulation of TLRs in M. mercenaria mantle following QPX challenge [15]. TLRs are among the most ancient and conserved PRRs. They are expressed by immune cells and interact with a large variety of MAMPs. Bivalve TLRs have been characterized in the oyster C. gigas and the scallop C. farreri where they exhibited significant response to LPS stimulation [45, 46]. Transcriptional modulation of TLRs has also been reported in Ruditapes philippinarum and Mya arenaria following MAMPs stimulation and bacterial challenge [47, 48]. Interestingly, Toll-8 (Tollo) has been shown to participate in Drosophila epithelial immunity where it mediates host cells communication that subsequently activates systemic immune responses [49]. This suggests that the Toll pathway could be one of the crucial pivoting links that allow coordination between focal and systemic immune components during infection.

The QPX nodules are formed as the result of granulomatous inflammation, which is a chronic inflammatory reaction characterized by focal accumulation of activated immune cells to isolate the invading agent [50, 51]. The formation of granuloma requires local recruitment of hemocytes at the site of infection to execute extracellular defense processes around the invaders [50]. An array of transcripts associated with cell migration, adhesion and proliferation was regulated in nodules, including G-protein coupled receptors (GPCRs) and integrins families (Table 2). GPCRs regulate inflammatory response via binding to chemokines and chemoattractants, thus activating pathways mediating hemocyte migration and adhesion [52]. They also activate transcription factors in immune cells, thus modulating the synthesis and secretion of certain pro- or anti-inflammatory substances [53]. On the other hand, integrins represent a major group of cell adhesion mediators [54]. They not only modulate the cell-cell and cell-extracellular matrix adhesion, but also affect multiple signal transduction cascades regulating cell survival and proliferation [54]. Overexpression of GPCRs and integrins in nodules suggests their role in hemocytes adhesion and aggregation associated with the formation of granuloma [50, 55].

Several enzymes regulating ROS production were also over-expressed during focal response (Table 3). These included a dual oxidase, which is a key component mediating host-microbe interactions in mucosa [56, 57]. Dual oxidase regulates oxidative burst and ROS production in the gill muscosa of the shrimp Marsupenaeus japonicus, favoring shrimp survivorship during viral infections [58]. Interestingly, transcripts of dual oxidase were only expressed in nodules, suggesting this enzyme was induced upon direct clam-QPX interactions as a part of the mantle mucosa-related immune response. Other transcripts associated with oxidation-reduction processes also exhibited somewhat nodule-exclusive pattern, including the allene oxide synthase-lipoxygenase (AOSL), lysyl oxidase-like protein (LOXL), ww domain-containing oxidoreductase (WWOX), c-terminal binding protein (CtBP), isocitrate dehydrogenase (ICD) and methylenetetrahydrofolate reductase (MTHFR). These molecules are important for maintaining the redox homeostasis of extracellular environment as they are key regulators for oxi-reduction reactions. Over-expression of these transcripts in nodules suggests the need for the host to timely balance out excessive ROS and other toxic intermediates produced during interaction with QPX. In addition to redox-regulation, many of these molecules also take part in the immune modulation indirectly. For example, AOSL play a role in coral immunity by controlling the production of the inflammation regulator arachidonic acid during apoptosis [59], and LOXL acts both as a scavenger receptor and regulator for extracellular matrix remodeling that initiate hemocyte migration and tissue regeneration [60], while WWOX was shown to promote proliferation of immune cells through inhibition of their apoptosis [61, 62]. In addition, ICD, MTHFR and cytochrome p450 are major detoxification enzymes [63, 64]. In fact, immune cells and their secreted effectors require the proper redox state in the extracellular environments to exert their immune functions, which makes the maintenance of redox homeostasis essential for persistent and effective host defense [65, 66]. This is particularly true in the case of QPX disease where the neutralization of parasites depends on extracellular killing pathways [26].

Table 3 Transcripts with putative functions (GO terms) related to immune effectors that were differentially expressed during M. mercenaria focal response against QPX. Additional information on these transcripts is given in Additional file 2. “Inf” designates an infinite fold change calculated for focal response as the expression of that transcript in non-nodule tissue was equal to 0

Apoptosis is an essential host mechanism to effectively remove damaged and infected cells without causing inflammatory destructions to surrounding tissues [67, 68]. Interestingly, apoptosis seems to be largely inhibited during M. mercenaria focal response, as shown by the under-expression of pro-apoptosis transcripts (Table 2). For example, the tumor necrosis factor (TNF)-like protein and 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase-like protein were significantly under-expressed in nodules. Similarly, the pro-apoptotic p53-induced protein and solute carrier family 25 member protein were also under-expressed in nodules. Meanwhile, inhibitor of apoptosis protein (IAP) was over-expressed in nodules. IAPs regulate immune cell expansion and survival in highly inflammatory environments in mammals [69] and they may share similar function in clams by preventing hemocytes from death during interaction with QPX. In fact, ROS production during parasite killing may trigger apoptotic cell death in molluscs [67], and proper control of apoptosis mechanisms is required to maintain cellular homeostasis during immune response. This suspected inhibition of host apoptosis is supported by the above-mentioned over-expression of integrins, as these were shown to protect cells from apoptosis and induce anti-apoptotic pathways during cell adhesion and spreading in the snails Lymnaea stagnalis [70] and B. glabrata [71].

Infection and tissue injury trigger host immune responses via immune signaling pathways [72], by activating transcription factors and initiating the production of immune effectors and regulators. Immune signaling pathways identified in mollusks include Toll, MAPK/JNK and JAK/STAT signaling pathways [27, 47, 73, 74]. During M. mercenaria focal response to QPX, a variety of transcripts encoding kinases and phosphatases were over-expressed (Table 2), suggesting the involvement of MAPK and other kinase-mediated cascades in regulating the focal inflammatory response [75], whereas the under-expression of EF-hand domain containing protein and calmodulin may indicate the suppression of calcium-regulated pathways [76]. Overexpression of E3 ubiquitin-protein ligase and its upstream regulator COP9 signalosome suggests the activation of the damage surveillance ubiquitin/proteasome pathway [7779]. In parallel, the over-expression of rho GTPase and rho kinase suggests the induction of anti-apoptotic Rho-mediated signaling pathway [80, 81] and reinforce the idea that apoptotic inhibition is extensively initiated by M. mercenaria to help fight QPX. However, a very limited number of DE transcripts was detected in relation to those conventional signaling pathways of innate immunity, such as the complement pathway and the Toll/TLR pathway. Only 1 transcript encoding macrophage-expressed gene protein 1 (MPEG1), a putative member of the complement pathway [82, 83], was differentially expressed in QPX nodules. As for the Toll/TLR pathway, only a few receptors were identified (Table 2) but none of their downstream components.

Interestingly, components of the Notch and the Wnt signaling pathways were over-expressed during focal response (Table 2). These included two putative Notch family members, the neurogenic locus Notch protein and the mediator protein nicastrin, and the tyrosin-protein kinase RYK which belongs to the Wnt pathway. Wnt signaling pathway regulates many cellular immune processes and is evolutionarily conserved across taxa [84, 85]. Wnt signaling has been intensively exploited for its regulatory functions during wound healing and tissue regeneration [86, 87], so its over-expression may be related to wound healing to repair damage resulting from tissue digestion by the parasite or tissue necrosis. The Notch signaling network was scarcely explored in bivalves even though it has been reported to be associated with several aspects of immune response in mammals [88], especially in regulating granulomatous reactions to foreign bodies [89]. As an evolutionarily conserved pathway involved in modulating the intercellular signaling, the Notch pathway presumably shares an equally important role in M. mercenaria by modulating the formation of granuloma. Notch signaling triggers macrophage expression of genes involved in pro-inflammatory responses [90], but can suppress inflammation responses triggered by canonical TLR cascade [91], in agreement with our observations. These tightly regulated mechanisms ensure tailored immune responses against different pathogens and are crucial for the host to achieve high immune efficiency while avoiding excessive immune activation and self-inflicted damages.

Several proteases (A.K.A proteinases, peptidases) were also differentially regulated during M. mercenaria focal response against QPX. These proteases mostly belong to the serine and metallo protease families and were generally over-expressed in QPX nodules (Table 3). Proteases serve as key immune modulators partially through their ability to digest and remodel the extracellular matrix and tissues associated with hemocyte activation [92]. Commonly associated with lysosomes and granules of inflammatory cells, serine proteases participate in immune regulation either directly by degradation of pathogens or indirectly through activation of cell surface receptors and signal molecules [9294]. The function of metalloproteases in immune regulation is even more diverse, acting as immune effectors, signal transducers, and mediators of immune cell development and migration [95]. Metalloproteases are also known to be involved in many pro-inflammatory pathways, particularly in the Notch pathway where they act as a type of downstream element to Notch [96].

At the same time, immune effectors with known universal protease inhibitor activities, such as alpha2-macroglobulin (α 2 M) and thioester-containing protein (TEP, a subfamily of α 2 M), were also collectively overexpressed in QPX nodules (Table 3). The α 2 M superfamily inhibits peptidases of diverse origins [97, 98]. The simultaneous over-expression of proteases and protease inhibitors may reflect a finely adjusted defense response of M. mercenaria to maintain homeostasis and regulate self- and pathogen-derived proteases, as shown in other host-pathogen systems [92, 99], including bivalves [100]. Proteases have been identified as major virulence factors of QPX [101, 102], and are thought to degrade host proteinaceous immune effectors and hydrolyze host tissues to fulfill nutritional requirements. Therefore, inhibition of pathogen proteases contribute to host protection, and was shown to represent a determinant factor for resistance against infectious diseases in bivalves [103106]. In parallel, TEPs have been extensively studied in mollusk immunity [99, 107], and beside their function as protease inhibitors, they also play a role as PRRs or opsonins to facilitate microbial phagocytosis and encapsulation. Consistent over-expression of TEPs was noted in this study in agreement with findings following experimental infection with QPX [15]. These results support a critical role of TEPs in clam immune response against QPX either via the protease-inhibitor activity of these proteins, or by mediating parasite encapsulation, or both.

In addition, several metal ion transporters were over-expressed in nodules (Table 3), including the putative copper ion binding protein ceruloplasmin precursor, the transferrin enzyme ferric-chelate reductase and the divalent metal transporter (A.K.A natural resistance-associated macrophage protein; 328 fold increase). These molecules contribute to host defense by controlling the supply of essential micronutrients in the vicinity of infection sites thus reducing parasite survival [108, 109] and favoring the production of antimicrobial factors [110]. It is noteworthy to point out that some of the focally over-expressed transcripts might be partly driven by the dramatic increase of hemocyte proportion within tissues in the vicinity of infection foci as compared to the surrounding host tissues [26, 111]. This would be especially the case for transcripts known to be highly expressed in hemocytes, such as the cell surface PRRs, secreted humoral immune effectors, cell signal transducers and enzymes associated with ROS production.

Systemic response

Significant transcriptomic regulations were observed during M. mercenaria systemic immune response against QPX, with a total of 1,681 DE transcripts, which is about two times the number of focal DE transcripts (829). However, only about one third of the DE transcripts (513) were over-expressed in response to QPX infection, the larger remaining part (1,168) represented significantly under-expressed transcripts (Fig. 4, Additional file 3), possibly due to the chronic stress imposed by the infection. Transcriptome-wide depression has been demonstrated in many marine invertebrates as the result of pathogenic or environmental stress [18, 59, 112114]. A considerable number of systemically under-expressed transcripts was related to metabolism and biomineralization, which could be the result of host resource allocation during on-going infection. The chronic inflammation induced by infection likely created extra energy demands, which require resources being allocated from other physiological processes, such as growth and reproduction, to immune processes which is critical for the survivorship of the host. In fact, slow growth and lower tissue conditions are frequently observed in QPX-infected M. mercenaria [2], and similar energy trade-offs existed between immune defense and other energy expenditure pathways [115118]. Interestingly, a suite of transcripts over-expressed during focal response was significantly under-expressed in non-nodule tissues as compared to healthy clams (Additional file 3), which included several immune effectors and mediators associated with nodule formation and focal inflammation (e.g., integrins, notch proteins and peroxidases). In fact, maintaining high levels of these focally-induced molecules could be costly and dangerous as some are toxic to both the parasite and the host, so their production must be restrained within areas where they can directly exert defense function, and reduced outside the infection foci to minimize risks of undesirable effects on the host [116, 117].

On the other hand, systemically over-expressed transcripts included stress proteins and other soluble immune factors such as lysozyme (c-type lysozyme 2), lectins (C1q domain containing protein, macrophage mannose receptor 1, low affinity immunoglobulin epsilon fc receptor), AMP (hemocyte defensin), proteases (cathepsin K, calpain 11, isoaspartyl peptidase/L-asparaginase, ASRGL, counting factor associated protein d) and ferric-chelate reductase (Table 4, Additional file 3). Over-expression of host stress proteins, such as heat shock proteins (HSP 70, HSP 90) and universal stress protein (USP) was also noted, in agreement with observations made during infection in other bivalve species [119121]. Increased levels of stress proteins provide host cells with protection against incorrect protein folding caused by infection, inflammation, oxidative stress and other destructive events [122, 123]. The systemic over-expression of soluble immune effectors (e.g., humoral proteins) may help maintaining comparatively high immune capacity to prevent the spread of QPX (or secondary pathogens) throughout the host. In addition, transcripts of anti-apoptotic factors (IL17, deoxyguanosine, baculoviral map repeat-containing proteins) were also over-expressed during the systemic response, indicating that anti-apoptotic processes noted during focal response are not limited to the infection foci.

Table 4 Selected transcripts with annotated functions (GO terms) related to immune response that were over-expressed during M. mercenaria systemic response against QPX. Additional information on these transcripts is given in Additional file 3. “Inf” designates an infinite fold change calculated for focal response as the expression of that transcript in non-nodule tissue was equal to 0

Pathway alterations during M. mercenaria’s response to QPX

Transcriptomic alterations during both focal and systemic response discussed above were also highlighted in the pathway enrichment analysis of the DE transcripts. This analysis aims at extracting an overview of phenotypic changes on the underlying functional level, to reduce the complexity of biological information given by the long lists of DE genes/transcripts [124]. The KEGG pathways of focal adhesion (04510), ECM-receptor interaction (04512), Notch signaling pathway (04330) and apoptosis (0421) were significantly over-represented during both focal and systemic response (Fig. 5), even though fold enrichment were generally higher during the focal response. Other immune-related pathways particularly enriched during the focal response included regulation of actin cytoskeleton (04810), cell adhesion molecules (CAMs, 04514), the leukocyte transendothelial migration (04670), complement and coagulation cascade (04610) and Wnt signaling pathway (04310). These pathways are critically involved in the immune cells activation during migration, attachment and parasite encapsulation, which serve as the underlying mechanisms for nodule formation and QPX killing. On the other hand, basic metabolic pathways such as the citrate cycle (TCA) and pyruvate metabolism were specially enriched during the systemic response. These alterations were largely in accordance with the under-expression of metabolism-associated DE transcripts in infected tissue compared to the healthy tissue (Fig. 4), possibly reflecting changes in the energy allocation strategy during infection as discussed above.

Fig. 5
figure 5

Significantly enriched KEGG pathways in M. mercenaria derived from the differentially expressed (DE) genes during focal (a) and systemic (b) response against QPX. The KEGG pathways having significant enrichment (P < 0.01) are presented, and the bar shows the x-fold enrichment of each KEGG pathway. Pathways involved in immune response are marked with red dots while metabolism-related pathways are marked with green dots

Interdependence of KEGG pathways widely exist and most of these are interrelated with each other via shared components, forming a signaling network to allow for pathway crosstalk. To investigate these interactions, we extracted the DE transcripts shared by multiple enriched pathways and constructed a sketch of the hypothetical pathways network that are significantly altered by QPX infection (Fig. 6). In this framework, M. mercenaria response to QPX infection was initiated upon the sensing of danger signals via cell membrane receptors. The signals subsequently transmitted down through the MAPK, Wnt and Notch pathways and triggered the production of a series of host defense factors as the end results. In parallel, activation of pathways regulating actin cytoskeleton and leukocyte transendothelial migration facilitated the recruitment of hemocytes to the infection area to build a barrier of cellular defense against the parasite. Recruited hemocytes then attached to and encapsulated QPX cells as suggested by the modulation of focal adhesion and ECM receptor interaction pathways. These cellular activities were performed under a tight regulation of the apoptosis pathway to determine cell fates, resulting in either the survival or death of M. mercenaria cells.

Fig. 6
figure 6

Overview of immune-related enriched pathways of differentially expressed (DE) transcripts during M. mercenaria response to QPX. a Schematic diagram of enriched pathways and their interactions. Only pathway components encoded by DE transcripts (shown in text boxes) are presented within each enriched KEGG pathway (framed in orthogonal polygons). Arrows display possible interactions (e.g., activation, inhibition) between pathway components. b Overview of DE transcripts expression associated with enriched pathways in nodule, non-nodule and healthy M. mercenaria tissues. The red and green heatmap values indicate log2 fold of relative expression levels for individual transcripts. Arp2/3: actin related protein 2/3 complex; CtBP: C-terminal binding protein; CAMK2: calcium/calmodulin-dependent protein kinase (CaM kinase) II; CycD: cyclin D1(CCND1); ECM: von Willebrand factor; ENDO-G: endonuclease G; FAK: focal adhesion kinase; FYN: tyrosine-protein kinase; HES1: hairy and enhancer of split 1; IAP: baculoviral IAP repeat-containing protein, cIAPs; ITGA: intergrin alpha; ITGB: intergrin beta; MLCP: serine/threonine-protein phosphatase PP1 catalytic subunit; NRXN: neurexin; NCSTN: nicastrin; Rac: Ras-related C3 botulinum toxin substrate; ROCK: Rho-associated protein kinase; RTK: proto-oncogene tyrosine-protein kinase; SELE: selectin, endothelial cell; SELP: selectin, platelet; TRAIL: tumor necrosis factor ligand superfamily member 10); WNT: wingless-type MMTV integration site family

Distinctive transcriptomic pattern of healthy clams

A suite of transcripts (407, Fig. 4, Additional file 4 and Table 5) exhibited higher transcription levels in healthy clams as compared to diseased clams (considering both nodule and non-nodule tissues). A considerable fraction of these transcripts were related to metabolic processes, nucleic acids binding and transcriptional regulation. The over-expressed immune-related transcripts identified in healthy clams are of particular interest as they may be involved in M. mercenaria resistance towards QPX. For example, an antimicrobial protein (aplysianin A [125, 126]) was exclusively identified in healthy clams with almost no detection in diseased clams. In addition, the highest expression levels of a serine protease inhibitor were also observed in healthy clams. A serine protease inhibitor has been linked to oyster (Crassostrea virginica) resistance against the protozoan parasite Perkinsus marinus [103106, 127]. Therefore, the high expression of the serine protease inhibitor in healthy clams supports its involvement in M. mercenaria’s resistance against QPX, likely by inhibiting the activity of parasite proteases. Moreover, a pathogen recognition protein (c-type lectin domain family 10 member A-like) was significantly higher in healthy clams as compared to diseased animals, which may also contribute to clam resistance against QPX by promoting microbial recognition and encapsulation [128]. Previous studies demonstrated that clam genetic background affects M. mercenaria resistance toward QPX [4, 129, 130], therefore immune-related transcripts specifically associated with healthy clams represent excellent candidates as molecular markers for further research on screening and breeding QPX-resistant M. mercenaria strains.

Table 5 Selected transcripts with annotated functions (GO terms) related to immune response that were differentially expressed in naïve M. mercenaria as compared to QPX infected individuals. Additional information on these transcripts is given in Additional file 4. “Inf” designates an infinite fold change calculated for focal response as the expression of that transcript in non-nodule tissue was equal to 0

On the other hand, a subset of transcripts (126, Fig. 4 and Table 5) exhibited lowest expression levels in healthy clams. These include stress protein HSP 70 and USP, as well as the detoxification molecule cytochrome p450, which together highlight the stress experienced by clams as the result of QPX infection [63, 131135]. In addition, other immune related transcripts including a protease (isoaspartyl peptidase/L-asparaginase, ASRGL), a protease inhibitor (GTP-binding protein yptV4), a C-type lectin (MRC1), and molecules involved in tissue regeneration and cell signaling, were also under-expressed in healthy clams as compared to infected individuals, suggesting their role in fighting the infection.

Conclusions

This is one of the first studies contrasting focal and systemic immune responses to infections in invertebrates using high-throughput sequencing. Resulting transcriptome represents a significant addition to the so far limited public genomic information available for M. mercenaria. The transcriptomic profiles of healthy and infected clams reflected complex interactions between the host immune system and the pathogen leading to molecular changes at both the infection foci and the systemic level. In general, the systemic responses of M. mercenaria reflected prevailing transcriptomic suppression accompanied with a contrasting over-expression of stress proteins and soluble antimicrobial effectors; whereas the focal response highlighted cell-cell interactions between hemocytes and the parasite that typically result in local inflammation, extracellular degradation, encapsulation, granuloma formation, and wound repair. What needs to be kept in mind is that the regulation of these genes can be the result of either an effective immune response or a symptom of a future death. In fact, several apoptotic and anti-apoptotic genes were regulated. This highlights a fine adjustment of M. mercenaria defense mechanisms to precisely adapt to the infection development (e.g., through the modulation of energy allocation, apoptotic and anti-apoptotic processes and mobilization of different signaling pathways). In addition, the identification of immune-related transcripts that were particularly associated with healthy clams offered new perspectives on the molecular features putatively involved in clam resistance against QPX.

Methods

Clam tissue and RNA samples preparation

Adult hard clams (54 ± 5 mm in length, mean ± standard deviation) were collected from a QPX enzootic area in Massachusetts. Clams were shucked and grossly examined for the presence of nodules along the mantle rim. Nodule tissues were dissected and a small piece of each nodule was microscopically examined (fresh mount) for the identification of QPX cells. Aliquot tissues from positive nodules were submitted to RNA extraction. Meanwhile, a piece of seemingly healthy tissue that is anatomically symmetrical to the nodule was collected from the same clam (e.g., “healthy” tissue from a diseased clam) and divided into 2 aliquots, with the first aliquot used for QPX detection by qPCR [10] and the second used for RNA extraction. Mantle tissues were also collected from “seemingly” healthy clams (no visible nodules) and used for QPX detection and RNA extraction. Following confirmation of disease status by qPCR, the samples were divided into the following 3 categories: (1) infected tissue from a diseased clam, (2) “healthy”(non-nodule) tissue from a diseased clam, and (3) “healthy” tissue from a healthy clam. Total RNA was individually extracted using TRIzol ® Reagent (Invitrogen, Carlsbad, CA, USA). Further RNA clean-up and on-column DNase digestion were performed with RNeasy Mini Kit (Qiagen Ltd., Crawley, UK) according to the manufacturer’s guideline. RNA quantity and quality were analyzed on Nanodrop ND-1000 (Thermo Scientific, Wilmington, USA). Only RNA samples with absorption ratios of A260/A280 close to 2.0 were used for RNA-seq analysis. A total amount of 3 μg RNAs per sample was used and pooled into a total of 3 pools (3 clams per pool) representing each tissue category (9 pools in all). The pooling strategy is shown in Table 1.

RNA sequencing, de novo assembly and annotation

The sequencing of each pooled RNA sample as a paired end (PE) reads library (100 pb) was performed on Illumina HiSeq 2000 platform at the McGill University and Genome Quebec Innovation Center (Montreal, Canada), producing from 27.1 up to 46.8 millions of PE reads per sample (Table 1). Raw reads were filtered and trimmed according to length and quality score (min length 60 nt, end trimming quality 25, min quality filtering: 20 on 75 % of the read length) using the FASTX-Toolkit software v 0.0.13 (http://hannonlab.cshl.edu/fastx_toolkit/). rRNA cleaning was performed using the riboPicker software v 4.0.3 (http://ribopicker.sourceforge.net) [136] against SILVA database v111. And finally pair retrieval was performed using a homemade python script accessible through GitHub (https://github.com/ppericard/bioinfo-toolkit). High quality filtered sequence reads from all libraries were combined and subsequently used for de novo assembly. Assembly was based on the de Bruijn graph assembler Trinity (http://trinityrnaseq.sourceforge.net/) [137] using the default parameters. Assembly quality was controlled by re-mapping raw reads back to the transcripts using bowtie 2 and RSEM within Trinity package scripts. The assembled transcripts with sequence length longer than 200 bp, re-mapping FPKM (fragments per kilobase of transcript per million mapped reads) greater than 1 and isoform discovery level greater than 1 % were then considered for the following annotation and transcript abundance quantitation. Annotation of this de novo assembled transcriptome was performed using Blast2GO (http://www.blast2go.com) with a semi-automated functional annotation based on sequence homology search. Putative gene identities was obtained by Blastx search against National Center for Biotechnology Information (NCBI) non-redundant sequences (nr) database with the E-value threshold setting at 1E-03. Putative gene functions were predicted by sequence similarity search against Gene Ontology (GO, http://www.geneontology.org/) database and assigning GO annotation terms to each mapped transcript. Protein domain search and enzyme annotation were also performed using InterPro scan and the Kyoto Encyclopedia of Genes and Genomes (KEGG). KEGG Orthology (KO) terms and KEGG pathways were also assigned to the assembled sequences using the online KEGG Automatic Annotation Server (KAAS, http://www.genome.jp/kegg/kaas/) using bi-directional best-hit methods [138]. This server provides KO annotation and pathway mapping.

Differential gene expression analysis

The consensus transcriptome generated in the previous steps was used as a reference for transcript abundance analysis. Reads from each single library were separately mapped to this reference transcriptome using the “RNA-Seq by Expectation-Maximization (RSEM)” method that is bundled within the Trinity package. The expression level of each transcript was determined as the total mapped reads count. The differences in gene expression between clam tissue samples (nodule, non-nodule and healthy tissues) were estimated using the DESeq Bioconductor package (https://github.com/Bioconductor-mirror/SESeq/tree/release-3.2) [139] in R statistical software (R Development Core Team, 2010; http://www.R-project.org). The threshold for defining significant differentially expressed (DE) transcripts between two different conditions (3 replications in each condition) was set as adjusted p-value smaller than 0.001 and absolute log2 (fold change) values greater than 2. Expression patterns of DE transcripts were also analyzed by a K-means clustering method using Euclidean distance based on expression levels over all input samples. For further analysis, only those DE transcripts with annotation were considered as candidates of interest and were subsequently divided into over- and under-expressed groups.

Availability of supporting data

The data sets supporting the results of this study are available at the NCBI short Read Archive database under the SRA accession number SRP068241.

References

  1. Whyte SK, Cawthorn RJ, McGladdery SE. QPX (Quahaug Parasite X), a pathogen of northern quahaug Mercenaria mercenaria from the gulf of St-Lawrence, Canada. Dis Aquat Organ. 1994;19(2):129–36.

    Article  Google Scholar 

  2. Smolowitz R, Leavitt D, Perkins F. Observations of a protistan disease similar to QPX in Mercenaria mercenaria (hard clams) from the coast of Massachusetts. J Invertebr Pathol. 1998;71(1):9–25.

    Article  PubMed  Google Scholar 

  3. Anderson RS, Kraus BS, McGladdery SE, Reece KS, Stokes NA. A thraustochytrid protist isolated from Mercenaria mercenaria: molecular characterization and host defense responses. Fish Shellfish Immun. 2003;15(3):183–94.

    Article  CAS  Google Scholar 

  4. Ford SE, Kraeuter JN, Barber RD, Mathis G. Aquaculture-associated factors in QPX disease of hard clams: density and seed source. Aquaculture. 2002;208(1–2):23–38.

    Article  Google Scholar 

  5. Ragone Calvo LM, Walker JG, Burreson EM. Prevalence and distribution of QPX, Quahog Parasite Unknown, in hard clams Mercenaria mercenaria in Virginia, USA. Dis Aquat Organ. 1998;33(3):209–19.

    Article  PubMed  CAS  Google Scholar 

  6. MacCallum GS, McGladdery SE. Quahog Parasite Unknown (QPX) in the northern quahog Mercenaria mercenaria (Linnaeus, 1758) and M. mercenaria var. notata from Atlantic Canada, survey results from three maritime provinces. J Shellfish Res. 2000;19(1):43–50.

    Google Scholar 

  7. Dove ADM, Bowser PR, Cerrato RM. Histological Analysis of an Outbreak of QPX Disease in Wild Hard Clams Mercenaria mercenaria in New York. J Aquat Anim Health. 2004;16(4):246–50.

    Article  Google Scholar 

  8. Lyons MM, Ward JE, Smolowitz R, Uhlinger KR, Gast RJ. Lethal marine snow: Pathogen of bivalve mollusc concealed in marine aggregates. Limnol Oceanogr. 2005;50(6):1983–8.

    Article  Google Scholar 

  9. Gast RJ, Moran DM, Audemard C, Lyons MM, DeFavari J, Reece KS, et al. Environmental distribution and persistence of Quahog Parasite Unknown (QPX). Dis Aquat Organ. 2008;81(3):219–29.

    Article  PubMed  Google Scholar 

  10. Liu QQ, Allam B, Collier JL. Quantitative Real-Time PCR Assay for QPX (Thraustochytriidae), a Parasite of the Hard Clam (Mercenaria mercenaria). Appl Environ Microb. 2009;75(14):4913–8.

    Article  CAS  Google Scholar 

  11. Perrigault M, Dahl SF, Espinosa EP, Gambino L, Allam B. Effects of temperature on hard clam (Mercenaria mercenaria) immunity and QPX (Quahog Parasite Unknown) disease development: II. Defense parameters. J Invertebr Pathol. 2011;106(2):322–32.

    Article  PubMed  Google Scholar 

  12. Bayne CJ. Phagocytosis and non-self recognition in invertebrates. Bioscience. 1990;40(10):723–31.

    Article  Google Scholar 

  13. Canesi L, Gallo G, Gavioli M, Pruzzo C. Bacteria-hemocyte interactions and phagocytosis in marine bivalves. Microsc Res Tech. 2002;57(6):469–76.

    Article  PubMed  Google Scholar 

  14. Perrigault M, Dahl SF, Espinosa EP, Allam B. Effects of salinity on hard clam (Mercenaria mercenaria) defense parameters and QPX disease dynamics. J Invertebr Pathol. 2012;110(1):73–82.

    Article  PubMed  CAS  Google Scholar 

  15. Perrigault M, Tanguy A, Allam B. Identification and expression of differentially expressed genes in the hard clam, Mercenaria mercenaria, in response to quahog parasite unknown (QPX). Bmc Genomics. 2009;10:377.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  16. Perrigault M, Allam B. Differential immune response in the hard clam (Mercenaria mercenaria) against bacteria and the protistan pathogen QPX (quahog parasite unknown). Fish Shellfish Immun. 2012;32(6):1124–34.

    Article  CAS  Google Scholar 

  17. Leite RB, Milan M, Coppe A, Bortoluzzi S, dos Anjos A, Reinhardt R, et al. mRNA-Seq and microarray development for the Grooved carpet shell clam, Ruditapes decussatus: a functional approach to unravel host -parasite interaction. Bmc Genomics. 2013;14(1):741.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  18. Zhao X, Yu H, Kong L, Li Q. Transcriptomic Responses to Salinity Stress in the Pacific Oyster Crassostrea gigas. PLoS ONE. 2012;7(9):e46244.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  19. Teaniniuraitemoana V, Huvet A, Levy P, Klopp C, Lhuillier E, Gaertner-Mazouni N, et al. Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes. Bmc Genomics. 2014;15(1):491.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  20. Adema CM, Hanington PC, Lun C-M, Rosenberg GH, Aragon AD, Stout BA, et al. Differential transcriptomic responses of Biomphalaria glabrata (Gastropoda, Mollusca) to bacteria and metazoan parasites, Schistosoma mansoni and Echinostoma paraensei (Digenea, Platyhelminthes). Mol Immunol. 2010;47(4):849–60.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  21. Wang JP, Lindsay BG, Leebens-Mack J, Cui L, Wall K, Miller WC, et al. EST clustering error evaluation and correction. Bioinformatics. 2004;20(17):2973–84.

    Article  PubMed  CAS  Google Scholar 

  22. Mittapalli O, Bai X, Mamidala P, Rajarapu SP, Bonello P, Herms DA. Tissue-Specific Transcriptomics of the Exotic Invasive Insect Pest Emerald Ash Borer (Agrilus planipennis). PLoS ONE. 2010;5(10):e13708.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  23. Pelechano V, Wei W, Steinmetz LM. Extensive transcriptional heterogeneity revealed by isoform profiling. Nature. 2013;497(7447):127–31.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  24. Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10(1):161.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  25. Humphries JE, Yoshino TP. Cellular receptors and signal transduction in molluscan hemocytes: connections with the innate immune system of vertebrates. Integr Comp Biol. 2003;43(2):305–12.

    Article  PubMed  CAS  Google Scholar 

  26. Soudant P, E Chu FL, Volety A. Host-parasite interactions: Marine bivalve molluscs and protozoan parasites, Perkinsus species. J Invertebr Pathol. 2013;114(2):196–216.

    Article  PubMed  Google Scholar 

  27. Canesi L, Betti M, Ciacci C, Lorusso LC, Pruzzo C, Gallo G. Cell signalling in the immune response of mussel hemocytes. Invertebrate Survival Journal: ISJ. 2006;3:40–9.

    Google Scholar 

  28. Hatanaka R, Sekine Y, Hayakawa T, Takeda K, Ichijo H. Signaling pathways in invertebrate immune and stress response. Invertebrate Survival Journal: ISJ. 2009;6:32–43.

    Google Scholar 

  29. Jeannin P, Jaillon S, Delneste Y. Pattern recognition receptors in the immune response against dying cells. Curr Opin Immunol. 2008;20(5):530–7.

    Article  PubMed  CAS  Google Scholar 

  30. Wang L, Wang L, Zhang D, Li F, Wang M, Huang M, et al. A novel C-type lectin from crab Eriocheir sinensis functions as pattern recognition receptor enhancing cellular encapsulation. Fish Shellfish Immun. 2013;34(3):832–42.

    Article  CAS  Google Scholar 

  31. Vazquez-Mendoza A, Carrero JC, Rodriguez-Sosa M. Parasitic infections: a role for C-type lectins receptors. Biomed Res Int. 2013;2013:456352.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  32. Yang J, Qiu L, Wei X, Wang L, Wang L, Zhou Z, et al. An ancient C-type lectin in Chlamys farreri (CfLec-2) that mediate pathogen recognition and cellular adhesion. Dev Comp Immunol. 2010;34(12):1274–82.

    Article  PubMed  CAS  Google Scholar 

  33. Wang L, Wang L, Yang J, Zhang H, Huang M, Kong P, et al. A multi-CRD C-type lectin with broad recognition spectrum and cellular adhesion from Argopecten irradians. Dev Comp Immunol. 2012;36(3):591–601.

    Article  PubMed  CAS  Google Scholar 

  34. Ao J, Ling E, Yu XQ. Drosophila C-type lectins enhance cellular encapsulation. Mol Immunol. 2007;44(10):2541–8.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  35. Lasky LA. Selectin-Carbohydrate Interactions and the Initiation of the Inflammatory Response. Annu Rev Biochem. 1995;64(1):113–40.

    Article  PubMed  CAS  Google Scholar 

  36. Gazi U, Martinez-Pomares L. Influence of the mannose receptor in host immune responses. Immunobiology. 2009;214(7):554–61.

    Article  PubMed  CAS  Google Scholar 

  37. Moreira R, Milan M, Balseiro P, Romero A, Babbucci M, Figueras A, et al. Gene expression profile analysis of Manila clam (Ruditapes philippinarum) hemocytes after a Vibrio alginolyticus challenge using an immune-enriched oligo-microarray. Bmc Genomics. 2014;15(1):267.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  38. Gordon S. Pattern recognition receptors: doubling up for the innate immune response. Cell. 2002;111(7):927–30.

    Article  PubMed  CAS  Google Scholar 

  39. Liu L, Yang J, Qiu L, Wang L, Zhang H, Wang M, et al. A novel scavenger receptor-cysteine-rich (SRCR) domain containing scavenger receptor identified from mollusk mediated PAMP recognition and binding. Dev Comp Immunol. 2011;35(2):227–39.

    Article  PubMed  CAS  Google Scholar 

  40. Mukhopadhyay S, Gordon S. The role of scavenger receptors in pathogen recognition and innate immunity. Immunobiology. 2004;209(1–2):39–49.

    Article  PubMed  CAS  Google Scholar 

  41. Blumbach B, Pancer Z, Diehl-Seifert B, Steffen R, Munkner J, Muller I, et al. The putative sponge aggregation receptor. Isolation and characterization of a molecule composed of scavenger receptor cysteine-rich domains and short consensus repeats. J Cell Sci. 1998;111(Pt 17):2635–44.

    PubMed  CAS  Google Scholar 

  42. van der Laan LJ, Dopp EA, Haworth R, Pikkarainen T, Kangas M, Elomaa O, et al. Regulation and functional involvement of macrophage scavenger receptor MARCO in clearance of bacteria in vivo. J Immunol. 1999;162(2):939–47.

    PubMed  Google Scholar 

  43. Bowdish DM, Gordon S. Conserved domains of the class A scavenger receptors: evolution and function. Immunol Rev. 2009;227(1):19–31.

    Article  PubMed  CAS  Google Scholar 

  44. Rast JP, Smith LC, Loza-Coll M, Hibino T, Litman GW. Genomic insights into the immune system of the sea urchin. Science. 2006;314(5801):952–6.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  45. 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.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  46. 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;22(5):451–66.

    Article  PubMed  CAS  Google Scholar 

  47. Moreira R, Balseiro P, Planas JV, Fuste B, Beltran S, Novoa B, et al. Transcriptomics of in vitro immune-stimulated hemocytes from the Manila clam Ruditapes philippinarum using high-throughput sequencing. PLoS One. 2012;7(4):e35009.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  48. 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;34(7):710–4.

    Article  PubMed  CAS  Google Scholar 

  49. Akhouayri I, Turc C, Royet J, Charroux B. Toll-8/Tollo Negatively Regulates Antimicrobial Response in the Drosophila Respiratory Epithelium. Plos Pathog. 2011;7(10):e1002319.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  50. Williams GT, Williams WJ. Granulomatous inflammation--a review. J Clin Pathol. 1983;36(7):723–33.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  51. de Brito T, Franco MF. Granulomatous inflammation. Revista do Instituto de Medicina Tropical de Sao Paulo. 1994;36(2):185–92.

    Article  PubMed  Google Scholar 

  52. Sun L, Ye RD. Role of G protein-coupled receptors in inflammation. Acta pharmacologica Sinica. 2012;33(3):342–50.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  53. Cho H, Kehrl JH. Regulation of immune function by G protein-coupled receptors, trimeric G proteins, and RGS proteins. Prog Mol Biol Transl Sci. 2009;86:249–98.

    Article  PubMed  CAS  Google Scholar 

  54. Zhang Y, Wang H. Integrin signalling and function in immune cells. Immunology. 2012;135(4):268–75.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  55. Stavitsky AB. Regulation of granulomatous inflammation in experimental models of schistosomiasis. Infect Immun. 2004;72(1):1–12.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  56. Ha EM, Oh CT, Bae YS, Lee WJ. A direct role for dual oxidase in Drosophila gut immunity. Science. 2005;310(5749):847–50.

    Article  PubMed  CAS  Google Scholar 

  57. Bae YS, Choi MK, Lee WJ. Dual oxidase in mucosal immunity and host-microbe homeostasis. Trends Immunol. 2010;31(7):278–87.

    Article  PubMed  CAS  Google Scholar 

  58. Inada M, Kihara K, Kono T, Sudhakaran R, Mekata T, Sakai M, et al. Deciphering of the Dual oxidase (Nox family) gene from kuruma shrimp, Marsupenaeus japonicus: Full-length cDNA cloning and characterization. Fish Shellfish Immun. 2013;34(2):471–85.

    Article  CAS  Google Scholar 

  59. Libro S, Kaluziak ST, Vollmer SV. RNA-seq Profiles of Immune Related Genes in the Staghorn Coral Acropora cervicornis Infected with White Band Disease. PLoS ONE. 2013;8(11):e81821.

    Article  PubMed Central  PubMed  Google Scholar 

  60. Bignon M, Pichol-Thievend C, Hardouin J, Malbouyres M, Brechot N, Nasciutti L, et al. Lysyl oxidase-like protein-2 regulates sprouting angiogenesis and type IV collagen assembly in the endothelial basement membrane. Blood. 2011;118(14):3979–89.

    Article  PubMed  CAS  Google Scholar 

  61. Chang NS, Hsu LJ, Lin YS, Lai FJ, Sheu HM. WW domain-containing oxidoreductase: a candidate tumor suppressor. Trends Mol Med. 2007;13(1):12–22.

    Article  PubMed  CAS  Google Scholar 

  62. Yang Z, Zhao T, Liu Y. Upregulation of tumor suppressor WWOX promotes immune response in glioma. Cell Immunol. 2013;285(1–2):1–5.

    Article  PubMed  CAS  Google Scholar 

  63. Behrendt L, Jonsson ME, Goldstone JV, Stegeman JJ. Induction of cytochrome P450 1 genes and stress response genes in developing zebrafish exposed to ultraviolet radiation. Aquat Toxicol. 2010;98(1):74–82.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  64. Pan L, Liu N, Xu C, Miao J. Identification of a novel P450 gene belonging to the CYP4 family in the clam Ruditapes philippinarum, and analysis of basal- and benzo(a)pyrene-induced mRNA expression levels in selected tissues. Environ Toxicol Pharmacol. 2011;32(3):390–8.

    Article  PubMed  CAS  Google Scholar 

  65. Lotze MT, Zeh HJ, Rubartelli A, Sparvero LJ, Amoscato AA, Washburn NR, et al. The grateful dead: damage-associated molecular pattern molecules and reduction/oxidation regulate immunity. Immunol Rev. 2007;220:60–81.

    Article  PubMed  CAS  Google Scholar 

  66. Mone Y, Monnin D, Kremer N. The oxidative environment: a mediator of interspecies communication that drives symbiosis evolution. Proc Biol Sci/The Royal Society. 2014;281(1785):20133112.

    Article  CAS  Google Scholar 

  67. Terahara K, Takahashi KG. Mechanisms and immunological roles of apoptosis in molluscs. Curr Pharm Des. 2008;14(2):131–7.

    Article  PubMed  CAS  Google Scholar 

  68. Sokolova IM. Apoptosis in molluscan immune defense. Invertebrate Survival Journal: ISJ. 2009;6:49–58.

    Google Scholar 

  69. Gentle IE, Moelter I, Lechler N, Bambach S, Vucikuja S, Hacker G, et al. Inhibitors of apoptosis proteins (IAPs) are required for effective T-cell expansion/survival during antiviral immunity in mice. Blood. 2014;123(5):659–68.

    Article  PubMed  CAS  Google Scholar 

  70. Plows LD, Cook RT, Davies AJ, Walker AJ. Phagocytosis by Lymnaea stagnalis haemocytes: A potential role for phosphatidylinositol 3-kinase but not protein kinase A. J Invertebr Pathol. 2006;91(1):74–7.

    Article  PubMed  CAS  Google Scholar 

  71. Yoshino TP, Laursen JR. Production of Schistosoma mansoni daughter sporocysts from mother sporocysts maintained in synxenic culture with Biomphalaria glabrata embryonic (Bge) cells. J Parasitol. 1995;81(5):714–22.

    Article  PubMed  CAS  Google Scholar 

  72. Newton K, Dixit VM. Signaling in innate immunity and inflammation. Cold Spring Harb Perspect Biol. 2012;4(3). doi: 10.1101/cshperspect.a006049.

  73. Toubiana M, Rosani U, Giambelluca S, Cammarata M, Gerdol M, Pallavicini A, et al. Toll signal transduction pathway in bivalves: complete cds of intermediate elements and related gene transcription levels in hemocytes of immune stimulated Mytilus galloprovincialis. Dev Comp Immunol. 2014;45(2):300–12.

    Article  PubMed  CAS  Google Scholar 

  74. Philipp EER, Kraemer L, Melzner F, Poustka AJ, Thieme S, Findeisen U, et al. Massively Parallel RNA Sequencing Identifies a Complex Immune Gene Repertoire in the lophotrochozoan Mytilus edulis. PLoS ONE. 2012;7(3):e33091.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  75. Theodosiou A, Ashworth A. MAP kinase phosphatases. Genome Biol. 2002;3(7):REVIEWS3009.

    Article  PubMed Central  PubMed  Google Scholar 

  76. Grabarek Z. Insights into modulation of calcium signaling by magnesium in calmodulin, troponin C and related EF-hand proteins. Biochim Biophys Acta (BBA) - Molecular Cell Research. 2011;1813(5):913–21.

    Article  CAS  Google Scholar 

  77. Liu Y, Schiff M, Serino G, Deng XW, Dinesh-Kumar SP. Role of SCF ubiquitin-ligase and the COP9 signalosome in the N gene-mediated resistance response to Tobacco mosaic virus. The Plant cell. 2002;14(7):1483–96.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  78. Bech-Otschir D, Seeger M, Dubiel W. The COP9 signalosome: at the interface between signal transduction and ubiquitin-dependent proteolysis. J Cell Sci. 2002;115(Pt 3):467–73.

    PubMed  CAS  Google Scholar 

  79. Nezames CD, Deng XW. The COP9 signalosome: its regulation of cullin-based E3 ubiquitin ligases and role in photomorphogenesis. Plant Physiol. 2012;160(1):38–46.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  80. Tosello-Trampont AC, Nakada-Tsukui K, Ravichandran KS. Engulfment of apoptotic cells is negatively regulated by Rho-mediated signaling. J Biol Chem. 2003;278(50):49911–9.

    Article  PubMed  CAS  Google Scholar 

  81. Krijnen PA, Sipkens JA, Molling JW, Rauwerda JA, Stehouwer CD, Muller A, et al. Inhibition of Rho-ROCK signaling induces apoptotic and non-apoptotic PS exposure in cardiomyocytes via inhibition of flippase. J Mol Cell Cardiol. 2010;49(5):781–90.

    Article  PubMed  CAS  Google Scholar 

  82. He X, Zhang Y, Yu Z. An Mpeg (macrophage expressed gene) from the Pacific oyster Crassostrea gigas: molecular characterization and gene expression. Fish Shellfish Immunol. 2011;30(3):870–6.

    Article  PubMed  CAS  Google Scholar 

  83. D’Angelo ME, Dunstone MA, Whisstock JC, Trapani JA, Bird PI. Perforin evolved from a gene duplication of MPEG1, followed by a complex pattern of gene gain and loss within Euteleostomi. Bmc Evol Biol. 2012;12(1):59.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  84. Staal FJ, Luis TC, Tiemessen MM. WNT signalling in the immune system: WNT is spreading its wings. Nat Rev Immunol. 2008;8(8):581–93.

    Article  PubMed  CAS  Google Scholar 

  85. Holstein TW. The evolution of the Wnt pathway. Cold Spring Harb Perspect Biol. 2012;4(7):a007922.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  86. Zhang Z, Deb A, Zhang Z, Pachori A, He W, Guo J, et al. Secreted frizzled related protein 2 protects cells from apoptosis by blocking the effect of canonical Wnt3a. J Mol Cell Cardiol. 2009;46(3):370–7.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  87. Fathke C, Wilson L, Shah K, Kim B, Hocking A, Moon R, et al. Wnt signaling induces epithelial differentiation during cutaneous wound healing. Bmc Cell Biol. 2006;7(1):4.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  88. Yuan JS, Kousis PC, Suliman S, Visan I, Guidos CJ. Functions of notch signaling in the immune system: consensus and controversies. Annu Rev Immunol. 2010;28:343–65.

    Article  PubMed  CAS  Google Scholar 

  89. Rangel JR, Chung Y, Rosenbach M, Ziober AF, Zhang PJ, Schaffer A. Expression of Notch signaling components in cutaneous foreign body and sarcoidal granulomas and fusing macrophages. Am J Dermatopathol. 2014;36(5):409–13.

    Article  PubMed Central  PubMed  Google Scholar 

  90. Palaga T, Buranaruk C, Rengpipat S, Fauq AH, Golde TE, Kaufmann SH, et al. Notch signaling is activated by TLR stimulation and regulates macrophage functions. Eur J Immunol. 2008;38(1):174–83.

    Article  PubMed  CAS  Google Scholar 

  91. Zhang Q, Wang C, Liu Z, Liu X, Han C, Cao X, et al. Notch signal suppresses Toll-like receptor-triggered inflammatory responses in macrophages by inhibiting extracellular signal-regulated kinase 1/2-mediated nuclear factor kappaB activation. J Biol Chem. 2012;287(9):6208–17.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  92. Manoury B, Roghanian A, Sallenave JM. Serine and cysteine proteases and their inhibitors as antimicrobial agents and immune modulators. In: Vergnolle N, Chignard M, editors. Proteases and their receptors in inflammation. Progress in Inflammation Research. Basel: Springer; 2011. p. 27–50.

  93. Shpacovitch V, Feld M, Hollenberg MD, Luger TA, Steinhoff M. Role of protease-activated receptors in inflammatory responses, innate and adaptive immunity. J Leukoc Biol. 2008;83(6):1309–22.

    Article  PubMed  CAS  Google Scholar 

  94. Meyer-Hoffert U, Wiedow O. Neutrophil serine proteases: mediators of innate immune responses. Curr Opin Hematol. 2011;18(1):19–24.

    Article  PubMed  CAS  Google Scholar 

  95. Khokha R, Murthy A, Weiss A. Metalloproteinases and their natural inhibitors in inflammation and immunity. Nat Rev Immunol. 2013;13(9):649–65.

    Article  PubMed  CAS  Google Scholar 

  96. Murthy A, Shao YW, Narala SR, Molyneux SD, Zuniga-Pflucker JC, Khokha R. Notch activation by the metalloproteinase ADAM17 regulates myeloproliferation and atopic barrier immunity by suppressing epithelial cytokine synthesis. Immunity. 2012;36(1):105–19.

    Article  PubMed  CAS  Google Scholar 

  97. Sottrup-Jensen L. Alpha-macroglobulins: structure, shape, and mechanism of proteinase complex formation. J Biol Chem. 1989;264(20):11539–42.

    PubMed  CAS  Google Scholar 

  98. Borth W. Alpha 2-macroglobulin, a multifunctional binding protein with targeting characteristics. Faseb J. 1992;6(15):3345–53.

    PubMed  CAS  Google Scholar 

  99. Armstrong PB. Proteases and protease inhibitors: a balance of activities in host–pathogen interaction. Immunobiology. 2006;211(4):263–81.

    Article  PubMed  CAS  Google Scholar 

  100. Allam B, Pales Espinosa E, Tanguy A, Jeffroy F, Le Bris C, Paillard C. Transcriptional changes in Manila clam (Ruditapes philippinarum) in response to Brown Ring Disease. Fish Shellfish Immunol. 2014;41(1):2–11.

    Article  PubMed  CAS  Google Scholar 

  101. Anderson RS, Kraus BS, McGladdery S, Smolowitz R. QPX, a pathogen of quahogs (hard clams), employs mucoid secretions to resist host antimicrobial agents. J Shellfish Res. 2003;22(1):205–8.

    Google Scholar 

  102. Rubin E, Tanguy A, Perrigault M, Pales Espinosa E, Allam B. Characterization of the transcriptome and temperature-induced differential gene expression in QPX, the thraustochytrid parasite of hard clams. Bmc Genomics. 2014;15(1):245.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  103. He Y, Yu H, Bao Z, Zhang Q, Guo X. Mutation in promoter region of a serine protease inhibitor confers Perkinsus marinus resistance in the eastern oyster (Crassostrea virginica). Fish Shellfish Immunol. 2012;33(2):411–7.

    Article  PubMed  CAS  Google Scholar 

  104. La Peyre JF, Xue QG, Itoh N, Li Y, Cooper RK. Serine protease inhibitor cvSI-1 potential role in the eastern oyster host defense against the protozoan parasite Perkinsus marinus. Dev Comp Immunol. 2010;34(1):84–92.

    Article  PubMed  CAS  Google Scholar 

  105. Yu H, He Y, Wang X, Zhang Q, Bao Z, Guo X. Polymorphism in a serine protease inhibitor gene and its association with disease resistance in the eastern oyster (Crassostrea virginica Gmelin). Fish Shellfish Immunol. 2011;30(3):757–62.

    Article  PubMed  CAS  Google Scholar 

  106. Xue QG, Waldrop GL, Schey KL, Itoh N, Ogawa M, Cooper RK, et al. A novel slow-tight binding serine protease inhibitor from eastern oyster (Crassostrea virginica) plasma inhibits perkinsin, the major extracellular protease of the oyster protozoan parasite Perkinsus marinus. Comp Biochem Physiol B Biochem Mol Biol. 2006;145(1):16–26.

    Article  PubMed  CAS  Google Scholar 

  107. 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;4(9):e813.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  108. Ong ST, Shan Ho JZ, Ho B, Ding JL. Iron-withholding strategy in innate immunity. Immunobiology. 2006;211(4):295–314.

    Article  PubMed  CAS  Google Scholar 

  109. Stafford SL, Bokil NJ, Achard ME, Kapetanovic R, Schembri MA, McEwan AG, Sweet MJ. Metal ions in macrophage antimicrobial pathways: emerging roles for zinc and copper. Biosci Rep. 2013;33(4). doi: 10.1042/BSR20130014.

  110. Sun Y, Zhang Y, Fu X, Zhang R, Zou J, Wang S, et al. Identification of two secreted ferritin subunits involved in immune defense of Yesso scallop Patinopecten yessoensis. Fish Shellfish Immunol. 2014;37(1):53–9.

    Article  PubMed  CAS  Google Scholar 

  111. Torreilles J, Guerin MC, Roch P. Reactive oxygen species and defense mechanisms in marine bivalves. Cr Acad Sci Iii-Vie. 1996;319(3):209–18.

    CAS  Google Scholar 

  112. Sun Y, Li F, Sun Z, Zhang X, Li S, Zhang C, et al. Transcriptome Analysis of the Initial Stage of Acute WSSV Infection Caused by Temperature Change. PLoS One. 2014;9(3):e90732.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  113. Pauletto M, Milan M, Moreira R, Novoa B, Figueras A, Babbucci M, et al. Deep transcriptome sequencing of Pecten maximus hemocytes: a genomic resource for bivalve immunology. Fish Shellfish Immunol. 2014;37(1):154–65.

    Article  PubMed  CAS  Google Scholar 

  114. Moya A, Huisman L, Ball EE, Hayward DC, Grasso LC, Chua CM, et al. Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to CO(2)-driven acidification during the initiation of calcification. Mol Ecol. 2012;21(10):2440–54.

    Article  PubMed  CAS  Google Scholar 

  115. Lourenco AP, Martins JR, Bitondi MM, Simoes ZL. Trade-off between immune stimulation and expression of storage protein genes. Arch Insect Biochem Physiol. 2009;71(2):70–87.

    Article  PubMed  CAS  Google Scholar 

  116. Rauw WM. Immune response from a resource allocation perspective. Front Genet. 2012;3:267.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  117. Hegemann A, Matson KD, Versteegh MA, Tieleman BI. Wild Skylarks Seasonally Modulate Energy Budgets but Maintain Energetically Costly Inflammatory Immune Responses throughout the Annual Cycle. PLoS ONE. 2012;7(5):e36358.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  118. Simmons LW. Resource allocation trade-off between sperm quality and immunity in the field cricket, Teleogryllus oceanicus. Behav Ecol. 2011;23(1):168–73.

    Article  Google Scholar 

  119. Fabbri E, Valbonesi P, Franzellitti S. HSP expression in bivalves. Invertebrate Survival Journal: ISJ. 2008;5:135–61.

    Google Scholar 

  120. Moreira R, Balseiro P, Romero A, Dios S, Posada D, Novoa B, et al. Gene expression analysis of clams Ruditapes philippinarum and Ruditapes decussatus following bacterial infection yields molecular insights into pathogen resistance and immunity. Dev Comp Immunol. 2012;36(1):140–9.

    Article  PubMed  CAS  Google Scholar 

  121. Genard B, Miner P, Nicolas J-L, Moraga D, Boudry P, Pernet F, et al. Integrative Study of Physiological Changes Associated with Bacterial Infection in Pacific Oyster Larvae. PLoS ONE. 2013;8(5):e64534.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  122. Zugel U, Kaufmann SH. Role of heat shock proteins in protection from and pathogenesis of infectious diseases. Clin Microbiol Rev. 1999;12(1):19–39.

    PubMed Central  PubMed  CAS  Google Scholar 

  123. Osterloh A, Breloer M. Heat shock proteins: linking danger and pathogen recognition. Med Microbiol Immunol. 2008;197(1):1–8.

    Article  PubMed  CAS  Google Scholar 

  124. Ouzounis CA, Khatri P, Sirota M, Butte AJ. Ten Years of Pathway Analysis: Current Approaches and Outstanding Challenges. Plos Comput Biol. 2012;8(2):e1002375.

    Article  CAS  Google Scholar 

  125. Kamiya H, Muramoto K, Yamazaki M. Aplysianin-A, an antibacterial and antineoplastic glycoprotein in the albumen gland of a sea hare, Aplysia kurodai. Experientia. 1986;42(9):1065–7.

    Article  PubMed  CAS  Google Scholar 

  126. Iijima R, Kisugi J, Yamazaki M. Antifungal activity of Aplysianin E, a cytotoxic protein of sea hare (Aplysia kurodai) eggs. Dev Comp Immunol. 1995;19(1):13–9.

    Article  PubMed  CAS  Google Scholar 

  127. Oliver JL, Gaffney PM, Allen SK, Faisal M, Kaattari SL. Protease Inhibitory Activity in Selectively Bred Families of Eastern Oysters. J Aquat Anim Health. 2000;12(2):136–45.

    Article  Google Scholar 

  128. Chernikov OV, Molchanova VI, Chikalovets IV, Kondrashina AS, Li W, Lukyanov PA. Lectins of marine hydrobionts. Biochemistry (Moscow). 2013;78(7):760–70.

    Article  CAS  Google Scholar 

  129. Smolowitz R, Defaveri J, Walton W, Murphy D, Leavitt D. The interaction of temperature and hard clam (Mercenaria mercenaria) strain on the occurrence of QPX disease in the laboratory and in the field. J Shellfish Res. 2008;27(4):1052–3.

    Google Scholar 

  130. Dahl SF, Perrigault M, Allam B. Laboratory transmission studies of QPX disease in the hard clam: Interactions between different host strains and pathogen isolates. Aquaculture. 2008;280(1–4):64–70.

    Article  Google Scholar 

  131. Monari M, Foschi J, Rosmini R, Marin MG, Serrazanetti GP. Heat shock protein 70 response to physical and chemical stress in Chamelea gallina. J Exp Mar Biol Ecol. 2011;397(2):71–8.

    Article  CAS  Google Scholar 

  132. Metzger DC, Pratt P, Roberts SB. Characterizing the Effects of Heavy Metal and Vibrio Exposure on Hsp70 Expression in Crassostrea Gigas Gill Tissue. J Shellfish Res. 2012;31(3):627–30.

    Article  Google Scholar 

  133. Roberts RJ, Agius C, Saliba C, Bossier P, Sung YY. Heat shock proteins (chaperones) in fish and shellfish and their potential role in relation to fish health: a review. J Fish Dis. 2010;33(10):789–801.

    Article  PubMed  CAS  Google Scholar 

  134. Narusaka Y, Narusaka M, Seki M, Umezawa T, Ishida J, Nakajima M, et al. Crosstalk in the responses to abiotic and biotic stresses in Arabidopsis: Analysis of gene expression in cytochrome P450 gene superfamily by cDNA microarray. Plant Mol Biol. 2004;55(3):327–42.

    Article  PubMed  CAS  Google Scholar 

  135. Cardona P-J, Seifart Gomes C, Izar B, Pazan F, Mohamed W, Mraheil MA, et al. Universal Stress Proteins Are Important for Oxidative and Acid Stress Resistance and Growth of Listeria monocytogenes EGD-e In Vitro and In Vivo. PLoS ONE. 2011;6(9):e24965.

    Article  CAS  Google Scholar 

  136. Schmieder R, Lim YW, Edwards R. Identification and removal of ribosomal RNA sequences from metatranscriptomes. Bioinformatics. 2012;28(3):433–5.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  137. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  138. Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35(Web Server issue):W182–5.

    Article  PubMed Central  PubMed  Google Scholar 

  139. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

Download references

Acknowledgments

This research was supported by projects R/XG-19 and R/FBM-34 (to BA and EPE), funded under award NA07OAR4170010 from the National Sea Grant College Program of NOAA to the Research Foundation of State University of New York on behalf of New York Sea Grant. The study was also partially supported by project IOS 1050596 (to BA and EPE) supported by the National Science Foundation. The statements, findings, conclusions, views and recommendations are those of the authors and do not necessarily reflect the views of any of those organizations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bassem Allam.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BA and EPE designed the study, coordinated experiments and supervised data analysis and manuscript drafting. KW carried out the experiments, analyzed the data and drafted the manuscript. EC participated in data analysis and helped in drafting the methodology section of the manuscript. CDC helped with data analysis. All authors read and approved the final manuscript.

Additional files

Additional file 1:

Sequence summary of M. mercenaria de novo assembled transcriptome. This is a Microsoft Excel worksheet that contains descriptions and annotation information of individual transcript sequences. (XLSX 6600 kb)

Additional file 2:

Differentially expressed (DE) sequences associated with M. mercenaria focal response against QPX. This file is a Microsoft Excel worksheet that contains separate lists of over- and under-expressed transcripts with their expression levels and fold change values given. (XLSX 143 kb)

Additional file 3:

Differentially expressed (DE) sequences associated with M. mercenaria systemic response against QPX. This file is a Microsoft Excel worksheet that contains separate lists of over- and under-expressed transcripts with their expression levels and fold change values given. (XLSX 327 kb)

Additional file 4:

Transcripts that were differentially expressed in naïve M. mercenaria as compared to QPX infected individuals. This file is a Microsoft Excel worksheet that contains separate lists of over- and under-expressed transcripts with their expression levels and fold change values given. (XLSX 128 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, K., del Castillo, C., Corre, E. et al. Clam focal and systemic immune responses to QPX infection revealed by RNA-seq technology. BMC Genomics 17, 146 (2016). https://doi.org/10.1186/s12864-016-2493-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-016-2493-9

Keywords