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.

Comparative transcriptomic and proteomic analysis of yellow shell and black shell pearl oysters, Pinctada fucata martensii

BACKGROUND

The pearl oyster Pinctada fucata martensii (Pfu.), widely cultured in the South China Sea, is a precious source of sea pearls and calcifying materials. A yellow shell variety of Pfu. was obtained after years of artificial breeding. To identify differentially expressed genes between yellow shell and normal black shell pearl oysters, we performed transcriptomic sequencing and proteomic analyses using mantle edge tissues.

RESULTS

A total of 56,969 unigenes were obtained from transcriptomic, of which 21,610 were annotated, including 385 annotated significant up-regulated genes and 227 significant down-regulated genes in yellow shell oysters (| log2 (fold change) | ≥2 and false discovery rate < 0.001). Tyrosine metabolism, calcium signalling pathway, phototransduction, melanogenesis pathways and rhodopsin related Gene Ontology (GO) terms were enriched with significant differentially expressed genes (DEGs) in transcriptomic. Proteomic sequencing identified 1769 proteins, of which 51 were significantly differentially expressed in yellow shell oysters. Calmodulin, N66 matrix protein, nacre protein and Kazal-type serine protease inhibitor were up-regulated in yellow shell oysters at both mRNA and protein levels, while glycine-rich protein shematrin-2, mantle gene 4, and sulphide: quinone oxidoreductase were down-regulated at two omics levels. Particularly, calmodulin, nacre protein N16.3, mantle gene 4, sulphide: quinone oxidoreductase, tyrosinase-like protein 3, cytochrome P450 3A were confirmed by quantitative real-time PCR. Yellow shell oysters possessed higher total carotenoid content (TCC) compared than black shell oyster based on spectrophotography.

CONCLUSIONS

The yellow phenotype of pearl oysters, characterised by higher total carotenoids content, may reflect differences in retinal and rhodopsin metabolism, melanogenesis, calcium signalling pathway and biomineralisation. These results provide insights for exploring the relationships between calcium regulation, biomineralisation and yellow shell colour pigmentation.

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.

Colour plasticity in the shells and pearls of animal graft model Pinctada margaritifera assessed by HSV colour quantification

The bivalve Pinctada margaritifera has the capacity to produce the most varied and colourful pearls in the world. Colour expression in the inner shell is under combined genetic and environmental control and is correlated with the colour of pearls produced when the same individual is used as a graft donor. One major limitation when studying colour phenotypes is grader subjectivity, which leads to inconsistent colour qualification and quantification. Through the use of HSV (Hue Saturation Value) colour space, we created an R package named 'ImaginR' to characterise inner shell colour variations in P. margaritifera. Using a machine-learning protocol with a training dataset, ImaginR was able to reassign individual oysters and pearls to predefined human-based phenotype categories. We then tested the package on samples obtained in an experiment testing the effects of donor conditioning depth on the colour of the donor inner shell and colour of the pearls harvested from recipients following grafting and 20 months of culture in situ. These analyses successfully detected donor shell colour modifications due to depth-related plasticity and the maintenance of these modifications through to the harvested pearls. Besides its potential interest for standardization in the pearl industry, this new method is relevant to other research projects using biological models.

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.

Regulation of IL-17 by lncRNA of IRF-2 in the pearl oyster

Long noncoding RNAs (lncRNAs), once thought to be nonfunctional, have recently been shown to participate in the multilevel regulation of transcriptional, posttranscriptional and epigenetic modifications and to play important roles in various biological processes, including immune responses. However, the expression and roles of lncRNAs in invertebrates, especially nonmodel organisms, remain poorly understood. In this study, by comparing a transcriptome to the PfIRF-2 genomic structure, we identified lncIRF-2 in the PfIRF-2 genomic intron. The results of the RNA interference (RNAi) and the nucleus grafting experiments indicated that PfIRF-2 might have a negative regulatory effect on lncIRF-2, and PfIRF-2 and lncIRF-2 may have a positive regulatory effect on PfIL-17. Additionally, lncIRF-2, PfIRF-2 and PfIL-17 were involved in responses to the nucleus graft. These results will enhance the knowledge of lncIRF-2, IRF-2, and IL-17 functions in both pearl oysters and other invertebrates.

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.

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.

Cultured Pearl Surface Quality Profiling by the Shell Matrix Protein Gene Expression in the Biomineralised Pearl Sac Tissue of Pinctada margaritifera

Nucleated pearls are produced by molluscs of the Pinctada genus through the biomineralisation activity of the pearl sac tissue within the recipient oyster. The pearl sac originates from graft tissue taken from the donor oyster mantle and its functioning is crucial in determining key factors that impact pearl quality surface characteristics. The specific role of related gene regulation during gem biogenesis was unknown, so we analysed the expression profiles of eight genes encoding nacreous (PIF, MSI60, PERL1) or prismatic (SHEM5, PRISM, ASP, SHEM9) shell matrix proteins or both (CALC1) in the pearl sac (N = 211) of Pinctada margaritifera during pearl biogenesis. The pearls and pearl sacs analysed were from a uniform experimental graft with sequential harvests at 3, 6 and 9 months post-grafting. Quality traits of the corresponding pearls were recorded: surface defects, surface deposits and overall quality grade. Results showed that (1) the first 3 months of culture seem crucial for pearl quality surface determination and (2) all the genes (SHEM5, PRISM, ASP, SHEM9) encoding proteins related to calcite layer formation were over-expressed in the pearl sacs that produced low pearl surface quality. Multivariate regression tree building clearly identified three genes implicated in pearl surface quality, SHEM9, ASP and PIF. SHEM9 and ASP were clearly implicated in low pearl quality, whereas PIF was implicated in high quality. Results could be used as biomarkers for genetic improvement of P. margaritifera pearl quality and constitute a novel perspective to understanding the molecular mechanism of pearl formation.

Adaptive response of pearl oyster Pinctada fucata martensii to low water temperature stress

The pearl oyster Pinctada fucata martensii is a warm-water shellfish that is sensitive to cold environments. To investigate its potential adaptation to low-temperature stress, the selected line (SL) and based population (BP) were sampled to undergo transcriptome sequence. Results of transcriptome analysis showed 572 significant differentially expressed genes. The typical HSP70 and HSP40 exhibited the polar expression model in the two groups. Meanwhile, the related genes that involved in energy release mediated by oxidative phosphorylation and the biosynthesis of unsaturated fatty acid were increased in the SL. The apparent enrichment of different expressed genes in amino acid metabolism indicated that the small molecule system with amino acids was one of the main regulator for low-temperature stress. The different expressions of immune-related and lysosome protein encoding genes also reflected the variation of immunity in the two groups and indicated that it could affect the adaptation ability in different temperature. In addition, the similar trends of different expression of typical genes between two groups were obtained by using RNA-seq and qRT-PCR. These results suggested that multi-system adjustments are involved in the processes of low water temperature stress in pearl oyster, providing insights into the response systems of shellfish to acclimatise with ambient environment change.

Optimal age of the donor graft tissue in relation to cultured pearl phenotypes in the mollusc, Pinctada margaritifera

Ageing is defined as the progressive decline in tissue and organ functions over time. This study aims to evaluate the ageing effect on cultured pearl quality phenotypes (including size and quality traits) in the graft-recipient animal model: Pinctada margaritifera. For this, eight uniform grafting experiments were designed using two hatchery-produced pearl oyster families as donors, which were followed through time, between 7 and 30 months in age. For each age category, 20 donors were studied for each culture site giving a total of 2400 grafted oysters. Several phenotypic measurements were made: 1) donor family growth performance from shell size records, 2) pearl size and corresponding quality traits, and 3) expression of some genes related to biomineralization processes on both the mantle graft and on pearl sac tissues. Results showed that: 1) donor age has an impact on pearl size, with grafts coming from the youngest donors yielding the biggest pearls; and 2) grafts from donors between 12 and 18 months in age produced pearls of the highest quality (grade and surface quality), a result supported by an analysis where the level of expression for a panel of genes associated with biomineralization was greatest in donors within the 12 to 18 months age group. These results indicate that donors aged between 12 and 18 months have high potential for biomineralisation and nacre deposition, and likely produce larger and higher quality cultured pearls than older donors.

Novel Isoforms of N16 and N19 Families Implicated for the Nacreous Layer Formation in the Pearl Oyster Pinctada fucata

Although a wide variety of proteins and genes possibly related to the shell formation in bivalve have been identified, their functions have been only partially approved. We have recently performed deep sequencing of expressed sequence tags (ESTs) from the pearl oyster Pinctada fucata using a next-generation sequencer, identifying a dozen of novel gene candidates which are possibly associated with the nacreous layer formation. Among the ESTs, we focused on three novel isoforms (N16-6, N16-7, and N19-2) of N16 and N19 families with reference to five known genes in the families and determined the full-length cDNA sequences of these isoforms. Reverse transcription-polymerase chain reaction revealed that N16-6 was expressed in gill, gonad, adductor muscle, and mantle, whereas N16-7 exclusively in mantle. N19-2 was expressed in all tissues examined. In situ hybridization demonstrated their regional expression in mantle and pearl sac, which well corresponded to those shown by EST analysis previously reported. Shells in the pearl oyster injected with dsRNAs of N16-7 and N19-2 showed abnormal surface appearance in the nacreous layer. Taken together, novel isoforms in N16 and N19 families shown in this study are essential to form the nacreous layer.

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.

Preparative expression and purification of a nacreous protein N16 and testing its effect on osteoporosis rat model

N16, a nacreous protein isolated from Pinctada martensii, is related to nacreous layer formation. Our previous study indicated that N16 showed dual regulatory effects by inducing osteoblast biomineralization as well as inhibiting osteoclast formation. In order to obtain large quantity of N16 for animal experiment and clinical trial, a fermentation and preparative purification method was established. The N16 cDNA was cloned to a BL21(DE3)plysE-pET32a vector and grown in a 20 L fermenter. The medium, temperature, pH and dissolved oxygen (DO) were optimized. N16 was expressed in inclusion bodies. It was denatured and refolded in 8 M urea buffer and purified to 97% purity by passing through a gel filtration column. The glucocorticoid induced osteoporosis (GIO) rat model was used to investigate the anti-osteoporosis activity of N16 in vivo. Results showed that the decrease of the bone mineral density (BMD) and the ultimate load was significantly relieved after N16 treatment. N16 displayed dual regulatory effects by promoting osteogenesis as well as inhibiting bone resorption in vivo. Our work will contribute to further clinical studies on N16 for osteoporosis treatment.

First evidence of establishment of the rayed pearl oyster, Pinctada imbricata radiata (Leach, 1814), in the eastern Adriatic Sea

The Mediterranean Sea is increasingly under threat from invasive species that may negatively affect biodiversity and/or modify ecosystem structure and function. The bivalve mollusc Pinctada imbricata radiata is listed among the 100 most invasive species in the Mediterranean. A first finding of an established population of P. imbricata radiata in the coastal waters of the eastern Adriatic Sea, is presented in this paper. Six and then 30 live specimens were collected in 2015 and in 2017, respectively, at depths of 5 to 15m from the island of Mljet, Croatia. DNA sequencing of the mitochondrial cytochrome c oxidase I gene (COI) revealed three different haplotypes. All samples showed greatest similarity (98 to >99%) to P. radiata COI sequence records in GenBank (=P. imbricata radiata as used in this paper). A Neighbour Joining tree placed all Croatian samples within the 100% bootstrap supported clade for P. imbricata radiata.

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.

A Spirochaete is suggested as the causative agent of Akoya oyster disease by metagenomic analysis

Mass mortality that is acompanied by reddish browning of the soft tissues has been occurring in cultured pearl oyster, Pinctada fucata martensii. The disease is called Akoya oyster disease (AOD). Although spreading pattern of the disease and transmission experiments suggest that the disease is infectious, the causative agent has not yet been identified. We used shotgun and 16S rRNA-based metagenomic analysis to identify genes that are present specifically in affected oysters. The genes found only in diseased oysters were mostly bacterial origin, suggesting that the causative agent was a bacterial pathogen. This hypothesis was supported by the inhibition of AOD development in naïve oysters injected with the hemolymph of diseased animals followed immediately with penicillin bath-administration. Further analyses of the hemolymph and mantle specifically and universally detected genes of bacteria that belong to phylum Spirochaetes in diseased pearl oysters but not in healthy oysters. By in situ hybridization or immunostaining, a Brachyspira-like bacterium was observed in the smears of hemolymph from affected oysters, but not from healthy oysters. Phylogenetic analysis using 16S rRNA sequences showed that the presumptive causative bacterium was outside of but most closely related to family Brachyspiraceae. We propose ‘Candidatus Maribrachyspira akoyae’ gen. nov, sp nov., for this bacterium.

Transcriptome analysis of the immune reaction of the pearl oyster Pinctada fucata to xenograft from Pinctada maxima

The pearl oyster Pinctada maxima exhibits great difficulty to culture pearls through nuclear insertion with an allograft, but it is easy for P. fucata to culture pearls after allografting. If P. fucata could be used as a surrogate mother to culture P. maxima pearls, it would benefit the pearl culture industry of P. maxima. However, this is blocked by the immune rejection of P. fucata against P. maxima mantle grafts. In this study, the immune responses of P. fucata hemocyte to allograft and xenograft were investigated after transplantation by transcriptome analysis. In total, 107.93 Gb clean reads were produced and assembled using the reference genome of P. fucata. Gene Ontology Term enrichment and KEGG enrichment analyses indicated that apoptosis, hippo signaling pathway, oxidation-reduction, MAPK signaling pathway, ribosome, protein processing in endoplasmic reticulum, purine metabolism, NF-kappa B signaling pathway, oxidative phosphorylation, Ras signaling pathway, and ubiquitin mediated proteolysis were involved in response to transplantation. Many genes related to oxidation-reduction reactions, the MAPK signaling pathway, and apoptosis were identified by comparison of the allograft group and the xenograft group at 0 h, 6 h, 12 h, 24 h, 48 h, 72 h, and 96 h post-transplantation. Among them, the expression levels of NADH dehydrogenase, succinate dehydrogenase and other dehydrogenases were increased significantly in the xenograft groups compared with allograft groups at 0 h post transplantation, indicating that a respiratory burst of neutrophils occurred immediately after xenograft transplantation. Additionally, HSP70 was highly expressed from 0 h to 96 h in the xenograft groups, indicating an oyster immune response to the xenograft. The genes enriched in the ribosome and hippo-signaling pathways were also identified, and expression patterns of these DEGs were different as compared between transplantation and control groups. Finally, altered expression levels of 10 randomly selected immune-related DEGs were confirmed by quantitative real-time PCR. These results indicated that oxidation-reduction is likely the key factor responsible for immune rejection to transplantation. The findings should provide some new insight into the molecular mechanism of immune rejection of the host against xenograft, and thus benefit to development of immunosuppressive reagents to facilitate effective xenograft pearling.

Characterization of the distal-less homologue gene, PfDlx, involved in regulating the expression of Pif in the pearl oyster, Pinctada fucata

Distal-less (Dlx) homeobox transcription factors play an important role in regulating various aspects of vertebrate biology. In vertebrates and invertebrates, distal-less is a highly conserved and well-studied transcription factor. In pearl oyster, we have identified a homologue of this gene, Dlx, and cloned the full-length cDNA. Genomic structure analysis revealed that PfDlx genomic DNA contained three exons and two introns. Their deduced amino acid sequences all showed the highest identity with homologues in Crassostrea gigas. Analyses of PfDlx mRNA in tissues and developmental stages showed high expressions in gonad, polar body stage, 2-4 cells and 32 cells. After shell notching, the changes in expression of Dlx shows that it reached a maximum at 24h. In co-transfection experiments, PfDlx significantly activates reporter constructs containing a Pif promoter. Through using RNAi techniques, we demonstrated that down-regulation of Dlx in P. fucata did not significantly disrupt the development of the nacreous layer in scanning electron microscopy, but it significantly down-regulated the expression of Pif gene. Thus, our work suggests that PfDlx might participate in regulating the expression of the Pif gene in the pearl oyster.

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.

Response of the pearl oyster Pinctada margaritifera to cadmium and chromium: Identification of molecular biomarkers

This study was designed to identify in the pearl oyster Pinctada margaritifera, used as a bio-accumulator, molecular biomarkers for the presence of heavy metals in the lagoon environment. Pearl oysters were exposed to 2 concentrations (1 and 10μgL^-1) of cadmium (Cd) and chromium (Cr) compared to a control. Twelve target genes encoding proteins potentially involved in the response to heavy metal contamination with antioxidant, detoxification or apoptosis activities were selected. P. margaritifera accumulated Cd but not Cr, and mortality was related to the amount of Cd accumulated in tissues. In response to Cd-Cr contamination, metallothionein (MT) was significantly up-regulated by Cd-Cr at both concentrations, while 7 others (SOD, CAT, GPX, GSTO, GSTM, CASP, MDR) were down-regulated. Based on the development of these molecular tools, we propose that the pearl oyster, P. margaritifera, could be used as a sentinel species for heavy metal contamination in the lagoons of tropical ecosystems.

Co-Option and De Novo Gene Evolution Underlie Molluscan Shell Diversity

Molluscs fabricate shells of incredible diversity and complexity by localized secretions from the dorsal epithelium of the mantle. Although distantly related molluscs express remarkably different secreted gene products, it remains unclear if the evolution of shell structure and pattern is underpinned by the differential co-option of conserved genes or the integration of lineage-specific genes into the mantle regulatory program. To address this, we compare the mantle transcriptomes of 11 bivalves and gastropods of varying relatedness. We find that each species, including four Pinctada (pearl oyster) species that diverged within the last 20 Ma, expresses a unique mantle secretome. Lineage- or species-specific genes comprise a large proportion of each species' mantle secretome. A majority of these secreted proteins have unique domain architectures that include repetitive, low complexity domains (RLCDs), which evolve rapidly, and have a proclivity to expand, contract and rearrange in the genome. There are also a large number of secretome genes expressed in the mantle that arose before the origin of gastropods and bivalves. Each species expresses a unique set of these more ancient genes consistent with their independent co-option into these mantle gene regulatory networks. From this analysis, we infer lineage-specific secretomes underlie shell diversity, and include both rapidly evolving RLCD-containing proteins, and the continual recruitment and loss of both ancient and recently evolved genes into the periphery of the regulatory network controlling gene expression in the mantle epithelium.

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.

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.