大豆NLR基因GmDBA1与GmDBA2调控矮化表型并增强对大豆疫霉（Phytophthora sojae）抗性的候选基因鉴定

《Plants》：GmDBA1 and GmDBA2 Are Candidate NLR Genes Associated with Dwarfism and Enhanced Resistance to Phytophthora sojae in Soybean

　　

由大豆疫霉（Phytophthora sojae）引起的疫霉根腐病（Phytophthora root and stem rot, PRR）严重降低大豆产量。研究人员从大豆品种合豆12（Hedou 12, H12）的甲基磺酸乙酯（ethyl methanesulfonate, EMS）诱变群体中分离得到隐性突变体gmdba1。该突变体表现出矮化、节数减少、分枝增多及叶片皱缩等表型。通过遗传定位，将与表型相关的位点界定于18号染色体上约154 kb区间内，该区间包含22个注释基因。经重测序和靶向序列分析，研究人员优先锁定两个NLR候选基因，即Glyma.18g226500（GmDBA1，矮化与基础免疫激活1，Dwarfism and Basal Activation 1）和Glyma.18g226800（GmDBA2），二者共含有四个突变特异性单核苷酸多态性（single nucleotide polymorphism, SNP）位点（其中三个为非 synonymous 替换）。gmdba1突变体以及一株在两个基因上同时编辑的纯合T2株系，在接种大豆疫霉P6497菌株后均表现出减轻的病害症状。上述结果支持GmDBA1/GmDBA2所在位点的扰动与发育改变及在测试条件下抗性增强之间存在关联。转录组和激素分析显示，gmdba1中防御相关基因、茉莉酸（jasmonic acid, JA）和水杨酸（salicylic acid, SA）含量呈组成性激活，这可能解释所观察到的生长-防御权衡现象。综上所述，研究结果表明GmDBA1和GmDBA2抑制对大豆疫霉的基础抗性，为培育抗病大豆品种提供了有价值的遗传资源。

一、研究背景与科学问题

大豆（Glycine max）是全球重要的经济作物，其产量长期受到由卵菌、真菌、细菌、病毒及线虫等多种病原侵害的严重制约。其中，由大豆疫霉（Phytophthora sojae）引起的疫霉根腐病（Phytophthora root and stem rot, PRR）是极具破坏性的土传卵菌病害，可在大豆从苗期到成熟期的整个发育阶段造成侵染。由于大豆疫霉的土传特性及其快速克服寄主抗性的能力，该病害的防控面临重大挑战。寄主遗传抗性仍是PRR防治的核心策略，但大多数经典的大豆疫霉抗性（Resistance to P. sojae, Rps）基因仅赋予小种特异性抗性，易因病原群体致病性变异而失效。因此，大豆对疫霉的抗性日益被视为由主效Rps位点与数量抗性组分共同构成的复合性状，这为发掘不依赖特定病原基因型的抗性机制提供了契机。尽管大豆中已定位了大量赋予数量抗性的数量性状位点（quantitative trait locus, QTL），但成功克隆并完成功能鉴定的抗性基因仍然有限。

核苷酸结合富亮氨酸重复受体（nucleotide-binding leucine-rich repeat receptor, NLR）是激活防御反应的重要胞内免疫受体，可单独、以遗传连锁配对或作为受体网络组分发挥功能。NLR信号调控失常可导致组成性防御相关表型及生长代价。植物免疫与发育密切相关，持续性的防御信号往往伴随营养生长减弱和株型改变。茉莉酸（jasmonic acid, JA）与水杨酸（salicylic acid, SA）是兼具免疫与形态建成调控功能的双功能植物激素，二者在不同组织、发育阶段及病原情境下可呈现拮抗或协同效应。基于上述背景，本研究旨在通过正向遗传学手段鉴定参与大豆生长与免疫调控的新基因，并阐明其潜在机制。

二、研究内容与主要结论

研究人员从高产耐逆中国大豆品种合豆12（H12）的EMS诱变群体中筛选获得隐性突变体gmdba1。该突变体在整个营养与生殖生长阶段均表现出多效性矮化表型，株高较野生型降低超过60%，节数减少、分枝数增加，且新出叶片呈现明显皱缩形态。遗传分析表明，gmdba1×Williams 82（W82）杂交产生的F₂群体中，野生型表型与矮化植株数分别为161株和44株，符合单隐性基因座3:1分离模式（χ² = 1.367, p = 0.242）。通过全基因组InDel标记初筛将目标位点初步定位于18号染色体，随后利用94株和131株突变体F₂单株逐步缩小定位区间，最终将候选区间界定于AGm18-17与AGm18-24之间约154 kb的物理区域内，该区间包含22个注释基因。对全部22个基因进行Sanger测序后排除20个与野生型序列一致的基因，剩余两个基因即Glyma.18g226500（GmDBA1）和Glyma.18g226800（GmDBA2）被优先锁定为候选基因。池化全基因组重测序鉴定出四个突变特异性SNP，其中三个为非同义替换。GmDBA1全长43,083 bp，含9个外显子和8个内含子，编码2882个氨基酸；GmDBA2全长8610 bp，含6个外显子和5个内含子，编码2082个氨基酸。二者编码序列一致性达69%，氨基酸一致性达87%，且均含有预测的NB-ARC和LRR结构域。系统发育分析表明二者为近期共同祖先衍生的大豆旁系同源基因。

组织特异性表达分析显示，GmDBA1/2在根组织中表达量最高，约为其他检测组织（顶端分生组织、花、嫩荚、叶片、根瘤）的4倍，与大豆疫霉主要侵染根组织的生物学特性相符。接种致病性大豆疫霉P6497菌株后，gmdba1突变体离体叶片和胚轴均表现出显著减小的病斑面积和病原生物量，呈现典型的超敏反应特征；扫描电镜观察显示突变体叶片表皮蜡质沉积不规则且表皮细胞形态改变。为提供直接遗传证据，研究人员利用CRISPR/Cas9系统在W82背景下同时编辑GmDBA1和GmDBA2，获得纯合T₂双敲除株系（GmDBA1/2-CR）。测序确认GmDBA1发生216 bp缺失（造成72个氨基酸的框内缺失），GmDBA2发生160 bp缺失（导致翻译提前终止）。该双编辑株系精确重现了gmdba1突变体的矮化、节间缩短和叶片皱缩表型，并在接种试验中获得显著抗性增强，与原始突变体表型高度一致。

转录组分析（RNA-seq）在V₆期顶端组织中共鉴定出大量差异表达基因（differentially expressed gene, DEG），KEGG和GO富集分析显示这些基因显著富集于植物-病原互作、MAPK信号级联以及JA和SA激素生物合成与信号转导等通路。qRT-PCR验证表明，SA防御标志基因GmPR1、JA-Ile合成基因GmJAR1以及WRKY和MYB转录因子等关键免疫调控因子在gmdba1中均显著上调。激素定量检测进一步证实，gmdba1中JA、JA-Ile、JA-Val、MeJA以及SA及其储存形式SAG的含量均显著高于野生型。上述结果共同表明，GmDBA1/GmDBA2位点在基础条件下参与抑制SA和JA依赖的免疫输出，该模块的突变导致激素过度积累和下游免疫信号激活，伴随组成性抗性和严重发育缺陷。

三、主要技术方法

本研究综合运用了多种关键技术手段。遗传定位方面，研究人员以W82×gmdba1构建F₂作图群体，采用全基因组InDel标记初筛结合高分辨率重组单株精细定位，将目标位点锁定于18号染色体约154 kb区间。候选基因鉴定采用池化全基因组重测序（Illumina NovaSeq平台，PE150双端测序，平均30×深度）结合ΔSNP-index分析及Sanger测序验证。功能验证采用农杆菌介导的CRISPR/Cas9基因组编辑技术，在W82背景下获得GmDBA1/GmDBA2双基因编辑纯合株系。转录组分析基于BGI测序平台完成RNA-seq，经HISAT2比对、RSEM定量和DESeq差异分析；激素检测采用超高效液相色谱-串联质谱（UPLC-MS/MS）技术；系统发育分析基于Phytozome同源蛋白序列，以MUSCLE比对和MEGA11邻接法构建进化树。样本来源包括EMS诱变的H12突变体、W82×gmdba1杂交群体以及W82背景的CRISPR编辑株系，田间试验于中国山东省青岛市（35.35° N, 119.3° E）夏季生长季进行。

四、研究结果分述

2.1 gmdba1表型鉴定：通过田间和温室对比观察，确认gmdba1在营养和生殖阶段均表现多效性矮化，株高降低超60%，节数减少、分枝数增加，V₄至V₁₀期新生叶片明显皱缩。

2.2 遗传分析与目标位点精细定位：F₂群体分离比符合单隐性基因座模型；全基因组标记筛选指向18号染色体，经两轮重组单株筛选将区间缩小至154 kb，含22个注释基因；序列分析锁定GmDBA1和GmDBA2两个NLR候选基因，二者含四个突变特异性SNP（三个非同义替换），编码蛋白均含NB-ARC和LRR结构域，且为高相似性旁系同源基因。

2.3 组织特异性表达：qRT-PCR显示GmDBA1/2在根组织中表达量约为其他组织的4倍，与其作为土传根病原侵染靶标的生物学背景相符。

2.4 gmdba1接种P6497后病害症状减轻：离体叶片和胚轴接种试验均显示gmdba1病斑面积和病原生物量显著小于野生型；扫描电镜观察发现突变体叶片表皮蜡质沉积不规则、表皮细胞形态改变。

2.5 CRISPR/Cas9双编辑株系表型：W82背景下的GmDBA1/2双敲除纯合T₂株系精确重现gmdba1的矮化、节间缩短和叶片皱缩表型，并在接种P6497后获得显著抗性增强，确认该NLR模块同时调控发育与免疫。

2.6 gmdba1中防御相关转录变化：RNA-seq鉴定出大量DEG，显著富集于植物-病原互作、MAPK信号及JA/SA激素通路；qRT-PCR确认GmPR1、GmJAR1及WRKY/MYB转录因子显著上调，表明突变驱动不依赖效应子的组成性防御开启状态。

2.7 gmdba1内源JA和SA衍生防御激素积累增强：LC-MS/MS检测显示JA、JA-Ile、JA-Val、MeJA、SA和SAG含量均显著升高，与转录组数据共同表明GmDBA1/2在基础条件下抑制JA和SA依赖的免疫输出。

五、讨论与结论

讨论部分指出，遗传定位将隐性突变表型相关位点限定于18号染色体154 kb区间内的两个串联NLR基因，CRISPR双敲除株系的表型重现提供了该位点调控大豆发育和疫霉抗性的汇聚遗传证据。但研究亦明确承认若干待解问题：双敲除策略无法区分GmDBA1和GmDBA2各自的独立贡献，也不能确定二者是否构成经典的传感器-辅助NLR配对或协同作用于汇聚信号模块；单基因敲除、互补验证及异源二聚化检测是阐明其分子层级关系的必要后续工作。此外，抗性评价仅针对P6497单一菌株，该位点是否赋予广谱抗性尚待系统验证。转录组和激素数据表明gmdba1中JA和SA信号通路持续激活，与矮化表型共同构成生长-防御权衡的生理基础，但多组学结果仅建立相关性而非因果证据。育种应用方面，GmDBA1/2位点的价值在于其赋予增强抗性的同时伴随可量化的发育代价，全功能缺失等位基因因严重自免疫代价而不适合直接农艺应用，通过等位基因挖掘、基因编辑或调控工程来解耦抗性与发育输出是更现实的目标。

研究结论部分总结认为，遗传定位和序列分析将GmDBA1和GmDBA2确定为与隐性gmdba1表型相关的18号染色体位点候选基因；原始EMS突变体和W82背景的CRISPR双编辑株系均表现出矮化及接种大豆疫霉P6497后病害症状减轻；转录组分析揭示广泛的防御相关转录重编程，并伴随突变体中JA和SA衍生代谢物积累升高。该研究将这两个NLR基因鉴定为解析大豆生长与免疫分子互作的有前景候选基因，所获得的突变遗传材料为进一步研究大豆免疫机制提供了宝贵资源，GmDBA1/GmDBA2亦为未来改良大豆疫霉抗性的分子育种提供了潜在靶基因。