Next Article in Journal
Mitochondrial Dysfunction in Aging and Diseases of Aging
Previous Article in Journal
Plasticity in Standard and Maximum Aerobic Metabolic Rates in Two Populations of an Estuarine Dependent Teleost, Spotted Seatrout (Cynoscion nebulosus)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 8
Issue 2
10.3390/biology8020047
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata)

by
Songqian Huang
1,
Yuki Ichikawa
1,
Kazutoshi Yoshitake
1,
Shigeharu Kinoshita
1,
Yoji Igarashi
1,
Fumito Omori
2,
Kaoru Maeyama
2,
Kiyohito Nagai
3,
Shugo Watabe
4 and
Shuichi Asakawa
1,*
1
Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo 113-8657, Japan
2
Mikimoto Pharmaceutical CO., LTD., Kurose 1425, Ise, Mie 516-8581, Japan
3
Pearl Research Laboratory, K. MIKIMOTO & CO., LTD., Osaki Hazako 923, Hamajima, Shima, Mie 517-0403, Japan
4
School of Marine Biosciences, Kitasato University, Minami-ku, Sagamihara, Kanagawa 252-0313, Japan
*
Author to whom correspondence should be addressed.
Biology 2019, 8(2), 47; https://doi.org/10.3390/biology8020047
Submission received: 18 April 2019 / Revised: 12 June 2019 / Accepted: 13 June 2019 / Published: 17 June 2019
(This article belongs to the Section Marine Biology)
Browse Figures
Versions Notes

Abstract

:
The biological process of pearl formation is an ongoing research topic, and a number of genes associated with this process have been identified. However, the involvement of microRNAs (miRNAs) in biomineralization in the pearl oyster, Pinctada fucata, is not well understood. In order to investigate the divergence and function of miRNAs in P. fucata, we performed a transcriptome analysis of small RNA libraries prepared from adductor muscle, gill, ovary, and mantle tissues. We identified 186 known and 42 novel miRNAs in these tissues. Clustering analysis showed that the expression patterns of miRNAs were similar among the somatic tissues, but they differed significantly between the somatic and ovary tissues. To validate the existence of the identified miRNAs, nine known and three novel miRNAs were verified by stem-loop qRT-PCR using U6 snRNA as an internal reference. The expression abundance and target prediction between miRNAs and biomineralization-related genes indicated that miR-1990c-3p, miR-876, miR-9a-3p, and novel-3 may be key factors in the regulatory network that act by controlling the formation of matrix proteins or the differentiation of mineralogenic cells during shell formation in mantle tissue. Our findings serve to further clarify the processes underlying biomineralization in P. fucata.

1. Introduction

The pearl oyster (Pinctada fucata) is a well-studied organism, owing to the economic potential of pearl production as well as the fascinating biology of mollusks. P. fucata is also a representative experimental model for biomineralization analysis [1,2]. A number of genes involved in biomineralization have been identified, and their functions in pearl and shell formation have been clarified [3,4,5,6,7]. Recently, transcriptomics [8,9], proteomics [10,11,12], genomics [1,2,13,14], and gene interference techniques [15,16] have been used to investigate genetic components of shell and pearl formation in mollusks.
MicroRNAs (miRNAs) are a class of small endogenous noncoding RNAs, 20–25 nucleotides (nt) in length. They are embedded within the stem regions of hairpin transcripts that exist in a wide range of invertebrates and vertebrates [17]. A single miRNA can regulate hundreds of target genes. Up to 30% of human protein coding genes may be regulated by miRNAs [18]. miRNAs play a vital role in the regulation of gene expression at the post-transcriptional level, especially for signaling pathways involved in cellular development, proliferation, apoptosis, oncogenesis, and differentiation [19,20]. miRNAs negatively regulate gene expression through sequence-specific interactions with the 3′ untranslated regions (UTRs) of target genes, and thereby cause translational repression or mRNA destabilization [21,22]. Prior research has shown that some miRNAs, such as miR-223, miR-125b and, miR-302a, participate in the control of biomineralization during bone formation in animals [23,24,25].
Mollusks are the second largest phylum in Metazoa, yet very few studies have addressed the diversity and function of miRNAs in mollusks. So far, five miRNAs have been identified from Holiotis rufescens, 60 miRNAs from Lottia gigantean [26], 258 miRNAs from Pinctada martensii [27], 199 miRNAs from Crassostrea gigas [28], and 289 miRNAs from Crassostrea hongkongensis [29]. Moreover, a computational prediction of candidate miRNAs was performed for the Patella vulgata [30] and P. martensii [31] genomes. In addition, prior work revealed that miR-2305 participates in nacre formation by targeting pearlin in P. martensii [32]. However, data regarding the involvement of miRNAs in pearl formation remain limited.
To obtain various miRNAs and explore their potential functions in biomineralization, we constructed and sequenced eight small RNA libraries prepared from adductor muscle, gill, mantle, and ovary tissues of two P. fucata. The genomic and transcript sequences of P. fucata are very helpful for performing functional analysis of miRNAs in this organism. The clarification of the functions of miRNAs in P. fucata will provide new insights to further understand the molecular mechanisms of biomineralization. In addition, this will provide useful information for further research into various biological processes in P. fucata. The aims of this study were (1) to identify known and novel miRNAs of P. fucata, (2) to clarify the expression profile of identified miRNAs in the tissues examined, (3) to predict genes targeted by the miRNAs, and (4) to identify the miRNAs involved in biomineralization.

2. Materials and Methods

2.1. Ethics

This study was conducted in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the University of Tokyo. All efforts were made to minimize the suffering of the animals.

2.2. Small RNA Isolation and cDNA Library Construction

Two female pearl oysters (approximately two years old) were collected from Mikimoto Pearl Research Institution Base, Mie Prefecture, Japan in 2014. Adductor muscle (Ad for short), gill (Gi for short), mantle (Ma for short), and ovary (Ov for short) tissues were collected from the pearl oysters. Eight RNA samples in total were extracted from the tissues, and small RNA extraction was performed following the protocol of the mirVana miRNA Isolation Kit (Life Technologies, California, USA). In brief, the sequencing libraries were constructed as per the Ion Total RNA-seq Kit v2 small RNA libraries construction protocol (Life Technologies, California, USA). This was followed by adapter ligation, synthesis of cDNA by reverse transcription, purification of the small RNA fraction, and amplification by PCR to construct cDNA libraries. An Agilent 2200 TapeStation (Agilent Technologies, Waldbronn, Germany) was used to determine the cDNA concentration. The cDNA libraries were then used to prepare a template for loading onto the analysis chip, according to the Ion PI™ Template OT2 200 Kit v3 protocol (Life Technologies, California, USA). Sequencing was then performed following the manufacturer’s protocol for the Ion PITM Sequencing 200 Kit v3 (Ion Proton Sequencer and Ion Proton Server, Life Technologies).

2.3. Small RNA Sequence Analysis

The raw sequencing data were saved as FASTQ files. We first filtered out the reads with low quality, containing poly-A tails, and those outside the range of 15–31 nt length (Step 1). After filtering, quality and length distribution statistics for these sequencing data were examined. Then, the filtered reads were evaluated by a homology search using basic local alignment search tool (BLAST) against the NCBI non-coding RNA database [33] and Rfam database [34] to remove rRNA, tRNA, snRNA, and snoRNA sequences (Step 2). Subsequently, the filtered reads were mapped onto the P. fucata genome by bowtie (-v1) (Step 3) [35]. Then, miRNA identification was performed by comparing the sequences with the known mature miRNAs and miRNA precursors in miRBase 22.0 [36]. Unannotated sequences that failed to match with the sequences in the databases above were analyzed by miRDeep2 to predict novel miRNA candidates [37]. The base bias composition at each nucleotide position of all identified miRNAs was analyzed. The genome of P. fucata (version 2.0) [2] was used as the reference genome for miRNA precursor prediction.

2.4. Analysis of Differentially Expressed miRNAs among Tissues

The total miRNA reads of each library were used for normalization to analyze differentially expressed miRNAs between two tissues’ libraries. The miRNA counts were normalized according to the reads per million reads (RPM) method. The expression levels of known and novel miRNAs in the various tissues were compared to identify differentially expressed miRNAs. DESeq2 software (version 1.14.1) [38] was used to analyze the differences in the expression levels of the miRNAs between two tissues with the criteria: |log2(FoldChange)| > 1 and p-adj < 0.05. Six paired comparisons were evaluated by DESeq2, including Ad vs. Gi, Ad vs. Ov, Ad vs. Ma, Gi vs. Ov, Gi vs. Ma, and Ov vs. Ma. Each tissue had two individuals as replicates.

2.5. Target Analysis for miRNAs and Biomineralization-Related Genes

Complete sequences of 105 biomineralization-related genes from the Pinctada genus were downloaded from NCBI GenBank [39], and the 3′ UTRs were used for miRNA-target prediction. Target predictions for miRNAs and biomineralization-related genes were analyzed by RNAhybrid [40], miRanda [41] and RNA22 [42]. For the RNAhybrid algorithm, the free energy (e) and p-value (p) were used to estimate the interaction between miRNAs and mRNAs. The conditions for functional targets were e < −20 kcal/mol and p < 0.05. Functional target genes of miRNAs identified using the miRanda algorithm were restricted by the conditions of score > 100 and free energy < −10 kcal/mol. The RNA22 algorithm was used to determine the most favorable hybridization site between miRNA and mRNA with a sensitivity of 63%, specificity of 61%, and seed size of 7, allowing a maximum of one unpaired base in the seed.

2.6. Validation of P. fucata miRNAs by Stem-Loop qRT-PCR

Stem-loop qRT-PCR analysis [43] was employed to validate and determine the specific expression of miRNAs in P. fucata. Ten adult female P. fucata were obtained for tissue collection from Mikimoto Pearl Research Institution Base in May 2018. Total RNA was isolated from adductor muscle, gill, ovary, and mantle tissues using an RNeasy Mini Kit (QIAGEN, Maryland, USA). RNA quality and purity were assessed using RNA ScreenTape (Agilent, California, USA). Each RNA sample (0.5 μg) was reverse transcribed into cDNA using a cDNA Synthesis Kit (TaKaRa, Shiga, Japan) as per the manufacturer’s instructions. The stem-loop primer was used for small RNA reverse transcription; qRT_PCR forward primer and common reverse primer were used to amplify small RNA sequences. PCR assays were carried out in a quantitative thermal cycler (Applied Biosystems 7300 Real-Time PCR System) (Life Technologies) in a 20 μL reaction volume containing 10 μL SYBR Premix Ex Taq II (TaKaRa), 2 μL cDNA, 0.4 μL ROX Reference Dye (50×) (TaKaRa), 6 μL sterile distilled water, and 0.4 μM of each primer. Primer sequences are listed in Table 1. The thermal program included 95 °C for 30 s, 40 cycles of 95 °C for 5 s, 60 °C for 31 s, and a dissociation stage of 95 °C for 15 s, 60 °C for 60 s and 95 °C for 15 s. The expression of each miRNA quantified by stem-loop qRT-PCR was presented as n = 10 done for three replicates, and U6 snRNA was used as the internal reference for the stem-loop qRT-PCR. Ct values were represented by the mean values of three independent replicates, and the relative expression levels were calculated using the 2−∆∆Ct method [44].

2.7. Data Accessibility

The raw sequencing reads obtained in the present study are deposited at DNA Data Bank of Japan (DDBJ) Sequence Read Archive (DRA) under submission accession number DRA006953.

3. Results

3.1. Overview of Small RNA Library Sequencing

Eight small RNA libraries from the adductor muscle (Ad), gill (Gi), ovary (Ov), and mantle (Ma) tissues of two individual P. fucata were sequenced by Ion Proton sequencing to survey the miRNA diversity in P. fucata. In total, 50.32 million raw reads were obtained from all libraries (Table 2). After the removal of low-quality sequences, adaptor sequences, poly-A tail sequences, and sequences shorter than 15 nt or longer than 31 nt, 41.72 million reads remained for the statistical analysis of sequence length (Figure 1 and Table 2, Step 1). The length distribution of each tissue library is shown in Figure S1. Two obvious small RNA peaks with different lengths were found in the sequencing data. After removing known RNAs and genome unmapped reads, 33.26 million reads (66.1% of raw reads) were remaining for miRNA identification in this study.

3.2. Identification of Known and Novel miRNAs

To identify known miRNAs in P. fucata, we compared our dataset with the known miRNAs in miRBase 22.0. First, the remaining sequencing reads were mixed for miRNA identification. A read found more than five times was considered as a possible miRNA in P. fucata. Allowing no more than one mismatch between sequences, 186 known miRNAs were identified (Table S1). We then analyzed the abundance of the known miRNA, which revealed a large distribution in the expression levels. The read numbers of known miRNAs ranged from 1 to 1359669 in a single library, indicating that not only highly expressed miRNAs but also weakly expressed miRNAs were present in P. fucata.
One of the important characteristics that distinguish miRNAs from other endogenous small RNAs is the capacity of the precursor miRNA sequences to fold back into a canonical stem-loop hairpin structure [27]. The genomic data of P. fucata were used for miRNA precursor identification and were further parsed through the miRDeep2 software for the prediction of precursor sequences and secondary structures. A total of 280 sequences could form proper secondary hairpins and were considered miRNA precursors. Moreover, 186 known and 42 novel miRNAs were identified based on expression levels and precursor identification. The base bias at each position of all identified miRNAs is shown in Figure S2. The majority of miRNAs tended to start with 5′-U, which is consistent with typical miRNA sequence patterns [17]. Of these miRNAs, 38 (32 known and six novel miRNAs) have multiple (from two to five) precursor sequences (Table S1). For example, miR-9649 and novel-18 have five and three precursors, respectively.

3.3. Expression Levels of miRNAs in P. fucata

Among 228 miRNAs identified, 165, 198, 159 and 185 miRNAs were found in adductor muscle, gill, ovary and mantle tissues, respectively. Of these, 119 miRNAs (52.19% of total miRNA) were ubiquitously expressed in all examined tissues. On the contrary, 19 miRNAs (11.95% of ovary miRNAs) and eight miRNAs (4.04% of gill miRNAs) were uniquely expressed in ovary and gill tissues, respectively (Figure S3). One miRNA, miR-8871, was also uniquely expressed in mantle tissues, albeit with a low expression level. The expression levels of P. fucata miRNAs were calculated by reads per million reads (RPM), and these data are shown in Table S2. The novel miRNAs were expressed at lower levels than the previously known miRNAs (Figure 2).
Among the miRNAs identified in P. fucata, miR-279 showed the highest expression levels in gill and mantle tissues, whereas miR-100 and miRNA-549a showed the highest expression levels in adductor muscle and ovary tissues. Let-7a, miR-184-3p, miR-71, and miR-87a also demonstrated high expression in all examined tissues in P. fucata, while miR-133a, miR-549a, and miR-1990c-3p were highly expressed in adductor muscle, ovary and mantle tissues, respectively. The expression patterns of miRNAs in ovary tissue showed rather different features from those in the somatic tissues (Figure 3). In Figure 3a, four clusters of miRNAs are shown based on their expression patterns, in which each cluster corresponds to highly and specifically expressing miRNA groups in each tissue based on the average expression levels of two replicates. The expression pattern of ovary tissue differed from somatic tissues, and 24 (57.14% of novel miRNAs) novel miRNAs were highly expressed in ovary tissue. Essentially identical clustering patterns were observed based on individual expression levels (Figure 3b,c). However, when we examined clustering patterns of the eight samples (four tissues of the two individuals), mixed and complicated patterns were observed in somatic and ovary tissues (Figure 3d). Then, we examined clustering patterns of the six somatic samples, in which a pair for each tissue clustered together, indicating that the expression of miRNAs in the somatic tissues showed similar patterns between the two individuals (Figure 3e). The principal component analysis (PCA) also showed large differences in miRNA expressions between ovary and somatic tissues (Figure S4). Data regarding the miRNA expressions as obtained from high-throughput sequencing are shown in Table S2. For mantle tissues, eight miRNAs, including miR-1493, miR-1993, miR-9a-3p, miR-9b-3p, novel-3, miR-876, miR-1990c-3p, and miR-10a, were significantly highly expressed when compared with the other three tissues (Figure S5).

3.4. Validation of miRNA by Stem-Loop qRT-PCR

The total RNA extracted from adult P. fucata was subjected to stem-loop qRT-PCR to validate the authenticity of miRNA expression shown by the transcriptome analysis. Twelve miRNAs, including nine known miRNAs and three novel miRNAs, were selected for the checking of miRNA expression. Of those, six miRNAs, namely miR-1493, miR-1990c-3p, miR-1993, miR-876, miR-9a-3p, and novel-3, were selected because of their high expression in mantle tissues. The other miRNAs were selected randomly. All of these miRNAs could be readily detected by stem-loop qRT-PCR. Except for inconsistency in the expression of miR-1493 and miR-1993 in adductor muscle between the transcriptome and qRT-PCR data, the other ten miRNAs demonstrated consistent levels of relative expression by qRT-PCR when compared with those obtained from the high-throughput sequencing data (Figure 4). Let-7a, novel-1, miR-279, and miR-200a were ubiquitously expressed in all examined tissues, while miR-279 and miR-200a were highly expressed in gill, ovary, and mantle tissues. Among these miRNAs, miR-1990c-3p, miR-876, miR-9a-3p, and novel-3 were highly expressed in mantle tissues, while miR-183 and novel-10 were highly expressed in ovary tissue when compared with the other three tissues.

3.5. Functional Prediction of miRNA in Biomineralization

To elucidate the functions of the identified miRNAs in biomineralization, the putative target associations between miRNAs and the biomineralization-related genes were analyzed using three informatics software packages, RNAhybrid, miRanda, and RNA22, which were programmed to predict the miRNA:mRNA interaction sites. One hundred and five biomineralization-related genes from the Pinctada genus with completed 3′UTR sequences were downloaded from the NCBI database. As shown in Table S3, 248, 323, and 577 miRNA:mRNA interaction sites were predicted by RNAhybrid, miRanda, and RNA22, respectively. The results also demonstrated that the majority of these genes could be targeted by more than one miRNA. For example, Fam20c could be targeted by 27 miRNAs, of which twenty-five were known miRNAs and two were novel miRNAs. The predicted interaction sites for the same miRNA:mRNA pair varies widely among the software packages used. For example, when RNA22 was employed, the MPN88 family genes were predicted to be regulated by miR-872 (Figure 5a), while RNAhybrid predicted that MPN88 family genes were also regulated by miR-463-3p (Figure 5b).
Of these predictions, twenty-two genes are shown to be regulated by multiple miRNAs by all three programs (Table S3). On the other hand, no interaction sites were predicted in eight genes by all three programs. Moreover, some miRNAs are predicted to regulate multiple genes, whereas some genes are predicted to be regulated by multiple miRNAs. For example, while let-7a was predicted to regulate 13 genes, including BMP3/7, CaM, Cathepsin-B, GRP, Nacrein, PfCHS1, PFMG6, PFMG7, Pif177-like, and shematrin-5/6/8, the gene BMP3, a subgroup of the bone morphogenetic protein (BMP) family, could be regulated by either let-7a, miR-133a, or miR-9007. Among these potential target genes, most of them are shell matrix proteins, such as GRP, KRMP, Linkine, Pif177-like, PFMG, Prismalin-14, SGMP1, Shematrins etc., which constitute the basic calcium carbonate and organic matrices formed during the biomineralization process. Moreover, many protein kinases (ALP, CaM, and MMP) and transcription factors (OCT4) are also regulated by miRNA according to target site predictions (Table S3).
miR-1990c-3p, miR-876, miR-9a-3p, and novel-3 were highly expressed in mantle tissues, and each of these miRNAs has multiple target genes related to biomineralization. Therefore, these miRNAs are very likely to be involved in the biomineralization process. Indeed, miR-1990c-3p was predicted to regulate MMP and SGMP1; miR-9a-3p regulated Pif177-like; miR-876 regulated Fam20c, GRP, PMMG1, Prismalin-14, shematrin-2/2β, and Tyr-1; and novel-3 regulated GRP, Linkine, MMP, and OCT4 (Figure 6).

4. Discussion

Prior research has demonstrated that miRNAs are a major class of post-transcriptional regulatory molecules that have critical functions in diverse biological processes [45]. Abundant miRNAs have been identified in various organisms, while only a small number of miRNAs have been identified and undergone functional prediction in mollusks such as P. fucata [31]. In the present study, two main peaks of different sizes of small RNAs were observed in the distribution of small RNA lengths by high-throughput sequencing, with small RNAs with a length of 21–23 nt likely representing miRNAs, while those 29–31 nt in length mainly represented putative Piwi-interacting RNAs (piRNA) [46]. Separate analyses were conducted on these two different sets of small RNAs (we focused on the miRNA fraction in this study).
A total of 228 miRNAs, including 186 known and 42 novel miRNAs, were identified by high-throughput sequencing in P. fucata. To gain a better understanding of miRNA expression in P. fucata, the normalized RPM of each miRNA was used to compare miRNA expression among tissues. More than 50% of miRNAs were commonly expressed in all tissues, while a few miRNAs were specifically expressed in a single tissue. Clustering analysis and PCA of miRNA expression patterns, based on average expression levels, clustered the somatic tissues together, indicating that the miRNA expression patterns in P. fucata differ between the somatic and ovary tissues. These results demonstrate that the expression of miRNAs in P. fucata is both common and tissue-specific, which may play an important role in development, differentiation, and biological processes. The results of the eight samples in Figure 3d showed mixed and complicated clustering patterns. This result might be caused by the perturbation of relatively differential expression levels of miRNAs in ovaries and mantles between the two individuals.
Among these 228 miRNAs, many known miRNAs were highly expressed in P. fucata. miR-279 has been observed as being well preserved in arthropods, but is rarely observed in vertebrates [47]. Previous studies showed that the miR-279 family is involved in the proper development of CO2 sensory neurons and the maintenance of circadian rhythm in Drosophila [48,49]. miR-100 was also ubiquitously expressed in all examined tissues by high-throughput sequencing (Table S2). In previous studies, it was considered a tumor suppressor miRNA in patient tissue, and the decreased expression of miR-100 appears to play an important role in the development and progression of disease and may contribute to reduced sensitivity to ionizing radiation [50,51]. Since no predicted interaction site was identified between these two miRNAs and the biomineralization-related genes (Table S3), miR-279 and miR-100 may be involved outside of the biomineralization process in P. fucata. Let-7a was also found in most animals, suggesting that these miRNAs, which regulate cell proliferation and differentiation, are evolutionarily well conserved [52,53].
In addition, 42 novel miRNAs were identified in P. fucata. The expression density of novel miRNAs was lower than that of the majority of the known miRNAs. This result is consistent with previous observations suggesting that novel miRNAs are often expressed at lower levels than known miRNAs [27]. Many of the novel miRNAs (57.14%) are likely specialized for reproductive functions given their tissue-specific biases in expression. This finding indicates that these novel miRNAs may have not general but specific functions. Further experiments could provide further insight into the functions of these novel miRNAs in P. fucata.
To verify the authenticity of miRNAs identified by Ion Proton sequencing, we performed stem-loop qRT-PCR to detect and measure the relative expression levels of miRNAs in P. fucata. All of the selected miRNAs were easily detected by stem-loop qRT-PCR. Most of the miRNA expression levels obtained by stem-loop qRT-PCR are consistent with those measured by Ion Proton sequencing, except for miR-1493 and miR-1993 in adductor muscle. This inconsistency may be caused by stem-loop qRT_PCR primer amplification efficiency, sequencing errors during the complex sequencing protocols [27] or differences among individual P. fucata. Nevertheless, the overall results of sequencing technology and stem-loop qRT-PCR in our study are highly consistent, suggesting that expression analysis was largely able to reveal accurate miRNA expression levels in P. fucata.
The biomineralization process is very complex and precise, and the expression of each related protein is subject to fine regulation [31]. Therefore, it would be beneficial to gain a better understanding of the regulators involved in biomineralization. Generally, miRNAs are fine regulators that function by silencing the expression of endogenous genes. In the present study, hundreds of interaction sites between miRNAs and biomineralization-related genes were predicted by RNAhybrid, miRanda, and RNA22, which were widely used in many reports, especially in the initial research when there was little knowledge about the miRNAs in this species [54,55,56]. Many biomineralization-related genes, such as ALP, BMP, CaM, GRP, KRMP, MMP, N16, OCT4, Pif177, PFMG, Prismin, SGMP, shematrin, etc., are regulated by miRNA according to the target sites’ prediction. Shematrins are synthesized in the mantle edge and secreted into the prismatic layer of the shell, where the protein family is thought to provide a framework for calcification in P. fucata [5]. Matrix metalloproteinases (MMPs), which can be regulated by miR-1990c-3p, digest proteins to prepare the peptides that contribute to calcium carbonate crystallization [57]. The transcription factor Pf-POU3F4 (OCT4), which is predicted to be regulated by novel-3, regulates the expression of the matrix protein genes Aspein and Prismalin-14 [58]. These results suggest that miRNAs are involved in the biomineralization process, not only through targeting the shell matrix protein, but also through the regulation of protein kinases and transcription factors. Here, we provided an important set of resources for future studies of the miRNA-dependent regulation of expression in the biomineralization process that forms pearls.

5. Conclusions

A total of 228 miRNAs were identified in P. fucata by high-throughput sequencing in the present study. Novel miRNAs demonstrate weaker expression than previously known miRNAs. Clustering analysis showed that the expression patterns of miRNAs were similar within the somatic tissues, but differed significantly between the somatic and ovary tissues. Hundreds of potential target sites of miRNAs and target genes were detected using multiple software tools. miR-1990c-3p, miR-876, miR-9a-3p, and novel-3, which were highly expressed in mantle tissues, may play a core role in biomineralization by regulating the formation of matrix proteins or protein kinase and transcription factor genes. These results further clarify the range of miRNAs found in P. fucata and their possible functions in biomineralization in mollusks.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-7737/8/2/47/s1, Figure S1: Length distribution of filtered reads in P. fucata. Figure S2: The sequence composition of miRNAs in P. fucata displays a first position nucleotide bias for uracil (U). Figure S3: Venn diagram of the number of identified miRNAs in P. fucata. Figure S4: PCA of miRNA expression patterns in P. fucata. Figure S5. Highly expressed miRNAs in P. fucata mantle tissues identified by high-throughput sequencing. Table S1: Detailed information miRNAs identified in P. fucata by high-throughput sequencing. Table S2: Read counts and normalized RPM of miRNA expression in P. fucata by high-throughput sequencing. Table S3: miRNA functioning in biomineralization as predicted by RNAhybrid, miRanda, and RNA22.

Author Contributions

Conceptualization, S.W., S.A. and S.H.; Methodology, S.H. and S.A.; Software, K.Y. and Y.I. (Yoji Igarashi); Formal analysis, Y.I. (Yuki Ichikawa), S.K. and F.O.; Investigation, K.M. and K.N.; Resources, K.M. and K.N.; Data curation, K.Y., S.K. and S.H.; Writing—original draft preparation, S.H.; Writing—review & editing, S.H., Y.I. (Yuki Ichikawa), K.Y., S.K., Y.I. (Yoji Igarashi), F.O., K.M., K.N., S.W. and S.A.; Supervision and project administration, S.A.; Funding acquisition, S.H. and S.A.

Funding

This work was supported by the Japan Society for the Promotion of Science [Project number JP24248034] and a Research Fellowship of the Japan Society for the Promotion of Science for Young Scientists [Project number 18J13176].

Acknowledgments

The authors thank the Mikimoto Pearl Research Institution Base for rearing and supporting the experiment animals.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Takeuchi, T.; Kawashima, T.; Koyanagi, R.; Gyoja, F.; Tanaka, M.; Ikuta, T.; Shoguchi, E.; Fujiwara, M.; Shinzato, C.; Hisata, K.; et al. Draft genome of the pearl oyster Pinctada fucata: A platform for understanding bivalve biology. DNA Res. 2012, 19, 117–130. [Google Scholar] [CrossRef] [PubMed]
  2. Takeuchi, T.; Koyanagi, R.; Gyoja, F.; Kanda, M.; Hisata, K.; Fujie, M.; Goto, H.; Yamasaki, S.; Nagai, K.; Morino, Y.; et al. Bivalve-specific gene expansion in the pearl oyster genome: Implications of adaptation to a sessile lifestyle. Zool. Lett. 2016, 2, 3. [Google Scholar] [CrossRef] [PubMed]
  3. Samata, T.; Hayashi, N.; Kono, M.; Hasegawac, K.; Horitad, C.; Akerad, S. A new matrix protein family related to the nacreous layer formation of Pinctada fucata. FEBS Lett. 1999, 462, 225–229. [Google Scholar] [CrossRef]
  4. Suzuki, M.; Murayama, E.; Inoue, H.; Ozaki, N.; Tohse, H.; Kogure, T.; Nagasawa, H. Characterization of Prismalin-14, a novel matrix protein from the prismatic layer of the Japanese pearl oyster (Pinctada fucata). Biochem. J. 2004, 382, 205–213. [Google Scholar] [CrossRef] [PubMed]
  5. Yano, M.; Nagai, K.; Morimoto, K.; Miyamoto, H. Shematrin: A family of glycine-rich structural proteins in the shell of the pearl oyster Pinctada fucata. Comp. Biochem. Phys. B 2006, 144, 254–262. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, H.L.; Liu, S.F.; Ge, Y.J.; Liu, J.; Wang, X.Y.; Xie, L.P.; Zhang, R.Q.; Wang, Z. Identification and characterization of a biomineralization related gene PFMG1 highly expressed in the mantle of Pinctada fucata. Biochemistry 2007, 46, 844–851. [Google Scholar] [CrossRef] [PubMed]
  7. Takagi, R.; Miyashita, T. Prismin: A new matrix protein family in the Japanese pearl oyster (Pinctada fucata) involved in prismatic layer formation. Zool. Sci. 2010, 27, 416–426. [Google Scholar] [CrossRef]
  8. Fang, D.; Xu, G.; Hu, Y.; Pan, C.; Xie, L.P.; Zhang, R.Q. Identification of genes directly involved in shell formation and their functions in pearl oyster, Pinctada fucata. PLoS ONE 2011, 6, e21860. [Google Scholar] [CrossRef]
  9. Kinoshita, S.; Wang, N.; Inoue, H.; Maeyama, K.; Okamoto, K.; Nagai, K.; Kondo, H.; Hirono, I.; Asakawa, S.; Watabe, S. Deep sequencing of ESTs from nacreous and prismatic layer producing tissues and a screen for novel shell formation-related genes in the pearl oyster. PLoS ONE 2011, 6, e21238. [Google Scholar] [CrossRef]
  10. Joubert, C.; Piquemal, D.; Marie, B.; Manchon, L.; Pierrat, F.; Zanella-Cleon, I.; Cochennec-Laureau, N.; Gueguen, Y.; Montagnani, C. Transcriptome and proteome analysis of Pinctada margaritifera calcifying mantle and shell: Focus on biomineralization. BMC Genom. 2010, 11, 613. [Google Scholar] [CrossRef]
  11. Berland, S.; Marie, A.; Duplat, D.; Milet, C.; Sire, J.Y.; Bedouet, L. Coupling proteomics and transcriptomics for the identification of novel and variant forms of mollusk shell proteins: A study with P. margaritifera. Chembiochem 2011, 12, 950–961. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, C.; Li, S.; Kong, J.; Liu, Y.; Wang, T.; Xie, L.; Zhang, R.Q. In-depth proteomic analysis of shell matrix proteins of Pinctada fucata. Sci. Rep. 2015, 5, 17269. [Google Scholar] [CrossRef] [PubMed]
  13. Miyamoto, H.; Endo, H.; Hashimoto, N.; Limura, K.; Isowa, Y.; Kinoshita, S.; Kotaki, T.; Masaoka, T.; Miki, T.; Nakayama, S.; et al. The diversity of shell matrix proteins: Genome-wide investigation of the pearl oyster, Pinctada fucata. Zool. Sci. 2013, 30, 801–816. [Google Scholar] [CrossRef] [PubMed]
  14. Du, X.D.; Fan, G.Y.; Jiao, Y.; Zhang, H.; Guo, X.M.; Huang, R.L.; Zheng, Z.; Bian, C.; Deng, Y.W.; Wang, Q.H.; et al. The pearl oyster Pinctada fucata martensii genome and multi-omic analyses provide insights into biomineralization. Gigascience 2017, 6, gix059. [Google Scholar] [CrossRef] [PubMed]
  15. Suzuki, M.; Saruwatari, K.; Kogure, T.; Yamamoto, Y.; Nishimura, T.; Kato, T.; Nagasawa, H. An acidic matrix protein, Pif, is a key macromolecule for nacre formation. Science 2009, 325, 1388–1390. [Google Scholar] [CrossRef]
  16. Fang, D.; Pan, C.; Lin, H.; Lin, Y.; Zhang, G.Y.; Wang, H.Z.; He, M.X.; Xie, L.P.; Zhang, R.Q. Novel basic protein, PfN23, functions as key macromolecule during nacre formation. J. Biol. Chem. 2012, 287, 15776–15785. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, S.Q.; Cao, X.J.; Tian, X.C.; Wang, W.M. High-throughput sequencing identifies microRNAs from posterior intestine of loach (Misgurnus anguillicaudatus) and their response to intestinal air-breathing inhibition. PLoS ONE 2016, 11, e0149123. [Google Scholar] [CrossRef]
  18. Di, L.G.; Calin, G.A.; Croce, C.M. MicroRNAs: Fundamental facts and involvement in human diseases. Birth Defects Res. C 2006, 78, 180–189. [Google Scholar]
  19. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  20. Urbich, C.; Kuehbacher, A.; Dimmeler, S. Role of microRNAs in vascular diseases, inflammation, and angiogenesis. Cardiovasc. Res. 2008, 79, 581–588. [Google Scholar] [CrossRef] [Green Version]
  21. Allen, E.; Xie, Z.; Gustafson, A.M.; Carrington, J.C. MicroRNA-directed phasing during transacting siRNA biogenesis in plants. Cell 2005, 121, 207–221. [Google Scholar] [CrossRef] [PubMed]
  22. Bushati, N.; Cohen, S.M. MicroRNA functions. Annu. Rev. Cell Dev. Biol. 2007, 23, 175–205. [Google Scholar] [CrossRef] [PubMed]
  23. Sugatani, T.; Hruska, K.A. MicroRNA-223 is a key factor in osteoclast differentiation. J. Cell. Biochem. 2007, 101, 996–999. [Google Scholar] [CrossRef] [PubMed]
  24. Xia, Z.; Chen, C.; Chen, P.; Xie, H.; Luo, X. MicroRNAs and their roles in osteoclast differentiation. Front. Med. 2011, 5, 414–419. [Google Scholar] [CrossRef] [PubMed]
  25. Kang, I.H.; Jeong, B.C.; Hur, S.W.; Choi, H.; Choi, S.H.; Ryu, J.H.; Hwang, Y.C.; Koh, J.T. MicroRNA-302a stimulates osteoblastic differentiation by repressing COUP-TFII expression. J. Cell. Physiol. 2015, 230, 911–921. [Google Scholar] [CrossRef] [PubMed]
  26. Wheeler, B.M.; Heimberg, A.M.; Moy, V.N.; Sperling, E.A.; Holstein, T.W.; Heber, S.; Peterson, K.J. The deep evolution of metazoan microRNAs. Evol. Dev. 2009, 11, 50–68. [Google Scholar] [CrossRef]
  27. Jiao, Y.; Zheng, Z.; Du, X.D.; Wang, Q.H.; Huang, R.L.; Deng, Y.W.; Shi, S.L.; Zhao, X.X. Identification and characterization of miRNAs in pearl oyster Pinctada martensii by solexa sequencing. Mar. Biotechnol. 2014, 16, 54–62. [Google Scholar] [CrossRef]
  28. Zhou, Z.; Wang, L.; Song, L.; Liu, R.; Zhang, H.; Huang, M.; Chen, H. The identification and characteristics of immune-related microRNAs in haemocytes of oyster Crassostrea gigas. PLoS ONE 2014, 9, e88397. [Google Scholar] [CrossRef]
  29. Zhao, X.L.; Yu, H.; Kong, L.F.; Liu, S.K.; Li, Q. High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response. Sci. Rep. 2016, 6, 22687. [Google Scholar] [CrossRef]
  30. Kenny, N.J.; Namigai, E.K.; Marletaz, F.; Hui, J.H.; Shimeld, S.M. Draft genome assemblies and predicted microRNA complements of the intertidal lophotrochozoans Patella vulgata (Mollusca, Patellogastropoda) and Spirobranchus (Pomatoceros) lamarcki (Annelida, Serpulida). Mar. Genom. 2015, 24, 139–146. [Google Scholar] [CrossRef]
  31. Zheng, Z.; Jiao, Y.; Du, X.D.; Tian, Q.L.; Wang, Q.H.; Huang, R.L.; Deng, Y.W. Computational prediction of candidate miRNAs and their potential functions in biomineralization in pearl oyster Pinctada martensii. Saudi J. Biol. Sci. 2016, 23, 372–378. [Google Scholar] [CrossRef] [PubMed]
  32. Jiao, Y.; Zheng, Z.; Tian, R.R.; Du, X.D.; Wang, Q.H.; Huang, R.L. Pm-miR-2305, participates in nacre formation by targeting pearlin in pearl oyster Pinctada martensii. Int. J. Mol. Sci. 2015, 16, 21442–21453. [Google Scholar] [CrossRef] [PubMed]
  33. Szymanski, M.; Erdmann, V.A.; Barciszewski, J. Noncoding RNAs database (ncRNAdb). Nucleic Acids Res. 2007, 35, D162–D164. [Google Scholar] [CrossRef] [PubMed]
  34. Kalvari, I.; Argasinska, J.; Quinones-Olvera, N.; Nawrocki, E.P.; Rivas, E.; Eddy, S.R.; Bateman, A.; Finn, R.D.; Petrov, A.I. Rfam 13.0: Shifting to a genome-centric resource for non-coding RNA families. Nucleic Acids Res. 2018, 46, D335–D342. [Google Scholar] [CrossRef] [PubMed]
  35. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, R25. [Google Scholar] [CrossRef] [PubMed]
  36. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. MiRBase: From microRNA sequences to function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef] [PubMed]
  37. Friedlander, M.R.; Mackowiak, S.D.; Li, N.; Chen, W.; Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40, 37–52. [Google Scholar] [CrossRef]
  38. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  39. Benson, D.A.; Cavanaugh, M.; Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2017, 41, D36–D42. [Google Scholar] [CrossRef]
  40. Rehmsmeier, M.; Steffen, P.; Hochsmann, M.; Giegerich, R. Fast and effective prediction of microRNA/target duplexes. RNA 2014, 10, 1507–1517. [Google Scholar] [CrossRef]
  41. Enright, A.J.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D.S. MicroRNA targets in Drosophila. Genome Biol. 2003, 5, R1. [Google Scholar] [CrossRef] [PubMed]
  42. Loher, P.; Rigoutsos, I. Interactive exploration of RNA22 microRNA target predictions. Bioinformatics 2012, 28, 3322–3323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Chen, C.; Ridzon, D.A.; Broomer, A.J.; Zhou, Z.; Lee, D.H.; Nguyen, J.T.; Barbisin, M.; Xu, N.L.; Mahuvakar, V.R.; Andersen, M.R.; et al. Real-time quantification of micro-RNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33, e179. [Google Scholar] [CrossRef] [PubMed]
  44. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2011, 29, e45. [Google Scholar] [CrossRef] [PubMed]
  45. Obernosterer, G.; Leuschner, P.J.; Alenius, M.; Martinez, J. Post-transcriptional regulation of microRNA expression. RNA 2006, 12, 1161–1167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Huang, S.Q.; Ichikawa, Y.; Igarashi, Y.; Yoshitake, K.; Kinoshita, S.; Omori, F.; Maeyama, K.; Nagai, K.; Watabe, S.; Asakawa, S. Piwi-interacting RNA (piRNA) expression patterns in pearl oyster (Pinctada fucata) somatic tissues. Sci. Rep. 2019, 9, 247. [Google Scholar] [CrossRef] [PubMed]
  47. Sun, K.; Jee, D.; de Navas, L.F.; Duan, H.; Lai, E.C. Multiple in vivo biological processes are mediated by functionally redundant activities of drosophila mir-279 and mir-996. PLoS Genet. 2015, 11, e1005245. [Google Scholar] [CrossRef] [PubMed]
  48. Cayirlioglu, P.; Kadow, I.G.; Zhan, X.; Okamura, K.; Suh, G.S.; Gunning, D.; Lai, E.C.; Zipursky, S.L. Hybrid neurons in a microRNA mutant are putative evolutionary intermediates in insect CO2 sensory systems. Science 2008, 319, 1256–1260. [Google Scholar] [CrossRef]
  49. Luo, W.; Sehgal, A. Regulation of circadian behavioral output via a microRNA-JAK/STAT circuit. Cell 2012, 148, 765–779. [Google Scholar] [CrossRef]
  50. Henson, B.J.; Bhattacharjee, S.; O’Dee, D.M.; Feingold, E.; Gollin, S.M. Decreased expression of miR-125b and miR-100 in oral cancer cells contributes to malignancy. Gene Chromosome Cancer 2009, 48, 569–582. [Google Scholar] [CrossRef] [Green Version]
  51. Giangreco, A.A.; Vaishnav, A.; Wagner, D.; Finelli, A.; Fleshner, N.; Van der Kwast, T.; Vieth, R.; Nonn, L. Tumor suppressor microRNAs, miR-100 and -125b, are regulated by 1,25-dihydroxyvitamin D in primary prostate cells and in patient tissue. Cancer Prev. Res. 2013, 6, 483–494. [Google Scholar] [CrossRef] [PubMed]
  52. Johnson, C.D.; Esquela-Kerscher, A.; Stefani, G.; Byrom, M.; Kelnar, K.; Ovcharenko, D.; Wilson, M.; Wang, X.; Shelton, J.; Shingara, J.; et al. The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 2007, 67, 7713–7722. [Google Scholar] [CrossRef] [PubMed]
  53. Iwasaki, T.; Tanaka, K.; Kawano, M.; Itonaga, I.; Tsumura, H. Tumor-suppressive microRNA-let-7a inhibits cell proliferation via targeting of E2F2 in osteosarcoma cells. Int. J. Oncol. 2015, 46, 1543–1550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Rajewsky, N. MicroRNA target predictions in animals. Nat. Genet. 2006, 38, S8–S13. [Google Scholar] [CrossRef] [PubMed]
  55. Witkos, T.M.; Koscianska, E.; Krzyzosiak, W.J. Practical aspects of microRNA target prediction. Curr. Mol. Med. 2011, 1, 93–109. [Google Scholar] [CrossRef]
  56. Peterson, S.M.; Thompson, J.A.; Ufkin, M.L.; Sathyanarayana, P.; Liaw, L.; Congdon, C.B. Common features of microRNA target prediction tools. Front. Genet. 2014, 5, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kubota, K.; Tsuchihashi, Y.; Kogure, T.; Maeyama, K.; Hattori, F.; Kinoshita, S.; Sakuda, S.; Nagasawa, H.; Yoshimura, E.; Suzuki, M. Structural and functional analyses of a TIMP and MMP in the ligament of Pinctada fucata. J. Struct. Biol. 2017, 199, 216–224. [Google Scholar] [CrossRef]
  58. Gao, J.; Chen, Y.; Yang, Y.; Liang, J.; Xie, J.; Liu, J.; Li, S.G.; Zheng, G.L.; Xie, L.P.; Zhang, R.Q. The transcription factor Pf-POU3F4 regulates the expression of the matrix protein genes Aspein and Prismalin-14 in pearl oyster (Pinctada fucata). FEBS J. 2016, 283, 1962–1978. [Google Scholar] [CrossRef]
Biology 08 00047 g001 550
Figure 1. Length distribution of total filtered reads in P. fucata.
Figure 1. Length distribution of total filtered reads in P. fucata.
Biology 08 00047 g001
Biology 08 00047 g002 550
Figure 2. Expression density of known and novel microRNAs (miRNAs) in P. fucata. Ad: adductor muscle; Gi: gill; Ma: mantle; Ov: ovary.
Figure 2. Expression density of known and novel microRNAs (miRNAs) in P. fucata. Ad: adductor muscle; Gi: gill; Ma: mantle; Ov: ovary.
Biology 08 00047 g002
Biology 08 00047 g003 550
Figure 3. Heat map of expression levels of miRNAs in P. fucata somatic and ovary tissues, based on (a) average expression levels of the two individuals, (b) individual expression levels of individual 1, (c) individual expression levels of individual 2, (d) all examined samples of the two individuals, and (e) six somatic samples of the two individuals. The details regarding miRNA expression are shown in Table S2. RPM: reads per million reads. Ad: adductor muscle; Gi: gill; Ov: ovary; Ma: mantle.
Figure 3. Heat map of expression levels of miRNAs in P. fucata somatic and ovary tissues, based on (a) average expression levels of the two individuals, (b) individual expression levels of individual 1, (c) individual expression levels of individual 2, (d) all examined samples of the two individuals, and (e) six somatic samples of the two individuals. The details regarding miRNA expression are shown in Table S2. RPM: reads per million reads. Ad: adductor muscle; Gi: gill; Ov: ovary; Ma: mantle.
Biology 08 00047 g003
Biology 08 00047 g004 550
Figure 4. qRT-PCR validation of miRNAs expressed in P. fucata somatic and ovary tissues. Ad: adductor muscle; Gi: gill; Ov: ovary; Ma: mantle; RPM: reads per million reads; Seq: high-throughput sequencing.
Figure 4. qRT-PCR validation of miRNAs expressed in P. fucata somatic and ovary tissues. Ad: adductor muscle; Gi: gill; Ov: ovary; Ma: mantle; RPM: reads per million reads; Seq: high-throughput sequencing.
Biology 08 00047 g004
Biology 08 00047 g005 550
Figure 5. The predicted interactions of miR-872 and miR-463-3p with the 3′ untranslated regions (UTRs) of MPN88 family genes. The target sites of miR-872 (a) and miR-463-3p (b) in the 3′UTRs of eight MPN88 family genes identified by RNA22 and RNAhybrid, respectively.
Figure 5. The predicted interactions of miR-872 and miR-463-3p with the 3′ untranslated regions (UTRs) of MPN88 family genes. The target sites of miR-872 (a) and miR-463-3p (b) in the 3′UTRs of eight MPN88 family genes identified by RNA22 and RNAhybrid, respectively.
Biology 08 00047 g005
Biology 08 00047 g006 550
Figure 6. Hairpin structure of miRNAs and the potential target sites. (a) The potential target sites between miRNA-1990c-3p and MMP by miRanda, and between miRNA-1990c-3p and SGMP1 by RNAhybrid. (b) The potential target site between miR-9a-3p and Pif177-like protein by miRanda. (c) The potential target sites between novel-3 and GRP, between novel-3 and Linkine by RNAhybrid, and between novel-3 and OCT4 by miRanda, and between novel-3 and MMP by RNA22. (d) The potential target sites between miR-876 and Prismalin-14 by RNAhybrid and RNA22, and between miR-876 and Fam20c, between miR-876 and PMMG1 by miRanda, between miR-876 and GRP, between miR-876 and shematrin-2/2β, and between miR-876 and Tyr-1 by RNA22. The position of target genes is numbered from the 3′UTR starting position.
Figure 6. Hairpin structure of miRNAs and the potential target sites. (a) The potential target sites between miRNA-1990c-3p and MMP by miRanda, and between miRNA-1990c-3p and SGMP1 by RNAhybrid. (b) The potential target site between miR-9a-3p and Pif177-like protein by miRanda. (c) The potential target sites between novel-3 and GRP, between novel-3 and Linkine by RNAhybrid, and between novel-3 and OCT4 by miRanda, and between novel-3 and MMP by RNA22. (d) The potential target sites between miR-876 and Prismalin-14 by RNAhybrid and RNA22, and between miR-876 and Fam20c, between miR-876 and PMMG1 by miRanda, between miR-876 and GRP, between miR-876 and shematrin-2/2β, and between miR-876 and Tyr-1 by RNA22. The position of target genes is numbered from the 3′UTR starting position.
Biology 08 00047 g006
Table
Table 1. Primers used in the present study.
Table 1. Primers used in the present study.
Small RNAStem-Loop PrimerqRT-PCR Forward Primer
let-7aCTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGAACTATACACACTCCAGCTGGGTGAGGTAGTAGGTTGT
miR-1493CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGACTGATGTACACTCCAGCTGGGAGAACTGTGTATGGAC
miR-1990c-3pCTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGCAAGTAGACACTCCAGCTGGGCGGGACTACGTCAACT
miR-1993CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGTCTCGTGAACACTCCAGCTGGGTATTATGCTGTTATTC
miR-279CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGGATGAGTACACTCCAGCTGGGTGACTAGATCCACAC
miR-876CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGTCACGGATACACTCCAGCTGGGTGGATTTCCCAAGAT
miR-9a-3pCTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGCCTCCGGTACACTCCAGCTGGGATAAAGCTAGGTTAC
miR-183CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGCCGTGAATACACTCCAGCTGGGAATGGCACTGGTAGAAT
miR-200aCTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGACATCTTACACTCCAGCTGGGTAATACTGTCAGGTAAA
novel-1CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGCGGAATCACACTCCAGCTGGGAGGCGAGCCTAAACGA
novel-3CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGCAGGTAGACACTCCAGCTGGGTGCCGTCACAAGGACT
novel-10CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGATACTGGACACTCCAGCTGGGTGCACCAAACTAATGCC
miRNA reverse primerTGGTGTCGTGGAGTCG
U6 forward primerTTGCTTCGGCGGTACATATA
U6 reverse primerATTTGCGTGTCATCCTTGC
Table
Table 2. Counts of filtered sequencing reads (in millions) for the small RNAs from different tissues.
Table 2. Counts of filtered sequencing reads (in millions) for the small RNAs from different tissues.
LibrariesAd1Ad2Gi1Gi2Ov1Ov2Ma1Ma2Total
Raw reads1.965.1413.595.746.743.045.578.5450.32
Step 11.584.6111.365.385.632.73.646.8241.72
Step 21.574.6111.355.375.632.73.646.8141.68
Step 31.333.599.424.313.971.923.075.6533.26
Step 1: Remove adapter sequences and remaining reads of 15–31 nt in length. Step 2: Remove known RNAs. Step 3: Remove unmapped reads from the reference genome.

Share and Cite

MDPI and ACS Style

Huang, S.; Ichikawa, Y.; Yoshitake, K.; Kinoshita, S.; Igarashi, Y.; Omori, F.; Maeyama, K.; Nagai, K.; Watabe, S.; Asakawa, S. Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata). Biology 2019, 8, 47. https://doi.org/10.3390/biology8020047

AMA Style

Huang S, Ichikawa Y, Yoshitake K, Kinoshita S, Igarashi Y, Omori F, Maeyama K, Nagai K, Watabe S, Asakawa S. Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata). Biology. 2019; 8(2):47. https://doi.org/10.3390/biology8020047

Chicago/Turabian Style

Huang, Songqian, Yuki Ichikawa, Kazutoshi Yoshitake, Shigeharu Kinoshita, Yoji Igarashi, Fumito Omori, Kaoru Maeyama, Kiyohito Nagai, Shugo Watabe, and Shuichi Asakawa. 2019. "Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata)" Biology 8, no. 2: 47. https://doi.org/10.3390/biology8020047

APA Style

Huang, S., Ichikawa, Y., Yoshitake, K., Kinoshita, S., Igarashi, Y., Omori, F., Maeyama, K., Nagai, K., Watabe, S., & Asakawa, S. (2019). Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata). Biology, 8(2), 47. https://doi.org/10.3390/biology8020047

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop