Low temperature stress represents a major threat to the lives of both farmed and wild fish species. However, biological pathways determining the development of cold resistance in fish remain largely unknown. Zebrafish larvae at 96 hpf were exposed to lethal cold stress (10 °C) for different time periods to evaluate the adverse effects at organism, tissue and cell levels. Time series RNA sequencing (RNA-seq) experiments were performed to delineate the transcriptomic landscape of zebrafish larvae under cold stress and during the subsequent rewarming phase. The genes regulated by cold stress were characterized by progressively enhanced or decreased expression, whereas the genes associated with rewarming were characterized by rapid upregulation upon return to normal temperature (28 °C). Genes such as trib3, dusp5 and otud1 were identified as the representative molecular markers of cold-induced damages through network analysis. Biological pathways involved in cold stress responses were mined from the transcriptomic data and their functions in regulating cold resistance were validated using specific inhibitors. The autophagy, FoxO and MAPK (mitogen-activated protein kinase) signaling pathways were revealed to be survival pathways for enhancing cold resistance, while apoptosis and necroptosis were the death pathways responsible for cold-induced mortality. Functional mechanisms of the survival-enhancing factors Foxo1, ERK (extracellular signal-regulated kinase) and p38 MAPK were further characterized by inhibiting their activities upon cold stress and analyzing gene expression though RNA-seq. These factors were demonstrated to determine the cold resistance of zebrafish through regulating apoptosis and p53 signaling pathway. These findings have provided novel insights into the stress responses elicited by lethal cold and shed new light on the molecular mechanisms underlying cold resistance of fish.

Properly regulated transcriptional responses to environmental perturbations are critical for the fitness of fish. Although gene expression profiles in the tissues of common carp upon cold stress were previously characterized, the transcriptional programs underlying cold acclimation are still not well known. In this study, the ability of three common carp strains including Hebao red carp (HB), Songpu mirror carp (SPM) and Yellow river carp (YR) to establish cold resistance after acclimation to a mild hypothermia stress at 18°C for 24 h was confirmed by measurements of the critical thermal minimums (CTMin). The gene expression profiles of the brain and the heart from these strains under both control and cold-acclimated conditions were characterized with RNA-sequencing. The data of the three common carp strains with different genetic background were combined in the differential gene expression analyses to balance the effects of genetic diversity on gene expression. Marked effects of tissue origins on the cold-induced transcriptional responses were revealed by comparing the differentially expressed gene (DEG) lists of the two tissues. Functional categories including spliceosome and RNA splicing were highly enriched in the DEGs of both tissues. However, steroid biosynthesis was specifically enriched in DEGs of the brain and response to unfolded protein was solely enriched in DEGs of the heart. Consistent with the up-regulation of the genes involved in cholesterol biosynthesis, total cholesterol content of the brain was significantly increased upon cold stress. Moreover, cold-induced alternative splicing (AS) events were explored and AS of the rbmx (RNA-binding motif protein, X chromosome) gene was confirmed by real-time quantitative PCR. Finally, a core set of cold responsive genes (CRGs) were defined by comparative transcriptomic analyses. Our data provide insights into the transcriptional programs underlying cold acclimation of common carp and offer valuable clues for further investigating the genetic determinants for cold resistance of farmed fish.