Toll signal transduction pathway in bivalves: complete cds of intermediate elements and related gene transcription levels in hemocytes of immune stimulated Mytilus galloprovincialis

Transcriptomic response of Mytilus coruscus mantle to acute sea water acidification and shell damage

Mytilus coruscus is an economically important marine calcifier living in the Yangtze River estuary sea area, where seasonal fluctuations in natural pH occur owing to freshwater input, resulting in a rapid reduction in seawater pH. In addition, Mytilus constantly suffers from shell fracture or injury in the natural environment, and the shell repair mechanisms in mussels have evolved to counteract shell injury. Therefore, we utilized shell-complete and shell-damaged Mytilus coruscus in this study and performed transcriptomic analysis of the mantle to investigate whether the expression of mantle-specific genes can be induced by acute seawater acidification and how the mantle responds to acute acidification during the shell repair process. We found that acute acidification induced more differentially expressed genes than shell damage in the mantle, and the biomineralization-related Gene Ontology terms and KEGG pathways were significantly enriched by these DEGs. Most DEGs were upregulated in enriched pathways, indicating the activation of biomineralization-related processes in the mussel mantle under acute acidification. The expression levels of some shell matrix proteins and antimicrobial peptides increased under acute acidification and/or shell damage, suggesting the molecular modulation of the mantle for the preparation and activation of the shell repairing and anti-infection under adverse environmental conditions. In addition, morphological and microstructural analyses were performed for the mantle edge and shell cross-section, and changes in the mantle secretory capacity and shell inner film system induced by the two stressors were observed. Our findings highlight the adaptation of M. coruscus in estuarine areas with dramatic fluctuations in pH and may prove instrumental in its ability to survive ocean acidification.

TLR4 involved in immune response against Vibrio Parahaemolyticus by MyD88-dependent pathway in Crassostrea hongkongensis

Vibrio parahaemolyticus (V. parahaemolyticus) is a salt-loving gram-negative bacterium, and is the leading cause of mortality in cultured shellfish in recent years. Toll-like Receptor 4 (TLR4) is a classical pattern recognition receptor (PRRs) that recognizes pathogen-associated molecular patterns (PAMPs) of pathogenic microorganism and activates the immune response. However, the function and signal pathway of TLR4 in oyster are still unknown. In this study, a new TLR4 gene was identified from the Crassostrea hongkongensis (C. hongkongensis). The ChTLR4 contained an open reading frame of 2643 bp, encoding 880 amino acids with seven leucine-rich repeat (LRR) domains and a Toll/IL-1R (TIR) domain. The ChTLR4 shared the highest sequence identity (83.0%) with TLR4 of Crassostrea gigas. Tissue expression analysis revealed that ChTLR4 showed the highest constitutive expression in the gill and hepatopancreas, and was significantly upregulated in immune tissues post V. parahaemolyticus infection, especially in gill and hemocytes. Moreover, TLR4 silencing significantly inhibited the immune-enzyme activities, including SOD, CAT, ACP, AKP in gill and LZM in hemolymph supernatant, and increased MDA content in hemolymph supernatant. Meanwhile, the antimicrobial activities of the hemolymph supernatant were also significantly inhibited by TLR4 silencing. These data demonstrated that the ChTLR4 involved in innate immune response of C. hongkongensis against V. parahaemolyticus challenge. Finally, qRT-PCR analysis showed that ChTLR4 silencing clearly inhibited the expression of genes in TLR4-MyD88 pathway, indicating that MyD88-dependent pathway played a crucial role in ChTLR4-mediated immune response against V. parahaemolyticus.

Comparative transcriptomic analysis revealed changes in multiple signaling pathways involved in protein degradation in the digestive gland of Mytilus coruscus during high-temperatures

As a result of global warming, the Mytilus coruscus living attached in the intertidal zone experience extreme and fluctuating changes in temperature, and extreme temperature changes are causing mass mortality of intertidal species. This study explores the transcriptional response of M. coruscus at different temperatures (18 °C, 26 °C, and 33 °C) and different times (0, 12, and 24 h) of action by analyzing the potential temperature of the intertidal zone. In response to high temperatures, several signaling pathways in M. coruscus, ribosome, endocytosis, endoplasmic reticulum stress, protein degradation, and lysosomes, interact to counter the adverse effects of high temperatures on protein homeostasis. Increased expression of key genes, including heat shock proteins (Hsp70, Hsp20, and Hsp110), Lysosome-associated membrane glycoprotein (LAMP), endoplasmic reticulum chaperone (BiP), and baculoviral IAP repeat-containing protein 7 (BIRC7), may further mitigate the effects of heat stress and delay mortality in M. coruscus. These results reveal changes in multiple signaling pathways involved in protein degradation during high-temperature stress, which will contribute to our overall understanding of the molecular mechanisms underlying the response of M. coruscus to high-temperature stress.

Understanding Crassostrea virginica tolerance of Perkinsus marinus through global gene expression analysis

Disease tolerance, a host's ability to limit damage from a given parasite burden, is quantified by the relationship between pathogen load and host survival or reproduction. Dermo disease, caused by the protozoan parasite P. marinus, negatively impacts survival in both wild and cultured eastern oyster (C. virginica) populations. Resistance to P. marinus has been the focus of previous studies, but tolerance also has important consequences for disease management in cultured and wild populations. In this study we measured dermo tolerance and evaluated global expression patterns of two sensitive and two tolerant eastern oyster families experimentally challenged with distinct doses of P. marinus (0, 10⁶, 10⁷, and 10⁸ parasite spores per gram wet weight, n = 3-5 individuals per family per dose). Weighted Gene Correlation Network Analysis (WGCNA) identified several modules correlated with increasing parasite dose/infection intensity, as well as phenotype. Modules positively correlated with dose included transcripts and enriched GO terms related to hemocyte activation and cell cycle activity. Additionally, these modules included G-protein coupled receptor, toll-like receptor, and tumor necrosis factor pathways, which are important for immune effector molecule and apoptosis activation. Increased metabolic activity was also positively correlated with treatment. The module negatively correlated with infection intensity was enriched with GO terms associated with normal cellular activity and growth, indicating a trade-off with increased immune response. The module positively correlated with the tolerant phenotype was enriched for transcripts associated with "programmed cell death" and contained a large number of tripartite motif-containing proteins. Differential expression analysis was also performed on the 10⁸ dosed group using the most sensitive family as the comparison reference. Results were consistent with the network analysis, but signals for "programmed cell death" and serine protease inhibitors were stronger in one tolerant family than the other, suggesting that there are multiple avenues for disease tolerance. These results provide new insight for defining dermo response traits and have important implications for applying selective breeding for disease management.

The localization of Toll and Imd pathway and complement system components and their response to Vibrio infection in the nemertean Lineus ruber

BACKGROUND

Innate immunity is the first line of defense against pathogens. In animals, the Toll pathway, the Imd pathway, the complement system, and lectins are well-known mechanisms involved in innate immunity. Although these pathways and systems are well understood in vertebrates and arthropods, they are understudied in other invertebrates.

RESULTS

To shed light on immunity in the nemertean Lineus ruber, we performed a transcriptomic survey and identified the main components of the Toll pathway (e.g., myD88, dorsal/dif/NFκB-p65), the Imd pathway (e.g., imd, relish/NFκB-p105/100), the complement system (e.g., C3, cfb), and some lectins (FreD-Cs and C-lectins). In situ hybridization showed that TLRβ1, TLRβ2, and imd are expressed in the nervous system; the complement gene C3-1 is expressed in the gut; and the lectins are expressed in the nervous system, the blood, and the gut. To reveal their potential role in defense mechanisms, we performed immune challenge experiments, in which Lineus ruber specimens were exposed to the gram-negative bacteria Vibrio diazotrophicus. Our results show the upregulation of specific components of the Toll pathway (TLRα3, TLRβ1, and TLRβ2), the complement system (C3-1), and lectins (c-lectin2 and fred-c5).

CONCLUSIONS

Therefore, similarly to what occurs in other invertebrates, our study shows that components of the Toll pathway, the complement system, and lectins are involved in the immune response in the nemertean Lineus ruber. The presence of these pathways and systems in Lineus ruber, but also in other spiralians; in ecdysozoans; and in deuterostomes suggests that these pathways and systems were involved in the immune response in the stem species of Bilateria.

Transcriptomic profiling of immune-associated molecules in the coelomocytes of lugworm Arenicola marina (Linnaeus, 1758)

Organization and functioning of immune system remain unevenly studied in different taxa of lophotrochozoan animals. We analyzed transcriptomic data on coelomocytes of the lugworm Arenicola marina (Linnaeus, 1758; Annelida, Polychaeta) to gain insights into the molecular mechanisms involved in polychaete immunity. Coelomocytes are specialized motile cells populating coelomic fluid of annelids, responsible for cellular defense reactions and providing humoral immune factors. The transcriptome was enriched with immune-related transcripts by challenging the cells in vitro with lipopolysaccharides of Escherichia coli and Zymosan from Saccharomyces cerevisiae. Our analysis revealed a multifaceted and complex internal defense system of the lugworm. A. marina possesses orthologs of proto-complement-like factors: six thioester-containing proteins, a complement-like receptor, and a MASP-related serine protease (MReM2). A. marina coelomocytes employ pattern-recognition receptors to detect pathogens and regulate immune responses. Among them, there are 18 Toll-like receptors and various putative lectin-like proteins with evolutionary conserved and taxa-specific domains. C-type lectins and a novel family of Gal-binding and CUB domains containing receptors were the most abundant in the transcriptome. The array of pore-forming proteins in the coelomocytes was surprisingly reduced compared to that of other invertebrate species. We characterized a set of conserved proteins metabolizing reactive oxygen species and nitric oxide and expanded the arsenal of potential antimicrobial peptides. Phenoloxidase activity in immune cells of lugworm is mediated only by laccase enzyme. The described repertoire of immune-associated molecules provides valuable candidates for further functional and comparative research on the immunity of annelids.

AIF-1 and RNASET2 are involved in the inflammatory response in the Mediterranean mussel Mytilus galloprovincialis following Vibrio infection

Filter-feeding bivalves, such as the Mytilus species, are exposed to different types of bacteria in the surrounding waters, in particular of the Vibrio genus. Mussels lack an adaptive immune system and hemocytes can recognize pathogen-associated molecular patterns (PAMPs) via pattern recognition receptors (PRRs) to activate intracellular signaling pathways to trigger the antimicrobial effectors synthesis. Among the areas of bivalve immunity that deserve study include the role of hemocyte subpopulations. Since little information are available on immune responses at the tissue level to human pathogenic vibrios commonly detected in coastal waters involved in seafood-borne diseases, in this work, immunological parameters of the hemocytes from the Mediterranean mussel M. galloprovincialis were evaluated in response to in vivo challenge with Vibrio splendidus. The histological approach has been first used in order to identify the hemocytes recruitment at the infection site and the morphological change of muscular fibers. In addition, using immunolabeling with specific antibody we detected the production of molecules involved in the inflammatory activated cascade: Toll-like receptors 4 (TLR4), the myeloid differentiation factor 88 (MyD88), the Allograft inflammatory factor-1 (AIF-1) and the ribonucleases RNASET2, belonging to the T2 family, that in vertebrates are involved in the recruitment and activation of macrophages. Our results indicate the activation of TLR4 during bacterial infection preparatory to the recruitment of the MyD88 adapter with a putative role in recognition and intracellular signalling. Furthermore, the data presented in this work suggest that challenging with Gram-negative bacteria causes a massive migration of AIF-1⁺ hemocytes and that the ribonuclease RNASET2 could play a key role in the recruitment of these activated hemocytes. Our approach is useful for further understanding the complex molecular defence mechanisms of the host in invertebrates, especially in relation to the need to develop methods to evaluate the immunological response of bivalve molluscs used in aquaculture.

The Evolution and Diversity of Interleukin-17 Highlight an Expansion in Marine Invertebrates and Its Conserved Role in Mucosal Immunity

The interleukin-17 (IL-17) family consists of proinflammatory cytokines conserved during evolution. A comparative genomics approach was applied to examine IL-17 throughout evolution from poriferans to higher vertebrates. Cnidaria was highlighted as the most ancient diverged phylum, and several evolutionary patterns were revealed. Large expansions of the IL-17 repertoire were observed in marine molluscs and echinoderm species. We further studied this expansion in filter-fed Mytilus galloprovincialis, which is a bivalve with a highly effective innate immune system supported by a variable pangenome. We recovered 379 unique IL-17 sequences and 96 receptors from individual genomes that were classified into 23 and 6 isoforms after phylogenetic analyses. Mussel IL-17 isoforms were conserved among individuals and shared between closely related Mytilidae species. Certain isoforms were specifically implicated in the response to a waterborne infection with Vibrio splendidus in mussel gills. The involvement of IL-17 in mucosal immune responses could be conserved in higher vertebrates from these ancestral lineages.

Silver nanoparticles' impact on the gene expression of the cytosolic adaptor MyD-88 and the interferon regulatory factor IRF in the gills and digestive gland of mytilus galloprovincialis

Silver nanoparticles (AgNPs) have been reported as stressors for the bivalves' immune system at different regulatory levels, impacting the detection step and receptors, and other mediators, as well as effector molecules. However, studies on how AgNPs impact the transmission of signals from receptors and whether they have an effect on mediators and transcription factors are still scarce. This study aims to investigate the effect of 12 hours of in vivo exposure to 100 µg/L of AgNPs on the gene expression of the cytosolic adaptor Myeloid, the differentiation protein 88 (MgMyD88-b), and the interferon regulatory factor (Me4-IRF) in the gills and digestive gland of Mytilus galloprovincialis, before and after blocking two major uptake pathways of nanoparticles (clathrin- and caveolae-mediated endocytosis). The results illustrate a tissue-specific gene expression of the MgMyD88-b and the Me4-IRF in the gills and digestive gland of M. galloprovincialis. In the gills, AgNPs did not significantly impact the expression of the two genes. However, blocking the caveolae-mediated endocytosis decreased the expression of Me4-IRF. However, inhibition of clathrin-mediated endocytosis in the digestive gland recorded a significant decrease in the expression of MgMyD88-b. Overall, the inhibition of the AgNPs' uptake routes have highlighted their potential interference with the immune response through the studied mediators' genes, which need to be studied further in future investigations.

Transcriptomic Response of Mussel Gills After a Vibrio splendidus Infection Demonstrates Their Role in the Immune Response

Mussels (Mytilus galloprovincialis) are filter feeder bivalves that are constantly in contact with a wide range of microorganisms, some of which are potentially pathogenic. How mussels recognize and respond to pathogens has not been fully elucidated to date; therefore, we investigated the immune mechanisms that these animals employ in response to a bacterial bath infection from the surrounding water, mimicking the response that mussels mount under natural conditions. After the bath infection, mussels were able to remove the bacteria from their bodies and from the water tank. Accordingly, antibacterial activity was detected in gill extracts, demonstrating that this tissue plays a central role in removing and clearing potential pathogens. A transcriptomic study performed after a bath infection with Vibrio splendidus identified a total of 1,156 differentially expressed genes. The expression levels of genes contributing to a number of biological processes, such as immune response activation pathways and their regulation with cytokines, cell recognition, adhesion and apoptosis, were significantly modulated after infection, suggesting that the gills play important roles in pathogen recognition, as well as being activators and regulators of the mussel innate immune response. In addition to RNA-seq analysis, long non-coding RNAs and their neighboring genes were also analyzed and exhibited modulation after the bacterial challenge. The response of gills against bath infection was compared with the findings of a previous transcriptomic study on hemocytes responding to systemic infection, demonstrating the different and specific functions of gills. The results of this study indicate that recognition processes occur in the gill, thereby activating the effector agents of the immune response to overcome bacterial infection.

Toll-like signaling pathway in the transcriptome of Littorina littorea

The transcriptome of the caenogastropod mollusk Littorina littorea was scanned for the presence of sequences encoding Toll-like receptors (TLRs) and corresponding proteins involved in downstream TLR signaling pathway. In the transcriptomic snapshots of hemocytes and kidney tissues, 45 complete TLRs encoded by 35 genes were identified. Out of the 59 non-TLR molecules involved in a canonical TLR signaling pathway, 35 genes were classified as homologous and could be placed within the TLR-mediated MyD88-and MAPK-dependent circuitries. No reference vertebrate adapters TIRAP, TRIF and TRAM were identified in the transcriptome. The results of RNA-seq experiments with an immune challenge (rediae of the digenean Himasthla elongata) indicate that four TLRs (LlTLR1, 3, 5 and 8) and a set of upregulated genes involved in signal transduction (LlMyd88, LlTNFα, LlCASP8, LlFADD, LlNFKBIA (IkBα), LlIRAK1, LlSTAT1, LlMAPK14 (P38), LlMAP2K1 (MEK1/2), LlIRF3 and LlIRF5) may participate in the anti-digenean immune response of L. littorea.

Contrasting Immunomodulatory Effects of Probiotic and Pathogenic Bacteria on Eastern Oyster, Crassostrea Virginica, Larvae

Several Vibrio spp. cause acute and severe mortality events in hatcheries where larvae of bivalve mollusks are reared, potentially leading to subsequent shortage of bivalve seed for the grow-out industry. In particular, strains of Vibrio coralliilyticus have been identified as a major cause of disease in Pacific, Crassostrea gigas, and eastern, C. virginica, oyster hatcheries in the United States of America. Probiotic bacteria are an inexpensive, practical, and natural method of disease control. Previous research shows that pretreatment of larval oysters with probiotic bacteria Bacillus pumilus RI06-95 (RI) and Phaeobacter inhibens S4 (S4) significantly decreases mortality caused by experimental challenge with the bacterial pathogen V. coralliilyticus RE22 (RE22). This study aims to characterize the immune response of 6-10-day-old eastern oyster larvae to experimental challenge with pathogen V. coralliilyticus RE22 and probionts RI and S4. Treatments included (a) pathogen and probiont exposure at a concentration of 5 × 10⁴ CFU per mL (~2500 bacterial cells per larva) for a duration of 6 h, (b) probiont exposure at the same concentration for a duration of 24 h, and (c) probiont RI daily treatment of larvae in the hatchery for 4, 11, and 15 days. Differential gene expression analysis compared pathogen or probiotic-treated transcriptomes to unexposed controls. Probiotic and pathogen treatment led to upregulation of transcripts coding for several immune pattern recognition receptors (PRRs) involved in environmental sensing and detection of microbes in oyster larvae. Larval oyster responses to pathogen RE22 suggested suppression of expression of genes in immune signaling pathways (myd88, tak1, nkap), failure in upregulation of immune effector genes, high metabolic demand, and oxidative stress that potentially contributed to mortality. On the other hand, the transcriptomic response to probiotic bacteria RI and S4 suggested activation of immune signaling pathways and expression of immune effectors (e.g., Cv-spi2, mucins and perforin-2). These key features of the host immune response to probiotic bacteria were shared despite the length of probiotic exposure, probiotic species, and the type of environment in which exposures were conducted. This study suggests that pre-exposure of eastern oyster larvae to probiotics for 6-24 h prior to pathogenic challenge leads to a robust and effective immune response that may contribute to protecting larvae from subsequent challenge with V. coralliilyticus RE22. This research provides new insights into host-microbe interactions in larval oysters that could be applied in the management of vibriosis in bivalve hatcheries.

In-Depth In Silico Search for Cuttlefish (Sepia officinalis) Antimicrobial Peptides Following Bacterial Challenge of Haemocytes

Cuttlefish (Sepia officinalis) haemocytes are potential sources of antimicrobial peptides (AMPs). To study the immune response to Vibrio splendidus and identify new AMPs, an original approach was developed based on a differential transcriptomic study and an in-depth in silico analysis using multiple tools. Two de novo transcriptomes were retrieved from cuttlefish haemocytes following challenge by V. splendidus or not. A first analysis of the annotated transcripts revealed the presence of Toll/NF-κB pathway members, including newly identified factors such as So-TLR-h, So-IKK-h and So-Rel/NF-κB-h. Out of the eight Toll/NF-κB pathway members, seven were found up-regulated following V. splendidus challenge. Besides, immune factors involved in the immune response were also identified and up-regulated. However, no AMP was identified based on annotation or conserved pattern searches. We therefore performed an in-depth in silico analysis of unannotated transcripts based on differential expression and sequence characteristics, using several tools available like PepTraq, a homemade software program. Finally, five AMP candidates were synthesized. Among them, NF19, AV19 and GK28 displayed antibacterial activity against Gram-negative bacteria. Each peptide had a different spectrum of activity, notably against Vibrio species. GK28—the most active peptide—was not haemolytic, whereas NF19 and AV19 were haemolytic at concentrations between 50 and 100 µM, 5 to 10 times higher than their minimum inhibitory concentration.

Molecular features of interleukin-1 receptor-associated kinase-b and a in Mytilus coruscus, regulating their function by infection of aquatic pathogens and the expression of their serine/threonine protein kinase functional domains

Interleukin-1 receptor-associated kinases (IRAKs) play important roles in the innate immune system of TLR (Toll-like receptor) signaling pathway. In this paper, interleukin-1 receptor-associated kinase-b (designated as McIRAK-b) and interleukin-1 receptor-associated kinase-a (named as McIRAK-a) were obtained based on the transcriptome data, the full length of McIRAK-b was 1815 bp and McIRAK-a was 3168bp, encoding 532 and 978 amino acids, respectively. BLASTp analysis and phylogenetic relationship strongly suggested that the deduced amino acid sequence of McIRAK-b had high homology with IRAK-4 and McIRAK-a was similar to IRAK-1 of other mollusks, especially at their function domains. The expressions of McIRAK-b and McIRAK-a were detected in six tissues including adductor muscle, hemocyte, gills, gonad and hepatopancreas, and the highest expressions appeared both in gills. The expressions of McIRAK-b and McIRAK-a in gills were observed with time-dependent manners after bacterial infections. After being challenged with Vibrio alginolyticus, McIRAK-b expressed significantly and got the peak at 8 h (9.47 times compared with the control group), but the peak appeared at 4 h by being infected with Vibrio parahaemolyticus (12.02 times compared with the control group). The highest point of McIRAK-a mRNA appeared at 12 h (5.16 times) after being challenged with V.alginolyticus and 8 h (4.21 times) for V.parahaemolyticus challenge. The results suggested that IRAK-b and IRAK-a might be important in immune signaling pathway of mussels. The kinase functional domain sequences (S_TKc) of McIRAK-b and McIRAK-a expressed in BL21(DE3) and purified by Ni-NAT Superflow resin conforming to the expected molecular weight with many active sites for their conferring protein-protein interaction functions. This study may provide some further understandings of the regulatory mechanisms in the bivalve innate immune system for IRAKs family.

Two toll-like receptors identified in the mantle of Mytilus coruscus are abundant in haemocytes

Toll-like receptors (TLRs) are a large family of pattern recognition receptors (PRRs) that play a critical role in innate immunity. TLRs are activated when they recognize microbial associated molecular patterns (MAMPs) of bacteria, viruses, or fungus. In the present study, two TLRs were isolated from the mantle of the hard-shelled mussel (Mytilus coruscus) and designated McTLR2 and McTLR3 based on their sequence similarity and phylogenetic clustering with Crassostrea gigas, CgiTLR2 and CgiTLR3, respectively. Quantitative RT-PCR analysis demonstrated that McTLR2 and McTLR3 were constitutively expressed in many tissues but at low abundance.

Immunity in mussels: An overview of molecular components and mechanisms with a focus on the functional defenses

Bivalves' immunity has received much more attention in the last decade, which resulted to a valuable growth in the availability of its molecular components. Such data availability coupled with the economical importance of these organisms aimed to shift the increase in the number of immunological and stress-related studies. Unfortunately, the crowd of generated data deciphering the involved physiological processes, investigators' differential conceptualization and the aimed objectives, has complicated the sensu stricto outlining of immune-related mechanisms. Overall, this review tried to compiles a summary about the molecular components of the mussels' immune response, surveying an overview of the mussels' functional immunity through gathering the most recent-related topics of bivalves' immunity as apoptosis and autophagy which deserves a great attention as stress-related mechanisms, the disseminated neoplasia as outbreak transmissible disease, not only within the same specie but also among different species, the hematopoiesis as topic that still generating interesting debate in the scientific community, the mucosal immunity described as the interface where host-pathogen interactions would occurs and determinate the late immune response, and innate immune memory and transgenerational priming, which described as very recent research topic with extensive applications in shellfish farming industry.

Rhodnius prolixus: Identification of missing components of the IMD immune signaling pathway and functional characterization of its role in eliminating bacteria

The innate immune system in insects is regulated by specific signalling pathways. Most immune related pathways were identified and characterized in holometabolous insects such as Drosophila melanogaster, and it was assumed they would be highly conserved in all insects. The hemimetabolous insect, Rhodnius prolixus, has served as a model to study basic insect physiology, but also is a major vector of the human parasite, Trypanosoma cruzi, that causes 10,000 deaths annually. The publication of the R. prolixus genome revealed that one of the main immune pathways, the Immune-deficiency pathway (IMD), was incomplete and probably non-functional, an observation shared with other hemimetabolous insects including the pea aphid (Acyrthosiphon pisum) and the bedbug (Cimex lectularius). It was proposed that the IMD pathway is inactive in R. prolixus as an adaptation to prevent eliminating beneficial symbiont gut bacteria. We used bioinformatic analyses based on reciprocal BLAST and HMM-profile searches to find orthologs for most of the "missing" elements of the IMD pathway and provide data that these are regulated in response to infection with Gram-negative bacteria. We used RNAi strategies to demonstrate the role of the IMD pathway in regulating the expression of specific antimicrobial peptides (AMPs) in the fat body of R. prolixus. The data indicate that the IMD pathway is present and active in R. prolixus, which opens up new avenues of research on R. prolixus-T. cruzi interactions.

The GRP94 gene of Yesso scallop (Patinopecten yessoensis): Characterization and expression regulation in response to thermal and bacterial stresses

The 94-kDa glucose-regulated protein (GRP94) belonging to the HSP90 family is an endoplasmic reticulum (ER) chaperone. It plays critical roles in ER quality control, and has been implicated as a specialized immune chaperone to regulate both innate and adaptive immunity. In this study, we identified and characterized a GRP94 gene (PyGRP94) from Yesso scallop (Patinopecten yessoensis). The protein sequence of PyGRP94 is highly conserved with its homologs in vertebrates, with a signal sequence in N-terminal, an ER retrieval signal sequence in C-terminal and a HATPase_c domain. Expression analysis suggests that PyGRP94 transcripts in early embryos are maternally derived and the zygotic expression is started from D-shaped larvae. This gene is also expressed in almost all the adult tissues examined except smooth muscle, with the highest expression level in hemocytes. Besides, PyGRP94 was demonstrated to be induced by heat shock and both Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacterial infection, with much more dramatic changes being observed after V. anguillarum challenge. Our results suggest the involvement of PyGRP94 in response to thermal stress, and that it might play an important role in the innate immune defense of scallop.

A novel toll-like receptor from Mytilus coruscus is induced in response to stress

Toll-like receptor (TLR) is considered to be an evolutionarily conserved transmembrane protein which promotes the Toll signal pathway to active the expression of transcription factors in the innate immunity of the organism. In this study, a full length of TLR homologue of 2525bp in Mytilus coruscus (named as McTLR-a, GenBank accession no: KY940571) was characterized. Its ORF was 1815 bp with a 5'untranslated region (UTR) of 128 bp and a 3'UTR of 582 bp, encoding 602 amino acid residues with a calculated molecular weight of 70.870 kDa (pI = 6.10). BLASTn analysis and phylogenetic relationship strongly suggested that this cDNA sequence was a member of TLR family. Quantitative real time RT-PCR showed that constitutive expression of McTLR-a was occurred, with increasing order in hemocyte, gonad, mantle, adducter, gill and hepatopancreas. Bacterial infection and heavy metals stimulation up-regulated the expression of McTLR-a mRNA in hepatopancreas with time-dependent manners. The maximum expression appeared at 12 h after pathogenic bacteria injection, with approximately 22-fold in Aeromonas hydrophila and 17-fold in Vibrio parahemolyticus higher than that of the blank group. In heavy metals stress group, they all reached peaks at 3d, while the diverse concentration caused the maximum expression were different. The highest expression reached approximately 7-fold higher than the blank in low concentration of Pb²⁺ exposure. In Cu²⁺ treated group, it reached the peak (approximately 12-fold higher than the blank)in middle concentration. These results indicated that McTLR-a might be involved in the defense response and had a significant role in mediating the environmental stress.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

Transcriptome analysis of Mytilus coruscus hemocytes in response to Vibrio alginnolyficus infection

As an economically important bivalve, the Mytilus coruscus is cultured widely in the eastern coast of China. In recent years, this bivalve has been seriously affected by the pathogenic infections. To elucidate the host defense mechanisms of M. coruscus against pathogenic challenge, the hemocyte transcriptomes of M. coruscus before and after Vibrio alginnolyficus infection were analyzed using the deep-sequencing platform Illumina/HiSeq-2500, meanwhile the differentially expressed genes (DEGs) were investigated. In total, 130,031,083 clean reads were obtained and then assembled into 63,942 unigenes with an average length of 810 bp and an N50 of 1056 bp. Unigenes were annotated by comparing against nr, Swiss-Prot, KEGG, COG, KOG, GO, and Pfam databases, and 27,345 unigenes (42.77%) were annotated in at least one database. After bacterial challenge, 1270 and 265 genes were identified as remarkably up- or down-regulated, respectively, amongst 1154 were associated with 122 pathways, including classical immune-related pathways, such as 'Toll-like receptor signaling', 'the complement cascades', 'MAPK signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. Besides, nine genes which were differently-expressed immuno-related were confirmed by using quantitative real-time PCR. These findings would provide new insights on the M. coruscus innate immunity, based on which, some novel strategies for management of diseases and long-term sustainability of M. coruscus culture could be developed.

Molecular identification of disk abalone (Haliotis discus discus) tetraspanin 33 and CD63: Insights into potent players in the disk abalone host defense system

Tetraspanins are a superfamily of transmembrane proteins involved in a diverse range of physiological processes including differentiation, adhesion, signal transduction, cell motility, and immune responses. In the present study, two tetraspanins, CD63 and tetraspanin 33 (TSPAN33) from disk abalone (AbCD63 and AbTSPAN33), were identified and characterized at the molecular level. The coding sequences for AbCD63 and AbTSPAN33 encoded polypeptides of 234 and 290 amino acids (aa) with predicted molecular mass of 25.3 and 32.5 kDa, respectively. The deduced AbCD63 and AbTSPAN33 protein sequences were also predicted to have a typical tetraspanin domain architecture, including four transmembrane domains (TM), short N- and C- terminal regions, a short intracellular loop, as well as a large and small extracellular loop. A characteristic CCG motif and cysteine residues, which are highly conserved across CD63 and TSPAN33 proteins of different species, were present in the large extracellular loop of both abalone tetraspanins. Phylogenetic analysis revealed that the AbCD63 and AbTSPAN33 clustered in the invertebrate subclade of tetraspanins, thus exhibiting a close relationship with tetraspanins of other mollusks. The AbCD63 and AbTSPAN33 mRNA transcripts were detected at early embryonic development stages of disk abalone with significantly higher amounts at the trochophore stage, suggesting the involvement of these proteins in embryonic development. Both AbCD63 and AbTSPAN33 were ubiquitously expressed in all the tissues of unchallenged abalones analyzed, with the highest expression levels found in hemocytes. Moreover, significant induction of AbCD63 and AbTSPAN33 mRNA expression was observed in immunologically important tissues, such as hemocytes and gills, upon stimulation with live bacteria (Vibrio parahaemolyticus and Listeria monocytogenes), virus (viral hemorrhagic septicemia virus), and two potent immune stimulators [polyinosinic:polycytidylic acid (poly I:C) and lipopolysaccharide (LPS)]. Collectively, these findings suggest that AbCD63 and AbTSPAN33 are involved in innate immune responses in disk abalone during pathogenic stress.

Microbial Diseases of Bivalve Mollusks: Infections, Immunology and Antimicrobial Defense

A variety of bivalve mollusks (phylum Mollusca, class Bivalvia) constitute a prominent commodity in fisheries and aquacultures, but are also crucial in order to preserve our ecosystem’s complexity and function. Bivalve mollusks, such as clams, mussels, oysters and scallops, are relevant bred species, and their global farming maintains a high incremental annual growth rate, representing a considerable proportion of the overall fishery activities. Bivalve mollusks are filter feeders; therefore by filtering a great quantity of water, they may bioaccumulate in their tissues a high number of microorganisms that can be considered infectious for humans and higher vertebrates. Moreover, since some pathogens are also able to infect bivalve mollusks, they are a threat for the entire mollusk farming industry. In consideration of the leading role in aquaculture and the growing financial importance of bivalve farming, much interest has been recently devoted to investigate the pathogenesis of infectious diseases of these mollusks in order to be prepared for public health emergencies and to avoid dreadful income losses. Several bacterial and viral pathogens will be described herein. Despite the minor complexity of the organization of the immune system of bivalves, compared to mammalian immune systems, a precise description of the different mechanisms that induce its activation and functioning is still missing. In the present review, a substantial consideration will be devoted in outlining the immune responses of bivalves and their repertoire of immune cells. Finally, we will focus on the description of antimicrobial peptides that have been identified and characterized in bivalve mollusks. Their structural and antimicrobial features are also of great interest for the biotechnology sector as antimicrobial templates to combat the increasing antibiotic-resistance of different pathogenic bacteria that plague the human population all over the world.

Diversity and evolution of TIR-domain-containing proteins in bivalves and Metazoa: New insights from comparative genomics

The Toll/interleukin-1 receptor (TIR) domain has a fundamental role in the innate defence response of plants, vertebrate and invertebrate animals. Mostly found in the cytosolic side of membrane-bound receptor proteins, it mediates the intracellular signalling upon pathogen recognition via heterotypic interactions. Although a number of TIR-domain-containing (TIR-DC) proteins have been characterized in vertebrates, their evolutionary relationships and functional role in protostomes are still largely unknown. Due to the high abundance and diversity of TIR-DC proteins in bivalve molluscs, we investigated this class of marine invertebrates as a case study. The analysis of the available genomic and transcriptomic data allowed the identification of over 400 full-length sequences and their classification in protein families based on sequence homology and domain organization. In addition to TLRs and MyD88 adaptors, bivalves possess a surprisingly large repertoire of intracellular TIR-DC proteins, which are conserved across a broad range of metazoan taxa. Overall, we report the expansion and diversification of TIR-DC proteins in several invertebrate lineages and the identification of many novel protein families possibly involved in both immune-related signalling and embryonic development.

Characterization of the TRAF3IP1 gene in Yesso scallop (Patinopecten yessoensis) and its expression in response to bacterial challenge

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is an important adaptor that transmits upstream activation signals to induce innate immune responses. TRAF3 interacting protein 1 (TRAF3IP1) interacts specifically with TRAF3, but its function in innate immunity remains unclear, especially in marine invertebrates. In this study, to better understand the functions of TRAFs in innate immune responses, we identified and characterized the first bivalve TRAF3IP1 gene, PyTRAF3IP1, from Yesso scallop (Patinopecten yessoensis), one of the most important mollusk species for aquaculture. The PyTRAF3IP1 cDNA is 2,367 bp, with an open reading frame of 1,629 bp encoding 542 amino acids. Phylogenetic and protein structural analysis confirmed the gene's identity and revealed that PyTRAF3IP1 was more similar to vertebrate TRAF3IP1s than to those of invertebrates. PyTRAF3IP1 was expressed in all the adult tissues and developmental stages sampled, implying that it plays versatile roles in many biological processes. Furthermore, PyTRAF3IP1 expression was dramatically induced in the acute phase (3-6 h) after infection with both Gram-positive (Micrococcus luteus) and Gram-negative (Vibrio anguillarum) bacteria, even stronger induction being observed after V. anguillarum challenge. This is the first report of the characterization and immune response involvement of TRAF3IP1 in marine invertebrates, and suggests that TRAF3IP1 contributes to innate immunity in bivalves.

Identification and characterization of two LBP/BPI genes involved in innate immunity from Hyriopsis cumingii

Lipopolysaccharide-binding protein and bactericidal permeability-increasing protein (LBP/BPI) play crucial role in modulating cellular signals in response to Gram-negative bacteria infection. In the present study, two isoforms of LBP/BPI genes, designated as HcLBP/BPI1 and HcLBP/BPI2, respectively, were cloned from the mussel Hyriopsis cumingii by RACE approach. The full-length cDNA sequences of HcLBP/BPI1 and HcLBP/BPI2 were 1887 and 2227 bp and encoded two secreted proteins of 501 and 518 amino acid residues, respectively. The deduced amino acid of HcLBP/BPI1 and HcLBP/BPI2 contained several conserved domains, such as signal peptide, two BPI/LBP and one central domain. Phylogentic analysis further supported that HcLBP/BPI1 and HcLBP/BPI2 belonged to new members of invertebrate LBP/BPI family. The mRNA transcripts of HcLBP/BPI1 and HcLBP/BPI2 were ubiquitously expressed in all examined tissues, and the expression level of HcLBP/BPI1 was higher than that of HcLBP/BPI2. The mRNA expression of HcLBP/BPI1 in hepatopancreas and hemocytes was significantly up-regulate after Aeromonas hydrophila and LPS challenge, and HcLBP/BPI2 in hepatopancreas was only up-regulated at 6 and 12 h after LPS challenge and at 12 h after A. hydrophila challenge. In addition, the recombinant HcLBP/BPIs displayed antibacterial activity against Gram-negative bacteria, and the antibacterial index of HcLBP/BPI1 was higher than that of HcLBP/BPI2. These results indicated that HcLBP/BPI1 and HcLBP/BPI2 probably played distinct roles in bacterial mediating immune response in Mollusca.

Uncovering the Complex Transcriptome Response of Mytilus chilensis against Saxitoxin: Implications of Harmful Algal Blooms on Mussel Populations

Saxitoxin (STX), a principal phycotoxin contributing to paralytic shellfish poisoning, is largely produced by marine microalgae of the genus Alexandrium. This toxin affects a wide range of species, inducing massive deaths in fish and other marine species. However, marine bivalves can resist and accumulate paralytic shellfish poisons. Despite numerous studies on the impact of STX in marine bivalves, knowledge regarding STX recognition at molecular level by benthic species remains scarce. Therefore, the aim of this study was to identify novel genes that interact with STX in the Chilean mussel Mytilus chilensis. For this, RNA-seq and RT-qPCR approaches were used to evaluate the transcriptomic response of M. chilensis to a purified STX as well as in vivo Alexandrium catenella exposure. Approximately 800 million reads were assembled, generating 138,883 contigs that were blasted against the UniProt Mollusca database. Pattern Recognition Receptors (PRRs) involved in mussel immunity, such as Toll-like receptors, tumor necrosis factor receptors, and scavenger-like receptors were found to be strongly upregulated at 8 and 16 h post-STX injection. These results suggest an involvement of PRRs in the response to STX, as well as identifying potential, novel STX-interacting receptors in this Chilean mussel. This study is the first transcriptomic overview of the STX-response in the edible species M. chilensis. However, the most significant contribution of this work is the identification of immune receptors and pathways potentially involved in the recognition and defense against STX's toxicity and its impact of harmful algae blooms on wild and cultivated mussel populations.

Transcriptome Analysis of the Sydney Rock Oyster, Saccostrea glomerata: Insights into Molluscan Immunity

BACKGROUND

Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic.

RESULTS

Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism.

CONCLUSIONS

This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

Exploring the immune signalling pathway-related genes of the cattle tick Rhipicephalus microplus: From molecular characterization to transcriptional profile upon microbial challenge

In dipteran insects, invading pathogens are selectively recognized by four major pathways, namely Toll, IMD, JNK, and JAK/STAT, and trigger the activation of several immune effectors. Although substantial advances have been made in understanding the immunity of model insects such as Drosophila melanogaster, knowledge on the activation of immune responses in other arthropods such as ticks remains limited. Herein, we have deepened our understanding of the intracellular signalling pathways likely to be involved in tick immunity by combining a large-scale in silico approach with high-throughput gene expression analysis. Data from in silico analysis revealed that although both the Toll and JAK/STAT signalling pathways are evolutionarily conserved across arthropods, ticks lack central components of the D. melanogaster IMD pathway. Moreover, we show that tick immune signalling-associated genes are constitutively transcribed in BME26 cells (a cell lineage derived from embryos of the cattle tick Rhipicephalus microplus) and exhibit different transcriptional patterns in response to microbial challenge. Interestingly, Anaplasma marginale, a pathogen that is naturally transmitted by R. microplus, causes downregulation of immune-related genes, suggesting that this pathogen may manipulate the tick immune system, favouring its survival and vector colonization.

cDNA and Gene Structure of MytiLec-1, A Bacteriostatic R-Type Lectin from the Mediterranean Mussel (Mytilus galloprovincialis)

MytiLec is an α-d-galactose-binding lectin with a unique primary structure isolated from the Mediterranean mussel (Mytilus galloprovincialis). The lectin adopts a β-trefoil fold that is also found in the B-sub-unit of ricin and other ricin-type (R-type) lectins. We are introducing MytiLec(-1) and its two variants (MytiLec-2 and -3), which both possess an additional pore-forming aerolysin-like domain, as members of a novel multi-genic “mytilectin family” in bivalve mollusks. Based on the full length mRNA sequence (911 bps), it was possible to elucidate the coding sequence of MytiLec-1, which displays an extended open reading frame (ORF) at the 5′ end of the sequence, confirmed both at the mRNA and at the genomic DNA sequence level. While this extension could potentially produce a polypeptide significantly longer than previously reported, this has not been confirmed yet at the protein level. MytiLec-1 was revealed to be encoded by a gene consisting of two exons and a single intron. The first exon comprised the 5′UTR and the initial ATG codon and it was possible to detect a putative promoter region immediately ahead of the transcription start site in the MytiLec-1 genomic locus. The remaining part of the MytiLec-1 coding sequence (including the three sub-domains, the 3′UTR and the poly-A signal) was included in the second exon. The bacteriostatic activity of MytiLec-1 was determined by the agglutination of both Gram-positive and Gram-negative bacteria, which was reversed by the co-presence of α-galactoside. Altogether, these data support the classification of MytiLec-1 as a member of the novel mytilectin family and suggest that this lectin may play an important role as a pattern recognition receptor in the innate immunity of mussels.

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

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2493-9) contains supplementary material, which is available to authorized users.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

MytiLec, a Mussel R-Type Lectin, Interacts with Surface Glycan Gb3 on Burkitt's Lymphoma Cells to Trigger Apoptosis through Multiple Pathways

MytiLec; a novel lectin isolated from the Mediterranean mussel (Mytilus galloprovincialis); shows strong binding affinity to globotriose (Gb3: Galα1-4Galβ1-4Glc). MytiLec revealed β-trefoil folding as also found in the ricin B-subunit type (R-type) lectin family, although the amino acid sequences were quite different. Classification of R-type lectin family members therefore needs to be based on conformation as well as on primary structure. MytiLec specifically killed Burkitt's lymphoma Ramos cells, which express Gb3. Fluorescein-labeling assay revealed that MytiLec was incorporated inside the cells. MytiLec treatment of Ramos cells resulted in activation of both classical MAPK/ extracellular signal-regulated kinase and extracellular signal-regulated kinase (MEK-ERK) and stress-activated (p38 kinase and JNK) Mitogen-activated protein kinases (MAPK) pathways. In the cells, MytiLec treatment triggered expression of tumor necrosis factor (TNF)-α (a ligand of death receptor-dependent apoptosis) and activation of mitochondria-controlling caspase-9 (initiator caspase) and caspase-3 (activator caspase). Experiments using the specific MEK inhibitor U0126 showed that MytiLec-induced phosphorylation of the MEK-ERK pathway up-regulated expression of the cyclin-dependent kinase inhibitor p21, leading to cell cycle arrest and TNF-α production. Activation of caspase-3 by MytiLec appeared to be regulated by multiple different pathways. Our findings, taken together, indicate that the novel R-type lectin MytiLec initiates programmed cell death of Burkitt's lymphoma cells through multiple pathways (MAPK cascade, death receptor signaling; caspase activation) based on interaction of the lectin with Gb3-containing glycosphingolipid-enriched microdomains on the cell surface.

IL-17 signaling components in bivalves: Comparative sequence analysis and involvement in the immune responses

The recent discovery of soluble immune-regulatory molecules in invertebrates takes advantage of the rapid growth of next generation sequencing datasets. Following protein domain searches in the transcriptomes of 31 bivalve spp. and in few available mollusk genomes, we retrieved 59 domains uniquely identifying interleukin 17 (IL-17) and 96 SEFIR domains typical of IL-17 receptors and CIKS/ACT1 proteins acting downstream in the IL-17 signaling pathway. Compared to the Chordata IL-17 family members, we confirm a separate clustering of the bivalve domain sequences and a consistent conservation pattern of amino acid residues. Analysis performed at transcript and genome level allowed us to propose an updated view of the components outlining the IL-17 signaling pathway in Mytilus galloprovincialis and Crassostrea gigas (in both species, homology modeling reduced the variety of IL-17 domains to only two 3D structures). Digital expression analysis indicated more heterogeneous expression levels for the mussel and oyster IL-17 ligands than for IL-17 receptors and CIKS/CIKSL proteins. Besides, new qPCR analyses confirmed such gene expression trends in hemocytes and gills of mussels challenged with heat-killed bacteria. These results uphold the involvement of an ancient IL-17 signaling pathway in the bivalve immune responses and, likewise in humans, suggest the possibility of distinctive modulatory roles of individual IL-17s/IL-17 receptors. Overall, the common evidence of pro-inflammatory cytokines and inter-related intracellular signaling pathways in bivalves definitely adds complexity to the invertebrate immunity.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.

Identification and Characterization of a Novel Family of Cysteine-Rich Peptides (MgCRP-I) from Mytilus galloprovincialis

We report the identification of a novel gene family (named MgCRP-I) encoding short secreted cysteine-rich peptides in the Mediterranean mussel Mytilus galloprovincialis. These peptides display a highly conserved pre-pro region and a hypervariable mature peptide comprising six invariant cysteine residues arranged in three intramolecular disulfide bridges. Although their cysteine pattern is similar to cysteines-rich neurotoxic peptides of distantly related protostomes such as cone snails and arachnids, the different organization of the disulfide bridges observed in synthetic peptides and phylogenetic analyses revealed MgCRP-I as a novel protein family. Genome- and transcriptome-wide searches for orthologous sequences in other bivalve species indicated the unique presence of this gene family in Mytilus spp. Like many antimicrobial peptides and neurotoxins, MgCRP-I peptides are produced as pre-propeptides, usually have a net positive charge and likely derive from similar evolutionary mechanisms, that is, gene duplication and positive selection within the mature peptide region; however, synthetic MgCRP-I peptides did not display significant toxicity in cultured mammalian cells, insecticidal, antimicrobial, or antifungal activities. The functional role of MgCRP-I peptides in mussel physiology still remains puzzling.

Characterizations and expression analyses of NF-κB and Rel genes in the Yesso scallop (Patinopecten yessoensis) suggest specific response patterns against Gram-negative infection in bivalves

Rel/NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) genes are evolutionarily conserved and play a pivotal role in several physiological events. They have been extensively studied from various species, including both vertebrates and invertebrates. However, the Rel/NF-κB genes have not been systematically characterized in bivalves. In this study, we identified and characterized PyNF-κB and PyRel in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine the identities and evolutionary relationships of Rel/NF-κB genes in Yesso scallop. Compared with the Rel/NF-κB genes from vertebrate species, the PyNF-κB and PyRel are relatively conserved in their structural features, but there were no paralogs found in P. yessoensis or other invertebrates. To gain insights into the roles of Rel/NF-κB genes during the innate immune response in scallop, quantitative real-time PCR was used to investigate the expression profiles of these genes at different developmental stages, in healthy adult tissues and in the hemolymph after bacterial infection with Micrococcus luteus and Vibrio anguillarum. The real-time PCR results indicated the abundance of PyNF-κB in the first four embryonic stages, including oocytes, fertilized eggs, morulae and blastulae. By contrast, PyRel was abundantly expressed in blastulae, trochophores and D-shaped larvae. In adult scallops, PyNF-κB and PyRel were ubiquitously expressed in most healthy tissues and highly expressed in most of the immune related tissues. Both genes were significantly up-regulated during the acute phase (3 h) after infection with Gram-positive (M. luteus) and negative (V. anguillarum) bacteria, while the much higher expression level of PyNF-κB suggested the involvement of the extra immune deficiency (IMD)-like pathway against the Gram-negative bacterial infection. The complex pattern of Rel/NF-κB induced expression suggested that PyNF-κB and PyRel both have specific and cooperative roles in the acute immune responses to bacterial infection.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

De novo transcriptome sequencing of the Octopus vulgaris hemocytes using Illumina RNA-Seq technology: response to the infection by the gastrointestinal parasite Aggregata octopiana

Background

Octopus vulgaris is a highly valuable species of great commercial interest and excellent candidate for aquaculture diversification; however, the octopus’ well-being is impaired by pathogens, of which the gastrointestinal coccidian parasite Aggregata octopiana is one of the most important. The knowledge of the molecular mechanisms of the immune response in cephalopods, especially in octopus is scarce. The transcriptome of the hemocytes of O. vulgaris was de novo sequenced using the high-throughput paired-end Illumina technology to identify genes involved in immune defense and to understand the molecular basis of octopus tolerance/resistance to coccidiosis.

Results

A bi-directional mRNA library was constructed from hemocytes of two groups of octopus according to the infection by A. octopiana, sick octopus, suffering coccidiosis, and healthy octopus, and reads were de novo assembled together. The differential expression of transcripts was analysed using the general assembly as a reference for mapping the reads from each condition. After sequencing, a total of 75,571,280 high quality reads were obtained from the sick octopus group and 74,731,646 from the healthy group. The general transcriptome of the O. vulgaris hemocytes was assembled in 254,506 contigs. A total of 48,225 contigs were successfully identified, and 538 transcripts exhibited differential expression between groups of infection. The general transcriptome revealed genes involved in pathways like NF-kB, TLR and Complement. Differential expression of TLR-2, PGRP, C1q and PRDX genes due to infection was validated using RT-qPCR. In sick octopuses, only TLR-2 was up-regulated in hemocytes, but all of them were up-regulated in caecum and gills.

Conclusion

The transcriptome reported here de novo establishes the first molecular clues to understand how the octopus immune system works and interacts with a highly pathogenic coccidian. The data provided here will contribute to identification of biomarkers for octopus resistance against pathogens, which could improve octopus farming in the near future.

Comparative study of immune responses in the deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the shallow-water mussel Mytilus galloprovincialis challenged with Vibrio bacteria

The deep-sea hydrothermal vent mussel Bathymodiolus azoricus and the continental European coast Mytilus galloprovincialis are two bivalves species living in highly distinct marine habitats. Mussels are filter-feeding animals that may accumulate rapidly bacteria from the environment. Contact with microorganism is thus inevitable during feeding processes where gill tissues assume a strategic importance at the interface between the external milieu and the internal body cavities promoting interactions with potential pathogens during normal filtration and a constant challenge to their immune system. In the present study B. azoricus and M. galloprovincialis were exposed to Vibrio alginolyticus, Vibrio anguillarum and Vibrio splendidus suspensions and to a mixture of these Vibrio suspensions, in order to ascertain the expression level of immune genes in gill samples, from both mussel species. The immune gene expressions were analyzed by means of quantitative-Polymerase Chain Reaction (qPCR). The gene expression results revealed that these bivalve species exhibit significant expression differences between 12 h and 24 h post-challenge times, and between the Vibrio strains used. V. splendidus induced the strongest gene expression level in the two bivalve species whereas the NF-κB and Aggrecan were the most significantly differentially expressed between the two mussel species. When comparing exposure times, both B. azoricus and M. galloprovincialis showed similar percentage of up-regulated genes at 12 h while a marked increased of gene expression was observed at 24 h for the majority of the immune genes in M. galloprovincialis. This contrasts with B. azoricus where the majority of the immune genes were down-regulated at 24 h. The 24 h post-challenge gene expression results clearly bring new evidence supporting time-dependent transcriptional activities resembling acute phase-like responses and different immune responses build-up in these two mussel species when challenged with Vibrio bacteria. High Pressure Liquid Chromatography (HPLC)-Electrospray ionization mass spectrometry (ESI-MS/MS) analyses resulted in different peptide sequences from B. azoricus and M. galloprovincialis gill tissues suggesting that naïve animals present differences, at the protein synthesis level, in their natural environment. B. azoricus proteins sequences, mostly of endosymbiont origin, were related to metabolic, energy production, protein synthesis processes and nutritional demands whereas in M. galloprovincialis putative protein functions were assumed to be related to structural and cellular integrity and signaling functions.

Target capture and massive sequencing of genes transcribed in Mytilus galloprovincialis

Next generation sequencing (NGS) allows fast and massive production of both genome and transcriptome sequence datasets. As the genome of the Mediterranean mussel Mytilus galloprovincialis is not available at present, we have explored the possibility of reducing the whole genome sequencing efforts by using capture probes coupled with PCR amplification and high-throughput 454-sequencing to enrich selected genomic regions. The enrichment of DNA target sequences was validated by real-time PCR, whereas the efficacy of the applied strategy was evaluated by mapping the 454-output reads against reference transcript data already available for M. galloprovincialis and by measuring coverage, SNPs, number of de novo sequenced introns, and complete gene sequences. Focusing on a target size of nearly 1.5 Mbp, we obtained a target coverage which allowed the identification of more than 250 complete introns, 10,741 SNPs, and also complete gene sequences. This study confirms the transcriptome-based enrichment of gDNA regions as a good strategy to expand knowledge on specific subsets of genes also in nonmodel organisms.