Take home message

（1）在植物叶片上，特定的有益细菌（共生假单胞菌）能够“训练”或“激活”植物自身的防御系统，从而选择性地抑制某些病原细菌（致病假单胞菌）的生长，保护植物健康。但这个过程的效果，取决于植物和细菌双方具体的基因型。

（2）本研究以拟南芥叶片为场景，研究了其中自然存在的共生假单胞菌群体与致病假单胞菌群体之间的相互作用。研究发现，有益菌（共生假单胞菌）并不会直接杀死病原菌。它们的作用方式是先与植物互动，触发植物自身的免疫或防御反应。

（3）植物被“激活”后，会产生一种选择性的抑制环境，能够特异性地针对和抑制某个特定的致病菌谱系，从而达到保护植物的目的。这种保护作用的强弱因植物基因型不同而异。这说明保护效果不是由细菌单方面决定的，而是由植物介导的。

（4）研究中也发现有个别致病假单胞菌菌株能够“逃脱”或“规避”由共生菌引发的这种植物保护。这再次证明了微生物菌株间微小的遗传差异会导致巨大的功能差异。

摘要

（1）The community structure in the plant-associated microbiome depends collectively on host–microbe, microbe–microbe and host–microbe–microbe interactions. 植物相关微生物群落的群落结构共同依赖于宿主-微生物、微生物-微生物和宿主-微生物-微生物的相互作用。 

（2） The ensemble of interactions between the host and microbial consortia may lead to outcomes that are not easily predicted from pairwise interactions 宿主和微生物菌群之间的整体相互作用可能导致不容易从成对相互作用中预测的结果 .

（3）Plant–microbe–microbe interactions are important to plant health but could depend on both host and microbe strain variation. 植物-微生物-微生物之间的相互作用对植物健康很重要，但可能取决于宿主和微生物菌株的变异。

（4）Here we study interactions between groups of naturally co-existing commensal and pathogenic Pseudomonas strains in the Arabidopsis thaliana phyllosphere 在这里，我们研究了拟南芥叶层中自然共存的共生和致病假单胞菌菌株之间的相互作用 .

（5）We find that commensal Pseudomonas prompt a host response that leads to selective inhibition of a specific pathogenic lineage, resulting in plant protection. 我们发现共生假单胞菌引起宿主反应，导致特定病原谱系的选择性抑制，从而产生植物保护。

（6）The extent of protection depends on plant genotype, supporting that these effects are host-mediated 保护的程度取决于植物基因型，支持这些效应是宿主介导的 .

（7）Strain-specific effects are also demonstrated by one individual Pseudomonas isolate eluding the plant protection provided by commensals. Our work highlights how within-species genetic differences in both hosts and microbes can affect host–microbe–microbe dynamics. 一个单独的假单胞菌分离株逃避了共生体提供的植物保护，这也证明了菌株特异性效应。我们的工作强调了宿主和微生物的种内遗传差异如何影响宿主-微生物-微生物动力学。

（8）Result

（9）Barcoding of Pseudomonas isolates and experimental design 假单胞菌分离株的条形码及实验设计

子图a  拟南芥和假单胞菌采样点。这些点有两类采集目标：拟南芥：即本研究的宿主植物。假单胞菌：即从这些地方的拟南芥叶片上分离得到的细菌。子图b， 这是一个系统发育树加上条形图的组合图。左侧 系统发育树：这是一棵基于细菌核心基因组序列构建的“家族树”。树干与分支：代表了14株本研究所使用的假单胞菌分离株之间的亲缘关系。树枝长短：下方标尺（0.1 nucleotide substitutions per site）表明，树枝长度代表了遗传距离（即序列差异的程度）。树枝越长，表示该菌株与其它菌株的遗传差异越大。右侧 条形图：数据来源：来自一个名为 “Karasov collection” 的、包含1524个假单胞菌分离株的大型公共菌种库。条形含义：每个条形对应左边系统发育树上的一个ATUE组（即一个特定的OTU）。条形的长度表示，在Karasov收藏集的1524个分离株中，与本研究使用的该组菌株亲缘关系极近的菌株的数量。“亲缘关系极近”的定义：在核心基因组上，序列差异小于 0.0001。这几乎是同一个克隆或非常近缘的菌株。本论文精心挑选的这14个实验菌株，并非稀有怪菌。在右侧条形图中，像ATUE5（P. viridiflava）这样的组别拥有很长的条形，意味着在自然界采集的1524个菌株中，有大量（可能上百个）与它几乎相同的菌株。这说明本研究关注的细菌互作现象，很可能在自然界中广泛存在且频繁发生。

（10） Host effects on composition of synthetic communities 宿主对合成群落组成的影响

（11） Host-dependent pathogenicity, growth promotion or protection 依赖宿主的致病性、生长促进或保护

 有益菌（共生菌）提供的保护作用因拟南芥品种不同而异，这直接证明了保护是由植物介导的，而不是细菌之间的直接对抗。顶部面板：原始数据与直观趋势。Y轴（左侧）： 鲜重。这是衡量植物生长的关键指标，鲜重越高代表植物长得越好，受病害影响越小。X轴（底部）： 六种不同的拟南芥生态型/基因型。CommonCom (共生菌组合)： 接种了多种共生假单胞菌的混合群体。这是为了测试有益菌群本身是否对植物有直接影响（如促进生长）。MixedCom (混合菌组合)： 接种了PathoCom和CommonCom的混合群体。这是实验的核心处理，测试有益菌群能否在病原菌存在时保护植物。

底部面板：效应大小与统计推断，Y轴（左侧）： 与对照组相比的平均差异。这是一个更精确的统计呈现，它量化了“保护效果”到底有多大。

（12） Load-dependent impact of pathogens and commensals 病原体和共生体的负荷依赖性影响

这张图旨在回答两个关键问题：不同处理下，植物叶片上的细菌总载量是多少？细菌载量与植物的健康状况（鲜重）有什么关系？子图a：细菌载量的分布，标题： 三种合成菌群的细菌载量（以10为底的对数）密度图。X轴： log10(bacterial load)，即细菌载量（取以10为底的对数）。这是为了将可能跨越几个数量级的数据压缩到更直观的尺度上。Y轴： Density（密度），表示在某个载量值附近数据点的相对集中程度。PathoCom（红色）： 仅病原菌群。CommenCom（绿色）： 仅共生菌群。MixedCom（蓝色）： 病原菌与共生菌的混合菌群。子图b：细菌载量与植物健康的相关性，X轴： log10(bacterial load)，Y轴： Rosette fresh weight (mg)，植物鲜重，是衡量植物生长和健康的指标。同样用红、绿、蓝三色区分PathoCom、CommenCom、MixedCom三种处理。PathoCom（红色）： 预期会有强烈的负相关。即病原菌载量越高，植物长得越差（鲜重越低）。改变病原菌的致病状态或增强植物耐受性，使得植物即使在有一定病原菌负载的情况下也能健康生长

（13） Protection by commensal members and pathogen suppression 共生成员的保护和病原体抑制

共生菌在植物体内（in planta）对病原菌的抑制模式，是否与在培养皿上（in vitro）的直接对抗模式一致？图 a：微生物群落整体结构变化（NMDS分析），直观展示三种处理（PathoCom， CommenCom， MixedCom）下，植物叶片上整个细菌群落的组成是否不同。PathoCom（红点）和CommenCom（绿点） 的样本点清晰地聚集在不同的区域，这说明仅病原菌和仅共生菌的群落结构截然不同，符合预期。MixedCom（蓝点） 的样本点并没有与PathoCom（红点）重叠，而是形成了一个独立的、居中的集群。这说明，当共生菌与病原菌共存时，形成的群落既不是简单的病原群落，也不是共生群落，而是一个全新的、独特的群落状态。这强烈暗示有更高层级的调控（即植物宿主）在影响群落的组装。子图 b & e：特定病原菌在植物体内的丰度变化，精确量化在MixedCom处理下，每一个单独的细菌菌株（特别是病原菌）的丰度，相比其单独存在时是增加还是减少了。P1, P2, P6, P7 这几个病原菌（都属于ATUE5谱系，紫色点）的可信区间完全位于负值区域（不跨越0）。这意味着在共生菌存在的情况下，它们在植物体内的数量被显著、特异性地抑制了。而其他非ATUE5的病原菌（如P3， P4， P5）没有被显著抑制。共生菌触发的保护，不是“无差别攻击”，而是植物发起的一场“精准打击”，主要针对ATUE5谱系的病原菌。子图 c：培养皿上的直接抑制网络（in vitro），展示在没有植物参与的情况下，这些细菌菌株两两之间直接的抑制关系。箭头从“抑制菌”指向“敏感菌”。在平板上，共生菌（C系列）对病原菌（P系列，包括ATUE5和非ATUE5）存在大量、广泛的直接抑制作用（黄色箭头很多）。这与图b/e在植物体内的选择性抑制形成鲜明对比。子图 d：植物体内菌群的相关性网络（in planta），揭示在真实的植物环境中，不同菌株之间的共生或竞争关系。连线表示菌株间的丰度相关性（正相关或负相关），线越粗、颜色越深（紫色为正，黄色为负），相关性越强。强烈的负相关（竞争）：主要发生在ATUE5病原菌内部（紫色节点之间），以及它们与部分共生菌之间。强烈的正相关（共生）：主要发生在非ATUE5的共生菌之间（绿色节点之间）。在植物这个特殊环境中，细菌之间的互作关系被重塑了。最强的“战争”不是发生在共生菌 vs. 所有病原菌之间，而是植物环境加剧了ATUE5病原菌内部的竞争，同时促进了有益菌之间的合作。这再次支持了植物环境是关键调节者。

（14） Defensive response elicited by non-ATUE5 非ATUE5引起的防御反应

从 “表型观察” 深入到 “分子机制” 的关键一步，它直接揭示了植物在转录组水平对不同微生物处理的响应。特异性激活：共生菌（CommenCom）和混合菌群（MixedCom）都能强烈且独特地激活植物的防御相关基因表达程序。叠加与协同：混合处理（MixedCom）激活的基因，不是病原菌处理和共生菌处理的简单叠加。它包含大量独有的防御响应基因，这对应于植物被“训练”后产生的增强的、准备就绪的免疫状态。分子验证：关键的防御标志基因在共生菌和混合菌处理下被显著诱导，从分子层面证实了宿主免疫系统被激活。

子图 a：差异表达基因的整体模式，目的：直观展示所有处理中，植物基因表达变化的全局图谱。这是一张热图。每一行代表一个差异表达基因，共2727个。每一列代表一个处理组与对照组的比较（如 Com vs Ctrl， P vs Ctrl， Mix vs Ctrl）。颜色表示相对表达水平：红色代表上调（表达量高于对照），蓝色代表下调（表达量低于对照）。基因根据表达模式的相似性进行了聚类（纵向分组）。共生菌和混合处理激活了一套相似但又有所区别的、大规模的基因表达程序。子图 b：差异表达基因的韦恩图，目的：精确量化并比较不同处理所特异性激活的基因集合。三个椭圆分别代表由 PathoCom、CommenCom、MixedCom 处理引起的差异表达基因集合。重叠部分代表两种或三种处理共同激活的基因。CommenCom 和 MixedCom 共享核心：两者有189个共同上调的基因。这很可能是共生菌触发的 “基础防御预备程序”。MixedCom 的特有响应：MixedCom 处理独有1370个上调基因，这个数量非常庞大。这说明，当植物在共生菌“预警”下遭遇病原菌时，它启动了一套更强大、更特异的防御反应程序，而不仅仅是把两套反应加起来。PathoCom 的有限响应：病原菌单独处理激活的基因最少，且大部分（409个）是其独有的，与共生菌诱导的程序重叠很少。这可能是一种被动的、滞后的病害反应，而非有效的主动防御。子图 c：基因功能（GO富集分析）：揭示上述不同基因集合主要参与哪些生物学过程，从而推断植物的生理状态。列出了三个关键基因集（b图中的189交集、630独有、1370独有）中最显著富集的前10个非冗余GO条目。GO条目是对基因功能的描述。CommenCom & MixedCom 交集： 很可能富集与 “免疫信号感知”、“早期防御信号转导” 相关的功能。这是共生菌触发的初始警报状态。CommenCom 独有： 可能富集与 “生长发育调节”、“基础代谢调整” 等相关的功能，反映了植物与共生菌的友好互动。MixedCom 独有（最重要）： 这里预期会强烈富集 “系统性获得性抗性”、“水杨酸介导的信号通路”、“病程相关蛋白活性”、“细胞壁加固”、“氧化爆发” 等典型且强烈的免疫防卫功能。这从功能上证实，混合处理将植物推入了一种全面备战的高免疫状态。子图 d：关键防御标志基因的表达验证，通过几个已知的、经典的防御基因，来直接验证并量化免疫激活的程度。选取了6个已知的防御标志基因（如 PR1， PDF1.2 等）。用柱状图展示它们在 Control、PathoCom、CommenCom、MixedCom 处理下的实际表达量。使用字母标注进行统计学差异比较（标有相同字母的组间无显著差异）。MixedCom： 关键防御基因的表达量最高，显著高于其他所有处理。这提供了分子水平的直接证据：正是共生菌与病原菌共存时，植物的防御系统才被最大化地激活，从而解释了表型上的保护效果。

（15） Lack of protection explained by a single pathogenic isolate 用单一病原分离株解释缺乏保护

通过机器学习、统计学和基因工程验证，精准地找到了整个保护效应的关键“开关”菌株， 在混合菌群（MixedCom）中，究竟是哪一个（或哪几个）具体的细菌菌株，对预测和保护植物健康起着最关键的作用，核心发现：关键病原菌株 P6 是决定植物在 MixedCom 处理下能否健康生长的最强预测因子和必要条件。P6 的丰度与植物鲜重呈强烈的负相关：P6 越多，植物长得越差。共生菌群提供的保护作用，其核心机制之一就是特异性地强烈抑制 P6 的生长。如果从混合菌群中移除 P6，保护作用就消失了。

子图 a：随机森林分析找到“关键先生”，在纷杂的数据中，用机器学习算法客观地找出影响植物健康的最重要变量，输入20个变量：包括14个单一菌株的丰度、病原菌组总丰度、共生菌组总丰度、细菌总丰度等。输出目标：预测植物的鲜重。算法：随机森林。它会评估每个变量对预测准确性的贡献度，贡献越大，重要性得分（Mean Decrease Accuracy）越高。最长的条形属于 P6，其重要性得分远高于其他任何变量，包括细菌总丰度。结论：在所有细菌中，P6 这个具体菌株的丰度，是决定植物生死存亡的最关键因素。这完美呼应了图4中发现的植物选择性抑制ATUE5谱系的现象，并将范围缩小到了一个具体的靶标——P6。子图 b：关键菌株 P6 的丰度具有基因型依赖性，目的：验证 P6 的丰度是否在不同植物基因型间存在差异，从而解释保护效果的基因型依赖性。横轴：六种拟南芥基因型。纵轴：细菌丰度（对数刻度）。展示了在每种基因型中，所有14个菌株的丰度估计（点与误差线）。灰色阴影区域 特别突出了 P6 的丰度范围。在受到保护的基因型（如 Lu3-30, kus3-1）中，P6 的丰度（灰色区域）非常低，处于所有菌株的底端。在未受保护的基因型（如 He-1, Ey15-2）中，P6 的丰度（灰色区域）非常高，甚至可能是最高的菌株之一。子图 c：移除验证实验（“分子手术刀”），目的：进行最严格的因果验证——如果 P6 是关键，那么把它从混合菌群中拿掉，保护作用是否还会存在？结论：P6 是共生菌群产生保护作用的必要条件。共生菌的保护机制必须通过抑制或抵消P6来实现。没有P6这个“靶子”，共生菌的“警报”和植物的“防御程序”即使被激活，也无法转化为实际的保护表型。这就像一场防御演习，如果假想敌（P6）不存在，演习就无法检验和产生真正的防护力。

The community structure in the plant-associated microbiome depends collectively on host–microbe, microbe–microbe and host–microbe–microbe interactions. 植物相关微生物群落的群落结构共同依赖于宿主-微生物、微生物-微生物和宿主-微生物-微生物的相互作用。