Transcriptomic profiling revealed the regulatory pathways and key genes associated with cold tolerance in two eel gobies

Closely related species of the eel goby family (Gobiidae) have evolved divergent resistance to low temperatures, but the molecular mechanisms remain poorly understood. This study used a comparative transcriptomic approach to identify key pathways and genes associated with cold tolerance in two eel goby species. Expression profiles of the cold-tolerant O. lacepedii and the cold-sensitive O. rebecca in control (23 °C) and cold stress groups (15 °C and 11 °C) were analyzed. Differentially expressed genes closely linked to interspecific cold tolerance divergence were identified through transcriptome profiling and Venn diagram analysis. GO and KEGG enrichment analyses revealed that processes related to cellular homeostasis, the PPAR signaling pathway, cellular respiration, and oxidative phosphorylation were activated during the cold tolerance response of eel gobies. WGCNA analysis indicated that the hub genes related to thermogenesis and microtubular stability, specifically PPARGC1A and α-tubulin, may contribute to the high cold tolerance in O. lacepedii. These findings provide key clues for dissection of the molecular mechanisms behind the formation of cold tolerance in eel gobies.