在马铃薯种植区，一种名为“早疫病”的病害正悄然蔓延，由链格孢菌（Alternaria）引起，它能导致叶片出现同心圆状病斑，严重时造成20-50%的产量损失，甚至绝收。尽管其危害性常被晚疫病的“光环”所掩盖，但在北美、北欧和亚洲部分冬季种植区，早疫病已跃升为主要生产限制因素。然而，与病原菌致病机理研究的热度相比，马铃薯自身抗早疫病（Early Blight Resistance, EBR）基因的挖掘却进展迟缓。过去几十年间，仅有零星几个耐受基因型被报道，明确功能的EBR基因屈指可数。这种“知己不知彼”的困境，使得培育抗病品种缺乏关键的“基因武器”，生产上过度依赖化学药剂，不仅增加成本和环境压力，更威胁全球粮食安全与农业可持续发展。因此，系统鉴定并验证马铃薯的EBR基因，成为一项紧迫而重要的科学任务。
在此背景下，一项发表于《Horticulturae》的研究为我们带来了突破。为了攻克这一难题，研究团队设计了一项精巧的实验。他们选取了地理来源遥远、抗病性迥异的两个马铃薯栽培品种——来自美国的相对抗病品种‘Atlantic’（Atl）和来自中国的感病品种‘Jinshu 16’（J16）——进行杂交，构建了一个F1分离群体。这个策略巧妙地利用了栽培种内部的遗传多样性，旨在打破四倍体基因组中不利的性状连锁，创造新的抗病基因组合。研究人员综合运用了抗性表型精准评价、组织病理学分析、混池分组测序（Bulked Segregant Analysis sequencing, BSA-seq）、转录组测序（RNA sequencing, RNA-seq）以及分子动力学模拟（Molecular dynamics simulation）等多学科技术手段，形成了一套强大的“组合拳”。
本研究主要采用了以下几项关键技术方法：首先，利用离体叶片接种和整株调查对F1群体进行抗性表型鉴定，并基于极端表型（各21株）构建抗、感混池用于BSA-seq。其次，对亲本Atl和J16在接种病原菌前后的叶片进行RNA-seq，分析差异表达基因。接着，通过BSA-seq结合六种算法（SNP-index, ED, GPS, QTGseq, Gprime, OcBSA）定位抗性相关数量性状位点（QTL），并与转录组数据整合筛选候选基因。然后，对候选基因进行突变位点分析和分子动力学模拟，预测蛋白质结构变化。最后，通过农杆菌介导的瞬时过表达技术，在马铃薯叶片中验证候选基因的功能。
3.1. 病原菌分离与鉴定
从马铃薯病叶中分离得到的病原菌，经形态学和分子系统发育分析，被鉴定为细极链格孢（Alternaria tenuissima），并完成了柯赫氏法则验证，确认为本研究的致病菌。
3.2. 抗、感品种的病害进程与细胞学防御表征
研究发现，感病品种J16在接种后24小时内即可观察到病原菌定殖，病斑迅速扩大；而抗病品种Atl能有效抑制病斑扩展长达72小时，最终病斑直径显著小于J16。细胞学观察进一步揭示，Atl能快速激活局部过敏反应和协调的胼胝质沉积，成功限制病原菌扩展；而J16则表现为大面积组织坏死和极少的胼胝质沉积，防御机制失效。
3.3. 抗、感品种的差异表达基因
RNA-seq分析显示，接种病原菌后，抗病品种Atl中显著上调的基因在几丁质酶活性、谷胱甘肽代谢过程、对真菌的防御反应等通路中富集，表明其防御反应更具协调性和特异性。值得注意的是，即使在健康条件下，Atl在植物激素信号转导、植物-病原互作等防御相关通路上也已呈现出更高的基础表达水平，体现出更强的“基础防御准备状态”。
3.4. BSA-seq鉴定候选区域
通过对极端抗、感混池的BSA-seq分析，并综合六种算法的结果，研究将主要的抗性QTL定位在马铃薯第8号染色体的三个重叠区域：25.07–29.20 Mb， 38.05–38.80 Mb， 以及 39.40–40.78 Mb。
3.5. 整合BSA-seq与RNA-seq鉴定CND41为马铃薯EBR关键候选基因
在定位区间内，通过整合表达和突变分析，天冬氨酸蛋白酶基因CND41（PGSC0003DMG400023144）被确定为核心候选基因。该基因在抗病亲本Atl和抗病混池中，存在六个位于TAXI-N结构域的非同义突变（G213R, I227V, I283L, Q261K, R302G, V251G）。表达分析表明，CND41在Atl中具有更高的基础表达水平，且在病原菌诱导后上调幅度更大。
3.6. 野生型与突变体蛋白的三维结构预测与分子动力学模拟
分子动力学模拟揭示了单点突变与六点联合突变导致抗性的不同机制。大多数单点突变（如G213R, Q261K）主要通过破坏蛋白内部关键相互作用，使整体结构变得松散、稳定性下降，从而破坏了原有的感病相关构象。而六个突变同时存在时，则协同作用，使蛋白整体结构变得更加紧凑、稳定，氢键网络增强，两个关键区域（100-125和400-425位氨基酸）的柔性显著降低，可能形成了一个新的、稳定的构象，从而阻断病原菌效应子结合或改变下游信号。
3.7. 马铃薯中CND41基因的瞬时过表达分析
功能验证表明，在感病品种J16中瞬时过表达CND41能显著减小病原菌接种后的病斑直径。在本身已具抗性的Atl品种（携带Q261K单点突变）中过表达该基因，也能观察到病斑减小的趋势。这为CND41正向调控马铃薯早疫病抗性提供了功能证据。
结论与讨论 本研究成功鉴定出天冬氨酸蛋白酶基因CND41是调控马铃薯早疫病抗性的一个关键候选基因。其重要性体现在多个层面：遗传上，该基因所在的染色体8区域被多个BSA算法共同定位；序列上，其TAXI-N结构域存在多个与抗性表型共分离的错义突变；表达上，在抗病基因型中呈现更强的组成型和诱导型表达；结构上，分子模拟预测其突变能根本性改变蛋白构象稳定性；功能上，瞬时过表达能增强植物的抗病性。CND41所属的TAXI蛋白家族已知具有木聚糖酶抑制活性，这可能通过抑制病原菌分泌的细胞壁降解酶或保存寡糖激发子来增强抗性。该研究不仅为理解四倍体马铃薯复杂数量抗性的遗传基础提供了范例，证明了整合多组学与计算生物学在挖掘作物抗病基因中的强大效力，更重要的是为马铃薯抗早疫病分子育种提供了一个经过多维度证据支持的强有力候选基因。未来通过稳定遗传转化或基因编辑对该基因进行功能确证，并将其与其它抗性位点聚合，有望培育出具有持久、广谱抗性的马铃薯新品种，从而减少对化学药剂的依赖，推动马铃薯产业的绿色可持续发展。
