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.     

cDNA Microarray Analysis Revealing Candidate Biomineralization Genes of the Pearl Oyster, <Emphasis Type="Italic">Pinctada fucata martensii</Emphasis>

Biomineralization is a common biological phenomenon resulting in strong tissue, such as bone, tooth, and shell. Pinctada fucata martensii is an ideal animal for the study of biomineralization. Here, microarray technique was used to identify biomineralization gene in mantle edge (ME), mantle center (MC), and both ME and MC (ME-MC) for this pearl oyster. Results revealed that 804, 306, and 1127 contigs expressed at least three times higher in ME, MC, and ME-MC as those in other tissues. Blast against non-redundant database showed that 130 contigs (16.17 %), 53 contigs (17.32 %), and 248 contigs (22.01 %) hit reference genes (E?≤??10), among which 91 contigs, 48 contigs, and 168 contigs could be assigned to 32, 26, and 63 biomineralization genes in tissue of ME, MC, and ME-MC at a threshold of 3 times upregulated expression level. The ratios of biomineralization contigs to homologous contigs were similar at 3 times, 10 times, and 100 times of upregulated expression level in either ME, MC, or ME-MC. Moreover, the ratio of biomineralization contigs was highest in MC. Although mRNA distribution characters were similar to those in other studies for eight biomineralization genes of PFMG3, Pif, nacrein, MSI7, mantle gene 6, Pfty1, prismin, and the shematrin, most biomineralization genes presented different expression profiles from existing reports. These results provided massive fundamental information for further study of biomineralization gene function, and it may be helpful for revealing gene nets of biomineralization and the molecular mechanisms underlining formation of shell and pearl for the oyster.     

Tissue inhibitor of metalloproteinase gene from pearl oyster Pinctada martensii participates in nacre formation

Tissue inhibitors of metalloproteinases (TIMPs) are nature inhibitors of matrix metalloproteinases and play a vital role in the regulation of extracellular matrix turnover, tissue remodeling and bone formation. In this study, the molecular characterization of TIMP and its potential function in nacre formation was described in pearl oyster Pinctada martensii. The cDNA of TIMP gene in P. martensii (Pm-TIMP) was 901 bp long, containing a 5′ untranslated region (UTR) of 51 bp, a 3′ UTR of 169 bp, and an open reading fragment (ORF) of 681 bp encoding 226 amino acids with an estimated molecular mass of 23.37 kDa and a theoretical isoelectric point of 5.42; The predicted amino acid sequence had a signal peptide, 13 cysteine residues, a N-terminal domain and a C-terminal domain, similar to that from other species. Amino acid multiple alignment showed Pm-TIMP had the highest (41%) identity to that from Crassostrea gigas. Tissue expression analysis indicated Pm-TIMP was highly expressed in nacre formation related-tissues, including mantle and pearl sac. After decreasing Pm-TIMP gene expression by RNA interference (RNAi) technology in the mantle pallium, the inner nacreous layer of the shells showed a disordered growth. These results indicated that the obtained Pm-TIMP in this study participated in nacre formation.     

Identification of microRNAs in Macaca fascicularis (Cynomolgus Monkey) by Homology Search and Experimental Validation by Small RNA-Seq and RT-qPCR Using Kidney Cortex Tissues

MicroRNAs (miRNAs) present in tissues and biofluids are emerging as sensitive and specific safety biomarkers. MiRNAs have not been thoroughly described in M. fascicularis, an animal model used in pharmaceutical industry especially in drug safety evaluation. Here we investigated the miRNAs in M. fascicularis. For Macaca mulatta, a closely related species of M. fascicularis, 619 stem-loop precursor miRNAs (pre-miRNAs) and 914 mature miRNAs are available in miRBase version 21. Using M. mulatta miRNAs as a reference list and homology search tools, we identified 604 pre-miRNAs and 913 mature miRNAs in the genome of M. fascicularis. In order to validate the miRNAs identified by homology search we attempted to sequence miRNAs expressed in kidney cortex from M. fascicularis. MiRNAs expressed in kidney cortex may indeed be released in urine upon kidney cortex damage and be potentially used to monitor drug induced kidney injury. Hence small RNA sequencing libraries were prepared using kidney cortex tissues obtained from three naive M. fascicularis and sequenced. Analysis of sequencing data indicated that 432 out of 913 mature miRNAs were expressed in kidney cortex tissues. Assigning these 432 miRNAs to pre-miRNAs revealed that 273 were expressed from both the -5p and -3p arms of 150 pre-miRNAs and 159 miRNAs expressed from either the -5p or -3p arm of 176 pre-miRNAs. Mapping sequencing reads to pre-miRNAs also facilitated the detection of twenty-two new miRNAs. To substantiate miRNAs identified by small RNA sequencing, 313 miRNAs were examined by RT-qPCR. Expression of 262 miRNAs in kidney cortex tissues ware confirmed by TaqMan microRNA RT-qPCR assays. Analysis of kidney cortex miRNA targeted genes suggested that they play important role in kidney development and function. Data presented in this study may serve as a valuable resource to assess the renal safety biomarker potential of miRNAs in Cynomolgus monkeys.     

Genome-wide identification and characterization of perirenal adipose tissue microRNAs in rabbits fed a high-fat diet

MicroRNAs (miRNAs) are a class of endogenous single-stranded RNA molecules that play an important role in gene regulation in animals by pairing with target gene mRNA. Extensive evidence shows that miRNAs are key players in metabolic regulation and the development of obesity. However, the systemic understanding of miRNAs in the adipogenesis of obese rabbits need further investigation. Here, seven small RNA libraries from rabbits fed either a standard normal diet (SND; n=3) or high-fat diet (HFD; n=4) were constructed and sequenced. Differentially expressed (DE) miRNAs were identified using the edgeR data analysis package from R. Software miRanda and RNAhybrid were used to predict the target genes of miRNAs. To further explore the functions of DE miRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. A total of 81449996 clean reads were obtained from the seven libraries, of which, 52 known DE miRNAs (24 up-regulated, 28 down-regulated) and 31 novel DE miRNAs (14 up-regulated, 17 down-regulated) were identified. GO enrichment analysis revealed that the DE miRNAs target genes were involved in intermediate filament cytoskeleton organization, intermediate filament-based process, and α-tubulin binding. DE miRNAs were involved in p53 signaling, linoleic acid metabolism, and other adipogenesis-related KEGG pathways. Our study further elucidates the possible functions of DE miRNAs in rabbit adipogenesis, contributing to the understanding of rabbit obesity.     

Highly efficient identification of thousands of microsatellite loci in the pearl oyster (Pinctada martensii) from RNA-Seq

High-throughput RNA-Seq affords a cost and time effective means of obtaining large numbers of genetic markers for aquatic genomics. Here, we present thousands of novel microsatellite loci developed for the pearl oyster, Pinctada martensii from the Illumina HiSeq™ 2000 library of the pearl sac. Free user-friendly bioinformatics tools were employed to screen for microsatellite loci and design appropriate primers in 102,762 unigenes with 7216 microsatellite loci identified in total, 4862 of which had flanking sequences suitable for polymerase chain reaction primer design. The 50 randomly chosen primer pairs were tested in two populations of pearl oyster (base population (POP1) and selected population (POP2), with 30 individuals of each population). All the primer pairs were amplified successfully in two populations. All loci were polymorphic in POP1, while there were 3 loci showing monomorphism in POP2. In POP1 and POP2, observed heterozygosity from 0.033 to 1.000 and 0.000 to 1.000, 19 and 16 microsatellite loci deviated significantly from Hardy–Weinberg expectations including a Bonferroni correction (P < 0.001). Thirteen loci were highly informative content (PIC ≥ 0.5) in both populations. These identified loci will be useful for potential application for evolutionary, population genetic and chromosome linkage mapping research on pearl oyster.