大量煤矸石长期露天堆存不仅侵占土地资源,还可能因自燃、淋溶等方式释放重金属和酸性物质污染周边水、土及大气环境,严重威胁矿区及周边生态安全
[1],实现煤矸石的无害化处理和资源化利用已成为矿区生态修复的关键任务。目前填埋复垦是矿区进行生态重建的主要手段,即直接在煤矸石表层覆盖一定厚度土壤,重构植物生长基质,以恢复矿区生态功能
[2]。覆土厚度作为复垦工程中的重要参数,直接影响土壤水分运移
[3]、土壤团聚体及养分
[4]、土壤微生物活性
[5]、土壤微生物量碳的积累
[6]、作物生物量
[7]及对水分的利用
[8],同时也关系到复垦的经济成本以及资源投入
[9]。因此,科学确定覆土厚度对协调生态修复效果与经济可行性具有重要意义。
土壤微生物多样性作为评价土壤质量与生态功能的重要指标,其变化与生境条件密切相关
[10]。在植物根区,因受根系分泌物影响,根际与非根际土壤微生物群落往往存在显著差异
[11]。目前对根际与非根际的研究大多集中于外源性物料施入对细菌群落的影响,如有研究表明生物炭施用对玉米非根际土壤细菌群落多样性的增加效果显著,而对根际土壤细菌多样性的影响不大
[12];土壤有机碳累积矿化量与玉米根际和非根际土壤细菌微多样性显著相关
[13]。相比之下,土壤真菌虽然在土壤微生物中数量占比少,却在土壤有机物质分解
[14]及土壤碳氮循环中发挥着重要作用,是土壤生态系统功能的重要驱动者。然而,目前关于覆土厚度如何影响煤矸石复垦区土壤真菌群落及其生态功能的研究仍较缺乏,尤其是不同覆土处理下真菌代谢功能与土壤养分的耦合机制尚不明确。
基于此,本研究提出如下假设:覆土厚度通过改变煤矸石填埋区复垦土壤的理化性质,进而影响土壤真菌群落的多样性、组成结构及其代谢功能,并存在一个最优厚度,能有效地促进根际微生态功能优化。为此,以山西省屯兰矿区煤矸石填埋复垦地为对象,分析不同覆土厚度下玉米根际和非根际土壤真菌群落结构、多样性及代谢功能变化,并结合土壤养分指标,揭示影响真菌群落构建的关键环境因子,以期为煤矸石复垦地土壤重构提供理论依据,促进矿区生态系统的可持续恢复。
1 材料与方法
1.1 研究区概况
本试验中矿区复垦地位于山西省古交市南梁上村(37°52'59"N,112°6'32"E),属温带大陆性气候,年均气温9.5 ℃,年均降水量460 mm。该矿曾经的排矸方式为倾倒于附近天然沟壑中,该区域煤矸石中重金属含量(
表1)均低于《GB15618-2018 土壤环境质量农用地土壤污染风险管控标准》,因此可利用附近山地黄土进行覆土复垦。该复垦地覆盖的黄土呈碱性,pH值为8.43,有机质含量低,仅为5.51 g/kg,氮磷钾全量分别为0.45、0.66和16.09 g/kg,有效磷、速效钾及碱解氮含量分别为17.35、92.32和9.68 mg/kg。本试验各覆土厚度(40、80和120 cm)下的试验样地均长5 m,宽4 m。本研究玉米种植试验播种与施肥分别采用人工点播与人工撒施方式,株距30 cm、行距55 cm,有机肥(7500 kg/hm²)配施无机肥(600 kg/hm²)。
1.2 样品采集
样品于玉米拔节期采集,在各覆土处理区选取3株长势一致的玉米,挖出植株根系后,采用抖根法轻抖去除附着在根系的大块土,用无菌软毛刷将附着在根系的0~2 mm土壤轻轻刷下后立即装入灭菌袋并密封,作为玉米根际土壤样品;采用对角线布点法,使用土钻获取株间0~20 cm土层土壤样品,放入灭菌袋并密封,作为非根际土壤样品。密封后的根际与非根际土壤样品均迅速放入4 ℃保温箱及时带回室内进行测定分析,所有的根际与非根际土壤样品均分为两部分,一部分用于真菌多样性分析,一部分用于土壤化学性质测定。未能及时做真菌试验的土样带回实验室后放在-80 ℃冰箱冷冻保存。用于测定化学性质的土壤样品需经自然风干,研磨并过筛后进行测定。每处理均设置3次重复。
1.3 实验方法与数据处理
土壤基础化学性质测定参照鲍士旦
[15]方法,真菌高通量测序试验包括真菌DNA提取、PCR扩增(ITS1区)、建库测序、质量控制、序列单元ASV(Amplicon Sequence Variant)测定等环节均在上海派森诺生物公司进行。
利用SPSS 26.0软件对不同处理下的土壤真菌多样性进行方差分析与LSD事后多重检验(P<0.05),利用Origin 2021软件进行绘图,使用R语言及其MetagenomeSeq包分析不同处理下有显著差异的真菌、利用Galaxy平台对不同处理下的真菌群落进行LEfSe(LAD Effect Size)分析,根据LAD默认值(LAD=2)确定有显著差异的标志性物种及数量;基于物种丰度矩阵(欧式距离)对不同处理下真菌群落差异与功能进行主成分分析(PCA,Principal component analysis);使用Canoco 5软件对真菌群落与土壤养分进行RDA(Redundancy Analysis)冗余分析。
2 结果与分析
2.1 玉米根际与非根际土壤养分
从
表2中可见,不同覆土厚度下玉米根际与非根际土壤pH值及全氮、有机质、有效磷和速效钾含量差异显著。对于根际土壤,80 cm覆土厚度下玉米根际土壤有机质、有效磷、速效钾含量明显高于40和120 cm,碱解氮与全氮含量与120 cm之间差异不大,但是显著高于40 cm。对于非根际土壤,40 cm覆土厚度下土壤中全氮、有机质、碱解氮、有效磷、速效钾含量均显著低于80和120 cm,而80 cm和120 cm覆土厚度之间全氮、有机质和速效钾含量差异不显著(
P>0.05)。
2.2 玉米根际与非根际土壤真菌Alpha多样性
从
图1a可见,不同覆土厚度下玉米根际与非根际土壤中独有的真菌ASV数量为83~470个,且根际数量均高于非根际。根际土壤中的ASV数量在80 cm覆土厚度下最多,为522,非根际土壤在40 cm的覆土厚度下数量最少,为135。玉米根际与非根际土壤中真菌ASV数量顺序分别为R80>R40>R120和S80>S120>S40。
从
图1b可以看出,对于玉米根际,80 cm厚度下Chao 1指数显著高于40和120 cm(
P<0.05),而Shannon指数在80与40 cm厚度间无显著差异(
P>0.05),但显著高于120 cm厚度(
P<0.05)。对于非根际,80 cm厚度下Chao1指数显著高于40 cm厚度(
P<0.05),但与120 cm厚度无显著差异(
P>0.05);而Shannon指数在不同覆土厚度下无显著差异(
P>0.05)。玉米根际土壤真菌群落在40和80 cm覆土厚度下Chao1指数相对于非根际增加115%和95%,Shannon多样性增加49%和53%,而120 cm厚度下根际与非根际间多样性指数差异不大(
P>0.05)。总体来看,80 cm的覆土厚度较40、120 cm更利于提高复垦土壤真菌多样性。
2.3 玉米根际与非根际土壤真菌群落组成
对不同覆土厚度下玉米根际与非根际真菌群落差异进行PCA主成分分析,从
图2可见,PC1与PC2分别解释48.8%和21.7%,总解释率为70.5%。对非根际土壤真菌群落,PC1将40 cm厚度与80、120 cm厚度处理分隔开,80和120 cm厚度处理间差异不大,均分布在第一象限,与40 cm厚度处理差异较大;根际土壤真菌群落80 cm和40、120 cm间差异较大,80 cm厚度处理位于第四象限,40 cm厚度处理位于二三象限,120 cm厚度处理位于第一二象限。由此可见,不同覆土厚度可以改变玉米根际与非根际真菌群落,尤其80与40 cm厚度间真菌群落组成差异显著。
从
图3a可以看出,玉米根际与非根际土壤中子囊菌门(Ascomycota,73.63%~95.84%)为该复垦地的优势真菌门,其相对丰度在同一覆土厚度下均表现为根际低于非根际,且无论根际还是非根际的相对丰度在80 cm覆土厚度下最低,分别为73.63%和83.87%;而被孢霉门(Mortierellomycota,0.88%~6.51%)的相对丰度在80 cm厚度下最高,且同一覆土厚度下玉米根际中被孢霉门的相对丰度高于非根际。从
图3b可发现,不同覆土厚度下玉米根际与非根际真菌中丰度较高的属有无茎真菌属(
Acaulium,1.6%~49.68%),瓶毛壳属(
Lophotrichus,1.62%~38.76%),被毛枝葡萄孢属(
Botryotrichum,0.41%~26.86%),赤霉属(
Gibberella,0.26%~21.98%),假裸囊菌属(
Pseudogymnoascus,1.24%~12.25%)等。在玉米根际土壤中,80 cm覆土厚度明显降低了瓶毛壳属的相对丰度,而提高了赤霉属的相对丰度;在非根际土壤中,80和120 cm覆土厚度相较于40 cm显著提高了无茎真菌属的丰度,降低了瓶毛壳属和赤霉属的丰度;各覆土厚度下玉米根际假裸囊菌属的相对丰度均大于非根际,而无茎真菌属却相反。总体来看,不同覆土厚度下玉米根际与非根际土壤中真菌门与属结构差异较大。
进一步分析造成不同覆土厚度下玉米根际与非根际土壤中真菌群落组成差异的原因(
图4),发现不同覆土厚度下玉米根际与非根际有显著组间差异的分类单元有31个,其中R80有24个,占总体的77%,有2个门,3个纲,3个目,5个科,11个属;R40有4个,分别为Bionectriaceae科、Helotiaceae科、Helotiales目和
Hydropisphaera属;S40有2个为
Gymnostellatospora和
Kernia,R120有1个为
Poaceascoma,而S80和S120均为0,没有表现显著差异的物种,说明覆土厚度对玉米根际真菌的影响较大,80 cm厚度下的根际土壤真菌群落的响应最强,可以使某些真菌的相对丰度显著增强。MetagenomeSeq分析表明各覆土厚度下根际相对于非根际上调的ASV数量在40、80和120 cm覆土厚度下分别有9、32和25个,这些ASV均以粪壳菌纲为主,在40、80和120 cm的厚度下其数量分别为5、14和8个。总体来看,80 cm覆土厚度下玉米根际与非根际土壤真菌群落差异最显著,且能显著增加粪壳菌纲丰度。
2.4 不同覆土厚度下根际与非根际土壤真菌代谢功能
不同覆土厚度下玉米根际与非根际土壤真菌群落一级功能单元表现为生物合成、降解、利用、同化和前体代谢物和能量的产生。二级功能主要为氨基酸的生物合成,碳水化合物生物合成,辅因子、假体基团、电子载体和维生素生物合成,核苷和核苷酸生物合成、脂肪酸和脂质生物合成、电子转移与呼吸(
图5)。不同样本共有84种代谢功能类群,对其进行因子分析后发现特征值大于1的有7个主成分,第1、2主成分总解释率为75.55%,大部分的代谢通路位于第一和第四象限,35条代谢通路对第1主成分的贡献较大,其中生物合成通路占比最大;对主成分2贡献率较大的为前体代谢物和能量的产生通路,如氮化合物代谢通路(D13)和辅酶Q-7(B21)等通路。从
表3可见,不同覆土厚度下根际真菌代谢功能均高于非根际,无论根际还是非根际,80 cm厚度下真菌代谢功能得分值最大,为0.71,S40得分最小,为-0.39,顺序为R80>R40>R120>S80>S120>S40,说明80 cm覆土厚度下土壤真菌代谢功能最强。
2.5 玉米根际与非根际土壤真菌群落与土壤养分
不同处理下的真菌群落与土壤基础养分RDA分析结果表明(
表4和
图6),玉米根际与非根际土壤真菌群落结构对覆土厚度的响应显著,其中贡献率较高的土壤有机质含量(36.8%,
P<0.01)、有效磷含量(27.3%,
P<0.05)和全氮含量(16.9%,
P<0.05)是主要驱动因子;根际与非根际间真菌群落在80和120 cm处理下差异显著;非根际土壤在40 cm与其它2个处理间存在显著差异,而80与120 cm处理间差异不显著;根际土壤则表现为80 cm与其他处理间差异显著。总体来看,无论根际还是非根际,80与40 cm厚度处理间差异较大。
3 讨论
3.1 覆土厚度与玉米根际与非根际土壤真菌多样性
根际是微生物、土壤与植物根系进行能量养分交换的场所,真菌是微生物的重要组成部分,研究根际真菌多样性对于复垦土壤质量的提升与植物生长有重要意义。对于根际与非根际真菌的研究报道结果不一
[16-17]。本研究表明各覆土厚度下玉米根际土壤真菌多样性均高于非根际,且80 cm的覆土厚度下二者间的差异最显著,与孙美美等
[18]和杨泽良等
[19]的研究结果相似。这是因为植物通过光合作用合成的同化产物会以有机与无机分泌物的形式向根际区域释放
[20],为根际微生物提供生长所需的物质与能量,作物根系活动也可能利于根际土壤中形成一定的特异性微生物群落,使得某些特定类群增殖和丰度提高
[21],导致作物根际与非根际土壤间的微生物多样性产生明显的差异。本研究也表明无论根际还是非根际,80 cm覆土厚度下土壤真菌ASV独有的数量与总数、物种多样性指数Shannon、丰富度指数chao1均高于40和120 cm,且玉米根际与非根际间差异性物种最多,说明该覆土厚度非常有利于根际真菌群落丰富度与物种多样性增加
[13],更利于真菌群落多样性恢复和生态环境可持续发展
[22]。
3.2 覆土厚度与玉米根际与非根际土壤真菌群落结构
有研究表明子囊菌门和担子菌门是矿区复垦地土壤真菌中的优势类群
[23]。本研究发现子囊菌门是该复垦区各覆土厚度下玉米根际与非根际土壤中的优势真菌门,这与子囊菌门可与植物形成共生关系
[24]、生存能力强的特点有关
[10]。该菌门下的无茎真菌属、瓶毛壳属和假裸囊菌属为相对丰度较高的真菌属。本研究发现80 cm覆土厚度下被孢霉门相对丰度最高,原因可能与该覆土厚度下土壤有机质含量最高有关,因为该真菌门对有机质丰富的腐败物适应能力强。同时发现该覆土厚度下根际与非根际土壤真菌群落结构差异最显著,根际相对于非根际显著增加了子囊菌门下粪壳菌纲的相对丰度,原因可能是作物根系分泌物促进该类真菌在根际富集并生长发育。粪壳菌纲对土壤中纤维素的降解功能较强
[25],由于土壤微生物群落多样性与土壤有机碳累积矿化量相关
[13],说明该覆土厚度利于复垦土壤中碳的转化与累积。我们也发现80和120 cm覆土厚度相较于40 cm非根际土壤中无茎真菌属的相对丰度显著提高,说明覆土越厚,越利于该真菌的生长,将促进其分解土壤有机物,形成多种次生代谢产物
[26]。总体来看,80 cm覆土厚度对玉米根际与非根际土壤真菌群落的影响较大,该厚度更利于复垦土壤质量的提升。
3.3 覆土厚度与玉米根际与非根际土壤真菌代谢功能
利用PICRUSt预测真菌群落功能,发现不同覆土厚度下真菌代谢通路主要表现为生物合成,因子分析表明玉米根际得分值高于非根际,且80 cm覆土厚度下真菌代谢功能得分值最高。这与80 cm的覆土厚度相较于40和120 cm增加了煤矸石填埋区根际土壤真菌多样性相关,因为大多数微生物在土壤中占据特有的生态位并执行不同生态功能,微生物多样性增加通常意味着生态功能增强
[27]。意味着80 cm的覆土厚度不仅可以为玉米生长提供良好的土壤环境,而且对于维持和增强该复垦区土壤生态系统功能有良好的作用。在玉米根际及非根际土壤样品中,unclassified真菌类群均表现出较高的相对丰度,这一现象提示在该复垦区域不同覆土厚度条件下,可能存在大量尚未被分类鉴定的真菌新类群,这些潜在的新类群在土壤生态系统中的功能角色及其生态学意义,仍需通过后续研究深入解析。
3.4 玉米根际与非根际土壤真菌群落与养分的相关性
土壤真菌群落组成随环境变化而改变
[10]。土壤养分如有机质、碳、氮、磷、钾等可能会通过对土壤真菌施加生理限制,进而改变土壤真菌群落结构
[28]。如陈建文等
[10]研究表明土壤有机质与土壤微生物群落组成密切相关,土壤有机质越多,微生物生物量越高,进而可以促进土壤微生物群落发展
[29];孙美美等
[18]的研究也表明土壤有机质是改变土壤真菌群落的关键因子;Chen等
[30]的研究表明土壤中除了有机质,硝态氮含量与真菌群落多样性的相关性也较强;谭小敏等
[31]的研究结果表明土壤有机质、全氮、速效钾含量及pH值是驱动农田玉米根际土壤真菌群落结构的关键环境因子。本研究表明土壤中有机质、全氮和有效磷含量对不同覆土厚度下玉米根际与非根际土壤真菌群落组成的贡献率较高,说明土壤环境资源可利用性的改变会直接导致土壤真菌群落的重组。
4 结论
各覆土厚度下玉米根际真菌ASV数量均高于非根际,子囊菌门(Ascomycota)为玉米根际与非根际土壤中的优势真菌门。80 cm的覆土厚度相较于40和120 cm增加了ASV数量,增加了被孢霉门的相对丰度,标志性差异物种最多,增加了土壤中的真菌群落物种多样性和丰富度,真菌代谢功能最强,得分最高。土壤全氮、有机质和有效磷含量是驱动该复垦地不同覆土厚度下根际与非根际土壤真菌群落变化的重要因素。未来的研究应该更加关注不同覆土厚度下的功能菌群对复垦地玉米产量、生物量等农艺指标的影响,切实提高复垦地地力水平。
山西省科技重大专项计划“揭榜挂帅”项目(202201140601028)
山西省基础研究计划资助项目(202203021221223)
忻州市科技局重点研发计划项目(20240301)