Comprehensive analysis of hub mRNA, lncRNA and miRNA, and associated ceRNA networks implicated in cobia (Rachycentron canadum) scales under hypoosmotic adaption

Abstract

Salinity plays a vital role in fish aquaculture, profoundly influencing the growth and development of fish. Scales, as the protective outer layer of fish, function as a critical defense against external factors. In this study, we employed transcriptome sequencing to analyze the ceRNA expression profiles to reveal the effect of salinity acclimation on transcriptional expression changes in the scales of cobia (Rachycentron canadum). The results revealed that after being exposed to a salinity level of 15 ‰ for just one day (1D), a total of 407 mRNAs/genes were significantly regulated; 66 miRNAs were respectively significantly regulated; and 109 target genes of the differentially expressed miRNAs were significantly regulated; a total of 185 differently expressed lncRNAs and 292 differently expressed target genes (DetGenes) of differently expressed lncRNAs were also identified. After 7 days (7D), a total of 2195 mRNAs/genes were found to be significantly regulated and 82 miRNAs were significantly regulated; among the target genes of the differentially expressed miRNAs, 245 were regulated. Moreover, 438 differently expressed lncRNAs and 681 DetGenes of these lncRNAs were identified. Subsequent analysis through GO, KEGG pathway, in 1D vs. CG (control group), DeGenes, which first respond to changes in salinity, are mainly involved in negative regulation of macrophage differentiation, negative regulation of granulocyte differentiation and negative regulation of phagocytosis, and are mainly related to biological processes related to the immune function of fish. After a 7-day process, DeGenes were enriched in the collagen fibril organization, regulation of nodal signaling pathway and cell recognition biology processes. These biological processes are not only related to the immune function of fish, but more importantly, to the physiological structure of fish. By analyzing the co down-regulated miRNAs of 1D vs. CG, as well as 7D vs. CG, the functions of these miRNAs are mainly related to bone differentiation and development. In addition，ceRNA network uncovered that the effect of salinity is temporal. The first competing lncRNAs mainly regulated genes related to physiological processes and biological development, while target genes related to immunity and body defense were less competitive. On the contrary, after a period of salinity treatment, the types of competing lncRNAs involved changed.