微生物作为生态系统的核心驱动者,在土壤环境中发挥着不可替代的作用。它们通过分解有机物质、调控养分循环,直接促进植物生长,同时维系着土壤碳功能的稳定性
[1-3]。土壤微生物群落结构(包括细菌、真菌、古菌等的组成、多样性和功能群分布)与土壤有机碳的形态特征之间存在密切的动态互馈关系
[2-3]。一方面,微生物作为有机碳分解与转化的核心驱动者,其群落结构直接调控着不同碳组分的周转速率与路径:例如,真菌及其菌丝网络往往促进大团聚体的形成,有利于颗粒态有机碳的物理保护;而某些细菌类群则倾向于利用易降解的可溶性有机碳或通过分泌胞外酶破解难降解的芳香族化合物,从而影响顽固性有机碳的积累
[2]。另一方面,土壤有机碳的不同形态(如颗粒态有机碳、矿物结合态有机碳、可溶性有机碳等)通过提供差异化的碳源与栖息微环境,反过来塑造了微生物的群落组成与功能
[3-4]。这种相互作用共同决定了土壤碳库的稳定性与周转动态,是理解土壤碳循环响应环境变化的关键机制。
海拔作为重要的环境梯度因子,通过影响温度与降水等气候要素,形成了显著的空间异质特征。海拔梯度变化不仅直接影响土壤理化性质(如pH值、养分有效性)和植被分布格局,还能够改变地表凋落物输入与根系分泌物组成等,进而驱动土壤微生物群落结构与功能的垂直分异
[5-6]。因此,海拔变化形成的天然环境梯度,为解析微生物群落演变规律提供了理想的研究平台。目前,土壤微生物群落多样性对海拔梯度的响应模式尚无统一定论。例如,日本富士山和秦岭东坡的细菌多样性均呈“单峰型”分布,而意大利阿尔卑斯山的真菌多样性则随海拔升高逐渐下降
[5-7]。相比之下,韩国汉拿山的细菌多样性在中海拔最低,呈现“U”型分布
[8]。除多样性分布模式外,微生物群落的海拔分布特征同样备受关注。如Islam等
[9]研究发现,喜马拉雅山濒危植物
Trillium govanianum根际微生物群落结构具有显著海拔差异:变形菌门(Proteobacteria)丰度随海拔升高而增加,放线菌门(Actinobacteriota)则显著下降,酸杆菌门(Acidobacteriota)在中海拔区占优势。林惠瑛等
[10]的研究表明武夷山黄山松林不同海拔细菌群落结构差异显著,其中酸杆菌门和绿弯菌门(Chloroflexi)丰度均随海拔升高递增。此外,土壤微生物群落的海拔分布格局受到多种因素的影响。例如,长白山森林土壤细菌群落不同海拔间的差异性是由于pH变化引起的
[11]。而Li等
[12]在北京东灵山的研究发现了土壤有效磷(Available Phosphorus, AP)和碳氮比(C/N)在影响土壤微生物分布方面的重要作用。然而,因不同研究区域的环境条件和生态系统类型差异显著,导致现有研究结果间也存在较大差异。因此,进一步探究不同海拔梯度下微生物群落的分布格局及其驱动机制,对于揭示微生物多样性变化的潜在规律具有重要意义。
秦岭不同海拔区域作为中国气候分界线和生态安全屏障,其独特的地理格局对区域降水分布和气候调节功能产生了深远的影响。位于秦岭不同海拔火地塘林区,以其超过90%的森林覆盖率和典型垂直带谱特征,成为研究土壤微生物群落海拔分布和有机碳结构组成的天然试验平台。如周煜杰等
[13]探讨了秦岭火地塘林区不同海拔土壤真菌群落的海拔分异规律,揭示土壤有效磷含量、pH和海拔是该地区真菌群落变化的核心驱动因素;王海华
[14]研究发现,秦岭火地塘地区不同海拔油松、华山松根际微生物群落结构的季节变化受土壤养分含量的显著影响。红桦(
Betula albosinensis)作为该区域的优势树种之一,不仅具有显著的经济价值,更是山地生态恢复的先锋树种,在涵养水源、碳汇固存以及生物多样性保护等关键生态过程中,发挥着不可替代的屏障作用。然而,目前针对秦岭不同海拔地区红桦林的研究多聚焦于土壤营养成分以及森林群落结构与稳定性,对红桦林土壤微生物群落沿海拔梯度的分布特征及其驱动因素的系统研究仍较匮乏。分子生物学技术的发展,如高通量测序和分子生态网络分析,为深入解析森林土壤细菌群落多样性及功能提供了有力工具。本研究选取秦岭不同海拔典型植被红桦林为研究对象,运用高通量测序技术,探究低海拔(1 500 m, LA)、中海拔(2 000 m, MA)、高海拔(2 500 m, HA)下土壤细菌群落结构组成和有机碳分子结构,并阐明其对海拔的响应机制,为揭示秦岭生态系统微生物地理分布格局的形成机制提供科学依据。
1 材料与方法
1.1 研究区概况
研究区位于陕西省宁陕县秦岭山区林场(33°18′—34°26′N,106°04′—110°40′E),海拔1 420~2 574 m,年平均温度8~10 ℃,年降水量900~1 200 mm,年日照数1 100~1 300 h,无霜期170 d,属北亚热带山地气候。土壤类型主要有山地棕壤、暗棕壤和山地草甸土,腐殖质层较厚,土壤肥力高,且结构良好。林区地形复杂,坡度一般在20°~50°。
1.2 样品采集
2025年6月于秦岭不同海拔林区采集土壤样本(
表1),采样地海拔高度分别设置为低海拔(1 500 m, LA)、中海拔(2 000 m, MA)、高海拔(2 500 m, HA),每个海拔梯度内设置4个10 m×10 m的重复样地,样地间隔≥15 m。所选样地的地形、坡度、坡向基本一致。在每个样地内,采用“S”型取样法选取5个采样位点,使用土壤取样器采集0—10 cm土层土壤,去除可见根系与碎石后,将同一样地内的土样充分混匀,过2 mm尼龙筛以剔除落叶及有机残体。过筛后土样分为3份:一份置于阴凉通风处自然晾干,风干后取其部分研磨并通过0.15 mm(100目)筛用于全量养分测定,剩余部分保持2 mm粒径用于pH与速效养分测定;另一份立即放入4 ℃冰箱保存,用于土壤酶活性测定;最后一份立即放入-80 ℃冰箱中保存,用于微生物分析。
1.3 土壤理化性质和有机碳的测定
土壤pH按水土比为2.5∶1使用pH计测定;土壤有机质(Soil Organic Matter, SOM)、全氮(Total Nitrogen, TN)含量分别采用重铬酸钾外加热法、凯氏定氮法;全磷(Total Phosphorus, TP)采用硫酸—高氯酸消煮—钼锑抗分光光度法测定。
易氧化有机碳(Easily Oxidizable Organic Carbon, EOC),轻组分有机碳(Light Fraction Organic Carbon, LOC)、重组分有机碳(Heavy Fraction Organic Carbon, HOC)、颗粒有机碳(Particulate Organic Carbon, POC)、矿物态结合碳(Mineral-Associated Organic Carbon, MOC)和可溶性有机碳(Dissolved Organic Carbon, DOC)参照文献
[15]的方法测定。
1.4 土壤细菌群落分析
1.4.1 土壤DNA提取和测序
本研究采用试剂盒法提取土壤总DNA,并通过琼脂糖凝胶电泳对其浓度及完整性进行定量分析。随后,基于蓝色荧光定量系统对PCR扩增产物进行检测,并将各样本扩增产物按等摩尔比混合。混合产物经磁珠纯化后,对目标片段进行验证。进一步借助MiSeq测序平台完成文库构建与双端测序,回收产物经洗脱纯化后通过电泳进行质检,然后对原始测序数据进行质控过滤,获得高质量序列。最后采用FLASH等工具对原始序列进行优化处理。
1.4.2 生物信息分析及数据处理
首先对每个样品的测序数据进行序列拆分、质量剪切与拼接处理。随后,基于QIIME平台对高通量序列进行进一步分析,具体包括:对有效序列进行聚类操作,将相似度≥97%的序列归并为同一分类单元(ASVs),并从中选取代表序列进行物种注释
[4]。
测序深度指数Coverage:该指数可用于评估测序数据的代表性,数值越高,表明测序的完整性和准确性越好,所反映的微生物群落真实情况也更为可靠。
细菌丰度和多样性指数选择Chao1指数、Ace指数、Simpson指数、Shannon指数。
1.5 数据分析
运用Excel 2019整理基础数据及绘图。采用SPSS 20.0统计软件进行试验数据处理,用方差分析法(Analysis of Variance, ANOVA)进行差异显著性分析,以及邓肯法进行组间多重分析。利用QIIME2.0计算微生物多样性指数,FastUniFrac进行NMDS排序,采用AMOS 21.0构建、验证并解析结构方程模型(Structural Equation Model, SEM)。
2 结果与分析
2.1 秦岭不同海拔土壤理化性质
由
表2可知,秦岭不同海拔土壤pH值变化范围为6.03~7.89,平均值为6.67,随着海拔梯度的增加,土壤pH值呈逐渐增加趋势;土壤有机碳SOC变化范围为9.11~12.11 g/kg,平均值为11.07 g/kg,随着海拔梯度的增加,SOC呈逐渐降低趋势;NH
4+-N和NO
3--N的平均值分别为5.27 mg/kg和5.00 mg/kg,最大值都出现在低海拔,随着海拔梯度的增加,NH
4+-N和NO
3--N呈逐渐降低趋势;对于TN,与SOC变化趋势相一致,随着海拔梯度的增加呈逐渐降低趋势,平均值为0.57 g/kg, TP的变化趋势并不明显,不同海拔之间差异均不显著(
p>0.05)。
2.2 秦岭不同海拔土壤有机碳形态结构
表3为秦岭不同海拔土壤有机碳形态结构,由表可知,海拔梯度对土壤碳形态产生了显著影响。土壤轻组有机碳LOC变化范围在0.38~0.59 g/kg之间,平均值为0.47 g/kg,随着海拔梯度的增加,LOC呈逐渐降低趋势;土壤重组有机碳HOC变化范围在2.12~3.67 g/kg之间,平均值为2.85 g/kg,随着海拔梯度的增加,HOC呈逐渐降低趋势;易氧化有机碳EOC和可溶性有机碳DOC平均值分别为249 mg/kg和19 mg/kg,随着海拔梯度的增加,EOC和DOC呈逐渐降低趋势。
2.3 秦岭不同海拔土壤细菌多样性分析
表4为秦岭不同海拔90个土壤样品中的优化序列、ASV数量及多样性指数,通过对秦岭不同海拔采集的90份土壤样品进行高通量测序,共获得35 864条高质量序列,平均长度为8 966 bp,在97%相似度水平下进行分析。各样品文库的Coverage指数均接近99%,表明测序深度充分,能够较好反映土壤细菌群落的实际组成。由
表4结果可知,不同海拔土壤细菌群落的丰度与多样性变化趋势一致:随着海拔升高,土壤细菌的Chao1指数、Ace指数、Shannon指数及ASV数量均呈逐渐下降趋势,其中土壤细菌盖度和Simpson指数变化幅度较小。
2.4 秦岭不同海拔土壤细菌门水平分类和纲水平分类
高通量测序结果显示,研究区土壤样品中共检出细菌32门、69纲、178目。经质控过滤后,从25 678条有效序列中获得19 871条高质量序列,其中超过99%的序列可鉴定至细菌门水平。如
图1所示,土壤细菌群落的优势菌门为变形菌门(相对丰度32%)、放线菌门(相对丰度25%)和酸杆菌门(相对丰度17%),其次为绿弯菌门、芽单胞菌门和厚壁菌门,这6类菌门约占总序列95%;拟杆菌门和硝化螺旋菌门也属于主要菌门。此外,还发现了27个稀有细菌门类,相对丰度低于1%。
在细菌纲分类水平上,α-变形菌纲(相对丰度28%)、放线菌纲(相对丰度22%)和酸杆菌纲(相对丰度14%)为优势类群,其次为β-变形菌纲、嗜热油菌纲和芽单胞菌纲,共占总序列约95%;此外也检出杆菌纲和δ-变形菌纲等,但相对丰度较低(
图2)。
2.5 秦岭不同海拔土壤微生物群落结构分布
非度量多维尺度分析(Non-Metric Multidimen-sional Scaling, NMDS)结果表明(
图3),细菌群落结构的胁迫系数(stress)为0.121。该系数反映实际生态距离与低维空间映射距离的偏差程度,本研究所得数值均低于0.2的可接受阈值,证实分析结果具有可靠性。基于ANOSIM分析,
R值(范围0~1)表征组内相似性与组间差异性,当
R值趋近1且
p<0.05时,表明组间差异具有统计学意义。本研究中,细菌群落的
R值达到0.812(
p=0.001),充分证实不同海拔对土壤细菌群落的分布影响较大。
2.6 土壤有机碳与细菌群落结构的相关性
土壤理化性质和细菌群落结构的相关性采用Mantel检验方法(Mantel tests)分析。相关性分析表明,Proteobacteria, Acidobacteria, Actinobacteria, Firm-icutes, Nitrospirae的相对丰富度与土壤pH值呈显著的负相关;Proteobacteria, Acidobacteria, Chlor-oflexi, Gemmatimonadetes, Firmicutes, Nitrospirae与碳形态呈显著的正相关;Proteobacteria, Actinob-acteria, Acidobacteria, Gemmatimonadetes, Firm-icutes, Nitrospirae与NH
4+-N呈显著的正相关;Proteo-bacteria, Acidobacteria, Chloroflexi, Gemmatim-onadetes, Bacteroidetes,与NO
3--N呈显著的正相关;Proteobacteria, Chloroflexi, Gemmatimonadetes, Firmicutes, Bacteroidetes与TN呈显著的正相关;不同细菌门相对丰富度与TP没有显著的相关性(
图4和
图5)。综上,与秦岭不同海拔土壤细菌群落相关性较大的土壤化学指标是pH, SOC和碳形态指标。
从
表5可知,土壤EOC和DOC与细菌门的优势菌Proteobacteria, Actinobacteria, Acidobacteria相对丰度呈显著的正相关,土壤EOC和DOC与细菌门的优势菌Chloroflexi, Gemmatimonadetes, Firmicutes, Bacteroidetes相对丰度呈显著的负相关。土壤HOC和LOC与细菌门的优势菌相对丰度均没有显著的相关性。
从
表6可知,土壤EOC和DOC与细菌纲的优势菌α-Proteobacteria, Actinobacteria, Acidobacteria、β-Proteobacteria, Thermoleophilia相对丰度呈显著的正相关,土壤EOC和DOC与细菌纲的优势菌Gemmatimonadetes, Bacilli、δ-Proteobacteria相对丰度呈显著的负相关。土壤HOC和LOC与细菌纲的优势菌相对丰度均没有显著的相关性。
2.7 秦岭不同海拔环境因素对土壤细菌群落结构的影响
采用结构方程模型探讨秦岭不同海拔环境因素对土壤细菌群落结构的影响,
图6中卡方自由度比值为0.51,小于3证明模型适配度较好。对于海拔梯度,对细菌群落多样性有显著负向影响,路径系数为-0.74(
p<0.01),对细菌群落多样性有正向影响,路径系数为0.59。土壤碳形态对细菌群落结构有显著负向影响,路径系数为-0.89(
p<0.01),对细菌群落多样性有正向影响,路径系数为0.61。结果表明,土壤碳组分对细菌群落多样性有负向影响,对细菌群落结构有正向影响。
3 讨 论
秦岭不同海拔土壤pH随海拔升高而显著增加,这与赵明松等
[16]在河南西部研究得出的海拔与土壤pH呈显著正相关结果相同。这种格局可能主要源于凋落物分解所驱动的土壤酸化强度的变化,低海拔区域气温较高,秦岭植被茂盛,凋落物输入量大且土壤有机质含量丰富,显著促进了微生物活动。正如Singh等
[8]的研究指出,凋落物分解过程会产生单宁、有机酸等酸性物质,在低海拔活跃的微生物代谢环境下,大量酸性中间产物持续积累,导致土壤pH值降低,酸化加剧。相反,高海拔地区气温低,凋落物输入少且分解慢,微生物活动相对较弱,酸性物质生成较少,土壤酸化较弱,pH相对低海拔较高。因此,低海拔强酸化与高海拔弱酸化的差异,最终使土壤pH随海拔升高而升高。此外,土壤TP含量随海拔升高而升高,这可能与高海拔红桦林土壤中更高的磷酸酶活性有关。磷酸酶在土壤代谢过程中发挥着重要作用,能够将土壤中的有机磷化合物水解为无机磷,从而增加土壤中可利用磷的含量
[17]。值得注意的是,本研究发现土壤NH
4+-N和NO
3--N含量相对较高,这可能源于土壤较高的脲酶活性和固氮功能,脲酶活性的增强可促进有机氮矿化,而固氮作用加强则有利于大气氮的生物固定,这共同导致了土壤氮素的相对富集。土壤EOC(易氧化碳)、LOC(轻组分有机碳)、HOC(重组分有机碳)的增加表明,有机质的加入不仅提高了土壤活性有机碳含量,还促进了较为稳定的碳组分累积,这有助于提高土壤碳库的稳定性。海拔梯度对DOC(可溶性有机碳)的影响较大,表明海拔越低,输入的碳源更容易被微生物利用。随着海拔的逐渐增加,土壤水分含量增加,土壤的通气性和保水率也上升,然而温度降低,抑制了土壤微生物的生理活性,减缓了土壤中的物质转化。
微生物群落多样性是维系生态系统物质循环的核心动力,其丰富的物种组成与功能差异,能显著促进碳、氮、磷等关键营养盐的生物地球化学循环过程
[3, 18]。海拔梯度改变了土壤环境参数(如pH值、有机质含量等),从而显著改变土壤微生物的群落结构特征和生物多样性水平
[3]。秦岭不同海拔土壤细菌多样性指数呈现降低趋势,这与前人的研究结果一致。低海拔地区土壤较高的NH
4+-N和NO
3--N含量和有机碳含量为细菌生长繁殖提供了优越的营养条件,从而提升了该区域土壤细菌的丰度、物种丰富度及群落均匀度。在微生物群落组成方面,变形菌门作为优势细菌菌群,这些优势菌门普遍具有分解和转化土壤有机质的功能,而林区土壤中丰富的有机质为其提供了适宜的生长环境。已有研究证实,该菌门通过参与氮磷循环的关键生化过程,显著促进速效养分的转化与释放,同时在有机物质降解中扮演核心角色
[19-20]。与之功能互补的拟杆菌门则以其突出的有机物分解能力著称,特别是在大分子有机物的矿化过程中发挥不可替代的作用
[21-22]。本研究结果显示,随着海拔的增加,变形菌门相对丰度显著提升,变形菌门作为土壤优势菌群,凭借其卓越的代谢可塑性(包括有机质降解、固氮作用及植物生长促进等多元功能),在驱动土壤养分循环和维持根际微生态平衡中发挥着核心作用
[22]。变形菌门细菌通过其特有的氮磷代谢途径(包括固氮、硝化、反硝化及有机磷矿化等关键过程),显著促进土壤养分的生物有效性转化,同时在复杂有机质的分解与矿化过程中扮演着不可替代的生态功能角色
[23-24]。本研究中低海拔土壤变形菌门丰度最高,说明土壤物理结构的优化显著改善了通气条件和水分保持能力,同时加速了有机物质的分解转化过程,这种微环境改良为变形菌门的增殖创造了有利条件,最终导致其群落相对丰度的显著提升。值得注意的是,不同海拔土壤微生物群落在属水平上表现出显著差异,这与β多样性分析结果一致,进一步证实了土壤微生物群落在海拔上展现出较强的空间异质性。
土壤作为土壤微生物生存和繁衍的生境,其理化特性通过直接调控微生物的生长繁殖、代谢活性和种间竞争关系,成为驱动土壤微生物群落结构动态演替的关键性环境决定因子。这些理化参数的变化会显著影响微生物的生态位分化,进而导致群落组成在时间和空间尺度上发生规律性的演替变化
[25-26]。植物根系分泌物通过调控土壤理化性质(如pH值、养分有效性等)和微环境条件,直接或间接地塑造了根际微生物群落的组成特征和多样性格局
[27]。本研究结果显示,根系分泌物有助于提升土壤微生物多样性,优化群落结构,进而促进养分循环与植物生长。冗余分析表明,秦岭山地土壤微生物多样性与养分含量随海拔变化表现出显著协同变化规律,这与多数山地生态研究结论相一致。其相互作用机制主要体现为海拔差异导致的土壤有机碳、氮等养分的空间分异,直接制约着微生物群落的代谢活性与组成结构;同时,微生物通过增强酶活性等途径,又反向驱动着有机质的转化与养分释放过程。海拔梯度导致的植物群落与土壤理化性质分异,是影响微生物群落空间格局的重要因素。土壤pH值是调控微生物群落结构的关键环境因子,本研究中,土壤pH与微生物群落组成呈显著负相关。此外,本研究发现土壤全磷与微生物群落结构无显著关联。这可能与研究区土壤类型及磷素形态有关,土壤中的磷主要来源于母岩风化,高海拔低温限制了风化与释放,同时植物吸收导致磷持续消耗、补充不足,使得土壤全磷普遍偏低。这种低含量且难利用的磷形态,削弱了磷与微生物群落间的耦合关系。相关性分析表明,环境因子可影响土壤细菌群落组成与多样性,土壤碳组分对细菌群落多样性有负向影响,对细菌群落结构有正向影响。马媛等
[28]在天山北坡雪岭云杉林土壤微生物群落的研究中,发现氮素、pH、年均降水量和年均温度是影响细菌群落的关键环境因子。万盼等
[29]在甘肃小陇山油松(
Pinus tabulaeformis)林和柴胡(
Bupleurum chinense)林土壤研究中发现影响土壤细菌群落结构的关键环境因子包括pH、有机碳和土壤湿度。秦岭不同海拔土壤细菌群落,有机碳结构组成和pH主要影响其分布。pH较低的低海拔地区主要富集了酸杆菌门如酸杆菌纲、酸杆菌目,酸杆菌门在低海拔的富集可能和其嗜酸性有关
[28]。这暗示了秦岭火地塘地区红桦土壤细菌与有机碳结构之间的密切关系,对植物适应高海拔环境具有重要意义。
4 结 论
本文研究了秦岭不同海拔土壤细菌群落特征与有机碳形态结构,结果显示:秦岭不同海拔土壤pH随海拔升高而显著增加,这种格局可能主要源于凋落物分解所驱动的土壤酸化强度的变化。土壤TP含量随海拔升高而升高,这可能与高海拔红桦林土壤中更高的磷酸酶活性有关。土壤NH4+-N和NO3--N含量相对较高,这可能源于土壤较高的脲酶活性和固氮功能,脲酶活性的增强可促进有机氮矿化,而固氮作用加强则有利于大气氮的生物固定,这共同导致了土壤氮素的相对富集。土壤EOC(易氧化碳)、LOC(轻组分有机碳)、HOC(重组分有机碳)在低海拔的增加表明,有机质的加入不仅提高了土壤活性有机碳含量,还促进了较为稳定的碳组分的累积,有助于提高土壤碳库的稳定性。海拔梯度对DOC(可溶性有机碳)的影响较大,表明海拔越低,输入的碳源更容易被微生物利用。随着海拔升高,土壤细菌的Chao1,Ace, Shannon指数以及ASV数量均呈显著递减趋势,而细菌盖度与Simpson指数则无明显变化。在秦岭不同海拔土壤中,变形菌门(Proteobacteria)、放线菌门(Actinobacteria)和酸杆菌门(Acidobacteria)为优势菌门;在菌纲水平上,放线杆菌纲(Actinobacteria)、α-变形菌纲(α- Proteobacteria)和酸杆菌纲(Acidobacteria)占主导地位。相关性分析表明,环境因子可影响土壤细菌群落组成与多样性,土壤碳组分对细菌群落多样性有负向影响,对细菌群落结构有正向影响。秦岭不同海拔土壤细菌群落,有机碳结构组成和pH主要影响其分布。pH较低的低海拔地区主要富集了酸杆菌门如酸杆菌纲、酸杆菌目,酸杆菌门在低海拔的富集可能和其嗜酸性有关。
国家重点研发计划“秦巴山区黑木耳等食用菌良种及先进培植技术研究与应用示范”(2021YFD1600401)