Zhe Zheng

Identification and Characterization of MicroRNAs in Pearl Oyster Pinctada martensii by Solexa Deep Sequencing

MicroRNAs (miRNAs) are short-nucleotide RNA molecules that function as negative regulators of gene expression in various organisms. However, miRNAs of Pinctada martensii have not been reported yet. P. martensii is one of the main species cultured for marine pearl production in China and Japan. In order to obtain the repertoire of miRNAs in P. martensii, we constructed and sequenced small RNA libraries prepared from P. martensii by Solexa deep sequencing technology and got a total of 27,479,838 reads representing 3,176,630 distinct sequences. After removing tRNAs, rRNAs, snRNAs, and snoRNAs, 10,596,306 miRNA reads representing 18,050 distinct miRNA reads were obtained. Based on sequence similarity and hairpin structure prediction, 258 P. martensii miRNAs (pm-miRNA) were identified. Among these pm-miRNAs, 205 were conserved across the species, whereas 53 were specific for P. martensii. The 3′ end sequence of U6 snRNA was obtained from P. martensii by 3′ rapid amplification of cDNA end PCR reaction and sequence-directed cloning. Eight conserved pm-miRNAs and two novel pm-miRNAs were validated by stem-loop quantitative real-time PCR with U6 snRNA as an internal reference gene. pm-miRNAs and the reported biomineralization-related genes were subjected to target analysis by using target prediction tools. Some of the pm-miRNAs, such as miR-2305 and miR-0046, were predicted to participate in biomineralization by regulating the biomineralization-related genes. Thus, this study demonstrated a large-scale characterization of pm-miRNAs and their potential function in biomineralization, providing a foundation to understand shell formation.

The pearl oyster Pinctada fucata martensii genome and multi-omic analyses provide insights into biomineralization

Abstract Nacre, the iridescent material found in pearls and shells of molluscs, is formed through an extraordinary process of matrix-assisted biomineralization. Despite recent advances, many aspects of the biomineralization process and its evolutionary origin remain unknown. The pearl oyster Pinctada fucata martensii is a well-known master of biomineralization, but the molecular mechanisms that underlie its production of shells and pearls are not fully understood. We sequenced the highly polymorphic genome of the pearl oyster and conducted multi-omic and biochemical studies to probe nacre formation. We identified a large set of novel proteins participating in matrix-framework formation, many in expanded families, including components similar to that found in vertebrate bones such as collagen-related VWA-containing proteins, chondroitin sulfotransferases, and regulatory elements. Considering that there are only collagen-based matrices in vertebrate bones and chitin-based matrices in most invertebrate skeletons, the presence of both chitin and elements of collagen-based matrices in nacre suggests that elements of chitin- and collagen-based matrices have deep roots and might be part of an ancient biomineralizing matrix. Our results expand the current shell matrix-framework model and provide new insights into the evolution of diverse biomineralization systems.

Computational prediction of candidate miRNAs and their potential functions in biomineralization in pearl oyster Pinctada martensii

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.

MicroRNA, Pm-miR-2305, Participates in Nacre Formation by Targeting Pearlin in Pearl Oyster Pinctada martensii

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.

Mantle Branch-Specific RNA Sequences of Moon Scallop Amusium pleuronectes to Identify Shell Color-Associated Genes.

Amusium pleuronectes (Linnaeus) that secretes red- and white-colored valves in two branches of mantle tissues is an excellent model for shell color research. High-throughput transcriptome sequencing and profiling were applied in this project to reveal the detailed molecular mechanism of this phenotype differentiation. In this study, 50,796,780 and 54,361,178 clean reads were generated from the left branch (secreting red valve, RS) and right branch (secreting white valve, WS) using the Illumina Hiseq 2000 platform. De novo assembly generated 149,375 and 176,652 unigenes with an average length of 764 bp and 698 bp in RS and WS, respectively. Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathway analysis indicated that the differentially expressed genes were involved in 228 signaling pathways, and 43 genes were significantly enriched (P<0.01). Nineteen of 20 differentially expressed vitellogenin genes showed significantly high expression in RS, which suggested that they probably played a crucial role in organic pigment assembly and transportation of the shell. Moreover, 687 crystal formation-related (or biomineralization-related) genes were detected in A. pleuronectes, among which 144 genes exhibited significant difference between the two branches. Those genes could be classified into shell matrix framework participants, crystal nucleation and growth-related elements, upstream regulation factors, Ca level regulators, and other classifications. We also identified putative SNP and SSR markers from these samples which provided the markers for genetic diversity analysis, genetic linkage, QTL analysis. These results provide insight into the complexity of shell color differentiation in A. pleuronectes so as valuable resources for further research.

miR-29a Participated in Nacre Formation and Immune Response by Targeting Y2R in Pinctada martensii

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.

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

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

PmRunt regulated by Pm-miR-183 participates in nacre formation possibly through promoting the expression of collagen VI-like and Nacrein in pearl oyster Pinctada martensii

Heterodimeric PEBP2/CBFs are key regulators in diverse biological processes, such as haematopoietic stem-cell generation, bone formation and cancers. In this work, we cloned runt-like transcriptional factor (designated as PmRunt) and CBF β (designated as PmCBF) gene, which comprise the heterodimeric transcriptional factor in Pinctada martensii. PmRunt was identified with an open reading frame that encodes 545 amino acids and has typical Runt domain. Phylogenetic analysis results speculated that runt-like transcriptional factors (RDs) in vertebrates and invertebrates are separated into two branches. In molluscs, PmRunt and other RDs are clustered in one of these branches. Direct interaction between PmRunt and PmCBF was evidenced by yeast two-hybrid assay results. Gene repression by RNA interference decreased the expression level of PmRunt, and subsequent observation of the inner surface of the nacre by scanning electron microscopy demonstrated disordered growth. The luciferase activities of reporters that contain promoter regions of Collagen VI-like (PmColVI) and PmNacrein were enhanced by PmRunt. Meanwhile, Pm-miR-183 apparently inhibited the relative luciferase activity of reporters containing the 3′-UTR of PmRunt. The expression level of PmRunt was repressed after Pm-miR-183 was overexpressed in the mantle tissue. Therefore, we proposed that PmRunt could be targeted by Pm-miR-183 and regulate the transcription of PmColVI and PmNacrein by increasing their transcriptional activity, thereby governing nacre formation.

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

ABSTRACT 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. HIGHLIGHTSHSP70 and HSP40 exhibited the polar expression model in the two populations.Biosynthesis of unsaturated fatty acid and energy release was increased in the SL.The small molecular system with amino acids might play a key role in cold stress.Lysosome and immune function could influence the low temperature resistance.

Transcriptome analysis of growth heterosis in pearl oyster Pinctada fucata martensii

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