As one of the critical abiotic factors, temperature controls fish development and reproduction. However, the effects of low temperature on the transcriptional regulation of zebrafish reproduction remain largely unclear. In this study, the fecundity of zebrafish was examined after exposure to cold temperatures at 19.5 °C, 19 °C, 18.5 °C, or 18 °C. The temperature at 19 °C showed no significant influence on the fecundity of zebrafish, but temperature at 18.5 °C or 18 °C significantly blocked the spawning of females, suggesting the existence of a low temperature critical point for the spawning of zebrafish females. Based on these observations, the brains of anesthetized fish under cold stress at different cold temperatures were collected for high-throughput RNA-seq assays. Key genes, hub pathways and important biological processes responding to cold temperatures during the spawning of zebrafish were identified through bioinformatic analysis. The number of down-regulated and up-regulated genes during the temperature reduction from egg-spawning temperatures at 19.5 °C and 19 °C to non-spawning temperatures at 18.5 °C and 18 °C were 2588 and 2527 (fold change ≥ 1.5 and p-value ≤ 0.01), respectively. Venn analysis was performed to identify up- and down-regulated key genes. KEGG enrichment analysis indicated that the hub pathways overrepresented among down-regulated key genes included the GnRH signaling pathway, vascular smooth muscle contraction, C-type lectin receptor signaling pathway, phosphatidylinositol signaling system and insulin signaling pathway. GO enrichment analysis of down-regulated key genes revealed the most important biological processes inhibited under non-spawning temperatures at 18.5 °C and 18 °C were photoreceptor cell outer segment organization, circadian regulation of gene expression and photoreceptor cell maintenance. Furthermore, 99 hormone-related genes were found in the brain tissues of non-spawning and spawning groups, and GnRH signaling pathway and insulin signaling pathway were enriched from down-regulated genes related to hormones at 18.5 °C and 18 °C. Thus, these findings uncovered crucial hormone-related genes and signaling pathways controlling the spawning of female zebrafish under cold stress.

Autophagy and endoplasmic reticulum (ER) stress response are among the key pathways regulating cold resistance of fish through eliminating damaged cellular components and facilitating the restoration of cell homeostasis upon exposure to acute cold stress. The transmembrane protein 39A (TMEM39A) was reported to regulate both autophagy and ER stress response, but its vertebrate-specific paralog, the transmembrane protein 39B (TMEM39B), has not been characterized. In the current study, we generate tmem39b-knockout zebrafish lines and characterize their survival ability under acute cold stress. We observed that the dysfunction of Tmem39b remarkably decreased the cold resilience of both the larval and adult zebrafish. Gene transcription in the larvae exposed to cold stress and rewarming were characterized by RNA sequencing (RNA-seq) to explore the mechanisms underlying functions of Tmem39b in regulating cold resistance. The results indicate that the deficiency of Tmem39b attenuates the up-regulation of both cold- and rewarming-induced genes. The cold-induced transcription factor genes bif1.2, fosab, and egr1, and the rewarming-activated immune genes c3a.3, il11a, and sting1 are the representatives influenced by Tmem39b dysfunction. However, the loss of tmem39b has little effect on the transcription of the ER stress response- and autophagy-related genes. The measurements of the phosphorylated H2A histone family member X (at Ser 139, abbreviated as γH2AX) demonstrate that zebrafish Tmem39b protects the cells against DNA damage caused by exposure to the cold-warming stress and facilitates tissue damage repair during the recovery phase. The gene modules underlying the functions of Tmem39b in zebrafish are highly enriched in biological processes associated with immune response. The dysfunction of Tmem39b also attenuates the up-regulation of tissue C-reactive protein (CRP) content upon rewarming. Together, our data shed new light on the function and mechanism of Tmem39b in regulating the cold resistance of fish.

The liver plays an essential role in multiple biological functions including metabolism, detoxification, digestion, coagulation, and homeostasis in vertebrates. The specification and differentiation of embryonic hepatoblasts, the proliferation of hepatocytes, and the hepatic tissue architecture are well documented, but molecular events governing the maturation of hepatocytes during liver development remain largely unclear. In this study, we performed a comparative transcriptome analysis of hepatocytes that were sorted by flow cytometry from developing zebrafish embryos at 60, 72, and 96 hpf. We identified 667 up-regulated and 3640 down-regulated genes in hepatocytes between 60 and 72 hpf, 606 up-regulated and 3924 down-regulated genes between 60 and 96 hpf, and 1693 up-regulated genes and 1508 down-regulated genes between 72 and 96 hpf. GO enrichment analysis revealed that key biological processes, cellular components, and molecular functions in hepatocytes between 60 to 72 hpf, such as cell cycle, DNA replication, DNA repair, RNA processing, and transcription regulation, are mainly associated with the proliferation of hepatocytes. In addition to biological processes, cellular components, and molecular functions for cell proliferation, molecular functions for carbohydrate metabolism were enriched in hepatocytes during 72 to 96 hpf. KEGG enrichment analysis identified key signaling pathways, such as cell cycle, RNA degradation, ubiquitin-mediated proteolysis, ErbB and Hedgehog signaling, basal transcription factors, Wnt signaling, and glycan degradation, which are closely associated with cell proliferation or carbohydrate metabolism in hepatocytes between 60 to 72 hpf. Newly enriched signaling pathways in hepatocytes during 72 to 96 hpf include metabolisms of pyrimidine, purine, nicotinate and nicotinamide, caffeine, glycine, serine and threonine, ABC transporters, and p53 signaling that function in metabolisms of lipid, protein and energy, cellular secretion, or detoxification, indicating the functional maturation of hepatocytes between 72 to 96 hpf. These findings provide novel clues for further understanding the functional differentiation and maturation of hepatocytes during liver development.