《The Plant Genome》：Genome-wide association study reveals the genetic architecture of drought tolerance in maize using yield-based resistance indices and transcriptomic evidence
编辑推荐：
干旱是全球范围内限制玉米(Zea mays L.)生产的主要环境因子。解析耐旱性的遗传基础并精准定位关键位点与候选基因，是培育耐旱性增强玉米品种分子育种工作的核心基础。本研究对200份玉米自交系在2022与2023两个连续年份中，分别于水分胁迫(WS)与充分灌
干旱是全球范围内限制玉米(Zea mays L.)生产的主要环境因子。解析耐旱性的遗传基础并精准定位关键位点与候选基因，是培育耐旱性增强玉米品种分子育种工作的核心基础。本研究对200份玉米自交系在2022与2023两个连续年份中，分别于水分胁迫(WS)与充分灌溉(WW)两种对立水分条件下进行表型评价。研究人员测定了成熟期籽粒产量及抗旱指数(DRI)，其中DRI作为基于产量的指标，能够更全面地表征耐旱性。基于FarmCPU模型开展的基因组关联分析共检测到126个显著单核苷酸多态性(SNP)位点，其中43个位点与WS条件下的产量相关（表型方差解释率[PVE] = 4.64%–10.67%），55个位点与DRI相关（PVE = 1.47%–10.67%）。结合两套独立的RNA-seq数据集，研究人员鉴定出43个核心候选基因，其中29个获得功能注释。基因本体论(GO)富集分析与蛋白互作(PPI)分析表明，这些基因主要参与代谢调控、信号转导及环境胁迫响应相关的生物学过程。值得注意的是，包括NAC35 (GRMZM5G813651)、MADS-box蛋白 (GRMZM2G148693)及C3HC4型RING锌指蛋白 (GRMZM2G147319)在内的多个转录因子，此前已被证实参与干旱响应通路。该研究通过采用基于成熟期植株的产量抗旱指数，更直接地反映了田间条件下的耐旱性，为解析玉米耐旱性的遗传结构提供了新的见解。
研究背景
玉米(Zea mays L.)是全球重要的粮食、饲料及工业原料作物，其产量稳定性对国家粮食安全至关重要。干旱是限制玉米生产的主要非生物胁迫因子，可导致产量下降1%至76%。虽然传统育种与生物技术是应对干旱的有效手段，但干旱胁迫的复杂性导致常规育种周期长、效率低。尽管全基因组关联研究(GWAS)已广泛应用于玉米复杂性状的遗传解析，但既往研究多聚焦于苗期性状或间接指标（如散粉吐丝间隔、根系构型等），针对大田条件下成熟期最终产量水平的耐旱性遗传基础研究仍显不足。为此，研究人员采用基于籽粒产量的抗旱指数(DRI)，旨在更直接地揭示玉米耐旱性的遗传架构。
技术方法
研究人员构建了包含200份玉米自交系的自然群体，于2022至2023年在山西晋中和忻州两地进行大田试验，设置充分灌溉(WW)与水分胁迫(WS)处理。在成熟期测定籽粒产量(GY)，并计算抗旱指数(DRI)。利用48K SNP芯片进行基因型分型，使用FarmCPU模型进行全基因组关联分析，结合R语言进行最佳线性无偏预测(BLUP)及表型方差分析。通过整合NCBI GEO数据库中的GSE71723和GSE132113两套转录组数据，筛选差异表达基因(DEG)，并利用DAVID和STRING数据库进行基因本体(GO)富集及蛋白互作网络分析。
研究结果
3.1 表型变异评估
两年间，水分胁迫导致籽粒产量显著降低，2022年降幅达27.58%，2023年为14.15%。方差分析显示基因型是导致产量变异的主导因素，其次是环境因素。籽粒产量的广义遗传力(H2)在充分灌溉条件下略高于胁迫条件，分别为0.63和0.62。
3.2 抗旱指数与胁迫下产量的GWAS
基因组膨胀因子(λ)介于0.97至1.04之间，表明模型有效控制了假阳性。研究共检测到126个显著SNP位点。其中55个位点与DRI相关，PVE范围为1.47%–10.67%；43个位点与WS下的产量相关，PVE范围为4.64%–10.67%。有6个SNP位点在两年中均稳定与DRI关联，位于第1、4、8和10号染色体上。
3.3 DRI关联位点的等位变异效应
对4个稳定SNP进行等位效应分析发现，有利等位基因均显著提高DRI值。例如，Affx-291407911位点的C等位基因相较于T等位基因，在2022年和2023年分别使DRI提高了0.13和0.16。群体中优势等位基因的频率平均为46.49%，且DRI值与优势等位基因的累加呈正相关，携带近50个优势等位基因的自交系DRI值可达1.3。
3.4 干旱相关候选基因鉴定
在DRI关联位点上下游200 kb范围内共鉴定出751个候选基因，结合转录组数据筛选出43个核心候选基因。功能注释显示，GRMZM5G813651编码NAC结构域蛋白35，GRMZM2G148693编码MADS-box蛋白，GRMZM2G147319编码C3HC4型RING锌指蛋白。GO富集分析表明，这些基因显著富集于有机物质代谢、信号转导（如生长素介导的信号通路）及逆境胁迫响应等生物学过程。
3.5 候选基因蛋白互作分析
蛋白互作网络分析显示，7个候选基因与其他功能蛋白存在相互作用。例如GRMZM2G056786（编码RING-H2 finger ATL16）、GRMZM2G031615（编码生长素响应因子ARF）及GRMZM2G144668（编码ACC氧化酶）等，这些基因在植物耐热性、根系发育及激素信号传导中发挥重要作用。
讨论与结论
讨论部分指出，以往研究多采用苗期或生理指标，而本研究采用的DRI直接反映了全生育期的耐旱性，具有更高的育种参考价值。研究发现的6个稳定位点（如Affx-291389821和Affx-88998430）代表了稳健的遗传信号。此外，第3、4、5、10号染色体的特定基因组区域被确定为耐旱热点区域。研究人员强调，NAC35、MADS-box及RING finger家族基因在转录水平调控干旱响应，是关键的候选靶点。
结论部分总结，本研究通过两年的水分梯度试验，利用FarmCPU模型定位了126个与产量及DRI相关的显著SNP。其中6个位点在不同年份重复检出，结合转录组分析锁定了43个核心候选基因，主要参与代谢调控、信号转导及胁迫响应。这些稳定的位点与基因为后续的标记辅助选择及基因组选择育种提供了重要的基因资源，将加速耐旱玉米品种的培育进程。
