Trachinotus ovatus is an economically important mariculture fish, and hypoxia has become a critical threat to this hypoxia-sensitive species. However, the molecular adaptation mechanism of T. ovatus liver to hypoxia remains unclear. In this study, we investigated the effects of acute hypoxic stress (1.5 ± 0.1 mg·L⁻¹ for 6 h) and re-oxygenation (5.8 ± 0.3 mg·L⁻¹ for 12 h) in T. ovatus liver at both the transcriptomic and metabolic levels to elucidate hypoxia adaptation mechanism. Integrated transcriptomics and metabolomics analyses identified 36 genes and seven metabolites as key molecules that were highly related to signal transduction, cell growth and death, carbohydrate metabolism, amino acid metabolism, and lipid metabolism, and all played key roles in hypoxia adaptation. Of these, the hub genes FOS and JUN were pivotal hypoxia adaptation biomarkers for regulating cell growth and death. During hypoxia, up-regulation of GADD45B and CDKN1A genes induced cell cycle arrest. Enhancing intrinsic and extrinsic pathways in combination with glutathione metabolism triggered apoptosis; meanwhile, anti-apoptosis mechanism was activated after hypoxia. Expression of genes related to glycolysis, gluconeogenesis, amino acid metabolism, fat mobilization, and fatty acid biosynthesis were up-regulated after acute hypoxic stress, promoting energy supply. After re-oxygenation for 12 h, continuous apoptosis favored cellular function and tissue repair. Shifting from anaerobic metabolism (glycolysis) during hypoxia to aerobic metabolism (fatty acid β-oxidation and TCA cycle) after re-oxygenation was an important energy metabolism adaptation mechanism. Hypoxia 6 h was a critical period for metabolism alteration and cellular homeostasis, and re-oxygenation intervention should be implemented in a timely way. This study thoroughly examined the molecular response mechanism of T. ovatus under acute hypoxic stress, which contributes to the molecular breeding of hypoxia-tolerant cultivars.

Trachinotus ovatus is a major economically important cultured marine fish in the South China Sea. However, extreme weather and increased culture density result in uncontrollable problems, such as increases in water temperature and a decline in dissolved oxygen (DO), hindering the high-quality development of aquaculture. In this study, liver transcriptional profiles of T. ovatus were investigated under acute high-temperature stress (31 °C and 34 °C) and normal water temperature (27 °C) using RNA sequencing (RNA-Seq) technology. Differential expression analysis and STEM analysis showed that 1347 differentially expressed genes (DEGs) and four significant profiles (profiles 0, 3, 4, and 7) were screened, respectively. Of these DEGs, some genes involved in heat shock protein (HSPs), hypoxic adaptation, and glycolysis were up-regulated, while some genes involved in the ubiquitin-proteasome system (UPS) and fatty acid metabolism were down-regulated. Our results suggest that protein dynamic balance and function, hypoxia adaptation, and energy metabolism transformation are crucial in response to acute high-temperature stress. Our findings contribute to understanding the molecular response mechanism of T. ovatus under acute heat stress, which may provide some reference for studying the molecular mechanisms of other fish in response to heat stress.

Savalani hairtail Lepturacanthus savala is a widely distributed fish along the Indo-Western Pacific coast, and contributes substantially to trichiurid fishery resources worldwide. In this study, the first chromosome-level genome assembly of L. savala was obtained by PacBio SMRT-Seq, Illumina HiSeq, and Hi-C technologies. The final assembled L. savala genome was 790.02 Mb with contig N50 and scaffold N50 values of 19.01 Mb and 32.77 Mb, respectively. The assembled sequences were anchored to 24 chromosomes by using Hi-C data. Combined with RNA sequencing data, 23,625 protein-coding genes were predicted, of which 96.0% were successfully annotated. In total, 67 gene family expansions and 93 gene family contractions were detected in the L. savala genome. Additionally, 1825 positively selected genes were identified. Based on a comparative genomic analysis, we screened a number of candidate genes associated with the specific morphology, behaviour-related immune system, and DNA repair mechanisms in L. savala. Our results preliminarily revealed mechanisms underlying the special morphological and behavioural characteristics of L. savala from a genomic perspective. Furthermore, this study provides valuable reference data for subsequent molecular ecology studies of L. savala and whole-genome analyses of other trichiurid fishes.

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-β superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.