检索项 检索词
  土壤学报  2021, Vol. 58 Issue (4): 1028-1039  DOI: 10.11766/trxb202005110021
0

引用本文  

刘茗, 曹林桦, 刘彩霞, 等. 亚热带4种典型森林植被土壤固碳细菌群落结构及数量特征. 土壤学报, 2021, 58(4): 1028-1039.
LIU Ming, CAO Linhua, LIU Caixia, et al. Characterization of Population and Community Structure of Carbon-Sequestration Bacteria in Soils under Four Types of Forest Vegetations Typical of Subtropical Zone. Acta Pedologica Sinica, 2021, 58(4): 1028-1039.

基金项目

浙江省自然科学基金委(LZ16C160002)资助

通讯作者Corresponding author

徐秋芳, E-mail: xuqiufang@zafu.edu.cn

作者简介

刘茗(1992—), 女, 安徽宿州人, 硕士研究生, 主要从事土壤生物与生物化学研究。E-mail: liumingsunrise@icloud.com
亚热带4种典型森林植被土壤固碳细菌群落结构及数量特征
刘茗1,2, 曹林桦1,2, 刘彩霞1,2, 梁辰飞1,2, 秦华1,2, 陈俊辉1,2, 邵帅1,2, 徐秋芳1,2,3    
1. 浙江省森林生态系统碳循环与固碳减排重点实验室, 浙江农林大学, 浙江临安 311300;
2. 浙江农林大学环境与资源学院, 浙江临安 311300;
3. 浙江省竹资源与高效利用协同创新中心, 浙江农林大学, 浙江临安 311300
摘要:为研究亚热带不同森林植被类型土壤固碳微生物特征及其影响因子,选取毛竹林(Moso banboo groves)、阔叶林(Broad-leaved forest)、杉木林(Chinese fir forest)和马尾松林(Masson pine forest)等4种森林植被为研究对象,以cbbL为固碳细菌指示基因,利用荧光定量PCR(real-time PCR)和MiSeq高通量测序为研究手段。结果表明,4种林分土壤的细菌16S rRNA基因和固碳细菌cbbL基因丰度范围分别是5.40×1010~2.81×1011 copies·g-1干土和4.55×108~3.53×109 copies·g-1干土,其中毛竹林显著高于其它三种林分(P < 0.05);基因丰度显著关联的环境因子是阔叶林土壤的有效磷、不同土层的pH(P < 0.05)。杉木林土壤固碳细菌多样性显著低于其他3种林分(P < 0.05),其亚表层土壤高于表层(P < 0.05);双因子分析表明,林型、土层之间土壤固碳细菌多样性均存在显著或极显著差异。所有土壤具有相似的优势属但相对丰度不同,其中毛竹林和杉木林土壤的甲基化石油杆菌属(Methylibium)和诺卡菌属(Nocardia)占比明显高于阔叶林和马尾松林。冗余分析结果显示,不同林分土壤pH、土壤有机碳、有效磷、全氮差异是影响土壤固碳细菌群落特征形成的主要因素。综上,4种植被对土壤固碳微生物数量及群落结构多样性影响明显,从土壤理化性质、固碳细菌基因丰度、多样性以及结构特征等多维度结果证明,毛竹林对土壤肥力以及固碳细菌影响效果最好,固碳微生物对毛竹林土壤有机质积累贡献大于阔叶林,定量结论有待进一步研究。
关键词固碳细菌    cbbL基因    亚热带森林    实时荧光定量PCR    高通量测序    
Characterization of Population and Community Structure of Carbon-Sequestration Bacteria in Soils under Four Types of Forest Vegetations Typical of Subtropical Zone
LIU Ming1,2, CAO Linhua1,2, LIU Caixia1,2, LIANG Chenfei1,2, QIN Hua1,2, CHEN Junhui1,2, SHAO Shuai1,2, XU Qiufang1,2,3    
1. Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon sequestration, Lin'an, Zhejiang 311300, China;
2. School of Environmental and Resources, Zhejiang Agriculture and Forestry University, Lin'an, Zhejiang 311300, China;
3. Zhejiang Provincial Collaborative Innovation Center for High-Efficiency Utilization of Bamboo Resources, Zhejiang Agriculture and Forestry University, Lin'an, Zhejiang 311300, China
Abstract: 【Objective】To alleviate the global climate change is one of the most important environmental challenges facing mankind. Autotrophic microorganisms, especially those in forest ecosystems, have been reported to have a strong ability to adapt to environmental changes and a high potential to sequestrate carbon. However, how their carbon-sequestrating effects vary with type of forests vegetation remains unclear. To explore this complex relationship, investigations were made of populations and community structures of soil carbon-sequestrating bacteria in the soils of four typical subtropical forests.【Method】Soil samples were collected from the surface (0-20 cm) and subsurface (20-40 cm) soil layers in the four types of forests, that is Moso banboo, Broad-leaved trees, Chinese fir and Masson pine for analysis of cbbL as indicator gene with the aid of quantitative PCR and MiSeq high-throughput sequencing technology.【Result】Abundance of the bacteria with 16S rRNA genes and those with cbbL genes varied in the range of 5.40×1010-2.81×1011 copies·g-1 dry soil and of 4.55×108to 3.53×109copies·g-1dry soil, respectively, in the four types of forest soils, and was the highest in the soil under bamboo groves (P < 0.05); among soil physico-chemical properties, soil available P and soil pH was significantly related to the abundance of cbbL genes in the two soil layers under the broad-leaved forest (P < 0.05). The soil under the Chinese fir forest was significantly lower than the others, and the subsurface soil was higher than the surface soil in diversity of carbon-sequestrating bacteria (P < 0.05). Double factor variance analysis shows that the diversity of soil carbon-sequestrating bacteria varied significantly or ultra-significantly with forest type and soil layer. All the soils had similar dominant species of cbbL-carrying bacteria, which, however, differed in relative abundance beteween the four types of vegetations. The soil of the Moso bamboo groves and the Chinese fir forest was signiticantly higher than that of the Broad-leaved forest and Masson pine forest in proportion of Methylibium and Nocardia. Redundancy analysis (RDA) shows that soil pH, organic carbon, available P and total N were the main factors affecting the formation of community structure of the soil carbon-sequestrating bacteria in the soils.【Conclusion】To sum up, all the findings in this study indicate that population and community structure diversity of the soil carbonsequestrating microbes varied with type of the vegetation. Comparison of the four forest soils in soil physico-chemical property, and genetic abundance, diversity and community structure of the carbon-sequestrating bacteria shows Moso bamboo groves are the best type of vegetation benefitting soil fertility and carbon-sequestrating bacteria. Contribution of the bacteria to accumulation of organic matter was higher in the soil under bamboo groves than under broad-leaved forests, but as to how much, further investigations should be done.
Key words: Carbon-sequestrating bacteria    cbbL gene    Subtropical forests    Real-time quantitative PCR    High-throughput sequencing    

工业飞速发展以及人类活动加剧了全球性温室效应[1],CO2是影响“温室效应”最主要的气体,降低大气CO2浓度有利于缓解温室效应。土壤微生物参与的土壤碳循环与温室气体排放密相关[2],一方面土壤微生物分解有机质和自身呼吸释放CO2,另一方面,自养微生物生长繁殖固定CO2[3]。通过植物或自养微生物将CO2转化成有机化合物的生物固碳过程、提高生态系统的碳固定和储存能力,是减少大气CO2浓度的有效途径[4-5]。研究证明,自养微生物同化CO2是全球生态系统碳循环中的关键微生物过程,对调节大气中二氧化碳浓度具有极重要作用[6],全球陆地土壤的自养微生物每年可捕获大气0.5%~4.1%的CO2、总量达0.6~4.9 Gt C[7]。土壤碳同化量主要来自于土壤自养微生物固定[8],自养微生物广泛分布于不同的生态系统中,具有很强的环境适应能力,可以在多种环境条件下如火山口[9]、海洋深处[10]、极地湖泊[11]等植物无法生存的生境中参与CO2的同化固定。迄今为止发现的5条固碳途径中卡尔文循环是自养微生物固定大气CO2的主要途径,其中1, 5-二磷酸核酮糖羧化酶(RubisCO)是其关键酶[12],研究者们将RubisCO I大亚基的编码基因cbbL基因作为固碳微生物标记物,对不同生态系统土壤固碳微生物开展了大量研究,固碳细菌数量及多样性受土地利用方式、海拔梯度、植被类型等影响[13-16]

综上所述,从整个生物圈的物质、能量流动来看,研究自养微生物的固碳功能具有重要意义,但影响固碳微生物因素却较为复杂。森林是陆地重要的碳汇,全球森林生态系统存储碳量高达650 Gt,其中地下土壤占45%[17]。中国亚热带森林面积大(2.5×108 hm2)、植被类型多样,在全球森林碳汇中发挥着独特作用[18]。因凋落物质量、根系分泌物和养分吸收的差异[19],地上植被种类和多样性的变化能直接或间接影响碳在地下生态系统的分配,并调节影响参与这些过程的土壤微生物[20]。土壤微生物最常利用的碳源是碳水化合物和羧酸,其次是氨基酸,酚酸和聚合物,而胺的利用率最低,土壤微生物群落利用各种碳源的数量及其比例会因植被类型而不同;森林土壤中细菌群落结构与功能多样性与地上植被组成密切相关[21],而土壤固碳细菌占总细菌的0.1%~10%不等[22-23],群落结构和丰度也会受植被类型影响,然而,森林植被类型对土壤固碳微生物影响研究却鲜见报道。

本研究以植物多样性丰富的亚热带不同森林植被类型为对象,选取亚热带典型森林植被毛竹林(Moso banboo groves)、阔叶林(Broad-leaved forest)、杉木林(Chinese fir forest)和马尾松林(Masson pine forest)等4种林分,以卡尔文循环功能基因cbbL为标记基因,利用qPCR和Illumina MiSeq高通量测序技术,探明4种林分土壤固碳细菌丰度、多样性和群落结构差异,为揭示亚热带4种林分对土壤固碳细菌影响机理提供依据。

1 材料与方法 1.1 研究区概况

研究样地位于浙江省杭州市临安区玲珑山(30°14′N,119°42′E),属中亚热带季风气候,热量丰富,雨量丰沛,光照充足,四季分明,具有春多雨、夏湿热、秋气爽、冬干冷的气候特征。多年平均降雨量1 628.6 mm,降水日158 d,全年日照时数1 847 h,无霜期237 d,全年平均气温16.4 ℃。土壤为发育于凝灰岩的黄红壤,具有典型的中亚热带的森林生态系统和森林景观,植被主要包括毛竹、阔叶、杉木、马尾松等典型亚热带植被。其中,阔叶林、杉木林和马尾松林下植被均包含丛生灌木,有蕨类植物,毛竹林下植被很少,仅有稀疏的草本植物。

1.2 样品采集与处理

于2017年10月按照生态控制原则,选择同一成土母质的毛竹林(MB)、阔叶林(BL)、杉木林(CF)和马尾松林(MP)。分别于各个林分的下坡位、中坡位、上坡位确定采样区,在采样区周围10 m2范围内按照三角形确定3个取样点,分别采集表层(0~20 cm)及亚表层(20~40 cm)土壤混匀得到一个土样,每种林分选择4个重复,4种林分共得32个(4×4×2)土样。新鲜土样去除植物残体和大的石头过筛(2 mm)后装入密封袋,将样品放入冰盒带回实验室。样品共分为2份,1份测定含水率后鲜土提取土壤细菌总DNA,并冷冻干燥存放于-70 ℃冰箱,用于其他分子生物学分析;另一份于室内自然风干,充分研磨过筛后用于土壤基本理化性质的测定。

1.3 分析方法

土壤理化性质分析:分析方法参考文献[24]。土壤pH采用酸度计测定,水土比为: 1(m/v);土壤有机碳(SOC)及全氮(TN)使用元素分析仪(ELEMENTAR Vario El Cube)测定;土壤碱解氮(AN)含量采用碱解扩散法测定;有效磷(AP)含量采用盐酸-氟化铵溶液浸提,钼锑抗比色法测定;速效钾(AK)含量采用醋酸铵提取,火焰光度计测定。

土壤总DNA提取:采用Power Soil DNA Isolation Kit试剂盒提取土壤总DNA。测定鲜土的含水率后,换算成0.35 g干土重,称取相应的新鲜土样,按试剂盒说明书进行提取。DNA提取成功后经1%(m/v)的琼脂糖凝胶电泳检测DNA片段大小,并用微量分光光度计检测其浓度和纯度。提取后的DNA样品保存于-40 ℃冰箱。

实时荧光定量PCR(Real-time Quantutative PCR,qPCR),测定固碳功能菌cbbL基因拷贝数引物对为K2f(ACCAYCAAGCCSAAGCTSGG),V2r(GCCTTCSAGCTTGCCSACCRC),得到492 bp-495 bp的扩增产物[25]。测定土壤细菌16S rRNA基因拷贝数通用引物对为338F/518R,16S-F(AGAGT TTGATCMTGGCTCAG),16S-R(GCTGCCTCCCGT AGGAGT)[26]。使用CFX96TM Real-Time System(Bio-Rad,USA)仪器对cbbL基因进行荧光定量PCR扩增。反应体系及程序详见文献[25-26]。

高通量测序:采用同上K2f/V2r引物对,对土壤DNA进行扩增,将PCR产物纯化后送至上海派森诺公司,采用Illumina MiSeq平台对群落DNA片段进行双端(Paired-end)测序。下机数据选择NCBI(National Center for Biotechnology Information)数据库的氨基酸序列进行对比,进行物种注释得到分类鉴定结果。

1.4 数据处理

利用Microsoft Excel 2010软件对数据进行初处理,采用SPSS 21.0软件对数据进行和双因素方差分析(Two-Way ANOVA)和最小显著差数(LSD)法进行显著性和多重比较(α=0.05),采用皮尔逊(Pearson)进行相关性分析。根据OTU表计算α多样性指数。利用Origin软件绘制基因丰度图和群落结构丰度图。使用Canoco 5.0软件对固碳细菌群落结构与环境因子的相关性进行冗余分析(Redundancy analysis,RDA)。

2 结果 2.1 不同林分土壤理化性质

4种林分虽然发育于同一母岩,但土壤理化性质却存在差异(表 1)。毛竹林和阔叶林表层土壤各项指标总体均较高,而杉木林和马尾松林土壤总体较低,尤其是杉木林土壤除pH和碱解氮外的其他指标均为最低。亚表层土壤2组林分之间的差异不如表层明显,而杉木林土壤的各项指标(pH除外)仍然为最低。表层土壤碳氮比总体高于亚表层,但没有显著差异,其中毛竹林两层土壤碳氮比均为最低。双因素检验结果表明,林分对5个土壤理化性质均产生显著影响(P < 0.05),而对碳氮比和有效磷则无影响;土层对土壤有机碳、全氮、碱解氮、速效钾具有显著影响(P < 0.05);林分和土层二者交互作用仅土壤有机质达到显著影响(P < 0.05),说明两者对土壤性质的影响可能不同步、甚至相反。

表 1 不同林分土壤理化性质(均值±标准差) Table 1 Soil physico-chemical propert ties relative to type of forest(mean ± SD)
2.2 不同林分土壤细菌16S rRNA基因与固碳细菌cbbL基因丰度特征

采用荧光实时定量PCR分析,4种林分土壤细菌16S rRNA基因和固碳细菌cbbL基因丰度结果表明(图 1)细菌丰度在5.40×1010~2.81×1011 copies·g-1干土之间(图 1A),固碳细菌丰度在4.55×108~3.53×109 copies·g-1干土之间(图 1B)。毛竹林土壤2种基因丰度在表层、亚表层土壤中均位于四种林分之首,且显著高于其他三种林分(P < 0.05),而其他三种林分细菌数量无显著差异,固碳细菌马尾松林显著高于其他两种林分(P < 0.05)。4种林分土壤cbbL与16S rRNA基因丰度的比值结果为亚表层略高于表层(图 1C);两层土壤比值最高的是马尾松林,其次为毛竹林,马尾松林显著高于阔叶林和杉木林(P < 0.05)。

注:不同小写字母表示差异显著(P<0.05)。Note:Different letters in the same column meant significant difference at 0.05 level. 图 1 不同林分土壤细菌16S rRNA(A)、固碳细菌cbbL基因(B)丰度及其比值(C) Fig. 1 Abundance of 16S rRNA(A)and cbbL genes(B)and their ratio relative to forest stand(C)
2.3 不同林分土壤固碳细菌基因丰度与理化因子的关系

为揭示土壤理化性质对不同林分土壤细菌16S rRNA和固碳细菌cbbL基因丰度的影响,分别以林分和土层为单位,将基因丰度与对应土壤理化性质进行相关性分析,结果表明,细菌和固碳细菌丰度与理化性质相关性总体较差(表 2)。在不同林分中,只有阔叶林土壤的2种基因丰度与有效磷呈显著正相关(P < 0.05),马尾松林土壤的cbbL基因丰度与pH、速效钾呈显著正相关(P < 0.05);在不同土层中,表层及亚表层土壤的两种基因丰度均与pH呈显著正相关(P < 0.05)。

表 2 土壤细菌16S rRNA和固碳细菌cbbL基因功能基因与环境参数之间的皮尔逊相关系数 Table 2 Pearson's correlations(r)of 16S rRNA and cbbL genes abundance with environmental parameters
2.4 不同林分土壤固碳微生物α多样性比较

对4种林分土壤固碳细菌进行高通量测序,毛竹林、阔叶林、杉木林和马尾松林表层土壤分别测得36 598、33 599、35 093和35 129条序列,亚表层土壤测得32 552、33 783、30 381和38 502条有效序列,将测序结果在97%相似度下进行聚类得到OTU的代表序列,并与NCBI数据库比对进行物种注释,置信度阈值设为0.8,认为测序深度已经基本覆盖到样品中所有的物种。经比对鉴定得到4个门、6个纲、12个目、20个科、35个属、40个种的土壤固碳细菌群落信息。

基于OTU的α多样性指数可反应土壤固碳细菌群落结构的多样性。从表层土壤固碳细菌α多样性4个指标看(表 3),杉木林所有指标最低,且均显著低于其他3种林分(P < 0.05),说明杉木林土壤固碳细菌的物种丰富度及多样性极差。亚表层土壤不同林分之间仅香农指数存在差异,其中马尾松林最高、杉木林最低,两者之间存在显著差异(P < 0.05)。马尾松林的所有指标值均为最高,说明马尾松林亚表层土壤固碳细菌的物种丰富度及多样性最高。同一林分不同土层之间的差异结果表明,杉木林土壤亚表层的所有指标都显著大于表层(P < 0.05),双因素检验结果表明,林分及土层对4种α多样性指标均达到了显著影响(P < 0.05),而两因子交互作用仅对香农指数与辛普森指数有显著影响(P < 0.05),且数值小于各单因子,说明两个因子对土壤固碳细菌α多样性影响不同步、甚至相反。

表 3 不同林分土壤固碳细菌cbbL基因多样性指数统计表(均值±标准差) Table 3 CbbL gene diversity of carbon sequestrating bacteria relative to forest stand(mean ± SD)
2.5 不同林分土壤固碳细菌群落相似性比较

为分析4种不同林分土壤固碳细菌群落组成的差异,以固碳细菌OTU种类数为依据,对其分组构建韦恩图分析不同处理之间的“相似关系”(图 2)。图 2A结果表明,天然林(阔叶林和马尾松林)土壤固碳细菌明显高于人工林(毛竹林和杉木林),但毛竹林土壤固碳细菌群落组成特异性强,与其他几种林分的差异最大;4种林分在表层土壤中共有OTU只有28个,占固碳细菌总OTU比例仅为0.5%,说明不同林分之间固碳细菌结构差异较大。亚表层土壤固碳细菌特有OTU数量明显低于表层土壤,不同林分之间的差异没有表层明显。这一结果与4种林分亚表层高达151个共有OTU、以及高占比(2.8%)结果吻合。综上所述,杉木林表层和亚表层土壤固碳细菌的物种数最低,亚表层土壤固碳细菌特有OTU数量大于表层,但表层土壤特有OTU数量及其占比则明显高于亚表层,尤其是毛竹林和杉木林土壤,说明不同林分表层土壤固碳细菌特异性强。

图 2 不同林分土壤固碳细菌OTUs韦恩图 Fig. 2 OTUs Veen of soil carbon sequestrating bacteria relative to forest stand
2.6 不同林分土壤固碳细菌群落物种组成

在属分类水平上(图 3),表层土壤中的毛竹林和杉木林,分属于变形菌门(Proteobacteria)的甲基化石油杆菌属(Methylibium)占比分别高达55.1%、70.85%,为绝对优势菌属;分属于放线菌门(Actinobacteria)的诺卡菌属(Nocardia)为次要优势菌属,分别占21.6%、7.75%;与之相反,阔叶林和马尾松林的诺卡菌属比例高于毛竹林和杉木林,分别占29.65%、26.25%,而甲基化石油杆菌属分别占28.15%、5.15%;变形菌门的红长命菌属(Rubrivivax)占比仅次于诺卡菌属也是优势菌属。在亚表层土壤中,毛竹林土壤诺卡菌属为优势菌属(33.80%),其次为放线菌门的假诺卡氏菌属(Pseudonocardia)(19.86%);杉木林土壤甲基化石油杆菌仍属优势菌群(38.6%)、其次诺卡菌属(20.1%);阔叶林与马尾松林的诺卡菌属分别占28.86%、28.26%,甲基化石油杆菌占28.76%、6.16%;4种林分土壤的假诺卡氏菌属均高于表层土壤。结果进一步表明土壤固碳细菌在人工林与天然林上具有较大差异。

图 3 不同林分土壤固碳细菌在属水平上的群落结构组成分布 Fig. 3 Composition of the soil carbon sequestrating bacteria community structure at genus level relative to forest stand
2.7 固碳细菌群落结构与环境因子的冗余分析

为分析4种林分土壤理化因子对土壤固碳细菌群落结构的影响,选取表层及亚表层土壤具有代表性的优势菌属(聚类大于1%)为物种变量、土壤理化性质为环境变量进行冗余分析,根据蒙特卡洛检验,选取影响较大的环境因子分析。结果表明,表层土壤第一排序轴和第二排序轴分别解释了34.4%和20.7%的变异(图 4A),4种林分样品分散于坐标的4个象限中,说明4种林分土壤固碳细菌的结构差异明显,分布在第一排序轴左侧的毛竹林和杉木林与右侧的阔叶林和马尾松林之间差异明显。第一排序轴与土壤有效磷呈正相关,与pH和全氮呈负相关;第二排序轴与pH、全氮、有效磷成呈正相关。土壤pH(F=4.152,P=0.048)、全氮(F=9.084,P=0.01)有效磷(F=7.317,P=0.01)对土壤固碳细菌属水平达到显著影响。

注:AP:有效磷;TN:全氮;SOC:有机碳。Note:AP:Available P;TN:Totle N;SOC:Soil organic carbon. 图 4 土壤固碳细菌在属水平上与环境因子的冗余分析 Fig. 4 Redundancy analysis of soil carbon sequestrating bacteria at the genus level with environmental parameters relative to forest stand

根据蒙特卡洛检验,选取亚表层土壤中(图 4B)影响较大的3个环境因子分析,结果表明,第一排序轴和第二排序轴分别解释了77.7%和3.7%的变异。第一排序轴上的群落结构与pH呈正相关,与有机碳、有效磷呈负相关;第二排序轴与pH呈负相关,与有机碳、有效磷呈负相关;土壤pH(F=2.574,P=0.044)达到显著水平,但有效磷(F=2.558,P=0.0560)、有机碳(F=1.882,P=0.154)未达到显著水平。

3 讨论 3.1 不同林分土壤细菌与固碳细菌基因丰度差异及其原因分析

通过qPCR分析亚热带典型森林植被土壤微生物16S rRNA基因与cbbL基因拷贝数得知,亚热带森林土壤中存在相当数量的土壤细菌与固碳细菌,对土壤碳汇贡献较大。本研究中4种林分土壤16S rRNA基因与cbbL基因数量普遍高于其他学者在湿地[15]、稻田[22]、旱地[23]、盐沼地[27]的研究结果。一般而言,土壤有机质与细菌的数量呈正相关[28],森林中凋落物的分解、积累与补充速率远高于农田生态系统,森林土壤丰富的有机质为细菌(包括固碳细菌)的生长繁殖提供能量和养分。土壤固碳细菌是一类利用CO2而非有机碳为碳源的自养微生物[29],土壤固碳细菌以兼性自养微生物为主[30],说明其还可同时利用有机碳作为碳源,有机质分解产生的CO2亦可为其提供底物。4种林分土壤固碳细菌占总细菌的比值(0.91%~2.24%)与相同气候带稻田[22]、旱地[23]相差不大,说明土地利用方式、生态类型、气候区域等对土壤总细菌和固碳细菌丰度影响是同步的。

不同生态系统因为土壤环境因子的巨大差异导致土壤固碳细菌丰度差异[22, 27],本研究4种林分土壤发育于同一母质,土壤基础理化性质一致,但由于不同植物对养分的选择性吸收、根系分泌物与凋落物分解差异进而影响土壤性质[16, 19]。因此,4种植被长期生长导致其林下土壤性质差异显著(表 1)。4种林分之间土壤固碳细菌基因丰度和总细菌丰度存在差异,毛竹林最高、阔叶林最低,毛竹林土壤显著高于其他3种林分,推测毛竹林土壤具有较强的固碳潜力,对土壤碳汇贡献较大。在天目山毛竹入侵阔叶林的研究中同样观察到阔叶林土壤固碳细菌基因丰度显著低于毛竹林的结果[31]。毛竹林地下竹鞭系统极其发达、根系分泌物丰富、老鞭死亡腐烂、凋落物易降解等为土壤提供大量碳源和养分。通常认为,土壤有机碳水平与固碳细菌丰度正相关[13],土壤养分水平高对微生物生长有利[32],本研究中毛竹、杉木与马尾松林均符合此规律;然而,阔叶林表现异常,土壤有机碳含量和养分水平最高、固碳细菌丰度却最低,说明本研究的4种土壤有机质含量不是影响16S rRNA与cbbL基因的最主要因素。虽然土壤固碳基因丰度与有机碳非正相关的结果也有报道,如西北干旱地区不同植被土壤有机质差异显著却并未引起土壤固碳细菌丰度明显差异[23],可能是因其土壤缺水限制其他因子作用,但本研究区降水充沛,为何土壤有机碳和养分含量丰富的阔叶林却出现基因丰度最低现象似乎难以解释。基因丰度与土壤理化性质的相关性分析(表 2)表明阔叶林土壤有效磷是影响其基因丰度的重要因子,但尚无法解释阔叶林土壤基因丰度最低的结果;表层及亚表层土壤的两种基因丰度均与pH呈显著正相关(P < 0.05)的结果也许可解释阔叶林的异常现象。大空间尺度下土壤pH是影响微生物多样性的关键因子,在一定范围内pH与细菌多样性呈正比[19, 22]。pH会间接影响有机质分解和养分释放,从而导致对生境变化较敏感的微生物种群丰度发生变化,同时提高某些有害金属元素如铝、锰等活性,对生物造成危害[32];高浓度H+会抑制细菌生长[33],较低的pH还会改变营养元素的形态等[34]。综观文献报道可知,土壤理化性质对基因丰度的影响比较复杂[30-31],有时候几种因子的作用可能相互抵消从而掩盖某项或几项土壤理化性质与基因丰度的密切关系,尤其是对于土壤理化性质差异不悬殊、且样本数量不大的情况,基因丰度与理化性质之间的相关性更加难以达到显著水平。因此,不同研究者所选择的对象不同,难以获得驱动土壤固碳细菌变化因子的一致规律。

3.2 不同林分土壤固碳细菌群落多样性和结构差异及其环境驱动因子分析

对固碳细菌功能基因cbbL进行高通量测序,全面分析了4森林植被土壤固碳细菌的群落结构与多样性特征。α多样性主要关注局域均匀生境下的物种数目,适合本研究固碳细菌多样性的描述。由α多样性指数表(表 3)可以推测,杉木和毛竹林土壤固碳细菌存在着明显的优势种群,一方面可能与纯林植被多样性低有关;另一方面是杉木林长期种植分泌出的酚类物对土壤理化性质和微生物产生不良影响[35],而毛竹林的特殊性可能也与其发达的地下竹鞭系统有关。4种林分亚表层土壤多样性均高于表层,类似的结果在棉花和花生等农田土壤也有报道[16],可能与亚表层根系密集分布、丰富的根系分泌物和呼吸产生大量CO2有利于固碳细菌生长繁殖有关。

不同林分之间固碳细菌丰度和多样性差异的根本原因是群落结构不同,从韦恩图结果显示的林分之间差异(图 3),毛竹林土壤固碳细菌群落组成特异性强、与其他几种林分的差异最大;杉木林土壤固碳细菌丰富度最低;而阔叶和马尾松林土壤虽然OTU总数多,但特有OTU比例却最少,说明两者土壤固碳细菌种类丰富。李永春等[36]利用DGGE对4种林分土壤氨氧化细菌菌群落结构研究发现,林分之间的相似关系与本研究完全不同。虽然固碳细菌和氨氧化细菌均属于细菌,但固碳细菌是指所有自养细菌的总称,种类多,而氨氧化细菌是特殊生理功能的群体,种类少,因此其对环境变化的响应不同。

本研究中4种林分表层土壤占绝对优势的甲基化石油杆菌属,在桂林毛村岩溶的岩溶区农田土壤中为第二优势属[37],甲基化石油杆菌属是杉木林的绝对优势种群,表层土壤中占比高达70.85%,马尾松林则占比最低(38.60%);优势属的明显差异可能受到植物凋落物以及根系分泌物的影响较大。红长命菌属在阔叶林、马尾松林中大量存在,但在毛竹林与杉木林却极少发现,可能该菌与自然林中的某些植物有关。土壤性质对固碳细菌群落结构差异的影响是主要原因[13],从RDA结果判断驱动土壤固碳细菌变化的因子发现(图 4),表层土壤固碳细菌群落结构受有效磷、全氮、pH影响;丰富的氮、磷营养塑造了阔叶林表层土壤固碳细菌群落结构特征,较高的氮素水平和pH则是形成毛竹林土壤固碳细菌群落结构主要因素,其他林分依次类推。亚表层主要受有机碳、有效磷、pH调控,虽然4种林分在坐标图中的分布区沿逆时针旋转90度,但驱动不同林分土壤固碳细菌群落结构形成的环境与表层土壤基本一致,土壤pH对毛竹和杉木土壤固碳细菌结构形成密切相关。虽然不同研究者得到的驱动土壤固碳细菌群落结构变化的主要因子不尽相同,但土壤pH是较为一致的驱动因子[14, 30]

4 结论

毛竹林、阔叶林、杉木林和马尾松林4种森林植被土壤中均存在相当数量的细菌及固碳细菌,其中毛竹林显著高于其他3种林分,阔叶林土壤有效磷以及不同土层的pH分别与对应的细菌与固碳细菌基因丰度显著相关。杉木土壤固碳细菌多样性显著低于其他3种林分;双因子分析显示,林型、土层之间均存在显著差异。4种林分土壤固碳细菌群落结构与农田、草地和海洋等生态系统有较大差异,甲基化石油杆菌属和诺卡菌属为4种林分共同的优势菌属;毛竹和杉木两种人工林土壤固碳细菌群落结构组成较为相似。RDA结果显示,不同林分土壤pH、有机碳、有效磷、全氮差异是驱动土壤固碳细菌群落特征形成的主要因素。毛竹林土壤肥力最好,固碳细菌丰度最高,不同林分之间土壤微生物对土壤有机质的积累的贡献存在差异,确切定量结论有待进一步研究。

参考文献
[1]
Adams J M, Piovesan G. Uncertainties in the role of land vegetation in the carbon cycle[J]. Chemosphere, 2002, 49(8): 805-819. DOI:10.1016/S0045-6535(02)00382-X (0)
[2]
Wu J S, Li Y, Tong C L, et al. The key geo-biochemical processes of the long-term carbon sequestration and its mechanisms in the subtropical paddy soils (In Chinese)[J]. Research of Agricultural Modernization, 2018, 39(6): 895-906. [吴金水, 李勇, 童成立, 等. 亚热带水稻土碳循环的生物地球化学特点与长期固碳效应[J]. 农业现代化研究, 2018, 39(6): 895-906.] (0)
[3]
Ge T D, Wu X H, Chen X J, et al. Microbial phototrophic fixation of atmospheric CO2 in China subtropical upland and paddy soils[J]. Geochimica et Cosmochimica Acta, 2013, 113: 70-78. DOI:10.1016/j.gca.2013.03.020 (0)
[4]
Zhang W, Li Y L. Research status and prospects of carbon dioxide sequestration technique (In Chinese)[J]. Environmental Pollution and Control, 2006, 28(12): 950-953. DOI:10.3969/j.issn.1001-3865.2006.12.021 [张炜, 李义连. 二氧化碳储存技术的研究现状和展望[J]. 环境污染与防治, 2006, 28(12): 950-953.] (0)
[5]
Parry M, Canziani O, Palutikof J, et al. Climate change 2007: Impacts, adaptation and vulnerability//Contributionof working group Ⅱ to the fourth assessment report of the Intergovernmental Panel on climate change[EB/OL]. 2007. (0)
[6]
Hügler M, Sievert S M. Beyond the Calvin cycle: Autotrophic carbon fixation in the ocean[J]. Annual Review of Marine Science, 2011, 3(1): 261-289. DOI:10.1146/annurev-marine-120709-142712 (0)
[7]
Falkowski P. The global carbon cycle: A test of our knowledge of earth as a system[J]. Science, 2000, 290(5490): 291-296. DOI:10.1126/science.290.5490.291 (0)
[8]
Jian Y, Zeng G J, Zhou P, et al. Input and distribution of autotrophic microbe-assimilated carbon in humus and aggregate fractions of soils (In Chinese)[J]. Research of Environmental Sciences, 2014, 27(12): 1499-1504. [简燕, 曾冠军, 周萍, 等. 自养微生物同化碳在土壤腐殖质组分及团聚体中的分配[J]. 环境科学研究, 2014, 27(12): 1499-1504.] (0)
[9]
Nanba K, King G M, Dunfield K. Analysis of facultative lithotroph distribution and diversity on volcanic deposits by use of the large subunit of ribulose 1, 5-bisphosphate carboxylase/oxygenase[J]. Applied and Environmental Microbiology, 2004, 70(4): 2245-2253. DOI:10.1128/AEM.70.4.2245-2253.2004 (0)
[10]
Nanba K, King G M, Dunfield K. Analysis of facultative lithotroph distribution and diversity on volcanic deposits by use of the large subunit of ribulose 1, 5-bisphosphate carboxylase/oxygenase[J]. Applied and Environmental Microbiology, 2004, 70(4): 2245-2253. DOI:10.1128/AEM.70.4.2245-2253.2004 (0)
[11]
Tourova T P, Kovaleva O L, Sorokin D Y, et al. Ribulose-1, 5-bisphosphate carboxylase/oxygenase genes as a functional marker for chemolithoautotrophic halophilic sulfur-oxidizing bacteria in hypersaline habitats[J]. Microbiology, 2010, 156(7): 2016-2025. DOI:10.1099/mic.0.034603-0 (0)
[12]
Kong W D, Ream D C, Priscu J C, et al. Diversity and expression of RubisCO genes in a perennially ice-covered Antarctic lake during the polar night transition[J]. Applied and Environmental Microbiology, 2012, 78(12): 4358-4366. DOI:10.1128/AEM.00029-12 (0)
[13]
Tabita F R. Molecular and cellular regulation of autotrophic carbon dioxide fixation in microorganisms[J]. Microbiological Reviews, 1988, 52(2): 155-189. DOI:10.1128/mr.52.2.155-189.1988 (0)
[14]
Yuan H Z, Qin H L, Liu S L, et al. Abundance and composition of CO2 fixating bacteria in relation to long-term fertilization of paddy soils (In Chinese)[J]. Acta Ecologica Sinica, 2012, 32(1): 183-189. [袁红朝, 秦红灵, 刘守龙, 等. 长期施肥对稻田土壤固碳功能菌群落结构和数量的影响[J]. 生态学报, 2012, 32(1): 183-189.] (0)
[15]
Gao J, Muhanmmad S, Yue L Y, et al. Changes in CO2-fixing microbial community characteristics with elevation and season in alpine meadow soils on the northern Tibetan Plateau (In Chinese)[J]. Acta Ecologica Sinica, 2018, 38(11): 3816-3824. [高静, Said Muhanmmad, 岳琳艳, 等. 藏北高原草甸土壤固碳微生物群落特征随海拔和季节的变化[J]. 生态学报, 2018, 38(11): 3816-3824.] (0)
[16]
Cao X B, Lin D, Cai L, et al. Effects of different vegetation communities on soil carbon fraction, RubisCO activity and cbbL genes in Nanjishan Wetland of Poyang Lake (In Chinese)[J]. Acta Pedologica Sinica, 2017, 54(5): 1269-1279. [曹煦彬, 林娣, 蔡璐, 等. 鄱阳湖南矶山湿地不同植被类型对土壤碳组分、羧化酶及cbbL基因的影响[J]. 土壤学报, 2017, 54(5): 1269-1279.] (0)
[17]
Tolli J, King G M. Diversity and structure of bacterial chemolithotrophic communities in pine forest and agroecosystem soils[J]. Applied and Environmental Microbiology, 2005, 71(12): 8411-8418. DOI:10.1128/AEM.71.12.8411-8418.2005 (0)
[18]
Wang S. Research progress of carbon cycling of forest ecosystem (In Chinese)[J]. Journal of Anhui Agricultural Sciences, 2013, 41(4): 1560-1563. DOI:10.3969/j.issn.0517-6611.2013.04.062 [汪森. 森林生态系统碳循环研究进展[J]. 安徽农业科学, 2013, 41(4): 1560-1563.] (0)
[19]
An X J, Li P, Dai W, et al. The variation characteristics of soil organic carbon and its relationship with soil properties in typical subtropical plantations (In Chinese)[J]. Chinese Agricultural Science Bulletin, 2012, 28(22): 53-58. [安晓娟, 李萍, 戴伟, 等. 亚热带几种林分类型土壤有机碳变化特征及与土壤性质的关系[J]. 中国农学通报, 2012, 28(22): 53-58.] (0)
[20]
Haichar F Z, Marol C, Berge O, et al. Plant host habitat and root exudates shape soil bacterial community structure[J]. The ISME Journal, 2008, 2(12): 1221-1230. DOI:10.1038/ismej.2008.80 (0)
[21]
Bardelli T, Gómez-Brandón M, Ascher-Jenull J, et al. Effects of slope exposure on soil physico-chemical and microbiological properties along an altitudinal climosequence in the Italian Alps[J]. Science of the Total Environment, 2017, 575: 1041-1055. DOI:10.1016/j.scitotenv.2016.09.176 (0)
[22]
Wu Z Y, Lin W X, Chen Z F, et al. Characteristics of soil microbial community under different vegetation types in Wuyishan National Nature Reserve, East China (In Chinese)[J]. Chinese Journal of Applied Ecology, 2013, 24(8): 2301-2309. [吴则焰, 林文雄, 陈志芳, 等. 武夷山国家自然保护区不同植被类型土壤微生物群落特征[J]. 应用生态学报, 2013, 24(8): 2301-2309.] (0)
[23]
Liu Q, Wei X M, Wu X H, et al. Characteristic of abundances and diversity of carbon dioxide fixation microbes in paddy soils (In Chinese)[J]. Environmental Science, 2017, 38(2): 760-768. [刘琼, 魏晓梦, 吴小红, 等. 稻田土壤固碳功能微生物群落结构和数量特征[J]. 环境科学, 2017, 38(2): 760-768.] (0)
[24]
Tang Z X, Fan F L, Wan Y F, et al. Abundance and diversity of RuBisCO genes responsible for CO2 fixation in arid soils of northwest China[J]. Pedosphere, 2015, 25(1): 150-159. DOI:10.1016/S1002-0160(14)60085-0 (0)
[25]
Lu R K. Analytical methods for soil and agro-chemistry (In Chinese). Beijing: China Agricultural Science and Technology Press, 2000. [鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000.] (0)
[26]
Yuan H Z, Ge T D, Zou S Y, et al. Effect of land use on the abundance and diversity of autotrophic bacteria as measured by ribulose-1, 5-biphosphate carboxylase/oxygenase(RubisCO)large subunit gene abundance in soils[J]. Biology and Fertility of Soils, 2013, 49(5): 609-616. DOI:10.1007/s00374-012-0750-x (0)
[27]
Fierer N, Jackson J A, Vilgalys R, et al. Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays[J]. Applied and Environmental Microbiology, 2005, 71(7): 4117-4120. DOI:10.1128/AEM.71.7.4117-4120.2005 (0)
[28]
Chaudhary D R, Gautam R K, Yousuf B, et al. Nutrients, microbial community structure and functional gene abundance of rhizosphere and bulk soils of halophytes[J]. Applied Soil Ecology, 2015, 91: 16-26. DOI:10.1016/j.apsoil.2015.02.003 (0)
[29]
刘达文. 长期不同施肥处理下南方典型水稻土有机碳积累的微生物学机制初探[D]. 南京: 南京农业大学, 2011.
Liu D W. A basic study on the microbial mechanism of soil organic carbon accumulation structures under different fertilization treatments of long-term Argo-ecosystem experiments[D]. Nanjing: Nanjing Agricultural University, 2011. (0)
[30]
Yuan H Z, Qin H L, Liu S L, et al. Advances in research of molecular ecology of carbon fixation microorganism (In Chinese)[J]. Scientia Agricultura Sinica, 2011, 44(14): 2951-2958. [袁红朝, 秦红灵, 刘守龙, 等. 固碳微生物分子生态学研究[J]. 中国农业科学, 2011, 44(14): 2951-2958.] (0)
[31]
Liu C X, Zhou Y, Xu Q F, et al. Effects of intensive management on the community structure and diversity of CO2-assimilating bacteria in a Phyllostachys pubescens stand (In Chinese)[J]. Acta Ecologica Sinica, 2018, 38(21): 7819-7829. [刘彩霞, 周燕, 徐秋芳, 等. 毛竹林集约经营对土壤固碳细菌群落结构和多样性的影响[J]. 生态学报, 2018, 38(21): 7819-7829.] (0)
[32]
梁雪. 毛竹入侵阔叶林对土壤固碳功能菌群落特征的影响及其机制[D]. 杭州: 浙江农林大学, 2017.
Liang X. Effects of Phyllostachys pubscens invasion of native broadleaf forest on community characteristics of soil CO2-fixing bacteria and its mechanism[D]. Hangzhou: Zhejiang A & F University, 2017. (0)
[33]
Shen P, Chen X D. Microbiology (In Chinese). Beijing: Higher Education Press, 2009. [沈萍, 陈向东. 微生物学[M]. 北京: 高等教育出版社, 2009.] (0)
[34]
Rousk J, Brookes P C, Bååth E. Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization[J]. Applied and Environmental Microbiology, 2009, 75(6): 1589-1596. DOI:10.1128/AEM.02775-08 (0)
[35]
Stockdale E A, Shepherd M A, Fortune S, et al. Soil fertility in organic farming systems - fundamentally different?[J]. Soil Use and Management, 2006, 18: 301-308. DOI:10.1111/j.1475-2743.2002.tb00272.x (0)
[36]
Yang Y S, Yu X T, Qiu R H, et al. Variation of rhizospheric soil fertility under different Chinese fir rotations (In Chinese)[J]. Chinese Journal of Applied and Environmental Biology, 1999, 5(3): 254-258. [杨玉盛, 俞新妥, 邱仁辉, 等. 不同栽杉代数根际土壤肥力及生物学特性变化[J]. 应用与环境生物学报, 1999, 5(3): 254-258.] (0)
[37]
Li Y C, Liu B R, Guo S, et al. Characteristics of soil mmonia-oxidation microbial communities in different subtropical forests, China (In Chinese)[J]. Chinese Journal of Applied Ecology, 2014, 25(1): 125-131. [李永春, 刘卜榕, 郭帅, 等. 亚热带不同林分土壤氨氧化菌群落特征[J]. 应用生态学报, 2014, 25(1): 125-131.] (0)
[38]
Zhang S S, Jin Z J, Jia Y H, et al. Community structure of CO2-fixing soil bacteria from different land use types in Karst areas (In Chinese)[J]. Environmental Science, 2019, 40(1): 412-420. [张双双, 靳振江, 贾远航, 等. 岩溶地区不同土地利用方式土壤固碳细菌群落结构特征[J]. 环境科学, 2019, 40(1): 412-420.] (0)