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

Melanism: Cryptic control by non-coding RNAs

Melanism drives both crypsis and mimicry in butterflies and moths. To date, melanism has been mapped to a structural gene called cortex, but now more detailed work shows that in fact it is controlled by non-coding RNAs at the same locus.

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.

Protoporphyrin IX metabolism mediated via translocator protein (CgTspO) involved in orange shell coloration of pacific oyster (Crassostrea gigas)

Mollusc shell color polymorphism is influenced by various factors. Pigments secreted in vivo by animals play a critical role in shell coloration. Among the different shell-color hues, orange pigmentation has been partially attributed to porphyrins. However, the detailed causal relationship between porphyrins and orange-shell phenotype in molluscs remains largely unexplored. The various strains of Pacific oyster (Crassostrea gigas) with different shell color provide useful models to study the molecular regulation of mollusc coloration. Accordingly, oysters with orange and gold-shells, exhibiting distinct porphyrin distributions, were selected for analysis of total metabolites and gene expression profile through mantle metabolomic and transcriptomic studies. Translocator protein (TspO) and protoporphyrin IX (PPIX) were identified as potential factors influencing oyster shell-color. The concentration of PPIX was measured using HPLC, while expression profiling of CgTspO was analyzed by qPCR, in situ hybridization, Western blotting, and immunofluorescence techniques. Moreover, the roles of CgTspO in regulating PPIX metabolism and affecting the orange-shell-coloration were investigated in vitro and in vivo. These studies indicate that PPIX and its associated metabolic protein, CgTspO may serve as new regulators of orange-shell-coloration in C. gigas. Data of this study offer new insights into oyster shell coloration and enhancing understandings of mollusc shell color polymorphism.

ARGONAUTE2 Localizes to Sites of Sporocysts in the Schistosome-Infected Snail, Biomphalaria glabrata

MicroRNAs (miRNAs) are a class of small regulatory RNA that are generated via core protein machinery. The miRNAs direct gene-silencing mechanisms to mediate an essential role in gene expression regulation. In mollusks, miRNAs have been demonstrated to be required to regulate gene expression in various biological processes, including normal development, immune responses, reproduction, and stress adaptation. In this study, we aimed to establishment the requirement of the miRNA pathway as part of the molecular response of exposure of Biomphalaria glabrata (snail host) to Schistosoma mansoni (trematode parasite). Initially, the core pieces of miRNA pathway protein machinery, i.e., Drosha, DGCR8, Exportin-5, Ran, and Dicer, together with the central RNA-induced silencing complex (RISC) effector protein Argonaute2 (Ago2) were elucidated from the B. glabrata genome. Following exposure of B. glabrata to S. mansoni miracidia, we identified significant expression up-regulation of all identified pieces of miRNA pathway protein machinery, except for Exportin-5, at 16 h post exposure. For Ago2, we went on to show that the Bgl-Ago2 protein was localized to regions surrounding the sporocysts in the digestive gland of infected snails 20 days post parasite exposure. In addition to documenting elevated miRNA pathway protein machinery expression at the early post-exposure time point, a total of 13 known B. glabrata miRNAs were significantly differentially expressed. Of these thirteen B. glabrata miRNAs responsive to S. mansoni miracidia exposure, five were significantly reduced in their abundance, and correspondingly, these five miRNAs were determined to putatively target six genes with significantly elevated expression and that have been previously associated with immune responses in other animal species, including humans. In conclusion, this study demonstrates the central importance of a functional miRNA pathway in snails, which potentially forms a critical component of the immune response of snails to parasite exposure. Further, the data reported in this study provide additional evidence of the complexity of the molecular response of B. glabrata to S. mansoni infection: a molecular response that could be targeted in the future to overcome parasite infection and, in turn, human schistosomiasis.

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.

Comprehensive analysis of differentially expressed mRNA, lncRNA and miRNA, and their ceRNA networks in the regulation of shell color in the Manila clam (Ruditapes philippinarum)

The regulatory mechanism of ceRNA network plays an important role in molecular function and biological processes, however, the molecular mechanism in the shell color of Ruditapes philippinarum has not yet been reported. In this study, we performed transcriptome sequencing on the mantle of R. philippinarum with different shell colors, and screened for mRNA, miRNA, and lncRNA. A total of 61 mRNAs, 3725 lncRNAs and 90 miRNAs were obtained from all the shell color comparison groups (all mRNAs, lncRNAs and miRNAs P < 0.05), and 7 mRNAs, 8 lncRNAs, and 4 miRNAs of the porphyrin pathway and melanin pathway were screened for competitive endogenous RNA (ceRNA) network construction. The results indicate that the ceRNA network composed of mRNA and lncRNA, centered around efu-miR-101, mle-bantam-3p, egr-miR-9-5p, and sma-miR-75p, may play a crucial regulatory role in shell color formation. This study reveals for the first time the mechanism of ceRNA regulatory networks in the shell color of R. philippinarum and providing important reference data for molecular breeding of shell color 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.

The mechanisms and factors that induce trained immunity in arthropods and mollusks

Besides dividing the organism's immune system into adaptive and innate immunity, it has long been thought that only adaptive immunity can establish immune memory. However, many studies have shown that innate immunity can also build immunological memory through epigenetic reprogramming and modifications to resist pathogens' reinfection, known as trained immunity. This paper reviews the role of mitochondrial metabolism and epigenetic modifications and describes the molecular foundation in the trained immunity of arthropods and mollusks. Mitochondrial metabolism and epigenetic modifications complement each other and play a key role in trained immunity.

Heme-Peroxidase 2 Modulated by POU2F1 and SOX5 is Involved in Pigmentation in Pacific Oyster (Crassostrea gigas)

Color polymorphism is frequently observed in molluscan shellfish, while the molecular regulation of shell pigmentation is not well understood. Peroxidase is a key enzyme involved in melanogenesis. Here, we identified a heme-peroxidase 2 gene (CgHPX2), and characterized the expression patterns and transcriptional regulation of CgHPX2 in the Pacific oyster Crassostrea gigas. Tissues expression analysis showed that CgHPX2 was a mantle-specific gene and primarily expressed in the edge mantle in black shell color oyster compared with white shell oyster. In situ hybridization showed that strong signals for CgHPX2 were detected in the both inner and outer surface of the outer fold of mantle in the black shell color oyster, whereas positive signals in white shell oyster were mainly localized in the outer surface of the outer fold of mantle. In the embryos and larvae, a high expression level of CgHPX2 was detected in the trochophore stage in both black and white shell color oysters. The temporal localization of CgHPX2 was mainly detected in the shell gland and edge mantle of trochophore and calcified shell larvae, respectively. In addition, a 2227 bp of 5' flanking region sequence of CgHPX2 was cloned, which contained a presumed core promoter region and many potential transcription factor binding sites. Further luciferase assay experiment confirmed that POU domain, class 2, transcription factor 1 (POU2F1), and SRY-box transcription factor 5 (SOX5) were involved in transcriptional regulation of CgHPX2 gene through binding to its specific promoter region. After CgPOU2F1 and CgSOX5 RNA interference, the CgHPX2 gene expression was significantly decreased. These results suggested that CgPOU2F1 and CgSOX5 might be two important transcription factors that positively regulated the expression of CgHPX2 gene, improving our understanding of the transcriptional regulation of molluscan shell pigmentation.

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.

Role of miRNAs in melanin metabolism: Implications in melanin-related diseases

BACKGROUND

MicroRNAs (miRNAs) are short single-stranded non-coding RNAs that regulate degradation and expression of messenger RNA (mRNA) and play a wide range of key roles in different biological processes. They mediate different stages of melanocyte differentiation, growth, and apoptosis through a variety of pathways and can mediate melanin production by targeting key enzymes.

AIMS

This article was aimed to review the role of miRNAs in melanin metabolism and to introduce the role and significance of miRNAs in melanin-related diseases.

MATERIALS & METHODS

Systematic search and retrospective review were performed on the published data.

RESULTS

This paper reviews the process of melanin synthesis and the regulatory mechanism, explores the miRNA expression profiles in different model organisms, and introduces the mechanisms of several key miRNAs participating in melanin metabolism through target genes. We also explore the potential role of miRNA as a new target for the treatment of melanin metabolism disease, including vitiligo, melanoma, and chloasma.

CONCLUSION

miRNAs play a key role in melanin-related diseases, and the miRNAs involved may be potential therapeutic targets.

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.

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.

Environmentally Driven Color Variation in the Pearl Oyster Pinctada margaritifera var. cumingii (Linnaeus, 1758) Is Associated With Differential Methylation of CpGs in Pigment- and Biomineralization-Related Genes

Today, it is common knowledge that environmental factors can change the color of many animals. Studies have shown that the molecular mechanisms underlying such modifications could involve epigenetic factors. Since 2013, the pearl oyster Pinctada margaritifera var. cumingii has become a biological model for questions on color expression and variation in Mollusca. A previous study reported color plasticity in response to water depth variation, specifically a general darkening of the nacre color at greater depth. However, the molecular mechanisms behind this plasticity are still unknown. In this paper, we investigate the possible implication of epigenetic factors controlling shell color variation through a depth variation experiment associated with a DNA methylation study performed at the whole genome level with a constant genetic background. Our results revealed six genes presenting differentially methylated CpGs in response to the environmental change, among which four are linked to pigmentation processes or regulations (GART, ABCC1, MAPKAP1, and GRL101), especially those leading to darker phenotypes. Interestingly, the genes perlucin and MGAT1, both involved in the biomineralization process (deposition of aragonite and calcite crystals), also showed differential methylation, suggesting that a possible difference in the physical/spatial organization of the crystals could cause darkening (iridescence or transparency modification of the biomineral). These findings are of great interest for the pearl production industry, since wholly black pearls and their opposite, the palest pearls, command a higher value on several markets. They also open the route of epigenetic improvement as a new means for pearl production improvement.