《Molecular Plant Pathology》：GWAS With a Global Population of Wheat Stripe Rust Identified High-Confidence Avirulence Gene Candidates
编辑推荐：
小麦条锈病由小麦条锈菌（Puccinia striiformis f. sp. tritici，Pst）引起，是一种全球性毁灭性病害，其特征是毒力快速进化。尽管其重要性，全球范围内无毒力变异的遗传基础仍知之甚少。在本研究中，研究人员将针对18个Yr抗性基因的毒
小麦条锈病由小麦条锈菌（Puccinia striiformis f. sp. tritici，Pst）引起，是一种全球性毁灭性病害，其特征是毒力快速进化。尽管其重要性，全球范围内无毒力变异的遗传基础仍知之甚少。在本研究中，研究人员将针对18个Yr抗性基因的毒力表型与来自全球范围Pst分离株的基因组变异相结合，以解析无毒力的群体基因组结构。研究人员鉴定出多个无毒力相关基因组区域，并定义了69个候选效应子，其中11个高置信候选效应子满足基于效应子特征和关联显著性的五项严格标准。等位基因频率分析揭示了显著的地理分化；许多变异在特定区域（尤其是东非、南美和中国内陆）呈高频或近乎固定状态，而在其他区域则保持稀有。这一模式表明，由宿主识别基因的区域性局限多态性驱动的强群体结构。此外，结构建模支持抗性蛋白Yr9与多个候选AvYr9效应子之间的直接分子识别，为特异性提供了机制学见解。总之，本研究建立了效应子发现的全球群体基因组框架，并阐明了Pst的毒力进化。这些候选效应子代表了理解锈病致病性和开发小麦持久抗性的宝贵资源。
小麦条锈病由小麦条锈菌（Puccinia striiformis f. sp. tritici，Pst）引起，是全球小麦生产最持久的生物威胁之一，反复造成亚洲、欧洲、非洲及美洲的严重产量损失。尽管利用抗病品种是最有效且环境可持续的防控策略，但抗性常在品种释放后数年内被克服，这反映了Pst群体强大的进化潜力。植物免疫多遵循基因对基因模型：宿主抗性（R）基因识别病原无毒（Avr）基因，触发效应子触发免疫（ETI）。Avr基因通常编码分泌效应蛋白，在抑制宿主基础防御的同时也暴露给宿主免疫识别。然而，与已鉴定超过80个小麦条锈病抗性（Yr）基因相比，Pst中对应Avr基因的认知极为有限，迄今仅AvrYr7被明确克隆。这种滞后源于Pst作为专性活体营养菌的生物学特性：基因组庞大且富含重复序列，侵染阶段为双核状态，难以进行遗传操作和经典图位克隆。近年来，群体规模基因组测序和统计基因组学为克服这些限制提供了新机会，全基因组关联分析（GWAS）无需杂交即可在自然病原群体中鉴定与表型相关的遗传变异。因此，研究人员整合了全球234个Pst分离株在18个Yr单基因系上的毒力表型与全基因组变异，利用GWAS解析无毒力变异的遗传基础，并进一步通过连锁不平衡（LD）区间定义、效应子候选优先化、比较基因组学和群体遗传学分析，系统鉴定候选Avr效应子及其地理分化，最后结合蛋白质结构预测和体内实验验证候选效应子与抗性蛋白Yr9的互作。研究建立了全球群体基因组框架，为理解Pst毒力进化和开发持久抗性提供了重要资源。
研究人员以全球六大洲234个Pst分离株为对象，利用18个Yr单基因系测定的侵染型（IT）表型（0-9标度，0-6为有毒，7-9为无毒），结合Illumina NovaSeq 6000全基因组重测序数据（比对至参考基因组CYR34，GATK变异检测），采用rMVP中的FarmCPU和混合线性模型进行全基因组关联分析（显著性p≤1.0×10
?8），用PopLDdecay估计LD衰减区间；通过SignalP/TMHMM、EffectorP、OrthoFinder、HMMSCAN等流程筛选效应子；用AlphaFold3预测蛋白互作，用农杆菌瞬时表达和Split-LUC验证互作。
2.1 表型统计与性状间相关性
对234个Pst分离株在18个Yr单基因系上的表型分析显示，各性状存在广泛变异。相关性分析揭示了模块化结构：Yr1、Yr9和Yr76呈强正相关；Yr24与Yr32相关性最高（r=0.84，p<0.001）；而Yr15与YrExp2等性状大多独立，与其余性状相关性弱且多不显著；少数存在负相关，如Yr1与Yr8。这表明不同Yr基因对应的毒力变异可能受不同的进化或功能约束。
2.2 与Yr基因对应的候选Avr基因的GWAS
采用FarmCPU模型鉴定出16个无毒力（Avr）位点（AvYr5和AvYr15因所有分离株表现一致无毒而被排除），MLM仅检测到8个位点。整合两种模型后，AvYr9关联基因数最多（134个），其次AvYr1（122个）、AvYr44（78个）。全基因组平均LD衰减距离为26.66 kb，以此定义关联SNP上下游置信区间，联合模型在LD区间内共鉴定出与AvYr1、AvYr9、AvYr44、AvYr6等位点关联的多个基因。
2.3 与Avr位点相关的候选信号肽基因鉴定
通过SignalP预测和TMHMM排除跨膜螺旋，在显著关联SNP所在基因和±26.66 kb区间内共鉴定出77个候选SP基因，去冗余后得到69个。其中Pst_SP12同时关联AvYr1/AvYr9，Pst_SP13关联AvYr1/AvYr44/AvYrExp2，Pst_SP14关联AvYr1/AvYr17，Pst_SP15关联AvYr1/AvYr9/AvYr44。
2.4 Pst效应子候选的特征
69个候选SP基因经五项标准评估：短氨基酸长度（<300 aa）、半胱氨酸丰富（≥4个Cys）、EffectorP预测、属/种特异性、无已知功能域及种内多态性。EffectorP预测其中37个为候选效应子；Pfam和InterPro注释显示55个缺乏保守结构域。比较基因组分析从14个真菌分离株中鉴定出51个Puccinia或P. striiformis特异蛋白。多态性分析显示66个候选基因在高质量基因组中存在多态。亚细胞定位预测提示31.9%的候选效应子可能靶向宿主细胞核。与127个已报道Pst效应子比对，13个有显著相似性，其中Pst_SP09与已知效应子KNE95698.1完全匹配；AvrYr7等位基因与本研究候选物无显著相似性。
2.5 全球群体中Avr相关多态性的地理分化
等位基因频率分析显示，Avr关联SNP存在强烈地理分化，约45.4%位于CM045852.1染色体上，多个变异在至少4个地理群体中高频或近固定，在其他群体中稀有。AvYr44关联SNP在所有8个群体中频率均高。东非（G3）和南美（G2）群体表现出与其他区域明显不同的频率谱。对AvYr1、AvYr6、AvYr9和AvYr44四个位点的详细分析揭示同一座位内不同变异存在区域特异的相反富集模式，提示局部选择或扩张事件。
2.6 毒性相关Avr位点表现出群体区分能力
基于所有Avr位点高置信SNP构建的系统发育树显示，东非、南美和中国内陆分离株分别形成优势聚类（MG2中57.1%为东非，MG3中79.2%为南美，MG7中88.9%为中国内陆），而欧洲与新月沃地等群体存在混合。判别主成分分析（DAPC）支持保留14个主成分，前三个判别轴捕获大部分群体间变异，东非、南美和中国内陆形成不重叠簇，欧洲分离株显示部分混合。这些结果表明Avr相关位点携带强烈的群体判别信号。
2.7 Pst_SP05作为与Yr9互作的AvYr9候选的体内验证和结构预测
针对AvYr9最大的功能簇（5个高置信候选效应子），在本氏烟草中瞬时表达后，Pst_SP05和Pst_SP08能完全抑制BAX蛋白诱导的细胞死亡，表明其具有免疫抑制活性。Split-LUC实验证实Pst_SP05与Yr9存在物理互作，而Pst_SP08无显著互作。AlphaFold3结构预测显示，全长Yr9预测置信度较低（ipTM=0.21），而Pst_SP05与Yr9的富亮氨酸重复（LRR）结构域组合达到高置信度（ipTM=0.73，pTM=0.53，排序得分=0.95），证明LRR结构域是识别Pst_SP05的主要结构决定簇，关键残基如Pst_SP05的Arg-17与Yr9的Glu-26/Ser-84形成盐桥和极性接触。
讨论部分强调，该研究的优势在于采用双模型整合策略（FarmCPU和MLM）并应用严格显著性阈值（p<1.0×10
?8），有效控制了假阳性。通过五维标准从69个候选SP基因中筛选出11个高置信候选，且它们与单个Avr位点呈一一对应关系，符合经典基因对基因模型。多个候选效应子与多个Avr位点关联（如Pst_SP13和Pst_SP15），提示单一效应子等位基因多态可能被多个R基因识别；而某些关联（如Yr1与Yr44表型负相关）暗示适应性代价可能限制多毒力基因型的积累。亚细胞定位预测提示效应子功能多样化，包括核靶向效应子。等位基因频率和群体结构分析表明，Avr位点的地理分化主要由区域部署的Yr基因施加的局部选择驱动，而非中性过程。系统发育和DAPC显示Avr位点比中性标记更能捕捉适应性群体结构。结构预测和实验验证共同支持Pst_SP05是AvYr9的候选效应子，并建立了基于AlphaFold3的NLR-效应子互作预测框架。
研究结论可概括为：区域抗性基因部署驱动Avr位点的局部选择，塑造等位基因频率和群体结构；效应子与含核苷酸结合域和富亮氨酸重复的受体蛋白（NLR）互作的分子多样化构成毒力进化的机制基础；Avr相关基因组区域是连接生态选择、群体历史和分子识别的进化热点，为全球监测和持久抗性育种提供了直接资源。
