在全球变暖持续加速的背景下,高海拔生态系统因其低温、强紫外辐射、季节性冻融和养分贫瘠等多重胁迫而表现出极端的脆弱性和高度的气候敏感性,被视为全球变化的“预警区”
[1,2]。随着全球变暖加剧,高山森林的物种组成、群落结构和生态功能正经历显著变化
[3],急剧变暖正在重塑高山森林的物种组成、群落结构和生态功能
[4]。外生菌根(Ectomycorrhizal, ECM)真菌作为植物-土壤互作的关键纽带,其多样性及生态适应性对维持高山森林的稳定性至关重要
[5]。
ECM真菌在森林生态系统中的功能虽然已被广泛研究,但多集中于低海拔或温带森林,而对高海拔地区,尤其是青藏高原边缘山地森林的ECM真菌群落特征及驱动因素仍缺乏系统性认知。海拔作为生态环境的一个重要复合因子,其变化可带来温度、降水、土壤理化性质、微生物活性及宿主植物生理特性的变化,显著影响ECM真菌的群落组成和多样性
[6-8]。相关研究表明,ECM真菌的丰富度多随海拔升高呈单峰模式或递减趋势
[9],但在阿尔卑斯山地区,某些耐寒真菌类群在高海拔环境下则表现出更强的竞争优势
[10],青藏高原作为“世界屋脊”,多样的生态系统和气候变化必然会引起ECM真菌多样的环境响应机制。急尖长苞冷杉作为青藏高原色季拉山优势建群树种,其分布范围狭窄,受气候变化和其他因素的影响,种群数量呈下降趋势,其共生ECM真菌的多样性可能进一步受到宿主特异性、生境片段化及环境过滤效应的调控
[8],但目前相关研究仍较为匮乏。
急尖长苞冷杉(
Abies georgei var
. smithii)是青藏高原东南缘色季拉山地区亚高山暗针叶林的建群种和优势树种,对维持区域水源涵养、碳固定及生物多样性具有不可替代的作用
[11]。高海拔环境下的低温、强紫外线辐射、季节性冻融及养分匮乏等胁迫条件严重限制了植物的生长和养分吸收
[12],但ECM真菌可通过与宿主植物形成共生关系,显著增强植物对氮、磷等矿物元素的吸收能力,并提高其抗逆性
[5,13-15],增强其对高海拔环境的适应能力。因此,研究急尖长苞冷杉ECM真菌多样性沿海拔梯度的分布格局及其环境驱动因子,对深入理解高山森林生态系统对未来气候变化的响应策略具有重要意义。
本研究以藏东南色季拉山急尖长苞冷杉林为研究对象,对不同海拔梯度ECM真菌进行采样与测定,结合DNA提取技术、PCR技术和Sanger测序方法,分析色季拉山急尖长苞冷杉ECM真菌的物种组成和生物多样性,同时分析不同海拔ECM群落的变化规律及差异,探讨关键环境因子对ECM真菌群落构建的影响,进一步完善高海拔冷杉林ECM真菌群落结构与多样性,以期为该区域的生物多样性保护、预测未来气候变化下植物-微生物共生关系的演变趋势及退化森林生态系统的恢复提供理论支撑。
1 材料与方法
1.1 研究区概况
色季拉山(29°30′ - 30°10′ N,94°20′ - 94°50′E)位于西藏自治区林芝市境内,属于念青唐古拉山脉东南段,是雅鲁藏布江下游的重要支流——尼洋河与帕隆藏布江的分水岭
[16]。该区域地处青藏高原东南缘,属于典型的高山峡谷地貌,海拔范围约2 800 ~ 5 300 m,地势起伏显著,坡度陡峭,形成明显的垂直气候带和植被带。受印度洋暖湿气流和高原寒流的共同影响,气候表现为温凉湿润的高原季风气候,年均气温-2 ℃ ~ 2 ℃,最热月(7月)均温9 ℃ ~ 12℃,最冷月(1月)均温-15 ℃ ~ -5℃。年降水量800 ~ 1 200 mm,主要集中在5 - 9月,占全年降水的80%以上。由于海拔梯度变化大,气候垂直分异明显,低海拔地区( < 3 000 m)气候较温和,而高海拔地区( > 4 000 m)常年低温,冬季积雪覆盖期长,季节性冻融作用显著。受气候、植被和地形的影响,土壤类型随海拔变化呈现明显的垂直地带性分布,低海拔地区( < 3 500 m):以山地棕壤为主,腐殖质层较厚,有机质含量较高(5% ~ 10%),pH值5.0 ~ 6.5,土壤肥力较好;中海拔地区(3 500 ~ 4 200 m):过渡为山地暗棕壤和亚高山草甸土,有机质积累更显著(10% ~ 15%),但随海拔升高,土壤风化减弱,砾石含量增加。高海拔地区( > 4 200 m):以高山草甸土和寒冻土为主,土层浅薄,有机质分解缓慢,土壤呈微酸性至中性(pH 6.0 ~ 7.0)
[17-19]。
1.2 样品采集
2024年7 - 9月,在色季拉山东坡,以200 m为一个海拔梯度,设置4个50 m×200 m的样带,4个样带海拔分别为3 700 m、3 900 m、4 100 m、4 300 m,每个样带内按左、中、右选取3个急尖长苞冷杉外生菌根采样点,即3个重复,每个采样点内随机选取5株标准木作为采样木,每株采样木之间间距大于10 m,共计60株,记录每株独立木的树高、胸径、冠幅、地理坐标等信息。以每株采样木为中心分别在其东南和西北两个方向对侧根进行外生菌根样品采集,每株采样木采集4 ~ 6根侧根(所有侧根数不少于150个根尖)。采集后装入塑料袋中并编号,在塑料袋中加入适量腐殖质,喷洒些许水分,防止菌根真菌样品失活,随后将当天采集的急尖长苞冷杉外生菌根放入实验室冰箱中进行冷藏,温度保持在4 ℃,以保持其活性。在每株采样木根系样品取样完毕后将根际周围土壤装入自封袋中作为土壤样品,将每个采样点的土壤样品混合均匀,剔除砾石杂物等,带回实验室阴干,过100目土壤筛后保存,用于测定土壤理化性质(样带信息见
表1)。
1.3 菌根样品测序
1.3.1 DNA提取
用于DNA提取的根尖全部取自色季拉山急尖长苞冷杉树种。采用EZup柱式真菌基因组DNA抽提试剂盒(上海生工生物股份有限公司)进行提取,操作步骤严格按照试剂盒操作说明进行。
1.3.2 PCR扩增
PCR处理以2x Taq PCR Mix为酶,ITS1-ITS4为引物进行扩增,各体系用量为模板DNA 10 μL、2x Taq PCR Mix 25 μL、ITS1 2 μL、ITS4 2 μL、Sterilized ddH
2O 11 μL,总体积为50 μL。反应程序见
表2。PCR扩增产物的定量测定、均一化、文库制备、测序及质量控制等工作均由上海生工生物股份有限公司完成。通过质量筛查获得有限的序列,对优化后的序列进行97%聚类分析。运用GenBank的序列局部相似性查询系统(BLAST),在GenBank、EMBL、DDBJ及PDB 4个DNA序列数据库中搜索与所测DNA序列相似的序列,并进行分类鉴定。
1.4 土壤理化性质测定
土壤理化性质测定:土壤pH值用电位法测定;土壤含水率(WC)采用烘干法计算;土壤电导率(EC)使用电导率仪测定;土壤样品进行风干后,经高氯酸-硫酸消煮后使用凯氏定氮法测定土壤全氮(TN);利用重铬酸钾氧化测定土壤有机碳(TOC);采用氢氧化钠熔融-钼锑抗比色法测定土壤全磷(TP)。
1.5 统计分析
土壤理化数据利用Excel软件进行预处理,采用SPSS 25.0软件进行统计分析。利用单因素方差分析方法分析不同海拔ECM真菌多样性指数,显著性检验采用Duncan′s法;采用主成分分析法(PCA)分析不同海拔ECM真菌组成差异性;使用R Studio软件进行多样性及丰富度指数计算;利用Canoco 5.0软件进行真菌多样性指数、优势菌属与环境因子的冗余分析(RDA);利用Origin 2021软件进行相关性分析绘图;并计算OUT相对丰度:
2 结果与分析
2.1 土壤理化性质
4个海拔样带土壤理化性质见
表3,不同海拔之间的电导率EC、pH、TN、TP、TOC和WC有显著性差异(
P < 0.05)。3 700 m的EC、pH和TN含量最高,且随着海拔的提升,EC、pH和TN这3个土壤理化性质会下降,TP和TOC在4 100 m最高,TP随海拔的提升而增大,TOC随海拔的提升呈现先增后减的规律,WC在3 900 m最大,与海拔变化未呈现明显的规律。
2.2 色季拉山不同海拔急尖长苞冷杉外生菌根真菌侵染率
图1显示,各海拔急尖长苞冷杉的ECM侵染率差异不显著(
P > 0.05),以4 100 m海拔带的侵染最高(41.29%),3 700 m海拔带的侵染率次之(41.17%),而后依次是3 900 m海拔带(41.12%)、4 300 m海拔带(39.73%)。
2.3 急尖长苞冷杉ECM真菌测序结果统计
对不同海拔急尖长苞冷杉ECM真菌
OTUs经BLAST进行DNA序列比对,经鉴定共有49个
OTUs与目前已知的种类相似度较高,隶属于担子菌门和子囊菌门(见
表4)。其中,担子菌门和子囊菌门的
OTU比例分别为92.59%和7.41%。
通过计算
OTUs的相对丰度和组成,不同海拔的急尖长苞冷杉林分的ECM真菌群落组成结构存在差异(
图2)。3 700 m急尖长苞冷杉林分中相对丰度较高的菌属是红菇属(
Russula)和丝膜菌属(
Cortinarius),其相对丰度为32.35%和23.53%,其次是乳菇属(
Lactarius),为14.71%。在3 900 m的急尖长苞冷杉林分中,相对丰度较高的菌属是红菇属,占25.71%;乳菇属占22.86%;丝膜菌属占20.67%。在4 100 m的急尖长苞冷杉林分中,优势菌属是红菇属,占27.50%;乳菇属和丝膜菌属分别占20.00%。在4 300 m的急尖长苞冷杉林分中,优势菌属是红菇属和乳菇属,分别占37.50%和29.17%;丝膜菌属占16.67%。可见,红菇属、乳菇属和丝膜菌属是急尖长苞冷杉林分ECM真菌群落中的优势菌属,但不同海拔其相对丰度比例不同,红菇属在不同海拔急尖长苞冷杉林分中相对丰度均最高。
2.4 不同海拔急尖长苞冷杉ECM真菌群落组成
根据测序鉴定结果表明,利用RStudio软件进行去重和排除,所有序列被分为49个
OTUs。不同海拔急尖长苞冷杉ECM共有和特有
OTU数为:Pop1(3 700 m)共有急尖长苞冷杉ECM的
OTU数为35个,特有2个;Pop2(3 900 m)共有36个,特有4个;Pop3(4 100 m)最多,共有42个,特有2个;Pop4(4 300 m)最少,共有24个,特有1个。各海拔共有急尖长苞冷杉ECM的OTU数为17个(见
图3)。
2.5 不同海拔ECM真菌多样性差异
通过对4个海拔样带的ECM真菌群落进行Simpson多样性指数、Shannon-Wiener指数及Pielou均匀度指数分析,4 100 m海拔的急尖长苞冷杉ECM真菌群落最为丰富,Simpson多样性指数最小,为0.53,与其他海拔呈现显著性差异(P < 0.05);Shannon-Wiener指数最大,为0.98,显著高于其他海拔。当海拔低于4 100 m时,Simpson多样性指数和Shannon-Wiener指数随海拔的提升呈现出上升的趋势,Pielou均匀度指数在不同海拔之间无显著差异。整体来看,仅4 100 m海拔样带的真菌多样性最为突出,其他海拔样带的真菌多样性变化差异不大。
2.6 海拔梯度及土壤理化性质对ECM真菌多样性的影响
根据不同海拔中的真菌信息进行主成分分析(PCA),结果如
图5所示,第一轴解释率(PCA1)为57.99%,第二轴解释率(PCA2)为32.70%,总解释率可达到90.69%,说明不同海拔样带急尖长苞冷杉ECM真菌样本相对分散,同一海拔样带ECM真菌样本相对聚集,这表明不同海拔急尖长苞冷杉ECM真菌群落组成存在明显差异。
对海拔及土壤理化性质与不同海拔急尖长苞冷杉ECM真菌优势菌属进行RDA分析(
图6A)表明,RDA1轴的解释量为64.02%,RDA2轴的解释量为29.90%,总解释量为93.92%,可以反映相应实际情况。海拔(ALT)、pH、TP以及含水率WC与红菇属和乳菇属呈负相关,TN、TOC、EC以及土壤温度(AST)与红菇属和乳菇属呈正相关;丝膜菌属与TN、TOC、EC以及AST呈负相关,与海拔ALT、pH、TP以及含水率呈正相关。对ECM真菌多样性与土壤理化性质进行RDA分析(
图6B)发现,TOC、TP、ALT以及pH与ECM真菌Shannon-Wiener指数、Simpson指数和Pielou指数呈正相关,与其余因子呈负相关。其中,ECM真菌Shannon-Wiener指数与TOC和TP正相关关系较强,与TN呈现较强的负相关;ECM真菌Simpson指数与TOC和TP呈现较强的正相关,与TN呈现较强的负相关;ECM真菌Pielou指数与ALT和pH呈现较强的正相关,与EC、WC以及AST呈现较强的负相关。
对海拔和土壤理化性质进行Pearson相关性分析(
图7)表明,土壤TP和WC会随着海拔的升高而增大,AST、EC和TN会随着海拔的升高而减小,土壤TP与EC、TN以及AST呈极显著负相关(
P < 0.01),与pH呈显著正相关(
P < 0.05),TOC与EC和TN呈极显著正相关(
P < 0.01),EC与TN和AST呈极显著正相关(
P < 0.01),全氮TN与AST呈极显著正相关(
P < 0.01)。
3 讨论
3.1 不同海拔ECM真菌群落结构
本研究中不同海拔的ECM真菌侵染率差异不大,可能是因为急尖长苞冷杉本身适应性较强,能在不同海拔广泛分布,其根系系统可为ECM真菌提供稳定的生态位。即使海拔变化,宿主植物的生物量或根系活性未显著改变,ECM真菌的侵染率则会保持稳定。通过Sanger测序共获得49个急尖长苞冷杉ECM真菌
OTUs,说明急尖长苞冷杉ECM真菌资源较为丰富。49个真菌OTUs绝大部分是担子菌门真菌,小部分是子囊菌门真菌。担子菌门真菌是急尖长苞冷杉ECM真菌的绝对优势真菌,占所有菌根ECM真菌的92.59%,这与冷杉属(
Abies)树种ECM真菌多样性的多数研究结果基本一致
[20,21],体现了担子菌门真菌在森林生态系统中的绝对重要地位。但在相关研究中,也有子囊菌门的相对丰度高于担子菌门的情况
[22],这可能与生境差异及ECM真菌的共性及专一性有关
[23],子囊菌门里的许多腐生性真菌可以通过分解土壤中难以降解的有机质和营养物质,提高土壤养分利用效率
[24],这可能导致子囊菌门真菌的进化及其物种多样性可能高于担子菌门。担子菌门真菌大多为腐生或者寄生菌,而子囊菌门真菌绝大多数为腐生菌
[25]。从属层面来看,优势菌属分别是红菇属、乳菇属和丝膜菌属,且都属于担子菌门。有研究表明,这3个属的外生菌根真菌不仅可以与大部分裸子植物如冷杉属、落叶松属、云杉属等广泛形成ECM共生系统,还能与部分被子植物形成互利共生系统
[26-28],如栎、松、桦、椴、杨、柳等
[29,30]。进一步说明担子菌门外生菌根真菌不仅是色季拉山急尖长苞冷杉种群的优势菌种,同时还能与其他多种被子植物形成共生关系,具有重要应用开发前景。
优势菌属在不同海拔急尖长苞冷杉林分中的分布和相对丰度有所差别,这是由林分生长状态决定的,在相同生态位上,ECM真菌会因为寄主的需求改变而产生直接替代或直接失去优势地位,不同海拔林分中ECM真菌群落动态的明显变化可能是优势菌群相互竞争相互制约的结果
[31]。尽管随海拔变化,急尖长苞冷杉林分的ECM真菌优势菌属分布及其相对丰度发生动态改变,但其整体多样性变化差异并不显著,这一特征也体现了真菌多样性在环境变化中的弹性及稳定性
[5]。
真菌多样性在环境变化中的稳定性并非偶然,而是ECM真菌群落对环境变化的功能性适应
[32]。在本研究中,4 100 m海拔的急尖长苞冷杉林ECM真菌的Shannon-Wiener指数最高,为0.98,显著高于其他海拔(
P < 0.05),该海拔具备相对稳定的水热条件和较高的根系分泌物输入,有利于维持多样化的ECM功能群;随着海拔的提升,ECM真菌的Shannon-Wiener指数呈现出先增后减的趋势,3 700 m海拔的急尖长苞冷杉林ECM真菌的Shannon-Wiener指数仅为0.52,显著低于其他海拔(
P < 0.05),这可能是因为海拔变化通过调控温度和水分可利用性间接影响土壤微生物活性,进而改变ECM真菌的多样性模式,中海拔(4 100 m)样带的ECM真菌功能群最为丰富,这可能与该海拔区间稳定的水热条件和较高的根系分泌物输入有关。而在高海拔(4 300 m)区域,低温限制土壤微生物的分解作用,导致养分矿化速率降低,从而可能迫使ECM真菌与宿主植物形成更紧密的共生关系以获取养分
[33],表现出功能冗余与生态位替代机制,从而维持ECM真菌群落结构的相对稳定,这一机制与鲁莹等
[34]在樟子松人工林研究中得出的结论一致:土壤真菌群落在不同林龄间α多样性与β多样性均无显著差异,群落功能趋于稳定,生态功能更趋有利于植物生长。
ECM真菌群落的稳定性还体现在其对环境变化扰动的缓冲能力
[32],高海拔季节性冻融导致的土壤扰动可能进一步减少高海拔ECM真菌的物种共存
[35,36]。但在本研究中,红菇属、乳菇属和丝膜菌属广谱性菌根真菌在不同海拔中仍保持优势地位,表明这些类群具备广泛的宿主适应性与环境可塑性,可在资源受限或扰动频繁的环境中维持其生态功能。这与“微生物群落构建的生态驱动框架”相呼应
[37],即微生物群落通过生态位分化与扩散限制机制,在环境变化中实现结构与功能的动态平衡。
3.2 ECM真菌多样性与土壤因子的关系
土壤作为影响植被生长的重要因素之一,其理化性质影响植物根系真菌群落组成
[38]。本研究结果显示,急尖长苞冷杉ECM真菌多样性指数与土壤因子呈较强的相关性。中海拔区域的ECM真菌多样性更高,这与该海拔范围内较高的土壤有机质含量和适宜的pH密切相关。已有研究证实,土壤pH是驱动ECM真菌群落的关键因素
[39],土壤pH可以影响ECM真菌的生长及其菌丝的代谢,同时,土壤pH可以通过改变土壤养分的可用性来间接改变ECM真菌群落种类
[40]。随着海拔的提升,pH呈降低趋势,不同海拔之间的pH值存在显著差异,虽然4个海拔之间的平均pH值变化相差不大,但是即使pH值变化很小,环境土壤中的ECM真菌群落也会发生明显变化。pH值还可以通过影响真菌胞外酶的活性来调控其养分获取能力
[41],高海拔样带的ECM真菌多样性较低,可能也与该区域土壤温度低、冻融作用频繁及养分可利用性较差有关
[36]。
基于RDA分析,TOC和TP是影响ECM真菌群落结构的关键因子。在高海拔样带,较高的磷含量与红菇属和乳菇属等腐生型ECM真菌的富集显著相关,这些类群可能更依赖易分解有机质。而在低海拔区域,较低的磷有效性更有利于丝膜菌属等耐胁迫真菌的生存,这些类群通常具有更强的养分捕获能力。这一结果支持了“环境过滤”假说
[42],即土壤养分状况通过选择适应性强的ECM真菌类群来塑造群落结构
[43]。
4 结论
(1)优势群落组成:色季拉山急尖长苞冷杉的ECM真菌群落中,占据优势地位的菌属主要隶属于担子菌门和子囊菌门,这些优势菌属包括红菇属、乳菇属和丝膜菌属。
(2)海拔梯度下的群落差异:结合ECM真菌物种丰富度和系统发育进化距离两个维度进行分析,发现不同海拔的急尖长苞冷杉ECM真菌群落在组成结构上差异变化不大。
(3)最适海拔与丰富度:色季拉山海拔4 100 m处的土壤理化性质最有利于急尖长苞冷杉ECM真菌群落的生长和发育,这直接导致了该海拔ECM真菌群落的物种丰富度显著高于其他研究海拔样带的ECM真菌群落的物种丰富度。
(4)关键影响环境因子:土壤pH值、有机碳含量和全磷含量是影响色季拉山急尖长苞冷杉ECM真菌群落结构和分布的主要环境驱动因子。
西藏高原森林生态教育部重点实验室开放课题(XZA-JYBSYS-2022-03)
西藏高原森林生态教育部重点实验室研究生创新计划项目(XZAJYBSYS-202505)