长期氮添加改变土壤微生物群落结构及潜在的碳-氮与多环芳烃代谢功能

《Emerging Contaminants》：Long-Term Nitrogen Addition Shifts Soil Microbial Community Structure and Potential Carbon-Nitrogen and PAH Metabolic Functions

　　

在一项为期15年的田间试验中，研究人员设置了4个氮（N）施肥水平，采集表层土壤（0–20 cm）以测定其对土壤理化性质、微生物群落结构、探索性共现网络、宏基因组功能谱以及残留多环芳烃（PAH）的影响。相对于对照（CK），高氮处理（N3）中总氮和钙浓度分别为对照的2.25倍和1.38倍。氮添加与探索性微生物网络拓扑差异相关：在N3中，细菌平均度增加15%，真菌增加31%；而在低氮处理（N1）中，细菌和真菌模块度分别降低12%和30%。高氮处理还表现出碳代谢相关基因相对占比升高，以及PAH降解途径（ko00624）中基因相对丰度的变化。Mantel检验表明微生物群落组成与残留PAH谱之间存在依赖于处理的关联，尤其在低氮下明显。总体而言，长期氮施肥与土壤微生物群落及其碳、氮和PAH代谢相关功能潜力的改变相关。这些基于关联的结果并未证明途径活性、PAH降解速率或环境风险；要得出此类推断需直接的过程速率与生物有效性测量。上述发现支持兼顾作物养分需求与遗留有机污染物持久性的位点特异性氮管理。

研究背景与立项依据：土壤微环境受大气沉降与农业活动等多重驱动，长期过量施氮可造成氮盈余、酸化及养分循环失衡，进而重塑微生物群落。华北平原某长期污水灌溉农田遗留有多环芳烃（PAH，polycyclic aromatic hydrocarbons），而氮施肥、微生物演替与残留PAH转化三者间的关系尚不清楚。此前同站点研究提示施氮可改变PAH在各密度组分的分布，但未从微生物功能潜力层面解析。为此，研究人员依托河北石家庄始于2004年、采样于2019年的15年尿素梯度试验（CK、N1=100 kg N ha^-1 yr^-1、N2=200、N3=300，每处理3个小区重复），联合扩增子测序、鸟枪法宏基因组与PAH残留化学分析，检验“长期氮添加以剂量依赖方式改变土壤化学与微生物结构，并与残留PAH谱及PAH降解功能潜力相关”的假设。论文发表于《Emerging Contaminants》。

主要技术方法概览：研究以12个小区级表层土（0–20 cm）为样本队列，分别进行土壤pH、总有机碳（TOC）、总氮（TN）等11项元素与阴/阳离子测定；PAH经冻干、索氏提取、硅胶柱净化后由气相色谱-质谱（GC-MS）定量16种美国环保署优先PAH；微生物端并行两条独立流程——细菌16S rRNA V3–V4区扩增子测序（Illumina MiSeq）与全宏基因组鸟枪测序（Illumina NovaSeq 6000），后者经DIAMOND/KofamScan注释KEGG同源组（KO）并做通路水平reporter-score分析；统计上采用Kruskal-Wallis检验、DESeq2、基于Bray-Curtis的PCoA与ANOSIM、Spearman（|r|>0.8）共现网络（Gephi计算度/模块度等）、Mantel检验关联微生物-环境-PAH矩阵，所有处理内网络与Mantel均因每处理仅3个小区样本而作探索性解读。

研究结果

3.1 氮改良土壤理化性质：随氮水平升高，TN在N3达CK的2.25倍，钙1.38倍，锰3.74倍，铁升高，钾降至CK的0.67倍；CK均pH 7.93、N3均pH 6.53但归档单变量分析未达显著。Mantel显示微生物组成与pH、N、Ca、Cu等协变，CK组群落与N/P/Ca/Mg/S/Zn正相关，N3转而与N/Cu相关，说明氮梯度作为环境过滤器重分配类群而非单纯毒效。

3.2 施氮对微生物群落结构的影响：细菌Shannon多样性在N2/N3降低；Proteobacteria（尤其Deltaproteobacteria、Gammaproteobacteria）随N升高占比升，Chloroflexi、Bacteroidota降，Acidobacteriota升。真菌N1与CK相似、N2/N3偏离。原生生物/真核数据中未分类真核在N1升至56.51%、N3降至34.88%。中性模型（NCM）估得细菌迁移率m=0.652、R²=0.8712，真菌m=0.0214、R²=0.5852；细菌生态位宽度随N收窄、真菌拓宽。OTU水平ANOSIM支持处理间差异（R=0.60802，p=0.001），PCoA轴1解释53.61%变异。

3.3 相关网络与KEGG富集分析：细菌共现网络在N3节点微减、边增12%、平均度增15%、密度增17%、模块度降12%；真菌边增33%、平均度增31%，N1模块度降30%。KEGG reporter-score显示N1分泌系统（3.69）富集，N3甲烷代谢（5.00）、碳代谢ko01200（4.33）、多样环境微生物代谢（3.80）、ABC转运器（3.80）富集。上述仅为功能潜力相对表征，非途径激活证据。

3.4 施氮后微生物群落对残留PAH的响应：1343个细菌OTU为四处理共有；Bacteroidota、Actinobacteriota与高环PAH负相关，Nitrospirota、Zixibacteria正相关；真菌仅Basidiomycota与6环PAH（r=0.72）等少数显著。宏基因组中PAH降解途径ko00624注释到pht3、pht5、ligA、pcaG、todD、tcbB等基因相对丰度变化。Mantel热图显示N1以正关联为主、CK/N2/N3多负关联。研究人员强调丰度≠表达≠降解速率≠风险。

讨论总结：全文反复界定“关联非因果”——网络拓扑不等同稳定性，KO丰度不等同途径流量，PAH残留谱与微生物协变不排除土壤有机质、pH、作物吸收等混杂。单站点、单氮肥（尿素）、每处理3重复限制了推广性；真菌/原生生物引物与中间组装参数未完整留存，属可重复性约束。结论段直译如下：在该15年田块，氮施肥与土壤化学、微生物组成、生态位宽度、共现拓扑及宏基因组功能谱变化相关，且呈剂量依赖——细菌生态位收窄、真菌拓宽，高氮下碳与甲烷代谢基因及PAH降解途径注释基因相对占比升高，残留PAH谱与处理依赖的群落/功能潜力协变；但因基于基因丰度、共现与Mantel关联，未证明途径活性、生态稳定性、污染物固持或实际PAH降解速率。位点特异性氮管理应同时考量作物养分需求与遗留污染物持久性；未来多站点研究须联合残留化学、PAH生物有效性、转录/酶活及直接降解 assay 以建立机制与风险导向施肥指南。