T. fasciatus (Takifugu fasciatus) faces the same problem as most warm water fish: the water temperature falls far below the optimal growth temperature in winter, causing a massive death of T. fasciatus and large economic losses. Understanding of the cold-tolerance mechanisms of this species is still limited. Integrated application of multi-omics research can provide a wealth of information to help us improve our understanding of low-temperature tolerance in fish. To gain a comprehensive and unbiased molecular understanding of cold-tolerance in T. fasciatus, we characterized mRNA-seq and metabolomics of T. fasciatus livers using Illumina HiSeq 2500 and UHPLC-Q-TOF MS. We identified 2544 up-regulated and 2622 down-regulated genes in the liver of T. fasciatus. A total of 40 differential metabolites were identified, including 9 down-regulated and 31 up-regulated metabolites. In combination with previous studies on proteomics, we have established an mRNA-protein-metabolite interaction network. There are 17 DEMs (differentially-expressed metabolites) and 14 DEGs-DEPs (differentially co-expressed genes and proteins) in the interaction network that are mainly involved in fatty acids metabolism, membrane transport, signal transduction, and DNA damage and defense. We then validated a number of genes in the interaction network by qRT-PCR. Additionally, a number of SNPs (single nucleotide polymorphisms) were revealed through the transcriptome data. These results provide key information for further understanding of the molecular mechanisms of T. fasciatus under cold stress. The data generated by integrated application of multi-omics can facilitate our understanding of the molecular mechanisms of fish response to low temperature stress. We have not only identified potential genes and SNPs involved in cold tolerance, but also show that some nutrient metabolites may be added to the diet to help the overwintering of T. fasciatus.

Hypoxia frequently occurs in aquatic ecosystem, which is influenced by salinity, water temperature, weather, and surface water runoff. In order to shed further light on the evolutionary and adaptive mechanisms in fish under hypoxic condition, the impact of acute hypoxia (1.63 ± 0.2 mg/L) and reoxygenation (7.0 ± 0.3 mg/L) on oxygen sensors, energy metabolism, and hematological indices was evaluated in Takifugu fasciatus. Data from transcriptional level analysis show that the expressions of genes related to oxygen sensors (HIF-1α, PHD2, and VHL) were upregulated in the brain and liver under hypoxia and recovered under reoxygenation. The upregulation of GLUT2, VEGF-A, and EPO in conjugation with VEGF-A protein and hematological indices conferred the rapid adjustments of cellular glucose uptake and blood oxygen-carrying capacities in pufferfish. Higher levels of glycolysis-related mRNAs (HK, PGK1, and PGAM2), HK activity, and proteins (PGK1 and PGAM2) were detected in the brain and liver under hypoxic condition compared with control. Interestingly, the expression of MDH1 at the mRNA, enzyme activity, and protein levels was significantly increased in the brain at 0 or 2 h and in the liver at 8 h under hypoxic condition. In addition, although the enzyme activity and mRNA expression of LDH in the brain were not significantly changed, a persistent upregulation was observed in the liver during hypoxia exposure. This study demonstrated that pufferfish could counterpoise the energetic demands and hematological functional properties evoked by oxygen sensors after hypoxia. Our findings provided new insights into the molecular regulatory mechanism of hypoxia in pufferfish.