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

Integrated mRNA and miRNA transcriptomic analysis reveals the response of Rapana venosa to the metamorphic inducer (juvenile oysters)

Metamorphosis, as a critical developmental event, controls the population dynamics of most marine invertebrates, especially some carnivorous gastropods that feed on bivalves, whose population dynamics not only affect the maintenance of the ecological balance but also impact the protection of bivalve resources; therefore, the metamorphosis of carnivorous gastropods deserve attention. Here, we investigated the mechanism underlying the response of the carnivorous gastropod Rapana venosa to its metamorphic inducer juvenile oysters through integrated analysis of miRNA and mRNA profiles. According to the results, we speculated that the AMPK signaling pathway may be the critical regulator in the response to juvenile oysters in R. venosa competent larvae. The NF-kB and JAK-STAT signaling pathways that regulated apoptosis were also activated by the metamorphic inducer, which may result in the degeneration of the velum. Additionally, the significant changes in the expression of the SARP-19 precursor gene and protein cibby homolog 1-like gene may indicate that these signaling pathways also regulate growth and development during metamorphosis. This study provides further evidence that juvenile oysters can induce metamorphosis of R. venosa at the transcriptional level, which expands our understanding of the metamorphosis mechanism in carnivorous gastropods.

Identification and Characterization of miRNAs and Their Predicted mRNAs in the Larval Development of Pearl Oyster Pinctada fucata

As an important economic shellfish, the pearl oyster, Pinctada fucata, and its larvae are an ideal model for studying molecular mechanisms of larval development in invertebrates. Larval development directly affects the quantity and quality of pearl oysters. MicroRNAs (miRNAs) may play important roles in development, but the effects of miRNA expression on P. fucata early development remain unknown. In this study, miRNA and mRNA transcriptomics of seven different P. fucata developmental stages were analyzed using Illumina RNA sequencing. A total of 329 miRNAs, including 87 known miRNAs and 242 novel miRNAs, and 33,550 unigenes, including 26,333 known genes and 7217 predicted new genes, were identified in these stages. A cluster analysis showed that the difference in the numbers of miRNAs was greatest between fertilized eggs and trochophores. In addition, the integrated mRNA transcriptome was used to predict target genes for differentially expressed miRNAs between adjacent developmental stages, and the target genes were subjected to a gene ontology enrichment analysis. Using the gene ontology annotation, 100 different expressed genes and 95 differentially expressed miRNAs were identified as part of larval development regulation. Real-time PCR was used to identify eight mRNAs and three miRNAs related to larval development. The present findings will be helpful for further clarifying the regulatory mechanisms of miRNA in invertebrate larval development.

Differential Expression of Long Non-Coding RNA (lncRNA) in Mediterranean Mussel (Mytilus galloprovincialis) Hemocytes under Immune Stimuli

The Mediterranean mussel is one of the most economically relevant bivalve mollusk species in Europe and China. The absence of massive mortalities and their resistance to pathogens affecting other cultured bivalves has been under study in recent years. The transcriptome response of this species to different immune stimuli has been extensively studied, and even the complexity of its genome, which has recently been sequenced, has been suggested as one of the factors contributing to this resistance. However, studies concerning the non-coding RNA profiles remain practically unexplored-especially those corresponding to the lncRNAs. To the best of our knowledge, this is the second characterization and study of lncRNAs in this bivalve species. In this work, we identified the potential repertoire of lncRNAs expressed in mussel hemocytes, and using RNA-Seq we analyzed the lncRNA profile of mussel hemocytes stimulated in vitro with three different immune stimuli: LPS, poly I:C, and β-glucans. Compared to unstimulated hemocytes, LPS induced the highest modulation of lncRNAs, whereas poly I:C and β-glucans induced a similar discrete response. Based on the potential cis-regulatory activity of the lncRNAs, we identified the neighboring protein-coding genes of the regulated lncRNAs to estimate-at least partially-the processes in which they are implicated. After applying correlation analyses, it seems that-especially for LPS-the lncRNAs could participate in the regulation of gene expression, and substantially contribute to the immune response.

Digging into bivalve miRNAomes: between conservation and innovation

Bivalves are a diverse mollusc group of economic and ecological importance. An evident resilience to pollution, parasites and extreme environments makes some bivalve species important models for studying adaptation and immunity. Despite substantial progress in sequencing projects of bivalves, information on non-coding genes and gene-regulatory aspects is still lacking. Here, we review the current repertoire of bivalve microRNAs (miRNAs), important regulators of gene expression in Metazoa. We exploited available short non-coding RNA (sncRNA) data for Pinctada martensii, Crassostrea gigas, Corbicula fluminea, Tegillarca granosa and Ruditapes philippinarum, and we produced new sncRNA data for two additional bivalves, the Mediterranean mussel Mytilus galloprovincialis and the blood clam Scapharca broughtonii. We found substantial heterogeneity and incorrect annotations of miRNAs; hence, we reannotated conserved miRNA families using recently established criteria for bona fide microRNA annotation. We found 106 miRNA families missing in the previously published bivalve datasets and 89 and 87 miRNA complements were identified in the two additional species. The overall results provide a homogeneous and evolutionarily consistent picture of miRNAs in bivalves and enable future comparative studies. The identification of two bivalve-specific miRNA families sheds further light on the complexity of transcription and its regulation in bivalve molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Identification of shell-color-related microRNAs in the Manila clam Ruditapes philippinarum using high-throughput sequencing of small RNA transcriptomes

Shell-color polymorphism is a common phenomenon in several mollusk species and has been associated with thermal capacity, developmental stability, shell strength, and immunity. Shell-color polymorphism has been related to the differential expression of genes in several signal transduction pathways; however, the functions of micro-RNAs (miRNAs) in shell-color formation remain unclear. In the present study, we compared high-quality, small-RNA transcriptomes in three strains of the Manila clam Ruditapes philippinarum with specific shell-color patterns, artificially selected for six generations. Totals of 114 known and 208 novel miRNAs were identified by high-throughput sequencing, of which nine known and one novel miRNA were verified by stem-loop quantitative real time-polymerase chain reaction. Predicted miRNA targets were subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. miR-137 and miR-216b and the Hedgehog signaling pathway and Wnt signaling pathway were identified as being potentially involved in pigment formation and regulation in R. philippinarum. These results may help to clarify the role of miRNAs in shell coloration and shed light on the mechanisms regulating color formation in bivalve shells.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

Comprehensive analysis of microRNAs in the mantle central and mantle edge provide insights into shell formation in pearl oyster Pinctada fucata martensii

MicroRNAs (miRNAs) are a class of non-coding RNA molecules with post-transcriptional regulatory activity in various biological processes. Pearl oyster Pinctada fucata martensii is one of the main species cultured for marine pearl production in China and Japan. In this study, we constructed two small RNA libraries of mantle central (MC) and mantle edge (ME) from P. f. martensii and obtained 24,175,537 and 21,593,898 clean reads, respectively. A total of 258 miRNAs of P. f. martensii (Pm-miRNA) were identified, and 93 differentially expressed miRNAs (DEMs) including 49 known Pm-miRNAs and 44 novel Pm-miRNAs were obtained from the MC and ME. The target transcripts of these DEMs were obviously enriched in neuroactive ligand-receptor interaction pathway, and others. After over-expression of Pm-miR-124 and Pm-miR-9a-5p in the MC by mimic injection into the muscle of P. f. martensii, nacre exhibited a disorderly growth as detected by scanning electron microscopy. Pm-nicotinic acetylcholine receptor alpha subunit, Pm-neuropeptide Y and Pm-chitin synthase were investigated as the targets of Pm-miR-124; and Pm-tumor necrosis factor receptor associated factor 2 and Pm-chitin synthase were investigated as the targets of Pm-miR-9a-5p. These predicted target transcripts were down-regulated after the over-expression of Pm-miR-124 and Pm-miR-9a-5p in MC. This study comprehensively analyzed the miRNAs in mantle tissues to enhance our understanding of the regulatory mechanism underlying shell formation.

Profiling of microRNAs and mRNAs in marine mussel Mytilus galloprovincialis

MicroRNAs (miRNAs) are a class of noncoding RNA molecules containing 18-24 nucleotides, and those with conserved structures are able to regulate the expression of eukaryotic genes by inhibition or enhancement of mRNA translation. However, miRNAs of the blue mussel, Mytilus galloprovincialis have not been reported. M. galloprovincialis is a primary species distributed along coastal zones worldwide. To reveal the repertoire of miRNAs in M. galloprovincialis, we constructed small RNA libraries prepared from three different mussels, which were then sequenced by Solexa deep sequencing technology. A total of 32,836,817, 33,359,113 and 33,093,562 clean reads from the tissues of the three M. galloprovincialis were obtained. Based on sequence similarities and hairpin structure predictions, 137 M. galloprovincialis miRNAs (mg-miRNA) were identified. Among the mg-miRNAs, 104 were conserved across species, whereas 33 might be novel and specific for M. galloprovincialis. Some of the mg-miRNAs, such as let-7 and the miR-100 family are playing key roles in many metabolic pathways and are worthy of further study. By performing a whole genome-scale characterization of mg-miRNAs and proposing their potential functions, these results provide a foundation for understanding the biological processes of the blue mussel, M. galloprovincialis.

PmCBP, a novel poly (chitin-binding domain) gene, participates in nacreous layer formation of Pinctada fucata martensii

Chitin participates in shell formation as the main component of an organic framework. Chitin-binding protein contains domains that can bind to chitin specifically. In this study, a novel chitin-binding protein from Pinctada fucata martensii (PmCBP) with poly (chitin-binding domain) was cloned, which contains a 5'-untranslated region (UTR) of 114 bp and 3'UTR of 116 bp, and encodes a putative protein of 2044 amino acids. The predicted PmCBP protein was structurally typical of the CBP family with 20 ChtBD2 domains. Phylogenetic and linear relation analyses showed that the ChtBD2 domain has a highly conserved structure among the three species of P. f. martensii, Crassostrea gigas, and Mizuhopecten yessoensis. qRT-PCR and in-situ hybridization analysis revealed that PmCBP was most abundant in the mantle pallium whose expression level was significantly correlated with the growth traits. After RNAi, PmCBP expression was significantly inhibited in the mantle pallium (P < 0.05) and the microstructure of nacreous layers showed a disordered growth in the experiment group. These results indicated that PmCBP may be involved in nacreous layer formation through participation in the process of binding chitin in pearl oyster P. f. martensii.

Improved annotation of Lutzomyia longipalpis genome using bioinformatics analysis

Lutzomyia longipalpis, a sand fly, is a vector-spreading pathogenic protozoan in the New World. MicroRNA (miRNA) is evolutionarily-conserved non-coding RNA, which plays critical roles in various biological processes. To date, the functions of most proteins in L. longipalpis are unknown, and few studies have addressed the roles of miRNAs in this species. In the present study, we re-annotated the protein-coding genes and identified several miRNAs using a set of comparative genomics tools. A large number of L. longipalpis proteins were found to be homologous with those in the mosquito genome, indicating that they may have experienced similar selective pressures. Among these proteins, a set of 19 putative salivary proteins were identified, which could be used for studying the transmission of Leishmania. Twenty-one novel miRNAs were characterized, including two miRNAs, miR-4113-5p and miR-5101, which are unique to L. longipalpis. Many of the targets of these two genes were found to be involved in ATP hydrolysis-coupled proton transport, suggesting that they may have important roles in the physiology of energy production. Topology analysis of the miRNA-gene network indicated that miR-9388-5p and miR-3871-5p regulate several critical genes in response to disease development. In conclusion, our work provides a basis for improving the genome annotation of L. longipalpis, and opens a new door to understanding the molecular regulatory mechanisms in this species.

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

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.

Circular RNA expression profile in gingival tissues identifies circ_0062491 and circ_0095812 as potential treatment targets

Circular RNAs (circRNAs) emerging as a novel class of endogenous, conserved noncoding RNAs, which have been reported to be participated in immune-inflammatory response recently, yet their function in periodontal inflammation has remained elusive. This study aimed to analyze the specific circRNAs expressed in periodontal inflammation process through Illumina (San Diego, CA) sequencing technology combining with experimental validation. The inflamed and healthy gingival tissues from patients were selected to explore the expression statues of circRNAs. In brief, 1304 dysregulated circRNAs were identified in inflamed gingival tissues. Besides, Gene Ontology (GO) enrichment analysis was conducted to investigate the functions of abnormally expressed circRNAs in periodontitis as well as their host-linear genes. Furthermore, the interaction network of circRNAs-miRNAs (microRNA) was constructed to reveal the key circRNAs (circ_0095812, circ_0120299, circ_0125699, circ_0062491, and circ_0043115) in the pathobiology of periodontitis. Subsequently, we have investigated the expression pattern of circ_0062491 (downregulated) and circ_0095812 (upregulated) among 30-paired periodontitis patients and healthy individuals by quantitative real-time polymerase chain reaction. Moreover, circ_0062491 was identified as the sponge of miR-584 which play a key role in periodontitis.

Transcriptome Dynamics During Turbot Spermatogenesis Predicting the Potential Key Genes Regulating Male Germ Cell Proliferation and Maturation

Spermatogenesis is a dynamic developmental process in which spermatogonial stem cells proliferate, differentiate and mature into functional spermatozoa. These processes require an accurate gene regulation network. Here, we investigated the dynamic changes that occur during spermatogenesis through a combination of histological and transcriptome analyses of different developmental stages of the testis. We constructed 18 testis transcriptome libraries, and the average length, N50, and GC content of the unigenes were 1,795 bp; 3,240 bp and 49.25%, respectively. Differentially expressed genes (DEGs) that were related to germ cell proliferation and maturation, such as NANOS3, RARs, KIFs, steroid hormone synthesis-related genes and receptor genes, were identified between pairs of testis at different developmental stages. Gene ontology annotation and pathway analyses were conducted on DEGs with specific expression patterns involved in the regulation of spermatogenesis. Nine important pathways such as steroid hormone biosynthesis related to spermatogenesis were identified. A total of 21 modules that ranged from 49 to 7,448 genes were designed by a weighted gene co-expression network analysis. Furthermore, a total of 83 candidate miRNA were identified by computational methods. Our study provides the first transcriptomic evidence for differences in gene expression between different developmental stages of spermatogenesis in turbot (Scophthalmus maximus).

Identification and characterization of microRNAs involved in scale biomineralization in the naked carp Gymnocypris przewalskii

The mineralized scale derived from skin plays a protective role for the fish body and also possesses important application values in the biomaterial field. However, little is known about fish scale biomineralization and related molecular regulatory mechanisms. Here, we used a comparative microRNA sequencing approach to identify and characterize differentially expressed microRNAs (DEMs) involved in scale biomineralization in the naked carp Gymnocypris przewalskii. A total of 18, 43, and 66 DEMs were obtained from skin tissues covered with initial, developing, and mature scales (IS, DS, and MS) compared with scale-uncovered skin. The target genes of these DEMs were significantly enriched in a sole biomineralization-related sphingolipid signaling pathway. Seven DEMs (dre-miR-124-3p, dre-miR-133a-2-5p, dre-miR-184, dre-miR-206-3p, novel_33, novel_56 and novel_75) were common in IS, DS, and MS. Dre-miR-124-3p, dre-miR-206-3p, and novel_33 were predicted to be able to target biomineralization-related genes. Stem-loop real-time quantitative PCR further confirmed that the common DEMs had higher expression levels in scale-covered skin tissues than that in the gill, intestine, and brain, except for dre-miR-133a-2-5p. Our results suggest that these identified microRNAs may play a role in scale biomineralization in G. przewalskii, and the obtained microRNAs are expected to be candidates in understanding the molecular mechanism of scale biomineralization in fish species.

Molecular modularity and asymmetry of the molluscan mantle revealed by a gene expression atlas

Background

Conchiferan molluscs construct a biocalcified shell that likely supported much of their evolutionary success. However, beyond broad proteomic and transcriptomic surveys of molluscan shells and the shell-forming mantle tissue, little is known of the spatial and ontogenetic regulation of shell fabrication. In addition, most efforts have been focused on species that deposit nacre, which is at odds with the majority of conchiferan species that fabricate shells using a crossed-lamellar microstructure, sensu lato.

Results

By combining proteomic and transcriptomic sequencing with in situ hybridization we have identified a suite of gene products associated with the production of the crossed-lamellar shell in Lymnaea stagnalis. With this spatial expression data we are able to generate novel hypotheses of how the adult mantle tissue coordinates the deposition of the calcified shell. These hypotheses include functional roles for unusual and otherwise difficult-to-study proteins such as those containing repetitive low-complexity domains. The spatial expression readouts of shell-forming genes also reveal cryptic patterns of asymmetry and modularity in the shell-forming cells of larvae and adult mantle tissue.

Conclusions

This molecular modularity of the shell-forming mantle tissue hints at intimate associations between structure, function, and evolvability and may provide an elegant explanation for the evolutionary success of the second largest phylum among the Metazoa.

Understanding microRNA Regulation Involved in the Metamorphosis of the Veined Rapa Whelk (Rapana venosa)

The veined rapa whelk (Rapana venosa) is widely consumed in China. Nevertheless, it preys on oceanic bivalves, thereby reducing this resource worldwide. Its larval metamorphosis comprises a transition from pelagic to benthic form, which involves considerable physiological and structural changes and has vital roles in its natural populations and commercial breeding. Thus, understanding the endogenous microRNAs (miRNAs) that drive metamorphosis is of great interest. This is the first study to use high-throughput sequencing to examine the alterations in miRNA expression that occur during metamorphosis in a marine gastropod. A total of 195 differentially expressed miRNAs were obtained. Sixty-five of these were expressed during the transition from precompetent to competent larvae. Thirty-three of these were upregulated and the others were downregulated. Another 123 miRNAs were expressed during the transition from competent to postlarvae. Ninety-six of these were upregulated and the remaining 27 were downregulated. The expression of miR-276-y, miR-100-x, miR-183-x, and miR-263-x showed a >100-fold change during development, while the miR-242-x and novel-m0052-3p expression levels changed over 3000-fold. Putative target gene coexpression, gene ontology, and pathway analyses suggest that these miRNAs play important parts in cell proliferation, migration, apoptosis, metabolic regulation, and energy absorption. Twenty miRNAs and their target genes involved in ingestion, digestion, cytoskeleton, cell adhesion, and apoptosis were identified. Nine of them were analyzed with real-time polymerase chain reaction (PCR), which showed an inverse correlation between the miRNAs and their relative expression levels. Our data elucidate the role of miRNAs in R. venosa metamorphic transition and serve as a solid basis for further investigations into regulatory mechanisms of gastropod metamorphosis.

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.