一个转座子的跳跃，给猕猴桃装上了抗病武器
一个转座子的跳跃，给猕猴桃装上了抗病武器
猕猴桃溃疡病是全球猕猴桃产业面临的最严重威胁之一。这种由丁香假单胞菌猕猴桃致病变种（Pseudomonas syringae pv. actinidiae，简称 Psa）引起的细菌性溃疡病，可以在几年内摧毁整片果园。自2008年以来，Psa 在意大利、新西兰、韩国、日本等主要猕猴桃产区引发了大规模疫情，造成了数十亿美元的经济损失。
栽培猕猴桃对 Psa 的抗性普遍较弱。但猕猴桃的野生近缘种——特别是毛花猕猴桃（Actinidia eriantha）——对 Psa 表现出明显更强的抗性。
问题是：这种抗性是怎么来的？
2026年4月发表在 Genome Biology 上的一项研究给出了一个有趣的答案：一个转座子的”跳跃”，给抗病基因装上了一个新的开关。
泛基因组：不只看一个参考基因组
传统基因组研究通常依赖一个”参考基因组”——选择一个代表性个体进行测序，作为后续分析的基准。但这种方法有一个根本性的局限：任何一个个体的基因组都无法代表整个物种的遗传多样性。
泛基因组（pangenome）方法试图弥补这一不足。它对一个物种的多个个体进行基因组测序和组装，然后把这些基因组整合成一个”图谱”，记录所有个体之间的共同序列和变异序列。
这项研究构建了猕猴桃的泛基因组。研究团队对11个野生毛花猕猴桃材料和1个种间杂交种进行了染色色体级别的基因组组装，加上之前已经发布的12个基因组（包括3个毛花猕猴桃和9个中华猕猴桃），构建了一个”参考无偏”的图结构泛基因组。
这个泛基因组揭示了什么？ 3100多万个 SNP（单核苷酸多态性）、1300多万个 InDel（插入缺失）、以及62万多个结构变异（structural variations）。这是一个巨大的遗传变异图谱，记录了栽培猕猴桃和野生猕猴桃之间的基因组差异。
结构变异：比 SNP 更大的遗传差异
在这62万多个结构变异中，有一种特别值得关注：由转座子插入导致的结构变异。
结构变异（SV）是指基因组中较大规模的序列变化——通常定义为长度超过50个碱基对的插入、缺失、倒位或易位。相比只有单个碱基变化的 SNP，结构变异对基因功能的影响往往更大——它们可以打断基因、改变基因的调控环境、或者引入全新的序列。
转座子插入是结构变异的重要来源之一。当一个活跃的转座子在基因组中跳跃时，它可以在新位置插入一段几百到几千个碱基对的序列。这段插入序列本身可能携带调控信息——比如转录因子结合位点。
这正是这项研究中发现的情况。
NLR25-1：一个抗病基因的”升级”
研究团队通过泛基因组分析，在野生毛花猕猴桃中发现了一个对 Psa 抗性有显著贡献的基因：AeNLR25-1。
NLR 基因（Nucleotide-binding Leucine-rich Repeat genes）编码的蛋白质是植物先天免疫系统的核心组件。它们可以识别病原菌释放的效应蛋白，触发植物的免疫反应——这种反应通常以局部细胞死亡（过敏性坏死反应）的形式表现，限制病原菌的扩散。
关键的发现是：AeNLR25-1 在野生和栽培猕猴桃之间存在等位基因特异性表达（allele-specific expression）差异。也就是说，来自野生亲本的等位基因比来自栽培亲本的等位基因表达水平更高。
为什么同一个基因在不同个体中的表达水平不同？
答案在启动子区域。研究发现，在野生猕猴桃的 AeNLR25-1 启动子区域，有一段转座子插入。这段转座子插入引入了一个新的 WRKY 转录因子结合位点。
WRKY 转录因子是植物免疫反应中的关键调控因子。当病原菌侵染时，WRKY 转录因子被激活，结合到目标基因的启动子上，增强这些基因的表达。在野生猕猴桃中，转座子插入带来的新 WRKY 结合位点，让 AeNLR25-1 在病原菌侵染时获得了更强的表达上调能力，从而产生了更强的免疫反应。
简单说：一个转座子的跳跃，给抗病基因装上了一个更强的”加速器”。
这不是”演化在设计”
需要明确的是：转座子插入产生新的 WRKY 结合位点，是一个随机事件，不是”演化在主动设计”抗性。
转座子在基因组中的跳跃是相对随机的——它们不”知道”应该插入到哪里对宿主有利。在大多数情况下，转座子插入可能是有害的（打断功能基因）或者中性的（插入到基因之间的非功能区域）。但在极少数情况下，转座子插入恰好落在了有利的位置，比如在一个免疫基因的启动子区域，引入了一个新的调控元件。
自然选择在这个过程中起作用的方式是：拥有这个有利转座子插入的个体在 Psa 流行的环境中更容易存活和繁殖，导致这个等位基因在种群中的频率逐渐升高。但选择筛选的是结果，不是原因——转座子的跳跃本身是随机的，选择只是保留了碰巧有利的那些插入。
泛基因组在作物改良中的价值
这项研究的意义不只是揭示了一个有趣的抗病机制，它还展示了泛基因组方法在作物改良中的实际价值。
传统育种在利用野生近缘种的有利基因时，通常依赖杂交和回交，通过表型筛选来选择具有目标性状的后代。这种方法有效，但效率有限——特别是当目标性状由多个基因控制时，筛选过程可能需要很多代。
泛基因组方法提供了一种更精确的工具：通过比较栽培种和野生种的基因组，直接定位产生表型差异的具体变异——就像这项研究中定位到 NLR25-1 启动子的转座子插入一样。这种精确的定位可以帮助育种者设计更高效的分子标记辅助选择策略，甚至直接利用基因编辑技术将有利变异引入栽培品种。
但需要指出的是：从基因发现到实际的品种改良还有很长的路。 将一个野生等位基因引入栽培品种需要考虑连锁累赘（不利基因随有利基因一起被引入）、基因表达在不同遗传背景中的稳定性、以及田间条件下的实际抗性表现。基因组数据提供了线索，但最终的验证仍然需要田间实验。
生信在这项工作中的角色
这项研究高度依赖生物信息学方法，贯穿了从数据产生到结论得出的整个流程：
基因组组装。 将测序产生的短读段组装成完整的染色体级别基因组，需要专门的组装算法和质量评估工具。
泛基因组构建。 将24个基因组整合成一个图结构泛基因组，需要检测所有个体之间的 SNP、InDel 和结构变异，然后以图的形式表示这些变异。
结构变异检测。 从基因组比对中识别大尺度的结构变异（包括转座子插入），比 SNP 检测更困难——因为结构变异往往发生在重复序列附近，比对的歧义性更大。
等位基因特异性表达分析。 判断来自父本和母本的等位基因是否表达水平不同，需要将 RNA-seq 数据精确地映射到包含已知变异的基因组上，然后分别统计两个等位基因的表达量。
顺式调控元件预测。 在转座子插入序列中识别 WRKY 转录因子结合位点，需要结合 motif 分析和转录因子结合数据。
每一步分析都需要考虑数据的质量、算法的参数选择、以及结果的生物学解释。这些分析的可靠性直接影响了研究结论的可信度。
回头看
这项研究是一个很好的例子，展示了现代生物学中”计算分析”和”实验验证”如何协同工作。
泛基因组分析发现了 NLR25-1 这个候选基因——这是计算分析提供的线索。遗传和分子实验验证了这个候选基因确实与抗性相关，并且揭示了转座子插入引入 WRKY 结合位点的机制——这是实验验证确认的因果关系。
如果只有计算分析，我们可能得到一个”看起来相关”的候选基因，但无法确定它是否真的有因果作用。如果只有实验验证，我们可能需要逐一测试成百上千个候选基因，效率极低。
计算提出假设，实验验证假设——这个循环是现代生物研究的核心模式。
对于非专业读者来说，这个故事也提供了一个有趣的视角：你吃的猕猴桃，它的抗病能力可能来自一段数百万年前”跳”进基因组的转座子。基因组不是一成不变的——它是一个动态演化的系统，充满了历史的痕迹和意外的创新。
参考文献
Wu, Y., Wang, Y., Lin, Y., et al. (2026). Allelic structural variation at the NLR25-1 locus enhances defense against Pseudomonas syringae in kiwifruit. Genome Biology. https://doi.org/10.1186/s13059-026-04068-0
Vanneste, J. L. (2017). The scientific, economic, and social impacts of the New Zealand outbreak of bacterial canker of kiwifruit (Pseudomonas syringae pv. actinidiae). Annual Review of Phytopathology, 55, 377-399.
Li, J., et al. (2022). Plant NLR immune receptors: from guarded to guard. Annual Review of Plant Biology, 73, 191-221.
Bourque, G., et al. (2018). Ten things you should know about transposable elements. Genome Biology, 19, 199.