免费文献传递   相关文献

Bacterial community structure and diversity in soils of different forest ages and types in Baotianman forest, Henan Province, China.

河南宝天曼不同林龄与林型森林土壤的细菌群落结构与多样性


采用高通量16S rRNA标签测序法,比较了地处北亚热带与暖温带过渡带的宝天曼自然保护区不同林龄与林分类型的土壤细菌群落结构及多样性.结果表明: 宝天曼森林土壤细菌以变形菌门(29%)、酸杆菌门(18.5%)、疣微菌门(10%)等为主,共检测到60门1209属,优势属主要有疣微菌门的DA101(6.3%)、酸杆菌门的Acidobacteria2(5.9%)和Candidatus Solibacter(2.9%)、泉古菌门的Candidatus Nitrososphaera(2.6%)等.不同林龄和林分类型土壤分别具有特有的种属组成及高丰度和低丰度种属.林龄与林分类型都对土壤微生物群落结构影响显著,且林分类型的影响大于林龄.80年林龄的锐齿栎土壤菌群多样性在不同林龄和林分类型中均最低.pH、土壤全氮、有机碳等是不同林龄及林分类型下土壤菌群结构变化的重要影响因子.

To compare the microbial compositions and diversities in soils of different forest ages and types in Baotianman forest, Henan Province, China, genomic DNA of forest soils was extracted for amplifying the 16S rRNA V4 hyper variable region by PCR and sequencing by Illumina MiSeq. The BIPES, UCHIME and QIIME were employed to analyze the soil bacterial community. It was shown that 60 phyla were identified, with Proteobacteria, Acidobacteria, and Verrucomicrobia representing the most dominant lineages and accounting for 29%, 18.5% and 10% of all sequences, respectively. At the genus level, 1209 genera were identified, the most abundant phylotypes were DA101(6.3%), Acidobacteria2 (5.9%), Candidatus Solibacter (2.9%) and Candidatus Nitrososphaera(2.6%). Different forest age and type soil samples had unique compositions and specific high and rare genus. Forest type and age both impacted the soil microbial community structure, and the influence of the former was stronger than the latter. The soil microbial diversity of the 80yearold  Quercus aliena forest was the lowest among all age and type forest soil samples. Soil pH, soil nitrogen and organic carbon contents were the most important factors affecting soil bacterial community structure.


全 文 :河南宝天曼不同林龄与林型森林土壤的
细菌群落结构与多样性∗
白晓旭1,2  史荣久1  尤业明3  盛华芳4  韩斯琴1  张  颖1∗∗
( 1中国科学院沈阳应用生态研究所污染生态与环境工程重点实验室, 沈阳 110016; 2中国科学院大学, 北京 100049; 3北京林
业大学森林培育与保护教育部重点实验室, 北京 100083; 4南方医科大学公共卫生与热带医学学院环境卫生学系, 广州
510515)
摘  要  采用高通量 16S rRNA标签测序法,比较了地处北亚热带与暖温带过渡带的宝天曼
自然保护区不同林龄与林分类型的土壤细菌群落结构及多样性.结果表明: 宝天曼森林土壤
细菌以变形菌门(29%)、酸杆菌门(18.5%)、疣微菌门(10%)等为主,共检测到 60 门 1209
属,优势属主要有疣微菌门的 DA101(6.3%)、酸杆菌门的 Acidobacteria⁃2(5.9%)和 Candidatus
Solibacter(2.9%)、泉古菌门的 Candidatus Nitrososphaera(2.6%)等.不同林龄和林分类型土壤
分别具有特有的种属组成及高丰度和低丰度种属.林龄与林分类型都对土壤微生物群落结构
影响显著,且林分类型的影响大于林龄.80 年林龄的锐齿栎土壤菌群多样性在不同林龄和林
分类型中均最低.pH、土壤全氮、有机碳等是不同林龄及林分类型下土壤菌群结构变化的重要
影响因子.
关键词  宝天曼暖温带森林; 土壤微生物群落结构; 16S rRNA; Illumina Miseq
文章编号  1001-9332(2015)08-2273-09  中图分类号  Q938.1  文献标识码  A
Bacterial community structure and diversity in soils of different forest ages and types in Bao⁃
tianman forest, Henan Province, China. BAI Xiao⁃xu1,2, SHI Rong⁃jiu1, YOU Ye⁃ming3,
SHENG Hua⁃fang4, HAN Si⁃qin1, ZHANG Ying1 ( 1Key Laboratory of Pollution Ecology and Envi⁃
ronmental Engineering, Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang
110016, China; 2University of Chinese Academy of Sciences, Beijing 100049, China; 3Ministry of
Education Key Laboratory for Silvia Culture and Conservation, College of Forest Science, Beijing
Forestry University, Beijing 100083, China; 4Department of Environmental Health, School of Public
Health and Tropical Medicine, Southern Medical University, Guangzhou 510515, China) . ⁃Chin. J.
Appl. Ecol., 2015, 26(8): 2273-2281.
Abstract: To compare the microbial compositions and diversities in soils of different forest ages and
types in Baotianman forest, Henan Province, China, genomic DNA of forest soils was extracted for
amplifying the 16S rRNA V4 hyper variable region by PCR and sequencing by Illumina MiSeq. The
BIPES, UCHIME and QIIME were employed to analyze the soil bacterial community. It was shown
that 60 phyla were identified, with Proteobacteria, Acidobacteria, and Verrucomicrobia represen⁃
ting the most dominant lineages and accounting for 29%, 18.5% and 10% of all sequences, respec⁃
tively. At the genus level, 1209 genera were identified, the most abundant phylotypes were DA101
(6.3%), Acidobacteria⁃2 (5. 9%), Candidatus Solibacter (2. 9%) and Candidatus Nitrososphaera
(2.6%). Different forest age and type soil samples had unique compositions and specific high and
rare genus. Forest type and age both impacted the soil microbial community structure, and the influ⁃
ence of the former was stronger than the latter. The soil microbial diversity of the 80⁃year⁃old Quercus
aliena forest was the lowest among all age and type forest soil samples. Soil pH, soil nitrogen and or⁃
ganic carbon contents were the most important factors affecting soil bacterial community structure.
Key words: Baotianman warm temperate forest; soil bacterial community; 16S rRNA; Illumina
Miseq.
∗国家重点基础研究发展规划项目(2011CB403205)和国家自然科学基金重点项目(41330530)资助.
∗∗通讯作者. E⁃mail: yzhang@ iae.ac.cn
2014⁃10⁃01收稿,2015⁃05⁃06接受.
应 用 生 态 学 报  2015年 8月  第 26卷  第 8期                                                         
Chinese Journal of Applied Ecology, Aug. 2015, 26(8): 2273-2281
    森林土壤微生物是森林生态系统的重要组成部
分,在土壤有机碳库储存、养分循环、缔结地上与地
下生物学过程间的交互作用、植物多样性维持等方
面发挥着至关重要的作用[1] .土壤微生物的种群结
构与功能是评估生态系统稳定性的重要指标[2-3],
在森林生态系统多样性与功能恢复中扮演着重要角
色[4] .已有研究显示,微生物的多样性在景观与区域
尺度上受环境因子,如森林植被类型、土壤类型、营
养成分、土壤 pH 等的影响显著[5-8],但在生态系统
恢复与重建过程中土壤微生物群落结构及演化规
律等研究仍十分欠缺[9] .目前研究主要集中在地
上生物多样性、植被覆盖度和土壤理化性质等方
面[10-12],而对地下土壤微生物结构与功能的研究不
多,尤其是森林生态系统恢复演化过程对微生物种
群结构组成的影响研究更为缺乏.目前,土壤微生物
研究手段多基于变性梯度凝胶电泳(DGGE)、磷脂
脂肪酸分析( PLFA)等传统微生物分子生态学方
法[13],这些手段的微生物通量和分辨率均不高,检
测的灵敏度也不够.近几年来,16S rRNA 基因高通
量测序技术在准确揭示微生物种群结构方面提供了
强有力的手段[14-16],其中 454、Illumina 和 Ion torrent
是目前应用最为广泛的高通量测序平台.454 和 ABI
测序技术受到高成本和高错误率的制约,而 Illumina
的 MiSeq测序系统采用成熟的 TruSeq 边合成边测
序技术,大大降低了错误的发生,其较低的测序成
本、简便的操作流程使其逐渐在微生物生态学研究
中被广泛应用[17-18] .
宝天曼自然保护区地处暖温带向北亚热带的过
渡带和第二阶梯向第三阶梯的过渡地带,属季风型
大陆性气候,是我国同纬度生态结构保存最完整和
植物群落最复杂的地区,拥有独特的过渡带生态系
统及丰富的物种多样性[19-20] .本研究选取河南伏牛
山地区宝天曼过渡带森林生态系统中具有不同林龄
及不同林分组成的林地土壤为研究对象,通过基于
16S rRNA基因的高通量标签测序法对表层土壤细
菌群落结构进行分析,旨在进一步揭示该区域森林
生态系统土壤微生物多样性特征,为森林生态系统
管理提供微生物生态学层面的理论依据.
1  研究地区与研究方法
1􀆰 1  研究区概况
研究样地位于伏牛山地区宝天曼自然保护区.
宝天曼地处河南省西南部,秦岭东段,伏牛山南坡,
地理坐标为 33°25′—33°33′ N,110°53′—112° E.属
季风型大陆性气候,四季分明.年均温 15.2 ℃,年降
雨量 885.6 mm,年蒸发量 991.6 mm,相对湿度 68%,
极端最高温 41.2 ℃,极端最低温-14.8 ℃,≥10 ℃
的年积温 4200~4900 ℃,无霜期高山区 160 d,低山
区 227 d.植被区划属暖温带落叶阔叶林地带,以栎
类为优势树种[13] .
选取 3个优势树种均为锐齿栎(Quercus aliena)
但林龄不同的样地(33°30′ N,111°56′ E),来比较
研究林龄对土壤微生物的影响,分别为:锐齿栎幼龄
林(林龄约 40 年,文中表示为 QAF40);锐齿栎中龄
林(林龄约 80年,文中表示为 QAF80)和锐齿栎老龄
林(林龄约 160 年,文中表示为 QAF160).同时,选取
了该地区与锐齿栎中龄林林龄相同的均为 80年林龄
的其他 3种森林类型,比较不同林型对土壤微生物的
影响,分别为:短柄枹林[优势树种为短柄枹(Q. glan⁃
dulifera),文中表示为 QGB],栓皮栎林[优势树种为
栓皮栎(Q. variabilis),文中表示为 QV],锐齿栎 /华山
松混交林[优势树种为锐齿栎和华山松(Pinus arman⁃
dii),文中表示为 PA].每个样地分别设置 3 个重复.
样地基本信息及土壤理化性质见表 1.
1􀆰 2  研究方法
1􀆰 2􀆰 1样品采集  在每个样地设 3个 20 m×20 m的
重复样方,使用内径为 5 cm的土钻在每个样方内从
“S”形随机设置的 5 个取样点中取 0 ~ 5 cm 表层土
壤样品,混合均匀后装入无菌自封袋,过 2 mm 土筛
后于-20 ℃冰箱中保存.
1􀆰 2􀆰 2 DNA 提取与 16S rRNA 基因扩增   选用
Powersoil DNA Isolation Kit(MoBio,USA)提取土壤
总 DNA,经琼脂糖凝胶电泳定性检测后,再采用
Nanodrop 2000(Thermo,USA)测定其浓度及纯度,最
后置于-20 ℃冰箱中保存备用.
对每份样品,采用引物 515F / 806R 扩增 16S
rRNA基因的 V4 高变区片段.引物序列[21]为 bar⁃
code⁃GT⁃515F ( GCCAGCMGCCGCGGTAA ) 和 bar⁃
code⁃806R (CCGGACTACHVGGGTWTCTAAT). PCR
条件如下:94 ℃预变性 2 min, 94 ℃变性 30 s,57 ℃
退火 30 s,72 ℃延伸 30 s,25次循环后最终在 72 ℃
终延伸 5 min.每个 25 μL 反应体系都包括 2.5 μL
TaKaRa 10×ExTaq 缓冲液(不含 Mg2+),2 μL dNTP
Mix(2.5 mmol·L-1),1.5 μL Mg2+(2.5 mmol·L-1),
0􀆰 25 μL TaKaRa Ex Taq DNA聚合酶(2.5 U),1 μL模
板 DNA,浓度为 10 μmol·L-1的上游 barcode 引物
515F和下游引物 806R 各 0.5 μL,超纯水(ddH2O)
补齐.
4722 应  用  生  态  学  报                                      26卷
表 1  宝天曼样地基本信息
Table 1  Basic information of Baotianman sampling sites
林分类型
Forest type
海拔
Altitude
(m)
坡度
Slope
(°)
林龄
Stand age
(a)
pH 全氮
TN
(g·kg-1)
有机碳
Organic C
(g·kg-1)
锐齿栎幼龄林
40⁃year⁃old Quercus aliena forest
1383~1416 15~20 40 4.5±0.2 4.2±0.0 57.5±1.0
锐齿栎中龄林
80⁃year⁃old Quercus aliena forest
1368~1412 14~16 80 4.6±0.1 4.2±0.0 63.4±2.8
锐齿栎老龄林
160⁃year⁃old Quercus aliena forest
1398~1411 15~20 160 4.5±0.2 4.5±0.0 60.9±2.5
短柄枹林
Quercus glandulifera forest
1392~1425 12~15 80 4.2±0.2 2.7±0.0 46.6±1.7
栓皮栎林
Quercus variabilis forest
1096~1129 12~16 80 4.1±0.1 2.4±0.1 42.0±1.4
锐齿栎⁃华山松混交林
Mixed Pinus armandii and Q. aliena forest
1300~1396 1~15 80 4.5±0.1 4.5±0.5 65.3±4.2
1􀆰 2􀆰 3高通量测序  依据 PCR 产物浓度,将所有扩
增成功的 PCR产物等量混合,再经定量等质量控制
后,建 Illumina 测序文库,采用 Illumina Miseq 平台
进行双末端 250 bp测序(委托诺禾致源生物信息科
技有限公司完成).
1􀆰 2􀆰 4 生物信息学分析   序列数据经 BIPES(bar⁃
coded Illumina paired⁃end sequencing) [22]流程进行初
步处理,去除低质量数据;反向互补组装 V4 序列;
根据 barcode将序列分拆到对应样品,去除 barcode;
通过对标准品计算测序错误率初步评估测序精确
率,评估合格后去除一个或多个模糊不清的序列
(Ns),去除引物上含有任何错误的、发生一个以上
碱基错配的及长度不及匹配序列 80%的序列.采用
UCHIME程序[23]筛除嵌合体序列后用 QIIME 流程
做聚类及 α、β多样性分析.其中,用 UCLUST方法进
行操作分类单元( operational taxonomic unit, OTU)
聚类,聚类采用的序列相似度为 97%.对每个 OTU,
提取其中频数最高的序列为代表序列,用 RDP [24]对
其进行系统分类.
多样性分析方法:α多样性采用 Chao1 和 Shan⁃
non等多样性指数来计算.β多样性是基于两两样品
间距离矩阵来表示样本间的差异,距离矩阵中的数
据通过 PCoA 来进行可视化展示.选用 UniFrac 和
Jaccard两种方法来计算不同林龄和林分类型下土
壤细菌群落结构差异[25-26] .
微生物群落分布影响因子分析:将样品 OTUs
标准化到 1000个后生成微生物物种矩阵,结合基本
理化性质进行影响因子的相关性分析.选用 R3.0.2
的 Vegan软件包,首先对物种数据矩阵作除趋势对
应分析(DCA)分析,根据梯度值确定相应的典范对
应分析(CCA)模型(>3 适合进行 CCA 分析[27-28] ).
经蒙特卡罗置换检验实现(Monte Carlo permutation
test)分析影响的显著性[28] .
2  结果与分析
2􀆰 1  宝天曼森林土壤微生物群落结构特点
将宝天曼 18 个土壤样品中的质控后序列按
97%相似性进行聚类,共得到 46580 个 OTUs,分属
于 60 个门,1209 个属.对宝天曼土壤全样品菌群门
水平的分布展示如图 1,其中相对丰度最高的门分
别为变形菌门 ( Proteobacteria, 29. 0%)、酸杆菌门
(Acidobacteria,18.5%)、疣微菌门(Verrucomicrobia,
10.0%)、拟杆菌门(Bacteroidetes,8.5%)、放线菌门
图 1  宝天曼森林土壤微生物门水平分布
Fig.1  Phylum distribution of Baotianman forest soil microbial
communities.
丰度最大的 10个门按相对丰度分布分别标注为 1 ~ 10 1-10 repre⁃
sented the relative abundance of top ten phyla: 1) 变形菌 Proteobacte⁃
ria; 2) 酸杆菌 Acidobacteria; 3) 疣微菌 Verrucomicrobia; 4) 拟杆菌
Bacteroidetes; 5) 未知菌 Unknown; 6) 放线菌 Actinobacteria; 7) 浮
霉菌 Planctomycetes; 8) 泉古菌 Crenarchaeota; 9) 厚壁菌 Firmicutes;
10) 绿弯菌 Chloroflexi.
57228期                  白晓旭等: 河南宝天曼不同林龄与林型森林土壤的细菌群落结构与多样性       
(Actinobacteria,6. 3%)、浮霉菌门 ( Planctomycetes,
5􀆰 9%)等,其中有 25 个低丰度门类 (相对丰度
<0.1%).从属水平上看,优势属主要有疣微菌门的
DA101(6.3%)、酸杆菌门的 Acidobacteria⁃2(5.9%)
和 Candidatus Solibacter(2.9%)、泉古菌门的 Candi⁃
datus Nitrososphaera(2.6%)等,此外还有十分丰富的
未分类种属(Other,6. 9%)和低丰度属(相对丰度
<0.1%的属,14.0%).
2􀆰 2  不同林龄锐齿栎林土壤微生物群落结构特征
与多样性
锐齿栎幼龄林共检测到土壤微生物 53个门,中
龄林 44个门,老龄林 48个门.不同林龄间土壤菌群
最高丰度的 10个门类相同,都为变形菌、酸杆菌、疣
微菌、拟杆菌、放线菌、浮霉菌、厚壁菌、绿弯菌、泉古
菌和未知菌等.由门水平菌群堆积图可以看出(图
2),疣微菌在幼龄林中分布较少,而放线菌和泉古
菌比较多.从属水平看,主要属丰度和排序也有明显
不同,幼龄林中最高丰度属为疣微菌门的 DA101,而
中龄林和老龄林中酸杆菌门的 Acidobacteria⁃2纲的
未知属丰度最高.
不同林龄土壤样品共得到 46627 个 OTUs,根据
林龄的不同对土壤样品序列分组后进行 ANOVA
分析,统计得到 157 个丰度差异显著的 OTUs ( P
<0􀆰 05).对差异显著的 OTUs 进行分析发现,不同林
龄锐齿栎林土壤菌群中差异性 OTUs 主要属于 20
个门类,101个属,其中有 50个差异 OTUs出现在变
图 2  不同林龄锐齿栎土壤菌群主要门水平分布
Fig.2  Phylum distribution of soil microbial communities from
Quercus aliena forests of different ages.
QAF40:锐齿栎幼龄林 40⁃year⁃old Quercus aliena forest; QAF80:锐齿栎
中龄林 80⁃year⁃old Quercus aliena forest; QAF160: 锐齿栎老龄林 160⁃
year⁃old Quercus aliena forest. 下同 The same below. 丰度最大的 10个门
按相对丰度分布分别标注为 1~ 10 1-10 represented the relative abun⁃
dance of top ten phyla: 1) 变形菌 Proteobacteria; 2) 酸杆菌 Acidobacte⁃
ria; 3) 放线菌 Actinobacteria; 4) 疣微菌 Verrucomicrobia; 5) 拟杆菌
Bacteroidetes; 6) 未知菌 Unknown; 7) 浮霉菌 Planctomycetes; 8) 泉古
菌 Crenarchaeota; 9) 厚壁菌 Firmicutes; 10) 绿弯菌 Chloroflexi.
形菌门,最高丰度的差异 OTU 注释到红游动菌属
(Rhodoplanes),主要富集在中龄林土壤菌群中(幼
龄林样本中得到 128个,中龄林样本中得到 260 个,
老龄林样品中得到 122 个);次高丰度 OTU 注释到
疣微菌门的 DA101(幼龄林样本中得到 130个,中龄
林样本中得到 147个,老龄林样品中得到 122个).
    将每个样品的 reads标准化到 10000条后,通过
3种微生物多样性指数比较不同林龄土壤的 α 多样
性(图 3),对 3 种多样性指数结果分别做 ANOVA
分析可知,不同林龄锐齿栎林土壤菌群多样性相差不
大(P>0.05).从 Shannon 指数分析来看,锐齿栎老龄
林(11.25)和幼龄林(11.21)土壤细菌多样性较高,中
龄林最低(10.83),Chao1指数结果显示锐齿栎林幼龄
林物种丰富度最高,老龄林居中,中龄林最低.
选用 UniFrac 和 Jaccard 两种指数进行不同林
龄土壤细菌群落结构相似性分析.采用 Jaccard 指数
表征菌群相似度时,一般认为 0.5 ~ 0.75 为中等相
似[26,29],3种林龄土壤细菌群落间 Jaccard 相似性系
数为 0.557~0.750,说明 40、80和 160年林龄土壤菌
群结构处于中等相似水平.而采用 UniFrac距离法表
征进化的差异,一般认为 UniFrac值在 0.5时群落结
构中等相似,UniFrac值越接近 1表明群落结构差异
越大[30],这 3个林龄土壤细菌群落间 UniFrac 距离
在 0.686~0.803,说明尽管 Jaccard 指数显示样品间
中等相似,但从进化的角度上看,不同林龄间却存在
着明显的群落结构差异.从图4可以看出,3组不同
图 3  不同林龄锐齿栎土壤微生物 α多样性
Fig.3  α⁃diversity of soil microbes from Quercus aliena forests
of different ages.
6722 应  用  生  态  学  报                                      26卷
图 4  不同林龄锐齿栎土壤微生物 β多样性
Fig.4  β⁃diversity of soil microbes from Quercus aliena forests
of different ages.
-1、-2、-3表示重复样品-1, -2 and -3 indicated repeated samples.
下同 The same below.
林龄土样按林龄不同区分开来,样品聚类树图显示,
幼龄林和老龄林土壤菌群聚在一起,与中龄林的差
异较大.
2􀆰 3  不同林分类型土壤微生物群落结构特征及多
样性
在林龄同为 80年的 4个林地样品中,锐齿栎林
土壤样品中共检测到 44 个门,栓皮栎林中共含 46
个门,短柄枹林中有 49个门,锐齿栎 /华山松混交林
有 48个门类.不同林型间土壤菌群最高丰度的 10
个门类为变形菌、酸杆菌、疣微菌、拟杆菌等,与前述
不同林龄锐齿栎林 10大主要门类相同.由不同林型
间门水平菌群堆积图(图 5)可以看出,锐齿栎和锐
齿栎 /华山松混交林主要门类相对丰度分布较为相
似,而短柄枹林则与其他林型相对丰度分布差异最
大,主要表现在短柄枹林的疣微菌最少而泉古菌含
量最大.从其他分类层级上看,4 类森林土壤细菌分
布也明显不同,如在属水平上,锐齿栎林和锐齿栎 /
华山松混交林中最高丰度属都为疣微菌门的
DA101,而栓皮栎林和短柄枹林中都显示未分类属
含量最高,短柄枹林中 Candidatus Nitrososphaera 相
对丰度远高于其他 3种林型,且在锐齿栎和锐齿栎 /
华山松混交林中分布最多的 DA101 在短柄枹林中
丰度明显降低.
不同林分类型土壤样品菌群序列统计共得到
60597个 OTUs,根据林型的不同对土壤样品序列进
行分组后进行ANOVA分析,发现其中229个OTUs
图 5  宝天曼不同林分类型土壤菌群主要门水平分布
Fig.5  Phylum distribution of soil microbial communities from
different type forests.
QGB: 短柄枹林 Quercus glandulifera forest; QV: 栓皮栎林 Q. variabi⁃
lis forest; PA: 锐齿栎 /华山松混交林 Mixed P. armandii and Q. aliena
forest. 下同 The same below. 丰度最大的 10 个门按相对丰度分布分
别标注为 1~10 1-10 represented the relative abundance of top ten phy⁃
la: 1) 变形菌 Proteobacteria; 2) 酸杆菌 Acidobacteria; 3) 疣微菌
Verrucomicrobia; 4) 拟杆菌 Bacteroidetes; 5) 未知菌 Unknown; 6) 浮
霉菌 Planctomycetes; 7) 放线菌 Actinobacteria; 8) 泉古菌 Crenarcha⁃
eota; 9) 厚壁菌 Firmicutes; 10) 绿弯菌 Chloroflexi.
在各样品间丰度差异显著(P<0.05),主要属于 19
个门类,127个属,其中有 50个差异 OTUs属于变形
菌门.不同林型间的显著差异 OTUs 中,锐齿栎林和
锐齿栎 /华山松混交林最高丰度 OTU 注释为红游动
菌属(Rhodoplanes),锐齿栎林样本中得到 260个,锐
齿栎 /华山松混交林样本中得到 255个,栓皮栎林样
本中得到 84个,短柄枹林样本中得到 91个;栓皮栎
林中最高丰度 OTUs注释为肠杆菌科(Enterobacteri⁃
aceae),锐齿栎林检测到63个,锐齿栎 / 华山松混
图 6  不同林型土壤微生物 α多样性
Fig.6  α⁃diversity of soil microbes from different type forests.
77228期                  白晓旭等: 河南宝天曼不同林龄与林型森林土壤的细菌群落结构与多样性       
图 7  不同林型土壤微生物 β多样性
Fig.7  β⁃diversity of soil microbes from different type forests.
交林147 个,栓皮栎 196 个,短柄枹 78 个;而短柄
枹林样本中最高丰度 OTUs 则注释为酸杆菌科
(Acidobacteriaceae),短柄枹林检测到 116 个,锐齿
栎林 14个,锐齿栎 /华山松混交林 30 个,栓皮栎林
107个.
将每个样品的 reads标准化到 10000条后,计算
不同样品的 α多样性,经 ANOVA分析可知,不同林
分类型锐齿栎土壤菌群多样性相差不大(P>0.05),
但变化趋势相近.从 Shannon 指数看,4 个样地土壤
细菌的多样性指数由大到小排序依次为:锐齿栎 /华
山松混交林(11􀆰 47)、栓皮栎林(11.24)和短柄枹林
(11􀆰 15)、锐齿栎(10.83),Chao1 指数显示短柄枹林
丰富度略高于栓皮栎林,但整体趋势上看,仍以锐齿
栎 /华山松混交林最高,锐齿栎林最低.
    不同林型土壤细菌群落间 Jaccard 相似性系数
为 0.538 ~ 0. 769,说明不同林分类型间土壤菌群
结构存在一定的相似度.而 Unifrac 距离在 0.697 ~
0􀆰 810,说明从进化的角度上看确实存在明显的群落
结构差异.从图 7 可知,不同林型土样被区分开来,
从样品聚类树图可见,锐齿栎 /华山松混交林和锐齿
栎林土壤样品群落结构比较接近,其次是栓皮栎林,
短柄枹距离最远,不同林龄间群落结构差异显著
(P<0.01).
2􀆰 4  不同林龄和不同林分类型下土壤菌群影响因
子典范对应分析
将样品 OTUs 标准化到 1000 条后,结合土壤
pH、全氮、有机碳 3种基本理化性质,首先用 R软件
图 8  环境因子典范对应分析
Fig.8  Canonical correspondence analysis of environmental fac⁃
tors.
SOC: 土壤有机碳含量 Content of soil organic carbon; TN: 土壤全氮
含量 Content of soil total nitrogen.
做 DCA分析来判别选择线性排序还是单峰排序,本
文中 DCA最大排序轴值为 4.10,因此选用单峰模型
CCA进行后续分析和讨论.依据 CCA 排序结果(图
8), 可将全部样品划分为较明显的两组,其中不同
林龄锐齿栎林和锐齿栎 /华山松混交林聚成一组,而
短柄枹和栓皮栎林聚成一组.由 CCA 图中箭头的连
线长度可以明显看出,选取的几种土壤理化性质
对菌群分布都有一定程度的影响,3 种环境因子
的合并解释量为 17. 2%.经蒙特卡罗秩和检验实
现(Monte Carlo permutation test)分析,可知 3 种环
境因子都对土壤细菌物种分布具有显著影响(P =
0􀆰 001),其中土壤有机碳解释了土壤菌群分布 9.8%
的变异,土壤全氮解释了菌群分布 9.7%的变异,土
壤 pH解释了 8.8%的变异.
3  讨    论
宝天曼地区森林土壤细菌主要以变形菌
(29􀆰 0%)、酸杆菌(18.5%)、疣微菌(10.0%)3 大门
类丰度最高,拟杆菌、放线菌、浮霉菌等也比较丰富.
据 Janssen[31]对森林、草地、牧场、高山草甸等不同
利用方式的土壤细菌分布研究显示,土壤中变形菌
门平均相对丰度为 39%(10% ~ 77%),酸杆菌门平
均相对丰度为 20%(5% ~46%)、疣微菌门平均丰度
为 7%(0~ 21%),3 大类主要土壤微生物相对丰度
依土壤类型不同而在一定范围内波动.Shen 等[32]采
用 454焦磷酸测序法对温带森林长白山土壤菌群结
构研究发现,长白山土壤主要微生物门类包括
变形菌(49.6%)、酸杆菌(18.7%)和放线菌(12.8%)
等,而疣微菌只占 0.17%.这与本研究中宝天曼的
土壤差别很大.其中,变形菌的差异主要来自于
8722 应  用  生  态  学  报                                      26卷
β⁃proteobacteria的丰度变化(长白山 18.6%,而宝天
曼仅为 5.6%),β⁃变形菌包括很多好氧或兼性细菌,
通常具有可变的降解能力,也有一些无机化能种类
如可以氧化氨的亚硝化单胞菌属(Nitrosomonas)和
光合种类如红环菌属(Rhodocyclus)等,丰度主要受
pH值和营养等环境因子影响[33],长白山森林土壤
平均 pH 在 5 以上[27],高于宝天曼土壤( pH 2.4 ~
4􀆰 5).疣微菌广泛分布于根际和团聚体土壤中[34],
其丰度与多种碳水化合物代谢相关的基因相关,可
能在碳生物地球化学循环中发挥着重要作用,同时
地理格局也是其丰度改变的主要驱动因子[35] .因此
这种变形菌和疣微菌的丰度差异可能与地理因素引
起的生态多样性和土壤营养情况不同有关.
Chao1和 Shannon 指数的计算公式不同,因而
在表征 α 多样性时,所反映的侧重点略有不同:
Chao1指数在观测 OTUs 基础上加入了低丰度序列
(只有一条和只有两条的 OTU数目)的计算,因此在
生态学中常用来估计物种总数;而 Shannon 指数同
时受到物种丰富度和均度两者的影响[23-24] .由于本
研究中不同样品的 Chao1指数在所有样品间的趋势
基本一致,而加入了均度影响的 Shannon 指数在所
有样品间的趋势略有不同,主要表现在菌群分布比
较接近的群落如幼龄林和老龄林间、栓皮栎和短柄
枹林间,多样性分析结果存在不一致性,但在菌群结
构差异较大的样品如中龄林和老龄林间、锐齿栎林
和锐齿栎 /华山松混交林间则无显著差别.Haegeman
等[36]指出,Chao1 指数往往会引起多样性的低估,
且这种低估会随着样本数和菌群物种数的增加而逐
渐加大,这可能是本研究中采用不同多样性指数得
到不一致结果的主要原因.
不同林龄锐齿栎林土壤样品各自存在丰度显著
差异门类和特有低丰度和高丰度类群,一些类群随
林龄增加也存在一些趋势性变化,如拟杆菌门和蓝
藻门随林龄增加相对丰度逐渐升高,而厚壁菌则相
反.不同林龄土壤中存在差异显著 OTUs,其中属于
变形菌门的 OTUs最多.变形菌是土壤细菌中最大的
门类,包括很多病原菌,也有自由生活的种类,包括
很多可以进行固氮的细菌(如根瘤菌) [37] .差异显著
OTUs中丰度最高的是红游动菌,红游动菌是兼性厌
氧光养菌,能利用各种有机化合物进行光合异养生
长[38],在营养丰富的中龄林土壤菌群中含量最高.
多样性指数分析显示中龄林多样性最低,Chao1 指
数呈现幼龄林略高于老龄林,而 Shannon 指数则是
老龄林高于幼龄林,其原因可能是幼龄林物种丰富
图 9  不同来源样品聚类树图
Fig. 9   Phylogenetic tree of all soil samples from different
sources.
度略高于老龄林,老龄林土壤菌群逐渐稳定,其土壤
菌群均度较好.结合 UniFrac的群落差异性分析结果
可知,中龄林和其他两种林型土壤菌群存在较为明
显的聚类差异.这种差异可能是由于:幼龄林植被定
植时间较短,生物量相对较低,养分积累较少;而随
林龄增加植物生长加速,植被群落逐渐完整,林分密
度和林下植被盖度增大,凋落物和根系分泌物增多
和植物间种间竞争更为激烈,因此对应的林下土壤
菌群种类和丰度也与幼龄林和老龄林差异显
著[39-40] .老龄林植物生长减慢,植被群落趋于稳定,
逐年累积的枯枝落叶和植物残体等凋落物分解和腐
殖质化,一些降解菌如拟杆菌等的相对丰度明显
增加.
4种不同林分类型菌群分布明显不同,存在很
多差异性类群,但锐齿栎和锐齿栎 /华山松混交林的
菌群分布较为相似,从聚类角度看,锐齿栎 /华山松
混交林和锐齿栎林群落结构比较接近,而与短柄枹
和栓皮栎林显著不同,可能是由于两者间有相交叉
的优势植被树种(锐齿栎),因此有相似的凋落物组
成,导致表层土壤微生物存在相似种群.但锐齿栎林
土壤菌群丰富度和均一度最低,而锐齿栎 /华山松混
交林的多样性指数最高,说明混交林能提高土壤微
生物的多样性,这与已有的研究一致[41] .短柄枹林
中氨氧化古细菌 Candidatus Nitrososphaera 相对丰度
远高于其他 3种林型.
Shen等[32]指出,土壤 pH、全氮、有机碳都是影
响土壤菌群结构的重要因子.分析发现,随 pH、土壤
全氮和有机碳的增加,土壤菌群物种多样性都呈增
加趋势.而不同林龄间 pH 值十分接近,全氮变化与
多样性同步,但随土壤有机碳增加丰富度指数和
Chao1指数逐渐减小.表明 pH、土壤全氮和有机碳都
97228期                  白晓旭等: 河南宝天曼不同林龄与林型森林土壤的细菌群落结构与多样性       
是与森林土壤菌群变化密切相关的影响因子.结合 3
种环境因子与微生物物种分布的典范对应分析,发
现 3种因子对菌群分布影响极为显著,分别解释了
土壤菌群 8.8%~9.8%的变异. 森林生态系统极其复
杂,其林下土壤菌群分布受地形、植被、土壤营养与
健康、人类干扰等多种因素影响[42],因此还存在很
多重要的环境影响因子,需要进一步分析.
    将不同林龄和不同林型的样品群落结构聚类结
果进行整合分析发现(图 9),锐齿栎林和其他林型
很好地分开,锐齿栎 /华山松混交林聚类关系处于锐
齿栎林和其他两种林型之间,但不同林龄的锐齿栎
则聚在一起.结合典范对应分析图可知,不同锐齿栎
林在排序图中聚在一起,而另几种不同林分类型则
远远分开.表明在 40 ~ 160 年林龄尺度上,林分类型
变化引起的菌群结构变化大于林龄,可能与地表凋
落物的组成成分、数量及土层中植物根系分布状况、
数量不同等有关[43] .
森林土壤中微生物的空间异质性高[44],因此同
一标准地中的 3个重复样品其微生物群落有时也可
看到明显差异[45],而宝天曼地处由南至北的过渡
带,环境异质性强,这导致在 PCoA分析中散点分布
较为分散.
参考文献
[1]  Gobat JM, Aragno M, Matthey W. The Living Soil:
Fundamentals of Soil Science and Soil Biology. New
York: Science Publishers, 2004
[2]  Luan J, Liu S, Wang J, et al. Rhizospheric and hetero⁃
trophic respiration of a warm⁃temperate oak chronose⁃
quence in China. Soil Biology and Biochemistry, 2011,
43: 503-512
[3]  Banning NC, Gleeson DB, Grigg AH, et al. Soil micro⁃
bial community successional patterns during forest eco⁃
system restoration. Applied and Environmental Microbio⁃
logy, 2011, 77: 6158-6164
[4]  Harris J. Soil microbial communities and restoration
ecology: Facilitators or followers? Science, 2009, 325:
573-574
[5]  Preem JK, Truu J, Truu M, et al. Bacterial community
structure and its relationship to soil physico⁃chemical
characteristics in alder stands with different management
histories. Ecological Engineering, 2012, 49: 10-17
[6]  Kohler I, Leehan P, Yebernetsky K. Assessing soil bac⁃
terial diversity in southwestern Pennsylvania. Procee⁃
dings of the National Conference on Undergraduate
Research. Ithaca, NY: Ithaca College, 2011
[7]  Lin YT, Huang YJ, Tang SL, et al. Bacterial communi⁃
ty diversity in undisturbed perhumid montane forest soils
in Taiwan. Microbial Ecology, 2010, 59: 369-378
[8]   Liu ZF, Fu BJ, Zheng XX, et al. Plant biomass, soil
water content and soil N:P ratio regulating soil microbial
functional diversity in a temperate steppe: A regional
scale study. Soil Biology and Biochemistry, 2010, 42:
445-450
[9]  Harris J. Measurements of the soil microbial community
for estimating the success of restoration. European Jour⁃
nal of Soil Science, 2003, 54: 801-808
[10]   Luan JW, Liu SR, Wang JX, et al. Factors affecting
spatial variation of annual apparent Q10 of soil respiration
in two warm temperate forests. PLoS One, 2013, 8:
e64167
[11]  Liu YC, Liu SR, Wang JX, et al. Variation in soil res⁃
piration under the tree canopy in a temperate mixed
forest, central China, under different soil water condi⁃
tions. Ecological Research, 2014, 29: 133-142
[12]  Lai JS, Coomes DA, Du XJ, et al. A general combined
model to describe tree⁃diameter distributions within sub⁃
tropical and temperate forest communities. Oikos, 2013,
122: 1636-1642
[13]  You YM, Wang J, Huang XM, et al. Relating microbial
community structure to functioning in forest soil organic
carbon transformation and turnover. Ecology and Evolu⁃
tion, 2014, 4: 633-647
[14]  Baldrian P, Kola rˇík M, Štursová M, et al. Active and
total microbial communities in forest soil are largely
different and highly stratified during decomposition. The
ISME Journal, 2012, 6: 248-258
[15]  Lin YT, Jangid K, Whitman WB, et al. Soil bacterial
communities in native and regenerated perhumid mon⁃
tane forests. Applied Soil Ecology, 2011, 47: 111-118
[16]  Van Elsas JD, Boersma FGH. A review of molecular
methods to study the microbiota of soil and the myco⁃
sphere. European Journal of Soil Biology, 2011, 47:
77-87
[17]  Quail MA, Smith M, Coupland P, et al. A tale of three
next generation sequencing platforms: Comparison of Ion
Torrent, Pacific Biosciences and Illumina MiSeq se⁃
quencers. BMC Genomics, 2012, 13: 341, doi:10.1186 /
1471⁃2164⁃13⁃341
[18]  Logares R, Haverkamp THA, Kumar S, et al. Environ⁃
mental microbiology through the lens of high⁃throughput
DNA sequencing: Synopsis of current platforms and
bioinformatics approaches. Journal of Microbiological
Methods, 2012, 91: 106-113
[19]   Yuan Z⁃L (袁志良), Chen Y (陈   云), Wei B⁃L
(韦博良), et al. Species habitat correlation analysis in
temperate⁃subtropical ecological transition zone. Acta
Ecologica Sinica (生态学报), 2013, 33(24): 7819-
7826 (in Chinese)
[20]  Yan D⁃F (闫东锋), Zhu Y (朱  滢), Yang X⁃T (杨
喜田). Species diversity and stability of a natural Quer⁃
cus forest in Baotianman, Henan Province. Journal of
Zhejiang A&F University (浙江农林大学学报), 2011,
28(4): 628-633 (in Chinese)
[21]  Caporaso JG, Lauber CL, Walters WA, et al. Global
patterns of 16S rRNA diversity at a depth of millions
of sequences per sample. Proceedings of the National
0822 应  用  生  态  学  报                                      26卷
Academy of Sciences of the United States of America,
2011, 108: 4516-4522
[22]  Zhou HW, Li DF, Tam NF, et al. BIPES, a cost⁃
effective high⁃throughput method for assessing microbial
diversity. The ISME Journal, 2011, 5: 741-749
[23]  Edgar RC, Haas BJ, Clemente JC, et al. UCHIME im⁃
proves sensitivity and speed of chimera detection. Bioin⁃
formatics, 2011, 27: 2194-2200
[24]   Maidak BL, Cole JR, Lilburn TG, et al. The RDP⁃Ⅱ
(ribosomal database project) . Nucleic Acids Research,
2001, 29: 173-174
[25]  Caporaso JG, Kuczynski J, Stombaugh J, et al. QIIME
allows analysis of high⁃throughput community sequencing
data. Nature Methods, 2010, 7: 335-336
[26]   Qian Y⁃Q (钱迎倩), Ma K⁃P (马克平). Principles
and Methodologies of Biodiversity Studies. Beijing: Chi⁃
nese Science and Technology Press, 1994 (in Chinese)
[27]  Lepš J, Milauer P. Multivariate Analysis of Ecological
Data Using CANOCO. Cambridge, UK: Cambridge
University Press, 2003
[28]  Lai J⁃S (赖江山), Mi X⁃C (米湘成). Ordination
Analysis of Ecological Data Using Vegan Package in R:
Advances in Biodiversity Conservation and Research in
China Ⅸ. Beijing: China Meteorological Press, 2011
(in Chinese)
[29]  Su Y⁃C (苏永春), Gou Y⁃B (勾影波), Yu D (郁
达), et al. Diversity of soil invertebrate communities at
Yushan Hill, Changshu, Jiangsu Province. Biodiversity
Science (生物多样性), 2004, 12(3): 333-338 ( in
Chinese)
[30]   Lozupone C, Knight R. UniFrac: A new phylogenetic
method for comparing microbial communities. Applied
and Environmental Microbiology, 2005, 71: 8228-8235
[31]  Janssen PH. Identifying the dominant soil bacterial taxa
in libraries of 16S rRNA and 16S rRNA genes. Applied
and Environmental Microbiology, 2006, 72: 1719-1728
[32]  Shen CC, Xiong JB, Zhang HY, et al. Soil pH drives
the spatial distribution of bacterial communities along
elevation on Changbai Mountain. Soil Biology and Bio⁃
chemistry, 2013, 57: 204-211
[33]  Nacke H, Thürmer A, Wollherr A, et al. Pyrosequen⁃
cing⁃based assessment of bacterial community structure
along different management types in German forest and
grassland soils. PLoS One, 2011, 6: e17000
[34]  Kielak A, Pijl AS, Van Veen JA, et al. Differences in
vegetation composition and plant species identity lead to
only minor changes in soil⁃borne microbial communities
in a former arable field. FEMS Microbial Ecology,
2008, 63: 372-382
[35]  Fierer N, Ladau J, Clemente JC, et al. Reconstructing
the microbial diversity and function of pre⁃agricultural
tallgrass prairie soils in the United States. Science,
2013, 342: 621-624
[36]  Haegeman B, Hamelin J, Moriarty J, et al. Robust esti⁃
mation of microbial diversity in theory and in practice.
The ISME Journal, 2013, 7: 1092-1101
[37]  Spain AM, Krumholz LR, Elshahed MS. Abundance,
composition, diversity and novelty of soil Proteobacteria.
The ISME Journal, 2009, 3: 992-1000
[38]  Sun Y (孙  勇), Xie S⁃T (谢数涛), Yang Z⁃M (杨
泽民), et al. Isolation and identification of a strain of
photosynthetic bacteria with high⁃yield production of
carotenoid pigments. Journal of Jinan University (Natu⁃
ral Science) (暨南大学学报: 自然科学版), 2006,
27(1): 146-150 (in Chinese)
[39]  Jiao R⁃Z (焦如珍), Yang C⁃D (杨承栋), Sun Q⁃W
(孙启武), et al. Changes in soil microbial amount and
biomass during the development of Chinese fir planta⁃
tion. Scientia Silvae Sinicae (林业科学), 2005, 41
(6): 163-165 (in Chinese)
[40]  Mao Z⁃J (毛钟警), Wang G⁃Y (王冠玉), Zhou Y⁃H
(周远和), et al. Progress on the distribution regularity
of soil microorganism in forestry and the relationship with
soil enzymes. Guangdong Agricultural Sciences (广东农
业科学), 2012, 39(12): 60-62 (in Chinese)
[41]  Qi Y⁃J (齐雅静), Bai S⁃L (白淑兰), Han S⁃L (韩胜
利), et al. Study on the microbial diversity of rhizo⁃
sphere soil in different forest types of Pinus tabulaeformis
and Ostryopsis davidaiana. Chinese Agricultural Science
Bulletin (中国农学通报), 2011, 27(2): 53-56 ( in
Chinese)
[42]  Zheng H (郑  华), Ouyang Z⁃Y (欧阳志云), Wang
X⁃K (王效科), et al. Effects of forest restoration pat⁃
terns on soil microbial communities. Chinese Journal of
Applies Ecology (应用生态学报), 2005, 15 ( 11):
2019-2024 (in Chinese)
[43]   Bi J⁃T (毕江涛), He B⁃H (贺达汉). Research ad⁃
vances in effects of plant on soil microbial diversity. Chi⁃
nese Agricultural Science Bulletin (中国农学通报),
2009, 25(9): 244-250 (in Chinese)
[44]  Han Y⁃Z (韩有志), Wang Z⁃Q (王政权). Spatial
heterogeneity and forest regeneration. Chinese Journal of
Applies Ecology (应用生态学报), 2002, 13 ( 5):
615-619 (in Chinese)
[45]  Lou X (娄   鑫), Cui X⁃Y (崔晓阳). The study of
soil microbial diversity in temperate forest secondary
succession process. Chinese Journal of Soil Science (土
壤通报), 2011, 42(3): 557-561 (in Chinese)
作者简介  白晓旭,女,1988年生,硕士. 主要从事环境微生
物学和分子生态学研究. E⁃mail: baixiaoxu11@ mails. gucas.
ac.cn
责任编辑  肖  红
18228期                  白晓旭等: 河南宝天曼不同林龄与林型森林土壤的细菌群落结构与多样性