Background: Temperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq. Results: Pre-exposure of 96 hpf zebrafish larvae to cold stress (16°C) for 24 h significantly increased their survival rates under severe cold stress (12°C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains. Conclusions: These findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

Temperature influences nearly all biochemical, physiological and life history activities of fish, but the molecular mechanisms underlying the temperature acclimation remains largely unknown. Previous studies have identified many temperature-regulated genes in adult tissues; however, the transcriptional responses of fish larvae to temperature stress are not well understood. In this study, we characterized the transcriptional responses in larval zebrafish exposed to cold or heat stress using microarray analysis. In comparison with genes expressed in the control at 28°C, a total of 2680 genes were found to be affected in 96 hpf larvae exposed to cold (16°C) or heat (34°C) for 2 and 48 h and most of these genes were expressed in a temperature-specific and temporally regulated manner. Bioinformatic analysis identified multiple temperature-regulated biological processes and pathways. Biological processes overrepresented among the earliest genes induced by temperature stress include regulation of transcription, nucleosome assembly, chromatin organization and protein folding. However, processes such as RNA processing, cellular metal ion homeostasis and protein transport and were enriched in genes up-regulated under cold exposure for 48 h. Pathways such as mTOR signalling, p53 signalling and circadian rhythm were enriched among cold-induced genes, while adipocytokine signalling, protein export and arginine and praline metabolism were enriched among heat-induced genes. Although most of these biological processes and pathways were specifically regulated by cold or heat, common responses to both cold and heat stresses were also found. Thus, these findings provide new interesting clues for elucidation of mechanisms underlying the temperature acclimation in fish.

Hypoxia and temperature stress are two major adverse environmental conditions often encountered by fishes. The interaction between hypoxia and temperature stresses has been well documented and oxygen is considered to be the limiting factor for the thermal tolerance of fish. Although both high and low temperature stresses can impair the cardiovascular function and the cross-resistance between hypoxia and heat stress has been found, it is not clear whether hypoxia acclimation can protect fish from cold injury. Pre-acclimation of 96-hpf zebrafish larvae to mild hypoxia (5% O2) significantly improved their resistance to lethal hypoxia (2.5% O2) and increased the survival rate of zebrafish larvae after lethal cold (10°C) exposure. However, pre-acclimation of 96-hpf larvae to cold (18°C) decreased their tolerance to lethal hypoxia although their ability to endure lethal cold increased. RNA-seq analysis identified 132 up-regulated and 41 down-regulated genes upon mild hypoxia exposure. Gene ontology enrichment analyses revealed that genes up-regulated by hypoxia are primarily involved in oxygen transport, oxidation-reduction process, hemoglobin biosynthetic process, erythrocyte development and cellular iron ion homeostasis. Hypoxia-inhibited genes are enriched in inorganic anion transport, sodium ion transport, very long-chain fatty acid biosynthetic process and cytidine deamination. A comparison with the dataset of cold-regulated gene expression identified 23 genes co-induced by hypoxia and cold and these genes are mainly associated with oxidation-reduction process, oxygen transport, hemopoiesis, hemoglobin biosynthetic process and cellular iron ion homeostasis. The alleviation of lipid peroxidation damage by both cold- and hypoxia-acclimation upon lethal cold stress suggests the association of these genes with cold resistance. Furthermore, the alternative promoter of hmbsb gene specifically activated by hypoxia and cold was identified and confirmed. Acclimation responses to mild hypoxia and cold stress were found in zebrafish larvae and pre-acclimation to hypoxia significantly improved the tolerance of larvae to lethal cold stress. RNA-seq and bioinformatics analyses revealed the biological processes associated with hypoxia acclimation. Transcriptional events co-induced by hypoxia and cold may represent the molecular basis underlying the protection of hypoxia-acclimation against cold injury.

Fish skin serves as the first line of defense against a wide variety of chemical, physical and biological stressors. Secretion of mucus is among the most prominent characteristics of fish skin and numerous innate immune factors have been identified in the epidermal mucus. However, molecular mechanisms underlying the mucus secretion and immune activities of fish skin remain largely unclear due to the lack of genomic and transcriptomic data for most economically important fish species. In this study, we characterized the skin transcriptome of mud loach using Illumia paired-end sequencing. A total of 40364 unigenes were assembled from 86.6 million (3.07 gigabases) filtered reads. The mean length, N50 size and maximum length of assembled transcripts were 387, 611 and 8670 bp, respectively. A total of 17336 (43.76%) unigenes were annotated by blast searches against the NCBI non-redundant protein database. Gene ontology mapping assigned a total of 108513 GO terms to 15369 (38.08%) unigenes. KEGG orthology mapping annotated 9337 (23.23%) unigenes. Among the identified KO categories, immune system is the largest category that contains various components of multiple immune pathways such as chemokine signaling, leukocyte transendothelial migration and T cell receptor signaling, suggesting the complexity of immune mechanisms in fish skin. As for mucin biosynthesis, 37 unigenes were mapped to 7 enzymes of the mucin type O-glycan biosynthesis pathway and 8 members of the polypeptide N-acetylgalactosaminyltransferase family were identified. Additionally, 38 unigenes were mapped to 23 factors of the SNARE interactions in vesicular transport pathway, indicating that the activity of this pathway is required for the processes of epidermal mucus storage and release. Moreover, 1754 simple sequence repeats (SSRs) were detected in 1564 unigenes and dinucleotide repeats represented the most abundant type. These findings have laid the foundation for further understanding the secretary processes and immune functions of loach skin mucus.

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disorder caused by abnormal lipid metabolisms, such as reduced hepatic fatty acid oxidation (FAO), but intracellular control of FAO under physio- and pathological conditions remains largely undefined. Here, we demonstrate that deprivation of Slc7a3a leads to hepatic steatosis in fasted zebrafish as a result of defects in arginine-dependent nitric oxide (NO) synthesis. Fast-induced hepatic steatosis in slc7a3a-null mutants can be rescued by treatments with NO donor, cyclic guanosine monophosphate analog, adenosine-monophosphate–activated protein kinase (AMPK) activator, or peroxisome proliferator-activated receptor alpha (PPAR-α) agonist. In contrast, inhibitors of NO synthases, AMPK, or soluble guanylate cyclase and liver-specifically expressed dominant negatives of peroxisome proliferator-activated receptor-gamma coactivator 1 alpha and PPAR-α are sufficient to induce hepatic steatosis in fasted wild-type larvae. Moreover, knockdown of Slc7a3 in mice or SLC7A3 in human liver cells impaired AMPK-PPAR-α signaling and resulted in lipid accumulation under fasting or glucose starvation, respectively. Conclusion: These findings have revealed a NO-AMPK-PPAR-α–signaling pathway that is crucial for the control of hepatic FAO in vertebrates. (Hepatology 2014;60:1928–1940)

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.