《New Phytologist》：GmSNAP14: a key contributor to soybean cyst nematode resistance in soybean
编辑推荐：
本研究发现位于14号染色体的GmSNAP14是控制大豆胞囊线虫（SCN， Heterodera glycines Ichinohe）抗性的关键数量性状位点（QTL）。携带缺失（GmSNAP14-del）或插入（GmSNAP14-ins）等位基因的抗性品种，其GmSNAP14全长转录本表达缺失或显著降低，并产生编码C端变体蛋白的剪接异构体。CRISPR/Cas9敲除GmSNAP14可增强对SCN的抗性，证实其表现为感病基因，是大豆抗性基因网络的重要补充，为通过基因编辑和育种策略培育抗性新品种提供了新靶点。
引言：囊泡运输与大豆抗线虫机制
囊泡运输在植物响应生物胁迫中扮演着关键角色，其过程被重编程以将胁迫相关蛋白和代谢物递送至特定细胞位置。可溶性N-乙基马来酰亚胺敏感因子附着蛋白受体（SNARE）复合物是驱动膜融合的核心机器，而该复合物的拆解由N-乙基马来酰亚胺敏感因子（NSF）和可溶性NSF附着蛋白（SNAP）介导，对维持囊泡运输的连续性至关重要。在大豆中，α-SNAP蛋白已成为赋予对主要病原体大豆胞囊线虫（SCN）抗性的关键角色。大豆基因组包含五个注释的α-SNAP基因，其中三个（GmSNAP18、GmSNAP11和GmSNAP02）已被证明在赋予SCN抗性中发挥作用。本研究的亲本材料“北京”（Peking）是SCN抗性基因的主要来源之一，而“皮克特”（Pickett）则是通过“北京”基因渗入培育出的第一代抗性品种。两者携带已知的主要抗性基因（GmSNAP18-a、GmSNAP11-sp、GmSHMT08和GmNSFRAN07）的相同等位基因，但对SCN HG类型1.2.5.7（小种2）分别表现出中度抗性和中度感病性。这提示“北京”可能携带额外的抗性基因，而本研究正是通过构建“皮克特”ד北京”的分离群体来寻找这个基因。
材料与方法
研究构建了一个由168个F4:5家系组成的重组自交系（RIL）群体，亲本为“皮克特”和“北京”。使用SCN近交系TN22（HG类型1.2.5.7）对群体进行抗性表型鉴定。利用Illumina Infinium BARCSoySNP3K BeadChip对群体进行基因分型，并使用MapQTL 6.0软件进行连锁作图和数量性状位点（QTL）定位。通过全基因组重测序（WGRS）和PacBio HiFi长读长测序技术，鉴定了GmSNAP14基因的等位基因变异。设计了TaqMan检测方法用于高通量基因分型。通过实时定量逆转录聚合酶链反应（qRT-PCR）分析了GmSNAP14在不同基因型及线虫侵染后的表达模式。利用CRISPR/Cas9基因编辑技术，在“福里斯特”（Forrest）及其Gmshmt08突变体F234的背景下，构建了针对GmSNAP02、GmSNAP14以及两者同时敲除的复合植株，并对转基因毛状根进行SCN接种表型分析，以验证基因功能。
结果：抗性QTL定位与GmSNAP14等位基因的发现
对“皮克特”ד北京”群体的表型评估显示，该群体存在超亲分离现象，超过30%的家系表现出比中度抗性亲本“北京”更低的雌性指数（FI）。连锁作图鉴定出一个主效QTL14和一个微效QTL12。主效QTL14定位于14号染色体，解释了62.4%的表型变异，其加性效应表明抗性等位基因来源于“北京”。在该QTL的置信区间内，含有GmSNAP14基因（Glyma.14G054900）。由于GmSNAP14是已知抗性基因GmSNAP18、GmSNAP11和GmSNAP02的旁系同源物，因此被列为优先候选基因进行深入分析。
通过对SCN抗性大豆种质“北京”、PI 437654和PI 90763进行全基因组重测序和PacBio HiFi测序分析，发现了GmSNAP14的两个变异等位基因。在“北京”中鉴定出一个约400 bp的缺失等位基因，命名为GmSNAP14-del，该缺失覆盖了部分内含子6、外显子7、内含子7、外显子8和部分内含子8。在PI 90763和PI 437654中，发现了一个约9725 bp的大片段插入等位基因，命名为GmSNAP14-ins，该插入位于第三个内含子中，经比对与copia/Ty1类逆转录转座子家族高度相似。研究为此开发了用于高通量检测这些等位基因的TaqMan检测方法。
GmSNAP14的表达模式与转录本分析
转录本分析揭示了不同等位基因背景下的表达差异。在携带野生型等位基因（GmSNAP14-WT）的感病品种“威廉斯82”（Williams 82）及中度感病品种“福里斯特”、“皮克特”中，主要检测到全长的GmSNAP14转录本（T1， 870 bp）。而在携带GmSNAP14-del等位基因的“北京”中，主要检测到缺失外显子7和8的T2转录本（741 bp）以及缺失外显子6、7、8的T3转录本（581 bp），全长T1转录本缺失。携带GmSNAP14-ins等位基因的PI 90763和PI 437654，GmSNAP14的总体表达水平显著低于其他基因型，并且能够检测到少量T1、T2和T3转录本的混合。这些剪接变异体的产生与基因型相关，但与SCN侵染无关。GmSNAP14 T2和T3转录本分别编码246和174个氨基酸的蛋白质变体。T2变体保留了全部四个TPR结构域，而T3变体由于移码导致第四个TPR结构域被截短。这些数据表明，SCN抗性可能源于GmSNAP14表达的减弱以及功能丧失的GmSNAP14蛋白变体的组合效应。
启动子分析揭示了顺式作用元件的差异
比较“北京”、PI 437654和“威廉斯82”的GmSNAP14启动子序列，发现了多个单核苷酸多态性（SNP）、插入和缺失。特别值得注意的是，在“威廉斯82”和“北京”的启动子中存在HD-ZIP转录因子结合元件，而在表达水平较低的PI 437654中该元件缺失，这提示HD-ZIP元件可能在调控GmSNAP14基因表达中扮演潜在角色。
基因编辑功能验证：GmSNAP14是SCN的感病基因
为了直接验证GmSNAP14功能缺失是否贡献于SCN抗性，研究利用CRISPR/Cas9技术在“福里斯特”背景下敲除了GmSNAP14。与转入空载体对照的植株相比，GmSNAP14敲除的复合植株根系上的胞囊数量显著减少。类似地，单独敲除GmSNAP02也能增强抗性。而同时敲除GmSNAP02和GmSNAP14，在部分实验中显示出比单独敲除任一基因更强的抗性趋势。这些结果直接证明，靶向敲除GmSNAP14能增强大豆对SCN HG类型1.2.5.7的抗性，从而确认GmSNAP14是一个感病基因。
GmSNAP02/14介导的抗性独立于GmSHMT08
为了探究GmSHMT08是否调控GmSNAP02和GmSNAP14介导的抗性，研究在“福里斯特”的Gmshmt08功能缺失突变体F234中进行了相同的基因编辑实验。结果显示，在F234背景中敲除GmSNAP02、GmSNAP14或两者，同样能显著降低胞囊数量，其抗性增强效果与在“福里斯特”野生型背景下一致。这表明GmSNAP02和GmSNAP14介导的对HG类型1.2.5.7的抗性，独立于GmSHMT08的功能。然而，在缺乏其他已知抗性基因（如GmSNAP18-a和GmSNAP11-sp）的完全感病品种“埃塞克斯”（Essex）中，同时敲除GmSNAP02和GmSNAP14并不能显著降低胞囊数，说明GmSNAP02/14的功能缺失所提供的增强抗性，需要GmSNAP18-a和GmSNAP11-sp等位基因的共存。
讨论与结论：一个扩展的SCN抗性基因网络模型
本研究确凿地证明了GmSNAP14-del/ins等位基因在赋予大豆对SCN抗性中的作用。与之前发现的通过大片段插入/缺失导致转录本完全缺失或严重截短的GmSNAP02等位基因不同，GmSNAP14的等位基因通过独特的机制导致功能丧失：GmSNAP14-del产生框内缺失，生成更小的蛋白质异构体；而GmSNAP14-ins则可能通过大型内含子插入干扰转录或RNA加工，从而显著降低基因的基础表达水平。这些机制共同导致GmSNAP14功能受损或缺失。
基于这些发现，研究提出了一个扩展的SCN与大豆互作模型。SCN通过口针分泌效应子，靶向宿主细胞正常的α-SNAP蛋白，以促进囊泡运输，从而成功建立取食细胞。而抗性大豆则通过进化出多种策略破坏这一过程：GmSNAP18-a编码非典型α-SNAP变体，能“诱骗”效应子但无法与典型NSF结合；GmSNAP11-sp产生截短蛋白，可能无法积累，从而完全“规避”效应子；GmSNAP02和GmSNAP14则通过等位基因变异导致功能丧失或表达降低，进一步削弱了线虫效应子可利用的宿主靶标。在“北京”中，GmSNAP14-del变体与GmSNAP11-sp、GmSNAP18-a共同作用，导致囊泡运输严重受损，从而表现出对特定SCN种群的中度抗性。而在PI 90763/PI 437654中，GmSNAP02功能完全丧失，结合GmSNAP14的低水平表达，提供了更全面的抗性。
总之，GmSNAP14是“北京”抗性遗传架构中的关键组成部分，与GmSNAP18-a、GmSNAP11-sp、GmSHMT08和GmNSFRAN07共同构成了一个复杂的抗性基因网络。该研究不仅深化了对大豆SCN抗性分子机制的理解，而且为大豆抗性育种提供了新的工具和策略。通过使用本研究开发的检测方法，可以将GmSNAP14-del或GmSNAP14-ins等位基因与现有抗性基因堆叠，培育具有更广谱和更持久抗性的大豆新品种。利用基因编辑技术靶向修饰GmSNAP14，也为未来设计新型抗性提供了可能。然而，在应用时需考虑基因功能冗余及对植株正常生长发育的潜在影响，采取精准的编辑策略。
