《Microorganisms》：Microbes Drive Straw Decomposition and Microbial Metabolism in Mollisols with Different Straw Return Rates
编辑推荐：
秸秆分解过程高度复杂，受多种相互作用因素调控。然而，秸秆还田比例如何影响秸秆分解过程中的代谢副产物、胞外酶（extracellular enzymes）及微生物群落，仍属未解问题。为此，研究人员整合了秸秆与土壤化学、扩增子测序、非靶向代谢组学及酶活测定，在30
秸秆分解过程高度复杂，受多种相互作用因素调控。然而，秸秆还田比例如何影响秸秆分解过程中的代谢副产物、胞外酶（extracellular enzymes）及微生物群落，仍属未解问题。为此，研究人员整合了秸秆与土壤化学、扩增子测序、非靶向代谢组学及酶活测定，在30天和90天两个时间点比较了四种秸秆还田比例（不还田、1/3、1/2及全量还田）。90天时1/2还田处理秸秆质量损失最大，达53.7%，而土壤有机碳（SOC）与全氮（N）在1/3还田下最高，较对照分别高出17.6%与17.2%，表明分解与短期土壤碳氮积累呈脱钩关系。微生物群落发生明显时序演替，水解酶从30天至90天上升且与质量损失正相关，氧化酶则与分解无正向关联。代谢组由早期易化物转向后期芳香族与酚类化合物，秸秆与土壤代谢组呈阶段特异性耦合（Procrustes M2=0.15，p=0.001；Mantel r=0.69，p=0.001）。Mantel检验进一步显示代谢组与还田比例及SOC显著相关。偏最小二乘路径模型（PLS-PM）揭示30天至90天间连接还田比例、代谢组、酶与分解的直接及间接路径被重新组织。总体而言，秸秆还田比例通过阶段依赖性生化与微生物变化影响分解，秸秆质量损失与短期土壤碳氮积累脱钩，呈非线性响应且在不同还田比例下达峰。
研究背景与立项依据
东北黑土区秸秆还田是维持土壤有机碳（SOC）与养分循环的基础措施，全球整合分析证实其可显著重建农田SOC库。秸秆分解是由底物化学、微生物群落演替与胞外酶生产共同调控的微生物驱动过程，通常分为早期细菌与水溶性水解酶矿化糖类和氨基酸、后期真菌与氧化酶降解芳香族和酚类残留的两个阶段。已有研究多基于单一固定还田量，关于还田比例如何改变微生物与酶演替时序及其后果仍不清楚。非靶向代谢组学可直视微生物与酶响应的底物环境，但多组学耦合尤其在黑土中还田比例依赖的研究鲜见。关键科学问题在于：提高还田比例是否以牺牲后期SOC积累为代价加速早期分解，这种权衡是线性还是存在剂量阈值。研究人员据此在黑龙江嫩江长期田间试验（49°33′35″N，125°27′05″E，典型黑土）开展90天分解实验，设不还田（CK）、1/3、1/2、全量还田四处理，整合化学、扩增子测序、非靶向代谢组学与酶活测定，旨在阐明还田比例如何通过微生物-代谢物-酶耦合调控分解与固碳平衡。
主要技术方法概括
研究人员依托黑龙江农科院自2012年建立的玉米连作长期田间试验（2018—2020年设四种还田梯度），采用原位网袋法（100目内袋装15g秸秆+20目外袋填200g对应 plot 土壤）埋于0—20cm土层，在30天与90天取样。样本队列来自同一站点四处理×四重复 plot 的秸秆残留与根际土壤。核心技术包括：高通量扩增子测序（细菌16S V3-V4、真菌ITS1，DADA2/ASV）；非靶向UHPLC-MS/MS代谢组（MS-DIAL注释，PCA/OPLS-DA筛选差异代谢物）；五种胞外酶比色/吸光法测定（β-葡萄糖苷酶BG、N-乙酰-β-D-葡萄糖胺酶NAG、亮氨酸氨基肽酶LAP、多酚氧化酶PPO、过氧化物酶POD）；统计上采用PERMANOVA、Procrustes、Mantel、Spearman网络及偏最小二乘路径模型（PLS-PM）解析多组学耦合与驱动路径。
研究结果
3.1 土壤碳氮含量：30天各处理SOC差异有限；90天时1/3还田SOC与增量最高，1/2还田SOC增量为负，TN亦以1/3最高。表明短期土壤C、N积累峰值不在质量损失最大处理，二者脱钩。
3.2 秸秆分解率与碳氮：30天质量损失21%—25%且无处理差异；90天达显著分化，1/2还田峰值53.7%，全量次之，CK与1/3最低。秸秆总C与总N随时间下降，30天C/N随还田比例升高而降低（CK约26降至全量约17.5）。分解对还田比例呈非线性。
3.3 微生物群落演替：PCoA与PERMANOVA显示细菌与真菌群落随时间与处理显著分离（30天细菌dispersion异质需谨慎）。细菌中Alphaproteobacteria在90天升高，Actinobacteriota下降；真菌以Sordariomycetes为主但Eurotiomycetes等比例变动。Full vs CK差异菌属阶段特异：30天细菌Massilia、真菌Talaromyces/Cladosporium富集；90天细菌Pantoea/Tardiphaga/Lysinimonas、真菌Cercospora富集。
3.4 胞外酶活性：BG、NAG、LAP从30至90天显著上升，PPO无变，POD略降。BG随还田量总体走高且在全量最高。水解酶尤其BG与质量损失正相关，PPO/POD弱或负相关。
3.5 秸秆代谢组重编程：PCA/OPLS-DA清晰分30天与90天。30天上调为糖类、氨基酸、有机酸、脂肪酸等易化物；90天转为木质素相关酚、酚酸、黄酮、香豆素等芳香酚类，两阶段共有代谢物仅3个。KEGG显示通路阶段转换。
3.6 多组学耦合：秸秆-土壤代谢组Procrustes M2=0.15（p=0.001）、Mantel r=0.69（p=0.001）紧密耦合。30天与90天微生物属-代谢功能组Spearman网络阶段特异，90天环境因子-真菌群落-代谢组Mantel链接显著。
3.7 分解微生物驱动（PLS-PM）：30天还田比例→代谢组→微生物路径弱，分解主要靠酶介导，氧化酶负向；90天还田比例直接正向作用于分解，并重塑代谢组与水解酶，水解酶正向、氧化酶仍负向，路径结构重组。
讨论部分总结
讨论指出分解与短期SOC积累脱钩源于质量损失仅反映残留消失，而SOC是输入、转化、矿化净结果，本研究未测CO2与SOC组分故机制不能直接解析。代谢组阶段转换反映残留化学改变但不等同矿化路径证明。微生物无简单细菌-真菌更替而是分类组成时空重塑。水解酶而非氧化酶主导质量损失，但酶活不随还田量线性缩放。PLS-PM显示早期酶主导、后期还田比例-代谢组-水解酶主导。局限为仅90天两时点、缺CO2与13C示踪。
结论翻译
秸秆还田比例通过微生物群落、代谢物与酶活性的阶段依赖性耦合塑造秸秆分解。质量损失对还田比例呈非线性响应，从30天21%—25%升至90天1/2还田最大53.7%。残留质量损失与短期土壤C、N积累脱钩：90天1/3还田SOC与TN最高，较对照高17.6%与17.2%，而非质量损失最大处理。这表明驱动分解的因素不同于调控土壤C、N留存的因素。秸秆代谢组由30天易化物转向90天芳香酚类，秸秆-土壤代谢组紧密耦合（Procrustes M2=0.15；Mantel r=0.69），伴随时序微生物组成与酶活变化。水解酶较氧化酶更关联质量损失，还田比例-代谢组-酶-分解关系在阶段间重组。综上，微生物群落通过阶段特异性代谢与酶响应驱动秸秆分解，残留崩解与土壤C积累为非线性脱钩过程，须由其底层机制理解。该文发表于《Microorganisms》。
