Identification of genes potentially related to biomineralization and immunity by transcriptome analysis of pearl sac in pearl oyster Pinctada martensii

Phosphoproteomics changes due to allograft-induced stress responses of Pinctada fucata martensii

Protein phosphorylation modifications are post-translational modifications (PTMs) that play important roles in signal transduction and immune regulation. Implanting a spherical nucleus into a recipient shellfish is critical in marine pearl aquaculture. Protein phosphorylation may be important in the immune responses of Pinctada fucata martensii after nucleus implantation, but their involvement in regulation remains unclear. Here, phosphoproteomics of P. f. martensii gill tissues was conducted 12 h after nuclear implantation using label-free data-independent acquisition (DIA) with LC-MS/MS. Among the 4024 phosphorylated peptides with quantitative information, 181 were up-regulated and 148 were down-regulated. Functional enrichment analysis of these differentially expressed phosphorylated proteins (DEPPs) revealed significant enrichment in functions related to membrane trafficking, exosomes, cytoskeleton, and signal transduction mechanisms. Further, 16 conserved motifs were identified among the DEPPs, including the RSphP, SphP, RSphA, RSphE, PTphP, and ATphP motifs that were significantly conserved, and which may be related to specific kinase recognition. Parallel response monitoring (PRM) analysis validated the abundances of 12 DEPPs from the proteomics, indicating that the phosphoproteomics analyses were robust. 12 DEPPs were selected from the proteomics results through Quantitative real-time PCR (qPCR) technology, and verification analysis was conducted at the gene level. The study suggests that kinases such as MAPKs, Akt, and CK2 may regulate the phosphorylation of related proteins following nuclear implantation. Furthermore, the important signaling pathways of Rap 1, IL-17A, and NF-κB, which are influenced by phosphorylated or dephosphorylated proteins, are found to be involved in this response. Overall, this study revealed the protein phosphorylation responses after nucleus implantation in P. f. martensii, helping to elucidate the characteristics and mechanisms of immune regulation responses in P. f. martensii, in addition to promoting a further understanding of protein phosphorylation modification functions in P. f. martensii.

Whole-genome bisulfite sequencing reveals the function of DNA methylation in the allotransplantation immunity of pearl oysters

Introduction

In the pearl culture industry, a major challenge is the overactive immunological response in pearl oysters resulting from allotransplantation, leading to shell-bead rejection and death. To better understand the molecular mechanisms of postoperative recovery and the regulatory role of DNA methylation in gene expression, we analyzed the changes in DNA methylation levels after allotransplantation in pearl oyster Pinctada fucata martensii, and elucidated the regulatory function of DNA methylation in promoter activity of nicotinic acetylcholine receptor (nAChR) gene.

Methods

We constructed nine DNA methylomes at different time points after allotransplantation and used bisulfite genomic sequencing PCR technology (BSP) to verify the methylation status in the promoter of nAChR. We performed Dual luciferase assays to determine the effect of the dense methylation region in the promoter on transcriptional activity and used DNA pull-down and mass spectrometry analysis to assess the capability of transcription factor binding with the dense methylation region.

Result

The DNA methylomes reveal that CG-type methylation is predominant, with a trend opposite to non-CG-type methylation. Promoters, particularly CpG island-rich regions, were less frequently methylated than gene function elements. We identified 5,679 to 7,945 differentially methylated genes (DMGs) in the gene body, and 2,146 to 3,385 DMGs in the promoter at each time point compared to the pre-grafting group. Gene ontology and pathway enrichment analyses showed that these DMGs were mainly associated with "cellular process", "Membrane", "Epstein-Barr virus infection", "Notch signaling pathway", "Fanconi anemia pathway", and "Nucleotide excision repair". Our study also found that the DNA methylation patterns of the promoter region of nAChR gene were consistent with the DNA methylomics data. We further demonstrated that the dense methylation region in the promoter of nAChR affects transcriptional activity, and that the methylation status in the promoter modulates the binding of different transcription factors, particularly transcriptional repressors.

Conclusion

These findings enhance our understanding of the immune response and regulation mechanism induced by DNA methylation in pearl oysters after allotransplantation.

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 Skeleton and Biomineralization Mechanism as Part of the Innate Immune System of Stony Corals

Stony corals are among the most important calcifiers in the marine ecosystem as they form the coral reefs. Coral reefs have huge ecological importance as they constitute the most diverse marine ecosystem, providing a home to roughly a quarter of all marine species. In recent years, many studies have shed light on the mechanisms underlying the biomineralization processes in corals, as characterizing the calicoblast cell layer and genes involved in the formation of the calcium carbonate skeleton. In addition, considerable advancements have been made in the research field of coral immunity as characterizing genes involved in the immune response to pathogens and stressors, and the revealing of specialized immune cells, including their gene expression profile and phagocytosis capabilities. Yet, these two fields of corals research have never been integrated. Here, we discuss how the coral skeleton plays a role as the first line of defense. We integrate the knowledge from both fields and highlight genes and proteins that are related to biomineralization and might be involved in the innate immune response and help the coral deal with pathogens that penetrate its skeleton. In many organisms, the immune system has been tied to calcification. In humans, immune factors enhance ectopic calcification which causes severe diseases. Further investigation of coral immune genes which are involved in skeleton defense as well as in biomineralization might shed light on our understanding of the correlation and the interaction of both processes as well as reveal novel comprehension of how immune factors enhance calcification.

Transcriptome profiling based on protein-protein interaction networks provides a set of core genes for understanding the immune response mechanisms of the egg-protecting behavior in Octopus ocellatus

Protection via of the immune system is indispensable to the life of organisms. Within an immune network, problems with a given link will affect the normal life activities of the organism. Octopus ocellatus is cephalopod widely distributed throughout the world's oceans. Because of its unique nervous system and locomotive organs, research on this species has gradually increased in recent years. Many immune response mechanisms associated with behaviors of O. ocellatus are still unclear. Moreover, as a factor affecting the normal growth of O. ocellatus, egg protection has rarely been considered in previous behavioral studies. In this study, we analyzed the transcriptome profile of gene expression in O. ocellatus larvae, and identified 5936 differentially expressed genes (DEGs). GO and KEGG enrichment analyses were used to search for immune-related DEGs. Protein-protein interaction networks were constructed to examine the interactions between immune-related genes. Fifteen hub genes involved in multiple KEGG signaling pathways or with multiple protein-protein interaction relationships were obtained and verified by quantitative RT-PCR. We first studied the effects of egg protection on the immunity of O. ocellatus larvae by means of protein-protein interaction networks, and the results provide valuable genetic resources for understanding the immunity of invertebrate larvae. The data serve as a foundation for further research on the egg-protecting behavior of invertebrates.

The role of Smad1/5 in mantle immunity of the pearl oyster Pinctada fucata martensii

The Smad protein family is an important medium for transducing BMP-Smads signals, and which have been proved that their important role in regulating shell biomineralization in Pinctada fucata martensii in our previous study. The members of TGF-β superfamily were involved in innate immunity in vertebrates and invertebrates, and Smad regulatory networks construct a balanced immune system. However, little is known about the role of Smad1/5 in immunity in P. f. martensii. The present study shows that the tissue distribution and the expression profiles of Smad1/5 at developmental stages suggested its wide distribution and crucial role in development at embryonic stages other than larval stage; the increased expression of bone morphogenetic proteins 2 (BMP2), Smad4, Smad1/5 and MSX mRNAs at mantle tissue after LPS and Poly (I:C) challenged implied the potential immune role of Smad1/5 and BMP2-Smad signals to defense against bacterial and virus infections; the reduced expression of immune gene nuclear factor kappa-B (NF-κB), matrix metalloproteinase (MMP), interleukin 17 (IL-17), CuZn-superoxide dismutase (CuZn-SOD), tissue inhibitors of metalloproteinase (TIMP) and lipopolysaccharide-induced TNF-α factor (LITAF) mRNA following knockdown of Smad1/5 indicated that Smad1/5 can regulate their expression via BMP2-Smads pathway in the immunity process; the up-regulated expression of Smad1/5 and BMP2-Smad signals genes, and immune genes during wound healing indicated that Smad1/5 and BMP2-Smad signals genes may be involved in wound healing collaborated with immune genes via a different and complex Smads signaling pathway. These results indicated Smad1/5 could regulate innate immunity via BMP2-Smads signal pathway, and which provided new insights into the relationship between BMP2-Smads signal pathway and mantle immunity.

Identification of novel adhesive proteins in pearl oyster by proteomic and bioinformatic analysis

Byssuses, which are proteinaceous fibers secreted by mollusks, are remarkable underwater adhesives. Although mussel adhesives are well known, much less is known about the byssal proteins of pearl oysters especially in the adhesive regions. In this study, adhesive proteins from the pearl oyster Pinctada fucata were studied in depth by transcriptomics and proteomics approaches. In total, 16 novel proteins were identified including a von Willebrand factor type A domain-containing protein, a thrombospondin-1-like protein, tyrosinase, mucin-like proteins, protease inhibitors, and Pinctada unannotated foot protein 3 (PUF3) to PUF6. Interestingly, PUF3-6 are enriched with glycine, serine, and PXG (X = F/Y/W/K/L) motifs and are highly expressed in the foot. The identification of byssal proteins of the pearl oyster is a key step for understanding byssus formation and may inspire the synthesis of novel adhesives for underwater use and the development of anti-biofouling strategies.

Molecular cloning, characterization, and expression of two TNFRs from the pearl oyster Pinctada fucata martensii

Proteins in the tumor necrosis factor receptor (TNFR) superfamily play significant roles in many physiological and pathological events, such as inflammation, apoptosis, autoimmunity, and organogenesis. Here, two TNFR gene homologs (PmTNFR1 and PmTNFR5) were identified in the pearl oyster Pinctada fucata martensii. The predicted PmTNFR1 and PmTNFR5 protein sequences were 406 and 533 amino acids long, respectively, and both possessed motifs characteristic of the TNFR family, including a TNFR homology domain (CRD), a transmembrane domain (TM), and death domains. However, the predicted amino acid sequences of PmTNFR1 and PmTNFR5 had low identity (~16-23%) with sequences of vertebrate TNFR family proteins. Furthermore, PmTNFR5 had a death domain at the C-terminal, indicating that this protein may be a novel member of the TNFR superfamily. Constitutive PmTNFR1 and PmTNFR5 mRNA expression was detected in all six pearl oyster tissues tested, with comparatively greater transcript abundance in the hepatopancreas and gill. The gene expression levels of PmTNFR1 and PmTNFR5, as well as those of downstream signaling molecules related to the NF-κB pathway (RIP, TRAF2, TRAF3, IKK, and NF-κB), were quantified in the gill after LPS challenge and in the hemocytes after nucleus insertion surgery using real-time PCR (qRT-PCR). We found that all genes were significantly upregulated at 6 h and 12 h post-injection, as well as at 15 d post-insertion. We used RNAi to inhibit the expression of the PmTNFR1 and PmTNFR5 genes. We then quantified the expression levels of PmTNFR1 and PmTNFR5, as well as downstream genes, using qRT-PCR. We found that RNAi inhibition of PmTNFR1 and PmTNFR5 downregulated the downstream genes (RIP, TRAF2, TRAF3, IKK, and NF-κB). Therefore, our results suggested that PmTNFR1 and PmTNFR5 mediate the NF-κB signaling pathway, and are closely related to immune defense, particularly allograft immunity, in the pearl oyster P. fucata martensii.

SEM observation of early shell formation and expression of biomineralization-related genes during larval development in the pearl oyster Pinctada fucata

Shell formation of Pinctada fucata in larval development stages plays a crucial role in their survival. Scanning electron microscopy (SEM) was used to observe the morphological changes during larval development. We found that the early shell forms soon after enlargement of the blastopore at the anterior end of the trochophore stage and the complete shell forms in the spats stage, required for metamorphosis of P. fucata. Based on our transcriptome data of trochophore, D-shaped, umbonal, eyespots and spats stages, including the whole process of shell formation, 93 differentially expressed biomineralization-related genes were identified, of which 25 genes were unique to P. fucata, 30 were identical to genes in pacific oyster, and the remaining genes were annotated to other species. Two-dimensional and three-dimensional principal components analysis (PCA) showed that different developmental stages were significantly different, with the early two stages exhibiting a larger difference compared with the next stages. The 93 genes were sorted into 20 trends with three trends being significantly enriched: an initial increase and then a decrease, a monotonic decrease, and a monotonic increase. Gene expression patterns changed with regulatory function during shell formation. Almost all the biomineralization-related genes were up-regulated in the D-shaped stage, but only five genes were up-regulated in that stage but down-regulated in the remaining stages. There were also 11 genes up-regulated in the last three stages, and a total of 24 genes showed high expression level during the last four stages. The 55 genes selected for shell incision experiment sorted into five trends and most genes presented differences in expression between 24 h and other time points. Considering all these results, there is a correlation with the morphological change and the expression of biomineralization-related genes during larval developmental stages, especially of differently expressed genes.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Functional annotation and identification of MADS-box transcription factors related to tuber dormancy in Helianthus tuberosus L

Dormancy-associated MADS-box (DAM) genes play an important role in plant dormancy and release phases. Little is known about the dormancy characteristics of Jerusalem artichoke tubers. Using bioinformatics, we identified and annotated 23 MADS-box gene sequences from the genome of the Jerusalem artichoke and we analyzed the differential expression of these genes at different developmental stages of tuber dormancy. The results show that all 23 genes encode basic proteins and most of the genes of the same subgroup have similar pI values. MADS-box genes from the Jerusalem artichoke and from other closely related species were divided into ten categories using phylogenetic analysis software. Based on the amino acid sequence of the MADS-domain proteins, ten highly conserved motifs were identified. Gene ontology annotation, InterProScan protein function prediction, and RT-PCR analysis showed that ten MADS-box genes play important roles in the dormancy process of Jerusalem artichoke tubers. Our work lays a foundation for further study of the role of MADS-box genes in the dormancy of the Jerusalem artichoke and other tuber crops.

Histological and Expression Differences Among Different Mantle Regions of the Yesso Scallop (Patinopecten yessoensis) Provide Insights into the Molecular Mechanisms of Biomineralization and Pigmentation

The molecular mechanisms of shell formation and pigmentation are issues of great interest in molluscan studies due to the unique physical and biological properties of shells. The Yesso scallop, Patinopecten yessoensis, is one of the most important maricultural bivalves in Asian countries, and its shell color shows polymorphism. To gain more information about the underlying mechanisms of shell formation and pigmentation, this study presents the first analyses of histological and transcriptional differences between different mantle regions of the Yesso scallop, which are thought to be responsible for the formation of different shell layers. The results showed major microstructural differences between the edge and central mantles, which were closely associated with their functions. Different biomineralization-related GO functions, which might participate in the formation of different shell layers, were significantly enriched in the different mantle regions, indicating the different molecular functions of the two mantle regions in shell formation. The melanogenesis pathway, which controls melanin biosynthesis, was the most significantly enriched pathway in the DEGs between the two mantle regions, indicating its important role in shell pigmentation. Tyr, the key and rate-limiting gene in melanogenesis, was expressed at a remarkably high level in the central mantle, while the upstream regulatory genes included in melanogenesis were mainly upregulated in the edge mantle, suggesting the different molecular functions of the two mantle regions in shell pigmentation.

Comprehensive transcriptome analysis reveal key molecular events in the pearl oyster after pre-grafting conditioning

Pre-grafting conditioning is a crucial procedure before transplant surgery during pearl production. To investigate the molecular response of the pearl oyster Pinctada fucata martensii to conditioning, we constructed two hemocyte transcriptomes from pearl oysters with and without conditioning. A total of 134,222,686 raw reads were generated and assembled using the reference genome of the pearl oyster. Transcriptome analysis revealed 3,074 differentially expressed genes (DEGs). Gene ontology and pathway enrichment analyses revealed that these DEGs were mainly associated with "microtubule-based process", "regulation of actin cytoskeleton", and "cell cycle". All related genes were over-expressed in pearl oysters after conditioning. Some nucleotide-binding oligomerization domain-like receptors (NLR), toll-like receptor, myd88, proinflammatory cytokine interleukin-17 (IL-17), and apoptosis-related genes were highly expressed in pearl oysters after conditioning, indicating that conditioning induced the immune response of pearl oysters. "Fatty acid biosynthesis" (FA biosynthesis) was included in the enriched terms, and all eight FA synthase genes in this pathway were highly induced after conditioning. Four tandemly duplicated arginine kinase genes (PmAK) were found in the genome of P. f. martensii, gene structure and sequence analysis indicated PmAK genes were more diverse compared with that from human and zebra fish. The four tandemly duplicated PmAKs were highly up-regulated after conditioning. These findings will help to elucidate the responding molecular events after conditioning and explain the high pearl oyster survival rate with conditioning after transplantation, thereby providing useful information in perfecting the conditioning method to improve pearl oyster survival rate after transplantation.

The global effects of PmRunt co-located and co-expressed with a lincRNA lncRunt in pearl oyster Pinctada fucata martensii

Runt related transcription factors as trans-acting elements play critical roles in the developmental control of cell fate, hematopoiesis, bone formation and cancers. In previous study, the homologue of runt related transcription factor PmRunt has been identified from pearl oyster Pinctada fucata martensii and considered to play an important role in nacre formation. In this study, we used the same samples to perform RNA-seq to detect the global effects after the decrease of PmRunt expression. The transcription levels of several nacre shell matrix protein (NSMP) genes were significantly changed and the potential compensatory effect could happen internal gene families. Downregulation of PmRunt could also influence the biosynthesis of NSMPs through affecting amino acid metabolism, translation, protein processing and export. The inhibition of PmRunt also possibly affected the expression of caspases, IAPs and C1qs that related to apoptosis and immune. In addition, PmRunt highly expressed at 12 h and 12 d after transplantation in hemolymph, which was corresponded to transplantation immunity immune response and the morphology of pearl sac, suggested the cross-talk of biomineralization-immune regulation in hemocytes. Furthermore, a lincRNA (LncRunt) that co-located with PmRunt was identified and showed a significantly relative expression with PmRunt, which suggested the potential regulation. Therefore, these findings provided new idea to find the regulation targets of runt-related transcription factors and offers evidence of lncRNAs in potential biomineralization-immune regulation.

Genome and transcriptome analyses providing insight into the immune response of pearl oysters after allograft and xenograft transplantations

The immune response after allograft or xenograft transplantation in the pearl oyster is a major factor that cause its nucleus rejection and death. To determine the mechanism underlying the immune response after allograft and xenograft transplantations in the pearl oyster Pinctada fucata martensii, we constructed two sets of transcriptomes of hemocytes at different times (6 and 12 h; 1, 3, 6, 12, and 30 d) after allograft and xenograft transplantations, in which the xenografted mantle tissue was from Pinctada maxima. The transcriptomic analysis reveals many genes are involved in the immune response to transplantation, such as transient receptor potential cation channel (TRP), calmodulin (CaM), DNA replication-related genes, and sugar and lipid metabolism-related genes. The expression of these identified genes was higher in the host pearl oyster transplanted with xenograft than that by allograft. The histological analysis of the pearl sac also confirmed that many hemocytes were still gathered around the transplanted nucleus, and no pearl sac was formed in the host pearl oysters at 30 d after xenograft transplantation. The genomic analysis indicated that pearl oysters evolved many copies of genes, such as TRP, CaM, and GST, to sense and cope with the immune response after transplantation. "Ribosome" and "Cytosolic DNA-sensing pathway" were specifically induced in the xenograft group, whereas "Notch signaling pathway" specifically responded to the allograft transplantation. These results can improve our understanding of the mechanism underlying the immune response of pearl oysters after allograft and xenograft transplantations.

Integrated application of transcriptomics and metabolomics provides insights into unsynchronized growth in pearl oyster Pinctada fucata martensii

Similar to other marine bivalves, Pinctada fucata martensii presents unsynchronized growth, which is one of the problems farmers currently face. However, the underlying mechanisms have not been studied. In the present study, pearl oyster P. f. martensii from cultured stocks were selected to produce a progeny stock. At 180 days, the stock was sorted by size, and fast- and slow-growing individuals were separately sampled. Then, metabolomic and transcriptomic approaches were applied to assess the metabolic and transcript changes between the fast- and slow-growing P. f. martensii groups and understand the mechanism underlying their unsynchronized growth. In the metabolomics assay, 30 metabolites were considered significantly different metabolites (SDMs) between the fast- and slow-growing groups and pathway analysis indicated that these SDMs were involved in 20 pathways, including glutathione metabolism; sulfur metabolism; valine, leucine, and isoleucine biosynthesis; and tryptophan metabolism. The transcriptome analysis of different growth groups showed 168 differentially expressed genes (DEGs) and pathway enrichment analysis indicated that DEGs were involved in extracellular matrix-receptor interaction, pentose phosphate pathway, aromatic compound degradation. Integrated transcriptome and metabolome analyses showed that fast-growing individuals exhibited higher biomineralization activity than the slow-growing group, which consumed more energy than the fast-growing group in response to environmental stress. Fast-growing group also exhibited higher digestion, anabolic ability, and osmotic regulation ability than the slow-growing group. This study is the first work involving the integrated metabolomic and transcriptomic analyses to identify the key pathways to understand the molecular and metabolic mechanisms underlying unsynchronized bivalve growth.

Evolution and function analysis of interleukin-17 gene from Pinctada fucata martensii

Interleukin-17 (IL-17) is a proinflammatory cytokine that plays an important role in immune responses. In this study, we identified 57 IL-17 genes from the genomes of six marine invertebrates, including Pinctada fucata martensii, Crassostrea gigas, Lottia gigantea, Capitella teleta, Mizuhopecten yessoensis, and Mytilus galloprovincialis. Phylogenetic analysis showed that all invertebrate IL-17 genes were clustered into one group, implying that invertebrate IL-17 evolved from one common ancestral gene. From the extron-intron analysis, we found many intronless IL-17 genes in mollusks, which may be caused by retroposition. Tissue and development transcriptomic analysis showed that the expression of PmIL-17 was tissue and developmental stage-specific. Moreover, we cloned the full length of the IL-17-2 gene from P. f. martensii (PmIL-17-2) and explored its function in the immune response. The full-length cDNA of PmIL-17-2 is 719 bp, containing an open reading frame of 564 bp, a 5' -untranslated region (UTR) of 31 bp, and a 3' -UTR of 124 bp with a 30 bp poly (A) tail. PmIL-17-2 had a strong response to lipopolysaccharide (LPS), indicating that the PmIL-17-2 participates in innate immune responses. In situ hybridization of hemocytes showed that PmIL-17-2 was mainly produced by granulosa cells, and the number of the stained granulosa increased after LPS stimulation. These results lay the foundation for the research of IL-17 family in marine invertebrates.

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.

Towards a better understanding of allograft-induced stress response in the pearl oyster Pinctada fucata martensii: Insights from iTRAQ-based comparative proteomic analysis

Implantation of a spherical nucleus into a recipient oyster is a critical step in artificial pearl production. The implanted nucleus is known to trigger cellular stress responses at several levels, yet the molecular mechanism underpinning physiological adaptation of the pearl oysters to nucleus implantation is still poorly understood. In this study, we took advantage of the iTRAQ-based proteomics and LC-MS/MS approach to look into allograft induced gene regulation at the protein expression level in the pearl oyster Pinctada fucata martensii, across a period of 30 days following nucleus implantation. A wide variety of proteins, including a group of immune-related proteins such as E3 ubiquitin-ligase and heat shock proteins, exhibited differential expression in response to the surgical operation. Further comparisons between different sampling points revealed that GO terms including "translation" and "oxidation-reduction process" and KEGG pathways including "glycolysis/gluconeogenesis" and "pyruvate metabolism" were significantly enriched at several time points, indicating the important roles of these molecular events in the stress response of pearl oysters to nucleus implantation. In addition, considerable discrepancy between protein expression level and gene transcript abundancy was identified, as only a few genes showed at least 2-fold expression changes at both proteomic and transcriptomic levels. The result implies that post-transcriptional gene regulation for the key proteins may represent an important aspect of allograft-induced stress response in the pearl oysters. Taken together, the data obtained would contribute to a deeper understanding of the molecular mechanisms enabling stress adaptation of the pearl oysters in response to nucleus implantation.

Characterization of the Tibet plateau Jerusalem artichoke (Helianthus tuberosus L.) transcriptome by de novo assembly to discover genes associated with fructan synthesis and SSR analysis

BACKGROUND

Jerusalem artichoke (Helianthus tuberosus L.) is a characteristic crop in the Qinghai-Tibet Plateau which has rapidly developed and gained socioeconomic importance in recent years. Fructans are abundant in tubers and represent the foundation for their formation, processing and utilization of yield; and are also widely used in new sugar-based materials, bioenergy processing, ecological management, and functional feed. To identify key genes in the metabolic pathway of fructans in Jerusalem artichoke, high-throughput sequencing was performed using Illumina Hi Seq™ 2500 equipment to construct a transcriptome library.

RESULTS

Qinghai-Tibet Plateau Jerusalem artichoke "Qingyu No.1" was used as the material; roots, stems, leaves, flowers and tubers of Jerusalem artichoke in its flowering stage were mixed into a mosaic of the Jerusalem artichoke transcriptome library, obtaining 63,089 unigenes with an average length of 713.6 bp. Gene annotation through the Nr, Swiss Prot, GO, KOG and KEGG databases revealed 34.95 and 46.91% of these unigenes had similar sequences in the Nr and Swiss Prot databases. The GO classification showed the Jerusalem artichoke unigenes were divided into three ontologies, with a total of 49 functional groups encompassing biological processes, cellular components, and molecular functions. Among them, there were more unigenes involved in the functional groups for cellular processes, metabolic processes, and single-organism processes. 38,999 unigenes were annotated by KOG and divided into 25 categories according to their functions; the most common annotation being general function prediction. A total of 13,878 unigenes (22%) were annotated in the KEGG database, with the largest proportion corresponding to pathways related to carbohydrate metabolism. A total of 12 unigenes were involved in the synthesis and degradation of fructan. Cluster analysis revealed the candidate 12 unigene proteins were dispersed in the 5 major families of proteins involved in fructan synthesis and degradation. The synergistic effect of INV gene is necessary during fructose synthesis and degradation in Jerusalem artichoke tuber development. The sequencing data from the transcriptome of this species can provide a reliable data basis for the identification and assessment of the expression of the members of the INV gene family.A simple sequence repeat (SSR) loci search was performed on the transcriptome data of Jerusalem artichoke, identifying 6635 eligible SSR loci with a large proportion of dinucleotide and trinucleotide repeats, and the most different motifs were repeated 5 times and 6 times. Dinucleotide and trinucleotide repeat motifs were the most frequent, with AG/CT and ACC/GGT repeat motifs accounting for the highest proportion.

CONCLUSIONS

In this study, a database search of the transcriptome of the Jerusalem artichoke from the Qinghai Tibet Plateau was conducted by high throughput sequencing technology to obtain important transcriptional and SSR loci information. This allowed characterization of the overall expression features of the Jerusalem artichoke transcriptome, identifying the key genes involved in metabolism in this species. In turn, this offers a foundation for further research on the regulatory mechanisms of fructan metabolism in Jerusalem artichoke.

Whole transcriptome sequencing and biomineralization gene architecture associated with cultured pearl quality traits in the pearl oyster, Pinctada margaritifera

Background

Cultured pearls are unique gems produced by living organisms, mainly molluscs of the Pinctada genus, through the biomineralization properties of pearl sac tissue. Improvement of P. margaritifera pearl quality is one of the biggest challenges that Polynesian research has faced to date. To achieve this goal, a better understanding of the complex mechanisms related to nacre and pearl formation is essential and can now be approached through the use of massive parallel sequencing technologies. The aim of this study was to use RNA-seq to compare whole transcriptome expression of pearl sacs that had producing pearls with high and low quality. For this purpose, a comprehensive reference transcriptome of P. margaritifera was built based on multi-tissue sampling (mantle, gonad, whole animal), including different living stages (juvenile, adults) and phenotypes (colour morphotypes, sex).

Results

Strikingly, few genes were found to be up-regulated for high quality pearls (n = 16) compared to the up-regulated genes in low quality pearls (n = 246). Biomineralization genes up-regulated in low quality pearls were specific to prismatic and prism-nacre layers. Alternative splicing was further identified in several key biomineralization genes based on a recent P. margaritifera draft genome.

Conclusion

This study lifts the veil on the multi-level regulation of biomineralization genes associated with pearl quality determination.

Electronic supplementary material

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

Transcriptional changes in the Japanese scallop (Mizuhopecten yessoensis) shellinfested by Polydora provide insights into the molecular mechanism of shell formation and immunomodulation

The Japanese scallop (Mizuhopecten yessoensis) is one of the most important aquaculture species in Asian countries; however, it has suffered severe infection by Polydora in northern China in recent years, causing great economic losses. The Polydora parasitizes the shell of scallops, badly destroying the shell's structure. To investigate the molecular response mechanism of M. yessoensis to Polydora infestion, a comprehensive and niche-targeted cDNA sequence database for diseased scallops was constructed. Additionally, the transcriptional changes in the edge mantle, central mantle and hemocytes, tissues directly related to the disease, were first described in this study. The results showed that genes involved in shell formation and immunomodulation were significantly differentially expressed due to the infestation. Different transcriptional changes existed between the two mantle regions, indicating the different molecular functions likely responsible for the formation of different shell layers. The differential expression of genes for immune recognition, signal transduction and pathogen elimination presented an integrated immune response process in scallops. Moreover, neuromodulation and glycometabolism involved in the regulation process with relevant function significantly enriched. The study provides valuable information for mechanism study of shell formation and immunomodulation in scallops.

Transcriptome analysis of mantle tissues reveals potential biomineralization-related genes in Tectus pyramis Born

The marine mollusk Tectus pyramis is a valuable shellfish primarily distributed in the tropical waters of the South China Sea, as well as in the Indo-Pacific Ocean and areas near the southern portion of the Japanese Peninsula. Despite major economic interest in this mollusk, limited genomic resources are available for this species, which has prevented studies of the molecular mechanism, such as biomineralization. Here, we report the first comprehensive transcript dataset of T. pyramis mantle tissue. From a total of 16,801,141 reads, 173,671 unique transcripts were assembled, which provides new genomic resources for the understanding of biomineralization in T. pyramis. The most abundant unique sequences of the top 30 most highly expressed genes were annotated as shematrin, while other highly expressed genes included glycine-rich protein and shematrin-1. Based on transcriptome annotation and Gene Ontology classification, 130 biomineralization-related genes were found including members of the BMP (bone morphogenetic proteins), calmodulin, perlucin, and shematrin families, as well as mantle genes, nacrein, and MSI60. The results of qPCR showed that 14 of 24 examined genes were highly expressed in the mantle. A phylogenetic tree of BMP, perlucin, shematrin proteins revealed conservation of their structure and functions and indicated that some members participated in biomineralization in T. pyramis. Taken together, the results presented herein will be useful in studies of molecular mechanisms and pathways of biomineralization in T. pyramis, as well as provide new insight into the mechanisms of biomineralization in gastropods.

Transcriptome analysis of growth heterosis in pearl oyster Pinctada fucata martensii

Heterosis improves growth and survival of shellfish species. Although breeders have widely exploited heterosis, its underlying molecular mechanisms remain unclear. In this study, a 2 × 2 complete diallel cross was facilitated between two full-sib families to produce two inbred families (A and D) and their reciprocal hybrid families (B and C) of pearl oyster Pinctada fucata martensii. Growth traits of the four families were compared at the adult stages. Transcriptome analysis was conducted on the four families using an Illumina sequencing platform. The results revealed that the growth traits of the four families significantly varied (P < 0.05). The mid-parent heterosis values of shell length, shell height, shell width, shell weight, and total weight were 12.9%, 14.9%, 18.2%, 17.2%, and 33.2%, respectively. The B- and C-inbred (A and D) triads had 79 and 68 differentially expressed genes (DEGs), respectively, which were dominantly nonadditive, including overdominance, underdominance, and low-parent dominance. Gene ontology term analysis showed that the DEGs in the B- and C-inbred triads were enriched for metabolic process, cellular process cell part, binding, and catalytic activity. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that the DEGs in the B- and C-inbred triads were involved in focal adhesion, the P13K-Akt signaling pathway, the mRNA surveillance pathway, and the focal adhesion pathway. The reliability of the sequencing data was confirmed by real-time polymerase chain reaction analysis of six growth-related genes. The findings of this study provide new insights into heterosis for growth traits and the design of genetic breeding programs for this species.

Comparison of Two Pearl Sacs Formed in the Same Recipient Oyster with Different Genetic Background Involved in Yellow Pigmentation in Pinctada fucata

Color is one of the most important factors determining the commercial value of pearls. Pinctada fucata is a well-known pearl oyster producing high-quality Akoya pearls. Phenotypic variation in amount of yellow pigmentation produces white and yellowish pearls. It has been reported that polymorphism of yellow pigmentation of Akoya pearls is genetically regulated, but the responsible gene(s) has remained unknown. Here, we prepared pearl sac pairs formed in the same recipient oyster but coming from donor oysters that differ in their color. These two pearl sacs produced pearls with different yellowness even in the same recipient oyster. Yellow tone of produced pearls was consistent with shell nacre color of donor oysters from which mantle grafts were prepared, indicating that donor oysters strongly contribute to the yellow coloration of Akoya pearls. We also conducted comparative RNA-seq analysis and retrieved several candidate genes involved in the pearl coloration. Whole gene expression patterns of pair sacs were not grouped by pearl color they produced, but grouped by recipient oysters in which they were grown, suggesting that the number of genes involved in the yellow coloration is quite small, and that recipient oyster affects gene expression of the majority of genes in the pearl sac.

Immune-Related Functional Differential Gene Expression in Koi Carp (Cyprinus carpio) after Challenge with Aeromonas sobria

In order to understand the molecular basis underlying the host immune response of koi carp (Cyprinus carpio), Illumina HiSeq^TM 2000 is used to analyze the muscle and spleen transcriptome of koi carp infected with Aeromonas sobria (A. sobria). De novo assembly of paired-end reads yielded 69,480 unigenes, of which the total length, average length, N50, and GC content are 70,120,028 bp, 1037 bp, 1793 bp, and 45.77%, respectively. Annotation is performed by comparison against various databases, yielding 42,229 (non-redundant protein sequence (NR): 60.78%), 59,255 (non-redundant nucleotide (NT): 85.28%), 35,900 (Swiss-Prot: 51.67%), 11,772 (clusters of orthologous groups (COG): 16.94%), 33,057 (Kyoto Encyclopedia of Genes and Genomes (KEGG): 47.58%), 18,764 (Gene Ontology (GO): 27.01%), and 32,085 (Interpro: 46.18%) unigenes. Comparative analysis of the expression profiles between bacterial challenge fish and control fish identifies 7749 differentially expressed genes (DEGs) from the muscle and 7846 DEGs from the spleen. These DEGs are further categorized with KEGG. Enrichment analysis of the DEGs and unigenes reveals major immune-related functions, including up-regulation of genes related with Toll-like receptor signaling, complement and coagulation cascades, and antigen processing and presentation. The results from RNA-Seq data are also validated and confirmed the consistency of the expression levels of seven immune-related genes after 24 h post infection with qPCR. Microsatellites (11,534), including di-to hexa nucleotide repeat motifs, are also identified. Altogether, this work provides valuable insights into the underlying immune mechanisms elicited during bacterial infection in koi carp that may aid in the future development of disease control measures in protection against A. sobria.

Sex determination and differentiation genes in a functional hermaphrodite scallop, Nodipecten subnodosus

The lion-paw, Nodipecten subnodosus is one of three scallop species commercially exploited on the west coast of the Peninsula of Baja California. Because nothing is known about sex determination and sexual differentiation in hermaphrodite scallops, in the present work, a global transcriptomic analysis was performed in two early developmental stages, settling eyed-larvae and spat, as well as in three tissues (undifferentiated gonad, digestive gland, and adductor muscle). Over 27 million Illumina paired-end reads were obtained through the MiSeq platform. After processing the reads a total of 243,774 transcripts were assembled with an N50 of 980 and an average length of 775nt. A total of 43,252 proteins were inferred and 36,103 transcripts had at least one homolog in the SwissProt database according to a blastx search. After differential expression analyses and GO annotations it was possible to identify several sex-related genes in the scallop, including one known to be involved in the sex determination pathway of the hermaphrodite model organism Caenorhabditis elegans, N. subnodosus-sex1 (Ns-sex1). Other interesting sex determination and differentiation genes were Ns-dmrta2, Ns-sox9, Ns-wnt4, Ns-doa, Ns-ovo, Ns-vir, among others. Most of these genes were mainly expressed in the testis region, suggesting their participation in male gonad region sex differentiation. These results represent the first available information on the genetics of sex determination and differentiation in a functional hermaphrodite scallop.

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

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

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.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

Identification and Differential Expression of Biomineralization Genes in the Mantle of Pearl Oyster Pinctada fucata

A series of proteins are involved in shell formation of the pearl oyster Pinctada fucata, but the involved mechanisms and the relative expression levels of these proteins have not been elucidated. In this study, we sequenced and characterized the transcriptome of P. fucata mantle tissue. A total of 100,679 unique transcripts were assembled, 43687 Unigenes were annotated, and 48654 CDSs were determined. Of these, GO annotated 16353 Unigenes, COG defined 11585 unigenes into 25 categories, and KEGG sorted 25053 unigenes into 258 pathways. In total, 67 biomineralization-related genes were identified, of which 23 genes were newly described in P. fucata. These genes included ones that expressed shell matrix proteins, regulatory factors, and uncharacterized genes. Differential expression of these 67 genes and 9 other biomineralization-related genes was confirmed using qPCR. Of the 8 nacreous layer-related genes, MSI60 (774.00) was expressed at a much higher level than the others. KRMP2-4 and MSI31 were the most highly expressed of the 13 prismatic layer-related genes and KRMP2 was expressed at nearly 10000 times of the level of the 18S gene. For genes related to both layers, shematrin 2 (3977.84), nacrein (2404.75), PFMG 10 (2113.93), and PFMG 4 (1015.89) were highly expressed, and ferritin-like protein (877.54) and PFMG 8 (516.48) were highly expressed among the 16 undefined genes. The expression levels of regulation factors were generally low, and the highest level was 324.09 (EF-hand) and the lowest occurred in the BMP and wnt families. The expression levels of the prismatic matrix proteins were much higher than those of nacreous ones, consistent with a thicker prismatic layer. MSI60 and nacrein are likely the main components of the nacreous layer, and KRMP2-4, MSI31, shematrin 2, and PFMG 10 gene products are the main components of the prismatic layer. This is the first report of transient expression levels of a large number of biomineralization-related genes at the same time in mantle tissue of P. fucata. These findings provide a novel perspective to understand the molecular mechanisms of shell formation and will be beneficial to genetic improvement of P. fucata for the production of high-quality pearls as well.

Differential mantle transcriptomics and characterization of growth-related genes in the diploid and triploid pearl oyster Pinctada fucata

To explore the molecular mechanism of triploidy effect in the pearl oyster Pinctada fucata, two RNA-seq libraries were constructed from the mantle tissue of diploids and triploids by Roche-454 massive parallel pyrosequencing. The identification of differential expressed genes (DEGs) between diploid and triploid may reveal the molecular mechanism of triploidy effect. In this study, 230 down-regulated and 259 up-regulated DEGs were obtained by comparison between diploid and triploid libraries. The gene ontology and KEGG pathway analysis revealed more functional activation in triploids and it may due to the duplicated gene expression in transcriptional level during whole genome duplication (WGD). To confirm the sequencing data, a set of 11 up-regulated genes related to growth and development control and regulation were analyzed by RT-qPCR in independent experiment. According to the validation and annotation of these genes, it is hypothesized that the set of up-regulated expressed genes had the correlated expression pattern involved in shell building or other interactive probable functions during triploidization. The up- regulation of growth-related genes may support the classic hypotheses of 'energy redistribution' from early research. The results provide valuable resources to understand the molecular mechanism of triploidy effect in both shell building and producing high-quality seawater pearls.

Characterization of transcriptome and identification of biomineralization genes in winged pearl oyster (Pteria penguin) mantle tissue

The winged pearl oyster Pteria penguin is a commercially important marine pearl oyster species, with pearls that are quite different from those of other pearl oysters. Among such species, mantle tissue is the main organ responsible for shell and pearl formation, a biomineralization process that is regulated by a series of genes, most of which remain unknown. In this study, we sequenced and characterized the transcriptome of P. penguin mantle tissue using the HiSeq 2000 sequencing platform. A total of 93,204 unique transcripts were assembled from 51,580,076 quality reads, with a mean length of 608bp, and 40,974 unigenes were annotated. The sequence data enabled the identification of 79,702 potential single nucleotide polymorphism loci and 4345 putative simple sequence repeat loci. A total of 71 unique transcripts were identified homologous to known biomineralization genes, including mantle gene, nacrein, pearlin, pif, chitinase, and shematrin, of which only 3 were previously reported in P. penguin. qPCR analysis indicated that 10 randomly selected biomineralization genes were much more highly expressed in mantle tissue than in the other tissues. In addition, 30 unique sequences were identified as highly expressed, with FPKM values of >3000, and most of these were biomineralization-related genes, including shematrin family genes, a jacalin-related lectin synthesis gene, calponin-2, and paramyosin. These findings will be useful for future studies of biomineralization in P. penguin, as well as in other Pteria species.

Differential Gene Expression during Larval Metamorphic Development in the Pearl Oyster, Pinctada fucata, Based on Transcriptome Analysis

P. fucata experiences a series of transformations in appearance, from swimming larvae to sessile juveniles, during which significant changes in gene expression likely occur. Thus, P. fucata could be an ideal model in which to study the molecular mechanisms of larval metamorphosis during development in invertebrates. To study the molecular driving force behind metamorphic development in larvae of P. fucata, transcriptomes of five larval stages (trochophore, D-shape, umbonal, eyespots, and spats) were sequenced using an Illumina HiSeq™ 2000 system and assembled and characterized with the transcripts of six tissues. As a result, a total of 174,126 unique transcripts were assembled and 60,999 were annotated. The number of unigenes varied among the five larval stages. Expression profiles were distinctly different between trochophore, D-shape, umbonal, eyespots, and spats larvae. As a result, 29 expression trends were sorted, of which eight were significant. Among others, 80 development-related, differentially expressed unigenes (DEGs) were identified, of which the majority were homeobox-containing genes. Most DEGs occurred among trochophore, D-shaped, and UES (umbonal, eyespots, and spats) larvae as verified by qPCR. Principal component analysis (PCA) also revealed significant differences in expression among trochophore, D-shaped, and UES larvae with ten transcripts identified but no matching annotations.

De Novo Transcriptome Analysis of Differential Functional Gene Expression in Largemouth Bass (Micropterus salmoides) after Challenge with Nocardia seriolae

Largemouth bass (Micropterus salmoides) are common hosts of an epizootic bacterial infection by Nocardia seriolae. We conducted transcriptome profiling of M. salmoides to understand the host immune response to N. seriolae infection, using the Illumina sequencing platform. De novo assembly of paired-end reads yielded 47,881 unigenes, the total length, average length, N50, and GC content of which were 49,734,288, 1038, 1983 bp, and 45.94%, respectively. Annotation was performed by comparison against non-redundant protein sequence (NR), non-redundant nucleotide (NT), Swiss-Prot, Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro databases, yielding 28,964 (NR: 60.49%), 36,686 (NT: 76.62%), 24,830 (Swissprot: 51.86%), 8913 (COG: 18.61%), 20,329 (KEGG: 42.46%), 835 (GO: 1.74%), and 22,194 (Interpro: 46.35%) unigenes. Additionally, 8913 unigenes were classified into 25 Clusters of Orthologous Groups (KOGs) categories, and 20,329 unigenes were assigned to 244 specific signalling pathways. RNA-Seq by Expectation Maximization (RSEM) and PossionDis were used to determine significantly differentially expressed genes (False Discovery Rate (FDR) < 0.05) and we found that 1384 were upregulated genes and 1542 were downregulated genes, and further confirmed their regulations using reverse transcription quantitative PCR (RT-qPCR). Altogether, these results provide information on immune mechanisms induced during bacterial infection in largemouth bass, which may facilitate the prevention of nocardiosis.

Pm-miR-133 hosting in one potential lncRNA regulates RhoA expression in pearl oyster Pinctada martensii

Long non-coding RNAs (LncRNAs) are abundant in the genome of higher forms of eukaryotes and implicated in regulating the diversity of biological processes partly because they host microRNAs (miRNAs), which are repressors of target gene expression. In vertebrates, miR-133 regulates the differentiation and proliferation of cardiac and skeletal muscles. Pinctada martensii miR-133 (pm-miR-133) was identified in our previous research through Solexa deep sequencing. In the present study, the precise sequence of mature pm-miR-133 was validated through miR-RACE. Stem loop qRT-PCR analysis demonstrated that mature pm-miR-133 was constitutively expressed in the adductor muscle, gonad, hepatopancreas, mantle, foot, and gill of P. martensii. Among these tissues, the adductor muscle exhibited the highest pm-miR-133 expression. Target analysis indicated that pm-RhoA was the potential regulatory target of pm-miR-133. Bioinformatics analyses revealed that a potential LncRNA (designated as Lnc133) with a mature pm-miR-133 could generate a hairpin structure that was highly homologous to that of Lottia gigantea. Lnc133 was also highly expressed in the adductor muscle, gill, hepatopancreas, and gonad. Phylogenetic analysis further showed that the miR-133s derived from chordate and achordate were separated into two classes. Therefore, Lnc133 hosting pm-miR-133 could be involved in regulating the cell proliferation of adductor muscles by targeting pm-RhoA.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

cDNA Microarray Analysis Revealing Candidate Biomineralization Genes of the Pearl Oyster, Pinctada fucata martensii

Biomineralization is a common biological phenomenon resulting in strong tissue, such as bone, tooth, and shell. Pinctada fucata martensii is an ideal animal for the study of biomineralization. Here, microarray technique was used to identify biomineralization gene in mantle edge (ME), mantle center (MC), and both ME and MC (ME-MC) for this pearl oyster. Results revealed that 804, 306, and 1127 contigs expressed at least three times higher in ME, MC, and ME-MC as those in other tissues. Blast against non-redundant database showed that 130 contigs (16.17 %), 53 contigs (17.32 %), and 248 contigs (22.01 %) hit reference genes (E ≤ -10), among which 91 contigs, 48 contigs, and 168 contigs could be assigned to 32, 26, and 63 biomineralization genes in tissue of ME, MC, and ME-MC at a threshold of 3 times upregulated expression level. The ratios of biomineralization contigs to homologous contigs were similar at 3 times, 10 times, and 100 times of upregulated expression level in either ME, MC, or ME-MC. Moreover, the ratio of biomineralization contigs was highest in MC. Although mRNA distribution characters were similar to those in other studies for eight biomineralization genes of PFMG3, Pif, nacrein, MSI7, mantle gene 6, Pfty1, prismin, and the shematrin, most biomineralization genes presented different expression profiles from existing reports. These results provided massive fundamental information for further study of biomineralization gene function, and it may be helpful for revealing gene nets of biomineralization and the molecular mechanisms underlining formation of shell and pearl for the oyster.

Interactive Effects of Seawater Acidification and Elevated Temperature on the Transcriptome and Biomineralization in the Pearl Oyster Pinctada fucata

Interactive effects of ocean acidification and ocean warming on marine calcifiers vary among species, but little is known about the underlying mechanisms. The present study investigated the combined effects of seawater acidification and elevated temperature (ambient condition: pH 8.1 × 23 °C, stress conditions: pH 7.8 × 23 °C, pH 8.1 × 28 °C, and pH 7.8 × 28 °C, exposure time: two months) on the transcriptome and biomineralization of the pearl oyster Pinctada fucata, which is an important marine calcifier. Transcriptome analyses indicated that P. fucata implemented a compensatory acid-base mechanism, metabolic depression and positive physiological responses to mitigate the effects of seawater acidification alone. These responses were energy-expensive processes, leading to decreases in the net calcification rate, shell surface calcium and carbon content, and changes in the shell ultrastructure. Elevated temperature (28 °C) within the thermal window of P. fucata did not induce significant enrichment of the sequenced genes and conversely facilitated calcification, which was detected to alleviate the negative effects of seawater acidification on biomineralization and the shell ultrastructure. Overall, this study will help elucidate the mechanisms by which pearl oysters respond to changing seawater conditions and predict the effects of global climate change on pearl aquaculture.

De Novo Assembly and Characterization of Early Embryonic Transcriptome of the Horseshoe Crab Tachypleus tridentatus

The horseshoe crab Tachypleus tridentatus is a unique marine species and a potential model for marine invertebrate. Limited genomic and transcriptional data are currently available to understand the molecular mechanisms underlying the embryonic development of T. tridentatus. Here, we reported for the first time the de novo transcriptome assembly for T. tridentatus at embryonic developmental stage using Illumina RNA-seq platform. Approximate 38 million reads were obtained and further assembled into 133,212 unigenes. Sequence homology analysis against public databases revealed that 33,796 unigenes could be annotated with gene descriptions. Of the annotated unigenes, we identified a number of key components of several conserved metazoan signaling pathways (Hedgehog, Wnt, TGF-beta and Notch pathways) and other important regulatory genes involved in embryonic development. Targeted searching of Pax family genes which play critical roles in the formation of tissue and organ during embryonic development identified a complete set of Pax family genes. Moreover, the full length T. tridentatus Pax1/9a (TtPax1/9a) and Pax1/9b (TtPax1/9b) cDNA sequences were determined based on the transcriptome, demonstrating the immediate application of our database. Using quantitative real time PCR, we analyzed the expression patterns of TtPax1/9a and TtPax1/9b in different tissues of horseshoe crab. Taking advantage of Drosophila model, we further found that TtPax1/9b, but not TtPax1/9a, can partly rescue the Drosophila homolog Poxm dysfunction-caused lethality at the larval stage. Our study provides the embryonic transcriptome of T. tridentatus which could be immediately used for gene discovery and characterization, functional genomics studies in T. tridentatus. This transcriptome database will also facilitate the investigations of molecular mechanisms underlying embryonic development of T. tridentatus and other marine arthropods as well.

Identification of functional genes involved in Cd(2+) response of Chinese surf clam (Mactra chinensis) through transcriptome sequencing

The Chinese surf clam Mactra chinensis is an economically important bivalve species in the coastal waters of Liaoning and Shandong Province, China. In this study, we carried out transcriptome sequencing to develop molecular resources for M. chinensis and conducted an acute test of Cd(2+) stimulation through quantitative real-time PCR (qRT-PCR) to analyze the relative expression of six functional genes. A total of 100,839 transcripts and 56,712 unigenes were obtained from 39.9 million filtered reads and 21,305 unigenes were annotated by hitting against NCBI database. According to the results of qRT-PCR, heat shock protein 22 (Hsp22) and cytochrome P450 (CYP450(2C31)) were inhibited in the low concentration, and induced in the high concentration of Cd(2+); thioredoxin peroxidase (TPx-A) was at normal level in low concentration, but induced in high concentration of Cd(2+); glutathione peroxidase A (GPA), glutathione peroxidase 1 (GPA1) and Mn superoxide dismutase gene (MnSOD) were down-regulated when exposed to any treatment groups. Expression levels of the six functional genes following Cd(2+) exposure indicated that these genes were linked to environmental stress. Moreover, the present work enriched the molecule genetic data of M. chinensis.

Transcriptional response to sulfide in the Echiuran Worm Urechis unicinctus by digital gene expression analysis

Background

Urechis unicinctus, an echiuran worm inhabiting the U-shaped burrows in the coastal mud flats, is an important commercial and ecological invertebrate in Northeast Asian countries, which has potential applications in the study of animal evolution, coastal sediment improvement and marine drug development. Furthermore, the worm can tolerate and utilize well-known toxicant-sulfide. However, knowledge is limited on the molecular mechanism of U. unicinctus responding to sulfide due to deficiency of its genetic information.

Methods

In this study, we performed Illumina sequencing to obtain the first Urechis unicinctus transcriptome data. Sequenced reads were assembled and then annotated using blast searches against Nr, Nt, Swiss-Prot, KEGG and COG. The clean tags from four digital gene expression (DGE) libraries were mapped to the U. unicinctus transcriptome. DGE analysis and functional annotation were then performed to reveal its response to sulfide. The expressions of 12 candidate genes were validated using quantitative real-time PCR. The results of qRT-PCR were regressed against the DGE analysis, with a correlation coefficient and p-value reported for each of them.

Results

Here we first present a draft of U. unicinctus transcriptome using the Illumina HiSeq^TM 2000 platform and 52,093 unique sequences were assembled with the average length of 738 bp and N50 of 1131 bp. About 51.6 % of the transcriptome were functionally annotated based on the databases of Nr, Nt, Swiss-Prot, KEGG and COG. Then based on the transcriptome, the digital gene expression analysis was conducted to examine the transcriptional response to sulfide during 6, 24 and 48 h exposure, and finally 1705, 1181 and 1494 tag-mapped genes were identified as differentially expressed genes in the 6-h, 24-h and 48-h libraries, then were further subjected to pathway analyses.

Conclusions

In the DGE database of U. unicinctus, the alterations in certain known sulfide-related pathways indicate similar changes in response to sulfide. For more than 80 % of the identified pathway members, this is the first report on their association with sulfide stress, among which glycolysis pathway and PIDD involving pathways were unique and discussed in details, and were thought to play important roles in the sulfide tolerance of U. unicinctus. All the results are helpful to explain the mechanism of sulfide tolerance and detoxification.

Electronic supplementary material

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

The AP-1 transcription factor homolog Pf-AP-1 activates transcription of multiple biomineral proteins and potentially participates in Pinctada fucata biomineralization

Activator protein-1 (AP-1) is an important bZIP transcription factor that regulates a series of physiological processes by specifically activating transcription of several genes, and one of its well-chartered functions in mammals is participating in bone mineralization. We isolated and cloned the complete cDNA of a Jun/AP-1 homolog from Pinctada fucata and called it Pf-AP-1. Pf-AP-1 had a highly conserved bZIP region and phosphorylation sites compared with those from mammals. A tissue distribution analysis showed that Pf-AP-1 was ubiquitously expressed in P. fucata and the mRNA level of Pf-AP-1 is extremely high in mantle. Pf-AP-1 expression was positively associated with multiple biomineral proteins in the mantle. The luciferase reporter assay in a mammalian cell line showed that Pf-AP-1 significantly up-regulates the transcriptional activity of the promoters of KRMP, Pearlin, and Prisilkin39. Inhibiting the activity of Pf-AP-1 depressed the expression of multiple matrix proteins. Pf-AP-1 showed a unique expression pattern during shell regeneration and pearl sac development, which was similar to the pattern observed for biomineral proteins. These results suggest that the Pf-AP-1 AP-1 homolog is an important transcription factor that regulates transcription of several biomineral proteins simultaneously and plays a role in P. fucata biomineralization, particularly during pearl and shell formation.