反硝化菌群落特征决定缺氧条件下农业土壤反硝化与N2O消耗潜力的上限

一氧化二氮（N₂O）是一种主要的长寿命温室气体，其大气浓度持续上升。当氧气受限时，反硝化微生物（denitrifiers）在土壤中既是N₂O的重要来源，也是重要汇。反硝化潜力高度依赖多种环境因子，而微生物群落特征可能是控制反硝化速率的关键决定因素。然而，跨不同系统适用的微生物群落结构与过程潜力之间的稳健联系仍然缺乏。研究人员提出假设：对比鲜明的土壤因子（edaphic factors）塑造了彼此不同的土壤反硝化菌群落，而这些群落具有不同的动力学特征，从而限制土壤在厌氧条件下的反硝化能力和N₂O消耗能力。因此，研究测定了3种差异显著土壤在缺氧微宇宙（anoxic microcosms）中的硝酸盐与N₂O消耗表观米氏动力学（Michaelis-Menten kinetics）。硝酸盐和N₂O的表观最大反应速率分别介于40–340与20–540 nM min^-1。氮素质量平衡表明，在NO₃^-还原过程中，反硝化而非异化性硝酸盐还原为铵（DNRA）作出了显著贡献。对16S rRNA及编码反硝化和DNRA相关酶的标志基因进行定量PCR（qPCR）和扩增子测序表明，3种土壤具有差异显著的反硝化菌与DNRA群落。Bradyrhizobium sp.是最具区分性的反硝化类群，而Opitutaceae和Gemmatimonadaceae则是决定N₂O还原群落差异的关键类群。估算的“细胞”特异性最大反应速率在土壤反硝化菌中为1–10 fmol gene^-1 h^-1，在N₂O还原菌中为5–90 fmol gene^-1 h^-1。编码细胞色素依赖型亚硝酸盐还原酶的基因nirS以及编码N₂O还原酶的基因nosZII的丰度与表观最大反应速率相关性良好。上述数据表明，反硝化活性和N₂O排放是在特定反硝化菌群落能力所界定的边界条件内被调节的，这一认识（i）对氮循环（N-cycle）建模具有意义，（ii）有助于揭示微生物功能群及其活性的重要生态驱动因素，并为开发与微生物组活性水平相匹配的管理策略提供基础。

该论文发表于《Soil Biology and Biochemistry》，聚焦农业土壤在缺氧条件下反硝化与N₂O还原潜力的群落生态学控制机制。研究背景在于，N₂O是一种具有极高增温效应的长寿命温室气体，农业施氮导致其排放持续增加，而土壤中的反硝化过程既可能产生N₂O，也可能进一步将其还原为N₂。因此，决定土壤究竟表现为N₂O“源”还是“汇”的关键，不仅在于外部环境条件，也在于土壤中承担相关代谢功能的微生物群落。然而，既往研究虽然已认识到pH、有机碳、C/N比、氧限制、土地利用和水分等土壤因子对反硝化具有重要影响，但不同土壤系统中“群落结构—动力学参数—过程潜力”之间仍缺乏稳健、可量化的联系。尤其是，真实土壤群落的表观米氏动力学参数、功能基因组成与关键类群如何共同限定反硝化及N₂O消耗能力上限，尚未得到系统阐明。基于此，研究人员提出，差异化的土壤理化因子将塑造不同的反硝化菌群落及其动力学性状，进而约束农业土壤在标准化缺氧条件下的反硝化和N₂O还原潜力。

为验证上述假设，研究人员选取德国3种性质差异显著的农业土壤与草地土壤作为模型系统，分别为Fuhrberg的gleyic podzol、Rotthalmünster的haplic luvisol以及Gießen的fluvic gleysol。研究在标准化缺氧微宇宙中分别加入NO₃^-或N₂O作为末端电子受体，并设置是否补充乙酸盐的处理，以解析电子供体可利用性对过程速率的影响。结果表明，乙酸补充对反硝化和N₂O消耗的动力学参数几乎没有影响，说明在本研究设定下三种土壤的过程潜力主要受NO₃^-或N₂O供应及群落能力限制，而非受碳源限制。与此同时，氮素质量平衡和NH₄⁺变化结果共同显示，NO₃^-还原主要流向反硝化，而不是异化性硝酸盐还原为铵（DNRA）。

研究采用的关键技术方法主要包括以下几类：首先，基于3个德国不同样点来源土壤构建缺氧土壤浆液微宇宙，通过机器人联用气相色谱系统实时监测O₂、CO₂、NO和N₂O，并依据初始速率拟合表观米氏方程，计算v_max,app与K_m,app。其次，对培养前土壤进行DNA提取，利用定量PCR（qPCR）定量16S rRNA、narG、napA、nirK、nirS、nosZI、nosZII和nrfA等功能标志基因。再次，采用Illumina MiSeq扩增子测序解析16S rRNA与多种氮循环功能基因群落结构，并结合QIIME2、DADA2、Aitchison距离主坐标分析（PCoA）开展生物信息学分析。最后，通过绝对丰度换算、相关性分析和群落排序，识别与过程参数密切相关的关键类群与指示基因。

在研究结果部分，论文首先以“Effect of supplemental organic carbon on anaerobic N-oxide dissimilation”说明补充有机碳对厌氧氮氧化物异化的影响。该部分通过NO₃^-和N₂O添加培养实验发现，3种土壤均可迅速启动反硝化和N₂O还原，但过程强度显著不同：Gießen土壤最高，Fuhrberg土壤最低。乙酸既未被耗尽，也未显著提升N₂O生成或消耗速率，表明体系并非碳限制。NO始终低于检测限，而NH₄⁺多表现为净消耗，进一步支持DNRA贡献有限、反硝化占主导的判断。

在“Apparent Michaelis-Menten-Kinetics for N₂O production and consumption”部分，研究人员通过初始速率拟合量化了三种土壤的表观米氏动力学特征。NO₃^-浓度超过100 μM后，反硝化N₂O生成速率趋于饱和；N₂O浓度达到10–20 μM以上时，N₂O还原速率达到最大。Gießen土壤无论在反硝化还是N₂O还原中均表现出最高的v_max,app和较高K_m,app，Fuhrberg最低，Rotthalmünster居中。这说明不同土壤微生物群落在底物亲和性与最大反应能力方面存在显著差异，且这种差异在统一培养条件下依然保留，提示其具有稳定的群落性状基础。

在“Genetic denitrification potentials”部分，研究通过qPCR比较了各土壤总菌量及氮循环功能基因丰度。Gießen土壤的narG、napA、nirK、nirS、nosZI、nosZII和nrfA绝对丰度整体最高，Fuhrberg次之，Rotthalmünster最低；但以16S rRNA标准化后，土壤间差异缩小，仅部分基因如napA、nirK和nrfA仍保持显著差别。这表明决定过程潜力的不仅是总微生物量，还包括功能群在群落中的相对构成和功能配置。

在“Correlation of apparent maximum velocities and gene abundances”部分，研究将v_max,app与基因丰度进行Pearson相关分析，发现nirS与nosZII丰度与N₂O生成和消耗最大速率的相关性最为突出，nrfA也表现出一定正相关趋势。虽然统计显著性受样本点数量限制，但这一结果提示，nirS可能是表征土壤最大反硝化活性的潜在分子指标，而nosZII可能是表征土壤N₂O汇能力的重要代理指标。与此同时，nir(K+S)/nrfA和nir(K+S)/nosZ(I+II)比值均高于1，说明这三类土壤在遗传潜力上均更偏向反硝化而非DNRA，并且N₂O产生潜力整体高于完全还原潜力。

在“Phylogenetic and functional genetic microbial community traits”部分，16S rRNA扩增子结果显示，Acidobacteriota、Actinomycetota、Chloroflexota、Bacillota、Gemmatimonadota、Pseudomonadota和Verrucomicrobiota为优势菌门，不同土壤在家族水平上具有显著分化。PCoA结果表明，三种土壤的总体细菌群落以及各功能基因群落均明显分离，提示土壤类型对应不同的群落组装模式。更重要的是，多标志基因分析确定Bradyrhizobium sp.为最具区分性的反硝化类群，而Opitutaceae和Gemmatimonadaceae则是N₂O还原群落分异的关键决定者。各土壤中优势功能类群亦不同，例如Fuhrberg中Bradyrhizobium sp.、Mesorhizobium sp.、Azoarcus sp.和Deinococcus sp.较为突出，Rotthalmünster中则有Mycolicibacterium sp.、Marinobacter sp.、Skermanella sp.和Gemmatirosa sp.，Gießen中则以Piscinibacter sp.、Variovorax sp.、Bradyrhizobium sp.及Opitutaceae为代表。该结果说明，不同土壤并非只是拥有不同数量的反硝化微生物，而是拥有不同组成、不同功能模块搭配的反硝化与N₂O还原联合体。

在“Estimation of ‘cell’-specific reaction rates”部分，研究基于最丰富的nirK和nosZII基因丰度，估算了“细胞”特异性反应速率。结果显示，反硝化相关“细胞”特异性最大速率约为1–10 fmol gene^-1 h^-1，而N₂O还原相关速率约为5–90 fmol gene^-1 h^-1，后者约高一个数量级。尽管这一估算受DNA提取和引物偏倚影响，仅代表数量级水平，但仍清晰揭示不同土壤中功能模块能力存在显著差别，而且过程上限并不完全由功能基因数量单独决定，群落结构本身同样重要。

讨论部分强调，本研究将群落水平米氏动力学与多标志基因群落结构分析结合起来，填补了“过程参数—群落组成—关键类群”联动研究的不足。研究表明，pH、C/N比、土壤质地等土壤因子长期塑造了差异化的反硝化菌与N₂O还原菌群落，并使其在标准化缺氧条件下仍表现出不同的底物亲和性和最大反应速率。研究还指出，一些传统上更多被视为固氮类群的自由生活α-变形菌，如Bradyrhizobium sp.，可能是农业土壤中被低估的重要反硝化者；而携带clade II nosZ的Gemmatimonas sp.与Opitutaceae，则可能在降低土壤N₂O排放中发挥关键作用。与此同时，DNRA总体潜力较弱，但部分DNRA相关类群可能仍参与N₂O消耗，提示土壤N氧化物呼吸网络具有更高的功能耦联性。研究同时谨慎指出，相关性分析仅基于3种土壤，因而更适合视为趋势性证据，仍需更多土壤样本加以验证。

研究结论可译为：研究结果强调，局地土壤因子如pH和C/N比通过塑造反硝化菌群落组成，进而约束了土壤即使在电子供体过量的“理想化”缺氧条件下的反硝化和N₂O消耗过程潜力。差异显著的表观米氏动力学参数与反硝化菌群落结构相关。编码细胞色素c依赖型NO₂^-还原酶的nirS基因和编码clade II N₂O还原酶的nosZII基因丰度，可能解释表观最大反应速率，提示其在反硝化和N₂O还原建模中的应用价值值得进一步研究。多标志基因分析提示Bradyrhizobium sp.是迄今可能被低估的土壤反硝化类群，而Opitutaceae和Gemmatimonas sp.则是具有重要N₂O还原能力的土壤微生物。总体而言，研究结果强调，不同微生物作用者之间的相互作用，包括仅编码部分反硝化通路的微生物，共同导致了不同的反硝化过程潜力。尽管本研究测得的是标准化条件下的过程潜力上限，原位土壤条件仍可能因土壤结构和异质性而产生不同结果，但本研究为未来将微生物群落结构及其动力学参数整合进反硝化模型提供了定量基础。