MicroRNAs (miRNAs) are a class of non-coding RNA molecules with presumed post-transcriptional regulatory activity in various biological processes, such as development and biomineralization. Pinctada martensii is one of the main species cultured for marine pearl production in China and Japan. In our previous research, 258 pm-miRNAs had been identified by solexa deep sequencing in P. martensii, while it is far from the number of miRNAs found in other species. In this study, based on the transcriptome database of pearl sac, we identified 30 candidate pm-miRNAs by computational prediction. Among the obtained 30 pm-miRNAs, 13 pm-miRNAs were generated from the complementary strand of protein-coding mRNAs, and 17 pm-miRNAs could not be annotated using blastx and tblastn analysis. Notably, 10 of the 30 pm-miRNAs, such as pm-miR-1b, pm-miR-205b and pm-miR-375b, were homologous with the reported pm-miRNAs, respectively. To validate the existence of the identified pm-miRNAs, eight randomly selected pm-miRNAs were tested by stem loop quantitative RT-PCR analyses using 5.8S as the internal reference gene. Target prediction between the obtained pm-miRNAs and biomineralization-related genes by microTar, miRanda and RNA22 indicated pm-miR-2386 and pm-miR-13b may be the key factors in the regulation network by regulating the formation of organic matrix or the differentiation of mineralogenic cell during shell formation. Thus, this study enriched miRNA databases of pearl oyster and provided a new way to understand biomineralization.

MiR-29a is a conserved miRNA that participates in bone formation and immune response in vertebrates. miR-29a of Pinctada martensii (Pm-miR-29a) was identified in the previous research though deep sequencing. In this report, the precise sequence of mature Pm-miR-29a was validated using miRNA rapid amplification of cDNA ends (miR-RACE) technology. The precursor sequence of Pm-miR-29a was predicted to have 87 bp. Stem loop qRT-PCR analysis showed that Pm-miR-29a was easily detected in all the tissues, although expressions in the mantle and gill were low. The microstructure showed the disrupted growth of the nacre after Pm-miR-29a over-expression, which was induced by mimic injection into P. martensii. Results of the target analysis indicated that neuropeptide Y receptor type 2 (Y2R) was the potential target of Pm-miR-29a. Meanwhile, Pm-miR-29a mimics could obviously inhibit the relative luciferase activity of the reporter containing 3′ UTR (Untranslated Regions) of the Y2R gene. Furthermore, the expression of Y2R was downregulated whereas expressions of interleukin 17 (IL-17) and nuclear factor κB (NF-κB) were upregulated after Pm-miR-29a over-expression in the mantle and gill, thereby suggesting that Pm-miR-29a could activate the immune response of the pearl oyster. Results showed that Pm-miR-29a was involved in nacre formation and immune response by regulating Y2R in pearl oyster P. martensii.

Bone morphogenetic protein 7 (BMP7), also called osteogenetic protein-1, can induce bone formation. In this study, the obtained full-length cDNA of BMP7 from Pinctada martensii (Pm-BMP7) was 2972 bp, including a 5'-untranslated region (UTR) of 294 bp, an open reading fragment of 1290 bp encoding a 429 amino acid polypeptide and a 3'-UTR of 1388 bp. The deduced protein sequence of Pm-BMP7 contained a signal peptide, a pro-domain and a mature peptide. The mature peptide consisted of 135 amino acids and included a transforming growth factor β family domain with six shared cysteine residues. The protein sequence of Pm-BMP7 showed 66% identity with that from Crassostrea gigas. Two unigenes encoding Pm-BMPRI (Pm-BMP receptor I) and Pm-BMPRII were obtained from the transcriptome database of P. martensii. Tissue expression analysis demonstrated Pm-BMP7 and Pm-BMPRI were highly expressed in the mantle (shell formation related-tissue), while Pm-BMPRII was highly expressed in the foot. After inhibiting Pm-BMP7 expression using RNA interference (RNAi) technology, Pm-BMP7 mRNA was significantly down-regulated (p < 0.05) in the mantle pallium (nacre formation related-tissue) and the mantle edge (prismatic layer formation related-tissue). The microstructure, observed using a scanning electron microscope, indicated a disordered growth status in the nacre and obvious holes in the prismatic layer in the dsRNA-Pm-BMP7 injected-group. These results suggest that Pm-BMP7 plays a crucial role in the nacre and prismatic layer formation process of the shell.

MicroRNAs (miRNAs) are noncoding RNA molecules that function as negative regulators of target genes. In our previous research, 258 pm-miRNAs were identified in Pinctada martensii by Solexa deep sequencing. Pm-miR-2305 was one of the identified pm-miRNAs with a potential function in biomineralization. In the present study, the precursor of pm-miR-2305 was predicted with 96 bp, containing a characteristic hairpin structure. Stem-loop qRT-PCR analysis indicated that pm-miR-2305 was constitutively expressed in all the tissues (adductor muscle, gill, mantle, hepatopancreas, foot, and gonad) of P. martensii and was highly expressed in the foot. After the over-expression of pm-miR-2305 in the mantle by mimics injection into the muscle of P. martensii, nacre demonstrated disorderly growth, as detected by scanning electron microscopy. Dual luciferase reporter gene assay indicated that pm-miR-2305 mimics could significantly inhibit the luciferase activity of the reporter containing the 3′UTR of the pearlin gene. Western blot analysis demonstrated that the protein expression of pearlin was down-regulated in the mantle tissue after the over-expression of pm-miR-2305. Therefore, our data showed that pm-miR-2305 participated in nacre formation by targeting pearlin in P. martensii.

The paired-box 3 (Pax3) is a transcription factor and it plays an important part in melanin synthesis. In this study, a new Pax3 gene was identified from Pteria penguin (Röding, 1798) (P. penguin) by RACE-PCR (rapid-amplification of cDNA ends-polymerase chain reaction) and its effect on melanin synthesis was deliberated by RNA interference (RNAi). The cDNA of PpPax3 was 2250 bp long, containing an open reading fragment of 1365 bp encoding 455 amino acids. Amino acid alignment and phylogenetic tree showed PpPax3 shared the highest (69.2%) identity with Pax3 of Mizuhopecten yessoensis. Tissue expression profile showed that PpPax3 had the highest expression in mantle, a nacre-formation related tissue. The PpPax3 silencing significantly inhibited the expression of PpPax3, PpMitf, PpTyr and PpCdk2, genes involved in Tyr-mediated melanin synthesis, but had no effect on PpCreb2 and an increase effect on PpBcl2. Furthermore, the PpPax3 knockdown obviously decreased the tyrosinase activity, the total content of eumelanin and the proportion of PDCA (pyrrole-2,3-dicarboxylic acid) in eumelanin, consistent with influence of tyrosinase (Tyr) knockdown. These data indicated that PpPax3 played an important regulating role in melanin synthesis by Tyr pathway in P. penguin.

Marine pearl production is directly influenced by the growth speed of Pinctada fucata martensii. However, the slow growth rate of this organism remains the main challenge in aquaculture production. Epidermal growth factor receptor (EGFR), an important receptor of tyrosine kinases in animals, plays versatile functions in development, growth and tissue regeneration. In this study, we described the characteristic and function of an EGFR gene identified from P. f. martensii (PmEGFR). PmEGFR possesses a typical EGFR structure and is expressed in all studied tissues, with the highest expression level in adductor muscle. PmEGFR expression level is significantly higher in the fast-growing group than that in the slow-growing one. Correlation analysis represents that shell height and shell weight show positive correlation with PmEGFR expression (p < 0.05), and total weight and tissue weight exhibit positive correlation with it (p < 0.01). This study indicates that PmEGFR is a valuable functional gene associated with growth traits.