《Applied Soil Ecology》：Beyond taxonomic richness: Network topology and functional gene abundance reveal tillage and nitrogen as physiological and topological filters across heterogeneous corn agroecosystems
编辑推荐：
土壤健康日益被理解为物理、化学和生物属性相互作用所涌现出的整体特性，而非肥力或养分供应的简单同义词。当代土壤健康框架强调，土壤作为复杂的生命系统运行，其性能取决于土壤结构、化学成分和生物群落之间的相互作用（Bünemann等，2018；Fine等，2017；L
土壤健康日益被理解为物理、化学和生物属性相互作用所涌现出的整体特性，而非肥力或养分供应的简单同义词。当代土壤健康框架强调，土壤作为复杂的生命系统运行，其性能取决于土壤结构、化学成分和生物群落之间的相互作用（Bünemann等，2018；Fine等，2017；Lehmann等，2020）。土壤微生物群落通过驱动有机质分解和养分矿化、调节温室气体通量，同时促进团聚体形成、孔隙度和保水性，支撑着关键的生态系统功能（Bardgett和van der Putten，2014；Francioli等，2016）。
微生物的高生理可塑性和短世代周期使其对环境变化和管理尤为敏感，因此微生物多样性、组成和活性的变化往往先于土壤理化指标的可检测变化（Semenov等，2025）。因此，微生物属性越来越多地被用作土壤退化或恢复的早期预警指标，以及评估土壤健康中生命组分的操作性指标。
在集约化种植体系中，管理措施通过不同的化学和物理途径充当强大的生态过滤器，塑造土壤微生物组的组成、结构和功能。施氮是一个主要驱动因素，因为它同时改变土壤pH值、碳氮比和底物供应。在不同生态系统中，实验性氮添加反复被证明会抑制微生物呼吸，尤其是在高施用量下（Ramirez等，2010a），并持续改变草地和耕地系统中的细菌群落组成（Ramirez等，2010b；Leff等，2015）。全球综合分析和长期实验表明，慢性或过量氮输入通常会降低细菌多样性并促进富营养型类群（通常在变形菌门Proteobacteria和厚壁菌门Firmicutes内），同时抑制寡营养型谱系如酸杆菌门Acidobacteria（Francioli等，2016；Hartmann等，2015；Leff等，2015）。近期研究还表明，细菌群落可能沿施肥梯度表现出非线性响应和阈值行为，超过临界氮负荷后发生突然的组成变化（Chen等，2024）。总体而言，这些研究强调，施氮对土壤细菌群落的重塑作用在组成和性状分布上比常用的α多样性指数更为显著。
补充这些化学驱动因素的是，耕作强度施加了一种不同的、主要是物理的选择压力，通过重组土壤基质、改变通气性和水分、破坏团聚体和真菌菌丝网络来实现。常规翻耕通常会使耕层均质化，使受保护的有机质暴露于氧化，并增加资源可用性的时间变异性。实证研究和荟萃分析表明，相对于少免耕或免耕系统，集约化耕作通常会减少微生物生物量和酶活性，尽管效应大小取决于情境，并受到残茬管理、土壤类型和气候的影响（Helgason等，2010；Lupwayi等，2012；Zuber和Villamil，2016）。保护性耕作如少耕或直播倾向于保留地表残茬，维持更强的有机质垂直梯度，并创造更多样化的栖息地，这通常与更高的细菌多样性和更明显的随深度分层的群落结构有关（Schmidt等，2018；Sui等，2022）。
微生物多样性不仅是群落结构的描述符，也是土壤多功能性的决定因素。全球综合分析表明，较高的微生物丰富度和均匀度通常能增强土壤维持多种功能的能力，包括固碳、养分循环和初级生产力（Delgado-Baquerizo等，2016；Wagg等，2014）。实验工作进一步表明，多样化的微生物群落往往对干扰更具抵抗力和恢复力，在胁迫下维持功能并在之后更快恢复（Griffiths和Philippot，2013）。这些关系将微生物多样性定位为连接土壤健康与生态系统稳定性的中心轴。氮（N）循环说明了微生物群落组成和功能如何共同控制土壤过程。关键的转化过程，如固氮、硝化和反硝化，由携带诸如nifH、amoA和nirK等基因的特化功能类群介导（Kuypers等，2018；Ward和Jensen，2014）。这些类群丰度和相对重要性的变化可以改变氮保留与损失的比例，从而影响植物可利用氮、硝酸盐淋溶和一氧化二氮排放。因此，改变氮循环群落结构的农业管理实践通过改变作物氮营养与环境氮损失之间的平衡，直接影响土壤健康。除了分类学组成和功能基因库之外，微生物相互作用的组织为土壤生态系统稳定性增加了另一个维度。基于类群共变模式的共现网络被广泛用作潜在生态相互作用（如竞争、互惠和互养）的代理（Faust和Raes，2012）。网络属性，包括连通性、模块性、聚类性和集中性，与功能冗余相关联，而功能冗余是微生物群落抵抗和恢复干扰的基础（Biggs等，2020；Louca等，2018；Puente-Sánchez等，2024）。沿着农业集约化梯度，实证研究表明，强化管理可以降低网络复杂性和关键类群的相对丰度，表明在更受干扰或资源更丰富的环境中微生物共现网络趋于简化（Banerjee等，2019；Karimi等，2019）。这些发现意味着网络拓扑可能提供关于不同管理制度如何影响土壤微生物组内部组织和稳健性的整合信息。
尽管取得了实质性进展，但在农业管理与微生物群落组装之间的界面上仍存在若干知识空白。关于施氮和耕作如何塑造土壤微生物组的详细机制见解主要来自单个研究站的长期试验，在这些试验中，肥料制度和土壤扰动的因子组合是在相对均质的土壤和气候条件下实施的（Duan等，2022；Mo等，2024；Srour等，2020）。这些试验提供了很高的内部有效性，但只采样了狭窄范围的土壤类型和管理历史，使得难以评估推断出的微生物组装“规则”在异质性农场环境中的普适性。景观尺度的调查和多地点实验表明，耕作强度、施肥制度和固有土壤属性（如质地和pH值）可以解释微生物群落结构的大部分变异（Kracmarova等，2022；Le Guillou等，2019），而基于网络的研究沿着土地利用和种植强度梯度揭示了微生物共现结构的强烈转变（Banerjee等，2019；Idbella等，2025；Karimi等，2019）。然而，相对较少的田间研究同时量化了许多农场的施氮梯度和对比耕作制度，并同时整合了细菌多样性、氮循环基因库和共现网络结构。这一差距在温带玉米生产系统中尤为明显，因为从受控的单站试验中推导出的预测规则在商业环境中常常失效，因为它们未能考虑商业农场固有的土壤气候异质性、土壤质地变化和历史管理遗留效应的混杂背景噪声。
本研究评估了施氮量和耕作强度对加拿大魁北克省南部84块商业化玉米田土壤细菌群落的结构、功能潜力和共现网络的影响。该区域网络涵盖了主要的耕作制度（常规翻耕CP、少耕RT和直播DD）、高度变化的土壤质地以及代表北美温带农业生态系统当前生产实践的施氮量。通过使用高通量16S rRNA基因测序、关键氮循环基因（nifH、amoA、nirK）的定量PCR以及基于重复共现模式的网络推断，研究人员整合了对土壤微生物组的分类学、功能和拓扑视角。研究人员并非假设源自受控站试验的组装规则直接转化为复杂的农场环境，而是明确验证了嵌套过滤模型在商业农场网络中的适用性。具体来说，研究人员评估了施氮量和耕作制度在多大程度上作为次级确定性过滤器作用于细菌多样性、功能基因库和共现架构，而其主要背景过滤器是土壤质地和pH值。研究人员假设，虽然内在土壤属性建立了基础的生态模板，但农业管理措施施加了次级生理和拓扑过滤器，这些过滤器主要通过功能基因重新分配和相互作用网络重组表现出来，而非广泛的分类学变化。通过在商业化生产田中测试这些多层级的组装规则，这项工作为农业生态系统诊断提供了一个操作框架，展示了拓扑和功能指标如何能够检测到仅依赖传统分类学指标无法发现的管理压力。
研究背景与问题
土壤健康是物理、化学和生物属性交互作用的涌现特性，而土壤微生物群落作为关键驱动者，其组成和功能对管理措施极为敏感。尽管已有大量研究探讨施氮和耕作对微生物的影响，但这些认知多源于单一研究站的长期定位试验，缺乏在异质性商业农场环境下的验证。现有知识空白在于：在土壤质地、pH等本底属性差异巨大的实际生产场景中，管理措施（施氮、耕作）作为“二级过滤器”究竟如何重塑微生物群落？传统的分类学多样性指标是否足以捕捉这些变化？为此，本研究在加拿大魁北克省南部的84块商业化玉米田开展大规模田间试验，旨在揭示土壤本底属性与管理措施对细菌群落的多层级过滤效应，特别关注功能基因丰度和共现网络拓扑结构是否能提供超越分类学丰富度的敏感诊断指标。
研究方法概述
研究于2021至2023年在加拿大魁北克省南部84个商业化玉米田进行，涵盖三种耕作制度（常规翻耕CP、少耕RT、直播DD）和五个施氮水平（0、50、100、150、200 kg N/ha）。采用16S rRNA基因扩增子测序获取细菌群落组成；利用定量PCR（qPCR）测定氮循环关键基因（nifH、amoA、nirK）的绝对丰度；通过SpiecEasi方法构建微生物共现网络，并采用自举重抽样策略进行统计检验；结合距离冗余分析（dbRDA）和Mantel检验评估环境变量对群落结构的解释度。
研究结果
3.1 测序产出与群落概况
共获得62.7百万条高质量双端序列，经质控后保留28.6百万条，鉴定出16,643个扩增子序列变异（ASV）。优势菌门为放线菌门Actinobacteriota（28.63%）和变形菌门Proteobacteria（21.02%），群落组成在不同处理间保持保守。
3.2 氮、耕作与土壤质地对比α多样性的影响
α多样性整体稳定。少耕（RT）的香农指数最高（5.72），显著高于直播（DD，5.66），常规翻耕（CP）居中（5.70）。施氮引起轻微非线性波动，未产生显著差异。土壤质地显著影响均匀度，中粗质地土壤均匀度高于细质地土壤。
3.3 β多样性与群落组成
β多样性分析显示，土壤质地是群落组成的最强驱动力（ANOSIM R = 0.2407），其次是耕作制度（R = 0.0828），而施氮无显著效应（R = 0.0014）。群落变异主要由物种周转（turnover，占95.53%）驱动，而非嵌套性（nestedness）。
3.4 差异分类组成
ANCOM-BC2分析未检测到耕作制度间的差异丰度类群。施氮仅导致少数孤立类群变化。土壤质地则引起多个分类水平的显著差异，其中粗质地（G3）相对于细质地（G1）富集了Oligoflexia，而Rubrobacteria、Chloroflexia等显著减少。
3.5 功能潜力与氮循环基因丰度
功能预测（Tax4Fun2）显示，土壤质地显著影响氮循环功能谱：粗质地土壤硝化和厌氧氨氧化潜力更高，细质地土壤则具有更高的异化硝酸盐还原（DNRA）和固氮潜力。耕作制度也改变了功能谱：常规翻耕（CP）的DNRA和反硝化潜力更高，而少耕（RT）和直播（DD）的硝化和厌氧氨氧化潜力更高。施氮无显著影响。功能冗余指数（rFRI）仅在土壤质地间存在显著差异。qPCR结果显示：高氮（200 kg N/ha）显著增加amoA丰度，降低nifH丰度，但不影响nirK；土壤质地仅显著影响amoA丰度（细质地最高）；耕作制度对所有三个基因均无显著影响。
3.6 微生物共现网络结构
网络分析揭示了管理措施的强烈拓扑过滤效应。耕作制度方面，直播（DD）网络最大、最模块化，少耕（RT）网络最小但集中性最高，常规翻耕（CP）居中但聚类系数最高。施氮方面，低氮（50 kg N/ha）导致网络更密集、更集中、模块性更低；高氮（150-200 kg N/ha）虽网络规模增大，但正关联边比例下降，负关联增加。土壤质地方面，细质地（G1）网络密度和集中性最高但规模最小，粗质地（G3）网络最大、模块性最高。
3.7 群落结构的环境驱动因素
dbRDA显示所选环境变量共同解释了27.9%的群落变异，其中土壤pH贡献最大（27.0%），其次为黏粒含量（18.8%）。施氮管理指标贡献小于0.5%。Mantel检验证实土壤pH与群落β多样性相关性最高（r = 0.385）。地理距离效应显著，但控制空间距离后环境因素仍独立解释群落变异。
讨论与结论
本研究在商业化农场网络中验证了嵌套过滤模型：土壤质地和pH构成初级环境过滤器，主导细菌群落组成；耕作制度作为次级物理过滤器，主要重塑共现网络拓扑而非分类学组成；施氮则作为生理过滤器，主要体现于功能基因丰度（amoA上调、nifH下调）和网络边缘特征（正关联减少）的变化，而对分类学组成影响微弱。研究表明，传统的α和β多样性指标不足以敏感反映短期管理压力，而功能基因丰度和网络拓扑结构提供了更灵敏的诊断工具。具体结论如下：土壤质地和pH是群落组成的主要决定因素，其影响覆盖了管理信号；耕作制度通过改变栖息地异质性，显著改变了微生物共现网络的规模和模块化程度，但对分类学组成影响有限；施氮通过改变资源供给，触发了功能基因的重新分配和网络内部关联模式的转变，表现为高氮下正关联减少和竞争性排斥加剧。这项研究为农业生态系统健康监测提供了操作框架，表明将拓扑和功能指标纳入常规土壤健康评估体系，能够更早、更准确地检测管理压力的生态后果。
