棉花细菌性角斑病（CBB）是由革兰氏阴性细菌Xanthomonas citripv. malvacearum（Xcm）引起的毁灭性病害，在全球范围内严重威胁棉花生产。在Xcm已知的22个生理小种中，小种18在美国尤为流行且毒性最强。历史上，CBB曾造成巨大产量损失，而2011年后该病害在美国东南部棉区重新流行，部分抗性品种丧失抗性，凸显了发掘和部署新的抗性来源的紧迫性。
1. 研究背景与抗性基因挖掘
传统育种已鉴定了至少10个来自不同棉属物种的CBB主效抗性基因（B1-B10），其中一些被导入陆地棉栽培种。著名的多逆境抗性（MAR）育种项目成功将B2、B3等基因导入，为美国育种提供了抗性基础。非洲陆地棉品种S295对包括小种18在内的所有已知Xcm小种均表现出抗性，先前研究将其抗性定位在染色体D02上。近期，BB13抗性位点被初步定位。然而，在全球陆地棉种质中CBB抗性的多样性以及BB13位点的精细图谱仍不清晰。
2. 材料、方法与表型分析
本研究利用一个包含661份陆地棉种质（包括美国优良种质、热带地方品种和阿肯色大学种质）的多样性群体（DIV）进行全基因组关联研究（GWAS）。同时，构建了六个以感病品种Acala Maxxa为共同母本，以不同抗病种质为父本的杂交重组自交系（RIL）群体，用于连锁作图和验证。所有材料均在苗期接种Xcm小种18，根据接种部位产生水浸状病斑（感病）或过敏性坏死反应（抗病）进行表型鉴定。在多样性群体中，512份材料表型明确，其中133份（25.97%）抗病，379份（74.02%）感病。六个RIL群体均以1:1的比例分离抗感表型，证实了抗性由单显性基因控制。
3. 群体结构与抗性分布
对多样性群体的主成分分析（PCA）和群体结构分析揭示了三个主要类群：热带地方种、美国优良种质和阿肯色大学种质。抗性分布极不均衡，主要集中在与MAR育种项目相关的K04亚群（96.3%抗性）和与美国中南部育种项目相关的K08亚群（50.8%抗性）。而在以热带地方种为主的K01、K07、K11等亚群中，几乎全部感病，仅有一份来自巴拉圭的地方品种TX-2367表现抗性。这表明有效的CBB抗性基因在陆地棉地方种中非常稀有，当前美国种质中的抗性主要源自有限的供体材料。
4. 全基因组关联分析锁定主效位点
利用20,681个高质量单核苷酸多态性（SNP）标记进行GWAS分析，在512份材料中鉴定出107个与CBB抗性显著相关的标记-性状关联位点（MTA）。所有显著关联信号都集中在染色体D02长臂上一个3.31 Mb的区域内，没有在其他染色体上发现新的主效抗性位点。其中，SNP标记i46775Gh和i25755Gh关联最显著，分别解释了26.88%和26.77%的表型变异。这强有力地证实了Bb13是控制美国陆地棉对小种18抗性的最主要基因。
5. 精细作图与候选基因分析
为了精确定位BB13基因，研究者设计了针对目标区域的PCR等位基因竞争延伸（PACE）标记，并在多个RIL群体中进行连锁作图。最终，将BB13位点精细定位到染色体D02上物理位置介于981,030 bp和1,135,315 bp之间的一个154.28 kb的区间。该区间在参考基因组中包含15个推定的基因，其中6个的同源基因在拟南芥、水稻等模式植物中被证实参与抗病反应：
- •
Gh_D02G012100 (RGA4)和 Gh_D02G012200 (At4g27190)：编码含有NB-LRR结构域的典型抗病蛋白，可直接识别病原效应因子并触发过敏性细胞死亡。
- •
Gh_D02G012500 (NAC014)：NAC结构域转录因子，调控防御相关基因的表达。
- •
Gh_D02G012800 (GAI)：DELLA蛋白，通过调节赤霉素信号和活性氧（ROS）平衡参与胁迫响应。
- •
Gh_D02G013400 (ULP1B)和 Gh_D02G013500 (ESD4)：类泛素特异性蛋白酶，可能通过SUMO化修饰途径调控防御信号。
6. 开发实用分子标记
研究成功开发了与Bb13基因完全共分离的PACE标记CBB16，可用于准确追踪纯合抗性等位基因。同时，基于候选基因ESD4开发的CBB15标记，能有效区分杂合子与纯合子，在F2群体中呈现预期的1:2:1分离比，错误率低，非常适合用于分子标记辅助选择（MAS），加速抗病品种选育。
7. 连锁不平衡与单倍型分析
对目标区域进行连锁不平衡（LD）和局部单倍型分析，揭示了该区域的复杂遗传结构。单倍型B在119个抗病材料中高度富集，其单倍型组合（在三个标记组MG1、MG2、MG3上均携带替代等位基因）与抗性表型强烈相关。分析指出，MG2是决定抗性的关键区域。单倍型D（携带MG1和MG2的替代等位基因，但缺少MG3的）也表现抗性，表明MG1和MG2的组合足以在某些材料中赋予抗性。这些单倍型信息为理解抗性的等位变异和指导育种选择提供了更深层次的见解。
8. 抗性起源与等位性测试
为验证唯一抗病的热带地方种TX-2367是否携带新的抗性基因，将其与已知携带Bb13的抗病品种Arkot 8102杂交。F2群体中未出现感病单株，证明TX-2367的抗性与Bb13是等位的。结合群体遗传和谱系分析，本研究推测Bb13基因可能最初源自非洲品种S295（该品种携带著名的B12抗性基因），后通过MAR育种项目被导入美国陆地棉种质，并进一步扩散到多个育种群体中。BB13与B12很可能是同一基因，或是紧密连锁的基因。
9. 讨论与展望
本研究通过整合GWAS、连锁作图和单倍型分析，系统解析了陆地棉对CBB小种18抗性的遗传基础，明确Bb13是主导抗性基因，并完成了其精细定位。所开发的CBB15和CBB16等PACE标记为抗病育种提供了高效的分子工具。研究揭示了当前美国棉花抗性种质遗传基础相对狭窄，抗性等位基因分布不均的问题。未来，需要利用已鉴定的候选基因进行功能验证和图位克隆，彻底解析Bb13的抗病分子机制。同时，应继续从野生棉种和地方品种中挖掘新的抗性等位基因，与Bb13进行基因聚合，以培育具有广谱、持久抗性的棉花新品种，应对不断进化的病原菌威胁。
