Cloning and gene expression of signal transducers and activators of transcription (STAT) homologue provide new insights into the immune response and nucleus graft of the pearl oyster Pinctada fucata

Molecular confirmation and functional study of signal transducer and activator of transcription genes in the Pacific oyster Crassostrea gigas

The JAK (Janus kinase)-STAT (Signal transducer and activator of transcription) is a well-known functional signaling pathway that plays a key role in several important biological activities such as apoptosis, cell proliferation, differentiation, and immunity. However, limited studies have explored the functions of STAT genes in invertebrates. In the present study, the gene sequences of two STAT genes from the Pacific oyster (Crassostrea gigas), termed CgSTAT-Like-1 (CgSTAT-L1) and CgSTAT-Like-2 (CgSTAT-L2), were obtained using polymerase chain reaction (PCR) amplification and cloning. Multiple sequence comparisons revealed that the sequences of crucial domains of these proteins were conserved, and the similarity with the protein sequence of other molluscan STAT is close to 90 %. The phylogenetic analyses indicated that CgSTAT-L1 and CgSTAT-L2 are novel members of the mollusk STAT family. Quantitative real-time PCR results implied that CgSTAT-L1 and CgSTAT-L2 mRNA expression was found in all tissues, and significantly induced after challenge with lipopolysaccharide (LPS), peptidoglycan (PGN), or poly(I:C). After that, dual-luciferase reporter assays denoted that overexpression of CgSTAT-L1 and CgSTAT-L2 significantly activated the NF-κB signaling, and, interestingly, the overexpressed CgSTAT proteins potentiated LPS-induced NF-κB activation. These results contributed a preliminary analysis of the immune-related function of STAT genes in oysters, laying the foundation for deeper understanding of the function of invertebrate STAT genes.

Exosome-derived small non-coding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca)

Exosomes, a subset of small extracellular vesicles, carry various nucleic acids, proteins, lipids, amino acids and metabolites. They function as a mode of intercellular communication and molecular transfer. Exosome cargo molecules, including small non-coding RNAs (sncRNAs), are involved in the immune response in various organisms. However, the role of exosome-derived sncRNAs in immune responses in molluscs remains unclear. Here, we aimed to reveal the sncRNAs involved in the immune response during grafting transplantation by the pearl oyster Pinctada fucata. Exosomes were successfully extracted from the P. fucata haemolymph during graft transplantation. Abundant microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) were simultaneously discovered in P. fucata exosomes by small RNA sequencing. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to immune response upon grafting transplantation, whereas piRNAs may also be associated with transposon silencing by targeting with genome transposon elements. This work provides the basis for a functional understanding of exosome-derived sncRNAs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in molluscs.

The primitive interferon-like system and its antiviral function in molluscs

The phylum mollusca is a very important group in the animal kingdom for the large number and diversified species. Recently, interest in molluscan immunity has increased due to their phylogenetic position and importance in worldwide aquaculture and aquatic environment. As the main aquaculture animal, most molluscs live in the water environment and they have to cope with many pathogen challenges, in which virus is one of the primary causes for the mass mortality. In vertebrates, interferon (IFN) system is generally recognized as the first line of defence against viral infection, while the antiviral mechanisms in molluscs remain to be clearly illuminated. Recently, some IFN-like proteins and IFN-related components have been characterized from molluscs, such as pattern recognition receptors (PRRs), interferon regulatory factors (IRFs), IFN-like receptors, JAK/STAT and IFN-stimulated genes (ISGs), which reinforce the existence of IFN-like system in molluscs. This system can be activated by virus or poly (I:C) challenges and further regulate the antiviral response of haemocytes in molluscs. This review summarizes the research progresses of IFN-like system in molluscs with the emphases on the uniformity and heterogeneity of IFN-like system of molluscs compared to that of other animals, which will be helpful for elucidating the antiviral modulation in molluscs and understanding the origin and evolution of IFN system.

Functional characterization of STATa/b genes encoding transcription factors from Branchiostoma belcheri

The signal transducer and activator of transcription (STAT), as an important transcription factor of the Janus kinase (JAK)-STAT signaling pathway, is pivotal for development and immunity and well documented in vertebrates. However, the STAT gene has not been reported in chordate amphioxus (Branchiostoma belcheri). In this study, we firstly identify and characterize two STAT genes from Branchiostoma belcheri (designed as AmphiSTATa and AmphiSTATb). Secondly, our results reveal that AmphiSTATa is clustered with vertebrate STAT1, STAT2, STAT3 and STAT4, whereas AmphiSTATb is grouped with STAT5 and STAT6 based on phylogenetic analysis. Thirdly, AmphiSTATa and AmphiSTATb are found to widely express in five representative tissues of amphioxus (gill, hepatic cecum, intestine, muscle and notochord) by RT-qPCR analysis. Importantly, both AmphiSTATa and AmphiSTATb can be involved in innate immune responses to LPS stimulation. Fourthly, we demonstrate that AmphiSTATa and AmphiSTATb can form homodimers or heterodimers by Co-IP and Native-PAGE assay, and that AmphiSTATa and AmphiSTATb proteins can also distribute in cytoplasm and nucleus by the subcellular localization. Taken together, our findings not only reveal the roles of AmphiSTATa and AmphiSTATb in amphioxus innate immune responses to LPS stimulation, but provide a new insight into further elucidating the evolution and function of STATs in animals.

-Omic Analysis of the Sepia officinalis White Body: New Insights into Multifunctionality and Haematopoiesis Regulation

Cephalopods, like other protostomes, lack an adaptive immune system and only rely on an innate immune system. The main immune cells are haemocytes (Hcts), which are able to respond to pathogens and external attacks. First reports based on morphological observations revealed that the white body (WB) located in the optic sinuses of cuttlefish was the origin of Hcts. Combining transcriptomic and proteomic analyses, we identified several factors known to be involved in haematopoiesis in vertebrate species in cuttlefish WB. Among these factors, members of the JAK-STAT signaling pathway were identified, some of them for the first time in a molluscan transcriptome and proteome. Immune factors, such as members of the Toll/NF-κB signaling pathway, pattern recognition proteins and receptors, and members of the oxidative stress responses, were also identified, and support an immune role of the WB. Both transcriptome and proteome analyses revealed that the WB harbors an intense metabolism concurrent with the haematopoietic function. Finally, a comparative analysis of the WB and Hct proteomes revealed many proteins in common, confirming previous morphological studies on the origin of Hcts in cuttlefish. This molecular work demonstrates that the WB is multifunctional and provides bases for haematopoiesis regulation in cuttlefish.

Characterization of cyclin dependent kinase-7 and its differential response to grafting challenge in the black shell colored selected line of pearl oyster Pinctada fucata martensii

Cyclin dependent kinase-7 (Cdk-7) is a protein kinase associated with regulating the cell cycle, cell differentiation and proliferation, apoptosis and inflammatory response. This study characterized the full cDNA sequence of Cdk-7 in Pinctada fucata martensii (PmCdk-7). A full length sequence of 1473bp with an open reading frame (ORF) of 915bp and encodes a 304aa, 5'-UTR of 58bp and a 3'-UTR of 500bp was obtained. The construed amino acid sequence of PmCdk-7 comprised of a Serine/Threonine protein kinases, catalytic (S_TKc) domain with a protein kinases ATP-binding region signature (14-38aa) and the serine/Threonine protein kinases active-site signature (129-141aa) within the domain. Tissue distribution analysis revealed a high relative mRNA expression of PmCdk-7 within haemocytes. Following the insertion operation (grafting), the relative expression levels of PmCdk-7 in the haemocyte was expressed differentially among the studied groups; the black shell colored selected line (BS) and the control group (CG). High expression was recorded between 12 h and 5d with a peak at 3d suggesting a heightened level of DNA replication and inflammatory response during the pearl-sac formation and this expression was higher in BS than CS showcasing, the heightened immune capacity of BS to grafting operation. Immune stimulation experiment with bacterial endotoxin and a viral mimic revealed PmCdk-7 response to pathogenic stress. The results from our study showed that PmCdk-7 performs a vital function during the cell cycle by aiding DNA replication and also aid response to inflammations generated due to the incision from the grafting operation and long exposure to immune-stimulants (pathogens).

Characterization of thioredoxin-like PROTEIN-5 (TRXLP-5) and its differential response to grafting challenge in the black coloured selected line and control stocks of Pinctada fucata martensii

Thioredoxin-like protein 5 (Trxlp-5) is a thioredoxin isoform associated with cellular redox homeostasis through the activity of thiol-disulfide reductase. In our study, Trxlp-5 was identified and characterized in Pinctada fucata martensii. The expression of PmTrxlp-5 was detected in response to polyinosinic: polycytidylic acid (poly I:C) and lipopolysaccharides (LPS) stimulation. The differences in PmTrxlp-5 expression were evaluated between the black coloured selected line and the control stock after grafting operation. The open reading frame (ORF) consisted of 1167bp encoding a 388 amino acid, 5'-UTR of 41bp and a 3'-UTR of 846bp. PmTrxlp-5 exhibited a conserved WCXXC functional motif similar to thioredoxins from other species. Tissue analysis showcased the highest relative mRNA expressions of PmTrxlp-5 in the haemocytes. Interestingly, after the grafting operation, mRNA expression of PmTrxlp-5 in the haemocytes was differentially expressed post grafting with a peak 6 h after grafting suggesting the high involvement of the gene in immune response in the early stage after grafting. The black coloured selected line group (BS) had significantly higher expression than the control group (CG) at 24 h, 6 d and 30 d after grafting operation. PmTrxlp-5 also showed a wave-like pattern in mRNA expression after bacterial endotoxin LPS and viral mimic poly I:C. These results suggested that PmTrxlp-5 plays a vital function in cellular redox homeostasis and immune response against grafting operation and pathogenic infections and can be used as a gene marker for selective breeding programs.

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.

Molecular cloning and characterization of a putative mitogen-activated protein kinase (Erk1/2) gene: Involvement in mantle immunity of Pinctada fucata

Extracellular signal-regulated kinases (ERKs) are conserved and related with protein-serine/threonine kinases that participate in the regulation of multiple biological processes, such as cell survival, cell differentiation, proliferation, metabolism, and inflammation. However, little is known about the roles of this kinase in the pearl oyster. In this study, we cloned and identified an ERK homolog from Pinctada fucata (PfErk). Furthermore, we have unraveled its expressional kinetics after lipopolysaccharide (LPS) and polyinosinic-epolycytidylic acid (poly I:C) immune challenge. Pferk harbored a 5' untranslated region (UTR) of 12 bp, a coding sequence of 1074 bp, and a 3' UTR of 882 bp. The putative peptide comprised a predicted molecular mass of 41.19 kDa, with a theoretical pI of 6.15. Sequence analysis showed that it possesses one STK catalytic domain and a conserved His-Arg-Asp (HRD) domain. In addition, a canonical Thr-Glu-Tyr (TEY) dual phosphorylation motif and an ATRW substrate binding site were also identified in the coding protein. Homology assessment of PfErk showed high similarity to Homo sapiens ERK. Phylogenetic analysis supported a close evolutionary relationship with molluscan orthologs. The expression patterns of Pferk were observed in seven different tissues of pearl oyster, with highest expression in the mantle and lowest expression in the digestive gland. Pferk mRNA expression levels were detected at developmental stages, with the highest expression in D-shaped larvae, followed by the 32-cell stage. The mRNA expression of Pferk was upregulated significantly in P. fucata mantle primary cells and mantle tissue after LPS and poly (I:C) treatment, and PfErk phosphorylation levels were activated by LPS and poly (I:C) challenges. Overall, our results suggested that PfErk may play important roles in pearl oyster innate immunity, and provided a new understanding of mantle immunity in the pearl oyster.

Identification and analysis of an MKK4 homologue in response to the nucleus grafting operation and antigens in the pearl oyster, Pinctada fucata

The mitogen-activated protein kinase kinase 4 (MKK4) is a key component of the c-Jun N-terminal kinase (JNK) signaling pathway and regulates multiple cellular activities. However, little is known about the roles of this kinase in pearl oyster. In this study, we identified an MKK4 homologue in Pinctada fucata by using a transcriptome database. Sequence analysis and protein structure prediction showed that PfMKK4 is highly conserved to MKK4 from other vertebrate and invertebrate species. Phylogenetic analysis revealed that PfMKK4 has the closest relationship with that from Crassostrea gigas. QPCR was used to investigate expression profiles in different healthy adult tissues and developmental stages of P. fucata. We found that PfMKK4 was ubiquitously expressed in all tissues and developmental stages examined except for in D-shaped larvae. Gene expression analysis suggested that PfMKK4 is involved in the response to the nucleus insertion operation. Lipopolysaccharide (LPS) and polyinosinic:polycytidylic acid [poly(I:C)] stimulation in vivo reduced PfMKK4 mRNA expression at 6 h, 48 h and 48 h, 72 h, respectively. LPS and poly(I:C) induced PfMKK4 phosphorylation in a primary mantle cell culture. These results contribute to better understanding of the potential role played by PfMKK4 in protecting the pearl oyster from injury caused by grafting or disease.

Deep transcriptome profiling sheds light on key players in nucleus implantation induced immune response in the pearl oyster Pinctada martensii

Immunological rejection of the pearl oysters following nucleus implantation is a major issue limiting the successful rate of cultured pearls. To date, the molecular mechanism of immune tolerance during pearl formation in the pearl oysters is still largely unknown. Through the RNA sequencing platform and comparative transcriptomic analysis, we investigated the chronic gene expression changes at seven time points (0, 5, 10, 15, 20, 30, 60 days post implantation or dpi) over a period of 60 days following nucleus implantation in the pearl oyster Pinctada martensii. A total of 81,390 unique transcripts (or unigenes) with a combined length of 96.8 million bp and a N50 value of 2227 bp were obtained. When compared with sequences in the nr, nt, Swiss-Prot, KEGG, COG and GO databases, 36,380 unigenes can find homologous genes. Pairwise comparison of gene expression among all the samples showed that the largest number (or 6846) of differentially expressed genes was observed at 10 dpi. The number then decreased to below 5000 at 15, 20 and 30 dpi and increased again to 6679 at 60 dpi. PCA analysis further showed that the seven time points can be roughly divided into four groups. Comparative transcriptomic analysis between the four groups identified a variety of genes showing differential expression at different time points, including many immune-related genes such as those encoding for toll-like receptor, lectin, scavenger receptor, and peroxidase. In addition, GO and KEGG enrichment analysis revealed that these differentially expressed genes were mainly associated with metabolism, ribosome function, immune response, signaling transduction, and cytoskeleton organization. Notably, two KEGG pathways, namely "cell adhesion molecules" and "primary immunodeficiency" were significantly enriched during the whole process. This finding indicates that genes in these pathways are likely to play critical roles in the immune tolerance of the pearl oysters. To conclude, the data obtained contribute to a better understanding of the molecular mechanisms of nucleus implantation induced immune response in the pearl oysters, and will facilitate the development of effective measures to improve the performance of pearl culture.

Three STATs are involved in the regulation of the expression of antimicrobial peptides in the triangle sail mussel, Hyriopsis cumingii

Janus kinase (Jak) and signal transducers and activators of transcription (STAT) signaling pathway is associated in antiviral and antibacterial immune response. Previous studies primarily investigated the function of STATs in mammals. For most invertebrates, only one STAT was found in each species, such as STAT92E was found in Drosophila melanogaster. The studies, which focus on the functional difference between various STATs in the same species of invertebrate, are limited. In the present study, three STATs (HcSTAT1, HcSTAT2 and HcSTAT3) were identified in triangle shell pearl mussel, Hyriopsis cumingii. Phylogenetic analysis showed that HcSTAT1 and HcSTAT3 were clustered with Homo sapiens STAT5, and HcSTAT2 was clustered with Pinctada fucata STAT and Crassostea gigas STAT6. All three STATs could be detected in all tested tissues (hemocytes, hepatopancreas, gill, mantle and foot), and were induced expression when challenged with Staphylococcus aureus or Aeromonas hydrophilia in hemocytes and hepatopancreas. HcSTAT1 regulated the expression of HcDef, HcWAP, HcThe and HcTNF. The expression of HcWAP and HcTNF was down-regulated in HcSTAT2-RNAi mussel. And HcSTAT3 affected the expression of HcTNF. The study is the first report of different functions in antibacterial immune responses between STATs in mollusks.

Differentially expressed immune-related genes in hemocytes of the pearl oyster Pinctada fucata against allograft identified by transcriptome analysis

The pearl oyster Pinctada fucata is commonly cultured for marine pearls in China. To culture pearls, a mantle piece from a donor pearl oyster is grafted with a nucleus into a receptor. This transplanted mantle piece may be rejected by the immune system of the recipient oyster, thus reducing the success of transplantation. However, there have been limited studies about the oyster's immune defense against allograft. In this study, hemocyte transcriptome analysis was performed to detect the immune responses to allograft in P. fucata at 0 h and 48 h after a transplant. The sequencing reaction produced 92.5 million reads that were mapped against the reference genome sequences of P. fucata. The Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to identify all immune-related differentially expressed genes (DEGs). Compared with patterns at 0 h, a total of 798 DEGs were identified, including 410 up-regulated and 388 down-regulated genes at 48 h. The expression levels of interleukin receptor and toll-like receptor in hemocytes were increased significantly 48 h post-transplant, indicating that the oyster immune response was induced. Finally, altered levels of 18 randomly selected immune-related DEGs were confirmed by quantitative real-time PCR (qRT-PCR). Our results provide the basis for further analysis of the immune rejection of allotransplantation.

Serum immune response of pearl oyster Pinctada fucata to xenografts and allografts

The mantle piece from the donor pearl oyster would be rejected by the immune system of recipient oyster in pearl culture practice, especially in the case that the donor and receptor are different species. Thus, investigation of the immune response of recipient oyster to grafted mantle pieces, particularly to xenografts, is of importance in creating xenograft transplantation technology for pearl culture industry. The humoral immune responses of P. fucata to allograft (mantle piece of P. fucata) and xenografts (mantle pieces of P. maxima and P. margaritifera, respectively) were studied in this paper. The oysters receiving no transplantations were served as the control group. The serum was collected from recipient P. fucata at 1 d, 2 d, 3 d, 4 d, 5 d, 7 d, 9 d, 11 d, 13 d, and 15 d, respectively after transplantation, and the serum antibacterial activity, lysozyme activity (LZM), alkaline phosphatase (AKP), acid phosphatase (ACP), total antioxidant capacity (TAC), and agglutination to rabbit red blood cells were investigated. The result indicated that serum of both the experimental groups and the control group can agglutinate rabbit red blood cells, with variation between groups and between time points, respectively. The antibacterial activity in the experimental group was significantly higher than that in the control group at 2-4 d, but lower at 5-11 d and returned back to normal at 15 d, with significant differences among experimental groups (P < 0.05). The LZM in the experimental group was significantly higher than that in the control group at 3-7 d, with significant differences in bacteriolytic activity among various groups (P < 0.05). Both the ACP and AKP activity levels in the experimental groups were higher than those in the control group at 2-9 d, with significant differences among various groups at 3-9 d (P < 0.05). The TAC level in the experimental groups was higher than that in the control group at 1-7 d, with significant differences among various groups at 4-7 d (P < 0.05). The above results indicated that all of the humoral immune factors investigated showed immune responses to both allografts and xenografts, with no specific to any of them. Thus, it is necessary to further screen immune rejection factors specific to xenografts, including cellular immune components.

Hemocytes participate in calcium carbonate crystal formation, transportation and shell regeneration in the pearl oyster Pinctada fucata

In this study, light microscope, scanning and transmission electron microscope, hematoxylin-eosin and fluorescent staining, and mass spectrometry methods were employed to observe the calcium carbonate (CaCO3) crystal formation, hemocyte release and transportation, and hemocyte distribution at the shell regeneration area and to analyse the proteome of hemocytes in the pearl oyster, Pinctada fucata. The results indicated that intracellular CaCO3 crystals were observed in circulating hemocytes in P. fucata, implying that there was a suitable microenvironment for crystal formation in the hemocytes. This conclusion was further supported by the proteome analysis, in which various biomineralization-related proteins were detected. The crystal-bearing hemocytes, mainly granulocytes, may be released to extrapallial fluid (EPF) by the secretory cavities distributed on the outer surface of the mantle centre. These granulocytes in the EPF and between the regenerated shells were abundant and free. In the regenerated prismatic layer, the granulocytes were fused into each column and fragmented with the duration of shell maturation, suggesting the direct involvement of hemocytes in shell regeneration. Overall, this study provided evidence that hemocytes participated in CaCO3 crystal formation, transportation and shell regeneration in the pearl oyster. These results are helpful to further understand the exact mechanism of hemocyte-mediated biomineralization in shelled molluscs.