凋落物是木本植物残体的统称,主要由地上部分的枯枝落叶组成,是森林生态系统中有机质的重要来源
[1]。凋落物分解是一个复杂的过程,在非生物和生物因素驱动下发生物理、化学和生物变化,经过淋溶、粉碎、微生物分解后,凋落物中的营养物质逐步归还到土壤中
[2]。此外,凋落物还能防止地面侵蚀,减少雨滴对土壤团聚体的破坏,在保持生态平衡方面起到重要作用
[3]。微生物群落不仅能通过分解作用促进凋落物降解为自身生长发育繁殖提供营养物质、为植物提供可吸收利用的养分,在维持全球生态系统稳定性上也发挥着十分重要的作用
[4]。随着全球气候的持续变化,凋落物的分解过程及微生物群落结构都会受到影响
[5]。高寒森林作为全球气候变化的敏感地区,被认为是研究自然生态过程的天然实验室
[6]。然而,目前对于微生物群落结构和功能在全球气候变化背景下如何影响凋落物的分解尚不清楚。
环境变化对凋落物分解的影响仍需进一步的研究来量化。在海拔梯度间,温度、湿度和土壤理化性质等因素均存在较大差异
[7],因此不同海拔常被用来模拟全球气候变化所带来的影响
[8]。海拔不仅会改变环境中的物理和化学条件,生物性质也会发生改变
[9]。随着海拔的升高,森林树种趋向单一、凋落物种类变少
[10],土壤微生物群落结构也因此受到影响
[11]。细菌群落易利用高质量、低碳氮比的底物,而真菌群落在低质量、高碳氮比底物的分解中占主导地位
[12]。目前已有对色季拉山不同海拔微生物群落组成的研究发现微生物群落的分布可能主要由植被类型和土壤性质所决定
[13],但不同海拔梯度上的微生物群落会如何影响森林凋落物的分解仍未可知。
色季拉山位于青藏高原东南部,海拔梯度为2 200~5 300 m、植被类型丰富、以高寒森林为主,是自然环境极为独特的区域。森林凋落物分解是色季拉山高寒森林生态系统物质循环和能量流动的关键,由于色季拉山地形复杂且凋落物分解受多种因素影响,目前对该过程的研究较少。因此,本研究拟以西藏色季拉山海拔3 500~4 300 m的高寒森林为研究对象,采用凋落袋掩埋法和高通量测序技术,探究不同海拔下凋落物的降解规律及微生物群落对凋落物分解的响应,阐明凋落物分解过程中微生物群落结构和多样性的变化,揭示影响凋落物降解率和微生物群落结构的土壤理化因素,以期为色季拉山高寒森林不同海拔微生物群落结构变化和凋落物降解规律提供科学依据。
1 材料与方法
1.1 试验设计
研究区域位于西藏自治区林芝地区的色季拉山(29°10′~30°15′ N,93°12′~95°35′ E),属于念青唐古拉山余脉,整体呈西南-东北走向。根据气象资料显示,色季拉山全年平均气温为5 ℃,相对湿度为73%,年降水量为866 mm,降水多集中在4—10月。植被类型丰富,主要有急尖长苞冷杉林(Abies georgei)、康定柳(Salix paraplesia)、高山杜鹃(Rhododendron simsii)和方枝柏(Juniperus saltuaria)等,土壤多为黄棕壤和山地棕壤。
2020年8月进行踏勘选点、凋落物取样,收集代表高山杜鹃当年最新的自然凋落叶,进行清洗和烘干处理。准确称取10.00 g烘干凋落物,装入孔径为60目的20 cm×20 cm尼龙网袋中。样地分别设置在色季拉山海拔3 500、3 700、3 900、4 100、4 300 m处杜鹃林下,每个样地设置1个20 m×20 m的样方,具体信息见
表1。采用五点取样法在每个样方内确定5个样点,于2020年8月在每个样点的0~10 cm表层土壤掩埋3个凋落袋,分别于2020年12月(雪被覆盖期)、2021年4月(雪被融化期)、2021年8月(植被生长期)在每个样方回收1个凋落袋并采集0~10 cm表层土壤,将每个样地的5份土样混合为1个土壤样品,用冷藏箱转移至实验室。土壤样品一部分风干后常温保存,另一部分放入-80 ℃冰箱保存,用于土壤微生物测序。
1.2 试验方法
1.2.1 凋落物降解率及土壤理化性质测定
凋落袋回收后去除杂质,于65 ℃下烘干至质量不变,计算凋落物损失质量/10.00 g初始质量,即为凋落物降解率;利用烘干法测定土壤含水率,pH计测定土壤酸碱度(pH),元素分析仪测定土壤有机碳(soil organic carbon,SOC)、全氮(total nitrogen,TN)并计算土壤碳氮比(carbon nitrogen ratio,C/N),碳酸氢钠提取后用抗坏血酸/钼酸盐试剂法测定土壤有效磷(available phosphorus,AP),乙酸铵萃取-火焰光度法测定土壤速效钾(available potassium,AK),吸管法测定土壤机械组成中黏粒和粉粒的含量(clay and silt content,clay+silt)以及砂粒的含量(sand content,sand),通过测定对土壤进行HCl处理后释放的CO2体积来确定碳酸盐含量(carbonate content,CO32-)。
1.2.2 土壤微生物
利用FastDNA Spin Kit for Soil自旋试剂盒提取土壤DNA。在ABI GeneAmp 9700 PCR热循环仪(ABI,CA,United States)上,用细菌16S rRNA基因的引物对338F(5′-ATCCCTACGGGGGGGGGAGGCAG-3′)和806R(5′-GGATTACHVGGGTWTCTAAT-3′)以及真菌ITS1区域的引物ITS1F(5′-CTGGTCATTTAGGGAAGTAA-3′)和ITS2R(5′-GTGCGTTCTTCATCGATGC-3′)对每个样品的纯化DNA进行PCR扩增。使用AxyPrep DNA凝胶提取试剂盒(Axygen Biosciences,Union City,CA,United States)收集随后的PCR产物,测序工作由上海美吉生物医药科技有限公司完成。然后在Illumina MiSeq PE300平台/NovaSeq PE250平台(Illuminia,San Diego,United States)上进行配对测序,这是根据标准协议完成的。测序后,USEARCH(版本8.0.1623)使用>97%作为共享同一性阈值进行了操作分类单元(operational taxonomic unit,OTU)分配,OTU的代表性序列被分配用于分类分析。
1.3 统计分析
使用Microsoft Excel 2019初步整理数据,R 4.0.2对数据进行分析及作图,微生物共现网络分析参考Zheng等
[14]的方法。
2 结果与分析
2.1 不同海拔土壤微生物群落组成
不同海拔土壤中共有可检测出细菌门为43个,真菌门为17个。根据物种注释结果,不同海拔土壤细菌和真菌群落在门水平相对丰度排名的前10位的物种如
图1所示。可知:在细菌群落中,优势菌门为变形菌门(Proteobacteria)、放线菌门(Actinobacteriota)、酸杆菌门(Acidobacteria)和绿弯菌门(Chloroflexi),其中变形菌门相对丰度最大,在不同海拔均超过25%;细菌群落组成相对稳定,优势菌门的种类不变,而相对丰度较易受到海拔影响。在真菌群落中,优势菌门为担子菌门(Basidiomycota)、子囊菌门(Ascomycota)和被孢霉门(Mortierellomycota);真菌群落组成易受凋落物降解影响,优势菌门的种类和相对丰度的变化较大,在海拔4 300 m处差异尤为显著。
土壤理化性质对微生物群落具有重要影响。分别对细菌和真菌群落的优势菌门与土壤理化性质进行Mantel检验,结果见
表2。可知:在细菌群落中,变形菌门(Proteobacteria)受
C/N显著影响(
P<0.05),放线菌门(Actinobacteriota)受总氮显著影响(
P<0.05);在真菌群落中,担子菌门(Basidiomycota)和子囊菌门(Ascomycota)受土壤含水率显著影响(
P<0.05),被孢霉门(Mortierellomycota)受
C/N显著影响(
P<0.05)。
2.2 不同海拔土壤微生物群落多样性
在凋落物分解过程中,不同海拔土壤细菌和真菌群落的丰富度和多样性指数如
图2所示。可知:除海拔4 300 m外,细菌和真菌群落的Observed OTUs和Chao1指数均在雪被融化期,即凋落物降解8个月时达到最大,而海拔4 300 m处在雪被覆盖期,即凋落物降解4个月后达到最大。在凋落物分解的3个时期,Shannon和Simpson指数随海拔的变化趋势在细菌群落中相同,在真菌群落中也相同;在雪被覆盖期,细菌和真菌群落的Shannon和Simpson指数随海拔变化最为剧烈。如
表3所示,土壤细菌和真菌群落的α-多样性均受到凋落物降解时间、海拔及二者交互作用的显著影响(
P<0.001)。
通过NMDS分析进一步探究土壤微生物群落β-多样性在不同海拔间的差异如
图3所示。可知:在细菌群落中,NMDS2轴可将高海拔4 100、4 300 m处与低海拔3 500、3 700、3 900 m的细菌群落明显区分开,且高海拔处细菌群落结构在凋落物降解的不同时期差异较低海拔处更大;在真菌群落中,具有类似的分布特征。在凋落物分解过程中,细菌和真菌群落在低海拔3 500、3 700、3 900 m和高海拔4 100、4 300 m均具有明显的差异;海拔越高,微生物群落结构的差异越大。
2.3 不同海拔土壤微生物群落共现网络
对5个海拔高度处的土壤微生物群落进行网络分析结果见
图4,可见不同海拔间网络复杂性差异很大。网络参数结果见
表4。可知:海拔3 700 m处的细菌和真菌共生网络具有最大的平均度和图密度,网络的复杂程度最高,模块化程度最低;海拔4 300 m处,网络节点和边的数量最大,模块化程度较高;海拔3 500、3 900、4 100 m处,节点和边的数量、模块化程度、平均度均随海拔升高而增加,网络结构随海拔升高而更加复杂。在5个海拔的共生网络中,真菌群落均占主导地位,以担子菌门(Basidiomycota)和子囊菌门(Ascomycota)为主,且节点之间均以正相关关系为主,说明细菌和真菌群落合作关系大于竞争关系。
2.4 不同海拔凋落物降解率及土壤理化性质
不同海拔凋落物降解率和土壤理化性质结果见图
5和
6。可知:在凋落物分解过程中,不同海拔的凋落物降解率随时间变化而逐渐增大,海拔3 700>4 300>3 500>3 900>4 100 m,其中海拔3 500和3 900 m降解率相似程度较大。除海拔3 700和4 300 m外,凋落物降解率随海拔升高而下降。受雪被覆盖和融化的影响,凋落物的降解速率逐渐下降;经过一年的时间,凋落物在不同海拔的降解率为28.2%~45.2%。在凋落物分解过程中,不同海拔间的土壤理化性质差异较大;经历冻融现象后,各海拔的有效磷含量达到最大,黏粒和粉粒、总氮含量最低,在植被生长期时有较大提升;凋落物分解的不同时期,碳氮比随海拔升高总体呈下降趋势,土壤含水率随海拔升高总体呈上升趋势。
2.5 土壤微生物群落结构和凋落物降解率的影响因素
微生物群落与凋落物降解率及影响因素的结构方程模型如
图7所示。可知:海拔、土壤含水率、
C/
N分别对土壤微生物群落和凋落物降解率产生不同程度的影响:在细菌群落中,海拔不仅对群落结构(
P<0.001)和凋落物降解率(
P<0.05)有显著负的直接影响,还对土壤含水率和
C/
N有显著负的直接影响(
P<0.01),进而间接影响群落结构和凋落物降解率;此外,凋落物降解率还受到
C/
N显著负的直接影响(
P<0.01)。在真菌群落中,海拔也对群落结构有显著负的直接影响(
P<0.05),对土壤含水率和
C/
N有显著负的直接影响(
P<0.01),进而间接影响群落结构和凋落物降解率,凋落物降解率还受到真菌群落显著正的直接影响(
P<0.001)和
C/
N显著负的直接影响(
P<0.001)。
3 讨 论
目前普遍认同微生物群落结构是影响生态系统过程的关键生物因素
[15]。在凋落物分解过程中,细菌和真菌群落的协同作用不仅能够促进分解,微生物的群落结构也会发生改变
[16]。在门水平上,细菌群落的优势种类不随海拔和凋落物降解时间而发生改变,其相对丰度在不同海拔间有所差异,主要受总氮、
C/
N显著影响(
P<0.05);真菌群落的优势种及相对丰度随着凋落物的分解变化较大,主要受土壤含水率、
C/
N显著影响(
P<0.05)。真菌群落较细菌群落对凋落物分解的响应更大,且凋落物的分解过程是由真菌主导的
[17]。在凋落物分解早期,变形菌门(Proteobacteria)和子囊菌门(Ascomycota)是最常见的门
[18],分别在细菌和真菌群落中占主导地位。随着凋落物分解时间的推移,易降解的化合物率先被分解,Bhatnagar等
[19]研究发现,酸杆菌门(Acidobacteria)在木质素的降解中发挥着重要作用,但在本研究中并未发现不同海拔细菌群落优势菌群的显著变化,可能是因为变形菌门、放线菌门和酸杆菌门的生态幅较宽,对环境干扰有一定的抵抗能力
[20]。然而,真菌群落似乎发生了演替,不同海拔处子囊菌门(Ascomycota)相对丰度下降、担子菌门(Basidiomycota)相对丰度逐渐增加,担子菌门比子囊菌门对难降解化合物的分解转化能力更强。此外,微生物群落多样性也受到凋落物分解、海拔及二者的交互影响。除海拔4 300 m外,其余各处的细菌和真菌群落丰富度指数均在雪被融化期达到最大,说明冻融限制了微生物群落对凋落物分解的响应
[21];由于高海拔4 300 m处温度低、雪被融化慢,微生物群落活性的恢复落后于其他海拔。随着凋落物的分解,微生物群落多样性指数在不同海拔间的差异逐渐降低,这与底物基质和微生物群落的代谢有关
[22]。
NMDS分析表明不同海拔的细菌和真菌群落结构存在差异,高海拔4 100、4 300 m处、低海拔3 500、3 700、3 900 m处的微生物群落分别聚集在NMDS2轴的两侧,这是因为海拔会影响太阳辐射和降水等环境因子,进而导致微生物群落结构的差异
[23]。在凋落物分解的不同阶段,微生物群落结构存在明显差异。此外,不同海拔间细菌和真菌群落的共生网络复杂程度也存在差异,海拔3 700 m处网络最为复杂,海拔4 300 m处细菌和真菌群落OTUs联系最为密切、节点和边数量最大、模块化程度最高,海拔3 500、3 900、4 100 m处网络复杂程度随海拔升高而增加。这似乎与凋落物的降解存在一定的联系,凋落物在不同海拔处降解率为3 700 >4 300 >3 500>3 900>4 100 m,低海拔3 700 m处海拔温度、湿度均适宜微生物活动和代谢,因而凋落物降解率最大;高海拔4 300 m处光辐射强度大,光降解作用加速了凋落物分解,促进了微生物间的相互作用;海拔3 500、3 900、4 100 m处的凋落物降解率呈现出随海拔升高而下降的趋势,这与Liu等
[24]的研究结果相同,但与细菌和真菌共生网络的复杂程度相反。然而,微生物群落间的相互作用程度并不仅受海拔和凋落物分解影响而产生差异,也存在一定的共性。不同海拔的细菌和真菌共生网络中,主要类群均为真菌群落的担子菌门(Basidiomycota)和子囊菌门(Ascomycota),且二者在凋落物分解中占主导地位
[22],进而影响凋落物的降解率。
凋落物是高寒森林生态系统中物质与能量的重要来源,其分解是生物地球化学循环的关键过程。经过一年的时间,凋落物有近一半的质量损失,这与华西雨屏区常绿阔叶林的凋落物降解率较为相似
[25]。随着时间的推移,不同海拔的凋落物降解率逐渐增加,但凋落物降解速率逐渐降低。在凋落物分解过程中,易溶化合物如淀粉、氨基酸、多糖最先被快速分解,接着纤维素和木质素成为主要分解成分,部分半纤维素与木质素共价结合,其结构复杂因而分解速度较慢。在高寒森林生态系统中,凋落物分解是一个复杂的过程,除了受微生物因素的影响,也会受到非生物因素影响:1)气候因素,如温度、湿度;2)土壤理化性质,如pH、
C/
N等,且本研究发现凋落物降解率与
C/
N呈显著负相关(
P<0.001)。高寒森林所分布的海拔跨度较大,温度、湿度和土壤理化性质均受海拔梯度的影响,进而影响凋落物的降解率
[26]。本研究发现海拔对
C/
N和土壤含水率产生显著负的直接影响(
P<0.01),对细菌(
P<0.001)和真菌(
P<0.05)群落也分别产生显著负的直接影响,凋落物降解率也受到上述因素直接或间接的影响。在细菌群落中,海拔可以直接影响凋落物降解率(
R2=-0.87,
P<0.05);在真菌群落中,海拔通过影响真菌群落来间接影响凋落物降解率(
R2=0.67,
P<0.001),由此可见真菌群落能够直接影响凋落物降解率,并且在高寒森林凋落物的分解过程中占主导地位
[20]。
4 结 论
本研究为高寒森林生态系统中不同海拔凋落物降解规律和微生物群落的响应提供了重要依据,主要结论如下:
1)在凋落物分解过程中,细菌群落组成相对稳定,而不同海拔的真菌群落则发生了演替;细菌和真菌群落的丰富度和多样性指数不仅受凋落物分解和海拔影响,还受到季节性冻融影响;真菌群落不仅在微生物共现网络中占主导地位,还是凋落物分解的主要贡献者;
2)凋落物降解率在不同海拔间存在较大差异,海拔、C/N和土壤含水率是影响凋落物降解率、细菌和真菌群落的关键因素,真菌群落对凋落物分解的影响更大。基于全球气候变化的背景下,进一步阐明高寒森林等气候敏感地区的土壤微生物群落结构和功能对自然生态过程的响应具有十分重要的意义。
西藏高原森林生态教育部重点实验室开放课题(XZA-JYBSYS-2020-02)
西藏自治区科技创新基地自主研究项目(XZ2022JR0007G)