生物技术进展 ›› 2025, Vol. 15 ›› Issue (3): 432-445.DOI: 10.19586/j.2095-2341.2024.0184

• 研究论文 • 上一篇    下一篇

联合GWAS和转录组分析挖掘与白斑狗鱼耐高温性状关联的SNP

海萨·艾也力汗null1(), 焦飞1, 刘鸿1, 呼德拉提·阿那斯null1, 沈玉帮2()   

  1. 1.新疆维吾尔自治区水产科学研究所,乌鲁木齐 830000
    2.上海海洋大学水产与生命学院,上海 201306
  • 收稿日期:2024-11-24 接受日期:2025-02-17 出版日期:2025-05-25 发布日期:2025-07-01
  • 通讯作者: 沈玉帮
  • 作者简介:海萨·艾也力汗 E-mail:hbahjan@126.com
  • 基金资助:
    国家自然科学基金项目(32060826)

Integrating GWAS and Transcriptome Profiling to Identify SNP Markers Linked to High-temperature Tolerance in Esox lucius

AYELHAN·Haysa1(), Fei JIAO1, Hong LIU1, ANASI·Hudelati1, Yubang SHEN2()   

  1. 1.Xinjiang Fishery Research Institute, Urumqi 830000, China
    2.College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306, China
  • Received:2024-11-24 Accepted:2025-02-17 Online:2025-05-25 Published:2025-07-01
  • Contact: Yubang SHEN

摘要:

研究旨在提升白斑狗鱼对极端高温的适应能力和满足向南部温暖地区养殖发展的需求,进一步发掘与其耐高温性状显著关联的基因位点和候选基因。首先基于前期对热敏感与耐高温群体的简化基因组测序数据,利用混合线性模型(mixed-linear model, MLM)进行全基因组关联分析(genome-wide association study, GWAS)筛选出与耐高温性状显著关联的单核苷酸多态性(single nucleotide polymorphism, SNP)位点;继而采用转录组测序技术对4组不同热处理幼鱼的脑组织进行基因差异表达分析;随后联合GWAS和转录组分析筛选出候选基因和SNP位点;最后利用Sanger测序技术,以耐高温组和参考组为材料,对显著关联的SNP位点加以验证。结果显示,GWAS分析获得了471个与耐高温性状显著关联的SNP位点。GO和KEGG分析显示,DEGs显著富集于内质网应激和免疫系统相关通路。GSEA分析显示,热胁迫84 h后先天免疫反应被激活,并且与NAD(+)ADP-核糖基转移酶活性相关的基因被诱导。先天免疫反应启动和调控相关基因的PPI分析显示,parp14parp9stat1stat3是核心基因,同时也是互作基因。热应激反应中parp14可能通过调控stat1stat3等下游效应基因来发挥保护细胞的作用。GWAS获得的SNP位点中,有2个SNP位点位于parp14基因的第三外显子。验证结果显示,parp14基因第三外显子中的2个SNP位点与耐高温性状显著相关,且为完全连锁不平衡,其中位点rsc.646T>C为错义突变,位点rsc.777G>A为同义突变。研究筛选出与白斑狗鱼耐高温相关的候选基因及与耐高温性状显著关联的2个SNP位点,可为白斑狗鱼耐高温性状的遗传机制和分子标记辅助育种提供理论依据。

关键词: 白斑狗鱼, 耐高温性状, GWAS, 转录组, parp14, SNP

Abstract:

To enhance the extreme high-temperature adaptability of Esox lucius for aquaculture expansion to southern warm regions and identify genetic loci associated with thermal tolerance, we integrated genome-wide association study (GWAS) and transcriptome analysis. Using reduced-representation genome sequencing data from heat-sensitive and heat-tolerant populations, GWAS analysis with an MLM model identified 471 SNPs significantly associated with thermal tolerance traits. Concurrently, transcriptome sequencing of juvenile brain tissues exposed to varying heat stress durations revealed differential gene expression patterns. Functional enrichment analyses demonstrated significant enrichment of differentially expressed genes (DEGs) in endoplasmic reticulum stress and immune-related pathways. GSEA analysis indicated activation of innate immune responses at 84 h of heat stress, with specific induction of NAD(+) ADP-ribosyltransferase activity-related genes. Protein-protein interaction (PPI) network analysis identified parp14parp9stat1, and stat3 as core regulators of innate immune responses. Notably, two GWAS-derived SNPs located in exon 3 of the parp14 gene—rsc.646T>C (missense mutation) and rsc.777G>A (synonymous mutation)—exhibited complete linkage disequilibrium and were significantly correlated with thermal tolerance traits. Sanger sequencing validation in heat-stressed and reference populations confirmed these associations. Our findings reveal candidate genes and functional SNP markers for high-temperature tolerance in E. lucius, providing insights into its genetic mechanisms and supporting marker-assisted breeding strategies.

Key words: Esox lucius, high temperature tolerance, GWAS, transcriptome, parp14, SNP

中图分类号: