环境科学  2024, Vol. 45 Issue (6): 3605-3613   PDF    
宁夏黄土丘陵区农田土壤细菌海拔分布特征
吴霞1,2, 蔡进军1,2, 王长军1, 李维倩1, 陈刚1,2, 白阳阳1     
1. 宁夏农林科学院农业资源与环境研究所, 银川 750002;
2. 国家农业环境银川观测实验站, 银川 750002
摘要: 研究农田土壤微生物多样性、结构与功能的海拔分布特征及其影响因子, 对保护农田生态系统生物多样性具有重要意义. 利用高通量测序技术对宁夏黄土丘陵区农田土壤细菌多样性、群落结构和代谢功能沿海拔分布特征及其对土壤理化性质的响应进行分析. 结果表明:①土壤细菌Alpha多样性指数与海拔呈显著负相关, 沿海拔呈现先降低再略升高的变化趋势. ②变形菌(Proteobacteria)、放线菌(Actinobacteria)和酸杆菌(Acidobacteria)等7个菌门为优势菌群, 其中5个菌群在海拔间具有极显著性差异. ③二级分类水平上包括膜运输、碳水化合物代谢和氨基酸代谢等36个细菌代谢功能, 其中22个在不同海拔间具有显著性差异, 12个具有极显著性差异. ④Pearson相关分析表明, 土壤含水量、容重、pH和碳氮比对细菌Alpha多样性影响最为显著, 总有机碳、全氮和全磷等土壤养分对细菌Beta多样性具有显著影响. ⑤Mantel分析表明, 门水平细菌群落结构受土壤含水量、总有机碳和碳氮比影响, 细菌代谢功能与土壤pH、总有机碳、全氮、全磷和碳氮比均有显著相关性. 方差分解分析显示土壤含水量对群落结构解释度最高, pH对代谢功能解释度最高. 由此可知, 土壤含水量和pH是影响宁夏黄土丘陵区农田土壤细菌多样性、群落组成和代谢功能的主要土壤因子.
关键词: 细菌多样性      群落结构      代谢功能      海拔分布      农田土壤      黄土丘陵区     
Altitude Distribution Characteristics of Farmland Soil Bacteria in Loess Hilly Region of Ningxia
WU Xia1,2 , CAI Jin-jun1,2 , WANG Zhang-jun1 , LI Wei-qian1 , CHEN Gang1,2 , BAI Yang-yang1     
1. Institute of Agricultural Resources and Environment, Ningxia Academy of Agriculture and Forestry Science, Yinchuan 750002, China;
2. Station of Observation and Experiment National Agricultural Environment in Yinchuan, Yinchuan 750002, China
Abstract: It is of great significance for the conservation of biodiversity in farmland ecosystems to study the diversity, structure, functions, and biogeographical distribution of soil microbes in farmland and their influencing factors. High-throughput sequencing technology was used to analyze the distribution characteristics of soil bacterial diversity, community structure, and metabolic function along elevation and their responses to soil physicochemical properties in farmland in the loess hilly areas of Ningxia. The results showed that: ① The Alpha diversity index of soil bacterial was significantly negatively correlated with elevation (P < 0.05) and showed a trend of decreasing and then slightly increasing along the elevation. ② Seven phyla, including Proteobacteria, Actinobacteria, and Acidobacteria, were the dominant groups, and five of them showed highly significant differences between altitudes (P < 0.01). ③ At the secondary classification level, there were 36 metabolic functions of bacteria, including membrane transport, carbohydrate metabolism, and amino acid metabolism, of which 22 showed significant differences, and 12 showed extremely significant differences among different altitudes. ④ Pearson correlation analysis showed that soil water content, bulk density, pH, and carbon-nitrogen ratio had the most significant effects on bacterial Alpha diversity, whereas soil nutrients such as total organic carbon, total nitrogen, and total phosphorus had significant effects on bacterial Beta diversity. ⑤ Mantel test analysis showed that the soil water content, total organic carbon, and carbon-nitrogen ratio affected bacterial community structure at the phylum level, and soil pH, total organic carbon, total nitrogen, total phosphorus, and carbon-nitrogen ratio were significantly correlated with bacterial metabolic function. Variance partitioning analysis showed that soil water content had the highest explanation for the community structure of soil bacteria, whereas soil pH had the highest explanation for metabolic function. In conclusion, soil water content and pH were the main factors affecting the diversity, community composition, and metabolic function of soil bacteria in farmland in the loess hilly region of Ningxia.
Key words: bacterial diversity      community structure      metabolic function      altitude distribution      farmland soil      loess hilly region     

土壤微生物受到多个因素的影响, 其海拔分布存在一定的复杂性和异质性[1 ~ 4]. 一般来讲, 随着海拔的升高, 气温逐渐降低, 同时降雨量和空气湿度也随之减少, 以上自然环境的变化对土壤微生物群落产生了影响[5 ~ 7]. 有研究发现土壤微生物数量和多样性随海拔升高而递减, 低海拔地区的微生物群落比高海拔地区更为多样化和复杂化, 这可能是由于低海拔地区气候条件更适宜微生物生长和繁殖导致的[8 ~ 10]. 在高海拔地区, 由于环境寒冷和干旱, 土壤微生物数量相对较少, 但在某些情况下, 高海拔地区的土壤微生物可能具有更高的适应性和抗逆能力[11 ~ 13]. 近年来的研究发现, 在不同山区中, 土壤微生物的分布情况还会受到其他因素的影响, 例如植被种类、土壤类型和地形地貌等[12 ~ 15]. 在农田生态系统中, 土壤微生物参与了养分循环、作物生长调节、病虫害防治和污染物降解等多个过程, 是农业生产中不可或缺的组成部分, 对土壤健康和生态系统稳定有着重要作用, 其分布特征不仅受自然环境影响, 还受到人为干扰和制约[16 ~ 19]. 有研究发现, 机械扰动、密集灌溉、过度施肥和农药不合理使用等导致土壤微生物数量和多样性降低[20 ~ 23], 而少耕、免耕、间作和覆盖等保护性耕作措施可对土壤微生物多样性和生态系统功能产生积极影响[24 ~ 27]. 综上, 研究农田土壤微生物的分布格局及其驱动因子, 不仅对保护农田土壤微生物多样性和土壤可持续利用具有重要意义, 而且对于更好地保护和解决生态环境失衡问题也十分必要.

黄土丘陵区是中国重要的农业区之一, 海拔高度差异大, 因此不同海拔高度下的土壤微生物数量和种类也会有所不同. 近年来, 已有众多学者针对黄土丘陵区土壤微生物多样性和群落结构等开展了研究, 结果表明植物类型、种植结构、海拔高度和土壤性质等均对土壤微生物具有重要影响[28 ~ 32]. 但以上研究多以小尺度为主, 针对区域尺度上土壤微生物分布特征的研究较少. 本文以宁夏固原市黄土丘陵区农田土壤为研究对象, 分析了不同海拔土壤细菌多样性、群落结构、代谢功能及与土壤理化性质关系, 以期为黄土丘陵区土壤微生物多样性的认识和农田土壤健康保护提供依据.

1 材料与方法 1.1 研究区概况

固原市位于宁夏回族自治区南部, 地处黄土高原西北边缘, 地势起伏较大, 海拔多在1 320 ~ 2 928 m之间, 以黄土丘陵为主要地貌, 沟壑纵横, 地形破碎, 水蚀强烈. 该地区处于暖温带半干旱气候区, 年平均气温6 ~ 8℃, 年降水量400 ~ 650 mm, 干燥度1.55 ~ 1.30, 大于10℃积温2 200 ~ 3 000℃, 无霜期120 ~ 180 d[24]. 固原市是宁夏旱作耕地集中分布区, 截至2019年耕地面积3.302 × 105 hm2, 占固原市总面积的31.39%, 主要以坡耕地和旱作梯田为主, 土壤贫瘠, 耕地质量等级总体较低, 主要种植马铃薯、玉米、小麦和小杂粮等农作物.

1.2 样品采集

2022年8底至9月初, 在农作物收获期对固原市黄土丘陵区不同海拔农田土壤进行调查采样. 每个调查采样田块按照梅花法采集10个混样点进行混合, 采样深度0 ~ 20 cm. 将所有混样点剔除根系后充分混匀, 用无菌自封袋装取50 g左右土样立即放入冷藏箱保存, 运回实验室后存于-80℃冰箱, 用于细菌高通量测序. 剩余土样全部装入大自封袋, 带回实验室风干, 过2 mm筛后用于土壤理化性质测试, 采样同时用GPS记录田块中心点坐标和海拔. 研究区共调查采集121个农田土壤样点, 按100 m海拔间隔划分为6组, 样点分布如图 1所示.

图 1 样点分布示意 Fig. 1 Distribution of sampling sites

1.3 分析方法 1.3.1 土壤理化性质测定

土壤含水量(SWC)采用铝盒烘干法测定;容重(SBD)采用环刀法测定;土壤pH采用电位法(水土比为2.5∶1)测定;总有机碳(TOC)采用重铬酸钾容量法测定;全氮(TN)采用半微量凯氏定氮法测定;全磷(TP)采用氢氧化钠碱熔-钼锑抗比色法测定[33];碳氮比(C/N)为土壤总有机碳与全氮比值.

1.3.2 土壤细菌高通量测序

土壤样品DNA使用HiPure Soil DNA试剂盒(Magen, Guangzhou, China)提取, PCR扩增引物为341F(5′-CCTACGGGNGGCWGCAG-3′)和806R(5′-GGACTACHVGGGTATCTAAT-3′), 扩增区域为V3~V4区, 扩增产物分离纯化采用琼脂糖凝胶电泳法. 细菌高通量测序采用Illumina平台(PE250), 原始测序序列采用Usearch方法进行拼接、过滤, 按≥97%相似度聚类为操作分类单元(OTU), 比对NCBI 16S rDNA数据库进行物种分类注释, 按最小样本序列抽平至28 149条OTU进行后续土壤细菌多样性、群落组成和代谢功能分析.

1.4 数据处理

采用Excel进行基础数据整理, Alpha多样性指数采用QIIME软件(版本1.9.1)计算, 差异性检验、非度量多维度分析(NMDS分析)和方差分解分析(variance partitioning analysis, VPA)使用R语言的vegan包计算, 环境因子与多样性指数和物种之间的Pearson相关系数采用R语言的psych包计算, 绘图使用R语言的ggplot2包.

2 结果与分析 2.1 土壤理化性质海拔变化特征

单因素方差分析表明(表 1), 研究区农田土壤含水量随海拔升高先降低再升高, 低海拔农田(1 600 ~ 1 800 m)土壤含水量显著高于高海拔(1 900 ~ 2 100 m)农田土壤含水量, 最高为17.04%(1 600 m), 最低为11.54%(1 900 m). 不同海拔土壤容重差异不显著, 平均为1.29 g·cm-3, 无明显海拔变化规律. 土壤pH整体呈碱性, 1 700 m处最高(8.58), 2 100 m处最低(8.23), 沿海拔升高有明显的降低趋势. 土壤总有机碳、全氮和全磷含量总体上均随海拔升高而显著增大. 其中, 2 100 m处总有机碳和全氮含量最高, 分别为含量最低的1 600 m处的1.84倍和1.72倍, 2 000 m和2 100 m处全磷含量较最低的1 700 m处高27.59%. 土壤碳氮比在不同海拔间差异不显著, 且无明显的海拔分布规律.

表 1 不同海拔土壤理化性质1) Table 1 Soil physicochemical properties

2.2 土壤细菌海拔分布特征 2.2.1 土壤细菌多样性海拔分布特征

Alpha多样性反映单个样本内的物种多样性, 研究区农田土壤细菌Alpha多样性指数沿海拔分布状况如图 2所示, 包括反映物种丰富度的Sobs指数、反映物种均匀度的Pielou指数、综合反映物种丰富度与均匀度的Shannon指数和反映物种谱系多样性的PD指数. 整体来看, 研究区农田土壤细菌Sobs指数介于2 617 ~ 3 710之间, 均值为3 244, Pielou指数介于0.75 ~ 0.87之间, 均值为0.83, Shannon指数介于8.54 ~ 10.28之间, 均值为9.67, PD指数介于287.85 ~ 448.56之间, 均值为381.59. Pearson相关分析表明4种多样性指数均与海拔呈显著性负相关关系, 均沿海拔呈现先降低再略升高的变化趋势. Kruskal-Wallis检验显示研究区农田土壤细菌4种Alpha多样性指数在不同海拔间的差异显著. 其中, 1 600 ~ 2 000 m区间细菌多样性呈降低趋势, 1 600 m海拔土壤细菌丰富度、均匀度、物种多样性和谱系多样性最高, 显著或极显著高于高海拔地区, 2 000 m海拔的细菌均匀度和物种多样性最低, 1 900 m海拔的细菌谱系多样性最低, 均显著低于其他海拔, 2 100 m海拔细菌丰富度、均匀度和多样性较2 000 m均略有升高.

灰色虚线为多样性指数均值;红色曲线为样点拟合曲线, 灰色阴影为拟合曲线95%置信区间;*表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 2 不同海拔农田土壤细菌Alpha多样性指数 Fig. 2 Alpha diversity index of soil bacteria in farmland at different altitudes

Beta多样性指数反映样本间的物种组成差异, 本文基于细菌OTU丰度, 采用Bray-curtis距离的NMDS分析研究区农田土壤细菌群落结构间的差异性[34], 结果如图 3(a)所示. NDMS分析的Stress = 0.154(< 0.2), 说明海拔梯度对研究区农田土壤细菌群落结构具有一定的影响, 但从图中的样点分布和置信区间来看, 不同海拔之间的细菌群落结构相似性较大, 差异并不明显. 进一步采用Anosim法检验不同海拔间的细菌差异性, 结果表明细菌群落结构沿海拔梯度的差异性极为显著(P = 0.001), 但差异性较小(R = 0.105). 同时, 海拔与NDMS分析第一轴(NDMS1)的回归分析也显示[图 3(b)], 海拔与细菌群落结构之间的关系及其显著(P = 0.001), 但拟合度较低(R2 = 0.08).

(a)中不同颜色椭圆为不同海拔样点95%置信区间;(b)中黑色曲线为样点拟合曲线, 灰色阴影为拟合曲线95%置信区间 图 3 不同海拔农田土壤细菌Beta多样性指数 Fig. 3 Beta diversity index of soil bacteria in farmland at different altitudes

2.2.2 土壤细菌群落组成海拔分布特征

研究区农田土壤细菌在门分类水平上可划分为41门, 相对丰度平均值大于5%的细菌门包括:变形菌(Proteobacteria, 21.29%)、放线菌(Actinobacteria, 16.21%)、酸杆菌(Acidobacteria, 14.87%)、浮霉菌(Planctomycetes, 10.84%)、拟杆菌(Bacteroidetes, 10.00%)、绿弯菌(Chloroflexi, 7.01%)和芽单胞菌(Gemmatimonadota, 5.46%), 合计占比为85.69%, 是研究区农田土壤优势菌(图 4). Kruskal-Walli检验结果显示, 变形菌、放线菌、酸杆菌、浮霉菌和芽单胞菌在不同海拔间具有显著差异. 其中, 变形菌的相对丰度在1 600m最高(24.63%), 2 000 m最低(18.16%), 且沿着海拔呈先下降再升高的分布特征, 放线菌则随着海拔先升高再下降, 1 600 m最低(11.42%), 1 900 m最高(17.43%);酸杆菌、浮霉菌和芽单胞菌则沿海拔呈波动式分布, 酸杆菌在1 600 m(17.11%)最高, 浮霉菌在2 000 m(12.14%)最高, 芽单胞菌在1 800 m(6.04%)最高, 3种细菌均在1 700 m的相对丰度最低, 分别为13.22%、9.29%和4.77%(图 4).

E1、E2、E3、E4、E5和E6分别表示1 600、1 700、1 800、1 900、2 000和2 100m海拔;不同小写字母表示P < 0.05时不同海拔的差异显著 图 4 不同海拔农田土壤门水平细菌群落组成 Fig. 4 Composition of soil bacterial communities at the phylum level in farmland soil at different altitudes

2.2.3 土壤细菌功能海拔差异分析

采用Tax4Fun[35]工具对研究区农田土壤细菌功能进行预测, 二级分类水平上共鉴定出36种代谢功能, 其中相对丰度 > 1% 的代谢功能19个, 在不同海拔间的功能差异如图 5所示. 平均相对丰度较高的有:膜运输(membrane transport)(12.78%)、碳水化合物代谢(carbohydrate metabolism)(12.73%)、氨基酸代谢(amino acid metabolism)(12.59%)、能量代谢(energy metabolism)(7.13%)、辅因子和维生素代谢(metabolism of cofactors and vitamins)(7.05%)以及信号转导(signal transduction)(7.00%). Kruskal-Wallis检验显示36个代谢功能中在不同海拔间具有显著差异的功能达22个, 极显著差异的达12个, 其中, 碳水化合物代谢、聚糖生物合成和代谢(glycan biosynthesis and metabolism)、异生素生物降解和代谢(xenobiotics biodegradation and metabolism)等功能差异性最为显著.

*表示 P < 0. 05,**表示 P < 0. 01,***表示 P < 0. 001 图 5 不同海拔农田土壤细菌代谢功能差异 Fig. 5 Differences in metabolic function of soil bacteria in farmland at different altitudes

2.3 土壤细菌海拔分布格局影响因子分析 2.3.1 土壤理化性质对细菌多样性的影响

通过Pearson相关分析表明:土壤理化性质对细菌多样性具有显著性影响. 土壤含水量、容重、pH和碳氮比对细菌Alpha多样性影响最为显著. 其中, 土壤含水量与Pielou指数、Shannon指数和PD指数均呈显著正相关, pH与Sobs指数、Shannon指数和PD指数均呈显著正相关, 而土壤容重则与Pielou指数和Shannon指数呈显著负相关, 土壤碳氮比仅与PD指数呈显著负相关. 土壤总有机碳、全氮、全磷和碳氮比与NMDS分析第一轴(NMDS1)呈极显著正相关, pH则与NMDS1呈极显著负相关(表 2).

表 2 细菌多样性和土壤理化性质Pearson相关系数1) Table 2 Pearson correlation between bacterial diversity and soil physicochemical properties

2.3.2 土壤理化性质对细菌群落结构和功能的影响

基于研究区土壤细菌物种丰度和功能预测丰度, 采用Mantel检验探讨细菌群落结构和代谢功能与土壤理化性质之间的关联性. 结果显示[图 6(a)], 研究区农田土壤细菌门水平群落结构与土壤含水量呈极显著正相关, 与总有机碳和碳氮比呈显著正相关, 但相关性均较弱(r < 0.3), 细菌代谢功能与土壤pH、总有机碳和全氮呈极显著弱相关, 与全磷和碳氮比呈显著弱相关. 由此可知, 除土壤容重外所选土壤因子均对细菌群落结构和代谢功能具有显著影响. 进一步采用VPA法分析上述土壤因子对门类细菌群落结构和代谢功能差异的解释度, 结果如图 6(b)所示. 土壤含水量对细菌群落结构变化的解释度最高, 且对细菌群落结构的独立解释度和共同解释度均达到极显著水平(P < 0.001), 总有机碳和容重对细菌群落结构的解释度次之, 但远低于土壤含水量. pH对细菌代谢功能的解释度远高于其他因子, 其次为全磷、全氮和总有机碳, 上述因子除全磷外对细菌代谢功能的独立解释度和共同解释度均达到显著水平(P < 0.05).

*表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 6 土壤理化性质与细菌群落结构、代谢功能相关性网络和解释度 Fig. 6 Correlation network between soil physicochemical properties and bacterial community structure, metabolic function and explanation degree diagram

3 讨论 3.1 不同海拔农田土壤理化性质分布特征

海拔梯度是影响区域气候、植被、土壤和微生物变化的重要地形因子[3]. 研究区农田土壤含水量随海拔呈先下降再升高的分布特征, 1 900 m含水量最低, 仅为11.54%, 1 600 ~1 900 m径流较大且蒸发强烈, 导致土壤水分随海拔升高而降低, 这与Qiu等[36]对黄土高原土壤水分随海拔降低的研究结论相同. 而1 900~2 100 m温度较低且林草覆盖率较高, 蒸发有所减弱, 因此农田土壤含水量略有升高. 一般而言, 随着海拔高度的增加, 温度降低, 有机物分解较为缓慢, 与低海拔相比较易积累有机质, 导致表层土壤中有机碳随海拔升高而增加[37], 研究区农田土壤总有机碳含量随海拔升高而显著增加进一步证实了这一结论. 有机碳的变化会直接影响土壤pH和养分含量等[38], 相关分析显示[图 5(a)], 土壤总有机碳与pH显著负相关, 与全氮和全磷显著正相关, 因此研究区农田土壤呈现出pH随海拔升高降低、全氮和全磷含量随海拔升高增加的变化趋势.

3.2 不同海拔农田土壤细菌多样性、群落结构及代谢功能分布特征

本研究中低海拔农田土壤细菌Alpha多样性指数显著高于高海拔农田, 且随海拔升高呈先降低再略升高的变化趋势, 同时细菌群落结构间的Beta多样性亦表现出显著的海拔差异. 可能是由于海拔的升高导致温度降低, 土壤微生物生长和繁殖受到抑制, 其数量和种类逐渐减少, 因此1 600~1 900 m土壤细菌多样性随海拔升高而递减. 2 000~2 100 m海拔可能受土壤含水量和有机质增加影响, 细菌多样性略有升高. Li等[39]、Singh等[40]和Wang等[41]研究发现土壤微生物随海拔呈下凹趋势, 这与本研究的结果相似.

门分类水平上, 变形菌门、放线菌门、酸杆菌门和浮霉菌门相对丰度均大于10%, 是宁夏黄土丘陵区农田土壤的优势菌群, 且4种菌群在不同海拔间的差异极为显著. 土壤细菌在有机物分解、养分循环和土壤修复等生态过程中起着重要作用[18], 先前研究已证明变形菌、放线菌和酸杆菌的相对丰度在土壤变暖条件下增加[42]. 本研究中低海拔农田土壤变形菌丰度较高且随海拔升高显著递减, 可能正是由于低海拔温度较高, 有利于变形菌对土壤有机质矿化作用的发挥. 而放线菌随海拔的变化则与变形菌相反, 这与前人研究的结果并不一致. 放线菌对土壤中的各种难溶性惰性有机物具有更强的矿化分解作用[43], 而研究区高海拔农田土壤受低温影响有机碳活性降低, 从而需要更多放线菌来分解这部分有机物, 这可能是放线菌随海拔增加的因素之一, 具体机制有待进一步研究.

王颖等[44]对青藏高原高寒草甸土、吴则焰等[45]对中亚热带森林土壤和韩冬雪等[46]对长白山红松林土壤的研究均发现多种土壤微生物代谢功能具有垂直地带性差异. 本研究中膜运输、碳水化合物代谢、氨基酸代谢是研究区农田土壤的主要代谢途径, 差异性检验显示近2/3代谢功能在不同海拔间具有显著差异, 进一步说明海拔对土壤细菌功能的发挥至关重要, 这可能与海拔差异引起的温度、水分以及土壤养分变化等有关.

3.3 土壤细菌多样性、群落结构及代谢功能与理化性质的关系

土壤、气候和植被是驱动土壤微生物群落时空变化的主要因素[3], 在气候和植被差异较小的区域农田尺度上, 土壤理化性质对微生物多样性及群落结构分布的影响尤为显著. 有研究表明土壤pH是影响土壤微生物时空分布的重要因素[5, 20, 39, 47], Lauber等[48]研究发现pH与细菌谱系多样性和群落结构显著相关, Li等[34]研究证明pH是驱动贡嘎山土壤细菌海拔格局的决定因素. 本研究发现pH同时与土壤细菌Alpha多样性和Beta多样性显著相关(表 2), 且对细菌代谢功能的解释度最高[图 6(b)], 进一步证明了pH对细菌多样性和代谢功能的重要影响. 宋贤冲[49]和林春英等[50]研究发现土壤水分是微生物多样性的主要驱动因子, 本研究中VPA分析表明土壤含水量是细菌群落结构差异的主要驱动因素, 相关分析也显示土壤含水量与细菌Pielou指数和Shannon指数相关性最高. 由此可见, pH和土壤含水量是驱动研究区农田土壤细菌海拔分布差异的关键因子. 此外, 土壤总有机碳、全氮、全磷、容重和碳氮比亦对研究区农田土壤细菌群落具有一定影响.

4 结论

(1)研究区农田土壤理化性质与海拔呈显著相关, 含水量随海拔先降低再升高, pH随海拔升高逐渐降低, 总有机碳、全氮和全磷含量随海拔升高显著增大.

(2)研究区农田土壤细菌多样性随着海拔先降低再略有升高, 与海拔呈显著负相关, 其中1 600 m多样性指数最高, 1 900~2 000 m最低.

(3)变形菌、放线菌、酸杆菌和浮霉菌是研究区农田土壤优势菌群;细菌代谢功能相对丰度较高的有膜运输、碳水化合物代谢和氨基酸代谢, 近2/3的代谢功能在海拔间具有显著差异.

(4)土壤含水量和pH是研究区农田土壤细菌海拔分布差异的主要影响因子, 而有机碳、全氮、全磷和碳氮比对农田土壤细菌群落结构和代谢功能亦具有一定影响.

参考文献
[1] 褚海燕, 王艳芬, 时玉, 等. 土壤微生物生物地理学研究现状与发展态势[J]. 中国科学院院刊, 2017, 32(6): 585-592.
Chu H Y, Wang Y F, Shi Y, et al. Current status and development trend of soil microbial biogeography[J]. Bulletin of Chinese Academy of Sciences, 2017, 32(6): 585-592.
[2] Fierer N, Jackson R B. The diversity and biogeography of soil bacterial communities[J]. Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(3): 626-631.
[3] 厉桂香, 马克明. 土壤微生物多样性海拔格局研究进展[J]. 生态学报, 2018, 38(5): 1521-1529.
Li G X, Ma K M. Progress in the study of elevational patterns of soil microbial diversity[J]. Acta Ecologica Sinica, 2018, 38(5): 1521-1529.
[4] Chu H Y, Gao G F, Ma Y Y, et al. Soil microbial biogeography in a changing world: recent advances and future perspectives[J]. mSystems, 2020, 5(2). DOI:10.1128/mSystems.00803-19
[5] Bryant J A, Lamanna C, Morlon H, et al. Microbes on mountainsides: contrasting elevational patterns of bacterial and plant diversity[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(S1): 11505-11511.
[6] 王国敏, 曹嘉瑜, 倪健. 山地土壤微生物地理分布格局及其驱动机制[J]. 地球与环境, 2019, 47(4): 565-574.
Wang G M, Cao J Y, Ni J. Geographical distribution pattern of soil microorganisms in mountains and its driving mechanism[J]. Earth and Environment, 2019, 47(4): 565-574.
[7] Frindte K, Pape R, Werner K, et al. Temperature and soil moisture control microbial community composition in an arctic⁃alpine ecosystem along elevational and micro-topographic gradients[J]. The ISME Journal, 2019, 13(8): 2031-2043. DOI:10.1038/s41396-019-0409-9
[8] Tang M Z, Li L, Wang X L, et al. Elevational is the main factor controlling the soil microbial community structure in alpine tundra of the Changbai Mountain[J]. Scientific Reports, 2020, 10(1). DOI:10.1038/s41598-020-69441-w
[9] Nottingham A T, Fierer N, Turner B L, et al. Microbes follow Humboldt: temperature drives plant and soil microbial diversity patterns from the Amazon to the Andes[J]. Ecology, 2018, 99(11): 2455-2466. DOI:10.1002/ecy.2482
[10] 吴狄, 杨腾, 林臣峰, 等. 天目山银杏根区土壤固氮微生物群落沿海拔梯度的分布特征[J]. 土壤, 2022, 54(5): 958-967.
Wu D, Yang T, Lin C F, et al. Elevational distribution of soil diazotrophic community in root zone of ginkgo biloba in Tianmu Mountain[J]. Soils, 2022, 54(5): 958-967.
[11] 李宏林, 周浪敏, 薛静, 等. 高寒沼泽化草甸群落常见种对土壤水分条件的响应研究[J]. 生态科学, 2023, 42(4): 1-8.
Li H L, Zhou L M, Xue J, et al. Responses of common species from alpine swamp meadow to different soil water condition[J]. Ecological Science, 2023, 42(4): 1-8. DOI:10.3969/j.issn.2095-7300.2023.04.001
[12] Kumar S, Suyal D C, Yadav A, et al. Microbial diversity and soil physiochemical characteristic of higher altitude[J]. PLoS One, 2019, 14(3). DOI:10.1371/journal.pone.0213844
[13] Wang X J, Ren Y X, Yu Z Q, et al. Effects of environmental factors on the distribution of microbial communities across soils and lake sediments in the Hoh Xil Nature Reserve of the Qinghai-Tibetan Plateau[J]. Science of the Total Environment, 2022, 838. DOI:10.1016/j.scitotenv.2022.156148
[14] Bastida F, Eldridge D J, García C, et al. Soil microbial diversity-biomass relationships are driven by soil carbon content across global biomes[J]. The ISME Journal, 2021, 15(7): 2081-2091. DOI:10.1038/s41396-021-00906-0
[15] Cai Z Q, Zhang Y H, Yang C, et al. Land-use type strongly shapes community composition, but not always diversity of soil microbes in tropical China[J]. CATENA, 2018, 165: 369-380. DOI:10.1016/j.catena.2018.02.018
[16] Bahram M, Hildebrand F, Forslund S K, et al. Structure and function of the global topsoil microbiome[J]. Nature, 2018, 560(7717): 233-237. DOI:10.1038/s41586-018-0386-6
[17] Hartmann M, Six J. Soil structure and microbiome functions in agroecosystems[J]. Nature Reviews Earth & Environment, 2023, 4(1): 4-18.
[18] 朱永官, 彭静静, 韦中, 等. 土壤微生物组与土壤健康[J]. 中国科学: 生命科学, 2021, 51(1): 1-11.
Zhu Y G, Peng J J, Wei Z, et al. Linking the soil microbiome to soil health[J]. Scientia Sinica Vitae, 2021, 51(1): 1-11.
[19] Wilhelm R C, van Es H M, Buckley D H. Predicting measures of soil health using the microbiome and supervised machine learning[J]. Soil Biology and Biochemistry, 2022, 164. DOI:10.1016/j.soilbio.2021.108472
[20] Zabel F, Delzeit R, Schneider J M, et al. Global impacts of future cropland expansion and intensification on agricultural markets and biodiversity[J]. Nature Communications, 2019, 10(1). DOI:10.1038/s41467-019-10775-z
[21] 陈芬, 余高, 孙约兵, 等. 汞矿区周边农田土壤微生物群落结构特征及其环境驱动因子[J]. 环境科学, 2022, 43(8): 4342-4352.
Chen F, Yu G, Sun Y B, et al. Characteristics of microbial community structure in the surrounding farmlands of a mercury mining area and its environmental driving factors[J]. Environmental Science, 2022, 43(8): 4342-4352.
[22] 吴春晓, 高小峰, 闫本帅, 等. 长期施肥对黄土高原梯田土壤养分特征和微生物资源限制的影响[J]. 环境科学, 2022, 43(1): 521-529.
Wu C X, Gao X F, Yan B S, et al. Effects of long-term fertilization on soil nutrient characteristics and microbial resource restrictions in a terrace on the Loess Plateau[J]. Environmental Science, 2022, 43(1): 521-529.
[23] 孙子欣, 蔡柏岩. 连作对土壤微生物菌群影响及修复研究进展[J]. 作物杂志, 2022(6): 7-13.
Sun Z X, Cai B Y. Research advances on the effects of continuous cropping on soil microbial community and restoration techniques[J]. Crops, 2022(6): 7-13.
[24] Liu C, Li L L, Xie J H, et al. Soil bacterial diversity and potential functions are regulated by long-term conservation tillage and straw mulching[J]. Microorganisms, 2020, 8(6). DOI:10.3390/microorganisms8060836
[25] Hontoria C, García-González I, Quemada M, et al. The cover crop determines the AMF community composition in soil and in roots of maize after a ten-year continuous crop rotation[J]. Science of the Total Environment, 2019, 660: 913-922. DOI:10.1016/j.scitotenv.2019.01.095
[26] Koskey G, Avio L, Turrini A, et al. Durum wheat-lentil relay intercropping enhances soil mycorrhizal activity but does not alter structure of arbuscular mycorrhizal fungal community within roots[J]. Agriculture, Ecosystems & Environment, 2023, 357. DOI:10.1016/j.agee.2023.108696
[27] 王倩, 刘红文, 贾淑霞, 等. 保护性耕作对东北黑土微生物碳循环功能基因的影响[J]. 生态学报, 2023, 43(11): 4760-4771.
Wang Q, Liu H W, Jia S X, et al. Effect of conservation tillage on microbial functional genes related to carbon cycle of black soil[J]. Acta Ecologica Sinica, 2023, 43(11): 4760-4771.
[28] 张蕾, 王强, 杨新月, 等. 黄土丘陵区退耕还林对土壤真菌群落的影响[J]. 环境科学, 2023, 44(3): 1758-1767.
Zhang L, Wang Q, Yang X Y, et al. Effect of the process of returning farmland to forest in the loess hilly area on soil fungal communities[J]. Environmental Science, 2023, 44(3): 1758-1767.
[29] 梁香寒, 张克斌, 乔厦. 半干旱黄土区柠条林土壤水分和养分与群落多样性关系[J]. 生态环境学报, 2019, 28(9): 1748-1756.
Liang X H, Zhang K B, Qiao X. Relationship between soil moisture and nutrients and plant diversity of Caragana microphylla community in semi-arid loess region[J]. Ecology and Environmental Sciences, 2019, 28(9): 1748-1756.
[30] 谭雪莲, 郭天文, 胡新元, 等. 黄土高原旱作区马铃薯连作根际土壤微生物群落变化特征[J]. 作物学报, 2022, 48(3): 682-694.
Tan X L, Guo T W, Hu X Y, et al. Characteristics of microbial community in the rhizosphere soil of continuous potato cropping in arid regions of the Loess Plateau[J]. Acta Agronomica Sinica, 2022, 48(3): 682-694.
[31] 张健, 徐明. 黄土丘陵区不同植被类型土壤细菌群落多样性特征[J]. 生态与农村环境学报, 2022, 38(2): 225-235.
Zhang J, Xu M. Characteristic of soil bacterial community diversity among different vegetation types in the loess hilly region[J]. Journal of Ecology and Rural Environment, 2022, 38(2): 225-235.
[32] Xiao L, Huang Y M, Zhao J F, et al. Effects of planting structure on soil water-stable aggregates, microbial biomass and enzyme activity in a catchment of Loess Plateau terraces, China[J]. Applied Soil Ecology, 2021, 159. DOI:10.1016/j.apsoil.2020.103819
[33] 鲍 士旦. 土壤农化分析[M]. (第三版). 北京: 中国农业出版社, 2000.
[34] Li J B, Shen Z H, Li C N, et al. Stair-step pattern of soil bacterial diversity mainly driven by pH and vegetation types along the elevational gradients of Gongga Mountain, China[J]. Frontiers in Microbiology, 2018, 9. DOI:10.3389/fmicb.2018.00569
[35] Aßhauer K P, Wemheuer B, Daniel R, et al. Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data[J]. Bioinformatics, 2015, 31(17): 2882-2884. DOI:10.1093/bioinformatics/btv287
[36] Qiu Y, Fu B J, Wang J, et al. Spatial variability of soil moisture content and its relation to environmental indices in a semi-arid gully catchment of the Loess Plateau, China[J]. Journal of Arid Environments, 2001, 49(4): 723-750. DOI:10.1006/jare.2001.0828
[37] 周华坤, 李珊, 孙建, 等. 三江源区高寒草甸植物群落与土壤理化性质沿海拔梯度的变化特征[J]. 草地学报, 2023, 31(6): 1735-1743.
Zhou H K, Li S, Sun J, et al. Characteristics of plant community and soil physical and chemical properties in alpine meadow along altitude gradient in the Headwaters Region of Three-River on Tibetan Plateau[J]. Acta Agrestia Sinica, 2023, 31(6): 1735-1743.
[38] 黄斌, 王泉泉, 李定强, 等. 南岭山地土壤有机碳及组分海拔梯度变化特征[J]. 土壤通报, 2022, 53(2): 374-383.
Huang B, Wang Q Q, Li D Q, et al. Variation characteristics of organic carbon and fractions in soils along the altitude gradient in Nanling Mountains[J]. Chinese Journal of Soil Science, 2022, 53(2): 374-383.
[39] Li G X, Xu G R, Shen C C, et al. Contrasting elevational diversity patterns for soil bacteria between two ecosystems divided by the treeline[J]. Science China Life Sciences, 2016, 59(11): 1177-1186. DOI:10.1007/s11427-016-0072-6
[40] Singh D, Lee-Cruz L, Kim W S, et al. Strong elevational trends in soil bacterial community composition on Mt. Halla, South Korea[J]. Soil Biology and Biochemistry, 2014, 68: 140-149. DOI:10.1016/j.soilbio.2013.09.027
[41] Wang J T, Zheng Y M, Hu H W, et al. Soil pH determines the alpha diversity but not beta diversity of soil fungal community along altitude in a typical Tibetan forest ecosystem[J]. Journal of Soils and Sediments, 2015, 15(5): 1224-1232. DOI:10.1007/s11368-015-1070-1
[42] Li H, Yang S, Semenov M V, et al. Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community[J]. Global Change Biology, 2021, 27(12): 2763-2779. DOI:10.1111/gcb.15593
[43] 乔沙沙, 周永娜, 刘晋仙, 等. 关帝山针叶林土壤细菌群落结构特征[J]. 林业科学, 2017, 53(2): 89-99.
Qiao S S, Zhou Y N, Liu J X, et al. Characteristics of soil bacterial community structure in coniferous forests of Guandi Mountains, Shanxi Province[J]. Scientia Silvae Sinicae, 2017, 53(2): 89-99.
[44] 王颖, 宗宁, 何念鹏, 等. 青藏高原高寒草甸不同海拔梯度下土壤微生物群落碳代谢多样性[J]. 生态学报, 2018, 38(16): 5837-5845.
Wang Y, Zong N, He N P, et al. Soil microbial functional diversity patterns and drivers along an elevation gradient on Qinghai-Tibet, China[J]. Acta Ecologica Sinica, 2018, 38(16): 5837-5845.
[45] 吴则焰, 林文雄, 陈志芳, 等. 武夷山国家自然保护区不同植被类型土壤微生物群落特征[J]. 应用生态学报, 2013, 24(8): 2301-2309.
Wu Z Y, Lin W X, Chen Z F, et al. Characteristics of soil microbial community under different vegetation types in Wuyishan National Nature reserve, East China[J]. Chinese Journal of Applied Ecology, 2013, 24(8): 2301-2309.
[46] 韩冬雪, 王宁, 王楠楠, 等. 不同海拔红松林土壤微生物功能多样性[J]. 应用生态学报, 2015, 26(12): 3649-3656.
Han D X, Wang N, Wang N N, et al. Soil microbial functional diversity of different altitude Pinuskoraiensis forests[J]. Chinese Journal of Applied Ecology, 2015, 26(12): 3649-3656.
[47] 李巍, 刘洋, 罗钦, 等. 武夷山常绿阔叶林土壤微生物多样性的季节动态[J]. 热带亚热带植物学报, 2017, 25(2): 115-126.
Li W, Liu Y, Luo Q, et al. Seasonal dynamics in soil microorganisms diversity of evergreen broad-leaved forest in Wuyi Mountains, Southeastern China[J]. Journal of Tropical and Subtropical Botany, 2017, 25(2): 115-126.
[48] Lauber C L, Hamady M, Knight R, et al. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale[J]. Applied and Environmental Microbiology, 2009, 75(15): 5111-5120. DOI:10.1128/AEM.00335-09
[49] 宋贤冲. 广西猫儿山土壤微生物群落时空变化特征研究[D]. 北京: 中国林业科学研究院, 2019.
Song X C. Study on spatio-temporal variation of soil microbial community of Mao'er mountain in Guangxi[D]. Beijing: Chinese Academy of Forestry, 2019.
[50] 林春英, 李希来, 张玉欣, 等. 黄河源区高寒沼泽湿地土壤微生物群落结构对不同退化的响应[J]. 环境科学, 2021, 42(8): 3971-3984.
Lin C Y, Li X L, Zhang Y X, et al. Responses of different degradation stages of alpine wetland on soil microbial community in the Yellow River Source Zone[J]. Environmental Science, 2021, 42(8): 3971-3984.