Integrative Genomics and Multi-Tissue Transcriptomics Identify Key Loci and Pathways for Hypoxia Tolerance in Grass Carp

Hypoxia is a critical environmental stressor in aquaculture, significantly affecting the survival and growth performance of cultured fish. To explore the genetic basis of hypoxia tolerance in grass carp (Ctenopharyngodon idella), we integrated genome-wide association analysis (GWAS) and multi-tissue transcriptome profiling. A total of 2000 grass carp were subjected to hypoxic stress, from which the 150 most hypoxia-intolerant (HI) and 150 most hypoxia-tolerant (HT) individuals were selected based on the time to loss of equilibrium (LOE). GWAS using 3,730,919 SNPs and 851,595 InDels identified 21 SNPs and 6 InDels associated with hypoxia tolerance. Two SNPs on chromosomes 10 and 13 reached genome-wide significance, accounting for 2.7% and 4.8% of the phenotypic variance explained (PVE), respectively. Validation of identified SNPs was performed using kompetitive allele-specific PCR (KASP) analysis. Candidate genes within ±50 kb of these variants were enriched in steroid biosynthesis, insulin signaling, and glycosphingolipid biosynthesis pathways. Transcriptomic analysis of six tissues (brain, gill, intestine, kidney, liver, and spleen) revealed 1620, 1221, 796, 246, 210, and 58 differentially expressed genes (DEGs) in the HT group compared to the HI group, respectively. DEGs in the brain were primarily enriched in steroid metabolic processes and angiogenesis regulation, while those in kidney and spleen DEGs were associated with oxygen transport and erythrocyte development. Integrated analysis of GWAS and transcriptome data identified 16 shared genes, including usf1 and trpv4. These findings reveal key genomic loci and molecular pathways underlying hypoxia tolerance in grass carp, providing valuable markers for future selective breeding programs.