Insights into molecular mechanisms of skin color variation in hybrid F1 of wild-type rainbow trout ♂ × yellow mutant rainbow trout ♀ through whole transcriptome and ceRNA network analysis

Rainbow trout (Oncorhynchus mykiss) is a well-known economical fish globally, and its market value and acceptance are greatly influenced by skin color. However, little is known about the roles of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in rainbow trout skin color variation. In this study, we compared the expression profiles of whole transcriptome in dorsal skin from wild-type (WTrt) and yellow mutant rainbow trout (YMrt) that were F1 generation of WTrt ♂ × YMrt ♀. A total of 115 differentially expressed (DE) mRNAs, 126 lncRNAs, 18 circRNAs and 56 miRNAs were identified, and several core genes (gch1, pts, slc2all, pax7a/pax7b, dgat2, plin6 and pnpla2) involved in pteridine and carotenoid synthesis were significantly upregulated in YMrt. Interestingly, we found that melanogenesis-related genes tyr, dct, mitfa and pmel were also notably elevated in YMrt. Functional enrichment analysis revealed some pigmentation-related GO terms and KEGG pathways, including pteridine-containing compound biosynthetic process, pigment metabolic process, melanosome, melanosome transport, tyrosine metabolism and melanogenesis. Based on competing endogenous RNA analysis, novel_circ_000056 was found to regulate dgat2 by targeting novel-m0344-3p. This study provides a systematic perspective on the potential functions of these DEmRNAs and ncRNAs, and lays a foundation for further study on the molecular mechanisms governing rainbow trout skin color variation.