编辑推荐：
本研究针对马铃薯抗孢囊线虫（PCN）的经典H1抗性位点结构不明、克隆困难的问题，通过构建携带H1位点的双单倍体材料DH4_Athlete，结合牛津纳米孔（ONT）长读长测序进行单倍型分相，并整合基于抗性基因富集测序（RenSeq）的关联遗传学（AgRenSeq-d）分析，成功重构了包含5‘和3’端重组边界在内的完整H1位点结构，为解析复杂多倍体作物抗病位点提供了创新性研究框架。
马铃薯是全球最重要的非谷类粮食作物，养活了超过十亿人口。然而，其生产可持续性正日益受到生物胁迫的挑战，其中马铃薯孢囊线虫（PCN）是最重要的威胁之一，其在全球范围内造成约9%的产量损失。更棘手的是，线虫的孢囊可在土壤中休眠长达20年以上，使其极难根除，已被全球超100个国家列为检疫性有害生物，严重限制了国际种薯贸易。面对这一“沉睡的土壤杀手”，利用植物自身的抗性基因进行育种是最经济环保的策略。其中，源自野生马铃薯材料Solanum tuberosum ssp. andigena的H1抗性位点，在欧洲成功用于防控金色马铃薯孢囊线虫（Globodera rostochiensis）Ro1和Ro4致病型已超过60年，堪称作物抗病育种中经久不衰的明星基因。然而，尽管此前研究已将该抗性定位到5号染色体末端，但由于马铃薯基因组本身四体遗传的复杂性，以及抗病基因常成簇分布且序列高度相似，H1位点的完整结构及其内的关键基因始终扑朔迷离，如同一个未被完全打开的“黑匣子”。
为了揭开H1位点的神秘面纱，研究人员在本研究中开展了一项精细的解码工程。他们并未直接强攻复杂的四倍体栽培种基因组，而是巧妙地采用了“化繁为简”的策略。首先，他们从携带H1抗性的栽培品种‘Athlete’出发，通过诱导获得了其双单倍体克隆DH4_Athlete。双单倍体基因组复杂度减半，为后续的高质量组装扫清了障碍。利用牛津纳米孔技术对DH4_Athlete进行长读长测序，研究团队获得了覆盖度高、连续性好的基因组数据，并成功构建了部分分相的两个单倍型（单倍型A和单倍型B）基因组图谱。这好比获得了一张清晰度极高的局部基因地图。与此同时，他们构建了一个包含126个马铃薯品种（58个高抗，68个高感）的关联分析群体，并利用先前在抗性品种‘Buster’上通过PacBio HiFi测序组装的抗性基因富集序列作为参考，进行SMRT-AgRenSeq-d（即结合PacBio长读长、关联遗传学和诊断性检测的RenSeq流程）分析，初步筛选出一批与H1抗性显著相关的候选NLR（核苷酸结合富亮氨酸重复蛋白，植物细胞内主要的免疫受体）基因。
本研究的关键技术方法主要包括：1. 样本制备与测序：从抗性品种‘Athlete’诱导产生双单倍体DH4_Athlete，利用牛津纳米孔技术进行长读长基因组测序与分相组装。2. 抗性基因捕获与分析：使用第五版RenSeq探针库对包括‘Buster’和关联分析群体在内的材料进行抗性基因捕获测序，并通过SMRT-AgRenSeq-d和ONT-AgRenSeq-d流程进行候选基因关联分析和验证。3. 生物信息学分析：运用nfHISS、Helixer、Resistify等生信流程进行基因组组装、基因预测、NLR鉴定、关联作图和系统发育分析。
结果
SMRT-RenSeq与关联分析群体的建立
研究人员以抗性品种‘Buster’为参考，通过PacBio HiFi RenSeq测序组装了其NLR组。同时，他们建立了一个包含126个具有明确抗/感表型的马铃薯品种的关联分析群体，为后续的关联遗传学研究奠定了基础。
SMRT-AgRenSeq-d初步筛选候选基因
通过对上述关联群体应用SMRT-AgRenSeq-d分析，在‘Buster’的组装中鉴定出89个与H1抗性显著相关的重叠群。进一步利用诊断性RenSeq计算F1分数（一种结合精确度和召回率的评估指标）进行精炼，筛选出20个F1分数高于0.95的高置信度候选NLR基因。系统发育分析显示，这些高置信度候选基因紧密聚为一支，且与已知的抗根结线虫基因Mi-1.1/Mi-1.2及抗晚疫病基因Rpi-blb2同属一个进化枝，提示了功能的潜在相关性。
ONT测序与双单倍体DH4_Athlete的组装
对DH4_Athlete的ONT测序产生了高质量的基因组数据，并成功组装出两个部分分相的单倍型基因组（单倍型A和B）。利用Helixer进行基因预测，并通过Resistify鉴定出两个单倍型中共1356个NLR基因。
ONT-AgRenSeq-d精确定位H1位点
以DH4_Athlete单倍型中鉴定出的NLR基因为参考，重新进行关联分析。结果显示，所有F1分数为1的候选基因仅存在于单倍型A中，而单倍型B中没有。结合在特定品种（如‘Strachan’、‘Vitabella’和‘Alouette’）中候选基因的存在/缺失模式，最终将保守的H1抗性渗入片段界定在单倍型A上约270千碱基的区域内。
H1位点的结构特征与比较基因组学
与不携带H1抗性的参考基因组DM v6.1及DH4_Athlete的单倍型B相比，携带抗性的单倍型A在H1位点区域显示出显著的基因扩张和重排。共线性比较和系统发育分析均支持单倍型A中该位点经历了基因重复及后续的功能分化。
讨论与结论
本研究通过整合ONT长读长测序、分相基因组组装与RenSeq关联遗传学，成功克服了马铃薯复杂基因组对H1位点研究的阻碍，首次绘制了该重要抗性位点的高分辨率物理图谱。研究不仅将抗性区间精细界定至约270kb，还明确了所有高置信度候选基因均位于DH4_Athlete的单倍型A上，清晰区分了功能性等位基因与非功能性同源序列。通过比较不同品种中已知抗性基因同源序列（如SH1-SH10）的保留与丢失情况，研究还解析了该位点历史上的重组断点，为理解其进化与育种选择提供了线索。
该研究的成功，彰显了将长读长测序、单倍型解析与关联遗传学相结合的策略，在解析多倍体作物复杂数量性状位点，特别是富含重复序列的抗病基因簇方面的强大威力。所建立的ONT-AgRenSeq-d分析框架，不依赖于感病参考基因组，能提供更准确的基因组背景和等位基因分辨率，为未来克隆H1及其他抗性基因、深入解析植物与病原互作机制，以及指导马铃薯乃至其他多倍体作物的分子设计育种，提供了强有力的方法学支撑和关键的理论依据。
