Comparative Transcriptome Analysis of Gill Tissue in Response to Hypoxia in Silver Sillago (Sillago sihama)

Abstract

Silver sillago (Sillago sihama) is a commercially important marine fish species in East Asia. In this study, we compared the transcriptome response to hypoxia stress in the gill tissue of S. sihama. The fish were divided into four groups, such as 1 h of hypoxia (hypoxia1h, DO = 1.5 ± 0.1 mg/L), 4 h of hypoxia (hypoxia4h, DO = 1.5 ± 0.1 mg/L), 4 h of reoxygen (reoxygen4h, DO = 8.0 ± 0.2 mg/L) after 4 h of hypoxia (DO = 1.5 mg/L), and normoxia or control (DO = 8.0 ± 0.2 mg/L) groups. Compared to the normoxia group, a total of 3550 genes were identified as differentially expressed genes (DEGs) (log2foldchange > 1 and padj < 0.05), including 1103, 1451 and 996 genes in hypoxia1h, hypoxia4h and reoxygen4h groups, respectively. Only 247 DEGs were differentially co-expressed in all treatment groups. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, DEGs were significantly enriched in steroid biosynthesis, biosynthesis of amino acids, glutathione metabolism and metabolism of xenobiotics by cytochrome P450, ferroptosis and drug metabolism—cytochrome P450 pathways. Of these, the cytochrome P450 (CYP) and glutathione S-transferase (GST) gene families were widely expressed. Our study represents the insights into the underlying molecular mechanisms of hypoxia stress.