A primitive Toll-like receptor signaling pathway in mollusk Zhikong scallop Chlamys farreri

A newly identified scallop MyD88 interacts with TLR and functions in innate immunity

Myeloid differentiation factor-88 (MyD88) is a key adaptor of the toll-like receptor (TLR) signaling pathway and plays a crucial role in innate immune signal transduction in animals. However, the MyD88-mediated signal transduction mechanism in shellfish has not been well studied. In this study, a new MyD88 gene, CfMyD88-2, was identified in the Zhikong scallop, Chlamys farreri. The 1779 bp long open reading frame encodes 592 amino acids. The N-terminus of CfMyD88-2 contains a conserved death domain (DD), followed by a TIR (TLR/Interleukin-1 Receptor) domain. The results of the multi-sequence comparison showed that the TIR domain sequences were highly conserved. Phylogenetic analysis revealed that CfMyD88-2 was first associated with Mizuhopecten yessoensis MyD88-4 and Argopecten irradians MyD88-4. CfMyD88-2 mRNA was expressed in all scallop tissues, as detected by qRT-PCR, and the expression level was the highest in the mantle and hepatopancreas. In addition, CfMyD88-2 mRNA expression significantly increased after pathogen-associated molecular patterns (PAMPs, such as lipopolysaccharide, peptidoglycan, or polyinosinic-polycytidylic acid) stimulation. The results of the co-immunoprecipitation experiments in HEK293T cells showed that both CfMyD88-1 and CfMyD88-2 interacted with the TLR protein of scallops, suggesting the existence of more than one functional TLR-MyD88 signaling axis in scallops. Dual luciferase reporter gene assays indicated that the overexpressed CfMyD88-2 in HEK293T cells activated interferon (IFN) α, IFN-β, IFN-γ, and NF-κB reporter genes, indicating that the protein has multiple functions. The results of the subcellular localization experiment uncovered that CfMyD88-2 was mainly localized in the cytoplasm of human cells. In summary, the novel identified CfMyD88-2 can respond to the challenge of PAMPs, participate in TLR immune signaling, and may activate downstream effector genes such as NF-κB gene. These research results will be useful in advancing the theory of innate immunity in invertebrates and provide a reference for the selection of disease-resistant scallops in the future.

Glycosaminoglycan lyase: A new competition between bacteria and the pacific white shrimp Litopenaeus vannamei

Horizontal gene transfer (HGT) is an important evolutionary force in the formation of prokaryotic and eukaryotic genomes. In recent years, many HGT genes horizontally transferred from prokaryotes to eukaryotes have been reported, and most of them are present in arthropods. The Pacific white shrimp Litopenaeus vannamei, an important economic species of arthropod, has close relationships with bacteria, providing a platform for horizontal gene transfer (HGT). In this study, we analyzed bacteria-derived HGT based on a high-quality genome of L. vannamei via a homology search and phylogenetic analysis, and six HGT genes were identified. Among these six horizontally transferred genes, we found one gene (LOC113799989) that contains a bacterial chondroitinase AC structural domain and encodes an unknown glycosaminoglycan (GAG) lyase in L. vannamei. The real-time quantitative PCR results showed that the mRNA expression level of LOC113799989 was highest in the hepatopancreas and heart, and after stimulation by Vibrio parahaemolyticus, its mRNA expression level was rapidly up-regulated within 12 h. Furthermore, after injecting si-RNA and stimulation by V. parahaemolyticus, we found that the experimental group had a higher cumulative mortality rate in 48 h than the control group, indicating that the bacteria-derived GAG lyase can reduce the mortality of shrimp with respect to infection by V. parahaemolyticus and might be related to the resistance of shrimp to bacterial diseases. Our findings contribute to the study of the function of GAGs and provide new insights into GAG-related microbial pathogenesis and host defense mechanisms in arthropods.

A unique Ca2+-inhibited C-type lectin in shrimp Fenneropenaeuschinensis

C-type lectins (CTLs) are glycan-binding pattern recognition receptors (PRRs) that can bind to carbohydrates on pathogen surfaces, triggering immune responses in shrimp innate immunity. In this study, a unique Ca2+-inhibited CTL named FcLec was identified and characterized in Chinese shrimp Fenneropenaeus chinensis. The full-length cDNA sequence of FcLec was 976 bp (GenBank accession number KU361826), with a 615 bp open reading frame (ORF) encoding 204 amino acids. FcLec possesses a C-type lectin-like domain (CTLD) containing four conserved cysteines (Cys105, Cys174, Cys192, and Cys200) and two sugar-binding site structures (QPD and LNP). The tertiary structure of FcLec deduced revealed three α-helices and eight β-pleated sheets. The mRNA expression levels of FcLec in hemocytes and the hepatopancreas were markedly elevated after stimulation with Vibrio anguillarum and white spot syndrome virus (WSSV). The recombinant FcLec protein exhibited Ca2+-independent hemagglutination and bacterial agglutination, but these activities were observed only in the presence of EDTA to chelate metal ions. These findings suggest that FcLec plays important and functionally distinct roles in the shrimp's innate immune response to bacteria and viruses, enriching the current understanding of the relationship between CTL activity and Ca2+ in invertebrates.

The interactions between CpG oligodeoxynucleotides and Toll-like receptors in Pacific white shrimp Litopenaeus vannamei

CpG oligodeoxynucleotides (ODNs), as a novel type of adjuvant with immunomodulatory effects, are recognized by Toll-like receptors (TLRs) in Litopenaeus vannamei. In the present study, eleven LvTLRs-pCMV recombinants (rLvTLRs) were constructed to investigate the relationships between various CpG ODNs and different LvTLRs in human embryonic kidney 293T (HEK293T) cells, which was further confirmed by bio-layer interferometry (BLI) technique. The results of dual luciferase reporter assay showed that every LvTLR could activate multiple downstream genes, mainly including NF-κB, CREB, ISRE, IL-6-promoter, TNF-α-promoter and Myc, thereby inducing main signaling pathways in shrimps. Most CpG ODNs possessed affinities to more than one LvTLR, while each LvTLR could recognize multiple CpG ODNs, and the widely recognized ligands within CpG ODNs are A-class and B-class. Moreover, BLI analysis showed that CpG 2216, Cpg 2006, CpG 2143 and CpG 21425 exhibited dose-dependent affinity to the expressed TLR protein, which were consistent with the results in HEK293T cells. It suggested that the interactions of CpG ODNs with LvTLRs were indispensable for the immune regulation triggered by CpG ODNs, and these findings would lay foundations for studying the activations of LvTLRs to immune signaling pathways and shedding lights on the immune functions and mechanisms of CpG ODNs.

The Role of the TLR4-MyD88 Signaling Pathway in the Immune Response of the Selected Scallop Strain "Hongmo No. 1" to Heat Stress

The innate immunity of marine bivalves is challenged upon exposure to heat stress, especially with increases in the frequency and intensity of heat waves. TLR4 serves a classical pattern recognition receptor in recognizing pathogenic microorganisms and activating immune responses. In this study, three genes, HMTLR4, HMMyD88 and HMTRAF6, were characterized as homologs of genes in the TLR4-MyD88 signaling pathway in the selected scallop strain "Hongmo No. 1". According to RT-PCR, acute heat stress (32 °C) inhibited genes in the TLR4-MyD88 signaling pathway, and LPS stimulation-induced activation of TLR4-MyD88 signal transduction was also negatively affected at 32 °C. ELISA showed LPS-induced tumor necrosis factor alpha (TNF-α) or lysozyme (LZM) activity, but this was independent of temperature. RNA interference (RNAi) confirmed that HMTLR4 silencing suppressed the expression of its downstream gene, whether at 24 °C or at 32 °C. The level of TNF-α and the activity of LZM also decreased after injection with dsRNA, indicating a negative effect on the innate immunity of scallops. Additionally, acute heat stress affected the suppression of downstream gene expression when compared with that at 24 °C, which led us to the hypothesis that heat stress directly influences the downstream targets of HMTLR4. These results enrich the knowledge of scallop immunity under heat stress and can be beneficial for the genetic improvement of new scallop strains with higher thermotolerance.

A comprehensive review on the dynamic role of toll-like receptors (TLRs) in frontier aquaculture research and as a promising avenue for fish disease management

Toll-like receptors (TLRs) represent a conserved group of germline-encoded pattern recognition receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs) and play a crucial role in inducing the broadly acting innate immune response against pathogens. In recent years, the detection of 21 different TLR types in various fish species has sparked interest in exploring the potential of TLRs as targets for boosting immunity and disease resistance in fish. This comprehensive review offers the latest insights into the diverse facets of fish TLRs, highlighting their history, classification, architectural insights through 3D modelling, ligands recognition, signalling pathways, crosstalk, and expression patterns at various developmental stages. It provides an exhaustive account of the distinct TLRs induced during the invasion of specific pathogens in various fish species and delves into the disparities between fish TLRs and their mammalian counterparts, highlighting the specific contribution of TLRs to the immune response in fish. Although various facets of TLRs in some fish, shellfish, and molluscs have been described, the role of TLRs in several other aquatic organisms still remained as potential gaps. Overall, this article outlines frontier aquaculture research in advancing the knowledge of fish immune systems for the proper management of piscine maladies.

A metallothionein gene from hard clam Meretrix meretrix: Sequence features, expression patterns, and metal tolerance activities

Metallothioneins (MTs) are low-molecular weight cytoplasmic heavy metal binding proteins. MTs can regulate the concentration of essential or non-essential metals in organisms, and have many important biological functions, including detoxification, trace element metabolism, and anti-oxidation. In the present study, we cloned and characterized a metallothionein gene (designated as MmMT) from the hard clam Meretrix meretrix. The complete cDNA sequence of MmMT contained an open reading frame (ORF) of 629 bp, which encoded a protein of 76 amino acids with a predicted molecular mass of 7.66 kDa and a calculated theoretical isoelectric point of 7.24. MmMT is highly similar to previously identified MTs from other species, with typical metallothionein features such as a high cysteine residue content and the absence of histidine and aromatic residues. The mRNA transcripts of MmMT were prevalent in all the tested tissues, and the expression levels of MmMT were highest in the hepatopancreas and hemocytes. During the stimulation of Vibrio splendidus, the mRNA transcripts of MmMT in the hepatopancreas and hemocytes were significantly increased. The Escherichia coli overexpressing MmMT performed strong growth in the media supplemented with CdCl2 and CuSO4 compared to the control strains. These results provide useful information for further investigation of the functions of MmMT in metal detoxification and the innate immune system.

Identification and functional analysis of myeloid differentiation factor 88 (MyD88) in early development of Haliotis diversicolor

Myeloid differentiation factor 88 (MyD88) is a universal adaptor protein and plays an important role in the signal transduction of Toll like receptors (TLR) family. In this study, the MyD88 gene from the Haliotis diversicolor (hdMyD88) was identified. The full-length cDNA of hdMyD88 has a 1927 base pairs (bp), with an open reading frame of 1314 bp encoding 437 amino acids including a death domain (DD) at the N-terminus and TIR domain at the C-terminus which are typical features of MyD88 family proteins. Three conserved boxes are also found in the hdMyD88, which are similar to MyD88 in vertebrates. The expression levels of hdMyD88 mRNA at different early embryonic developmental stages of abalone were measured by qPCR revealed that their constitutive expression at all developmental stages analyzed with the considerably highest values at 8 cell stage and the lowest level at the trochosphere stage. Additionally, the mRNA expression of hdMyD88 decreased significantly (P < 0.05) after MyD88-dsRNA soak in the stage of trochosphere and veliger than EGFP-dsRNA group and blank control group. Whole embryo in situ hybridization showed that the positive signals of hdMyD88 were in visceral mass of trochophore larvae and veliger larvae. These results indicate hdMyD88 may could respond to pathogenic infection and may play an important role in early innate immunity in the process of abalone larval development.

Molecular characterization and expression analysis of a QM protein gene in Chinese mitten crab Eriocheir sinensis

QM protein was previously discovered as a tumor suppressor, and numerous studies have shown that QM protein also played important roles in the immune responses. To investigate the potential roles of the QM protein gene in Eriocheir sinensis, the QM protein gene (designated as EsQM) has been cloned from E. sinensis using the rapid amplification of cDNA ends (RACE) technique. The cDNA of EsQM is 781 bp in length, consisting of a 654 bp open reading frame (ORF), encoding 219 amino acids, a 27 bp 5' untranslated region (UTR) and a 94 bp 3' UTR. The EsQM protein has a calculated molecular weight of 25.4 kDa and a theoretical isoelectric point of 10.10. The deduced protein sequence of EsQM contains a Ribosomal_L16 domain, an SH3-binding motif, an N-acylation site, two putative antibiotic binding sites, two putative protein kinase C phosphorylation sites, and two amidation sites. EsQM is extremely conserved and exhibits more than 85% similarities to previously identified arthropod QM protein genes. By real-time quantitative PCR (qPCR) analysis, we found that EsQM mRNA transcripts were detectable in all the examined tissues, with the highest expression in hemocytes. The mRNA expression of EsQM in hemocytes was significantly upregulated after the stimulation of Aeromonas hydrophila or polybrominated diphenyl ether-47 (BDE-47). Moreover, EsQM mRNA expression in hemocytes responded more quickly and lasted longer when stimulated by A.hydrophila than BDE-47. Thus, EsQM can respond to bacterial infection and environmental pollution, and might be involved in the defense mechanism to both biological and non-biological stimulation of arthropods.

Genome-wide identification and characterization of Toll-like receptors (TLR) genes in Haliotis discus hannai, H. rufescens, and H. laevigata

Toll-like receptors (TLRs) play essential roles in the innate immune system and have been extensively studied in mollusks. In this study, through a genome-wide search, TLR genes were identified as 29 in Haliotis discus hannai, 33 in H. rufescens, and 16 in H. laevigata. Domain analysis indicated that these TLR genes contain leucine-rich repeat (LRR) and Toll/IL-1 receptor (TIR) domains and exons ranging from 1 to 5. Polymorphism analysis showed that the TLRs in abalones did not have high diversities with 143 SNPs and no Indel in H. discus hannai, 92 SNPs and 3 Indels together with 6 missense mutations in H. rufescens, and no SNP or Indel in H. laevigata. The expression of 8 TLR genes in H. discus hannai was confirmed in the hepatopancreas, gill, hemolymph, gonads, intestine, muscle, and mantle. The expression of five TLR genes (out of 8) in gills (p < 0.05), three in hepatopancreas (p < 0.05), and three in hemolymph (p < 0.05) was upregulated separately in response to the infection caused by Vibrio parahaemolyticus. The findings in this study would contribute to a better understanding of the molecular immune mechanism of H. discus hannai against stimulation by V. parahaemolyticus and provide a basis for the study of TLRs in abalones.

Regulation of innate immunity in marine mussel Mytilus coruscus: MicroRNA Mc-novel_miR_196 targets McTLR-like1 molecule to inhibit inflammatory response and apoptosis

Toll-like receptors (TLRs) are crucial players in immune recognition and regulation, with aberrant activation leading to autoimmune, chronic inflammatory, and infectious diseases. MicroRNAs (miRNAs) have been shown to regulate gene expression at transcriptional and post-transcriptional levels. While miRNA-mediated regulation of TLR signaling has been studied in mammals, the underlying mechanisms of TLR-miRNA interactions in molluscs remain unclear. In a previous study, one of the TLR genes potentially targeted by miRNAs was identified and named McTLR-like1. McTLR-like1 was later found to be targeted by miRNA Mc-novel_miR_196 through bioinformatic prediction. In this study, we aim to experimentally determine the interaction between McTLR-like1 and Mc-novel_miR_196, as well as their functional role in the innate immune response of molluscs. The results showed that the expression of Mc-novel_miR_196 was suppressed, while the expression of McTLR-like1 was enhanced in M. coruscus hemocytes treated with lipopolysaccharide (LPS). Moreover, in vitro assays demonstrated that Mc-novel_miR_196 directly targets the 5' UTR of McTLR-like1 and leads to the down-regulation of proinflammatory cytokines in hemocytes. In addition, co-transfection experiments confirmed that Mc-novel_miR_196 inhibits McTLR-like1 and inhibits the expression of proinflammatory cytokines. The Tunel assay also showed that Mc-novel_miR_196 inhibited apoptosis in hemocytes induced by LPS. Our findings suggest that microRNA Mc-novel_miR_196 acts as a regulator of innate immunity in M. coruscus by targeting McTLR-like1 and inhibiting inflammatory response and apoptosis. These results provide further insights into the complex molecular mechanisms underlying TLR signaling in molluscs.

Identification of olive flounder (Paralichthys olivaceus) toll-like receptor genes: Involvement in immune response to temperature stress and Edwardsiella tarda infection

Toll-like receptor (TLR) genes are best known for their roles in the innate immune defense. However, studies focusing on the reaction mechanisms of TLR genes in olive flounder (Paralichthys olivaceus) immune responses are still limited. In this study, 11 TLR family members (PoTLRs) were identified and classified from P. olivaceus genome. Phylogenetic analysis showed that PoTLRs were highly conserved in olive flounder. The analysis of motif prediction and gene structure indicated that TLRs had high sequence similarity. The expression patterns in developmental stages and different tissues showed that TLR members were spatially and temporally specific. RNA-Seq analysis of temperature stress and Edwardsiella tarda infection suggested that TLR members were involved in inflammatory responses, PoTLR5b and PoTLR22 showed significant differences in response to both temperature stress and E. tarda stress, indicating their potential immune functions. The results of this study suggested that TLR genes played important roles in the innate immune response of olive flounder, and would provide a solid basis for further study of their functions.

Seasonality in Synergism with Multi-Pathogen Presence Leads to Mass Mortalities of the Highly Endangered Pinna nobilis in Greek Coastlines: A Pathophysiological Approach

Mortalities of Pinna nobilis populations set at risk the survival of the species from many Mediterranean coastline habitats. In many cases, both Haplosporidium pinnae and Mycobacterium spp. are implicated in mass mortalities of P. nobilis populations, leading the species into extinction. In the context of the importance of these pathogens' role in P. nobilis mortalities, the present study investigated two Greek populations of the species hosting different microbial loads (one only H. pinnae and the second both pathogens) by the means of pathophysiological markers. More specifically, the populations from Kalloni Gulf (Lesvos Island) and from Maliakos Gulf (Fthiotis), seasonally sampled, were chosen based on the host pathogens in order to investigate physiological and immunological biomarkers to assess those pathogens' roles. In order to determine if the haplosporidian parasite possesses a major role in the mortalities or if both pathogens are involved in these phenomena, a variety of biomarkers, including apoptosis, autophagy, inflammation and heat shock response were applied. The results indicated a decreased physiological performance of individuals hosting both pathogens in comparison with those hosting only H. pinnae. Our findings provide evidence for the synergistic role of those pathogens in the mortality events, which is also enhanced by the influence of seasonality.

An LRR domain-containing membrane protein gene in rotifer Brachionus plicatilis: Sequence feature, expression pattern, and ligands binding activity

Leucine-rich repeat (LRR) domains mediate multiple innate immune responses via protein-ligand and protein-protein interactions, but their exact roles in invertebrates are poorly understood. Herein, an LRR domain-containing transmembrane protein (BpLRRm) was identified in the rotifer Brachionus plicatilis. The 1069 bp BpLRRm nucleotide sequence contains a 942 bp open reading frame (ORF) encoding a 313 amino acid polypeptide with four LRR motifs harbouring the LXXLXXLXLXXNXLXXL motif, and a transmembrane domain. Treatment with 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) decreased BpLRRm mRNA levels at 3 h, but they increased thereafter and peaked at 12 h. Lipopolysaccharide (LPS) treatment first increased BpLRRm mRNA levels at 3 h, but levels returned to normal at 12 h, then increased and peaked at 24 h. Recombinant BpLRRm protein bound pathogen-related molecular patterns (PAMPs), including LPS, peptidoglycan (PGN), glucan (GLU) and polyinosinic-polycytidylic acid (poly IC), in a dose-dependent manner. Thus, BpLRRm might function as a pattern recognition receptor (PRR) in the innate immunity of B. plicatilis, and mediate responses to environmental pollution.

Molecular cloning and characterization of a thioredoxin-like protein gene in rotifer Brachionus plicatilis

The thioredoxin-like protein exists widely, in various organisms, as a regulator of redox homeostasis. In this study, the full-length cDNA of a thioredoxin-like protein gene from rotifer Brachionus plicatilis (designated as BpTXNL) was obtained by 5' rapid amplification of cDNA end (RACE) technology. The complete cDNA of BpTXNL was 1111 bp, and contained a 5' untranslated region (UTR) of 69 bp, a 3' UTR of 163 bp with a polyadenylate additional signal and a polyadenylation site (PAS), and an open reading frame (ORF) of 878 bp, encoding 292 amino acids. The calculated molecular weight and the theoretical isoelectric point (pI) of the deduced BpTXNL peptide were 32.7 kDa and 4.97, respectively. The deduced protein sequence of BpTXNL contained a thioredoxin domain with the conserved redox-active site at ³³CGPC³⁶ and a proteasome-interacting thioredoxin (PITH) domain. Phylogenetic analysis demonstrated that BpTXNL was clustered with TXNLs of Strongyloides ratti and Caenorhabditis elegans. The temporal mRNA expression level of BpTXNL significantly decreased at 6 h, then increased to the peak 24h after the 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) challenge, while the mRNA transcripts of BpTXNL significantly increased and reached the peaks twice, at 6 h and 24 h after the lipopolysaccharide (LPS) challenge. The recombinant BpTXNL protein quickly exhibited a concentration-dependent antioxidant capacity and the peak occurred at 55 min in the 20 μM group. All these results showed that BpTXNL possesses an antioxidant capacity, and that it may be involved in the regulation of excessive reactive oxygen species (ROS) during environmental stress or pathogen invading.

Comparative study of five anti-lipopolysaccharide factor genes in Litopenaeus vannamei

Anti-lipopolysaccharide factors (ALFs) are a family of common innate immune effectors in crustaceans, and they exhibit broad spectrum antimicrobial activity. In this study, we identified and characterized five novel ALF genes (designated as LvALF1-5) from the Pacific white shrimp (Litopenaeus vannamei) to investigate their potential immune functions. The amino acid sequence alignments showed that LvALFs contained two conserved cysteine residues, a hydrophobic N-terminal region, and the conserved signature sequence W(T/K)CPG(S)WT(A). They all shared high similarity with previously reported ALFs and were clearly novel members of the ALF family. The mRNA transcripts of LvALFs were most highly expressed in hemocytes and the hepatopancreas. After shrimp were stimulated with Vibrio parahaemolyticus or white spot syndrome virus, expression of the LvALFs was significantly induced in hemocytes and the hepatopancreas with various expression profiles. Recombinant proteins of LvALFs exhibited potent bacteriostatic activity in vitro. Together, these results suggest that LvALF1-5 participate in the immune response of Pacific white shrimp against invading pathogens.

Comparative analysis of two cathepsin L genes in Asiatic hard clam (Meretrix meretrix): Similar in sequence features, different in expression profiles

Cathepsin L is widely found in eukaryotes and prokaryotes, and it plays important roles in innate immunity. In the present study, we cloned two cathepsin L genes (designated as MmCTSL1 and MmCTSL2, respectively) from Asiatic hard clam (Meretrix meretrix). The complete sequence of MmCTSL1 cDNA contained a 5' untranslated region (UTR) of 31 bp, a 3' UTR of 228 bp with a poly (A) tail, and an open reading frame (ORF) of 1005 bp encoding 334 amino acids with predicted molecular weight of 37.5 kDa and theoretical isoelectric point of 5.27, and contained a signal peptide (from M1 to A16), a protease inhibitor I29 family domain (from W27 to F87), and a papain family cysteine protease domain (from L118 to T333). The complete sequence of MmCTSL2 cDNA contained a 5' UTR of 50 bp, a 3' UTR of 162 bp with a poly (A) tail, and an ORF of 996 bp encoding a polypeptide of 331 amino acids with predicted molecular weight of 36.8 kDa and theoretical isoelectric point of 7.07. It contained a signal peptide (from M1 to A16), a protease inhibitor I29 family domain (from W30 to F89), and a papain family cysteine protease domain (from L115 to T330). Real-time quantitative PCR analysis demonstrated that MmCTSL1 and MmCTSL2 were widely expressed in all the tested tissues, including adductor muscle, foot, gill, hemocytes, hepatopancreas and mantle, with the highest mRNA expression level in hepatopancreas and hemocytes, respectively. After Vibrio splendidus challenge, the mRNA expression levels of MmCTSL1 and MmCTSL2 in hemocytes and hepatopancreas were both significantly up-regulated with different expression profiles. In hemocytes, the expression levels of MmCTSL1 and MmCTSL2 reached their respective peaks (3.4-fold and 13.0-fold compared with the control, respectively) at 12 h after bacterial challenge, and MmCTSL2 responds earlier than MmCTSL1. In hepatopancreas, the expression levels of MmCTSL1 and MmCTSL2 reached their respective peaks at 6 h (9.0-fold compared with the control) and 24 h (2.8-fold compared with the control) after bacterial challenge, meaning that MmCTSL1 responds earlier than MmCTSL2. At the same time, whether in hepatopancreas or hemocytes, MmCTSL1 persist for a while after the bacterial challenge peak, while MmCTSL2 would quickly return to the initial level after the bacterial challenge peak. These results indicate that cathepsin L may be involved in the immune process of hard clam against V. splendidus with different potential roles.

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.

CfIRF8-like interacts with the TBK1/IKKε family protein and regulates host antiviral innate immunity

The interferon regulatory factor (IRF) family, a class of transcription factors with key functions, are important in host innate immune defense and stress response. However, further research is required to determine the functions of IRFs in invertebrates. In this study, the coding sequence of an IRF gene was obtained from the Zhikong scallop (Chlamys farreri) and named CfIRF8-like. The open reading frame of CfIRF8-like was 1371 bp long and encoded 456 amino acids. Protein domain prediction revealed a typical IRF domain in the N-terminus of the CfIRF8-like protein and a typical IRF3 domain in the C-terminus. Multiple sequence alignment confirmed the conservation of the amino acid sequences of these two functional protein domains. Phylogenetic analysis showed that CfIRF8-like clustered with mollusk IRF8 proteins and then clustered with vertebrate IRF3, IRF4, and IRF5 subfamily proteins. Quantitative real-time PCR detected CfIRF8-like mRNA in all tested scallop tissues, with the highest expression in the gills. Simultaneously, the expression of CfIRF8-like transcripts in gills was significantly induced by polyinosinic-polycytidylic acid challenge. The results of protein interaction experiments showed that CfIRF8-like could directly bind the TBK1/IKKε family protein of scallop (CfIKK2) via its N-terminal IRF domain, revealing the presence of an ancient functional TBK1/IKKε-IRF signaling axis in scallops. Finally, dual-luciferase reporter assay results showed that the overexpression of CfIRF8-like in human embryonic kidney 293T cells could specifically activate the interferon β promoter of mammals and the interferon-stimulated response element promoter in dose-dependent manners. The findings of this preliminary analysis of the signal transduction and immune functions of scallop CfIRF8-like protein lay a foundation for an in-depth understanding of the innate immune function of invertebrate IRFs and the development of comparative immunology. The experimental results also provide theoretical support for the breeding of scallop disease-resistant strains.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.

A novel c-type lysozyme from Litopenaeus vannamei exhibits potent antimicrobial activity

Lysozyme is relevant to the innate immune system as a vital protein for crustaceans. In the present study, we cloned and characterized a novel c-type lysozyme gene (LvLYZ) from the Pacific white shrimp (Litopenaeus vannamei). The obtained full-length cDNA of LvLYZ was 990 bp and contained an open reading frame of 693 bp. Its deduced amino acid sequence consisted of 230 amino acids (aa) with a 17 aa signal peptide at the N-terminal and 130 aa functional domains. The multiple sequence alignment (MSA) indicated that the typical active sites in LvLYZ were similarly conserved as c-type lysozymes from other species. The transcription of LvLYZ appeared in all detected tissues and had relatively higher expression levels in hemocytes, hepatopancreas, gill and intestine. The mRNA expression profiles of LvLYZ were up-regulated in hemocyte and hepatopancreas post the stimulation of Vibrio parahaemolyticus or white spot syndrome virus (WSSV), respectively. The recombinant protein of LvLYZ (rLvLYZ) exhibited antibacterial activities against various microbes, including Escherichia coli, Vibrio splendidus, Micrococcaus luteus, Vibrio parahaemolyticus and Staphylococcus aureus. These results indicated that LvLYZ could cope with bacteria in L. vannamei and may play a significant role in immune response against invading pathogens.

A preliminary report of exploration of the exosomal shuttle protein in marine invertebrate Chlamys farreri

Exosomes are extracellular vesicles secreted by diverse cell under normal or abnormal physiological conditions, which could carry a range of bioactive molecules and play significant roles in biological processes, such as intercellular communication and immune response. In the current study, a preliminary study was performed to investigate the exosomal shuttle protein in Chlamys farreri (designated as CfesPro) and to predict the potential function of exosomes in scallop innate immunity. The serum derived exosomes (designated as CfEVs) were obtained from lipopolysaccharide (LPS)-stimulated C. farreri and untreated ones. After confirmation and characterization by transmission electron microscopy (TEM), nano-HPLC-MS/MS spectrometry was performed on CfEVs using a label-free quantitative method. Totally 2481 exosomal shuttle proteins were identified in CfEVs proteomic data, which included many innate immune related proteins. GO and KOG functional annotation showed that CfesPro participated in cellular processes, metabolism reactions, signaling transductions, immune responses and so on. Moreover, 1421 proteins in CfesPro were enriched to 324 pathways by KEGG analysis, including several immune-related pathways, such as autophagy, apoptosis and lysosome pathway. Meanwhile, eight autophagy-related proteins were initially identified in CfesPro, indicating that CfEVs had a potential role with autophagy. All these findings showed that CfEVs were involved in C. farreri innate immune defenses. This research would enrich the protein database of marine exosomes and provide a basis for the exploration of immune defense systems in marine invertebrates.

Effect of acute Cu exposure on immune response mechanisms of golden cuttlefish (Sepia esculenta)

Sepia esculenta is a common economic cephalopod that has received extensive attention due to the tender meat, rich protein content and certain medicinal value thereof. Over the past decade, multiple industries have discharged waste into the ocean in large quantities, thereby significantly increasing the concentration of heavy metals in the ocean. Copper (Cu) is a common heavy metal in the ocean. The increase of Cu content will affect numerous biological processes such as immunity and metabolism of marine organisms. High concentrations of Cu may inhibit S. esculenta growth, development, swimming, and other processes, which would significantly affect its culture. In this research, transcriptome analysis is used to initially explore Cu-exposed S. esculenta larval immune response mechanisms. And compared to control group with normally growing larvae, 2056 differentially expressed genes (DEGs) are identified in experimental group with Cu-exposed larvae. The results of DEGs functional enrichment analyses including GO and KEGG indicate that Cu exposure might promote inflammatory and innate immune responses in cuttlefish larvae. Then, 10 key genes that might regulate larval immunity are identified using a comprehensive analysis that combines protein-protein interaction (PPI) network and KEGG functional enrichment analyses, of which three genes with the highest number of protein interactions or involve in more KEGG signaling pathways are identified as hub genes that might significantly affect larval immune response processes. Comprehensive analysis of PPI network and KEGG signaling pathway are used for the first time to explore Cu-exposed S. esculenta larval immune response mechanisms. Our results preliminarily reveal immune response mechanisms of cephalopods exposed to heavy metals and provide valuable resources for further understanding mollusk immunity.

Integrated RNA-seq and RNAi Analysis of the Roles of the Hsp70 and SP Genes in Red-Shell Meretrix meretrix Tolerance to the Pathogen Vibrio parahaemolyticus

The "Wanlihong" Meretrix meretrix (WLH-M) clam is a new variety of this species that has a red shell and stronger Vibrio tolerance than ordinary M. meretrix (ORI-M). To investigate the molecular mechanisms responsible for the WLH-M strain's tolerance to Vibrio, we challenged clams with Vibrio parahaemolyticus and then assessed physiological indexes and conducted transcriptome analysis and RNA interference experiments. The mortality, tissue bacterial load, and hemocyte reactive oxygen species level of ORI-M were significantly higher than those of WLH-M, whereas the content and activity of lysozyme were significantly lower. Gene Ontology functional annotation analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that immune and metabolic pathways were enriched in Vibrio-challenged clams. The expressions of the heat shock protein 70 (Hsp70) and serine protease (SP) genes, which are involved in antibacterial immunity, were significantly upregulated in WLH-M but not in ORI-M, while the expression of the kynurenine 3-monooxygenase gene, a proinflammatory factor, was significantly downregulated in WLH-M. RNA interference experiments confirmed that Hsp70 and SP downregulation could result in increased mortality of WLH-M. Therefore, we speculate that Hsp70 and SP may be involved in the antibacterial immunity of WLH-M in vivo. Our data provided a valuable resource for further studies of the antibacterial mechanism of WLH-M and provided a foundation for the breeding of pathogen-resistant strains.

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.

Molecular and functional analysis of the microphthalmia-associated transcription factor (mitf) gene duplicates in red tilapia

In vertebrates, the microphthalmia-associated transcription factor (mitf) is at the hub of the melanin synthesis regulation network. However, little information is known about its molecular characterization, expression, location, or function in skin color differentiation and variation of red tilapia. The full-length cDNA sequences (1977 bp and 1999 bp) of mitfa and mitfb, encoding polypeptides of 491 and 514 amino acids, were effectively identified from red tilapia in this study. The Mitfa and Mitfb sequences of red tilapia clustered first with O. aureus, then with other teleost fish, according to phylogenetic analysis. Mitfa and mitfb mRNA were highly expressed in the brain, dorsal skin and eye tissues using quantitative real-time PCR. The mRNA expressions of mitfa and mitfb were the highest in the cleavage stage during the early development of red tilapia. Among three different colors of red tilapia, the expression levels of mitfa and mitfb were highest in the PB (pink with scattered black spots) dorsal skin. After overwintering, the mitfa and mitfb mRNA expressions were high in the dorsal skin of PB (color changed from pink to black). Mitfa and mitfb were mostly found in the epidermal layer of the dorsal skin, according to in situ hybridization (ISH) analysis. After injecting mitf-dsRNA duplicates along the tail vein of red tilapia, the activity of tyrosinase and the level of melanin in the dorsal skin both decreased significantly. The mRNA expressions of mitfa and its downstream genes (tyrb, tyrp1a and dct) decreased, whereas the mRNA expression of mitfb increased after mitfa-dsRNA injection. The mRNA expressions of mitfb, tyrb, tyrp1a and dct decreased, whereas the mRNA expression of mitfa increased after injecting mitfb-dsRNA. These findings suggest that mitf gene duplicates may play an important role in red tilapia skin color differentiation and variation via the melanogenesis pathway.

A novel molluscan TLR molecule engaged in inflammatory response through MyD88 adapter recruitment

Toll-like receptors (TLRs) mediated signaling plays a vital role in activating innate and adaptive immunity. Although TLR mediated signaling has been comprehensively investigated in mammalian species, the mechanisms underlying TLR signaling in molluscs remain obscure. In the present study, a novel TLR isoform namely McTLR-like1 was identified in the thick shell mussel Mytilus coruscus. McTLR-like1 was highly expressed in molluscan immune-related tissues, and its transcriptional levels in hemocytes were significantly increased when challenged by V. alginolyticus. McTLR-like1 activated nuclear factor κB (NF-κB) and strengthened the transcription and phosphorylation of NF-κB subunit P65 in mammalian cells. Upon the silencing of McTLR-like1, the mRNA expression levels of pro-inflammatory cytokines were down-regulated, and the animals exhibited higher levels of resistance when challenged with V. alginolyticus. McMyD88a mRNA expression was also downregulated alongside McTLR-like1. Furthermore, GST-pull down assays revealed a visible affinity between McTLR-like1 and McMyD88a. Collectively, these results demonstrated that the newly identified gene affiliated to the molluscan TLR family and plays a role in the TLR-mediated activation of inflammatory response via its affinity with MyD88. The present study enhances our knowledge of TLR signaling mechanisms in molluscs and provides new insights into the evolution of TLRs.

Scallop IKK1 Responds to Bacterial and Virus-Related Pathogen Stimulation and Interacts With MyD88 Adaptor of Toll-Like Receptor Pathway Signaling

IKK proteins are key signaling molecules in the innate immune system of animals, and act downstream of pattern recognition receptors. However, research on IKKs in invertebrates, especially marine mollusks, remains scarce. In this study, we cloned CfIKK1 gene from the Zhikong scallop (Chlamys farreri) and studied its function and the signaling it mediates. The open reading frame of CfIKK1 was 2190 bp and encoded 729 amino acids. Phylogenetic analysis showed that CfIKK1 belonged to the invertebrate IKKα/IKKβ family. Quantitative real-time PCR analysis revealed the ubiquitous expression of CfIKK1 mRNA in all scallop tissues and challenge with lipopolysaccharide, peptidoglycan, or poly(I:C) significantly upregulated the expression of CfIKK1. Co-immunoprecipitation assays confirmed the interaction of CfIKK1 with scallop MyD88 (Myeloid differentiation actor 88, the key adaptor of the TLR signaling pathway) via its N-terminal kinase domain. Additionally, CfIKK1 protein could form homodimers and even oligomers, with N-terminal kinase domain and C-terminal scaffold dimerization domain playing key roles in this process. Finally, the results of RNAi experiments showed that when the scallop IKK1 gene was suppressed, the expression of IRF genes also decreased significantly. In conclusion, CfIKK1 could respond to PAMPs challenge and interact with MyD88 protein of scallop TLR signaling, with the formation of CfIKK1 dimers or oligomers. At the same time, the results of RNAi experiments revealed the close regulatory relationship between IKK1 and IRF genes of scallop. Therefore, as a key signal transduction molecule and immune activity regulator, CfIKK1 plays important roles in the innate immune system of scallops.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Pathophysiological Responses of Pinna nobilis Individuals Enlightens the Etiology of Mass Mortality Situation in the Mediterranean Populations

Due to the rapid decrease of Pinna nobilis populations during the previous decades, this bivalve species, endemic in the Mediterranean Sea, is characterized as 'critically endangered'. In addition to human pressures, various pathogen infections have resulted in extended reduction, even population extinction. While Haplosporidium pinnae is characterized as one of the major causative agents, mass mortalities have also been attributed to Mycobacterium sp. and Vibrio spp. Due to limited knowledge concerning the physiological response of infected P. nobilis specimens against various pathogens, this study's aim was to investigate to pathophysiological response of P. nobilis individuals, originating from mortality events in the Thermaikos Gulf and Lesvos and Limnos islands (Greece), and their correlation to different potential pathogens detected in the diseased animals. In isolated tissues, several cellular stress indicators of the heat shock and immune response, apoptosis and autophagy, were examined. Despite the complexity and limitations in the study of P. nobilis mortality events, the present investigation demonstrates the cumulative negative effect of co-infection additionally with H. pinnae in comparison to the non-presence of haplosporidian parasite. In addition, impacts of global climate change affecting physiological performance and immune responses result in more vulnerable populations in infectious diseases, a phenomenon which may intensify in the future.

An Explanation Based on Energy-Related Changes for Blue Mussel Mytilus edulis Coping With Seawater Acidification

As ocean acidification (OA) is gradually increasing, concerns regarding its ecological impacts on marine organisms are growing. Our previous studies have shown that seawater acidification exerted adverse effects on physiological processes of the blue mussel Mytilus edulis, and the aim of the present study was to obtain energy-related evidence to verify and explain our previous findings. Thus, the same acidification system (pH: 7.7 or 7.1; acidification method: HCl addition or CO2 enrichment; experimental period: 21d) was set up, and the energy-related changes were assessed. The results showed that the energy charge (EC) and the gene expressions of cytochrome C oxidase (COX) reflecting the ATP synthesis rate increased significantly after acidification treatments. What's more, the mussels exposed to acidification allocated more energy to gills and hemocytes. However, the total adenylate pool (TAP) and the final adenosine triphosphate (ATP) in M. edulis decreased significantly, especially in CO2 treatment group at pH 7.1. It was interesting to note that, TAP, ATP, and COXs gene expressions in CO2 treatment groups were all significantly lower than that in HCl treatment groups at the same pH, verifying that CO2-induced acidification exhibited more deleterious impacts on M. edulis, and ions besides H+ produced by CO2 dissolution were possible causes. In conclusion, energy-related changes in M. edulis responded actively to seawater acidification and varied with different acidification conditions, while the constraints they had at higher acidification levels suggest that M. edulis will have a limited tolerance to increasing OA in the future.

Dorsal regulates the expression of two phage lysozymes acquired via horizontal gene transfer in triangle sail mussel Hyriopsis cumingii

Dorsal is a Rel/NF-κB transcription factor, which forms a key part of the Toll pathway. Lysozyme is a ubiquitous enzyme that degrades bacterial cell walls. In this study, a Dorsal homolog was cloned and characterized from triangle sail mussel Hyriopsis cumingii, namely, HcDorsal. Dorsal consisted of 3041 bp, including a 1938 bp open reading frame encoding a 645 amino acid protein. The deduced HcDorsal protein contained a Rel homology domain and an Ig-like, plexin, transcription factor domain. Analysis of expression patterns showed that HcDorsal was highly expressed in the hepatopancreas of H. cumingii. The expression level of HcDorsal continuously increased after Vibrio parahaemolyticus stimulation. When HcDorsal was knocked down by siRNA interference, two phage lysozyme genes (HcLyso1 and HcLyso2) obtained by horizontal gene transfer were significantly downregulated in hemocytes of mussels. Furthermore, knockdown of HcLyso1 and HcLyso2 could weaken V. parahaemolyticus clearance ability. Recombinant HcLyso1 and HcLyso2 proteins accelerated the bacterial clearance in vivo in mussels and evidently inhibited the growth of V. parahaemolyticus. These results suggested that HcDorsal could be activated after V. parahaemolyticus stimulation and then modulate the immune response through the transcriptional regulation of HcLyso1 and HcLyso2, thereby playing a protective role in mussels.

A Toll-like receptor identified in Gigantidas platifrons and its potential role in the immune recognition of endosymbiotic methane oxidation bacteria

Symbiosis with chemosynthetic bacteria is an important ecological strategy for the deep-sea megafaunas including mollusks, tubeworms and crustacean to obtain nutrients in hydrothermal vents and cold seeps. How the megafaunas recognize symbionts and establish the symbiosis has attracted much attention. Bathymodiolinae mussels are endemic species in both hydrothermal vents and cold seeps while the immune recognition mechanism underlying the symbiosis is not well understood due to the nonculturable symbionts. In previous study, a lipopolysaccharide (LPS) pull-down assay was conducted in Gigantidas platifrons to screen the pattern recognition receptors potentially involved in the recognition of symbiotic methane-oxidizing bacteria (MOB). Consequently, a total of 208 proteins including GpTLR13 were identified. Here the molecular structure, expression pattern and immune function of GpTLR13 were further analyzed. It was found that GpTLR13 could bind intensively with the lipid A structure of LPS through surface plasmon resonance analysis. The expression alternations of GpTLR13 transcripts during a 28-day of symbiont-depletion assay were investigated by real-time qPCR. As a result, a robust decrease of GpTLR13 transcripts was observed accompanying with the loss of symbionts, implying its participation in symbiosis. In addition, GpTLR13 transcripts were found expressed exclusively in the bacteriocytes of gills of G. platifrons by in situ hybridization. It was therefore speculated that GpTLR13 may be involved in the immune recognition of symbiotic methane-oxidizing bacteria by specifically recognizing the lipid A structure of LPS. However, the interaction between GpTLR13 and symbiotic MOB was failed to be addressed due to the nonculturable symbionts. Nevertheless, the present result has provided with a promising candidate as well as a new approach for the identification of symbiont-related genes in Bathymodiolinae mussels.

Characterization of a novel toll-like receptor and activation NF-κB signal pathway in triangle sail mussel Hyriopsis cumingii

Toll-like receptor (TLR) family plays an important role in innate immunity for detection of and defense against microbial pathogens. In this study, a novel toll-like receptor (HcTLRn) was characterized from freshwater pearl mussel H. cumingii. The complete sequence of HcTLRn was 3725 bp, and the open reading frame (ORF) encoded 718 amino acid residues. Predicted HcTLRn protein possessed seven atypical leucine-rich repeat (LRR) domains, two typical LRR subfamily domains, a C-terminal domain LRR, a transmembrane domain and an intracellular Toll/interleukin-1 (IL-1) receptor domain. Transcripts of HcTLRn were constitutive expressed in the tissues of healthy mussels and were markedly induced in hepatopancreas and gills after lipopolysaccharide (LPS), peptidoglycan (PGN) and polyinosinic polycytidylic acid (ploy I: C) stimulation. Knockdown of HcTLRn in vivo significantly decreased the mRNA levels of TLR pathway transcription factors p65 and p105 as well as antimicrobial peptides (AMPs) including lysozyme (HcLys), theromacin (HcTher), whey acidic protein (HcWAP), LPS-binding protein/bactericidal permeability increasing protein (HcLBP/BPI) 1 and 2 after mussels challenged by LPS. In situ hybridization results showed that HcTLRn mRNA was significantly increased in hemocytes after LPS, PGN and poly I:C stimulation. HcTLRn protein was mainly expressed in hepatopancreas and gills and was significantly increased after LPS stimulation. Moreover, recombinant extracellular domain of HcTLRn (HcTLRn-ECD) proteins could bind to a variety of bacterial and pathogen-associated molecular patterns such as LPS, PGN, and poly I:C in vitro. Subcellular localization results showed that HcTLRn was mainly distributed near the cell membrane and in cytoplasm. Over-expression of HcTLRn activated the NF-κB luciferase reporter in HEK293T cells. Collectively, these results suggested that HcTLRn was a TLR family member that might play an important role in activation of NF-κB signal pathway in Mollusca.

Transcriptome-wide identification and characterization of toll-like receptors response to Vibrio anguillarum infection in Manila clam (Ruditapes philippinarum)

The Manila clam (Ruditapes philippinarum), one of the major marine aquaculture species in China, is susceptible to infection with the pathogen Vibrio, which results in massive mortality and economic losses. Toll-like receptors (TLRs) are significant pattern recognition receptors (PRRs) of innate immunity that are involved in immune regulation against pathogenic invasion. Molecular characterization of Manila clam TLRs and investigations of their immune functions are essential to prevent and control Vibrio infection. In the present research, eight cDNA sequences of R. philippinarum TLRs (RpTLRs) were identified from previous transcriptome libraries and then classified into four groups, namely, P-TLR (one sequence), V-TLR (one sequence), Ls-TLR (two sequences) and sP-TLR (four sequences), based on the corresponding LRR domain arrangement of their protein structures within the typical TLR motifs. A selective pressure test firstly suggested that the molluscan P-TLR, V-TLR, Ls-TLR and sP-TLR families underwent positive selection, and different numbers of positive selection sites (PSSs) were identified in different domains of the four types of RpTLRs, as determined by PAML and analysis of website data. These findings indicated that the evolution of RpTLRs may be associated with their immune recognition and function. Furthermore, tissue-specific expression analysis showed that all RpTLRs were ubiquitously expressed in all test tissues and were dominant in hemocytes. Quantitative real-time PCR revealed that the cDNA expression of all eight RpTLRs was upregulated after injection with Vibrio anguillarum (P < 0.01) in R. philippinarum hemocytes, revealing that these RpTLRs play important roles in responding to pathogenic stimulation. In summary, these findings provide a foundation for future investigations of the molecular classification and evolutionary patterns of Toll-like receptors in invertebrates, and the innate immune responses of TLR signaling pathways in Mollusca.

Identification of a Novel Pattern Recognition Receptor DM9 Domain Containing Protein 4 as a Marker for Pro-Hemocyte of Pacific Oyster Crassostrea gigas

DM9 refers to an uncharacterized protein domain that is originally discovered in Drosophila melanogaster. Two proteins with DM9 repeats have been recently identified from Pacific oyster Crassostrea gigas as mannose-specific binding pattern-recognition receptors (PRRs). In the present study, a novel member of DM9 domain containing protein (designated as CgDM9CP-4) was identified from C. gigas. CgDM9CP-4, about 16 kDa with only two tandem DM9 domains, was highly enriched in hemocytes and gill. The transcripts level of CgDM9CP-4 in circulating hemocytes were decreased after LPS, PGN and Vibrio splendidus stimulations. The recombinant protein of CgDM9CP-4 (rCgDM9CP-4) displayed a broad binding spectrum towards various pathogen-associated molecular patterns (PAMPs) (LPS, PGN, β-glucan and Mannose) and microorganisms (Staphylococcus aureus, Micrococcus luteus, V. splendidus, V. anguillarum, Escherichia coli, Pichia pastoris and Yarrowia lipolytica). CgDM9CP-4 was mostly expressed in gill and some of the hemocytes. Flow cytometry analysis demonstrated that the CgDM9CP-4-positive hemocytes accounted for 7.3% of the total hemocytes, and they were small in size and less in granularity. CgDM9CP-4 was highly expressed in non-phagocytes (~82% of total hemocytes). The reactive oxygen species (ROS) and the expression levels of cytokines in CgDM9CP-4-positive hemocytes were much lower than that in CgDM9CP-4-negative hemocytes. The mRNA expression level of CgDM9CP-4 in hemocytes was decreased after RNAi of hematopoietic-related factors (CgGATA, CgRunt, CgSCL, and CgNotch). In addition, CgDM9CP-4-positive cells were found to be much more abundant in hemocytes from gill than that from hemolymph, with most of them located in the gill filament. All these results suggested that CgDM9CP-4 was a novel member of PRR that expressed in undifferentiated pro-hemocytes to mediate immune recognition of pathogens.

Identification, characterization and expression analyses of cholinesterases genes in Yesso scallop (Patinopecten yessoensis) reveal molecular function allocation in responses to ocean acidification

Cholinesterases are key enzymes in central and peripheral cholinergic nerve system functioning on nerve impulse transmission in animals. Though cholinesterases have been identified in most vertebrates, the knowledge about the variable numbers and multiple functions of the genes is still quite meagre in invertebrates, especially in scallops. In this study, the complete cholinesterase (ChE) family members have been systematically characterized in Yesso scallop (Patinopecten yessoensis) via whole-genome scanning through in silico analysis. Ten ChE family members in the genome of Yesso scallop (designated PyChEs) were identified and potentially acted to be the largest number of ChE in the reported species to date. Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of these genes. The expression profiles of PyChEs were determined in all developmental stages, in healthy adult tissues, and in mantles under low pH stress (pH 6.5 and 7.5). Spatiotemporal expression suggested the ubiquitous functional roles of PyChEs in all stages of development, as well as general and tissue-specific functions in scallop tissues. Regulation expressions revealed diverse up- and down-regulated expression patterns at most time points, suggesting different functional specialization of gene superfamily members in response to ocean acidification (OA). Evidences in gene number, phylogenetic relationships and expression patterns of PyChEs revealed that functional innovations and differentiations after gene duplication may result in altered functional constraints among PyChEs gene clusters. Collectively, our results provide the potential clues that the selection pressures coming from the environment were the potential inducement leading to function allocation of ChE family members in scallop.

NKL-24: A novel antimicrobial peptide derived from zebrafish NK-lysin that inhibits bacterial growth and enhances resistance against Vibrio parahaemolyticus infection in Yesso scallop, Patinopecten yessoensis

The extensive use of antibiotics in aquaculture has resulted in the prevalence of antibiotic-resistant bacteria and, consequently, new antibacterial strategies or drugs with clear modes of action are urgently needed. Antimicrobial peptides (AMPs) are currently widely considered as alternatives to antibiotics in the treatment of infections in aquatic animals. In this study, we aimed to evaluate the effects of NKL-24, a truncated peptide derived from zebrafish NK-lysin, against Yesso scallop (Patinopecten yessoensis) pathogen, Vibrio parahaemolyticus. The results showed that NKL-24 had a potent antibacterial effect against V. parahaemolyticus via a membrane active cell-killing mechanism. The in vitro study showed that sub-lethal levels of NKL-24 obviously reduced bacterial swimming ability and downregulated the transcription of the selected genes associated with V. parahaemolyticus virulence. Studies on NKL-24 biosafety in hemocytes and in Yesso scallop have shown no adverse effects from this peptide. Bacteria challenge test results demonstrated that NKL-24 significantly decreased the mortality and inhibited bacterial growth in the scallop infected with V. parahaemolyticus, while further in vivo examination revealed that NKL-24 could enhance non-specific immune parameters. Moreover, NKL-24 was capable of modulating a series of V. parahaemolyticus-responsive genes in the scallop. These results suggest the protective action of NKL-24 against V. parahaemolyticus and the potential of this peptide as a promising candidate for aquaculture applications.

Genome-wide characterization of Toll-like receptors in black rockfish Sebastes schlegelii: Evolution and response mechanisms following Edwardsiella tarda infection

As one of the key components of pattern recognition receptors, Toll-like receptors (TLRs) play pivotal roles in the innate immune system. However, little information is available about the TLR genes in Sebastes schlegelii. In the present study, 17 TLR genes were identified and classified based on the whole genome database. Tandem duplication events in TLR1, TLR2, TLR5 and TLR13 played major role in the expansion of S. schlegelii TLR genes; both TLR2-3 and TLR2-4 had the same largest number of introns/exons, 11 exons and 10 introns. The syntenic analysis showed neighboring genes of TLR genes were most conserved in S. schlegelii and in L. crocea. Phylogenetic and evolutionary analysis showed that these TLR genes were divided into five subfamilies and exhibited different selection pressures. Meanwhile, the expression patterns of TLR genes in the intestine after E. tarda infection were investigated by qRT-PCR. Finally, protein and protein interaction (PPI) network analysis indicated that TLR genes interacted with IFN-inducible genes, inflammatory cytokines, and participated in MyD88-dependent pathway. In summary, this study provided valuable information for further functional characterization of TLR genes in S. schlegelii.

Genome-wide identification, characterization of RLR genes in Yesso scallop (Patinopecten yessoensis) and functional regulations in responses to ocean acidification

Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), are crucial sensors with a conserved structure in cytoplasm, inducing the production of cytokines, chemokines and host restriction factors which mediate a variety of intracellular activities to interfere with distinct PAMPs (pathogen-associated molecular patterns) for eliminating pathogens in innate immune system. Although RLR genes have been investigated in most vertebrates and some invertebrates, the systematic identification and characterization of RLR genes have not been reported in scallops. In this study, four RLR genes (PY-10413.4, PY-10413.5, PY-443.7 and PY-443.8, designated PyRLRs) were identified in Yesso scallop (Patinopecten yessoensis) through whole-genome scanning through in silico analysis, including two pairs of tandem duplicate genes located on the same scaffold (PY-10413.4 and PY-10413.5, PY-443.7 and PY-443.8, respectively). Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of these genes. The expression profiles of PyRLRs were determined in all developmental stages, in healthy adult tissues, and in mantles that simulated ocean acidification (OA) exposure (pH = 6.5 and 7.5) at different time points (3, 6, 12 and 24 h). Spatiotemporal expression patterns suggested the functional roles of PyRLRs in all stages of development and growth of the scallop. Regulation expressions revealed PY-10413.4 and PY-10413.5 with one or two CARD(s) (caspase activation and recruitment domain) were up-regulated expressed at most time points, whereas PY-443.8 and PY-10413.4 without CARD were significantly down-regulated at each time points, suggesting functional differentiations in the two pairs of PyRLRs based on the structural differences in response to OA. Collectively, this study demonstrated gene duplication of RLR family genes and provide primary analysis for versatile roles in the response of the bivalve innate immune system to OA challenge.

Genome-wide identification and characterization of toll-like receptor genes in spotted sea bass (Lateolabrax maculatus) and their involvement in the host immune response to Vibrio harveyi infection

Toll-like receptor (TLR) genes are the earliest reported pathogen recognition receptors (PRRs) and have been extensively studied. These genes play pivotal roles in the innate immune defense against pathogen invasion. In this study, a total of 16 tlr genes were identified and characterized in spotted sea bass (Lateolabrax maculatus). The tlr genes of spotted sea bass were classified into five subfamilies (tlr1-subfamily, tlr3-subfamily, tlr5-subfamily, tlr7-subfamily, and tlr11-subfamily) according to the phylogenetic analysis, and their annotations were confirmed by a syntenic analysis. The protein domain analysis indicated that most tlr genes had the following three major TLR protein domains: a leucine-rich repeat (LRR) domain, a transmembrane region (TM) and a Toll/interleukin-1 receptor (TIR) domain. The tlr genes in spotted sea bass were distributed in 11 of 24 chromosomes. The mRNA expression levels of 16 tlr genes in response to Vibrio harveyi infection were quantified in the head kidney. Most genes were downregulated following V. harveyi infection, while only 5 tlr genes, including tlr1-1, tlr1-2, tlr2-2, tlr5, and tlr7, were significantly upregulated. Collectively, these results help elucidate the crucial roles of tlr genes in the immune response of spotted sea bass and may supply valuable genomic resources for future studies investigating fish disease management.

A novel invertebrate toll-like receptor is involved in TLR mediated signal pathway of thick shell mussel Mytilus coruscus

Toll-like receptors (TLRs) are the most well studied pattern recognition receptors (PRRs) that play a crucial role in both innate and adaptive immunity in animals. In the present study, a novel toll-like receptor (McTLRj) was identified and characterised in thick shell mussel Mytilus coruscus. McTLRj possessed a signal peptide, a transmembrane domain, leucine-rich repeats and an intracellular Toll/interleukin-1 receptor domain that were conserved in typical TLRs. McTLRj transcripts were constitutively expressed in all of the examined tissues with high expression level in immune-related tissues, and significantly induced in haemocytes upon live Vibrio alginolyticus, lipopolysaccharide, polyinosinic-polycytidylic acid and peptidoglycans challenge. The overexpression of the McTLRj TLR fragment in Drosophila S2 cells could induce the expression of Drosophila attacin A, drosomycin, cecropin A, and metchnikowin expression. The expression of McTLRj was obviously repressed by dsRNA-mediated RNA interference, and downstream TLR pathway factors, such as MyD88a, IRAK4, and TRAF6 were significantly repressed in McTLRj-silenced mussels upon LPS challenge. These results collectively indicated that McTLRj is a TLR family member that may play a potential PRR role in TLR-mediated signalling pathway. This research contributed to the clarification of innate immune response in molluscs.

Molecular Characterization of Two Toll Receptors in Hyriopsis cumingii and Their Potential Roles in Antibacterial Response

Tolls/Toll-like receptors (TLRs) play a key role in innate immunity by detecting the invading microbes and subsequently activating downstream signaling cascades. In this study, two new molluscan Toll members (designed as HcToll6 and HcToll7) were identified from triangle-shell pearl mussel (Hyriopsis cumingii). The obtained HcToll6 full-length cDNA was 3207 bp consisting of a 2223 bp open reading frame (ORF) that encoded a peptide of 740 amino acids. HcToll7 cDNA is a 3216 bp molecule that contains an ORF of 2139 bp encoding a protein of 712 amino acids. The deduced HcToll6 and HcToll7 proteins share two common structures: extracellular leucine-rich repeat (LRR) domains and intracellular Toll/interleukin-1 receptor (TIR) domain. Quantitative real-time PCR results showed that HcToll6 and HcToll7 were mainly expressed in the hepatopancreas and the gills, and they responded rapidly to bacterial stimulation. RNA interference by dsRNA results revealed that HcToll6 and HcToll7 RNAi strongly decreased the expression of lysozyme (HcLyso) and defensin (HcDef) in the gills of RNAi-treated mussels with Vibrio parahaemolyticus challenge. As a pattern recognition receptor, the prokaryotic expressed the recombinant LRR domains of HcToll6 and HcToll7 (rHcToll6-LRR and rHcToll7-LRR) could bind to Gram-positive and Gram-negative bacteria and had a strong tendency to recognize lipopolysaccharide (LPS) and peptidoglycan (PNG). rHcToll6-LRR and rHcToll7-LRR exhibited a significant in vitro bactericidal activity against V. parahaemolyticus and Staphylococcus aureus. These findings provide useful information to characterize Tolls in mussels.

Molecular cloning of the Rab7 effector RILP (Rab-interacting lysosomal protein) in Litopenaeus vannamei and preliminary analysis of its role in white spot syndrome virus infection

To investigate the role of the Rab7 effector RILP (Rab-interacting lysosomal protein) in white spot syndrome virus (WSSV) infection, the full-length cDNA of RILP (LvRILP) was cloned in Litopenaeus vannamei, which consists of 1595 bp and encodes a polypeptide of 411 amino acids. Sequence analysis and multiple sequence alignment displayed that LvRILP contained a conserved RILP region from 277 amino acid to 325 amino acid. Both the LvRILP and Rab7 mRNA were most highly expressed in stomach and most lowly expressed in hemocyte, which were significantly up-regulated and exhibited similar kinetics post WSSV infection. The interaction of Rab7 with LvRILP was verified by both GST Pull-down and ELISA. Meanwhile, the results of Pull-down assays showed that the GST-tagged VP28 (GST-VP28), His-tagged Rab7 (His-Rab7) and His-RILP formed a tripartite complex. After silencing by specific LvRILP dsRNA, the LvRILP mRNA level exhibited a significant reduction, and the expression levels of three WSSV genes ie1, wsv477 and vp28 all exhibited decreases at 24, 36 and 48 h post WSSV infection. These results suggested that the Rab7 effector RILP was involved in WSSV infection.

Evolutionary Origins of Toll-like Receptor Signaling

Toll-like receptors (TLRs) are transmembrane pattern recognition receptors that are best known for their roles in innate immunity for the detection of and defense against microbial pathogens. However, TLRs also have roles in many nonimmune processes, most notably development. TLRs direct both immune and developmental programs by activation of downstream signaling pathways, often by activation of the NF-κB pathway. There are two primary TLR subtypes: 1) TLRs with multiple cysteine clusters in their ectodomain (mccTLRs) and 2) TLRs with a single cysteine cluster in their ectodomain (sccTLRs). For some time, it has been known that TLRs and the biological processes that they control are conserved in organisms from insects to mammals. However, genome and transcriptome sequencing has revealed that many basal metazoans also have TLRs and downstream NF-κB signaling components. In this review, we discuss what is known about the structure, biological function, and downstream signaling pathways of TLRs found in phyla from Porifera through Annelida. From these analyses, we hypothesize that mccTLRs emerged in the phylum Cnidaria, that sccTLRs evolved in the phylum Mollusca, and that TLRs have dual immune and developmental biological functions in organisms as ancient as cnidarians.

A novel invertebrate toll-like receptor with broad recognition spectrum from thick shell mussel Mytilus coruscus

Toll-like receptors (TLRs) are a category of most well recognized pattern recognition molecules that act on a vital role in both innate and adaptive immunity. In the present study, a novel toll-like receptor (McTLRw) was identified and characterized in thick shell mussel Mytilus coruscus. McTLRw possesses one intracellular Toll/interleukin-1 (IL-1) receptor (TIR) domain, one transmembrane region (TM), one leucine rich repeat N-terminal domain (LRR_NT) and a few of leucine-rich repeats (LRRs), which all are common in typical TLRs. McTLRw transcripts were constitutively expressed in all examined tissues with higher expression levels in immune related tissues, and were significantly induced in haemocytes with the challenges of live Vibrio alginolyticus, lipopolysaccharide (LPS), peptidoglycans (PGN) and β-glucan (GLU), but not induced by polyinosinic-polycytidylic acid (poly I:C). rMcTLRw exhibited affinity to LPS, PGN and GLU while no affinity to poly I:C. Further, the downstream of TLR signaling pathway myeloid differentiation factor 88a (MyD88a), interleukin-1 receptor-associated kinase-4 (IRAK4) and tumor necrosis factor receptor-associated factor 6 (TRAF6) were significantly repressed in McTLRw silenced mussels while challenged with LPS. These results collectively indicated that McTLRw is one member of TLR family and involved in immune response to against invaders by taking participate in TLR mediated signaling pathway.

The transcriptional response of the Pacific oyster Crassostrea gigas under simultaneous bacterial and heat stresses

Bacterial infection and heat stress are considered as two major environmental threats for the aquaculture industry of oyster Crassostrea gigas. In the present study, the expression profiles of mRNA transcripts in the hemocytes of oysters under bacterial challenge and heat stress were examined by next-generation sequencing. There were 21,095, 21,957 and 21,141 transcripts identified in the hemocytes of oysters from three groups, respectively, including control group (designated as Con group), Vibrio splendidus challenge group (Bac group), and bacterial and heat stress combined treatment group (BacHeat group). There were 4610, 5093 and 5149 differentially expressed transcripts (DTs) in the three pairwise comparisons Con/Bac, Con/BacHeat and Bac/BacHeat, respectively. The main enriched GO terms in biological process category of the DTs included the metabolic processes, cellular process, response to stimulus and immune system process. The expression patterns of DTs involved in pattern recognition, immune signal transduction and energy metabolic indicated that the immune response to bacterial challenge was disturbed under acute heat stress, which was also confirmed by quantitative real-time PCR. The neuroendocrine immunomodulation, especially the catecholaminergic regulation, played indispensable roles in stress response. The total energy reserves as well as cellular energy allocation (CEA) in hepatopancreas of oysters decreased remarkably especially in BacHeat group, while the energy consumption generally increased, suggesting that the immune defense against the simultaneous stimulation of pathogen and heat stress imposed greater costs on oyster's energy expenditure than a single stressor. These results above indicated that, the heat stress disturbed the normal expression of genes involved in immune response and energy metabolism, accelerated energy consumption and broke the metabolic balance, leading to a decline in resilience to infection and mass mortality of oyster in summer.

Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata

Background

The most critical step in the pearl formation during aquaculture is issued to the proliferation and differentiation of outer epithelial cells of mantle graft into pearl sac. This pearl sac secretes various matrix proteins to produce pearls by a complex physiological process which has not been well-understood yet. Here, we aimed to unravel the genes involved in the development of pearl sac and pearl, and the sequential expression patterns of different shell matrix proteins secreted from the pearl sac during pearl formation by pearl oyster Pinctada fucata using high-throughput transcriptome profiling.

Results

Principal component analysis (PCA) showed clearly different gene expression profiles between earlier (before 1 week) and later stages (1 week to 3 months) of grafting. Immune-related genes were highly expressed between 0 h – 24 h (donor dependent) and 48 h – 1 w (host dependent), and in the course of wound healing process pearl sac was developed by two weeks of graft transplantation. Moreover, for the first time, we identified some stem cell marker genes including ABCG2, SOX2, MEF2A, HES1, MET, NRP1, ESR1, STAT6, PAX2, FZD1 and PROM1 that were expressed differentially during the formation of pearl sac. The expression profiling of 192 biomineralization-related genes demonstrated that most of the shell matrix proteins (SMPs) involved in prismatic layer formation were first up-regulated and then gradually down-regulated indicating their involvement in the development of pearl sac and the onset of pearl mineralization. Most of the nacreous layer forming SMPs were up-regulated at 2 weeks after the maturation of pearl sac. Nacrein, MSI7 and shematrin involved in both layer formation were highly expressed during 0 h – 24 h, down-regulated up to 1 week and then up-regulated again after accomplishment of pearl sac formation.

Conclusions

Using an RNA-seq approach we unraveled the expression pattern of the key genes involved in the development of pearl sac and pearl as a result of host immune response after grafting. These findings provide valuable information in understanding the molecular mechanism of pearl formation and immune response in P. fucata.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5579-3) contains supplementary material, which is available to authorized users.