Identification and characterization of microRNAs in pearl oyster Pinctada martensii by Solexa deep sequencing

The expression patterns of exosomal miRNAs in the Pacific oyster after high-temperature stress or Vibrio stimulation

The exosomal miRNA plays a crucial role in the intercellular communication response to environmental stress and pathogenic stimulation. In the present study, the expression of exosomal miRNAs in the Pacific oyster Crassostrea gigas after high-temperature stress or Vibrio splendidus stimulation was investigated through high-throughput sequencing. The exosomes were identified to be teardrop-like vesicles with the average size of 81.7 nm by transmission electron microscopy. There were 66 known miRNAs and 33 novel miRNAs identified, of which 10 miRNAs were differentially expressed after both high-temperature stress and Vibrio stimulation compared to the control group. A total of 1868 genes were predicted as the putative targets of miRNAs, of which threonine aspartase 1-like was targeted by the highest number of related miRNAs. The robustness and reliability of miRNA expression from the sRNA sequencing data were verified by employing eight miRNAs for qPCR. GO and KEGG clustering analyses revealed that apoptosis was significantly enriched by the target genes of differentially expressed exosomal miRNAs after high-temperature stress, and autophagy and cytokine activity were significantly enriched after Vibrio stimulation. Energy metabolism was found to be significantly shared in the target gene enrichments after both high-temperature stress and Vibrio stimulation. These findings would improve our understanding of the regulatory mechanisms of exosomal miRNAs in C. gigas after high-temperature stress or Vibrio stimulation.

The role of Pm-miR-184-3p in regulating the immune response in the pearl oyster Pinctada fucata martensii

microRNAs are a class of non-coding RNAs with post-transcriptional regulatory functions in eukaryotes. In our previous study, miR-184-3p was identified in the hemocyte transcriptome of Pinctada fucata martensii (Pm-miR-184-3p), and its expression was shown to be up-regulated following transplantation surgery; however, its role in regulating transplantation immunity has not yet been clarified. Here, the role of Pm-miR-184-3p in regulating the immune response of P. f. martensii was studied. The expression of Pm-miR-184-3p increased following the stimulation of pathogen-associated molecular patterns, and Pm-miR-184-3p overexpression increased the activity of antioxidant-related enzymes, such as superoxide dismutase and catalase. Transcriptome analysis obtained 1096 differentially expressed genes (DEGs) after overexpression of Pm-miR-184-3p, and these DEGs were significantly enriched in conserved pathways such as the Cell cycle pathway and NF-kappa B signaling pathway, as well as GO terms including base excision repair, cell cycle, and DNA replication, suggesting that Pm-miR-184-3p could enhance the inflammation process. Target prediction and dual luciferase analysis revealed that pro-inflammatory related genes Pm-TLR3 and Pm-FN were the potential target of Pm-miR-184-3p. We speculate that Pm-miR-184-3p may utilize negative regulation of target genes to delay the activation of corresponding immune pathways, potentially preventing excessive inflammatory responses and achieving a delicate balance within the organism. Overall, Pm-miR-184-3p play a key role in regulating cellular responses to transplantation. Our findings provide new insights into the immune response of P. f. martensii to transplantation.

miRNA-seq provides novel insight into the response to hyper- and hypo- salinity acclimation in Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis.

Mitochondrially mediated RNA interference, a retrograde signaling system affecting nuclear gene expression

Several functional classes of short noncoding RNAs are involved in manifold regulatory processes in eukaryotes, including, among the best characterized, miRNAs. One of the most intriguing regulatory networks in the eukaryotic cell is the mito-nuclear crosstalk: recently, miRNA-like elements of mitochondrial origin, called smithRNAs, were detected in a bivalve species, Ruditapes philippinarum. These RNA molecules originate in the organelle but were shown in vivo to regulate nuclear genes. Since miRNA genes evolve easily de novo with respect to protein-coding genes, in the present work we estimate the probability with which a newly arisen smithRNA finds a suitable target in the nuclear transcriptome. Simulations with transcriptomes of 12 bivalve species suggest that this probability is high and not species specific: one in a hundred million (1 × 10-8) if five mismatches between the smithRNA and the 3' mRNA are allowed, yet many more are allowed in animals. We propose that novel smithRNAs may easily evolve as exaptation of the pre-existing mitochondrial RNAs. In turn, the ability of evolving novel smithRNAs may have played a pivotal role in mito-nuclear interactions during animal evolution, including the intriguing possibility of acting as speciation trigger.

Immunomodulatory effects of one novel microRNA miR-63 in pearl oyster Pinctada fucata martensii

Novel microRNA miR-63 (novel-miR-63) from pearl oyster Pinctada fucata martensii (Pm-novel-miR-63) is a species-specific miRNA. Our previous research has shown that the expression of Pm-novel-miR-63 was significantly downregulated at 24 h after nucleus transplantation. In this study, we analyzed the function and regulatory role of Pm-novel-miR-63 in the immune response of pearl oysters. The results showed that Pm-novel-miR-63 expression increased after the stimulation of pathogen associated molecular patterns at 6-12 h, and the activity of immune and antioxidant enzymes in the serum decreased after Pm-novel-miR-63 overexpression. Transcriptome analysis revealed that Pm-novel-miR-63 participated in regulating transplantation immunity through the Notch and mRNA surveillance signaling pathways. Target prediction and dual luciferase analysis revealed that Pm-GDP-FucTP, Pm-CysLTR2, and Pm-RLR were the target genes of Pm-novel-miR-63. These results suggested that Pm-novel-miR-63 participated in regulating the immune response in pearl oysters and can serve as a new interference target to reasonably control excessive immune rejection in pearl culture.

miR-10a-3p Participates in Nacre Formation in the Pearl Oyster Pinctada fucata martensii by Targeting NPY

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression via the recognition of their target messenger RNAs. MiR-10a-3p plays an important role in the process of ossification. In this study, we obtained the precursor sequence of miR-10a-3p in the pearl oyster Pinctada fucata martensii (Pm-miR-10a-3p) and verified its sequence by miR-RACE technology, and detected its expression level in the mantle tissues of the pearl oyster P. f. martensii. Pm-nAChRsα and Pm-NPY were identified as the potential target genes of Pm-miR-10a-3p. After the over-expression of Pm-miR-10a-3p, the target genes Pm-nAChRsα and Pm-NPY were downregulated, and the nacre microstructure became disordered. The Pm-miR-10a-3p mimic obviously inhibited the luciferase activity of the 3' untranslated region of the Pm-NPY gene. When the interaction site was mutated, the inhibitory effect disappeared. Our results suggested that Pm-miR-10a-3p participates in nacre formation in P. f. martensii by targeting Pm-NPY. This study can expand our understanding of the mechanism of biomineralization in pearl oysters.

Identification and comparative analysis of miRNA transcriptomes after allograft and xenograft transplantation in Pinctada fucata martensii

Effective immune regulation after transplantation during pearl production is crucial for the cultivation of high-quality pearls. MicroRNAs (miRNAs) play an important role in a variety of physiological processes. To understand the regulatory rules of miRNAs after transplantation in Pinctada funcata martensii, we constructed 13 miRNA transcriptomes, including the control group (Con), allograft (Al), and xenograft (Xe) transplantation at six time points (6, 12, and 24 h and 3, 6, and 12 days), in which the xenografted mantle tissue was from Pinctada maxima. We identified 159 differentially expressed miRNAs (DEMs) and found that these DEMs showed high expression at 12 h, 24 h, and 3 days after transplantation. A total of 130 DEMs, such as Let-7, were present in the Al and Xe groups; miR-34 and 16 other DEMs were specifically present in the Al group; miR-216b and 13 other DEMs were specifically present in the Xe group. Compared with the Con group, the target genes of DEMs in the Al group were significantly enriched in protein complex, cytoskeleton, and macromolecular complex, and the Xe group was significantly enriched in ribonucleoside metabolic process, nucleoside binding, and cell division. Compared with the Al group, the target genes in the Xe group were significantly enriched in response to DNA damage stimulation. Overall, multiple pathways associated with cellular activity were enriched in higher numbers of genes in the Xe group than in the Al group. These findings enriched the information on immune regulatory mechanisms at the expression level of miRNAs in P. f. martensii after transplantation.

Histone deacetylase inhibitor butyrate inhibits the cellular immunity and increases the serum immunity of pearl oyster Pinctada fucata martensii

Histone acetylation is a dynamic epigenetic modification and sensitive to the changes in extracellular environment. Butyrate, a histone deacetylase inhibitor, can inhibit the deacetylation process of histones. In this study, we found that the acetylation level of H3 was enhanced at 12 h after lipopolysaccharide (LPS) stimulation and increased at 6 h after combining treatment with LPS and butyrate in pearl oyster Pinctada fucata martensii. Transcriptome analysis indicated that butyrate counter-regulated 29.95%-36.35% of the genes repressed by LPS, and these genes were mainly enriched in the "cell proliferation" and "Notch signaling pathway". Meanwhile, butyrate inhibited the up-regulation of 31.54%-54.96% of the genes induced by LPS, and these genes were mainly enriched in "Notch signaling pathway", "cell proliferation", "NF-kappa B signaling pathway", "TNF signaling pathway", "apoptosis", "NOD-like receptor signaling pathway", "RIG-I-like receptor signaling pathway" and "cytosolic DNA-sensing pathway". Gene expression analysis showed that butyrate downregulated most of cell proliferation, immune-related genes effected by LPS. The activities of LAP, LYS, ACP, ALP, and GSH-Px were up-regulated at 6 h after combining treatment with LPS and butyrate, suggesting that butyrate could activate serum immune-related enzymes in pearl oyster. These results can improve our understanding of the function of histone deacetylase in the immune response of pearl oyster and provide references for an in-depth study of the functions of histone deacetylase in mollusks.

Transcriptome analysis reveals the diverse response of pearl oyster Pinctada fucata martensii after different PAMP stimulation

Bivalves have evolved effective strategies to combat different pathogens in the environment. They rely on innate immunity to deal with the invasion of various bacteria, viruses, and other microorganisms. However, the molecular mechanisms underlying the responses remain largely unknown. Herein, we constructed 21 transcriptomes of the hemocytes after lipopolysaccharide (LPS), peptidoglycan (PGN) and polyinosinic-polycytidylic acid (poly(I:C)) stimulation to investigate the molecular mechanisms underlying adaptations and plastic responses to different pathogen-related molecular patterns (PAMPs) in pearl oyster Pinctada fucata martensii. Transcriptome analysis revealed 1986-3427 responsive genes enriched in the major immune and cell cycle-related pathways at different times after PAMP stimulation, and the expression patterns of genes under these pathways are complex and diverse. Moreover, "lysosomes" were enriched 6 h after LPS and PGN stimulation, while "peroxisomes" were only enriched in poly(I:C) group. These results suggest different response strategies of pearl oyster to different PAMPs. Furthermore, we identified 261 pattern-recognition receptors (PRRs) including 4 retinoic acid-inducible gene I-like receptors, 38 NOD-like receptors, 83 Toll-like receptors, and 136 C-type lectins in the genome of P. f. martensii. The diverse expression patterns of these PRRs after different PAMP stimulation indicated that pearl oyster evolved complex and specific recognition systems due to tandem repeat and diverse domain combination, which may help pearl oyster cope with the different pathogens in the environment. The present study improved our understanding of the molecular response of pearl oyster to different PAMP stimulation.

Comprehensive analysis of the lncRNAs, mRNAs, and miRNAs implicated in the immune response of Pinctada fucata martensii to Vibrio parahaemolyticus

Non-coding RNAs (ncRNAs) have been implicated in a variety of biological processes. However, most ncRNAs are of unknown function and are as-yet unannotated. The immune-related functions of ncRNAs in the pearl oyster Pinctada fucata martensii were explored based on transcriptomic differences in the expression levels of long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in the hemocytes of P.f. martensii after challenge by the pathogenic bacterium Vibrio parahaemolyticus. Across the challenged and control pearl oysters, 144 miRNAs and 14,571 lncRNAs were identified. In total, 13,375 ncRNAs were differentially expressed between the challenged and control pearl oysters; in the challenged pearl oysters as compared to the controls, 15 miRNAs and 5147 lncRNAs were upregulated, while 51 miRNAs and 8162 lncRNAs were downregulated. The sequencing results were validated using quantitative real-time polymerase chain reaction (qRT-PCR) analysis. GO and KEGG pathway analysis showed that genes targeted by the differentially expressed ncRNAs were associated with the vascular endothelial growth factor (VEGF) signaling pathway and the nuclear factor kappa-B (NF-κB) signaling pathway. An lncRNA-mRNA-miRNA network that was developed based on the transcriptomic results of this study suggested that lncRNAs may compete with miRNAs for mRNA binding sites. This study may provide a useful framework for the detection of additional novel ncRNAs, as well as new insights into the pathogenic mechanisms underlying the response of P.f. martensii to V. parahaemolyticus.

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.

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.

MiRNA-mRNA Integration Analysis Reveals the Regulatory Roles of MiRNAs in Shell Pigmentation of the Manila clam (Ruditapes philippinarum)

The shell color of the Manila clam (Ruditapes philippinarum) is an economically important trait. We used high-throughput sequencing and transcriptome analysis to study the molecular mechanisms that underlie shell color formation and regulation in this species. We constructed small RNA libraries from mantle tissues from four shell color strains of Manila clam, subjected them to high-throughput sequencing. Notably, the results suggested that a number of pigment-associated genes including Mitf, HERC2, were negatively regulated by nvi-miR-2a, tgu-miR-133-3p, respectively. They might be involved in melanin formation via the activation of the melanogenesis pathway. And aae-miR-71-5p and dme-miR-7-5p linked to shell formation-related genes such as Calmodulin and IMSP3 were considered to participate in the calcium signaling pathway. We then used quantitative PCR to verify the candidate miRNAs and target genes in different shell color groups. Our results indicated that miR-7, miR-71, and miR-133 may regulate target mRNAs to participate in shell color pigmentation. These results provide the foundation to further characterize miRNA effects on the regulation of shell color and have significant implications for the breeding of new varieties of clams.

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.

Characterization and expression analysis of tandemly duplicated nicotinic acetylcholine receptors in pearl oysters after stimulation of pathogen-related molecular patterns

Nicotinic acetylcholine receptors (nAChRs) are a class of ligand-gated ion channels that participate in signal transduction and are reported to play an important role in the immunomodulation of vertebrates and invertebrates. Previous studies have shown that the nAChRs in mollusks have undergone large-scale expansion after tandem repeats and retrotransposition, with the most expansion observed in bivalves. This study characterized the sequence of a tandem repeat nAChR unique to several bivalve mollusks and investigated its functions in Pinctada fucata martensii. Firstly, phylogenetic analysis revealed that the tandem arrays of nAChRs existed before bivalve differentiation and m ost tandem-replicated nAChR genes have a conserved genomic structure and domain combination. In present study, five tandemly duplicated nAChR genes were cloned from P. f. martensii and designated as PmnAChR-1 to PmnAChR-5. qRT-PCR analysis revealed that five PmnAChRs were specifically expressed in adult gills. In addition, after PAMP stimulation, the expression of PmnAChRs in hemocytes of P. f. martensii were strongly induced but exhibited different responses to different stimuli. PmnAChR-1, PmnAChR-4, and PmnAChR-5 exhibited strong and wide responsiveness to lipopolysaccharide (LPS) stimulation but had no response to peptidoglycan (PGN) stimulation. PmnAChR-2 expression was notably upregulated at 6 h after PGN challenge but had no response to LPS stimulation. Polyinosinic-polycytidylic acid challenge upregulated nearly all PmnAChRs, except for PmnAChR-5. Furthermore, Pm-miR-873-3p, Pm-miR-4577, Pm-miR-103a-3p, and Pm-miR-6753-3p were identified as the regulatory miRNA of PmnAChR-1, PmnAChR-3, PmnAChR-4, and PmnAChR-5, respectively. These findings suggested that these tandem arrays of nAChRs are unique to bivalves, and the tandem duplication of nAChR genes may be involved in the immune regulation process after pathogen stimulation.

MicroRNA-Mediated Regulation of Initial Host Responses in a Symbiotic Organ

One of the most important events in an animal's life history is the initial colonization by its microbial symbionts, yet little is known about this event's immediate impacts on the extent of host gene expression or the molecular mechanisms controlling it. MicroRNAs (miRNAs) are short, noncoding RNAs that bind to target mRNAs, rapidly shaping gene expression by posttranscriptional control of mRNA translation and decay. Here, we show that, in the experimentally tractable binary squid-vibrio symbiosis, colonization of the light organ induces extensive changes in the miRNA transcriptome. Examination of the squid genome revealed the presence of evolutionarily conserved genes encoding elements essential for the production and processing of miRNAs. At 24 h postcolonization, 215 host miRNAs were detected in the light organ, 26 of which were differentially expressed in response to the symbionts. A functional enrichment analysis of genes potentially targeted by downregulation of certain miRNAs at the initiation of symbiosis revealed two major gene ontology (GO) term categories, neurodevelopment and tissue remodeling. This symbiont-induced downregulation is predicted to promote these activities in host tissues and is consistent with the well-described tissue remodeling that occurs at the onset of the association. Conversely, predicted targets of upregulated miRNAs, including the production of mucus, are consistent with attenuation of immune responses by symbiosis. Taken together, our data provide evidence that, at the onset of symbiosis, host miRNAs in the light organ drive alterations in gene expression that (i) orchestrate the symbiont-induced development of host tissues, and (ii) facilitate the partnership by dampening the immune response.IMPORTANCE Animals often acquire their microbiome from the environment at each generation, making the initial interaction of the partners a critical event in the establishment and development of a stable, healthy symbiosis. However, the molecular nature of these earliest interactions is generally difficult to study and poorly understood. We report that, during the initial 24 h of the squid-vibrio association, a differential expression of host miRNAs is triggered by the presence of the microbial partner. Predicted mRNA targets of these miRNAs were associated with regulatory networks that drive tissue remodeling and immune suppression, two major symbiosis-induced developmental outcomes in this and many other associations. These results implicate regulation by miRNAs as key to orchestrating the critical transcriptional responses that occur very early during the establishment of a symbiosis. Animals with more complex microbiota may have similar miRNA-driven responses as their association is initiated, supporting an evolutionary conservation of symbiosis-induced developmental mechanisms.

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'.

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.

Integrated analysis of microRNA and mRNA expression profiles in Crassostrea gigas to reveal functional miRNA and miRNA-targets regulating shell pigmentation

MicroRNAs (miRNAs) regulate post-transcription gene expression by targeting genes and play crucial roles in diverse biological processes involving body color formation. However, miRNAs and miRNA-targets underlying shell color polymorphism remain largely unknown in mollusca. Using four shell colors full-sib families of the Pacific oyster Crassostrea gigas, we systematically identified miRNAs and miRNA-targets in the mantles, which organ could produce white, golden, black or partially pigmented shell. RNA sequencing and analysis identified a total of 53 known miRNA and 91 novel miRNAs, 47 of which were detected to differentially express among six pairwise groups. By integrating miRNA and mRNA expression profiles, a total of 870 genes were predicted as targets of differentially expressed miRNAs, mainly involving in biomineralization and pigmentation through functional enrichment. Furthermore, a total of four miRNAs and their target mRNAs were predicted to involve in synthesis of melanin, carotenoid or tetrapyrrole. Of them, lgi-miR-317 and its targets peroxidase and lncRNA TCONS_00951105 are implicated in acting as the competing endogenous RNA to regulate melanogenesis. Our studies revealed the systematic characterization of miRNAs profiles expressed in oyster mantle, which might facilitate understanding the intricate molecular regulation of shell color polymorphism and provide new insights into breeding research in oyster.

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.

Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in Crassostrea gigas infected by Ostreid herpesvirus 1

Background

Since 2008, the aquaculture production of Crassostrea gigas was heavily affected by mass mortalities associated to Ostreid herpesvirus 1 (OsHV-1) microvariants worldwide. Transcriptomic studies revealed the major antiviral pathways of the oyster immune response while other findings suggested that also small non-coding RNAs (sncRNA) such as microRNAs might act as key regulators of the oyster response against OsHV-1. To explore the explicit connection between small non-coding and protein-coding transcripts, we performed paired whole transcriptome analysis of sncRNA and messenger RNA (mRNA) in six oysters selected for different intensities of OsHV-1 infection.

Results

The mRNA profiles of the naturally infected oysters were mostly governed by the transcriptional activity of OsHV-1, with several differentially expressed genes mapping to the interferon, toll, apoptosis, and pro-PO pathways. In contrast, miRNA profiles suggested more complex regulatory mechanisms, with 15 differentially expressed miRNAs (DE-miRNA) pointing to a possible modulation of the host response during OsHV-1 infection. We predicted 68 interactions between DE-miRNAs and oyster 3′-UTRs, but only few of them involved antiviral genes. The sncRNA reads assigned to OsHV-1 rather resembled mRNA degradation products, suggesting the absence of genuine viral miRNAs.

Conclusions

We provided data describing the miRNAome during OsHV-1 infection in C. gigas. This information can be used to understand the role of miRNAs in healthy and diseased oysters, to identify new targets for functional studies and, eventually to disentangle cause and effect relationships during viral infections in marine mollusks.

A species-specific miRNA participates in biomineralization by targeting CDS regions of Prisilkin-39 and ACCBP in Pinctada fucata

Biomineralization is a sophisticated biological process precisely regulated by multiple molecules and pathways. Accumulating miRNAs have been identified in invertebrates but their functions in biomineralization are poorly studied. Here, an oyster species-specific miRNA, novel_miR_1 was found to regulate biomineralization in Pinctada fucata. Target prediction showed that novel_miR_1 could target Prisilkin-39 and ACCBP by binding to their coding sequences (CDS). Tissue distribution analysis revealed that the expression level of novel_miR_1 was highest in the mantle, which was a key tissue participating in biomineralization. Gain-of-function assay in vivo showed that biomineralization-related genes including Prisilkin-39 and ACCBP were down-regulated and shell inner surfaces of both prismatic and nacreous layer were disrupted after the over-expression of novel_miR_1, indicating its dual roles in biomineralization. Furthermore, the shell notching results indicated that novel_miR_1 was involved in shell regeneration. Dual-luciferase reporter assay in vitro demonstrated that novel_miR_1 directly suppressed Prisilkin-39 and ACCBP genes by binding to the CDS regions. Taken together, these results suggest that novel_miR_1 is a direct negative regulator to Prisilkin-39 and ACCBP and plays an indispensable and important role in biomineralization in both prismatic and nacreous layer of P. fucata.

Evolutionary and functional analysis of MyD88 genes in pearl oyster Pinctada fucata martensii

Myeloid differentiation factor 88 (MyD88) is an adapter protein that links toll-like receptor and interleukin 1 receptor-mediated signal transduction. In this study, we identified 20 MyD88 genes from eight mollusk genomes and found that MyD88 was expanded in bivalves. This expansion tends to be tandem duplication. Phylogenetic analysis suggested that the tandem duplication of MyD88 was formed before bivalve differentiation. All of the identified MyD88 contained both of death domain (DD) and toll/interleukin-1 receptor (TIR) domain, and 13 mollusks MyD88 have low complexity regions (LCRs), which were not found in the MyD88 from humans and zebrafish. The genomic structure showed that most of the mollusk MyD88 (14 of 19) contained five conserved introns, four of which were found in humans and zebrafish. Furthermore, the cDNA full length of PfmMyD88-2 (one of the two identified MyD88 in Pincatada fucata martensii) was obtained with 1591 bp, including 260 bp of 5'UTR, 257 bp of 3'UTR, and 1077 bp of open reading frame encoding 358 amino acids. Quantitative real-time PCR analysis demonstrated that PfmMyD88-2 mRNA was widely expressed in all detected tissues. The highest expression level was in the gills and followed by hepatopancreas and feet. After lipopolysaccharide stimulation, PfmMyD88-2 expression level increased and reached the highest level at 12 h and then gradually declined to the normal level. Over-expression of PfmMyD88-2 in HEK293T increased the luciferase activity of the pNF-κB-Luc reporter. We also identified that PfmmiR-4047 could regulate the expression of PfmMyD88-2. These results help us elucidate the mechanism underlying mollusk immune response.

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.

PfmPif97-like regulated by Pfm-miR-9b-5p participates in shell formation in Pinctada fucata martensii

Mollusk shell matrix proteins are important for the formation of organic frameworks, crystal nucleation, and crystal growth in Pinctada fucata martensii (P. f. martensii). MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play important roles in many biological processes, including shell formation. In this study, we obtained the full-length sequence of Pif97-like gene in P. f. martensii (PfmPif97-like). PfmPif97-like was mainly distributed in mantle pallial and mantle edge. Correlation analysis indicated that the average shell thickness and weight showed a positive correlation with PfmPif97-like expression (P < 0.05). The inner surface of the nacreous layer and prismatic layer showed atypical growth when we knocked down the expression of PfmPif97-like by RNA interference (RNAi). We used a luciferase reporter assay to identify that miR-9b-5p of P. f. martensii (Pfm-miR-9b-5p) downregulated the expression of PfmPif97-like by interacting with the 3′-untranslated region (UTR) while we obtained the same result by injecting the Pfm-miR-9b-5p mimics in vivo. After injecting the mimics, we also observed abnormal growth in nacre layer and prismatic layer which is consistent with the result of RNAi. We proposed that PfmPif97-like regulated by Pfm-miR-9b-5p participates in shell formation of P. f. martensii. These findings provide important clues about the molecular mechanisms that regulate biomineralization in P. f. martensii.

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.

Identification and characteristics of muscle growth-related microRNA in the Pacific abalone, Haliotis discus hannai

BACKGROUND

The Pacific abalone, Haliotis discus hannai, is the most important cultivated abalone in China. Improving abalone muscle growth and increasing the rate of growth are important genetic improvement programs in this industry. MicroRNAs are important small noncoding RNA molecules that regulate post-transcription gene expression. However, no miRNAs have been reported to regulate muscle growth in H. discus hannai.

RESULTS

we profiled six small RNA libraries for three large abalone individuals (L_HD group) and three small individuals (S_HD group) using RNA sequencing technology. A total of 205 miRNAs, including 200 novel and 5 known miRNAs, were identified. In the L_HD group, 3 miRNAs were up-regulated and 7 were down-regulated compared to the S_HD specimens. Bioinformatics analysis of miRNA target genes revealed that miRNAs participated in the regulation of cellular metabolic processes, the regulation of biological processes, the Wnt signaling pathway, ECM-receptor interaction, and the MAPK signaling pathway, which are associated with regulating growth. Bone morphogenetic protein 7 (BMP7) was verified as a target gene of hdh-miR-1984 by a luciferase reporter assay and we examined the expression pattern in different developmental stages.

CONCLUSION

This is the first study to demonstrate that miRNAs are related to the muscle growth of H. discus hannai. This information could be used to study the mechanisms of abalone muscle growth. These DE-miRNAs may be useful as molecular markers for functional genomics and breeding research in abalone and closely related species.

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.

SmithRNAs: Could Mitochondria "Bend" Nuclear Regulation?

Typically, animal mitochondria have very compact genomes, with few short intergenic regions, and no introns. Hence, it may seem that there is little space for unknown functions in mitochondrial DNA (mtDNA). However, mtDNA can also operate through RNA interference, as small non coding RNAs (sncRNAs) produced by mtDNA have already been proposed for humans. We sequenced sncRNA libraries from isolated mitochondria of Ruditapes philippinarum (Mollusca Bivalvia) gonads, a species with doubly uniparental inheritance of mitochondria, and identified several putative sncRNAs of mitochondrial origin. Some sncRNAs are transcribed by intergenic regions that form stable stem-hairpin structures, which makes them good miRNA-like candidates. We decided to name them small mitochondrial highly-transcribed RNAs (smithRNAs). Many concurrent data support that we have recovered sncRNAs of mitochondrial origin that might be involved in gonad formation and able to affect nuclear gene expression. This possibility has been never suggested before. If mtDNA can affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for it to interact with the nucleus, and makes metazoan mtDNA a much more complex genome than previously thought.

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.

An oyster species-specific miRNA scaffold42648_5080 modulates haemocyte migration by targeting integrin pathway

miRNAs are important gene regulators at post-transcriptional level and can modulate diverse biological processes, including immune response. Dozens of species-specific miRNAs have been identified in oyster Crassostrea gigas while their functions remain largely unknown. In the present study, an oyster species-specific miRNA scaffold42648_5080 was found responsive to LPS stimulation and might target a total of 31 oyster genes possibly involved in cell communication, cellular localization and cellular response to stimulus. Besides, in gain-of-function assay of scaffold42648_5080 in vivo, the phagocytosis (30.90% in miRNA group verse 23.20% in miRNA control group), apoptosis (3.10% in miRNA group verse 5.30% in miRNA control group) and migration rate (13.88% in miRNA group verse 21.03% in miRNA control group) of oyster haemocytes were found significantly altered after the injection of scaffold42648_5080 mimics. Among the target genes, integrin-linked kinase (CgILK) was considered crucial in cell migration and its interaction with scaffold42648_5080 was then verified both in vitro and in vivo. Consequently, a significant decrease of relative luciferase ratio was observed in CgILK 3'-UTR luciferase reporter assay after transfection of scaffold42648_5080 mimics (0.70-fold of that in blank group, p < 0.01). Meanwhile, when scaffold42648_5080 was overexpressed in vivo (5.41-fold of miRNA control group, p < 0.01), the expression of CgILK declined significantly to 0.25-fold of miRNA control group (p < 0.01). Comparatively, a significant decrease of the haemocyte migration rate (19.76% verse 34.82% in siEGFP control group, p < 0.01) was observed after knock-down of CgILK in vivo. The present study, as far as we know, for the first time revealed the immunomodulation role of an oyster species-specific miRNA, which might provide new insights into miRNA-mediated adaptation mechanism of oysters.

Long Non-Coding RNAs (lncRNAs) of Sea Cucumber: Large-Scale Prediction, Expression Profiling, Non-Coding Network Construction, and lncRNA-microRNA-Gene Interaction Analysis of lncRNAs in Apostichopus japonicus and Holothuria glaberrima During LPS Challenge and Radial Organ Complex Regeneration

Long non-coding RNA (lncRNA) structurally resembles mRNA but cannot be translated into protein. Although the systematic identification and characterization of lncRNAs have been increasingly reported in model species, information concerning non-model species is still lacking. Here, we report the first systematic identification and characterization of lncRNAs in two sea cucumber species: (1) Apostichopus japonicus during lipopolysaccharide (LPS) challenge and in heathy tissues and (2) Holothuria glaberrima during radial organ complex regeneration, using RNA-seq datasets and bioinformatics analysis. We identified A. japonicus and H. glaberrima lncRNAs that were differentially expressed during LPS challenge and radial organ complex regeneration, respectively. Notably, the predicted lncRNA-microRNA-gene trinities revealed that, in addition to targeting protein-coding transcripts, miRNAs might also target lncRNAs, thereby participating in a potential novel layer of regulatory interactions among non-coding RNA classes in echinoderms. Furthermore, the constructed coding-non-coding network implied the potential involvement of lncRNA-gene interactions during the regulation of several important genes (e.g., Toll-like receptor 1 [TLR1] and transglutaminase-1 [TGM1]) in response to LPS challenge and radial organ complex regeneration in sea cucumbers. Overall, this pioneer systematic identification, annotation, and characterization of lncRNAs in echinoderm pave the way for similar studies and future genetic, genomic, and evolutionary research in non-model species.

Identification of a novel miR-146a from Pinctada martensii involved in the regulation of the inflammatory response

Increasing evidence demonstrated that microRNAs (miRNAs) play critical roles in innate immunity in vertebrates and invertebrates. MiR-146a/b is reported as a key regulator of the immune response through mediating Toll-like receptor and cytokine signalling. In this study, a novel miR-146a was identified and characterised from Pinctada martensii (designated as pm-miR-146a), and its roles in modulating the inflammatory response after LPS stimulation were also investigated. Pm-miR-146a ubiquitously expressed in all examined tissues, with the highest level in the mantle and lowest expression in the haemolymph. Pm-miR-146a increased at 24 h after lipoplysaccharide injection, in union with up-regulated NF-κB (P < 0.05). The over-expression of pm-miR-146a in vivo could significantly inhibit the expression of macrophage migration inhibitory factor (MIF), the potential target gene predicted by miRanda, while enforcing pm-miR-146a involved in the down-regulation of NF-κB. Thus, we propose that pm-miR-146a plays a role of negative feedback regulation to the NF-κB signal by repressing the expression of the pro-inflammatory cytokine MIF. These findings revealed that miR-146a represents a critical role in inflammatory response and offers new evidence for miRNAs in the innate immunity of molluscs.

The miRNA biogenesis in marine bivalves

Small non-coding RNAs include powerful regulators of gene expression, transposon mobility and virus activity. Among the various categories, mature microRNAs (miRNAs) guide the translational repression and decay of several targeted mRNAs. The biogenesis of miRNAs depends on few gene products, essentially conserved from basal to higher metazoans, whose protein domains allow specific interactions with dsRNA. Here, we report the identification of key genes responsible of the miRNA biogenesis in 32 bivalves, with particular attention to the aquaculture species Mytilus galloprovincialis and Crassostrea gigas. In detail, we have identified and phylogenetically compared eight evolutionary conserved proteins: DROSHA, DGCR8, EXP5, RAN, DICER TARBP2, AGO and PIWI. In mussels, we recognized several other proteins participating in the miRNA biogenesis or in the subsequent RNA silencing. According to digital expression analysis, these genes display low and not inducible expression levels in adult mussels and oysters whereas they are considerably expressed during development. As miRNAs play an important role also in the antiviral responses, knowledge on their production and regulative effects can shed light on essential molecular processes and provide new hints for disease prevention in bivalves.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

Expression Profile Analysis of miR-221 and miR-222 in Different Tissues and Head Kidney Cells of Cynoglossus semilaevis, Following Pathogen Infection

Half-smooth tongue sole (Cynoglossus semilaevis) is an important marine commercial fish species in China, which suffers from widespread disease outbreaks. Recently, in this regard, our group identified immune-related microRNAs (miRNAs) of C. semilaevis following Vibrio anguillarum infection. Furthermore, miRNA microarray was utilized to characterize the immune roles of important miRNA candidates in response to bacterial infection. Therefore, in the present study, we characterized miR-221 and miR-222 and profiled their expression after challenge. Here, miR-221 and miR-222 precursors were predicted to have a typical hairpin structure. Both miRNAs were expressed in a broad range of tissues in C. semilaevis, while miR-221 and miR-222 were significantly differentially expressed in the immune tissues of C. semilaevis among three small RNA libraries [control group (CG), bacteria-challenged fish without obvious symptoms of infection (NOSG), and bacteria-challenged fish with obvious symptoms of infection (HOSG)]. In order to further characterize and understand the immune response of miR-221 and miR-222, therefore, we profiled miR-221 and miR-222 expression in selected immune tissues after challenge with V. anguillarum. Both miR-221 and miR-222 were upregulated in the liver and spleen, while different expression patterns were observed in the head kidney. In addition, in half-smooth tongue sole head kidney cell line after challenge with lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV), both miR-221 and miR-222 showed significant difference in expression response to pathogen. Meanwhile, the target gene of miR-221 and miR-222 was predicted, which indicated that tumor necrosis factor receptor-associated factor 6 (TRAF6) and interleukin-1 beta (IL-1β) were the target genes of miR-221 and miR-222, respectively. Collectively, these findings indicated that miR-221 and miR-222 have putative roles in innate immune response during C. semilaevis exposure to pathogens. Our findings could expand the knowledge of immune function of C. semilaevis miRNA and guide future studies on C. semilaevis immunity.

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.

The comprehensive immunomodulation of NeurimmiRs in haemocytes of oyster Crassostrea gigas after acetylcholine and norepinephrine stimulation

Background

Neural-endocrine-immune (NEI) system is a major modulation network among the nervous, endocrine and immune system and weights greatly in maintaining homeostasis of organisms during stress and infection. Some microRNAs are found interacting with NEI system (designated NeurimmiRs), addressing swift modulations on immune system. The oyster Crassostrea gigas, as an intertidal bivalve, has evolved a primary NEI system. However, the knowledge about NeurimmiRs in oysters remains largely unknown.

Results

Six small RNA libraries from haemocytes of oysters stimulated with acetylcholine (ACh) and norepinephrine (NE) were sequenced to identify neurotransmitter-responsive miRNAs and survey their immunomodulation roles. A total of 331 miRNAs (132 identified in the present study plus 199 identified previously) were subjected to expression analysis, and twenty-one and sixteen of them were found ACh- or NE-responsive, respectively (FDR < 0.05). Meanwhile, 21 miRNAs exhibited different expression pattern after ACh or NE stimulation. Consequently, 355 genes were predicted as putative targets of these neurotransmitter-responsive miRNAs in oyster. Through gene onthology analysis, multiple genes involved in death, immune system process and response to stimulus were annotated to be modulated by NeurimmiRs. Besides, a significant decrease in haemocyte phagocytosis and late-apoptosis or necrosis rate was observed after ACh and NE stimulation (p < 0.05) while early-apoptosis rate remained unchanged.

Conclusions

A comprehensive immune-related network involving PRRs, intracellular receptors, signaling transducers and immune effectors was proposed to be modulated by ACh- and NE-responsive NeurimmiRs, which would be indispensable for oyster haemocytes to respond against stress and infection. Characterization of the NeurimmiRs would be an essential step to understand the NEI system of invertebrate and the adaptation mechanism of oyster.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2150-8) contains supplementary material, which is available to authorized users.

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.

Expression profiling analysis of the microRNA response of Cynoglossus semilaevis to Vibrio anguillarum and other stimuli

To investigate the roles of microRNAs (miRNA) of Cynoglossus semilaevis in response to Vibrio anguillarum that were previously identified using high-throughput sequencing, microarray analyses was performed on three small RNA libraries (CG, NOSG, and HOSG) prepared from C. semilaevis immune tissues. In total, of 1279 designed probes, 739 (57.78 %) were detectable. The expression levels of these miRNAs were analyzed using pairwise comparisons among the three libraries, and a total of 99 miRNAs were observed to be significantly differentially expressed. The expression patterns of 10 differentially expressed miRNAs were validated by real-time quantitative PCR (RT-qPCR). In addition, expression of miR-142-5p, miR-223, and miR-181a in response to V. anguillarum at numerous time-points in four tissues, as well as the responses to lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV) in head kidney cells, were studied by qRT-PCR. Taken together, all of the expression profiles showed significant differences compared to the control group; both similarities and differences in the expression responses to the same pathogen were observed. Collectively, these findings highlighted the putative roles for miRNAs in the context of the innate immune response of C. semilaevis exposing to pathogens and that further studies are needed to understand the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.

Identification and profiling of microRNAs in goat endometrium during embryo implantation

Background

MicroRNAs (miRNAs) are short, highly conserved small noncoding RNAs that had fundamental roles in post-transcriptional gene expression, and they are crucial for proper control of biological processes and known to participate in embryo implantation. However, miRNA expression profiles in the pre-receptive and receptive phases of the goat endometrium during embryo implantation are unknown.

Results

A total of 1,069 and 847 miRNAs were expressed in receptive (R) and pre-receptive (P) goat endometrium, and 632 miRNAs were co-expressed in both phases. We identified 545 (50.98%) known miRNAs in the R library and 522 (61.63%) in the P library. There were 110 up-expressed miRNAs and 33 down-expressed miRNAs in receptive endometrium compared with the pre-receptive endometrium meeting the criteria of P-values< 0.05. Moreover, GO and KEGG analysis of the target genes of the differentially expressed miRNAs revealed some candidate miRNAs, genes and pathways that may involve in the formation of the receptive endometrium. Based on stem-loop RT-qPCR, 15 miRNAs were detected and the results suggested that the majority of the miRNA expression data measured by Solexa deep sequencing could represent actual miRNA expression levels.

Conclusions

Our data revealed the first miRNA profile related to the biology of the goat receptive endometrium during embryo implantation, and the results suggested that a subset of miRNAs might play important roles in the formation of endometrial receptivity. Thus, elucidating the physiological roles of endometrial miRNAs will help us better understand the genetic control of embryo implantation in goats.

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.

The microRNA signature in response to nutrient restriction and refeeding in skeletal muscle of Chinese perch (Siniperca chuatsi)

The Chinese perch (Siniperca chuatsi) is one of the most commercially important carnivorous fish species in aquaculture with its large-scale culture in China. Increasing evidence suggests that microRNAs (miRNAs) play an important role in muscle cell proliferation and differentiation. However, the knowledge of the identity of myogenic miRNAs and the effect of nutrient status on miRNA expression in teleost remains limited. In the present study, among the 21 miRNAs identified with high abundance in the fast muscle of adult Chinese perch, 19 miRNAs were differentially expressed in the adults and juveniles. The postprandial changes in the transcript abundance were determined for the 21 miRNAs following a single satiating meal in the juveniles after fasting for 1 week. The results showed that the seven miRNAs (miR-10c, miR-107a, miR-133a-3p, miR-140-3p, miR-181a-5p, miR-206, and miR-214) were sharply upregulated or downregulated within 1 h after refeeding. These miRNAs may be the promising candidate miRNAs involved in a fast-response signaling system that regulates fish skeletal muscle growth. Target prediction and expressional analysis suggested that four miRNAs (miR-10c, miR-107a, miR-140-3p, and miR-181a-5p) might play a role in regulating the translation of target gene transcripts such as myostatin following acute anabolic stimuli.

Identification of conserved and novel microRNAs in the Pacific oyster Crassostrea gigas by deep sequencing

MicroRNAs (miRNAs) play important roles in regulatory processes in various organisms. To date many studies have been performed in the investigation of miRNAs of numerous bilaterians, but limited numbers of miRNAs have been identified in the few species belonging to the clade Lophotrochozoa. In the current study, deep sequencing was conducted to identify the miRNAs of Crassostrea gigas (Lophotrochozoa) at a genomic scale, using 21 libraries that included different developmental stages and adult organs. A total of 100 hairpin precursor loci were predicted to encode miRNAs. Of these, 19 precursors (pre-miRNA) were novel in the oyster. As many as 53 (53%) miRNAs were distributed in clusters and 49 (49%) precursors were intragenic, which suggests two important biogenetic sources of miRNAs. Different developmental stages were characterized with specific miRNA expression patterns that highlighted regulatory variation along a temporal axis. Conserved miRNAs were expressed universally throughout different stages and organs, whereas novel miRNAs tended to be more specific and may be related to the determination of the novel body plan. Furthermore, we developed an index named the miRNA profile age index (miRPAI) to integrate the evolutionary age and expression levels of miRNAs during a particular developmental stage. We found that the swimming stages were characterized by the youngest miRPAIs. Indeed, the large-scale expression of novel miRNAs indicated the importance of these stages during development, particularly from organogenetic and evolutionary perspectives. Some potentially important miRNAs were identified for further study through significant changes between expression patterns in different developmental events, such as metamorphosis. This study broadened the knowledge of miRNAs in animals and indicated the presence of sophisticated miRNA regulatory networks related to the biological processes in lophotrochozoans.

Identification and comparative analysis of the Tegillarca granosa haemocytes microRNA transcriptome in response to Cd using a deep sequencing approach

Background

MicroRNAs (miRNAs) are endogenous non-coding small RNAs (sRNAs) that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. To identify miRNAs regulated by heavy metal stress, we constructed two sRNA libraries for the blood clam Tegillarca granosa: one for organisms exposed to toxic levels of cadmium (Cd) and one for a control group.

Results

Sequencing of the two libraries and subsequent analysis revealed 215 conserved and 39 new miRNAs. Most of the new miRNAs in T. granosa were up- or down-regulated in response to Cd exposure. There were significant differences in expression between the Cd and control groups for 16 miRNAs. Of these, five miRNAs were significantly up-regulated and 11 were significantly down-regulated in the Cd stress library. Potential targets were predicted for the 16 differential miRNAs in pre-miRNAs identified according to sequence homology. Some of the predicted miRNA targets are associated with regulation of the response to stress induced by heavy metals. Five differentially expressed miRNAs (Tgr-nmiR-8, Tgr-nmiR-21, Tgr-miR-2a, Tgr-miR-10a-5p, and Tgr-miR-184b) were validated by qRT-PCR.

Conclusion

Our study is the first large-scale identification of miRNAs in T. granosa haemocytes. Our findings suggest that some miRNAs and their target genes and pathways may play critical roles in the responses of this species to environmental heavy metal stresses.

Identification of microRNAs in the coral Stylophora pistillata

Coral reefs are major contributors to marine biodiversity. However, they are in rapid decline due to global environmental changes such as rising sea surface temperatures, ocean acidification, and pollution. Genomic and transcriptomic analyses have broadened our understanding of coral biology, but a study of the microRNA (miRNA) repertoire of corals is missing. miRNAs constitute a class of small non-coding RNAs of ∼22 nt in size that play crucial roles in development, metabolism, and stress response in plants and animals alike. In this study, we examined the coral Stylophora pistillata for the presence of miRNAs and the corresponding core protein machinery required for their processing and function. Based on small RNA sequencing, we present evidence for 31 bona fide microRNAs, 5 of which (miR-100, miR-2022, miR-2023, miR-2030, and miR-2036) are conserved in other metazoans. Homologues of Argonaute, Piwi, Dicer, Drosha, Pasha, and HEN1 were identified in the transcriptome of S. pistillata based on strong sequence conservation with known RNAi proteins, with additional support derived from phylogenetic trees. Examination of putative miRNA gene targets indicates potential roles in development, metabolism, immunity, and biomineralisation for several of the microRNAs. Here, we present first evidence of a functional RNAi machinery and five conserved miRNAs in S. pistillata, implying that miRNAs play a role in organismal biology of scleractinian corals. Analysis of predicted miRNA target genes in S. pistillata suggests potential roles of miRNAs in symbiosis and coral calcification. Given the importance of miRNAs in regulating gene expression in other metazoans, further expression analyses of small non-coding RNAs in transcriptional studies of corals should be informative about miRNA-affected processes and pathways.