两株芽孢杆菌PGPB菌株在小麦与大豆中的作用：小麦生长促进未伴随可检测的根际微生物组重构

《International Journal of Molecular Sciences》：Two Bacillus PGPB Strains in Wheat and Soybean: Wheat Growth Promotion Without Detectable Rhizosphere Microbiome Restructuring

　　

植物促生细菌（PGPB）作为生物肥料的应用日益广泛，然而接种菌的基因组潜力与其在植物上实现效应之间的联系，鲜少在联合捕获根际微生物组、植物表型与菌株基因组的整合框架内被评估。两株芽孢杆菌（Bacillus）菌株——盐耐受芽孢杆菌（B. halotolerans）1453与短小芽孢杆菌（B. pumilus）630——在因子盆栽实验中施用于小麦与大豆（2菌株×2施用方法×3频次+对照，3–4重复）。根际样本（过滤后n=67）通过16S rRNA测序结合PICRUSt2功能预测及组成学验证（Aitchison PERMANOVA、ALDEx2、ANCOM-BC2）进行剖析。PGPB基因库通过基因组挖掘（481个标记基因，14类）表征。小麦表型（六个性状）与大豆株高采用适用于计数数据的模型（负二项与二项广义线性模型）进行处理-对照比较，并用因子方差分析分解主效应与交互作用。作物种类是塑造微生物组结构与功能的主导因素（PERMANOVA分类R2=14.7%，功能R2=7.8%，均p<0.001），小麦与大豆间存在具生物学意义的分类差异；菌株、施用频次与方法对群落组成无显著影响（各R2<4%），共现网络在控制读长深度与样本量后作物间无可靠差异。尽管微生物组层面呈中性，接种显著提高了小麦穗数（NB-GLM，全部12处理vs对照，padj 0.0002–0.031）、穗重与茎数，其中施用频次为变异性最强来源，并存在显著的菌株×施用频次交互。菌株1453在穗数（+23.1%，p=0.012）与穗重（+20.4%，p=0.023）上优于630；研究人员假设这可能与其更完整的DNRA途径（narGHI+nirB-nirD）及生物防治基因（bacE、srfAA）有关。菌株630效应较弱但重复间波动更小（CV≈16–21% vs 1453的≈24–26%），可能反映其对环境波动的更好韧性，或因其已确认的rsbV/rsbW胁迫耐受调节子。根际微生物组组成依作物（小麦vs大豆）清晰分化，但对菌株、施用方法或频次无可见响应。尽管缺乏微生物组信号，接种仍显著增加了小麦穗数与穗重，效应大小与稳定性因菌株而异。研究人员假设这种菌株依赖性差异与潜在基因组差异——尤其氮代谢（DNRA途径）与胁迫耐受基因——有关，尽管此联系未经直接检验，仍为未来工作的假说。

研究背景方面，植物促生细菌（PGPB, plant growth-promoting bacteria）通过固氮、解磷、植物激素产生、病原生物防治及诱导胁迫耐受等机制提升作物生产力。芽孢杆菌（Bacillus）属因其产芽孢能力、抗菌脂肽产生及根际定殖能力成为商品化生物肥料中最常用PGPB类群之一。尽管研究数十年且生物肥料市场增长，田间效力仍高度不稳定：同株菌在不同条件下表现迥异，成功案例的机制常未在单次实验内直接确立。根际微生物组主要由寄主植物身份塑造，作物特异性根系分泌物筛选不同细菌群落，功能冗余（functional redundancy，不同分类单元执行等价代谢功能的能力）可能缓冲接种菌扰动。这引出核心问题：若根际微生物组具功能冗余且由寄主决定，PGPB接种菌如何发挥促生效应，何种竞争机制负责？提出的非互斥机制包括修改常驻微生物群落、接种菌与植物直接生理生化互作（氮磷供给、激素信号、铁载体（siderophore）递铁、ACC脱氨酶调节乙烯）、以及诱导系统抗性（ISR, induced systemic resistance）。厘清机制需联合捕获微生物组状态、植物表型与菌株基因组潜力的整合路径。研究人员选择芽孢杆菌为焦点属，选用盐耐受芽孢杆菌（B. halotolerans）1453与短小芽孢杆菌（B. pumilus）630，在小麦（Triticum aestivum）与大豆（Glycine max） factorial 盆栽实验中测试接种是否产生可检测根际微生物组结构/功能偏移，以及表型效应是否依赖此种偏移。

主要关键技术方法方面，研究采用温室因子盆栽设计（2菌株×2施用方法×3频次+对照，每处理3–4重复，样本来自小麦与大豆根际，n=67经滤后进入统计），根际DNA经MagBeads提取、16S rRNA V3–V4区扩增与Illumina MiSeq双端测序，DADA2去噪获ASV（amplicon sequence variant）；分类学以SILVA v138.2注释，β多样性用Bray–Curtis/Jaccard结合PERMANOVA（adonis2）、深度协变量校正与db-RDA，差异丰度用ANCOM-BC2、DESeq2、LEfSe三角验证；功能预测用PICRUSt2（NSTI评估）获MetaCyc/KO/EC，并以ALDEx2与ANCOM-BC2在CLR（中心对数比）变换下做组成学校验；菌株检测用100%一致ASV匹配与二项占用模型；基因组潜力用DIAMOND+HMMER3+Prokka对481个PGPB标记基因（14类）挖掘定级；表型用负二项GLM（NB-GLM）与二项GLM处理穗数等计数性状，因子ANOVA（2×2×3）分解主效与交互，CV%评估稳定性；体外生化面板测脯氨酸、渗透胁迫生长、解磷、铁载体、赤霉素、IAA。

研究结果部分，2.1节根际宏基因组分析显示，分类学上α多样性无处理间差异，β多样性作物间极显著（Bray–Curtis R²=14.7%，校正深度后9.5–10.7%仍显著），小麦富集链霉菌（Streptomyces）等，大豆富集Reyranella、Novosphingobium等豆科特征类群；功能层Aitchison PERMANOVA R²=7.8%，作物间三条通路（同型乳酸发酵、肽聚糖生物合成II、L-赖氨酸生物合成II）经三法共证富集于小麦，大豆两条（L-Phe/L-Tyr超通路、藻尿胆素合成）获两法确认，菌株/方法/频次无显著功能关联；接种菌检测中630在1/30处理以100% identity出现28条读长，1453未检出，二项模型表明中位小麦深度402下0.1%丰度检出率仅33.1%，不能确认或排除定殖。2.2节表型效应表明，播种10天苗高无差异，收获期穗数经NB-GLM全部12处理显著高于对照（p_adj 0.0002–0.031，CK穗数3.0 vs处理6.25–17.0），穗重10/12处理显著，有效分蘖比（穗/茎）9/12处理显著；因子ANOVA示施用频次为5/6性状主效（η² 0.27–0.44），菌株仅对穗数（η²=0.20）与穗重（η²=0.17）主效且1453占优，方法仅穗数与总重显著；菌株×频次交互对茎数（η²=0.27）与穗数（η²=0.24）显著，630在4×频次穗数由13.8降至6.8（p<0.0001）而1453由15.3至12.5不显著（p=0.158），2×为两株实用最优；1453均值更高但CV更大（穗重25.5% vs 630的16.4%），呈生产力–稳定性权衡；籽粒元素（P₂O₅、NH₄等）各处理低于对照，解释为生长稀释。2.3节基因组解读中，体外生化呈互补特化：630产脯氨酸70.0μg/mL、15% PEG下生长36.9%，1453解磷95.3μg/mL、铁载体35.6%、赤霉素1.56mg/mL；基因组挖掘两株共享270/481标记基因，1453独有bacE（杆菌溶素）、srfAA（表面活性素）、nprE，630独有rsbV/rsbW（σ^B调控子）；DNRA（dissimilatory nitrate reduction to ammonium，NO₃^?→NO₂^?→NH₄⁺）途径中1453具narG/H/J与nirD而630缺nirD，致630在亚硝酸盐步阻断，研究人员据此赋氮代谢假说高置信度，并假说4×频次下630亚硝酸盐累积致毒解释交互；生物防治与胁迫耐受基因差异对应1453高生产力高变异、630低变异稳植被。

讨论部分总结，作物身份以根分泌物筛选主导根际组装，本实验同土隔离土壤混杂使R²=14.7%为较纯寄主特异估计；分类作物效应（14.7%）大于功能（7.8%）印证功能冗余——不同分类单元承载重叠功能核心，与Castellano-Hinojosa & Strauss跨六作物“分类核心小、功能核心大”一致，也与Matthews等功能分组预测力优于分类分组吻合；PICRUSt2局限（NSTI中位0.0018但预测与实测KO相关仅~0.20）使功能层差异小既可为冗余亦可为预测天花板。微生物组无响应不排除根组织定殖，仅约束土壤池信号；接种效应不赖群落重构，支持直接植物–微生物互作或ISR而非“修改常驻群落”机制。DNRA完整性（nirD有无）解释菌株×频次交互为可行未证假说，需测根际NO₂^?/NH₄⁺与nirD回补；生产力–稳定性权衡指向1453适可控环境求峰值、630适田间变异求稳。局限含重复少、单时点、麦豆发育阶段混淆、缺灭活接种与载剂对照、籽粒元素为混样描述、大豆表型伪重复、两株自20候选预筛致胜者诅咒偏倚。

结论部分翻译：PGPB接种菌的基因组潜力可预测其植物表型效应，即便可用方法既不能确认也不能排除接种菌定殖。此种解离挑战了PGPB主要通过微生物群落重构起作用的观念，转而支持直接植物–微生物互作模型——这与小麦和大豆根际间稳健具生物学意义的分类差异本身独立于接种形成对照。DNRA途径是菌株特异性剂量反应的合理但未经验证决定因子：1453含nirD使硝酸还原至铵完整，而630缺nirD可能在高施用频次致亚硝酸盐累积与植物毒性。由生物防治（1453）介导的生产力与由胁迫耐受（630）介导的稳定性之间的权衡表明，PGPB产品设计应同时考虑促生效应幅度及其对条件变异的稳健性。该研究发表于《International Journal of Molecular Sciences》。