大蒜（Allium sativum）种植管理驱动的土壤微生物组变化及其与土壤健康和农业生产力的关联：基于全长16S rRNA基因纳米孔测序与生产力加权评分（PWS）的多农场田间研究

《Journal of Soil Science and Plant Nutrition》：Changes in Soil Microbiome Linked to Management of Allium sativum Cultivation: Impacts on Soil Health and Agricultural Productivity

　　

大蒜（Allium sativum）产生具有抗菌活性的有机硫化合物，但在田间条件下其种植如何改变土壤细菌群落、以及这些变化与产量如何关联，仍认识不足。研究人员在巴西南部四家土壤与管理条件各异的商业化大蒜农场，于播种前与收获后采样（每农场每时间点三个样本）。采用Oxford Nanopore技术对全长16S rRNA基因进行测序。通过排序分析、相似性分析（ANOSIM）与多样性指数比较群落，并以Welch t检验比较类群。以生产力加权评分（PWS）按产量加权的相对丰度变化对物种排序，并用PICRUSt2预测功能谱。种植后群落组成发生变化（ANOSIM R = 0.34，p < 0.001），而α多样性未变（p > 0.05）。Ehrlichia由15.2%降至5.6%，Legionella由7.3%降至1.4%，而Bacteroidota由23.3%升至39.4%，主要通过Pedobacter与Flavobacterium（均p < 0.05）。Flavobacterium buctense与Chryseolinea serpens具有最高的正向PWS，由高产农场中的增加所驱动，而Legionella spp.最为负向。19条预测的MetaCyc通路在采样时间之间存在差异（校正p < 0.05），主要涉及中心碳代谢、芳香族化合物与多胺代谢。大蒜种植及其管理措施与土壤细菌群落在条件各异的农场间的一致重组相关，且未伴随多样性丧失。PWS提供了一种可重复的方法，用于从田间数据中对候选类群排序，以供在受控的多点研究中评估；预测的功能变化是可供直接测定的假设。

研究背景与问题。土壤微生物组是地球上功能多样性最高的生态系统之一，每克土壤可含多达10⁹个微生物细胞，驱动支撑农业生产力的关键生物地球化学过程。随着全球农业追求更高产量与更低环境代价，利用土壤微生物群落已成为可持续农业的核心议题。植物被视为"全生物体"（holobiont），即宿主与微生物群及其复杂信号机制构成的整合系统，土壤微生物通过养分活化、病原抑制与植物激素产生影响植物营养、胁迫韧性与生长。大蒜（Allium sativum L.）因合成多种具强抗菌活性的有机硫化合物而成为研究作物—微生物组互作的理想案例：组织损伤时，蒜氨酸酶（alliinase）将蒜氨酸（alliin）转化为大蒜素（allicin），后者破坏微生物膜、抑制关键酶并干扰群体感应（quorum sensing）。原位条件下，这些生物活性物质施加的选择压力会重构微生物网络，而非简单清除敏感物种，其效应取决于土壤质地、管理历史与环境背景。尽管大蒜抗菌性的体外证据充分，但田间尺度大蒜种植如何改变土壤微生物组、以及这些变化如何与作物表现相关联，仍知之甚少。本研究即针对这一问题：在土壤与管理条件各异的农场中，大蒜种植及相应管理措施如何改变土壤细菌群落，这些变化又如何与各农场获得的产量相关联。研究假设为：大蒜种植及相应管理驱动土壤细菌群落分类组成与相对丰度的变化，且变化的方向与幅度在产量不同的农场间存在差异，从而可识别其变化与真实种植条件下较高或较低生产力相关的类群。研究目标为：（1）量化四个商业化农场一个蒜季中土壤细菌群落的重构；（2）利用考虑各农场产量的生产力加权排序，结合农场水平与合并检验，识别相对丰度变化与大蒜产量相关的细菌类群；（3）通过组成性微生物组数据的差异丰度分析，结合NanoCLUST分类、PICRUSt2的MetaCyc通路功能推断与LinDA线性模型框架，表征种植前后的功能谱；（4）评估农场特异性因素对微生物组响应的潜在调节作用。

关键技术方法（不超过250字）。研究在巴西南部圣卡塔琳娜州高海拔（881–986 m）四家商业化紫皮大蒜农场（Farm 3–6，土壤类型、管理历史与产量各异）采样，于2021年7月16日播种前与11月19日收获后各取样，每农场设4个采样点、处理其中3个，共24个样本，取0–20 cm耕层复合土样。提取土壤总DNA后，用携带条形码27F/1492R引物的16S条形码试剂盒扩增大肠杆菌全长16S rRNA基因（约1,500 bp），在MinION（Mk1C）测序仪R9.4.1流动池上测序。生物信息学采用Guppy碱基识别、Porechop去条形码与接头，NanoCLUST进行物种水平聚类（最小簇50条读长）并生成一致性序列，经BLASTn比对NCBI 16S核糖体RNA参考数据库完成分类。统计上采用Bray–Curtis非度量多维标度（NMDS）、ANOSIM、Kruskal–Wallis检验、Welch t检验，并以PICRUSt2结合LinDA进行通路差异丰度分析。

研究结果。3.1 土壤理化性质与大蒜产量：四农场理化背景差异显著，pH 6.05–6.81，质地从Farm 4的重黏土到Farm 6的壤砂土；有效磷将站点分为相差一个数量级的两组，Farm 4与6较高；有机碳与全氮以Farm 3最高、Farm 6最低。产量从Farm 3的3,500 kg ha^?1到Farm 6的11,000 kg ha^?1，随磷梯度而非碳氮梯度变化；Farm 3为唯一观察到镰刀菌（Fusarium）枯萎症状的站点。3.2 种植前后土壤微生物群落结构：基于Bray–Curtis相异度的NMDS显示种植前后微生物组分离（stress = 0.19），ANOSIM确认群落差异显著（R = 0.34，p < 0.001）。判别分析显示农场特异性关联：Farm 3与蓝细菌属Chamaesiphon、Farm 4与Pseudomonas、Farm 5与包括Flavobacterium在内的多个属、Farm 6与Ehrlichia相关联。α多样性指标在种植前后无显著变化（Kruskal–Wallis，p > 0.05）；Farm 4优势度增加约130%（5.53%→12.77%），其余农场优势度下降?29%至?52%；Shannon有效数由33.17 ± 14.16变为31.23 ± 16.39，Simpson有效数由18.03 ± 7.39变为17.26 ± 9.93。Bray–Curtis相似性热图显示种植前样本在农场内与农场间均紧密聚类，种植后样本离散度更大，Farm 6种植后样本与种植前保持高相似性。3.3 种植前后细菌相对丰度比较：四个门（Pseudomonadota、Bacteroidota、Cyanobacteriota、Acidobacteriota）合计约占全部农场土壤微生物组的87%，Pseudomonadota占优（平均37.55 ± 11.45%），门水平稳定性与属、种水平的剧烈变化形成对比。Pseudomonadota内五个属占读长>50%；Ehrlichia下降三倍，在Farm 5最显著（18.10 ± 6.74%→1.37 ± 1.09%；p < 0.05）；Legionella同样骤降，四个种（L. tunisiensis、L. dresdenensis、L. hackeliae、L. qingyii）贡献了大部分降幅（均p < 0.05）；Pseudomonas在Farm 4由0.36 ± 0.32%升至23.21 ± 17.33%（约64倍），与该农场优势度升高一致，但未达统计显著。种植后Bacteroidota显著上升（23.33 ± 9.99%→39.37 ± 14.04%；p < 0.01），除Farm 6外各农场一致；属水平上Pedobacter与Flavobacterium显著增加；新定殖者出现，包括P. terrae（种植后平均>2%）与F. beibuense（Farm 5由0%升至8.34 ± 3.68%），P. soli亦显著上升（0.07 ± 0.15%→1.74 ± 3.11%，约26倍）。Cyanobacteriota门水平变化不显著（14.88 ± 14.72%→10.30 ± 8.47%），但Chamaesiphon在Farm 3大幅下降（24.27 ± 19.57%→0.04 ± 0.07%）；Acidobacteriota显著下降（6.75 ± 3.75%→3.81 ± 2.54%；p < 0.05），由Stenotrophobacter下降驱动；Thermodesulfobacteriota由2.51 ± 2.01%降至0.88 ± 0.83%，40个已识别种中30个呈下降趋势；Mycoplasmatota骤降（2.75 ± 3.59%→0.30 ± 0.36%），以Spiroplasma急剧减少为主。3.4 按生产力加权评分排序的细菌物种：19个物种符合纳入标准（在任一空间尺度上相对丰度变化显著，或位列最正/最负30个评分且平均总丰度≥0.5%）。Flavobacterium buctense正向PWS最高（+5.15），丰度由1.18 ± 1.42%升至1.39 ± 1.57%；Chryseolinea serpens（PWS = +2.94）由0.44 ± 0.65%升至0.84 ± 0.90%；Leadbetterella byssophila（PWS = +1.09）由0 ± 0%升至1.36 ± 3.17%。Pedobacter中P. soli（PWS = +2.54）增加26倍，P. terrae（PWS = +1.21）新出现（0 ± 0%→2.10 ± 3.96%）；Occidentia massiliensis（PWS = +2.89）由0.14 ± 0.23%升至1.13 ± 1.39%。九个物种呈负PWS：Legionella dresdenensis（?1.04）下降91%，L. hackeliae（?0.94）下降83%，L. tunisiensis（?0.76）下降74%，Ehrlichia muris（?0.62）下降65%，Flavitalea flava（?0.59）下降62%，Moheibacter sediminis（?0.25）虽由0.01 ± 0.03%升至2.07 ± 3.92%但评分为负。逐农场数值表明多数评分由一两个农场驱动：F. buctense的评分几乎完全来自Farm 4（相对变化+13.96），而在Farm 5显著下降（3.29 ± 0.85%→0.18 ± 0.20%；p = 0.020）；C. serpens的评分主要来自最高产的Farm 6，其增加显著（0.25 ± 0.42%→1.53 ± 0.29%；p = 0.016）；P. terrae仅在Farm 4种植后检出；P. soli的正评分来自最低产Farm 3的大幅比例增加（0.11 ± 0.09%→4.52 ± 5.41%），而该种在Farm 6种植后不再检出；O. massiliensis的评分主要来自Farm 5。最负的评分反映在全部或多数农场中的下降，如三个Legionella种与E. muris；Niastella vici的负评分则源于其在两个最高产农场的消失，超过其在Farm 3与5的增加。Pedobacter soli（灰色）在农场或合并水平均无显著变化却评分较高，尽管Pedobacter属水平变化显著。3.5 大蒜种植后土壤微生物群落的预测功能谱：LinDA差异丰度分析识别出19条种植前后预测丰度存在差异的MetaCyc通路（校正p < 0.05），其中16条种植后预测丰度降低，3条略升高（log₂FC介于0.019与0.067）。样本加权NSTI为0.390 ± 0.127（范围0.224–0.654）。3.5.1 中心碳代谢的预测变化：种植后与好氧碳代谢相关的通路预测丰度降低，包括三羧酸（TCA）循环IV（log₂FC = ?0.139）、乙醛酸旁路通路（log₂FC = ?0.409至?0.498）以及糖酵解、丙酮酸脱氢酶、TCA与乙醛酸旁路的超级通路（log₂FC = ?0.246）；发酵相关通路中，乙酰辅酶A发酵为丁酸II预测丰度略高（log₂FC = 0.055），而丙酮酸发酵为丙酮与异丙醇生物合成预测丰度较低。3.5.2 芳香族化合物与多胺代谢的预测变化：最大预测差异涉及芳香族化合物，L-苯丙氨酸与L-酪氨酸生物合成超级通路种植后预测丰度降低（log₂FC = ?0.619与?0.542），邻苯二酚降解I与II（间位裂解途径；log₂FC = ?0.499与?0.690）以及苯乙酸降解I（好氧）（log₂FC = ?0.734）亦然；多胺生物合成的三条超级通路预测丰度亦降低（log₂FC为?0.443至?0.595），而L-缬氨酸降解I预测丰度略高（log₂FC = 0.067）。3.5.3 辅因子与胁迫相关通路的预测变化：S-腺苷-L-甲硫氨酸循环I预测丰度略高（log₂FC = 0.019），而NAD补救途径II预测丰度较低（log₂FC = ?0.504）。

讨论与结论。研究结果与大蒜种植及其相关管理在差异显著的土壤条件下作为土壤细菌群落生态过滤器的观点一致。NMDS与ANOSIM支持的群落变化表明作物周期伴随显著周转，尽管样本点间存在一定重叠，样本在排序空间不同象限的分散提示存在响应分化的亚群。这些结果与既往研究一致：Allium物种通过释放有机硫化合物（尤其大蒜素及其衍生物）影响土壤微生物组，其作为广谱抗菌剂抑制某些类群而利于另一些类群。值得注意的是α多样性未显著改变，因此观察到的是群落重组而非多样性丧失，这与微生物群落受扰动时抗性与韧性相区别的模式一致。该重组结合了选择性丧失与新生态位定殖：Bacteroidota的增加以及Pedobacter terrae、P. soli与Flavobacterium spp.在播种前未检出土壤中的出现，与竞争释放后的快速生态位定殖一致，也与中度干扰理论相符——适度扰动通过防止竞争排斥维持多样性。另一方面，降幅最大的Ehrlichia与Legionella所对应的谱系，其最近缘的可培养近亲依赖宿主或原生动物相关生活方式及对环境变化敏感的铁获取系统，这与作物周期中维持其生存的条件丧失相一致。由分类组成变化所预测的功能谱随之改变，研究人员将其解释为关于代谢潜力的假设而非实测活性。多数差异通路种植后预测丰度降低，包括TCA循环与乙醛酸旁路，而乙酰辅酶A发酵为丁酸II预测丰度略高；一种可能解释是群落转向具有更强厌氧碳加工潜力的类群，这可能源于作物周期中有机硫化合物与根系分泌物输入相关的土壤氧化还原条件变化；该解释仅部分得到预测支持，因为另两条发酵相关通路（丙酮酸发酵为丙酮与异丙醇生物合成）预测丰度亦较低，且本研究未测定氧化还原电位；但它产生可检验的预测：若成立，同一田块收获后土壤应表现出比播种前更低的氧化还原电位与更高的发酵终产物形成速率。芳香族化合物降解通路预测丰度降低若被证实，将与土壤有机质循环相关：邻苯二酚是木质素衍生化合物及其他植物芳香物生物降解的中心中间体，邻苯二酚与苯乙酸降解通路预测丰度降低提示处理这些难降解底物的潜力下降；由于芳香族化合物占土壤有机质相当比例，尤其来自富含木质素的作物残体，该潜力的持续降低可能影响大蒜种植土壤的分解与养分循环。另两项观察可作同样解读：L-缬氨酸降解预测丰度在总体下降背景下略升，可能反映能利用支链氨基酸作为替代碳源与能源的类群，这是功能冗余所允许的补偿；NAD补救途径II与多胺生物合成预测丰度降低、S-腺苷-L-甲硫氨酸循环I略升，可能反映氧化还原平衡、生长与甲基化活性潜力的差异，这与有机硫暴露下的细胞修复与膜重塑一致。总体而言，预测谱提示群落具有较低的好氧中心碳代谢与芳香族化合物降解潜力，并在胁迫相关通路上存在差异；这些是源自分类组成的假设而非活性测定。关于PWS，其追问上述变化中哪些在较高产农场中权重更大，为相对丰度的直接比较增加了另一层解释。Flavobacterium buctense评分最高，由其在第二高产Farm 4中的增加驱动；该属成员曾被报道产生生长素与细胞分裂素、合成铁载体并释放诱导系统抗性的挥发性化合物，这些特性与观察到的关联一致，尽管本研究未测定这些功能。两个Pedobacter种因不同原因获得正评分：P. terrae仅在Farm 4种植后检出，而P. soli的评分来自最低产Farm 3的大幅比例增加且种植后在Farm 6不再检出，因此它们与产量的关联弱于整个属的增加，这与生态演替为功能特化类群创造机会相一致；该属物种曾被报道通过分泌有机酸与产生磷酸酶溶解磷酸盐并产生生物活性代谢物，鉴于大蒜鳞茎发育期对磷的需求，这些功能值得关注，尽管本研究未测定溶磷作用。这些关联值得结合各农场土壤磷分布考量：两个最高产站点Farm 4与6也是收获后有效磷较高的站点；但关联类群是否贡献于磷有效性，抑或二者共同响应某一共同土壤条件，无法从观察性设计的四个农场中确定。Flavobacterium与Pedobacter在属水平的增加仍是促动在受控实验中检验基于微生物组的管理策略的那类关联。具复杂生物聚合物降解能力的细菌，包括Chryseolinea serpens与Leadbetterella byssophila，也出现在正评分类群中；C. serpens在最高产Farm 6显著增加，其被报道的有机质分解能力可促进作物残体养分释放，鉴于TCA循环与芳香族化合物通路预测丰度较低，这一点值得关注。Occidentia massiliensis为Rickettsiaceae成员，最初自蜱中描述，也获得高正评分，主要来自Farm 5；此类归属应与底层一致性序列的相似度一并解读。Legionella种的负评分反映其在全部或多数农场（包括两个最高产农场）中的下降；该属成员在土壤常见自由生活原生动物内胞内复制并形成生物膜以增强持久性与抗性，在土壤团聚体内形成受保护微环境并消耗氨基酸与低分子量有机化合物；这些特性为负关联提供了合理解读，此类消耗原则上可降低参与植物—微生物信号传导的化合物可利用性，但本研究未测定此类互作，该关联本身并不确立对植物生长的拮抗效应。排序提示大蒜生产力可能与分类多样性关系较小，而与能在有机硫压力下存续的特定类群的存在关系较大；检验这些类群是否贡献于生产力而非仅伴随其出现，是本研究自然的延续。农场间差异提示土壤理化性质与管理调节微生物组响应：Farm 6为壤砂质地且有效磷最高，兼具最高产量（11,000 kg ha^?1）与采样间群落组成最小变化；而Farm 4为重黏土，表现出最大的群落简化（种植后由Pseudomonas主导）并具较高产量（9,500 kg ha^?1）；因此不止一种群落轨迹与高生产力相容，土壤质地与管理可能参与决定其发生方式。最后，两方面限定了结果范围：其一，研究包含四个商业化农场、每个仅一个产量值，故微生物组变化与生产力的关系仅建立在四个独立的两变量组合之上，所识别类群因此是需在更大多站点集合中确认的候选；其二，各农场同时在土壤质地、前茬作物与管理上存在差异，这些因素彼此不独立，PWS按各农场产量对其变化加权而未分离这些贡献。结论为：大蒜种植及其管理措施与条件各异的农场间土壤细菌群落的一致重组相关，且未伴随多样性丧失；PWS提供了一种可重复的方法，用于从田间数据中对候选类群排序以供受控的多点研究评估；预测的功能变化是可供直接测定的假设。