Identification and functional characterization of a superoxide dismutase (CuZnSOD) from Pinctada fucata martensii

Copper/zinc superoxide dismutase (Cu/Zn-SOD) can effectively eliminate reactive oxygen species (ROS)，avoid damage from O2 to the body, and maintain O2 balance. In this study, multi-step high-performance liquid chromatography (HPLC), combined with Mass Spectrometry (MS), was used to isolate and identify Cu/Zn-SOD from the serum of Pinctada fucata martensii (P. f. martensii) and was designated as PmECSOD. With a length of 1864 bp and an open reading frame (ORF) of 1422 bp, the cDNA encodes a 473 amino acid protein. The PmECSOD transcript was detected in multiple tissues by quantitative real-time PCR (qRT-PCR), with its highest expression level being in the gills. Additionally, the temporal expression of PmECSOD mRNA in the hemolymph was highest at 48 h after in vivo stimulation with Escherichia coli and Micrococcus luteus. The results from this study provide a valuable base for further exploration of molluscan innate immunity and immune response.

Effects of titanium dioxide nanoparticle exposure on the gut microbiota of pearl oyster (Pinctada fucata martensii)

With the advancement of nanotechnology and the growing utilization of nanomaterials, titanium dioxide (TiO2) has been released into aquatic environments, posing potential ecotoxicological risks to aquatic organisms. In this study, the toxicological effects of TiO2 nanoparticles were investigated on the intestinal health of pearl oyster (Pinctada fucata martensii). The pearl oysters were subjected to a 14-day exposure to 5-mg/L TiO2 nanoparticle, followed by a 7-day recovery period. Subsequently, the intestinal tissues were analyzed using 16S rDNA high-throughput sequencing. The results from LEfSe analysis revealed that TiO2 nanoparticle increased the susceptibility of pearl oysters to potential pathogenic bacteria infections. Additionally, the TiO2 nanoparticles led to alterations in the abundance of microbial communities in the gut of pearl oysters. Notable changes included a decrease in the relative abundance of Phaeobacter and Nautella, and an increase in the Actinobacteria, which could potentially impact the immune function of pearl oysters. The abundance of Firmicutes and Bacteroidetes, as well as the expression of genes related to energy metabolism (AMPK, PK, SCS-1, SCS-2, SCS-3), were down-regulated, suggesting that TiO2 nanoparticles exposure may affect the digestive and energy metabolic functions of pearl oysters. Furthermore, the short-term recovery of seven days did not fully restore these levels to normal. These findings provide crucial insights and serve as an important reference for understanding the toxic effects of TiO2 nanoparticles on bivalves.

A novel birnavirus identified as the causative agent of summer atrophy of pearl oyster (Pinctada fucata (Gould))

The Akoya pearl oyster (Pinctada fucata (Gould)) is the most important species for pearl cultivation in Japan. Mass mortality of 0-year-old juvenile oysters and anomalies in adults, known as summer atrophy, have been observed in major pearl farming areas during the season when seawater temperatures exceed about 20 °C since 2019. In this study, we identified a novel birnavirus as the pathogen of summer atrophy and named it Pinctada birnavirus (PiBV). PiBV was first presumed to be the causative agent when it was detected specifically and frequently in the infected oysters in a comparative metatranscriptomics of experimentally infected and healthy pearl oysters. Subsequently, the symptoms of summer atrophy were reproduced by infection tests using purified PiBV. Infection of juvenile oysters with PiBV resulted in an increase in the PiBV genome followed by the atrophy of soft body and subsequent mortality. Immunostaining with a mouse antiserum against a recombinant PiBV protein showed that the virus antigen was localized mainly in the epithelial cells on the outer surface of the mantle. Although the phylogenetic analysis using maximum likelihood method placed PiBV at the root of the genus Entomobirnavirus, the identity of the bi-segmented, genomic RNA to that of known birnaviruses at the full-length amino acid level was low, suggesting that PiBV forms a new genus. The discovery of PiBV will be the basis for research to control this emerging disease.

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.

Transcriptome and exosome proteome analyses provide insights into the mantle exosome involved in nacre color formation of pearl oyster Pinctada fucata martensii

Color polymorphisms in molluscan shells play an important economic in the aquaculture industry. Among bivalves, shell color diversity can reflect properties such as growth rate and tolerance. In pearl oysters, the nacre color of the donor is closely related to the pearl color. Numerous genes and proteins involved in nacre color formation have been identified within the exosomes of the mantle. In this study, we analyzed the carotenoids present in the mantle of gold- and silver-lipped pearl oysters, identifying capsanthin and xanthophyll as crucial pigments contributing to coloration. Transcriptome analysis of the mantle revealed several differentially expressed genes (DEGs) involved in color formation, including ferric-chelate reductase, mantle genes, and larval shell matrix proteins. We also isolated and identified exosomes from the mantles of both gold- and silver-lipped strains of the pearl oyster Pinctada fucata martensii, revealing the extracellular transition mechanism of coloration-related proteins. From these exosomes, we obtained a total of 1223 proteins, with 126 differentially expressed proteins (DEPs) identified. These proteins include those associated with carotenoid metabolism and Fe(III) metabolism, such as apolipoproteins, scavenger receptor proteins, β,β-carotene-15,15'-dioxygenase, ferritin, and ferritin heavy chains. This study may provide a new perspective on the nacre color formation process and the pathways involved in deposition within the pearl oyster P. f. martensii.

Effect of elevated temperature, sea water acidification, and phenanthrene on the expression of genes involved in the shell and pearl formation of economic pearl oyster (Pinctada radiata)

Our study aims to examine the effect of some stressors on the gene expression levels of shell matrix proteins in a pearl oyster. Oysters were exposed to the different combinations of the temperature, pH, and phenanthrene concentration is currently measured in the Persian Gulf and the predicted ocean warming and acidification for 28 days. The expression of all the studied genes was significantly downregulated. Time and temperature had the greatest effects on the decreases in n19 and n16 genes expression, respectively. Aspein and msi60 genes expression were highly influenced by pH. Pearlin was affected by double interaction temperature and phenanthrene. Moreover, a correlation was observed among the expression levels of studied genes. This study represents basic data on the relationship between mRNA transcription genes involved in the shell and pearl formation and climate changes in pollutant presence conditions and acclimatizing mechanism of the oyster to the future scenario as well.

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.

Toll-like receptor may be involved in acquired immune response in pearl oyster Pinctada fucata

The increasing experimental evidence suggests that there are some forms of specific acquired immunity in invertebrates, in which Toll-like receptors (TLRs) play vital roles in activating innate and adaptive immunity and have been comprehensively investigated in mammalian species. Yet, the immune mechanisms underlying TLR mediation in mollusks remain obscure. In this study, we identified a TLR13 gene in the pearl oyster Pinctada fucata for the first time and named it PfTLR13 which consists of a 5'-untranslated terminal region (5'-UTR) of 543 bp, an open reading frame (ORF) of 2667 bp, and a 3'-UTR of 729 bp. We found that PfTLR13 mRNA was expressed in all tissues examined, with the highest level in the gills. The expression of PfTLR13 in the gills of oysters exposed to Vibrio alginolyticus or pathogen-associated molecular patterns (PAMPs) (including LPS, PGN, and poly(I:C)) was significantly higher than in the control group. Interestingly, the immune response to the first stimulation was weaker than the response to the second stimulation, suggesting that the primary stimulation may lead to immune priming of TLR in pearl oysters, similar to acquired immunity in vertebrates. Furthermore, we found that PfTLR13 expression was differentially associated with allograft and xenograft in the pearl oyster P. fucata, with the highest expression levels observed at 12 h post-allograft and 24 h post-xenograft. Overall, our findings provide new insights into the immune mechanisms underlying TLR mediation in mollusks and suggest that PfTLR13 may play a crucial role in the specific acquired immunity of pearl oysters.

The first reporting of prevalence Vibrio species and expression of HSP genes in rayed pearl oyster (Pinctada radiata) under thermal conditions

The main goal of the present study was to evaluate the influence of thermal exposure on Vibrio population and HSP genes expression (HSP 90, HSP70, and HSP20) in rayed pearl oyster (P. radiata). To this end, the oysters were reared for 30 days at temperatures of 22 °C (control), 25 °C, 27 °C, and 29 °C. The results showed that five dominate Vibrio strains including Vibrio hepatarius, V. harveyi, V. alginolyticus, V. parahaemolyticus, and V. rotiferianus were identified. The highest population of V. parahaemolyticus, V. alginolyticus, and V. harveyi, was found in 29οC group. According to real-time PCR, mantle exhibited the highest expression levels of HSP20, HSP70, and HSP90 genes. A higher level of HSP20 expression was observed at high temperatures (25 °C, 27 °C, and 29 °C) in the gonad and mantle compared to the control group (22 °C) while decrease in HSP90 expression level was recorded in 25 °C, 27 °C, and 29 °C groups. HSP20 expression level in adductor muscle was remarkably down-regulated in 27 °C and 29 °C groups. In this tissue, HSP70 was detected at highest levels in the 29οC group. In mantle, HSP90 gene expression was lowest at 22 °C water temperature. Several Vibrio strains have been identified from pearl Gulf oyster that haven't been previously reported. The identification of dominant Vibrio species is essential for epidemiological management strategies to control and prevent Vibrio outbreaks in pearl oyster farms. The expression pattern of HSP genes differs in rayed pearl oyster tissues due to differences in their thermal tolerance capability and physiological and biological characteristics. The present study provides useful molecular information for the ecological adaptation of rayed pearl oysters after exposure to different temperature levels.

Purification and Characterization of the Lecithin-Dependent Thermolabile Hemolysin Vhe1 from the Vibrio sp. Strain MA3 Associated with Mass Mortality of Pearl Oyster (Pinctada fucata)

In a previous study, we isolated a Vibrio sp. strain MA3 and its virulence factor, a hemolysin encoded by vhe1. This strain is associated with mass mortalities of the pearl oyster Pinctada fucata. In the present study, the vhe1 gene from strain MA3 was cloned and its encoded product was purified and characterized. Our results show that the vhe1 gene encodes a protein of 417 amino acids with an estimated molecular mass of 47.2 kDa and a pI of 5.14. The deduced protein, Vhe1, was found to contain the conserved amino acid sequence (GDSL motif) of the hydrolase/esterase superfamily and five conserved blocks characteristic of SGNH hydrolases. A BLAST homology search indicated that Vhe1 belongs the lecithin-dependent hemolysin/thermolabile hemolysin (LDH/TLH) family. In activity analyses, the optimal temperature for both the hemolytic and phospholipase activities of Vhe1 was 50 °C. Vhe1 hemolytic activity and phospholipase activity were highest at pH 8.5 and pH 8.0, respectively. However, both enzymatic activities sharply decreased at high temperature (> 50 °C) and pH < 7.0. Compared with previously reported hemolysins, Vhe1 appeared to be more thermal- and pH-labile. Both its hemolytic activity and phospholipase activity were significantly inhibited by CuCl2, CdCl2, ZnCl2, and NiCl2, and slightly inhibited by MnCl2 and CoCl2. Vhe1 showed higher phospholipase activity toward medium-chain fatty acids (C8-C12) than toward shorter- and longer-chain fatty acids. These results accumulate knowledge about the LDH/TLH of V. alginolyticus, which detailed characterization has not been reported, and contribute to solving of the mass mortality of pearl oyster.

Integrated transcriptomic and metabolomic analysis sheds new light on adaptation of Pinctada fucata martensii to short-term hypoxic stress

Analyses of the transcriptome and metabolome were conducted to clarify alterations of key genes and metabolites in pearl oysters following exposure to short-term hypoxic treatment. We totally detected 209 DEGs between the control and hypoxia groups. Enrichment analysis indicated the enrichment of GO terms including "oxidation-reduction process", "ECM organization", "chaperone cofactor-dependent protein refolding", and "ECM-receptor interaction" KEGG pathway by the DEGs. In addition, between the two groups, a total of 28 SDMs were identified, which were implicated in 13 metabolic pathways, such as "phenylalanine metabolism", "D-amino acid metabolism", and "aminoacyl-tRNA biosynthesis". Results suggest that pearl oysters are exposed to oxidative stress and apoptosis under short-term hypoxia. Also, pearl oysters might adapt to short-term hypoxic treatment by increasing antioxidant activity, modulating immune and biomineralization activities, maintaining protein homeostasis, and reorganizing the cytoskeleton. The results of our study help unveil the mechanisms by which pearl oysters respond adaptively to short-term hypoxia.

A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

Comprehensive analysis of the lncRNAs, mRNAs, and miRNAs implicated in the immune response of Pinctada fucata martensii to Vibrio parahaemolyticus

Non-coding RNAs (ncRNAs) have been implicated in a variety of biological processes. However, most ncRNAs are of unknown function and are as-yet unannotated. The immune-related functions of ncRNAs in the pearl oyster Pinctada fucata martensii were explored based on transcriptomic differences in the expression levels of long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in the hemocytes of P.f. martensii after challenge by the pathogenic bacterium Vibrio parahaemolyticus. Across the challenged and control pearl oysters, 144 miRNAs and 14,571 lncRNAs were identified. In total, 13,375 ncRNAs were differentially expressed between the challenged and control pearl oysters; in the challenged pearl oysters as compared to the controls, 15 miRNAs and 5147 lncRNAs were upregulated, while 51 miRNAs and 8162 lncRNAs were downregulated. The sequencing results were validated using quantitative real-time polymerase chain reaction (qRT-PCR) analysis. GO and KEGG pathway analysis showed that genes targeted by the differentially expressed ncRNAs were associated with the vascular endothelial growth factor (VEGF) signaling pathway and the nuclear factor kappa-B (NF-κB) signaling pathway. An lncRNA-mRNA-miRNA network that was developed based on the transcriptomic results of this study suggested that lncRNAs may compete with miRNAs for mRNA binding sites. This study may provide a useful framework for the detection of additional novel ncRNAs, as well as new insights into the pathogenic mechanisms underlying the response of P.f. martensii to V. parahaemolyticus.

Lysine acetylation plays a role in the allograft-induced stress response of the pearl oyster (Pinctada fucata martensii)

Implanting a spherical nucleus into a recipient oyster is a critical step in artificial pearl production using the pearl oyster Pinctada fucata martensii. However, little is known about the role of post-translational modifications (PTMs) in the response of the pearl oyster to this operation. Lysine acetylation, a highly conserved PTM, may be an essential adaptive strategy to manage multiple biotic or abiotic stresses. We conducted the first lysine acetylome analysis of the P. f. martensii gill 12 h after nucleus implantation, using tandem mass tags (TMT) labeling and K^ac affinity enrichment. We identified 2443 acetylated sites in 1301 proteins, and 1511 sites on 895 proteins were quantitatively informative. We found 25 conserved motifs from all of the identified lysine sites, particularly motifs K^ac H, K^ac S, and K^ac Y were strikingly conserved, of which K^ac Y, K^ac H, Y K^ac, K^ac K, K^ac *K, K^ac R, and K^ac F which have been observed in other species and are therefore highly conserved. We identified 58 sites that were significantly differently acetylated in P. f. martensii in response to allograft (|fold change|>1.2, P-value ≤ 0.05); 38 newly acetylated and 20 deacetylated. According to GO functional analysis, subcellar location, and KOG classIfication, these proteins were divided into four categories: cytoskeleton, response to stimulus, metabolism, and other. The differentially acetylated proteins (DAPs) enriched pathways include aminoacyl-tRNA biosynthesis, salmonella infection, and longevity regulating pathway-worm-Caenorhabditis elegans (nematode). Parallel reaction-monitoring (PRM) validation of the differential acetylation of 10 randomly selected differentially acetylated sites from the acetylome analysis. These results indicated that our acetylome analysis results were sufficiently reliable and reproducible. These results provide an essential resource for in-depth exploration of the stress responses and adaptation mechanisms associated with lysine acetylation in marine invertebrates and P. f. martensii.

Circulating exosome miRNA, is it the novel nutrient molecule through cross-kingdom regulation mediated by food chain transmission from microalgae to bivalve?

MicroRNAs (miRNAs) can efficiently regulate gene expression at intracellular and extracellular levels. Plant-derived miRNAs are highly enriched in animal haemolymph and regulate mammalian gene expression. However, evidence for food-derived miRNAs in Mollusca species is lacking. In this study, we fed the microalga Nannochloropsis oculata to the pearl oyster Pinctada fucata martensii and detected dietary miRNAs in exosomes isolated from the haemolymph by RNA-seq. In total, 273 endogenous miRNAs were identified in all biological replicates. We identified 23 microalgae-derived miRNAs in the exosomes of pearl oyster haemolymph. Most microalgae-derived miRNAs showed high expression levels in both exosomes and microalgae and exhibited apparent variation among individuals. These food-derived miRNAs were predicted to participate in endocytosis, apoptosis, signal transduction, energy metabolism, and biomineralization by targeting multiple genes. These findings demonstrated the cross-kingdom transport of miRNAs from microalgae to bivalves and provide insights into novel nutrient transmission through the food chain.

Cloning and functional analysis of a trypsin-like serine protease from Pinctada fucata martensii

Trypsin-like serine proteases (TLSs) play various roles in dietary protein digestion, hemolymph coagulation, antimicrobial peptide synthesis, and, in particular, the rapid immune pathways activated in response to pathogen detection. The cultured pearl industry, of which Pinctada fucata martensii is one of the most important species, is plagued by disease, thus leading to large economic losses. Herein, the molecular mechanisms underlying the innate immune response of P.f. martensii were explored. First, immune effector molecules from the P.f. martensii genome were screened and a TLS-like gene encoding a protein with a trypsin domain, herein designated as PmTLS, was identified. A multi-sequence alignment indicated a low sequence homology between PmTLS and other mollusk TLS-like proteins. Furthermore, a neighbor-joining phylogenetic analysis indicated that PmTLS has the closest genetic relationship to a Crassostrea gigas TLS. Additionally, real-time quantitative PCR (qPCR) analysis showed that PmTLS mRNA is constitutively expressed in all of the 6 examined P.f. martensii tissues, with significantly higher expression noted in hemocytes relative to the other tissues examined (p < 0.05). P.f. martensii samples were then challenged with various pathogen-associated molecular patterns (PAMPs), including lipopolysaccharide, peptidoglycan, and polyinosinic acid. In the challenge groups, PmTLS was significantly upregulated in hemocytes at 48 h post-challenge when compared to the unchallenged controls. Furthermore, treatment with recombinant PmTLS (rPmTLS) also significantly inhibited the growth of most of the examined gram-negative bacteria tested in vitro (p < 0.05), but it had little effect on the growth of the examined gram-positive bacteria. When examining morphological changes via transmission electron microscopy, rPmTLS treated bacteria exhibited morphological changes such as plasma wall separation. Thus, rPmTLS appears to play a bactericidal role by destroying bacterial cell membranes or cell walls, which subsequently leads to a release of the cellular contents and cell death. The findings presented herein have enabled further characterization of the immune defense mechanisms in P.f. martensii and may lead to improved disease control methods for the pearl cultivation industry.

High gene flow in the silverlip pearl oyster Pinctada maxima between inshore and offshore sites near Eighty Mile Beach in Western Australia

An understanding of stock recruitment dynamics in fisheries is fundamental to successful management. Pinctada maxima is a bivalve mollusc widely distributed throughout the Indo-Pacific and is the main species targeted for cultured pearl and pearl shell production in Australia. Pearl production in Australia relies heavily on wild-caught individuals, the majority of which come from the Eighty Mile Beach region near Broome in Western Australia. In this study, we used a genotyping by sequencing approach to explore fine-scale patterns of genetic connectivity among inshore shallow and offshore deep populations of P. maxima near Eighty Mile Beach. Our results revealed high-levels of gene flow among inshore and offshore sites and no differences in genetic diversity between depths. Global estimates of genetic differentiation were low (F ST = 0.006) but significantly different from zero, and pairwise estimates of genetic differentiation among sites were significant in only 3% of comparisons. Moreover, Bayesian clustering detected no separation of inshore and offshore sample sites, and instead showed all samples to be admixed among sites, locations and depths. Despite an absence of any clear spatial clustering among sites, we identified a significant pattern of isolation by distance. In a dynamic environment like Eighty Mile Beach, genetic structure can change from year-to-year and successive dispersal and recruitment events over generations likely act to homogenize the population. Although we cannot rule out the null hypothesis of panmixia, our data indicate high levels of dispersal and connectivity among inshore and offshore fishing grounds.

Characterization of transcription factor activator pretein-1 (AP-1) and its association with cold tolerance in Pinctada fucata martensii

AP-1 is an important transcription factor for cell proliferation/differentiation and animal immunity/development; however, its role in research in shellfish is poorly understood. Here, the cDNA of AP-1 gene from Pinctada fucata martensii was characterized. Its expression was detected in all six examined tissues, and a high level was observed in the gill and hepatopancreas. Analysis of the developmental transcriptomes showed that the PmAP-1 gene expression levels were high during D-stage larval and spat stages. The gene also exhibited a significantly high expression under cold tolerance stress. SNP analysis of the exon region and 5' flanking region of PmAP-1 revealed 19 SNPs of which 8 showed significant differences between cold tolerance selection line and base stock. Furthermore, three haplotypes generated by the SNPs of PmAP-1 were significantly associated with cold tolerance, respectively.These results suggest that the PmAP-1 gene plays an important role in the response of P. f. martensii to low temperature stress. These SNPs and haplotypes of PmAP-1 may be related to the cold tolerance of P. f. martensii, and could be candidate markers potentially for further selective breeding.

Exploring crucial molecular events in pearl oyster after pre-grafting conditioning by genome-wide bisulfite sequencing for DNA methylation analysis

Pre-grafting condition is an important method to promote recovery from transplant surgery during pearl production. In the present study, we constructed two DNA methylomes from pearl oysters with and without conditioning to investigate the molecular mechanism of the pearl oyster Pinctada fucata martensii underlying the pre-grafting condition. A total of 4,594,997 and 4,930,813 methyl CG in the control (Con) and pre-grafting group (PT) were detected, resulting in the whole genome methylation profile and methylation pattern in P. f. martensii. Results reveal that the promoter, especially the CpG island-rich region, was more infrequently methylated than the gene function elements in P. f. martensii. A total of 51,957 differently methylated regions (DMRs) between Con and PT were obtained, including 3789 DMR in the promoter and 16,021 in the gene body. Based on gene ontology and pathway enrichment analyses, these DMRs were mainly related to "cellular process", "metabolic process", "Epstein-Barr virus infection", and "Fanconi anemia pathway". The methylation site in the promoter region may be associated with the promoter activity and transcription factor binding. These results help our understanding of the mechanism of pre-grafting condition, thereby providing key information in guiding to improve the conditioning methods for enhanced pearl oyster survival rate after transplantation.

Ecdysone signal pathway participates in shell formation in pearl oysters Pinctada fucata martensii

Ecdysone exists in arthropods, Mollusca and other invertebrates and plays vital roles in exoskeleton formation of Ecdysozoa. However, little is known about its functions in bivalve species. Herein, we identified ecdysone from the serum of pearl oyster Pinctada fucata martensii and obtained the coding sequence of ecdysone receptor (PmEcR) and homologue of its heterodimer protein retinoid X receptor (PmRXR). The deduced amino acid sequences of PmEcR and PmRXR contained a DNA-binding and ligand-binding domain and were very similar to the orthologs of other species. Moreover, PmEcR and PmRXR were located in the nuclei and cytoplasm of HEK-293T cells. PmEcR and PmRXR were highly expressed in early embryos and biomineralized mantle tissue. Moreover, the serum concentration of ecdysone significantly increased at 2, 4, 6, and 8 h post-shell notching. The expression of PmEcR in the mantle tissue was significantly induced at the corresponding time points, while that of PmRXR was significantly induced at 6 h. Ecdysone stimulation remarkably induced the expression of growth factors (BMP2 and BMP7), transcription factors (PmRunt and AP-1), and shell matrix protein genes (chitinase, lysine-rich matrix protein (KRMP), TYR2, and PmCOLVI), which indicated that ecdysone signaling plays important roles in shell repair. However, yeast two-hybrid assay and bimolecular fluorescence complementation showed that PmEcR and PmRXR did not form dimers, suggesting the different molecular interactions of EcR in bivalves. These findings provide insights into the function of ecdysone and its regulation pathway in bivalve species.

Comparative proteomics and transcriptomics illustrate the allograft-induced stress response in the pearl oyster (Pinctada fucata martensii)

Implantation of a spherical nucleus into a recipient oyster is a critical step in artificial pearl production. However, the molecular mechanisms underlying the response of the pearl oyster to this operation are poorly understood. In this research, we used transcriptomic and proteomic analyses to examine allograft-induced changes in gene/protein expression patterns in Pinctada fucata martensii 12 h after nucleus implantation. Transcriptome analysis identified 688 differential expression genes (DEGs) (FDR<0.01 and |fold change) ＞ 2). Using a 1.2-fold increase or decrease in protein expression as a benchmark for differentially expressed proteins (DEPs), 108 DEPs were reliably quantified, including 71 up-regulated proteins (DUPs) and 37 down-regulated proteins (DDPs). Further analysis revealed that the GO terms, including "cellular process", "biological regulation" and "metabolic process" were considerably enriched. In addition, the transcriptomics analysis showed that "Neuroactive ligand-receptor interaction", "NF-kappa B signaling pathway", "MAPK signaling pathway", "PI3K-Akt signaling pathway', "Toll-like receptor signaling pathway", and "Notch signaling pathway" were significantly enriched in DEGs. The proteomics analysis showed that "ECM-receptor interaction", "Human papillomavirus infection", and "PI3K-Akt signaling pathway" were significantly enriched in DEPs. The results indicate that these functions could play an important role in response to pear oyster stress at nucleus implantation. To assess the potential relevance of quantitative information between mRNA and proteins, using Ward's hierarchical clustering analysis clustered the protein/gene expression patterns across the experimental and control samples into six groups. To investigate the biological processes associated with the protein in each cluster, we identified the significantly enriched GO terms and KEGG pathways in the proteins in each cluster. Gene set enrichment analysis (GSEA) was used to reveal the potential protein or transcription pathways associated with the response to nuclear implantation. Thus, the study of P. f. martensii is essential to enhance our understanding of the molecular mechanisms involved in pearl biosynthesis and the biology of bivalve molluscs.

The byssal-producing glands and proteins of the silverlip pearl oyster Pinctada maxima (Jameson, 1901)

Pinctada maxima are most well known for their production of high-quality natural pearls. They also generate another natural material, the byssus, an adhesive thread critical for steadfast attachment underwater. Herein, P. maxima byssal threads were analysed via proteotranscriptomics to reveal 49 proteins. Further characterisation was undertaken on five highly expressed genes: glycine-rich thread protein (GRT; also known as PUF3), apfp1/perlucin-like protein (Pmfp1); peroxidase; thrombospondin 1, and Balbiani ring 3 (BR3), which showed localised tissue expression. The spatial distribution of GRT and Pmfp1 via immunodetection combined with histology helped to identify glandular regions of the foot that contribute to byssal thread production: the byssal gland, the duct gland, and two thread-forming glands of basophilic and acidophilic serous-like cells. This work advanced primary knowledge on the glands involved in the creation of byssal threads and the protein composition of the byssus for P. maxima, providing a platform for the design of marine biopolymers.

Microplastics impact shell and pearl biomineralization of the pearl oyster Pinctada fucata

Microplastics are extremely widespread aquatic pollutants that severely detriment marine life. In this study, the influence of microplastics on biomineralization was investigated. For the first time, multiple forms and types of microplastics were detected and isolated from the shells and pearls of Pinctada fucata. According to the present study, the abundance of microplastics in shells and pearls was estimated at 1.95 ± 1.43 items/g and 0.53 ± 0.37 items/g respectively. Interestingly, microplastics were less abundant in high-quality round pearls. Microplastics may hinder the growth of calcite and aragonite crystals, which are crucial components required for shell formation. During the process of biomineralization microplastics became embedded in shells, suggesting the existence of a novel pathway by which microplastics accumulate in bivalves. After a 96-h exposure to microplastics, the expression level of typical biomineralization-related genes increased, including amorphous calcium carbonate binding protein (ACCBP) gene which experienced a significant increase. ACCBP promotes the formation of amorphous calcium carbonate (ACC), which is the pivotal precursor of shell formation-related biominerals. ACCBP is highly expressed during the developmental stage of juvenile oysters and the shell-damage repair process. The increased expression of ACCBP suggests biomineralization is enhanced as a result of microplastics exposure. These results provide important evidence that microplastics exposure may impact the appearance of biominerals and the expression of biomineralization-related genes, posing a new potential threat to aquatic organisms.

Understanding connectivity of pearl oyster populations within Tuamotu atoll semi-closed lagoons: Cumulative insight from genetics and biophysical modelling approaches

Connectivity affects species demography, (meta)population dynamics, evolution, phylogeny and biogeography. Various methodological approaches are applied to measure connectivity. Biophysical modelling can explore systematically the influence of atmospheric, oceanic and ecological forcing, while genetics measures connectivity patterns within the sampling strategy limit. In the Pacific Ocean pearl farming lagoons, the activity relies on spat collecting of the black lipped pearl oyster Pinctada margaritifera occurring after the larval dispersal phase, which follows spawning from wild or farmed populations. Biophysical 3D modelling and genomic studies have both separately brought insights on within-lagoon connectivity and on the origin of spats. Here, we combined previous genetics results with new realistic biophysical modelling scenarios to elucidate connectivity in Ahe Atoll lagoon. When combined, we identified the weather sequence likely explaining the realized connectivity observations. We discuss the strengths, weaknesses, opportunities and threats of combining these two approaches considering specific pearl farming demographic connectivity questions.

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.

PU14, a Novel Matrix Protein, Participates in Pearl Oyster, Pinctada Fucata, Shell Formation

Biomineralization is a widespread biological process, involved in the formation of shells, teeth, and bones. Shell matrix proteins have been widely studied for their importance during shell formation. In 2015, our group identified 72 unique shell matrix proteins in Pinctada fucata, among which PU14 is a matrix protein detected in the soluble fraction that solely exists in the prismatic layer. However, the function of PU14 is still unclear. In this study, the full-length cDNA sequence of PU14 was obtained and functional analyses of PU14 protein during shell formation were performed. The deduced protein has a molecular mass of 77.8 kDa and an isoelectric point of 11.34. The primary protein structure contains Gln-rich and random repeat units, which are typical characteristics of matrix protein and indicate its potential function during shell formation. In vivo and in vitro experiments indicated PU14 has prismatic layer functions during shell formation. The tissue expression patterns showed that PU14 was mainly expressed in the mantle tissue, which is consistent with prismatic layer formation. Notching experiments suggested that PU14 responded to repair and regenerate the injured shell. After inhibiting gene expression by injecting PU14-specific double-stranded RNA, the inner surface of the prismatic layer changed significantly and became rougher. Further, in vitro experiments showed that recombinant protein rPU14 impacted calcite crystal morphology. Taken together, characterization and functional analyses of a novel matrix protein, PU14, provide new insights about basic matrix proteins and their functions during shell formation.

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.

Novel globular C1q domain-containing protein (PmC1qDC-1) participates in shell formation and responses to pathogen-associated molecular patterns stimulation in Pinctada fucata martensii

The C1q protein, which contains the globular C1q (gC1q) domain, is involved in the innate immune response, and is found abundantly in the shell, and it participates in the shell formation. In this study, a novel gC1q domain-containing gene was identified from Pinctada fucata martensii (P. f. martensii) and designated as PmC1qDC-1. The full-length sequence of PmC1qDC-1 was 902 bp with a 534 bp open reading frame (ORF), encoding a polypeptide of 177 amino acids. Quantitative real-time PCR (qRT-PCR) result showed that PmC1qDC-1 was widely expressed in all tested tissues, including shell formation-associated tissue and immune-related tissue. PmC1qDC-1 expression was significantly high in the blastula and gastrula and especially among the juvenile stage, which is the most important stage of dissoconch shell formation. PmC1qDC-1 expression was located in the outer epithelial cells of mantle pallial and mantle edge and irregular crystal tablets were observed in the nacre upon knockdown of PmC1qDC-1 expression at mantle pallial. Moreover, the recombined protein PmC1qDC-1 increased the rate of calcium carbonate precipitation. Besides, PmC1qDC-1 expression was significantly up-regulated in the mantle pallial at 6 h and was significantly up-regulated in the mantle edge at 12 h and 24 h after shell notching. The expression level of PmC1qDC-1 in mantle edge was significantly up-regulated at 48 h after LPS stimulation and was significantly up-regulated at 12 h, 24 h and 48 h after poly I:C stimulation. Moreover, PmC1qDC-1 expression was significantly up-regulated in hemocytes at 6 h after lipopolysaccharide (LPS) and poly I:C challenge. These findings suggest that PmC1qDC-1 plays a crucial role both in the shell formation and the innate immune response in pearl oysters, providing new clues for understanding the shell formation and defense mechanism in mollusk.

Phylogenetic comparisons reveal mosaic histories of larval and adult shell matrix protein deployment in pteriomorph bivalves

Molluscan shells are organo-mineral composites, in which the dominant calcium carbonate is intimately associated with an organic matrix comprised mainly of proteins and polysaccharides. However, whether the various shell matrix proteins (SMPs) date to the origin of hard skeletons in the Cambrian, or whether they represent later deployment through adaptive evolution, is still debated. In order to address this issue and to better understand the origins and evolution of biomineralization, phylogenetic analyses have been performed on the three SMP families, Von Willebrand factor type A (VWA) and chitin-binding domain-containing protein (VWA-CB dcp), chitobiase, and carbonic anhydrase (CA), which exist in both larval and adult shell proteomes in the bivalves, Crassostrea gigas and Pinctada fucata. In VWA-CB dcp and chitobiase, paralogs for larval and adult SMPs evolved before the divergence of these species. CA-SMPs have been taken as evidence for ancient origins of SMPs by their presumed indispensable function in biomineralization and ubiquitous distribution in molluscs. However, our results indicate gene duplications that gave rise to separate deployments as larval and adult CA-SMPs occurred independently in each lineage after their divergence, which is considerably more recent than hitherto assumed, supporting the "recent heritage and fast evolution" scenario for SMP evolution.

Molecular and functional analysis of PmC1qDC in nacre formation of Pinctada fucata martensii

The C1q-domain-containing (C1qDC) proteins are a family of proteins characterized by a globular C1q (gC1q) domain in their C-terminus which hold the potential function in the shell formation as shell matrix proteins. In this study, a C1qDC protein was identified and characterized in pearl oyster (Pinctada fucata martensii) (PmC1qDC) to explore its function in nacre formation. The PmC1qDC-deduced protein sequence carried a typical globular C1q (gC1q) domain that possessed the typical 10-stranded β-sandwich fold with a jelly-roll topology common to all C1qDC family members and shared high homology with other gC1q domains. Homologous analysis of PmC1qDC presented it contained conserved secondary structure and Phe135, Phe155, Tyr166, Phe173, Tyr181, Phe183, and Phe256 amino acid residues. Expression pattern analysis showed that PmC1qDC expressed in all the detected tissues and exhibited a significantly higher expression level in nacre formation-associated tissues. After the shell notching, the expression level of PmC1qDC showed significantly up-regulation after 12 h in the central zone of mantle (MC). PmC1qDC expression significantly decreased in the MC after RNA interference (RNAi). Furthermore, disordered crystals with evident rough surface and irregular crystal tablets were observed in the nacre after RNAi. Results suggested that PmC1qDC affects the shell nacre formation, which is significant to improve the pearl production of pearl oyster.

Antimicrobial properties and immune-related gene expression of a C-type lectin isolated from Pinctada fucata martensii

C-type lectins are carbohydrate-binding proteins that play important roles in the innate immune response to pathogen infections. Here, multi-step high performance liquid chromatography (HPLC), combined with mass spectrometry (MS), was used to isolate and identify proteins with antibacterial activity from the serum of Pinctada fucata martensii. Using this method, we obtained a novel isoform of C-type lectin (PmCTL-1). PmCTL-1 strongly inhibited gram-positive bacteria. The complete cDNA sequence of PmCTL-1 was 636 bp in length, and encoded a protein 149 amino acids long, containing a typical carbohydrate recognition domain (CRD). A phylogenetic analysis based on a multiple sequence alignment indicated that PmCTL-1 was highly similar to C-type lectins from other mollusks. Fluorescent quantitative real-time PCR analysis showed that PmCTL-1 mRNA was strongly upregulated in the mantle of healthy P.f. martensii, but was expressed only at low levels in the gill, gonad, hepatopancreas, adductor muscle, and hemocytes. PmCTL-1 expression levels in the mantle and hemocytes increased significantly in response to bacterial stimulation. This study provides a valuable framework for further explorations of innate immunity and the immune response in mollusks.

Identification of a long non-coding RNA (LncMSEN2) from pearl oyster and its potential roles in exoskeleton formation and LPS stimulation

Long non-coding RNAs (lncRNAs) play regulatory roles in various biological processes, including exoskeleton formation and immune response. The exoskeleton-based mantle-shell defense system is an important defense mechanism in shellfish. In this study, we found a novel lncRNA, herein formally named, LncMSEN2, from the pearl oyster Pinctada fucuta martensii, and its sequence was validated via polymerase chain reaction (PCR). LncMSEN2 was highly expressed in mantle tissues, especially in the central region (P < 0.05), and was also expressed in the pearl sac as detected by quantitative real-time PCR. In situ hybridization experiments revealed that LncMSEN2 had a strong positive signal in the inner and outer epidermal cells of the mantle pallial and central regions. RNA interference experiments showed that interference of LncMSEN2 expression with dsRNA in mantle tissues led to an abnormal crystal structure of the nacre. In addition, LncMSEN2 expression significantly increased 6 h after lipopolysaccharide stimulation in mantle tissues (P < 0.05). These results indicated that LncMSEN2 may be a novel regulator of the mantle-shell defense system of pearl oyster.

A species-specific miRNA participates in biomineralization by targeting CDS regions of Prisilkin-39 and ACCBP in Pinctada fucata

Biomineralization is a sophisticated biological process precisely regulated by multiple molecules and pathways. Accumulating miRNAs have been identified in invertebrates but their functions in biomineralization are poorly studied. Here, an oyster species-specific miRNA, novel_miR_1 was found to regulate biomineralization in Pinctada fucata. Target prediction showed that novel_miR_1 could target Prisilkin-39 and ACCBP by binding to their coding sequences (CDS). Tissue distribution analysis revealed that the expression level of novel_miR_1 was highest in the mantle, which was a key tissue participating in biomineralization. Gain-of-function assay in vivo showed that biomineralization-related genes including Prisilkin-39 and ACCBP were down-regulated and shell inner surfaces of both prismatic and nacreous layer were disrupted after the over-expression of novel_miR_1, indicating its dual roles in biomineralization. Furthermore, the shell notching results indicated that novel_miR_1 was involved in shell regeneration. Dual-luciferase reporter assay in vitro demonstrated that novel_miR_1 directly suppressed Prisilkin-39 and ACCBP genes by binding to the CDS regions. Taken together, these results suggest that novel_miR_1 is a direct negative regulator to Prisilkin-39 and ACCBP and plays an indispensable and important role in biomineralization in both prismatic and nacreous layer of P. fucata.

Genome-wide DNA Methylation Analysis of Mantle Edge and Mantle Central from Pearl Oyster Pinctada fucata martensii

DNA methylation is a type of epigenetic modification that alters gene expression without changing the DNA sequence and mediates some cases of phenotypic plasticity. In this study, we identified six DNA methyltransferase (DNMT) genes and two methyl-CpG binding domain protein2 (MBD2) gene from Pinctada fucata martensii. We also analyzed the genome-wide DNA methylation levels of mantle edge (ME) and mantle central (MC) from P. f. martensii via methylated immunoprecipitation sequencing (MeDIP-Seq). Results revealed that both ME and MC had 122 million reads, and had 58,702 and 55,721 peaks, respectively. The obtained methylation patterns of gene elements and repeats showed that the methylation of the protein-coding genes, particularly intron and coding exons (CDSs), was more frequent than that of other genomic elements in the pearl oyster genome. We combined the methylation data with the RNA-seq data of the ME and MC of P. f. martensii and found that promoter, CDS, and intron methylation levels were positively correlated with gene expression levels except the highest gene expression level. We also identified 313 differential methylation genes (DMGs) and annotated 212 of them. These DMGs were significantly enriched in 30 pathways, such as amino acid and protein metabolism, energy metabolism, terpenoid synthesis, and immune-related pathways. This study comprehensively analyzed the methylomes of biomineralization-related tissues and helped enhance our understanding of the regulatory mechanism underlying shell formation.

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

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

Molecular characterization and functional analysis of a Rel gene in the Pacific oyster

The nuclear factor-κB (NF-κB) signaling pathway plays a crucial role in regulating many physiological processes such as development, inflammation, apoptosis, cell proliferation, differentiation and immune responses. And the NF-κB/Rel family members were considered as the most important transcription factors in the NF-κB signaling pathway. In this study, we cloned a Rel homolog gene (named as CgRel2) from the Pacific oyster, Crassostrea gigas. The 2115-bp open reading frame (ORF) encodes 704 amino acids and CgRel2 possesses a conserved Rel Homology Domain (RHD) at the N-terminus. Phylogenetic analysis revealed that CgRel2 is most closely related to Pinctada fucata dorsal protein. CgRel2 transcripts are widely expressed in all tested tissues, with the highest expression observed in the labial palp and the gill. Moreover, the expression of CgRel2 is significantly upregulated after lipopolysaccharide (LPS), peptidoglycan (PGN), and polyinosinic-polycytidylic acid [poly(I:C)] challenge. CgRel2 transfection into human cell lines activated NF-κB, TNFα and oyster IL-17 (CgIL-17) reporter genes in a dose-dependent manner, while CgRel2 overexpression cannot induce ISRE (Interferon stimulation response element) reporter gene's transcriptional activity. Additionally, the results of co-immunoprecipitation showed that CgRel2 or CgRel1 could interact with oyster IκB1, IκB2 and IκB3 proteins strongly, which may be critical for the immune signaling transduction and the regulation of its immune functions. Together, these results suggest that CgRel2 could respond to pathogenic infection, participate in the immune signal transduction and activate NF-κB, TNFα and CgIL-17 reporter genes. Thus, CgRel2 could play an important role in the oyster immune system.

Evolutionary and functional analysis of MyD88 genes in pearl oyster Pinctada fucata martensii

Myeloid differentiation factor 88 (MyD88) is an adapter protein that links toll-like receptor and interleukin 1 receptor-mediated signal transduction. In this study, we identified 20 MyD88 genes from eight mollusk genomes and found that MyD88 was expanded in bivalves. This expansion tends to be tandem duplication. Phylogenetic analysis suggested that the tandem duplication of MyD88 was formed before bivalve differentiation. All of the identified MyD88 contained both of death domain (DD) and toll/interleukin-1 receptor (TIR) domain, and 13 mollusks MyD88 have low complexity regions (LCRs), which were not found in the MyD88 from humans and zebrafish. The genomic structure showed that most of the mollusk MyD88 (14 of 19) contained five conserved introns, four of which were found in humans and zebrafish. Furthermore, the cDNA full length of PfmMyD88-2 (one of the two identified MyD88 in Pincatada fucata martensii) was obtained with 1591 bp, including 260 bp of 5'UTR, 257 bp of 3'UTR, and 1077 bp of open reading frame encoding 358 amino acids. Quantitative real-time PCR analysis demonstrated that PfmMyD88-2 mRNA was widely expressed in all detected tissues. The highest expression level was in the gills and followed by hepatopancreas and feet. After lipopolysaccharide stimulation, PfmMyD88-2 expression level increased and reached the highest level at 12 h and then gradually declined to the normal level. Over-expression of PfmMyD88-2 in HEK293T increased the luciferase activity of the pNF-κB-Luc reporter. We also identified that PfmmiR-4047 could regulate the expression of PfmMyD88-2. These results help us elucidate the mechanism underlying mollusk immune response.

PfSMAD1/5 Can Interact with PfSMAD4 to Inhibit PfMSX to Regulate Shell Biomineralization in Pinctada fucata martensii

The BMP2 signal transduced by SMAD1/5 plays an important role in osteoblast differentiation and bone formation. Shell formation of Pinctada fucata martensii is a typical biomineralization process that is similar to that of teeth/bone formation. However, whether the Pinctada fucata BMP2 (PfBMP2) signal transduced by PfSMAD1/5 occurs in P. f. martensii, how the PfBMP2 signal is transduced by PfSMAD1/5, and how PfSMAD1/5 regulates the biomineralization process in this species and other shellfish are poorly understood. Therefore, injection experiments of recombinant PfBMP2 and inhibitor dorsomorphin revealed that PfSMAD1/5 can transduce PfBMP2 signals. Subcellular localization and bimolecular fluorescence complementation assays indicated that PfSMAD1/5 phosphorylated by PfBMPR1b interacts with PfSMAD4 in the cytoplasm to form a complex, which translocates to the nucleus to transduce PfBMP2 signals. Co-immunoprecipitation and luciferase assays revealed that PfSMAD1/5 may interact with PfMSX to dislodge it from its binding element, resulting in initiation of mantle gene transcription. The in vivo functional assay showed that knockdown of PfMSAD1/5 decreased expression of shell matrix genes and disordered the nacreous layer, and the correlation assay of shell regeneration showed the concomitant expression pattern of PfSMAD1/5 and shell matrix genes. Together, these data showed that PfSMAD1/5 can transduce PfBMP2 signals to regulate shell biomineralization in P. f. martensii, which illustrated conservation of the BMP2-SMAD signal pathway among invertebrates. Particularly, the results suggest that there is only one PfMSX gene, which functions like the Hox gene in vertebrates, that interacts with PfSMAD1/5 in a protein-protein action form and plays the role of transcription repressor.

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.

Survival, retention rate and immunity of the black shell colored stocks of pearl oyster Pinctada fucata martensii after grafting operation

We have developed a black shell colored selected line observed to have higher survival ability. In this study, to understand its immune capacity, total carotenoid content (TCC) of the black shell colored line (BG) and the control group (CG) were compared. Survival and retention rates, immunity and antioxidant capacity of BG were compared relative to CG at different times after grafting operation. The results showed that BG had significantly larger TCC than CG (P < 0.05). BG had significantly higher survival and retention rates than CG on days 7, 30 and 360 after grafting (P < 0.05). On days 360, BG had significantly larger pearl thickness than CG (P < 0.05). BG exhibited increased ACP, AKP, SOD, CAT, TAOC and LZ activity than the CG on 0 h, 12 h, 1 d, 3 d, 5 d, 7 d and 30 d after grafting. BG had higher expression levels of Fascin, SOD, CDK-7, CDAP-1, IRAK-1, α2m, GST-1, TRAF-3 and Caspase-2 than CG. The results suggested that BG had higher immune competence and pearl production performances, which is promising to improve pearl quality and production.

The role of a new insulin-like peptide in the pearl oyster Pinctada fucata martensii

Pinctada fucata martensii, is an economically important marine bivalve species cultured for seawater pearls. At present, we know little about the molecular mechanisms of the insulin signalling pathway in this oyster. Herein, we cloned and analysed an insulin-like peptide (PfILP) and its signalling pathway-related genes. We detected their expression levels in different tissues and developmental stages. Recombinant PfILP protein was produced and found to significantly increase primary mantle cell activity and induce the expression of the proliferating cell nuclear antigen (PCNA) gene. PfILP could also regulate the 293T cell cycle by stimulating the S phase and inhibiting the G1 and G2 phases. Recombinant PfILP protein induced the expression of its signalling pathway-related genes in mantle cells. In vitro co-immunoprecipitation analysis showed that PfILP interacts with PfIRR. PfILP activated expression of the pfIRR protein, and also activated the mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase (PI3K) pathways by stimulating phosphorylation of MAPK and AKT. Further analysis showed that PfILP up-regulated glycogen synthesis-related genes glycogen synthase kinase-3 beta (GSK-3β), protein phosphatase 1 (PP1) and glucokinase (GK) at the mRNA level, as well as the expression of the PP1 protein, and phosphorylation of GSK-3β. These results confirmed the presence of a conserved insulin-like signalling pathway in pearl oyster that is involved in cell activity, glycogen metabolism, and other physiological processes.

LncMSEN1, a mantle-specific LncRNA participating in nacre formation and response to polyI:C stimulation in pearl oyster Pinctada fucata martensii

Long noncoding RNA (LncRNA) regulates various life processes, including biomineralization and innate immune response through complex mechanisms. In this research, we identified a LncRNA named LncMSEN1 from pearl oyster Pinctada fucata martensii. LncMSEN1 sequence was validated by PCR, and its expression was high in mantle tissues according to qRT-PCR. LncMSEN1 was co-located with the nacre matrix protein N-U8 and fibrinogen domain-containing protein. And LncMSEN1 and N-U8 expression levels in the mantle were positively correlated. RNA interference was used to detect its effect on nacre formation in shells. Results showed that the decreased LncMSEN1 expression in mantle can cause the disordered growth of crystals on the inner surface of nacre in the shells, as well as the decrease expression of N-U8. In addition, the LncMSEN1 expression level significantly increased at 24 h after polyI:C stimulation in the mantle (P < 0.05). These findings suggested the involvement of LncMSEN1 in the formation of nacre in shells and related to innate immune response in pearl oyster, which provided additional insights into the roles of LncRNAs in pearl oysters.

Separation, identification and gene expression analysis of PmAMP-1 from Pinctada fucata martensii

Antibacterial peptides (AMPs) constitute an important part of the body's innate immune system and are responsible for a wide range of inhibitory effects against pathogens such as bacteria, fungi, and viruses. In this study, multi-step high performance liquid chromatography (HPLC), combined with Mass Spectrometry (MS), was used to isolate and identify proteins with antibacterial activity from the serum of Pinctada fucata martensii (P.f. Martensii) and obtain a component named P.f. Martensii antimicrobial peptide-1 (PmAMP-1). PmAMP-1 cDNA was cloned and sequenced by rapid amplification of cDNA ends (RACE) and mRNA expression of was analyzed by quantitative real-time PCR (qRT-PCR). From the results of this study, full-length PmAMP-1 cDNA was shown to be 700 base pairs (bp) long with an open reading frame (ORF) of 294 bp, encoding 97 amino acids with a predicted structure that is mostly α-helices. PmAMP-1 mRNA was constitutively expressed in all tested tissues including the adductor muscle, mantle, hepatopancreas, gill, gonads and hemocytes. The highest level of PmAMP-1 transcription was observed at 8 h and 2 h after bacterial challenge in hemocytes and adductor muscle (p < 0.01), respectively. Furthermore, PmAMP-1 caused significant morphological alterations in E. coli, as shown by transmission electron microscopy (TEM). The results from this study provide a valuable base for further exploration of molluscan innate immunity and immune response.

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.

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.