《Agriculture》：Land Use Type Governs Vertical Distribution, Co-Occurrence Network Structure, and Functions of Soil Bacterial Communities in an Alpine Semi-Arid Basin
编辑推荐：
为揭示土地利用类型如何调控高寒半干旱湟水河流域0–100 cm剖面土壤细菌群落结构、共现网络及碳氮功能，研究人员采用16S rRNA测序对耕地、草地、林地、湿地和撂荒地五种土地利用类型的土壤进行分析。假单胞菌门（Pseudomonadota）、放线菌门（Act
为揭示土地利用类型如何调控高寒半干旱湟水河流域0–100 cm剖面土壤细菌群落结构、共现网络及碳氮功能，研究人员采用16S rRNA测序对耕地、草地、林地、湿地和撂荒地五种土地利用类型的土壤进行分析。假单胞菌门（Pseudomonadota）、放线菌门（Actinomycetota）、绿弯菌门（Chloroflexota）、酸杆菌门（Acidobacteriota）和芽孢杆菌门（Bacillota）在所有群落中占优势（77.25–87.04%）。各土地利用类型内均发生垂直分层：富营养型类群主导表层（49.77–70.46%），而耐胁迫厌氧菌随深度增加而增多。土地利用类型是群落分异的首要驱动因子（解释56.46%的变异），草地明显分离。土壤有机碳（SOC）、全氮（TN）、全磷（TP）、有效磷（AP）和钾为关键驱动因子。网络分析显示，草地正相关比例最高（53.3%），湿地模块化指数最高（0.878），耕地负相关比例最高（55.3%）且模块化指数最低（0.774）。功能预测表明，草地富集甲烷氧化、亚硝酸盐呼吸和固氮功能；耕地以甲醇氧化和尿素水解为特征；林地和湿地表现出较强的芳香族化合物降解潜力。上述发现阐明了土地利用对垂直细菌分化和功能的调控作用，可为高寒半干旱流域的土地优化与土壤保育提供支撑。
研究背景方面，土壤微生物是陆地生态系统生物地球化学循环的核心驱动者，在维持土壤肥力、调控养分循环和支撑生态系统多功能性方面发挥关键作用。土地利用变化直接改变土壤理化性质，进而重塑微生物群落组成、多样性及种间互作格局，最终影响生态功能表达。高寒半干旱区生态环境脆弱，土壤生态系统对人为干扰高度敏感。已有研究多集中于湿润低地、农业或森林生态系统，且多针对表层土壤或单一土地利用转变，未能解析土地利用类型与土壤深度如何共同调控高寒半干旱流域的群落组成、互作网络和碳氮功能潜力。湟水河流域是黄河上游重要支流流域，现有研究主要集中在水环境问题，土壤微生物群落对不同土地利用类型的响应特征尚未系统研究。基于此，研究人员提出假设：土地利用类型和土壤深度共同调控土壤细菌群落组成与多样性，土地利用类型为群落分异的首要驱动因子；不同土地利用类型形成独特的共现网络结构和互作模式；碳氮循环功能潜力在土地利用类型间存在显著差异并随剖面变化。
研究人员在湟水河流域西宁段选取耕地（CL）、草地（GL）、人工林地（FL）、河滨湿地（WL）和撂荒地（AL）五种典型土地利用类型，每种设置3个独立重复样地，样地间距不小于500 m，于2023年9月采样。每个样地按五点混合法采集0–20、20–40、40–60、60–80和80–100 cm五个层次土样；因撂荒地土壤发育浅，仅采集0–20 cm层。共采集分析63个土壤样品。主要关键技术方法包括：采用电位法、重铬酸钾外加热法、碳氮分析仪、钼锑抗比色法、火焰光度法和EDTA滴定法等测定土壤理化性质；采用FastDNA Spin Kit for Soil提取土壤总DNA，以338F/806R引物扩增细菌16S rRNA基因V3–V4高变区，在Illumina MiSeq PE300平台进行高通量测序；使用Fastp进行质控、DADA2降噪生成ASV，基于SILVA v138数据库在QIIME2中完成分类注释；采用PCoA分析β多样性，LEfSe鉴定指示类群，以|r|>0.6且p<0.001构建共现网络并用Gephi可视化，采用FAPROTAX预测碳氮循环功能，统计分析在R软件中完成。
研究结果方面，3.1不同土地利用类型下土壤理化性质的垂直分布：土地利用类型显著改变土壤碳、氮、磷养分的含量及垂直分布格局，草地0–80 cm层土壤有机碳（SOC）最高（28.91–24.67 g/kg），显著高于其他类型；全氮（TN）呈相同垂直格局，草地和耕地表层TN分别为2.33和2.12 g/kg，为湿地和撂荒地的3.5–5.55倍；全磷（TP）和有效磷（AP）在草地整个剖面显著较高；所有土壤均呈碱性（pH>8）。各土地利用类型均表现养分随深度耗竭，但下降速率不同。
3.2不同土地利用类型下细菌群落组成与多样性变异：细菌α多样性受土地利用类型和土壤深度共同调控，表层丰富度和多样性总体较高，随深度波动而无一致单调下降趋势；林地和湿地各层α多样性总体较低。门水平上群落呈显著水平分异和垂直分层，五个核心门占总序列77.25–87.04%；表层以假单胞菌门和放线菌门为主，二者合计占49.77–70.46%；随深度增加，绿弯菌门在地下层次富集，芽孢杆菌门在草地40–80 cm层显著富集并超过放线菌门成为第二优势门。纲水平上各优势门内部结构因土地利用类型而异，草地Clostridia明显高于其他类型。层次聚类显示草地20–100 cm样品独立成簇，凸显其群落结构的独特性；属水平上Sphingomonas为耕地、林地、湿地和撂荒地的核心优势属，草地优势属呈显著垂直分化，0–20 cm以Gaiella为主，20–80 cm以Chryseolinea、Aggregatilinea、Clostridium和Terrisporobacter为主。PCoA前两轴累计解释56.46%的群落变异（PC1为40.05%，PC2为16.41%），草地明显分离，耕地、林地和湿地部分重叠。Venn图显示表层五种类型仅共享22个ASV，全剖面四种类型共享245个ASV。LEfSe在LDA阈值3下共鉴定35个显著差异细菌类群，各土地利用类型具有独特生物标志类群，草地Clostridia的LDA得分超过4.5。
3.3共现网络结构与互作模式：各土地利用类型细菌共现网络结构和拓扑参数呈系统差异。林地节点数最多（224），草地最少（142），较林地减少36.6%；湿地边数最多（1006），草地最少（614）；湿地平均度最高（9.863），林地平均度最低（8.527）、平均路径长度最长（0.380）；聚类系数差异较小（0.951–0.958）。除草地外，耕地、林地和湿地负边比例均高于正边（分别为55.3%、52.8%和53.3%）；耕地模块化指数最小（0.774），湿地最大（0.878）；草地平均路径长度最短（0.323）。
3.4土壤细菌群落结构的环境驱动因子：门水平上绿弯菌门与SOC、TN、TP、AP和K+显著正相关（r=0.50至0.61），放线菌门和芽单胞菌门与这些养分指标显著负相关（r=?0.51至?0.62）；芽孢杆菌门仅与pH和溶解性有机碳（DOC）正相关，黏球菌门与多数养分指标显著正相关。属水平上优势属按养分响应分为两组：Clostridium、Aggregatilinea、Chryseolinea、Terrisporobacter和Devosia与养分显著正相关，其中Clostridium与氮磷比显著正相关（r=0.54）；Skermanella、Blastococcus、Solirubrobacter和Agromyces与养分显著负相关；Sphingomonas与TP、pH、DOC和AP负相关，Lysobacter与pH极显著负相关。
3.5土地利用驱动的细菌碳氮循环功能转变：基于FAPROTAX的功能预测显示四种土地利用类型功能谱差异显著，以碳氮循环相关功能为主，好氧化能异养和化能异养为最普遍的核心功能。碳循环方面，草地显著富集发酵、烃类降解、甲烷氧化和甲基营养，甲烷氧化与甲基营养相对丰度呈草地>耕地>林地>湿地的顺序；耕地甲醇氧化能力显著高于林地和湿地；芳香族化合物降解在林地和湿地整个剖面突出。氮循环方面，草地亚硝酸盐呼吸、固氮和氮呼吸显著高于其他类型；硝酸盐还原在各生境广泛分布；尿素水解主要富集于耕地各层和湿地表层。垂直格局上，甲烷氧化和甲基营养分别在草地和耕地0–40 cm显著富集并随深度下降；耕地烃类降解在全剖面维持较高丰度；草地表层亚硝酸盐还原显著较高并随深度明显下降；湿地硝酸盐还原沿剖面呈“先降后升”格局。
讨论部分，4.1土地利用与土壤深度共同调控土壤细菌群落的垂直分化：表层仅22个共享ASV表明表层细菌群落具有强土地利用特异性；草地全剖面ASV总数最高（4282），深层丰富度显著较高，可能与草地植物深根系通过根系周转和分泌物向深层输送有机碳有关；林地和湿地α多样性较低，分别与难降解凋落物限制生态位分化和长期淹水还原条件对好氧细菌的生理胁迫有关。门水平上假单胞菌门在耕地表层丰度最高（42.62%），随深度先降后升并在40–60 cm出现次峰；放线菌门在湿地和草地表层占优并主要富集于20–60 cm中层；绿弯菌门在地下40–80 cm增加；芽孢杆菌门在草地40–80 cm显著富集，归因于其产芽孢厌氧类群的强胁迫耐受性。属水平垂直更替揭示生态位分化，草地由表层Gaiella向深层厌氧、耐寡营养类群转变，反映从表层好氧高底物向深层缺氧低底物环境的连续过渡。研究指出，评估土地利用对土壤微生物的影响必须纳入垂直深度梯度。
4.2土地利用驱动的土壤细菌共现网络格局与结构分化：草地是唯一以协同共生为主的类型（正相关占53.3%），节点和边数最少，网络结构相对简单；耕地、林地和湿地以种间竞争为主，负相关比例均超过52%，耕地受长期耕作施肥影响负边比例最高（55.3%）且模块化最低（0.774）；湿地节点（204）、边（1006）和模块化指数（0.878）最高，网络更具模块化和结构稳健性；林地节点最多（224）、边955，物种组成更丰富。
4.3环境因子驱动土壤细菌群落结构分化：SOC、TN、TP、AP和K+是核心驱动因子，不同细菌类群响应模式不同。绿弯菌门与多数养分正相关而放线菌门负相关的对立格局，主要由底物质量（易分解与矿物结合态有机碳）、深度相关氧可利用性及富营养类群的竞争排斥共同决定，而非单纯由富营养/寡营养策略决定。属水平上，与养分正相关的多为厌氧或兼性厌氧发酵细菌，分布于耕地和草地表层等肥沃土壤；与养分负相关的多为寡营养或耐胁迫类群，可在贫瘠或碱性土壤中占据优势生态位。
4.4土地利用与土壤深度对土壤细菌碳氮循环功能的调控：土地利用类型与土壤深度通过协同改变底物供应、氧化还原条件和垂直养分分布，共同调控碳氮循环功能的组成与丰度，形成水平分异与垂直分层的双重功能格局。好氧化能异养和化能异养在所有类型中占优且随深度显著下降。碳循环方面，草地以大气低浓度甲烷好氧氧化为主且深层较稳定，湿地以厌氧甲烷氧化为主而FAPROTAX仅预测好氧途径故预测值较低；耕地甲醇氧化在表层最突出，与长期有机肥投入和秸秆还田释放的挥发性底物有关；林地和湿地芳香族化合物降解占优。氮循环方面，草地亚硝酸盐呼吸和固氮在中上层富集，与适中C/N比和微好氧—好氧交替环境有关；耕地高氮输入使尿素水解全层占优并显著抑制固氮；湿地长期厌氧及林地酸性条件和难降解底物限制关键氮循环功能表达。不同生态系统由此形成差异化生态服务能力：草地支撑甲烷氧化与生物固氮，耕地支撑甲醇氧化与尿素水解，湿地支撑硝酸盐还原，林地和湿地支撑芳香族化合物降解。
研究结论部分译为：本研究系统阐明了高寒半干旱湟水河流域五种土地利用类型下土壤细菌群落的垂直分化特征、共现网络拓扑及碳氮循环功能潜力。结果表明，土壤细菌群落结构受土地利用类型和土壤深度共同调控。假单胞菌门和放线菌门呈现显著分异的垂直分布格局，而草地深层富集芽孢杆菌门和Clostridium等厌氧类群，清晰呈现由表层富营养型类群向深层耐胁迫类群的连续垂直演替。土壤有机碳、全氮、全磷、有效磷和钾离子是驱动细菌群落分异的核心环境因子，其通过筛选具有不同养分适应策略的微生物类群，塑造了资源依赖型群落分布格局。共现网络分析显示各土地利用类型网络结构存在系统差异。草地是唯一以微生物协同为主的类型，湿地网络模块化（结构稳健性）最高，耕地以最强烈的种间竞争（负边比例最高）为特征。功能预测证实各土地利用类型功能谱显著分化：草地显著富集甲烷氧化、亚硝酸盐呼吸和生物固氮；耕地以甲醇氧化和尿素水解为主；林地和湿地表现出较强的芳香族化合物降解潜力。
