Comparative microRNA-seq Analysis Depicts Candidate miRNAs Involved in Skin Color Differentiation in Red Tilapia

The lgi-miR-2d is Potentially Involved in Shell Melanin Synthesis by Targeting mitf in Manila Clam Ruditapes philippinarum

Shell color as an important economic trait is also the crucial target trait for breeding and production. MicroRNA (miRNA) is an endogenous small non-coding RNA that can post-transcriptionally regulate the expression of target genes, it plays important roles in many life activities and physiological processes, such as shell color, stress response, and disease traits. In this study, we investigated the function of lgi-miR-2d in shell melanin formation and the expression patterns of lgi-miR-2d and target gene Rpmitf in Manila clam Ruditapes philippinarum. We further explored and verified the relationship between Rpmitf and lgi-miR-2d and identified the expression level of shell color-related gene changes by RNAi and injecting the antagomir of lgi-miR-2d, respectively. Our results indicated that lgi-miR-2d antagomir affected the expression of its target gene Rpmitf. In addition, the dual-luciferase reporter assay was conducted to confirm the direct interaction between lgi-miR-2d and Rpmitf. The results showed that the expression levels of melanin-related genes such as Rpmitf and tyr were significantly decreased in the positive treatment group compared with the blank control group after the Rpmitf dsRNA injection, indicating Rpmitf plays a crucial role in the melanin synthesis pathway. Taken together, we speculated that lgi-miR-2d might be negatively modulating Rpmitf, which might regulate other shell color-related genes, thereby affecting melanin synthesis in R. philippinarum.

Transcriptome analysis reveals candidate miRNAs involved in skin color differentiation of juvenile Plectropomus leopardus in response to different background colors

Red skin color in Plectropomus leopardus is important to its ornamental and economic value. However, the color of P. leopardus can change during the rearing process, darkening and turning black due to the influence of environmental background color. The underlying molecular mechanisms that regulate this phenomenon remain unclear. MicroRNAs (miRNAs) are endogenous, small non-coding RNAs that play important roles in numerous biological processes, such as skin differentiation and color formation in many animals. Therefore, we performed miRNA sequencing of P. leopardus skin before (initial) and after rearing with three different background colors (white, black, and blue) using Illumina sequencing to identify candidate miRNAs that may contribute to skin color differentiation. In total, 154,271,376 clean reads were obtained, with over 92 % of them successfully mapped to the P. leopardus reference genome. The miRNA length distributions of all samples displayed peaks around a typical length of 22 nt. Within these sequences, 243 known and 287 novel miRNAs were identified. A total of 65 significantly differentially expressed miRNAs (DEMs) were identified (P < 0.05), including 40 known DEMs and 25 novel DEMs. These DEMs included novel_561, miR-141-3p, and miR-129-5p, whose target genes were primarily associated with pigmentation related processes, including tyrosine metabolism, melanogenesis, and the Wnt signaling pathway. These findings shed light on the potential roles of miRNAs in the darkening of skin color in P. leopardus, thus enhancing our understanding of the molecular mechanisms involved in skin pigmentation differentiation in this species.

Genome-wide association study (GWAS) analysis of black color trait in the leopard coral grouper (Plectropomus leopardus) using whole genome resequencing

The leopard coral grouper (Plectropomus leopardus) is a coral reef fish species that exhibits rapid and diverse color variation. However, the presence of melanoma and the high proportion of individuals displaying black color in artificial breeding have led to reduced economic and ornamental value. To pinpoint single nucleotide polymorphisms (SNPs) and potential genes linked to the black pigmentation characteristic in this particular species, This study gathered a cohort of 360 specimens from diverse origins and conducted a comprehensive genome-wide association analysis (GWAS) employing whole-genome resequencing. As a result, 57 SNPs related to the black skin trait were identified, and a grand total of 158 genes were annotated within 50 kb of these SNPs. Subsequently, GWAS was applied to three populations (LED, QHH, and QHL), and the corresponding results were compared with the analysis results of the total population. The results of the four GWAS models showed significant enrichment in Rap1 signaling pathway, melanin biosynthesis, metabolic pathways, tyrosine metabolism, cAMP signaling pathway, AMPK signaling pathway, PI3K-Akt signaling pathway, EGFR tyrosine kinase inhibitor resistance, HIF-1 signaling pathway, Ras signaling pathway, MAPK signaling pathway, etc. (p < 0.05), which were mainly associated with eleven genes (POL4, MET, E2F2, COMT, ZBED1, TYRP2, FOXP2, THIKA, LORF2, MYH16 and SOX2). Significant differences (p < 0.05) were observed in the expression of all 11 genes in the dorsal skin tissue, in 10 genes except COMT in the ventral skin tissue, and in all 11 genes in the caudal fin tissue. These findings imply that the control of body color in the P. leopardus is the result of the joint action of multiple genes and signaling pathways. These findings will contribute to a more profound comprehension of the genetic attributes that underlie the development of black skin in the vibrant P. leopardus, thus furnishing a theoretical foundation for genetic enhancement.

microRNA regulation of skin pigmentation in golden-back mutant of crucian carp from a rice-fish integrated farming system

BACKGROUND

MicroRNAs (miRNAs) are endogenous small non-coding RNAs (21-25 nucleotides) that act as essential components of several biological processes. Golden-back crucian carp (GBCrC, Carassius auratus) is a naturally mutant species of carp that has two distinct body skin color types (golden and greenish-grey), making it an excellent model for research on the genetic basis of pigmentation. Here, we performed small RNA (sRNA) analysis on the two different skin colors via Illumina sequencing.

RESULTS

A total of 679 known miRNAs and 254 novel miRNAs were identified, of which 32 were detected as miRNAs with significant differential expression (DEMs). 23,577 genes were projected to be the targets of 32 DEMs, primarily those involved in melanogenesis, adrenergic signaling in cardiomyocytes, MAPK signaling pathway and wnt signaling pathway by functional enrichment. Furthermore, we built an interaction module of mRNAs, proteins and miRNAs based on 10 up-regulated and 13 down-regulated miRNAs in golden skin. In addition to transcriptional destabilization and translational suppression, we discovered that miRNAs and their target genes were expressed in the same trend at both the transcriptional and translational levels. Finally, we discovered that miR-196d could be indirectly implicated in regulating melanocyte synthesis and motility in the skin by targeting to myh7 (myosin-7) gene through the luciferase reporter assay, antagomir silencing in vivo and qRT-PCR techniques.

CONCLUSIONS

Our study gives a systematic examination of the miRNA profiles expressed in the skin of GBCrC, assisting in the comprehension of the intricate molecular regulation of body color polymorphism and providing insights for C. auratus breeding research.

Comparative transcriptome and miRNA analysis of skin pigmentation during embryonic development of Chinese soft-shelled turtle (Pelodiscus sinensis)

BACKGROUND

Aquatic animals show diverse body coloration, and the formation of animal body colour is a complicated process. Increasing evidence has shown that microRNAs (miRNAs) play important regulatory roles in many life processes. The role of miRNAs in pigmentation has been investigated in some species. However, the regulatory patterns of miRNAs in reptile pigmentation remain to be elucidated. In this study, we performed an integrated analysis of miRNA and mRNA expression profiles to explore corresponding regulatory patterns in embryonic body colour formation in the soft-shelled turtle Pelodiscus sinensis.

RESULTS

We identified 8 866 novel genes and 9 061 mature miRNAs in the skin of Chinese soft-shelled turtles in three embryonic stages (initial period: IP, middle period: MP, final period: FP). A total of 16 563 target genes of the miRNAs were identified. Furthermore, we identified 2 867, 1 840 and 4 290 different expression genes (DEGs) and 227, 158 and 678 different expression miRNAs (DEMs) in IP vs. MP, MP vs. FP, and IP vs. FP, respectively. Among which 72 genes and 25 miRNAs may be related to turtle pigmentation in embryonic development. Further analysis of the novel miRNA families revealed that some novel miRNAs related to pigmentation belong to the miR-7386, miR-138, miR-19 and miR-129 families. Novel_miR_2622 and novel_miR_2173 belong to the miR-19 family and target Kit and Gpnmb, respectively. The quantification of novel_miR_2622 and Kit revealed negative regulation, indicating that novel_miR_2622 may participate in embryonic pigmentation in P. sinensis by negatively regulating the expression of Kit.

CONCLUSIONS

miRNA act as master regulators of biological processes by controlling the expression of mRNAs. Considering their importance, the identified miRNAs and their target genes in Chinese soft-shelled turtle might be useful for investigating the molecular processes involved in pigmentation. All the results of this study may aid in the improvement of P. sinensis breeding traits for aquaculture.

MiR-382 Functions on the Regulation of Melanogenesis via Targeting dct in Rainbow Trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that are involved in a diverse collection of biological processes as important post-transcriptional regulators. However, little is known about the molecular regulatory mechanism of miRNAs in fish skin pigmentation. In this study, we first confirmed that dopachrome tautomerase (dct), a key gene of melanogenesis, is a target of miR-382 in rainbow trout (Oncorhynchus mykiss) using luciferase reporter assay. The analysis of different developmental stages and tissue expression patterns between wild-type and yellow mutant rainbow trout suggested that miR-382 is a potential regulator during the process of skin pigmentation. In vitro, miR-382 mimics in rainbow trout primary liver cells significantly downregulated dct expression and resulted in decreased expression of key melanogenic genes including tyrosine-related protein 1 (tyrp1) and premelanosome protein (pmel), whereas the expression level of dct was markedly increased after transfected with miR-382 inhibitor. In vivo, overexpression of miR-382 by injection of miR-382 agomir significantly depressed the expression of dct in dorsal skin, tail fin, and liver and then reduced the expression levels of tyrp1 and pmel. Furthermore, transfection of miR-382 mimics inhibited cell proliferation and induced apoptosis. Taken together, our results identified a functional role of miR-382 in rainbow trout skin pigmentation through targeting dct, which facilitate understanding the regulatory mechanism of rainbow trout skin color at the post-transcriptional level and provide a theoretical basis for molecular breeding with skin color as the target trait.

Drug response biomarkers of Pien Tze Huang treatment for hepatic fibrosis induced by carbon tetrachloride

OBJECTIVE

To explore biomarkers of Pien Tze Huang that ameliorated the symptoms of hepatic fibrosis.

METHODS

Two groups of carbon tetrachloride-induced hepatic fibrosis (HF) mice model were constructed in our study: one group received PZH treatment and another group received no treatment. We performed this study to investigate the role of PZH in the regulation process of hepatic fibrosis.

RESULTS

We identified 31 down-regulated and 39 up-regulated miRNAs using small RNA-seq analysis. Combining RNA-Seq data analysis, our study revealed 7 significant target genes (Sp4, Slc2a6, Tln2, Hmga2, Ank3, Pax9, Fgf9). The results of real-time quantitative polymerase chain reaction analysis suggested that the expression level of 6 genes (Sp4, Tln2, Hmga2, Ank3, Pax9, Fgf9) were down-regulated compared to control group. On the other hand, the expression level of Slc2a6 appeared to be up-regulated. The protein mass spectrometry showed that PZH group had lower protein expression of Tln2 compared to control group.

CONCLUSION

We identified 7 genes that were significantly related to PZH response in HF mice using multiple conjoint analysis methods. These genes could participate in underlying regulation mechanism of hepatic fibrosis during PZH treatment.

Comparing Transcriptomes Reveals Key Metabolic Mechanisms in Superior Growth Performance Nile Tilapia (Oreochromis niloticus)

Metabolic capacity is intrinsic to growth performance. To investigate superior growth performance in Nile tilapia, three full-sib families were bred and compared at the biochemical and transcriptome levels to determine metabolic mechanisms involved in significant growth differences between individuals under the same culture environment and feeding regime. Biochemical analysis showed that individuals in the higher growth group had significantly higher total protein, total triglyceride, total cholesterol, and high- and low-density lipoproteins, but significantly lower glucose, as compared with individuals in the lower growth group. Comparative transcriptome analysis showed 536 differentially expressed genes (DEGs) were upregulated, and 622 DEGs were downregulated. These genes were significantly enriched in three key pathways: the tricarboxylic acid cycle (TCA cycle), fatty acid biosynthesis and metabolism, and cholesterol biosynthesis and metabolism. Conjoint analysis of these key pathways and the biochemical parameters suggests that Nile tilapia with superior growth performance have higher ability to consume energy substrates (e.g., glucose), as well as higher ability to biosynthesize fatty acids and cholesterol. Additionally, the fatty acids biosynthesized by the superior growth performance individuals were less active in the catabolic pathway overall, but were more active in the anabolic pathway, and might be used for triglyceride biosynthesis to store excess energy in the form of fat. Furthermore, the tilapia with superior growth performance had lower ability to convert cholesterol into bile acids, but higher ability to convert it into sterols. We discuss the molecular mechanisms of the three key metabolic pathways, map the pathways, and note key factors that may impact the growth of Nile tilapia. The results provide an important guide for the artificial selection and quality enhancement of superior growth performance in tilapia.

Transcriptomic and proteomic analyses reveal the common and unique pathway(s) underlying different skin colors of leopard coral grouper (Plectropomus leopardus)

To gain a comprehensive and unbiased molecular understanding of the different skin colors of P. leopardus, we used Illumina HiSeq 2500 and TMT (Tandem Mass Tag) to compare transcription and protein levels between red and black skin of P. leopardus. We identified 797 upregulated and 314 downregulated genes (differentially expressed genes; DEGs) in red (RG) compared with black (BG) skin of P. leopardus. We also identified 377 differentially abundant proteins (DAPs), including 314 upregulated and 63 downregulated proteins. These DEGs and DAPs were significantly enriched in melanin synthesis (e.g., pyrimidine metabolism, Phenylalanine, tyrosine, and tryptophan biosynthesis, melanogenesis, phenylalanine metabolism, and tyrosine metabolism), oxidative phosphorylation (e.g., phosphonate and phosphinate metabolism, and oxidative phosphorylation), energy metabolism (e.g., HIF-1, glycolysis/gluconeogenesis, fatty acid biosynthesis, and fatty acid degradation), and signal transduction (e.g., Wnt, calcium, MAPK, and cGMP-PKG signaling pathways), etc. Further analysis of MAPKs showed that the activation levels of its main members JNK1 and ERK1/2 differed significantly between red and black skin colors. After RNAi was used to interfere with ERK1/2, it was found that the local skin of the tail of P. leopardus would turn black. Combined transcriptome and proteome analysis showed that most DEGs-DAPs in red skin were higher than in black skin (58 were upregulated, 1 was downregulated, and 4 were opposite). These DEGs-DAPs showed that the differences between red and black skin tissues of P. leopardus were related primarily to energy metabolism, signal transduction and cytoskeleton. These findings are not only conducive to understand the skin color regulation mechanism of P. leopardus and other coral reef fish, but also provide an important descriptive to the breeding of color strains. SIGNIFICANCE OF THE STUDY: The skin color of P. leopardus gradually darkens or blackens due to environmental factors such as changes in light intensity and human activities, and this directly affects its ornamental and economic value. In this study, RNAseq and TMT were used to conduct comparative quantitative transcriptomics and proteomics and analyze differences between red and black P. leopardus skin. The results showed that energy metabolism, signal transduction and cytoskeleton were the main metabolic pathways causing their skin color differences. These findings contribute to existing data describing fish skin color, and provide information about protein levels, which are of great significance to a deeper understanding of the skin color regulation mechanism in P. leopardus and other coral reef fishes.

Molecular and functional analysis of the microphthalmia-associated transcription factor (mitf) gene duplicates in red tilapia

In vertebrates, the microphthalmia-associated transcription factor (mitf) is at the hub of the melanin synthesis regulation network. However, little information is known about its molecular characterization, expression, location, or function in skin color differentiation and variation of red tilapia. The full-length cDNA sequences (1977 bp and 1999 bp) of mitfa and mitfb, encoding polypeptides of 491 and 514 amino acids, were effectively identified from red tilapia in this study. The Mitfa and Mitfb sequences of red tilapia clustered first with O. aureus, then with other teleost fish, according to phylogenetic analysis. Mitfa and mitfb mRNA were highly expressed in the brain, dorsal skin and eye tissues using quantitative real-time PCR. The mRNA expressions of mitfa and mitfb were the highest in the cleavage stage during the early development of red tilapia. Among three different colors of red tilapia, the expression levels of mitfa and mitfb were highest in the PB (pink with scattered black spots) dorsal skin. After overwintering, the mitfa and mitfb mRNA expressions were high in the dorsal skin of PB (color changed from pink to black). Mitfa and mitfb were mostly found in the epidermal layer of the dorsal skin, according to in situ hybridization (ISH) analysis. After injecting mitf-dsRNA duplicates along the tail vein of red tilapia, the activity of tyrosinase and the level of melanin in the dorsal skin both decreased significantly. The mRNA expressions of mitfa and its downstream genes (tyrb, tyrp1a and dct) decreased, whereas the mRNA expression of mitfb increased after mitfa-dsRNA injection. The mRNA expressions of mitfb, tyrb, tyrp1a and dct decreased, whereas the mRNA expression of mitfa increased after injecting mitfb-dsRNA. These findings suggest that mitf gene duplicates may play an important role in red tilapia skin color differentiation and variation via the melanogenesis pathway.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

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.

MicroRNAs May Play an Important Role in Sexual Reversal Process of Chinese Soft-Shelled Turtle, Pelodiscus sinensis

The Chinese soft-shelled (Pelodiscus sinensis) turtle exhibits obvious sex dimorphism, which leads to the higher economic and nutritional value of male individuals. Exogenous hormones can cause the transformation from male to female phenotype during gonadal differentiation. However, the molecular mechanism related to the sexual reversal process is unclear. In this study, we compared the difference between the small RNAs of male, female, and pseudo-female turtles by small RNA-seq to understand the sexual reversal process of Chinese soft-shelled turtles. A certain dose of estrogen can cause the transformation of Chinese soft-shelled turtles from male to female, which are called pseudo-female individuals. The result of small RNA-seq has revealed that the characteristics of pseudo-females are very similar to females, but are strikingly different from males. The number of the microRNAs (miRNAs) of male individuals was significantly less than the number of female individuals or pseudo-female individuals, while the expression level of miRNAs of male individuals were significantly higher than the other two types. Furthermore, we found 533 differentially expressed miRNAs, including 173 up-regulated miRNAs and 360 down-regulated miRNAs, in the process of transformation from male to female phenotype. Cluster analysis of the total 602 differential miRNAs among females, males, and pseudo-females showed that miRNAs played a crucial role during the sexual differentiation. Among these differential miRNAs, we found 12 miRNAs related to gonadal development and verified their expression by qPCR. The TR-qPCR results confirmed the differential expression of 6 of the 12 miRNAs: miR-26a-5p, miR-212-5p, miR-202-5p, miR-301a, miR-181b-3p and miR-96-5p were involved in sexual reversal process, which was consistent with the results of omics. Using these six miRNAs and some of their target genes, we constructed a network diagram related to gonadal development. We suggest that these miRNAs may play an important role in the process of effective sex reversal, which would contribute to the breeding of all male strains of Chinese soft-shelled turtles.

Expression profile of miRNAs involved in the hepatoprotective effects of curcumin against oxidative stress in Nile tilapia

Curcumin is a polyphenol with antioxidant activity that has been used to protect the health of fish livers. Our previous studies about comparative transcriptome have shown that curcumin can enhance the Nrf2-Keap1 signaling pathway and induce downstream anti-stress genes to maintain cell viability. However, the possible role of miRNAs in the protective mechanism of curcumin is not understood. In this study, the tilapia hepatocyte H₂O₂ stress model was used, and the miRNA expression profile for four groups (control group, curcumin group, H₂O₂ group, and protection group) were established by high-throughput sequencing. In our results, 278-333 types of Oreochromis niloticus miRNAs, 309-543 types of conserved miRNAs, and 535-746 types of novel miRNAs were identified in different samples. Differentially expressed miRNAs were identified by comparing miRNA expression profiles among the four groups. The expression levels were confirmed by q-PCR. The target genes of these differentially expressed miRNAs were predicted, and their functional annotations were enriched by GO and KEGG analysis, which revealed that many target genes were involved in "response to stimulus" and "antioxidant activity" in each pair of groups. Several miRNAs related to oxidative stress showed differential expression. For example, in the H₂O₂ group, the expression of miR-122 was decreased, and the expression of miR-21 and miR-489 increased significantly. In the curcumin group, the expression of miR-153b was decreased, and the expression of miR-200a and miR-29 was increased significantly. miR-153b, miR-200a, and miR-29 may be involved in the regulation of the Nrf2-Keap1 signaling pathway by curcumin. This work might provide insights into the molecular mechanisms of miRNA regulation of curcumin on the prevention and alleviation of liver diseases in fish.

Integrated Analysis of Coding Genes and Non-coding RNAs Associated with Shell Color in the Pacific Oyster (Crassostrea gigas)

Molluscan shell color polymorphism is important in genetic breeding, while the molecular information mechanism for shell coloring is unclear. Here, high-throughput RNA sequencing was used to compare expression profiles of coding and non-coding RNAs (ncRNAs) from Pacific oyster Crassostrea gigas with orange and black shell, which were from an F2 family constructed by crossing an orange shell male with a black shell female. First, 458, 13, and 8 differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified, respectively. Functional analysis suggested that the DEGs were significantly enriched in 9 pathways including tyrosine metabolism and oxidative phosphorylation pathways. Several genes related to melanin synthesis and biomineralization expressed higher whereas genes associated with carotenoid pigmentation or metabolism expressed lower in orange shell oyster. Then, based on the ncRNA analysis, 163 and 20 genes were targeted by 13 and 8 differentially expressed lncRNAs (DELs) and miRNAs (DEMs), severally. Potential DELs-DEMs-DEGs interactions were also examined. Seven DEMs-DEGs pairs were detected, in which tyrosinase-like protein 1 was targeted by lgi-miR-133-3p and lgi-miR-252a and cytochrome P450 was targeted by dme-miRNA-1-3p. These results revealed that melanin synthesis-related genes and miRNAs-mRNA interactions functioned on orange shell coloration, which shed light on the molecular regulation of shell coloration in marine shellfish.

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.

The stage-specific long non-coding RNAs and mRNAs identification and analysis during early development of common carp, Cyprinus carpio

Cyprinus carpio is considered an alternative vertebrate fish model to zebrafish. However, systemic times-series research on the lncRNAs and mRNAs during early development of C. carpio has not been reported yet. This study provides the first long non-coding RNA (lncRNA)-mRNA expression profiles during six main early development stages (2 h post-fertilization hpf, 6 hpf, 12 hpf, 20 hpf, 64 hpf and 1 day post-hatching). A total of 51,979 lncRNAs were identified. We screened the top 10 abundance lncRNAs and mRNAs and stage-specific lncRNAs and mRNAs (specificity measure SPM > 0.9). We identified significant differentially expressed lncRNAs and mRNAs (|log2 (fold change)| ≥ 1 and false discovery rate FDR of <0.05). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified numerous signaling pathways. Additionally, the lncRNA-mRNA co-regulated network analysis of two lncRNAs (lncrps25 and malat1) and two mRNAs (mitf and troponin T) were investigated. Our results provide new insight into the role of lncRNAs and mRNAs, and would advance the understanding of lncRNA-mediated mechanisms in early development of fish.

MicroRNA expression profile analysis of skin immune response in crucian carp (Carassius auratus) infected by Aeromonas hydrophila

MicroRNAs (miRNAs) are non-coding RNA molecules that regulate gene expression in fish, but its regulatory mechanism of the skin mucosal immune response remains poorly understood. In order to investigate the immunological role of miRNAs, three sRNA libraries (mSC, mST1, mST2) from skin samples of crucian carp (Carassiusauratus) infected with Aeromonas hydrophila at three time points (0, 6 and 12 hpi) were constructed and examined using Illumina Hiseq 2000 platform. All of the identified miRNA, rRNA and tRNA were 69444 (13.39%), 29550 (5.70%) and 10704 (2.06%) in skin, respectively. At 6 and 12 hpi, 829 and 856 miRNAs were differentially expressed, respectively. Among these DEMs, 53 known and 10 novel miRNAs were all significantly differentially expressed during early infection (p < 0.01). GO and KEGG enrichment analyses revealed that 118111 target-genes were primarily involved in cellular process, metabolic process, biological regulation and stress response, such as antigen processing and presentation, complement and coagulation cascades, phagosome, MAPK, TLR, NF-κB and JAK-STAT signaling pathways. These results will help to elucidate the mechanism of miRNAs involved in the skin mucosal immune response of crucian carp against Aeromonas hydrophila infection.

The Mantle Exosome and MicroRNAs of Hyriopsis cumingii Involved in Nacre Color Formation

The nacre color of shells has an effect on the pearl color in Hyriopsis cumingii and is an important indicator for its value. The nacre is part of the shell, and some studies have shown that exosomes of the mantle are involved in the formation of shells. Most of the RNA contained in exosomes are microRNAs (miRNAs); however, little information is available on the roles of exosomes and miRNAs on the formation of nacre color in mussels. In this study, exosomes of mantles were extracted from white and purple mussels. High-throughput Illumina sequencing was performed on the white and purple mussel mantle exosomes, and 7,665,167 and 10,994,115 reads were harvested. Using the standard of |log₂(Fold change)| ≥ 2, and a p value ≤ 0.05, a total of 54 differentially expressed miRNAs were identified. The miRNAs that regulated the target genes (hcApo, HcTyr, HcTyp-1, HcMitf, HcSRCR1, and HcSRCR2) involved in shell color formation were predicted. Moreover, miR-15b negatively regulated hcApo, which plays important roles in the absorption and transport of β-carotene in H. cumingii. These results improve our understanding of the molecular mechanisms of nacre color formation in H. cumingii.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.

Identification and characterization of skin color microRNAs in Koi carp (Cyprinus carpio L.) by Illumina sequencing

Background

MicroRNAs (miRNAs) are endogenous, small (21–25 nucleotide), non-coding RNAs that play important roles in numerous biological processes. Koi carp exhibit diverse color patterns, making it an ideal subject for studying the genetics of pigmentation. However, the influence of miRNAs on skin color regulation and variation in Koi carp is poorly understood.

Results

Herein, we performed small RNA (sRNA) analysis of the three main skin colors in Koi carp by Illumina sequencing. The results revealed 330, 397, and 335 conserved miRNAs (belonging to 81 families) and 340, 353, and 351 candidate miRNAs in black, red, and white libraries, respectively. A total of 164 differentially expressed miRNAs (DEMs) and 14 overlapping DEMs were identified, including miR-196a, miR-125b, miR-202, miR-205-5p, miR-200b, and etc. Target prediction and functional analysis of color-related miRNAs such as miR-200b, miR-206, and miR-196a highlighted putative target genes, including Mitf, Mc1r, Foxd3, and Sox10 that are potentially related to pigmentation. Determination of reference miRNAs for relative quantification showed that let-7a was the most abundant single reference gene, and let-7a and miR-26b was the most abundant combination.

Conclusions

The findings provide novel insight into the molecular mechanisms determining skin color differentiation in Koi carp, and serve as a valuable reference for future studies on tissue-specific miRNA abundance in Koi carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5189-5) contains supplementary material, which is available to authorized users.

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.