青藏高原东北部高寒草甸土壤细菌群落沿海拔梯度的空间分布及其与环境因子的关系

《Biology》：Spatial Distribution of Soil Bacterial Communities Along an Altitudinal Gradient in Alpine Meadows of the Northeastern Qinghai–Tibet Plateau and Their Relationship with Environmental Factors Wenfang Chen, Huichun Xie, Shuang Ji, Yue Zhang, Xunxun Qiu, Zhiqiang Dong and Jiaxiang Xu

　　

土壤微生物群落在驱动高寒草甸养分循环中扮演着关键角色，然而，其沿海拔梯度的分布格局及其对环境变化的响应在很大程度上仍未得到充分探索。为此，研究人员在青藏高原东北部沿海拔梯度（3300-4500 米）采集了土壤样本，以分析细菌群落组成和多样性及其与土壤理化性质和酶活性的关联。结果表明，不同海拔高度下细菌群落组成和多样性存在显著变异。在所有采样点，Actinomycetota、Pseudomonadota和 Acidobacteriota是优势菌门。以香农（Shannon）指数衡量的群落多样性总体上随海拔升高而增加，在 4500 米处达到峰值，在 3300 米处最低。皮尔逊（Pearson）相关分析和冗余分析（RDA）表明，土壤细菌群落结构与土壤养分因子和酶活性均显著相关。在这些变量中，全钾（TK）、有效磷（AP）、过氧化氢酶（CAT）和脲酶（URE）与细菌群落分异呈强相关。此外，PERMANOVA（基于距离的多元方差置换分析）结果显示，海拔是驱动群落变异的主要因子，在统计学显著水平上解释了群落组成的相当一部分变异。总体而言，本研究揭示了高寒草甸土壤中细菌群落沿海拔梯度的分布格局及其环境关联，为理解高寒生态系统中微生物群落对环境变化的响应提供了基础数据。

研究背景、目的与意义

青藏高原作为“世界屋脊”和“第三极”，其生态系统对区域乃至全球气候具有重要的调节作用。高寒草甸是该地区主要的草地生态系统，在维持生物多样性、涵养水源等方面发挥着关键的生态功能。由于其对环境变化高度敏感，已成为理解高原生态系统动态的研究重点。土壤微生物，特别是土壤细菌，是草地生态系统物质循环和能量流动的关键驱动者，在土壤有机质分解、养分转化和生物地球化学循环中扮演着不可或缺的角色。它们对环境变化高度敏感，其群落结构和多样性常被用作环境变化的生物指标。因此，研究土壤细菌群落沿海拔梯度的变化模式，对于理解高寒草甸生态系统对环境条件的响应至关重要。然而，已有关于土壤细菌群落沿海拔梯度变化的研究结果并不一致，这可能与研究区域、植被类型、海拔范围和环境背景的差异有关。为了在一致的区域内系统地阐明高寒草甸土壤细菌群落沿海拔梯度的变化模式并识别其潜在的驱动因素，Wenfang Chen, Huichun Xie, Shuang Ji, Yue Zhang, Xunxun Qiu, Zhiqiang Dong 和 Jiaxiang Xu 等研究人员在青藏高原东北部的典型高寒草甸地区，沿3300-4500米的海拔梯度设立了样地，通过采集土壤样本，分析了土壤细菌群落的组成、多样性及其与环境因子的关系。这项研究的目的是揭示细菌群落组成和多样性沿海拔梯度的变化，并识别塑造这些群落变异的关键环境驱动因素，为理解高寒草甸土壤细菌的空间格局及其对环境因子的响应提供科学依据，并为高原草地生态系统的生态保护和可持续管理提供见解。本研究发表在期刊《Biology》上。

关键研究方法概述

本研究在青藏高原东北部高寒草甸区，沿海拔梯度（从3300米到4500米，以300米为间隔）设置了五个海拔梯度（3300、3600、3900、4200、4500米），共计采集了60个复合土壤样品进行分析。研究方法主要涵盖以下几个方面：（1）土壤样品采集与理化及酶活性分析：采集0-20厘米表层土壤，测定包括全氮（TN）、全磷（TP）、全钾（TK）、碱解氮（AN）、有效磷（AP）、有效钾（AK）等土壤理化性质，以及脲酶（URE）、转化酶（INV）、过氧化氢酶（CAT）和碱性磷酸酶（AKP）的活性。（2）微生物DNA提取与高通量测序：使用E.Z.N.A.^?Soil DNA试剂盒提取土壤总DNA，对16S rRNA基因的V3-V4高变区进行PCR扩增，在Illumina测序平台进行高通量测序。利用UPARSE软件在97%相似度水平下聚类为操作分类单元（OTU），并在SILVA数据库进行物种注释。（3）生物信息学与统计分析：计算α多样性指数（如Shannon指数、Chao1指数、ACE指数和Coverage指数），并通过主坐标分析（PCoA）、线性判别分析效应量（LEfSe）等方法分析群落结构差异。使用PICRUSt工具对细菌群落功能进行预测注释。通过皮尔逊相关分析、冗余分析（RDA）等统计方法探究细菌群落与环境因子（理化性质和酶活性）之间的关系。

研究结果

3.1. 土壤理化性质沿海拔梯度的变化

土壤全养分（TN、TP、TK、AN）含量在不同海拔间无显著差异。而土壤速效养分表现出显著的海拔梯度差异，其中有效磷（AP）含量先升高后降低，在4500米处达到峰值；有效钾（AK）含量则呈现单峰模式，在3600米处显著高于其他海拔。

3.2. 不同海拔下的土壤酶活性

土壤酶活性对海拔的响应各异。参与氮循环的脲酶（URE）活性、与磷循环相关的碱性磷酸酶（AKP）活性和与碳循环相关的转化酶（INV）活性均在3900米处出现峰值，在4200米处显著下降，随后在4500米处再次上升。过氧化氢酶（CAT）活性沿海拔梯度呈波动趋势，在3600米处最高，在3900米处显著降低，在4200米处再次升高，最后在4500米处下降。

3.3. 土壤细菌群落组成随海拔的变化

在所有海拔点共检测到50个细菌门，其中14个为优势门。Actinomycetota在所有位点均为最丰富的门，其次是 Pseudomonadota、Acidobacteriota和 Chloroflexota。这些优势菌门的相对丰度沿海拔梯度发生显著变化。LEfSe分析揭示了不同海拔梯度下的特征性细菌类群，例如，高海拔位点（4500米）富含Elusimicrobiota、Patescibacteria等类群，表明高海拔土壤对微生物群落有更强的环境过滤作用。

3.4. 不同海拔的土壤细菌多样性分析

3.4.1. α多样性

Shannon指数随海拔升高呈整体增加趋势，在4500米处最高，表明高海拔地区的细菌多样性更高。ACE指数和Chao1指数在3900米处显著高于4200米，表明该海拔物种丰富度较高。Coverage指数在所有位点均接近1，表明测序深度足够。

3.4.2. β多样性

Venn图分析显示，不同海拔土壤细菌群落共享3658个OTU，占OTU总数的12.60%。每个海拔点都包含一组独特的OTU，其中3900米处的特有OTU数量最多。基于Bray-Curtis距离的主坐标分析（PCoA）也表明，不同海拔组的土壤细菌群落能明显分开，相同海拔的样本聚集紧密。

3.5. 不同海拔土壤细菌群落的功能多样性

基于PICRUSt的功能预测显示，在所有海拔点，代谢（Metabolism）均是土壤细菌群落中最主要的功能类别，占预测功能的78%以上。在KEGG（京都基因与基因组百科全书）二级功能层级上，“全局与概览图谱（Global and overview maps）”、碳水化合物代谢、氨基酸代谢和能量代谢是最丰富的亚类。COG（同源蛋白簇）功能注释同样显示，氨基酸转运与代谢、能量产生与转换、翻译/核糖体结构与生物发生、细胞壁/膜/被膜生物发生是主要的保守功能类别。这些主要功能类别的相对丰度在不同海拔间变化不大，表明高寒草甸土壤细菌群落的功能结构具有稳定性。

3.6. 影响土壤细菌群落组成和多样性的关键因素

Pearson相关分析表明，多个土壤养分、酶活性与主要细菌类群的相对丰度存在显著关联。冗余分析（RDA）进一步阐明了土壤细菌群落结构与环境因子之间的关系。RDA的前两轴累计解释了群落总变异的18.18%。排序结果显示，土壤养分因子（TN、TP、TK、AN、AP、AK）对观测到的群落结构变异贡献显著。不同海拔的样本在RDA排序空间中有部分分离，例如3300米的样本与CAT、URE和INV活性关联更强，而4500米的样本与N、P、K等养分因子的关联更强。

讨论与结论

讨论要点总结

研究结果表明，青藏高原东北部高寒草甸的土壤细菌群落在不同海拔存在明显差异。优势菌门 Actinomycetota、Pseudomonadota和 Acidobacteriota的丰度沿海拔梯度发生分异，表明细菌类群在高寒草甸生态系统中经历了沿环境梯度的生态位分化。LEfSe分析揭示了特定海拔的特征类群，表明海拔梯度通过环境过滤作用对微生物群落进行了选择性塑造。土壤细菌群落的α多样性（Shannon指数）随海拔升高而增加，这与部分类似高寒生态系统的研究结果一致。β多样性分析进一步显示，不同海拔的土壤细菌群落结构存在显著分化。功能预测分析揭示了一个“群落结构分化但功能保守”的模式：尽管细菌群落的物种组成沿海拔发生变化，但其核心代谢功能（如碳、氮、能量代谢）在不同海拔间保持高度一致，这表明功能冗余可能在维持高寒草甸土壤生态系统功能稳定性中发挥关键作用。细菌群落沿海拔梯度的分布格局是多种环境因子相互作用塑造的结果。除了海拔作为一个综合环境梯度外，土壤养分有效性（特别是AP和TK）和酶活性（如CAT和URE）的变化是驱动细菌群落分异的关键因素。冗余分析（RDA）结果也支持了这一观点，表明土壤养分状况是影响高寒草甸细菌群落空间分异的关键环境背景因子。研究也承认了其局限性，包括未全面考虑植被类型、放牧干扰等其他因素的交互作用，采样设计不平衡，以及环境因子覆盖不完整等。

研究结论（翻译自原文）

本研究揭示了青藏高原东北部高寒草甸土壤细菌群落沿海拔梯度（3300-4500 米）的一个独特生态模式，其特征是多样性增加、结构分异和功能稳定。尽管土壤理化性质在梯度上没有表现出显著的整体变化，但细菌群落表现出明显的海拔分异：群落多样性（香农指数）随海拔升高显著增加，表明高海拔环境拥有更丰富的细菌物种库。在门水平上，Pseudomonadota、Actinomycetota和 Acidobacteriota被确定为主导类群，它们的相对丰度沿海拔梯度发生了显著演替。值得注意的是，基于PICRUSt的功能预测揭示，尽管存在结构分异，但核心代谢功能（例如碳、氮和能量代谢）在不同海拔间保持高度一致，这表明“结构可塑，功能稳定”可能是高寒生态系统中微生物的一种潜在适应策略。冗余分析进一步将海拔、全钾（TK）、有效磷（AP）、过氧化氢酶（CAT）和脲酶（URE）确定为驱动群落结构分异的关键环境因子。

总而言之，本研究首次系统阐明了青藏高原东北部高寒草甸土壤细菌群落的海拔分布模式，揭示了细菌群落结构对海拔的敏感性及其功能谱系的保守性。这些发现为理解极端环境中微生物群落的生态适应策略提供了新的视角。此外，它们丰富了高寒生态系统微生物生物地理学的理论理解，并为预测土壤微生物群落对气候变化的响应建立了基础数据集。