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

cAMP-Mediated CREM-MITF-TYR Axis Regulates Melanin Synthesis in Pacific Oysters

Colorful shells in bivalves are mostly caused by the presence of biological pigments, among which melanin is a key component in the formation of shell colours. Cyclic adenosine monophosphate (cAMP) is an important messenger in the regulation of pigmentation in some species. However, the role of cAMP in bivalve melanogenesis has not yet been reported. In this study, we performed in vitro and in vivo experiments to determine the role of cAMP in regulating melanogenesis in Pacific oysters. Besides, the function of cAMP-responsive element modulator (CREM) and the interactions between CREM and melanogenic genes were investigated. Our results showed that a high level of cAMP promotes the expression of melanogenic genes in Pacific oysters. CREM controls the expression of the MITF gene under cAMP regulation. In addition, CREM can regulate melanogenic gene expression, tyrosine metabolism, and melanin synthesis. These results indicate that cAMP plays an important role in the regulation of melanogenesis in Pacific oysters. CREM is a key transcription factor in the oyster melanin synthesis pathway, which plays a crucial role in oyster melanin synthesis through a cAMP-mediated CREM-MITF-TYR axis.

Transcriptomic and Physiological Analysis Reveal Melanin Synthesis-Related Genes and Pathways in Pacific Oysters (Crassostrea gigas)

Shell color is one of the shell traits of molluscs, which has been regarded as an economic trait in some bivalves. Pacific oysters (Crassostrea gigas) are important aquaculture shellfish worldwide. In the past decade, several shell color strains of C. gigas were developed through selective breeding, which provides valuable materials for research on the inheritance pattern and regulation mechanisms of shell color. The inheritance patterns of different shell colors in C. gigas have been identified in certain research; however, the regulation mechanism of oyster pigmentation and shell color formation remains unclear. In this study, we performed transcriptomic and physiological analyses using black and white shell oysters to investigate the molecular mechanism of melanin synthesis in C. gigas. Several pigmentation-related pathways, such as cytochrome P450, melanogenesis, tyrosine metabolism, and the cAMP signaling pathway were found. The majority of differentially expressed genes and some signaling molecules from these pathways exhibited a higher level in the black shell oysters than in the white, especially after L-tyrosine feeding, suggesting that those differences may cause a variation of tyrosine metabolism and melanin synthesis. In addition, the in vitro assay using primary cells from mantle tissue showed that L-tyrosine incubation increased cAMP level, gene and protein expression, and melanin content. This study reveals the difference in tyrosine metabolism and melanin synthesis in black and white shell oysters and provides evidence for the potential regulatory mechanism of shell color in oysters.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Integrated analysis of mRNAs and lncRNAs reveals candidate marker genes and potential hub lncRNAs associated with growth regulation of the Pacific Oyster, Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is an economically important shellfish around the world. Great efforts have been made to improve its growth rate through genetic breeding. However, the candidate marker genes, pathways, and potential lncRNAs involved in oyster growth regulation remain largely unknown. To identify genes, lncRNAs, and pathways involved in growth regulation, C. gigas spat was cultured at a low temperature (15 ℃) to yield a growth-inhibited model, which was used to conduct comparative transcriptome analysis with spat cultured at normal temperature (25 ℃).

RESULTS

In total, 8627 differentially expressed genes (DEGs) and 1072 differentially expressed lncRNAs (DELs) were identified between the normal-growth oysters (cultured at 25 ℃, hereinafter referred to as NG) and slow-growth oysters (cultured at 15 ℃, hereinafter referred to as SG). Functional enrichment analysis showed that these DEGs were mostly enriched in the AMPK signaling pathway, MAPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. LncRNAs analysis identified 265 cis-acting pairs and 618 trans-acting pairs that might participate in oyster growth regulation. The expression levels of LNC_001270, LNC_003322, LNC_011563, LNC_006260, and LNC_012905 were inducible to the culture temperature and food abundance. These lncRNAs were located at the antisense, upstream, or downstream of the SREBP1/p62, CDC42, CaM, FAS, and PIK3CA genes, respectively. Furthermore, the expression of the trans-acting lncRNAs, including XR_9000022.2, LNC_008019, LNC_015817, LNC_000838, LNC_00839, LNC_011859, LNC_007294, LNC_006429, XR_002198885.1, and XR_902224.2 was also significantly associated with the expression of genes enriched in AMPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process.

CONCLUSIONS

In this study, we identified the critical growth-related genes and lncRNAs that could be utilized as candidate markers to illustrate the molecular mechanisms underlying the growth regulation of Pacific oysters.

Transcriptome Analysis Revealed That Hydrogen Peroxide-Regulated Oxidative Phosphorylation Plays an Important Role in the Formation of Pleurotus ostreatus Cap Color

Pleurotus ostreatus is widely cultivated in China. H2O2, as a signaling molecule, can regulate the formation of cap color, but its regulatory pathway is still unclear, severely inhibiting the breeding of dark-colored strains. In this study, 614 DEGs specifically regulated by H2O2 were identified by RNA-seq analysis. GO-enrichment analysis shows that DEGs can be significantly enriched in multiple pathways related to ATP synthesis, mainly including proton-transporting ATP synthesis complex, coupling factor F(o), ATP biosynthetic process, nucleoside triphosphate metabolic processes, ATP metabolic process, purine nucleoside triphosphate biosynthetic and metabolic processes, and purine ribonuclease triphosphate biosynthetic metabolic processes. Further KEGG analysis revealed that 23 DEGs were involved in cap color formation through the oxidative phosphorylation pathway. They were enriched in Complexes I, III, IV, and V in the respiratory chain. Further addition of exogenous uncoupling agents and ATP synthase inhibitors clarifies the important role of ATP synthesis in color formation. In summary, H2O2 may upregulate the expression of complex-encoding genes in the respiratory chain and promote ATP synthesis, thereby affecting the formation of cap color. The results of this study lay the foundation for the breeding of dark-colored strains of P. ostreatus and provide a basis for the color-formation mechanism of edible fungi.

Effects of Dietary Cystine and Tyrosine Supplementation on Melanin Synthesis in the Pacific Oyster (Crassostrea gigas)

Melanogenesis is a multistep process to produce melanin for dark pigmentation in skin coloration. Previous studies in vertebrates demonstrated that cystine and tyrosine amino acids are involved in the melanin synthesis. However, very little is known about the melanogenesis in bivalve. In this study, cystine supplementation for 30 days significantly upregulated the expression of CgB-aat1, CgCbs and CgTyr and pheomelanin content in the Pacific oyster Crassostrea gigas. Transmission electron microscope (TEM) results revealed more melanosomes in the connective tissue and melanin granules were secreted in epithelium of mantle. In contrast, tyrosine supplementation had no clear effect on melanogenesis except the gene expression changes of CgB-aat1 and CgCbs. In addition, prolonged supplementation of cystine or tyrosine for 60 days had a negative impact on melanogenesis. Indeed, after 60 days, expression of most of the melanin synthesis-related genes under study was decreased, and melanin content was significantly reduced, indicating that cystine and tyrosine might inhibit production of eumelanin and pheomelanin, respectively. In addition, in vitro analysis using primary cell culture from mantle tissue indicated that incubation with cystine, tyrosine, or B-AAT1 polypeptide, CBS/TYR recombinant proteins induced the increase of CgB-aat1 and CgCbs expression in a dose-dependent manner, suggesting the presence of a regulatory network in response to cystine and tyrosine amino acids intakes in pheomelanin synthesis-related gene expression. Taken together, these data indicate that cystine-CgB-aat1-CgCbs-CgTyr axis is a potential regulator of the pheomelanin biosynthesis pathway, and thus plays an important role in the mantle pigmentation in C. gigas. This work provides a new clue for selective cultivation of oyster strains with specific shell colors in bivalve breeding.

Transcription factor CgPOU3F4-like regulates expression of pheomelanin synthesis related gene CgB-aat1 in the Pacific oyster (Crassostrea gigas)

Previous study has found that b (0, +) -type amino acid transporter 1 (CgB-aat1) plays an essential role on mantle pigmentation in the Pacific oyster Crassostrea gigas. However, the molecular regulation of CgB-aat1 gene expression remains unclear. Herein, three POU domain family members, CgPOU2F1, CgPOU3F4-like and CgPOU4F3-X1 were characterized and they all had POUs and HOX domains, respectively, which were important in transcriptional regulation. CgPOU3F4-like gene expression was the highest in mantle edge. Subsequently, the dual-luciferase reporter result showed that the core regulatory region of CgB-aat1 gene was from -632 to -350 bp of promoter. In transient co-transfection assays, the strongest activity was activated only by CgPOU3F4-like, suggesting CgPOU3F4-like was a valid transcriptional activator of CgB-aat1 gene promoter. And the structural integrity of CgPOU3F4-like was essential for its activation function. In addition, site directed mutagenesis assay was applied to detect three key binding sites between CgPOU3F4-like and core region of CgB-aat1 gene promoter, and this interaction was verified by ChIP test. Furthermore, CgPOU3F4-like knockdown by RNA interference led to obvious decreases in CgB-aat1 and cystathionine beta-synthase (CgCbs) expressions at both mRNA and protein levels. Collectively, these results indicate that CgPOU3F4-like positively regulate CgB-aat1 gene expression and it may be a critical upstream transcriptional regulation factor in pheomelanin synthesis in C. gigas.

Histological, elemental, and ultrastructural analysis of melanin in mantle of Pacific oyster (Crassostrea gigas)

Colorful shell of bivalve is mainly because of the biological pigments, of which melanin plays an important role in shell color formation. More and more studies focus on the genes function involved in melanin synthesis, but relatively few studies address the biochemical character and ultrastructure of melanin in bivalve from microscopic perspective. Here, we investigated the histological structure of mantle of Crassostrea gigas with orange shell color. Distribution of melanin in mantle was verified with histochemical staining. In addition, immunofluorescence technique showed that strongly positive signal of CgTYR was specific to the mantle margin, which is consistence with the location of brown granules in H&E staining. The further result of elementary composition of melanin displayed that metal Ca, Fe, and Zn were detected using scanning transmission electron microscope and energy dispersive spectroscopy mapping methods. Next, based on TEM observations, it was speculated that the series of cellular events leading to the formation and release of melanin. Melanocyte in the primary stage showed many mitochondria and rough endoplasmic reticulum as well as an extensive Golgi complex with numerous vesicles intermingled with melanosome. Subsequently, melanosome was expended and their hue gradually intensified, and Golgi complex and mitochondria were still observed in the cytoplasm. Finally, after melanosome was discharged into intercellular spaces, the disintegration of membranes in some cells, and severe cellular vacuolization. These data enrich the understanding of ultrastructural characteristic and formation of melanin in mantle of bivalve and pave the way for further investigating shell coloration at the cellular level.

Involvement of B-aat1 and Cbs in regulating mantle pigmentation in the Pacific oyster (Crassostrea gigas)

BACKGROUND

Shell color formation is an important physiological process in bivalves, the molecular genetic basis has potential application in bivalve aquaculture, but there is still remaining unclear about this issue. The cystine/glutamate transporter (Slc7a11) and cystathionine beta-synthase (Cbs) are integral genes in pheomelanin synthesis pathway, which is vital to skin pigmentation.

METHODS AND RESULTS

Here, the sequences of b (0, +) -type amino acid transporter 1 (B-aat1) and Cbs in Pacific oyster (Crassostrea gigas) (CgB-aat1, CgCbs) were characterized. Phylogenetically, the deduced amino acid sequences of CgB-aat1 and CgCbs both possessed conserved features. Genes were both ubiquitously expressed in six tested tissues with more abundant expression level in central mantle. Besides, the polyclonal antibodies of CgB-aat1, CgCbs, CgTyr, and CgTyrp2 were successfully prepared. Immunofluorescence analysis revealed that CgB-aat1 and CgCbs proteins were both expressed in gill rudiments of eyed-larvae and concentrated mainly in cytoplasm of epithelial cell and nerve axons in mantle. Additionally, after CgB-aat1 or CgCbs silencing, expressions at mRNA and protein levels of CgB-aat1 and CgCbs involved in pheomelanin synthesis were significantly suppressed, and CgTyr, CgTyrp1 and CgTyrp2 related to eumelanin synthesis were also down-regulated but no apparent differences, respectively. Moreover, micrographic examination found less brown-granules at mantle edge in CgB-aat1 interference group.

CONCLUSION

These results implied that pheomelanin synthesis was possible induced by CgB-aat1-CgTyr-CgCbs axis, and it played an essential role on mantle pigmentation in the oysters. These findings provide the useful genetic knowledge and enrich the physiological information for the shell color formation in bivalve aquaculture.

Divergent Carry-Over Effects of Hypoxia during the Early Development of Abalone

After being exposed to environmental stimuli during early developmental stages, some organisms may gain or weaken physiological regulating abilities, which would have long-lasting effects on their performance. Environmental hypoxia events can have significant effects on marine organisms, but for breeding programs and other practical applications, it is important to further explore the long-term physiological effects of early hypoxia exposure in economically significant species. In this study, the Pacific abalone Haliotis discus hannai was exposed to moderate hypoxia (∼4 mg/L) from zygote to trochophora, and the assessments of hypoxia tolerance were conducted on the grow-out stage. The results revealed that juvenile abalones exposed to hypoxia at the early development stages were more hypoxia-tolerant but with slower weight growth, a phenomenon called the trade-off between growth and survival. These phenotypic effects driven by the hypoxia exposure were explained by strong selection of genes involved in signal transduction, autophagy, apoptosis, and hormone regulation. Moreover, long non-coding RNA regulation plays an important role modulating carry-over effects by controlling DNA replication and repair, signal transduction, myocardial activity, and hormone regulation. This study revealed that the ability to create favorable phenotypic differentiation through genetic selection and/or epigenetic regulation is important for the survival and development of aquatic animals in the face of rapidly changing environmental conditions.

Genome-Wide Identification of Long Non-coding RNAs in Crassostrea gigas and Their Association with Heat Stress

Oysters face a complex and changeable environment in the intertidal zone. Heat stress is the main cause of their mass summer deaths. Several important genes are identified to be associated with heat response for oysters. However, regulation of these heat response genes in oysters remains largely unknown. In this study, 27 RNA-seq datasets are used to give a relatively comprehensive of long noncoding RNA (lncRNA) sets for C. gigas. Then, the differential expressed genes and lncRNAs are identified under heat stress. Among all the heat shock proteins (HSPs) and inhibitor of apoptosis proteins (IAPs) in the C. gigas genome, 25 heat shock proteins and 14 IAPs are differential expressed. The Gene Ontology analysis reveals that differential expressed genes (DEGs) are enriched in 6, 7, and 7 GO terms in cellular components, molecular function, and biological process, respectively. Within these terms, cellular response to stimulus is the most abundant term. Furthermore, the potential cis target of the differential expressed lncRNAs (DELs) are predicted to investigate their functions. Of the 394 DELs, there are 80 DELs being found to be corresponded to 113 protein-coding genes. Of them, eight HSPs are found to be regulated by their lncRNA regulators under heat stress. This work provides valuable resource of lncRNA and their regulatory roles under heat stress in C. gigas and gives new insights into adaptive evolution in marine mollusks.

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.