全 文 :书犇犗犐:10.11686/犮狔狓犫20150115 犺狋狋狆://犮狔狓犫.犾狕狌.犲犱狌.犮狀
李刚,修伟明,王杰,吴元凤,赵建宁,宋晓龙,杨殿林.不同植被恢复模式下呼伦贝尔沙地土壤反硝化细菌狀犻狉犓基因组成结构和多样性研究.草
业学报,2015,24(1):115123.
LiG,XiuW M,WangJ,WuYF,ZhaoJN,SongXL,YangDL.Communitystructureanddiversityofsoildenitrifyingbacteriaofthe狀犻狉犓gene
typeunderdifferentvegetationrestorationpatternsintheHulunbeierSandyLand,InnerMongolia.ActaPrataculturaeSinica,2015,24(1):115123.
不同植被恢复模式下呼伦贝尔沙地土壤反硝化
细菌狀犻狉犓基因组成结构和多样性研究
李刚,修伟明,王杰,吴元凤,赵建宁,宋晓龙,杨殿林
(农业部产地环境质量重点实验室,农业部环境保护科研监测所,天津300191)
摘要:本研究选取呼伦贝尔沙地5种植被恢复模式(单播冰草、单播羊柴、单播柠条、冰草+柠条和冰草+柠条+羊
柴+披碱草)为研究对象,以完全退化的沙地为对照,采用分子生物学技术分析了土壤反硝化细菌狀犻狉犓 基因组成
结构和多样性,并探讨了与土壤理化因子的相关性。结果表明,土壤反硝化细菌狀犻狉犓 基因组成存在明显差异。单
播羊柴模式土壤反硝化细菌狀犻狉犓基因多样性指数最高,其次为单播冰草和单播柠条,而两种混播模式最低,这5
种植被恢复模式均显著高于对照裸地。序列和系统发育分析结果表明,所获得的狀犻狉犓基因片段与NCBI数据库中
的狀犻狉犓基因的相似度在81%~99%之间,分别归属于苍白杆菌属、根瘤菌属及不可培养的细菌,其中与不可培养
细菌相关的片段占到总数的50%以上。典范对应和相关分析结果表明,速效磷、pH、全氮、全磷和含水量对反硝化
细菌的影响均达到显著水平,且5种理化因子间均具有极显著相关性,说明这5种理化因子相互关联,共同对反硝
化细菌狀犻狉犓基因组成变化起到决定作用。
关键词:呼伦贝尔沙地;植被恢复;反硝化细菌;狀犻狉犓基因;PCR-DGGE
犆狅犿犿狌狀犻狋狔狊狋狉狌犮狋狌狉犲犪狀犱犱犻狏犲狉狊犻狋狔狅犳狊狅犻犾犱犲狀犻狋狉犻犳狔犻狀犵犫犪犮狋犲狉犻犪狅犳狋犺犲狀犻狉犓犵犲狀犲狋狔狆犲
狌狀犱犲狉犱犻犳犳犲狉犲狀狋狏犲犵犲狋犪狋犻狅狀狉犲狊狋狅狉犪狋犻狅狀狆犪狋狋犲狉狀狊犻狀狋犺犲犎狌犾狌狀犫犲犻犲狉犛犪狀犱狔犔犪狀犱,犐狀狀犲狉
犕狅狀犵狅犾犻犪
LIGang,XIU Weiming,WANGJie,WUYuanfeng,ZHAOJianning,SONGXiaolong,YANGDianlin
犓犲狔犔犪犫狅狉犪狋狅狉狔狅犳犗狉犻犵犻狀犪犾犃犵狉狅犲狀狏犻狉狅狀犿犲狀狋犙狌犪犾犻狋狔狅犳犕犻狀犻狊狋狉狔狅犳犃犵狉犻犮狌犾狋狌狉犲,犃犵狉狅犈狀狏犻狉狅狀犿犲狀狋犪犾犘狉狅狋犲犮狋犻狅狀犐狀狊狋犻狋狌狋犲,
犕犻狀犻狊狋狉狔狅犳犃犵狉犻犮狌犾狋狌狉犲,犜犻犪狀犼犻狀300191,犆犺犻狀犪
犃犫狊狋狉犪犮狋:Weselectedfivevegetationrestorationmethodsincludingmonoplantingof犆犪狉犪犵犪狀犪犽狅狉狊犺犻狀狊犽犻犻
(UC),犃犵狉狅狆狔狉狅狀犮狉犻狊狋犪狋狌犿 (UA),and犎犲犱狔狊犪狉狌犿犳狉狌狋犻犮狅狊狌犿 (UH),mixedplantingof犃犵狉狅狆狔狉狅狀犮狉犻狊狋犪
狋狌犿and犎犲犱狔狊犪狉狌犿犳狉狌狋犻犮狅狊狌犿 (AC)andof犃犵狉狅狆狔狉狅狀犮狉犻狊狋犪狋狌犿,犎犲犱狔狊犪狉狌犿犳狉狌狋犻犮狅狊狌犿,犆犪狉犪犵犪狀犪犽狅狉
狊犺犻狀狊犽犻犻,and犈犾狔犿狌狊狀狌狋犪狀狊(ACHE)toanalyzethecommunitystructureanddiversityof狀犻狉犓typesoildeni
trifyingbacteriaandtheirrelationshipswiththesoilphysicalandchemicalfactors.Barelandwasusedasacon
trol(CK).Therearesignificantdifferencesinthecompositionofthesoil狀犻狉犓typedenitrifyingbacteria.The
Shannonindexofthe狀犻狉犓geneisthehighestunderUH,folowedbyunderUA,UC,ACandACHE,andthe
lowestinCKwithsignificantdifferencesbetweenthefivemethods.Thesequenceandphylogeneticanalysis
第24卷 第1期
Vol.24,No.1
草 业 学 报
ACTAPRATACULTURAESINICA
2015年1月
Jan,2015
收稿日期:20131223;改回日期:20140225
基金项目:国家自然科学基金项目(31170435,31200424)和“十二五”国家科技计划项目(2012BAD13B07)资助。
作者简介:李刚(1981),男,辽宁辽阳人,助理研究员。Email:gangli1981@gmail.com
通讯作者Correspondingauthor.Email:yangdianlin@caas.cn
showedthatthesimilarityofthe狀犻狉犓fragmentstothe狀犻狉犓genesregisteredinGenBankwasbetween81%
and99%,andthattheyweremainlyfrom犗犮犺狉狅犫犪犮狋狉狌犿,犚犺犻狕狅犫犻狌犿,andsomeotherunidentifiedbacteria.
Thesebacterialspeciesaccountedformorethan50%ofthetotalfragments.Canonicalcorrespondenceandcor
relationanalysesshowedthatthesoilavailablephosphorus,pH,totalnitrogen,totalphosphorusandsoil
moisturehadsignificanteffectsonthedenitrifyingbacterialcommunities,andtherewereverysignificantcorre
lationsamongthesefivesoilphysicalandchemicalfactors.Itissuggestedthatinteractionsbetweenthefivesoil
physicalandchemicalfactorsplayedadecisiveroleinthepopulationchangesofthe狀犻狉犓typedenitrifyingbac
teria.
犓犲狔狑狅狉犱狊:Hulunbeiersandyland;vegetationrestoration;denitrifyingmicrobe;狀犻狉犓gene;PCRDGGE
反硝化作用是由反硝化细菌通过异化作用使硝态氮最终转化为N2 的生物过程,是生态系统氮循环的重要
组成部分[1]。土壤因反硝化作用而丧失了20%~30%的氮肥[2]。亚硝酸还原酶(nir)是反硝化细菌具有的关键
酶[3],是微生物驱动的反硝化作用中的限速步骤[4]。许多研究将狀犻狉犓 和狀犻狉犛 基因作为环境样品中反硝化细菌
群落的分子标记物,而且发现狀犻狉犓 基因对环境因子的响应较狀犻狉犛基因更敏感[5]。研究人员在多种环境中均发
现狀犻狉犓 型的反硝化细菌的存在,如农田、森林、草地、沉积物以及水体等[610]。同时还发现,土壤性质的改变可以
影响反硝化细菌的群落结构和种群丰度,如全氮含量、NO3-含量、含水量、pH等[11],而反硝化细菌的群落结构
和种群丰度又能反映这些调控因子的综合作用,因此研究土壤反硝化细菌群落和对环境的响应具有重要意义。
沙漠化是制约我国草原地区发展的重要环境问题,已严重威胁到我国经济的可持续发展[12]。植被恢复是治
理沙漠化的根本解决途径。我国对沙化草地恢复研究集中在植被恢复模式设计、恢复后地表植物群落结构、土壤
理化性质、地下生物量变化及土壤动物多样性变化等的变化特征[1317],而不同植被恢复模式下沙地土壤反硝化细
菌分子生态的相关研究未见报道。因此,通过研究不同植被恢复模式下沙地土壤反硝化细菌群落结构组成变化
及其与土壤理化因子的相互关联性,对认识沙地植被恢复过程反硝化细菌的响应是非常关键的。本研究运用
PCR-DGGE及克隆测序技术,通过对呼伦贝尔沙地不同植被恢复模式下土壤反硝化细菌群落多样性和系统发
育研究,探讨反硝化细菌多样性对沙地植被恢复模式和环境因子的响应,以期为呼伦贝尔沙地植被恢复和重建提
供理论依据。
1 材料与方法
1.1 研究区概况
本试验以内蒙古呼伦贝尔市陈巴尔虎旗长期定位试验站(49°12′N,118°54′E)为研究地点,该区域草原已显
著退化及沙化,具有呼伦贝尔地区沙地的典型特征。该区域海拔为618.0m,年均气温-2.6℃,年均降水量
307.7mm,年均风速3.5m/s,属中温带半温润和半干旱大陆性气候。沙地周边以大针茅(犛狋犻狆犪犵狉犪狀犱犻狊)和羊
草(犔犲狔犿狌狊犮犺犻狀犲狀狊犻狊)为主要自然植被,土壤类型为沙质栗钙土。
1.2 试验设计与样品采集
试验始于2006年7月,分别设置单播冰草(Unicast犃犵狉狅狆狔狉狅狀犮狉犻狊狋犪狋狌犿,UA)、单播羊柴(Unicast
犎犲犱狔狊犪狉狌犿犳狉狌狋犻犮狅狊狌犿,UH)、单播柠条(Unicast犆犪狉犪犵犪狀犪犽狅狉狊犺犻狀狊犽犻犻,UC)、冰草+柠条混播(犃犵狉狅狆狔狉狅狀
犮狉犻狊狋犪狋狌犿+犆犪狉犪犵犪狀犪犽狅狉狊犺犻狀狊犽犻犻,AC)、冰草+柠条+羊柴+披碱草混播(犃犵狉狅狆狔狉狅狀犮狉犻狊狋犪狋狌犿+犆犪狉犪犵犪狀犪
犽狅狉狊犺犻狀狊犽犻犻+犎犲犱狔狊犪狉狌犿犳狉狌狋犻犮狅狊狌犿+犈犾狔犿狌狊狀狌狋犪狀狊,ACHE)共5种植被恢复模式,并以完全退化的沙地为对
照(CK)。每种植被恢复模式小区面积均为1hm2,各3次重复,随机区组分布。于2010年7月在每个小区内沿
S形采集20点表层0~10cm土壤样品,将同一处理3个小区的土壤样品充分混匀,并去除石块、植物根茬等杂
物,作为一个混合样。采用四分法分出约1kg土样装入无菌袋内,低温条件下带回实验室。将土样分为两部分,
一部分用于测定土壤理化性质,一部分于-70℃条件下保存用于DNA提取。
1.3 土壤理化性质测定
采用烘干法测定土壤含水量,重铬酸钾容量法测定土壤有机质含量,凯氏定氮法测定土壤全氮,KCl浸提-
611 草 业 学 报 第24卷
紫外分光光度法测定土壤硝态氮,KCl浸提-靛酚蓝比色法测定土壤铵态氮,硫酸-高氯酸消煮法测定土壤全
磷,碳酸氢钠浸提-钼锑抗比色法测定土壤速效磷,2.5∶1.0水土比浸提后测定土壤pH。具体步骤见《土壤农
化分析》[18]。各土壤理化指标见李刚等[1920]的研究结果。
1.4 DNA提取、狀犻狉犓 基因扩增和变性梯度凝胶电泳
采用PowerSoil?DNAIsolationKit(MoBioLaboratories,SolanaBeach,CA,USA)按照操作说明从0.25
g土壤样品中提取总DNA,每个处理各3次重复。
狀犻狉犓 基因扩增采用半巢式PCR(seminestedPCR)方法,引物序列及扩增片段大小见表1。第1轮反应:50
μL的反应体系中含有1×Buffer,400μmol/LdNTPs,0.4μmol/L的每种引物,1.25U的TaqHotStart聚合
酶(TaKaRa),50ng的模板DNA。第2轮反应:50μL的反应体系中含有1×Buffer,400μmol/LdNTPs,0.4
μmol/L的每种引物,1.25U的TaqHotStart聚合酶(TaKaRa),以3μL第1轮PCR产物作为模板。PCR反应
均采用“touchdown”方法。第1轮反应条件:95℃,5min;95℃,30s,45℃,40s,72℃,40s,10个循环,每个循环
降低0.5℃;95℃,30s,43℃,40s,72℃,40s,20个循环;72℃,5min。第2轮反应条件:95℃,5min;95℃,30s,
55℃,40s,72℃,40s,10个循环,每个循环降低0.5℃;95℃,30s,53℃,40s,72℃,40s,20个循环;72℃,5min。
表1 狀犻狉犓基因扩增引物
犜犪犫犾犲1 犘狉犻犿犲狉狊狌狊犲犱犳狅狉狀犻狉犓犵犲狀犲犪犿狆犾犻犳犻犮犪狋犻狅狀
半巢式PCRSeminestedPCR 引物Primer 引物序列(5′3′)Primersequence(5′3′) 产物大小Productsize(bp)
第1轮Step1 nirK1F GGMATGGTKCCSTGGCA 515
nirK5R GCCTCGATCAGRTTRTGG
第2轮Step2 nirK1FGC GGMATGGTKCCSTGGCA 430
nirK3R GAACTTGCCGGTVGYCCAGAC
M=A/C;K=G/T;S=G/C;R=A/G;V=G/A/C;Y=C/T.GC 夹 子 GCclamp:CGCCCGCCGCGCCCCGCGCCCGGCCCGC
CGCCCCCGCCCC.
DGGE分析采用DcodeTM通用突变检测系统(BioRad,USA)进行。聚丙烯酰胺凝胶(37.5∶1)浓度为8%,
变性剂梯度为45%~65%[100%变性剂含有7mol/L尿素和40%(v/v)去离子甲酰胺]。将25μLPCR产物和
5μL6×loadingbuffer混合后共同上样,60℃、150V10h。SYBR
TM GreenI(1∶10000)(Invitrogen)染液中染
色25min。用GelDoxXR凝胶成像系统(BioRad,USA)观察并拍照。根据DGGE图谱,共回收26条DGGE
条带。
1.5 DGGE条带测序和系统发育分析
用不带GC夹的nirK1F和nirK3R引物扩增切割后的DGGE条带。将确定为单一条带的PCR产物进行电
泳纯化回收。用pGEMTeasyVector(Promega)对回收的目的片段克隆。由上海生工生物技术公司进行测序,
测序结果经Blast比对分析,获得相近典型菌株序列。将测序获得的狀犻狉犓 基因序列提交至NCBI,获得登录号。
用 Mega5.2.1中的邻接法(NeighborJoining)建立系统发育树。
1.6 数据分析方法
DGGE图谱采用QuantityOne软件(BioRad,USA)进行数字化处理。各样品采用UPGMA法进行聚类分
析,并用Shannon-Wiener指数(犎)和均匀度(犈犎)评价土壤反硝化细菌狀犻狉犓 基因多样性。计算公式如下:
犎=-∑犘犻ln犘犻
犈犎=犎/ln犛
式中,犎 为Shannon-Wiener指数,犈犎 为均匀度指数,犛为DGGE胶中条带数目,犘犻为第犻条带灰度占该样品
总灰度的比率。
采用SPSS16.0进行单因素方差分析(ANOVA)和相关分析(correlationanalysis)。采用Canoco4.5软件
包进行典范对应分析(canonicalcorrespondenceanalysis,CCA)。
711第1期 李刚 等:不同植被恢复模式下呼伦贝尔沙地土壤反硝化细菌狀犻狉犓基因组成结构和多样性研究
2 结果与分析
图1 土壤反硝化细菌狀犻狉犓基因犇犌犌犈图谱
犉犻犵.1 犇犌犌犈犳犻狀犵犲狉狆狉犻狀狋犻狀犵狅犳狊狅犻犾犱犲狀犻狋狉犻犳狔犻狀犵
犫犪犮狋犲狉犻犪犾狀犻狉犓犵犲狀犲
1-26为克隆条带。1-26werethenumbersoftheexcisedbands.
图2 土壤反硝化细菌狀犻狉犓基因聚类分析
犉犻犵.2 犆犾狌狊狋犲狉犪狀犪犾狔狊犻狊狅犳狊狅犻犾犱犲狀犻狋狉犻犳狔犻狀犵犫犪犮狋犲狉犻犪犾狀犻狉犓犵犲狀犲
图中数字表示相似度。Thenumbersinthe
Figurerepresentthesimilarity.
2.1 土壤反硝化细菌狀犻狉犓 基因DGGE图谱分析
DGGE(图1,3次重复之一)结果表明,不同植被
恢复模式土壤反硝化细菌狀犻狉犓 基因DGGE图谱存在
差异。采用聚类分析发现(图2),UA与 UC先聚,相
似性达到73%;再依次与AC、UH和ACHE聚类,相
似性分别为56%,50%和39%;最后再与CK相聚,相
似度仅为20%,这说明5种植被恢复模式土壤反硝化
细菌狀犻狉犓 基因结构组成较对照沙地发生了明显的变
化。王奇赞等[21]认为相似性较好的群体之间的相似
性应高于60%。由此可知,仅UA与UC的反硝化细
菌狀犻狉犓 基因结构组成较相似。条带4,5,6,7,17,18
和19为UA和UC所共有(图1),说明UA和UC恢
复模式下土壤性状相近,具有更多相似的反硝化细菌,
而条带7,18和19又存在于AC和ACHE中,说明这
几种反硝化细菌适应了这4种植被恢复模式所形成的
土壤环境。
2.2 土壤反硝化细菌狀犻狉犓 基因多样性分析
对土壤反硝化细菌狀犻狉犓 基因Shannon指数(犎)
和均匀度指数(犈犎)分析发现(表2),5种植被恢复模
式土壤反硝化细菌狀犻狉犓 基因Shannon指数均显著高
于对照沙地,其中 UH 最高,为1.370,其次为 UA、
UC、AC和ACHE,分别为1.090,1.070,0.900和
0.840。5种植被恢复模式土壤反硝化微生物均匀度
指数与Shannon指数呈相似的变化趋势,均显著高于
对照沙地,依次为 UH>UC>UA>ACHE>AC,其
中UC与UA之间无显著差异。
2.3 土壤反硝化细菌狀犻狉犓 基因测序结果及系统发
育分析
Blast比对结果表明(表3),26个阳性克隆与已知
序列的相似度较高,在81%~99%之间,在数据库中
没有找到完全相同的序列,这些序列可能是之前没有
描述过的基因序列。从系统发育树(图3)可以看出,
26个阳性克隆分属于布鲁氏菌科(Brucelaceae)的苍
白杆菌属(犗犮犺狉狅犫犪犮狋狉狌犿)、根瘤菌科(Rhizobiaceae)的
根瘤菌属(犚犺犻狕狅犫犻狌犿)以及不可培养的细菌。
2.4 土壤反硝化细菌狀犻狉犓 基因与土壤理化因子的
典范对应分析
在CCA分析结果中(表4,图4),前2个排序轴合
并解释了样本中47.2%的样本总变异,这一比例虽然
不高,但其中可能包含丰富信息[22]。在此排序图中,
表2 土壤反硝化细菌狀犻狉犓基因犇犌犌犈条带数、
犛犺犪狀狀狅狀指数和均匀度指数
犜犪犫犾犲2 犖狌犿犫犲狉狅犳犇犌犌犈犫犪狀犱狊,犛犺犪狀狀狅狀犻狀犱犲狓犪狀犱
犲狏犲狀狀犲狊狊犻狀犱犲狓狅犳狀犻狉犓犵犲狀犲狅犳犱犲狀犻狋狉犻犳狔犻狀犵犫犪犮狋犲狉犻犪
处理
Treatment
条带数
Numberofbands
Shannon指数
Shannonindex
均匀度指数
Evennessindex
CK 4.3±0.6a 0.500±0.010a 0.310±0.010a
UA 11.7±0.6d 1.090±0.030d 0.440±0.010d
UH 12.3±1.5d 1.370±0.020e 0.520±0.010e
UC 10.0±1.0c 1.070±0.010d 0.480±0.003d
AC 9.3±0.6c 0.900±0.010c 0.380±0.004b
ACHE 7.3±0.6b 0.840±0.003b 0.410±0.002c
注:同列不同字母表示差异显著(犘<0.05).狀=3.
Note:Differentletterswithinthesamecolumnindicatesignificant
differenceat犘<0.05.狀=3.
811 草 业 学 报 第24卷
表3 犇犌犌犈条带比对结果
犜犪犫犾犲3 犚犲狊狌犾狋狊狅犳狋犺犲犅犾犪狊狋犪狀犪犾狔狊犻狊狅狀狋犺犲犇犌犌犈犫犪狀犱狊犲狇狌犲狀犮犲狊
条带
Sourceofband
登录号
GenBankID
比对菌名称
Speciesnameofthesubject
比对菌登录号
GenBankIDofthesubject
相似度
Similarity(%)
1 KJ184635 不可培养苍白杆菌属Uncultured犗犮犺狉狅犫犪犮狋狉狌犿sp. GU136459.1 99
2 KJ184637 人苍白杆菌犗犮犺狉狅犫犪犮狋狉狌犿犪狀狋犺狉狅狆犻 GU207402.1 99
3 KJ184636 不可培养苍白杆菌属Uncultured犗犮犺狉狅犫犪犮狋狉狌犿sp. GU136459.1 97
4 KJ184638 不可培养细菌Unculturedbacterium AY683856.1 96
5 KJ184639 不可培养细菌Unculturedbacterium HM438925.1 97
6 KJ184640 人苍白杆菌犗犮犺狉狅犫犪犮狋狉狌犿犪狀狋犺狉狅狆犻 GU207402.1 97
7 KJ184641 不可培养苍白杆菌属Uncultured犗犮犺狉狅犫犪犮狋狉狌犿sp. GU136462.1 99
8 KJ184642 不可培养细菌Unculturedbacterium EU790859.1 84
9 KJ184643 不可培养细菌Unculturedbacterium AY121562.1 97
10 KJ184644 不可培养细菌Unculturedbacterium HM061091.1 96
11 KJ184645 不可培养细菌Unculturedbacterium HM061091.1 96
12 KJ184646 不可培养细菌Unculturedbacterium HM061091.1 96
13 KJ184647 不可培养细菌Unculturedbacterium EU790859.1 85
14 KJ184648 人苍白杆菌犗犮犺狉狅犫犪犮狋狉狌犿犪狀狋犺狉狅狆犻 GU207402.1 99
15 KJ184649 根瘤菌属犚犺犻狕狅犫犻狌犿sp. DQ096645.1 81
16 KJ184650 不可培养苍白杆菌属Uncultured犗犮犺狉狅犫犪犮狋狉狌犿sp. GU136459.1 99
17 KJ184651 不可培养细菌Unculturedbacterium GU270538.1 91
18 KJ184652 不可培养细菌Unculturedbacterium HM438925.1 96
19 KJ184653 不可培养苍白杆菌属Uncultured犗犮犺狉狅犫犪犮狋狉狌犿sp. GU207402.1 99
20 KJ184654 根瘤菌属犚犺犻狕狅犫犻狌犿sp. DQ096645.1 83
21 KJ184655 根瘤菌属犚犺犻狕狅犫犻狌犿sp. DQ096645.1 81
22 KJ184656 人苍白杆菌犗犮犺狉狅犫犪犮狋狉狌犿犪狀狋犺狉狅狆犻 GU207402.1 98
23 KJ184657 不可培养细菌Unculturedbacterium AY121567.1 95
24 KJ184658 不可培养细菌Unculturedbacterium JF772658.1 97
25 KJ184659 不可培养细菌Unculturedbacterium JF772658.1 97
26 KJ184660 不可培养细菌Unculturedbacterium HM061091.1 96
第1排序轴由速效磷(AP)、全磷(TP)、全氮(TN)、有机质(OM)、土壤含水量(SM)、pH(正相关)和硝态氮
(NO3-N)、铵态氮(NH4+N)(负相关)组成。第2排序轴由速效磷(AP)、全磷(TP)、全氮(TN)、有机质(OM)
(正相关)和土壤含水量(SM)、pH、硝态氮(NO3-N)、铵态氮(NH4+N)(负相关)组成。
各处理的土质状况以有机质(OM)、全氮(TN)、硝态氮(NO3-N)、铵态氮(NH4+N)、全磷(TP)、速效磷
(AP)、土壤含水量(SM)和pH共8项理化指标来表征。经蒙特卡洛检验(MonteCarlopermutationtest)发现
AP(犘:0.0020)、pH(犘:0.0020)、TN(犘:0.0020)、TP(犘:0.006)和SM(犘:0.006)可以解释反硝化细菌的种
群变化,置信度达到95%以上。如图4所示,AP(狉=0.7803)、TP(狉=0.5793)、TN(狉=0.4539)、SM(狉=
0.3733)、pH(狉=0.7398)和 OM(狉=0.4519)与第1排序轴正相关,NH4+N(狉=-0.0827)和 NO3-N(狉=
-0.0354)与第1排序轴负相关。TN(狉=0.0746)、TP(狉=0.0340)、AP(狉=0.0838)和OM(狉=0.0894)与第2
排序轴正相关,NH4+N(狉=-0.2237)、NO3-N(狉=-0.1241)、SM(狉=-0.0676)和pH(狉=-0.1885)与第2
排序轴负相关。第1排序轴将 UH 与 ACHE和 UA、UC与 AC分开,第2排序轴将 UH、UA和 UC与 AC、
ACHE分开。同时对具有显著影响的各理化因子之间进行相关分析发现,AP、pH、TN、TP和SM之间均极显著
相关(犘<0.01),说明土壤反硝化细菌的狀犻狉犓 基因组成变化是多种理化因子之间共同影响的结果。
911第1期 李刚 等:不同植被恢复模式下呼伦贝尔沙地土壤反硝化细菌狀犻狉犓基因组成结构和多样性研究
图3 土壤反硝化细菌狀犻狉犓基因系统发育树(邻接法)
犉犻犵.3 犖犲犻犵犺犫狅狌狉犼狅犻狀犻狀犵狋狉犲犲犱犲狆犻犮狋犻狀犵狋犺犲狆犺狔犾狅犵犲狀犲狋犻犮狉犲犾犪狋犻狅狀狊犺犻狆狊犪犿狅狀犵狋犺犲狀犻狉犓狊犲狇狌犲狀犮犲狊
Ⅰ:根瘤菌属犚犺犻狕狅犫犻狌犿sp.;Ⅱ:苍白杆菌属犗犮犺狉狅犫犪犮狋狉狌犿sp.;Ⅲ:不可培养细菌Unculturedbacterium;Ⅳ:不可培养细菌Unculturedbacterium.
表4 犆犆犃分析结果
犜犪犫犾犲4 犚犲狊狌犾狋狊狅犳犮犪狀狅狀犻犮犪犾犮狅狉狉犲狊狆狅狀犱犲狀犮犲犪狀犪犾狔狊犻狊(犆犆犃)
项目
Items
第1轴
Axes1
第2轴
Axes2
第3轴
Axes3
第4轴
Axes4
特征值Eigenvalues 0.709 0.598 0.554 0.332
种-环境相关系数Speciesenvironmentcorrelation 0.967 0.901 0.961 0.890
种累积百分比变化率Cumulativepercentagevarianceofspecies(%) 25.6 47.2 67.2 79.2
种-环境累积百分比变化率Cumulativepercentagevarianceofspeciesenvironmentrelation(%) 29.2 53.8 76.6 90.3
3 讨论
由微生物推动的反硝化使地球生物圈中被固定的氮素重新回到大气中去,从而实现整个自然界的氮素循
环[23]。由亚硝酸盐转化为一氧化氮的反应由亚硝酸还原酶催化实现,是反硝化这一生物过程中重要的限速步
021 草 业 学 报 第24卷
骤,因此亚硝酸还原酶基因(狀犻狉犓和狀犻狉犛)成为反硝化
图4 土壤理化因子与反硝化细菌群落结构组成的犆犆犃分析
犉犻犵.4 犆犪狀狅狀犻犮犪犾犮狅狉狉犲狊狆狅狀犱犲狀犮犲犪狀犪犾狔狊犻狊(犆犆犃)狅犳狋犺犲狊狅犻犾
狆犺狔狊犻犮狅犮犺犲犿犻犮犪犾犳犪犮狋狅狉狊犪狀犱犱犲狀犻狋狉犻犳狔犻狀犵犫犪犮狋犲狉犻犪犾
犮狅犿犿狌狀犻狋狔狊狋狉狌犮狋狌狉犲犪狀犱犮狅犿狆狅狊犻狋犻狅狀
OM:有机质 Organicmatter;TN:全氮 Totalnitrogen;NO3-N:
硝态氮Nitratenitrogen;NH4+N:铵态氮 Ammoniumnitrogen;TP:
全磷Totalphosphorus;AP:速效磷Availablephosphorus;SM:土壤含
水量Soilmoisture.
细菌功能基因中研究最多的基因,并作为反硝化的分
子标记用于研究其种群及多样性[24]。由于狀犻狉犓 基因
对环境的响应比狀犻狉犛基因更加敏感[2526]。因此,本文
以狀犻狉犓 基因作为研究对象,通过PCR-DGGE和克
隆测序分析探讨了呼伦贝尔沙地土壤中反硝化细菌
狀犻狉犓 基因多样性变化对不同植被恢复模式的响应。
研究结果表明土壤反硝化细菌狀犻狉犓 基因结构组
成存在明显差异,仅 UA和 UC间土壤反硝化细菌群
落结构组成相似度较高。通过对Shannon指数(犎)
和均匀度指数(犈犎)的分析发现,5种植被恢复模式土
壤反硝化细菌多样性均显著高于对照沙地,说明在植
被恢复过程中,土壤中的反硝化细菌种类明显增加,其
中混播模式下(AC和ACHE)最低。本实验室对土壤
固氮细菌多样性研究发现,混播模式下 (AC 和
ACHE)土壤固氮细菌的Shannon指数和均匀度最
高,而单播模式低于混播模式,且均高于对照沙地[19]。
对土壤氨氧化细菌多样性研究发现,5种植被恢复模
式均高于对照沙地,其中AC模式下Shannon指数和
均匀度最高,而另一种混播模式下Shannon指数和均
匀度仅排在第4位[20]。这些结果说明,土壤中不同的功能微生物对不同植被恢复模式产生了不同的响应。土壤
反硝化细菌多样性的变化与固氮细菌和氨氧化细菌呈现出一定的相反趋势,这一变化是否能够说明植被恢复后
沙地土壤向天然草原演变,仍需进一步研究。
土壤反硝化细菌的活性受到土壤理化性质的影响,典范对应分析表明,AP、TP、TN、SM 和pH对反硝化细
菌区系的影响均达到显著水平,并且这5种理化因子之间均极显著相关,说明反硝化细菌的区系组成受到土壤各
理化因子之间相互关联共同影响。Marschner等[27]认为多个相互联系的环境因子共同控制草原生态系统的微
生物群落结构。Pengthamkeerati等[28]认为土壤微生物的分布、活动以及多样性与土壤微生境的改变密切相关,
而这均是土壤物理性质改变的结果。Wolsing和Prieme[29]应用TRFLP技术对土壤中狀犻狉犓 基因片段的细菌群
落组成分析认为,土壤中狀犻狉犓 基因群落组成的差异不能归因于单一影响因素,反硝化细菌种群对氮含量增加的
影响可能是氮含量增加(直接效应)和其他环境因素(间接效应)共同作用的结果。尽管NO3-被认为是反硝化速
率的重要控制因素,而本研究并未发现狀犻狉犓 基因群落多样性与NO3-含量具有显著相关性,这与 Mergel等[30]
研究发现酸性森林土的土壤剖面中NO3-含量与反硝化细菌群落有很小的相关性具有一致性,表明反硝化细菌
的数量和多样性由多种因素共同调控。
土壤理化性质与地表植物的种类密切相关,同时对土壤微生物群落结构造成影响[3132]。本实验室研究发现,
植被恢复后植物物种数量均有所增加(结果未列出)。植被组成、生长状况、地表覆盖程度、根系分泌物的化学组
成及根系周转速率均存在差异[33],造成土壤养分含量的变化,从而导致土壤反硝化细菌多样性的不同。Baneras
等[34]发现植被类型的改变对狀犻狉犓 反硝化群落具有显著影响。Bremer等[35]也发现植物根系分泌物能够直接影
响土壤微生物的种群组成,反硝化细菌群落结构的形成是根系生物量、植被类型以及土壤理化性质等多种因素共
同作用的结果。
从系统发育树可以看出,部分条带与已知种类细菌相似,具有一定的亲缘关系,但同源性均小于85%;大多
数的条带与未知的反硝化细菌具有近的亲缘关系,这表明它们可能代表新的反硝化细菌基因序列。从系统发育
分析结果可知,不同植被恢复模式下土壤反硝化细菌群落结构主要由α变形菌纲(Alphaproteobacteria)的苍白
121第1期 李刚 等:不同植被恢复模式下呼伦贝尔沙地土壤反硝化细菌狀犻狉犓基因组成结构和多样性研究
杆菌属和根瘤菌属以及不可培养的细菌组成,其中与不可培养细菌相关的片段占到总数的50%以上。在已有的
关于研究不同生物地理条件下反硝化细菌多样性的报道中,多数检测到的片段均为不可培养的细菌。Yoshida
等[25]研究水稻土反硝化微生物群落组成时发现,大部分狀犻狉犛和狀犻狉犓 片段均与不可培养的细菌相关。Falk
等[36]研究海洋水体和沉积物中狀犻狉犛基因的群落多样性时,发现多数狀犻狉犛片段为不可培养细菌狀犻狉犛片段。从反
硝化细菌的结构组成上可以看到,对照沙地中土壤反硝化细菌主要属于苍白杆菌属;UA和UC模式中反硝化细
菌主要属于苍白杆菌属、根瘤菌属和不可培养细菌;AC模式中反硝化细菌主要属于根瘤菌属;而 UH和ACHE
模式中反硝化细菌全部为不可培养细菌。这再次说明不同植被恢复模式下,土壤中的反硝化细菌群落组成差异
明显,而产生这种结果的原因还有待进一步研究。近年来,好氧反硝化作用逐渐被认识,而本研究所发现的苍白
杆菌属和根瘤菌属均属于好氧反硝化细菌,且形成由苍白杆菌属(对照沙地、UA和UC)向根瘤菌属(UA、UC和
AC)一种潜在的过度趋势,这是否说明好氧反硝化细菌在植被恢复过程中发挥了重要作用,仍需要通过采用更为
先进的分子生物学技术进行长期跟踪研究。
本研究首次分析了呼伦贝尔沙地不同植被恢复模式土壤反硝化细菌狀犻狉犓 基因多样性及群落组成,为沙地
植被恢复模式的建立提供了一定的理论依据。由于反硝化细菌在陆地生态系统具有重要的生态位,因此开展此
研究对沙地的恢复和重建具有重要意义。土壤反硝化过程由4种反硝化酶逐步催化完成,而携带不同反硝化基
因的微生物种群对环境因子变化的响应并不一致[37]。因此,为进一步揭示呼伦贝尔沙地植被恢复过程中反硝化
细菌群落组成变化,有必要对不同类型的反硝化微生物开展研究。同时还要结合多种分子生物学技术,如实时定
量PCR、RT-PCR及高通量测序等[38],分别在DNA和RNA水平上对植被恢复过程中土壤反硝化细菌群落结
构、动态变化以及数量变化等进行分析,深入剖析土壤反硝化细菌种群动态变化对植被恢复的响应机制。
犚犲犳犲狉犲狀犮犲:
[1] ZhangYG,LiDQ,WangHM,犲狋犪犾.PreliminarystudyofalpinemeadowsoildenitrifyingbacteriadiversityintheThreeRiversareaofTibet
anPlateau.ChineseScienceBuletin,2006,51(6):715723.
[2] FirestoneMK.Biologicaldenitrification.In:StevensonFJ.NitrogeninAgriculturalSoils[M].Madison,Wisconsin:AmericanSocietyofAg
ronomy,CropScienceSocietyofAmerica,SoilScienceSocietyofAmerica,1982:289326.
[3] RoschC,MergelA,BotheH.Biodiversityofdenitrifyinganddinitrogenfixingbacteriainanacidforestsoil.AppliedandEnvironmentalMi
crobiology,2002,68(8):38183829.
[4] BrakerG,ZhouJ,WuL,犲狋犪犾.Nitritereductasegenes(狀犻狉犓and狀犻狉犛)asfunctionalmarkerstoinvestigatediversityofdenitrifyingbacteriain
Pacificnorthwestmarinesedimentcommunities.AppliedandEnvironmentalMicrobiology,2000,66(5):20962104.
[5] SmithJM,OgramA.Geneticandfunctionalvariationindenitriferpopulationsalongashorttermrestorationchronosequence.AppliedandEn
vironmentalMicrobiology,2008,74:56155620.
[6] WertzS,GoyerC,ZebarthBJ,犲狋犪犾.EffectsoftemperaturesnearthefreezingpointonN2Oemissions,denitrificationandontheabundance
andstructureofnitrifyinganddenitrifyingsoilcommunities.FEMSMicrobiologyandEcology,2013,83(1):242254.
[7] JungJ,YeomJ,HanJ,犲狋犪犾.Seasonalchangesinnitrogencyclegeneabundancesandinbacterialcommunitiesinacidicforestsoils.Journalof
Microbiology,2012,50(3):365373.
[8] HartmannAA,BarnardRL,MarhanS,犲狋犪犾.EffectsofdroughtandNfertilizationonNcyclingintwograsslandsoils.Oecologia,2013,171
(3):705717.
[9] HuangS,ChenC,WuQ,犲狋犪犾.Distributionoftypicaldenitrifyingfunctionalgenesanddiversityofthe狀犻狉犛encodingbacterialcommunityre
latedtoenvironmentalcharacteristicsofriversediments.BiogeosciencesDiscussions,2011,8(3):52515280.
[10] StevenAW,PaulNN,JohnDA,犲狋犪犾.Bacterialcommunitystructureanddenitrifier(狀犻狉gene)abundanceinsoilwaterandgroundwater
beneathagriculturalland.SoilResearch,2011,49:6576.
[11] RobertsonGP,PaulEA,HarwoodRR.Greenhousegasesinintensiveagriculture:contributionsofindividualgasestotheradiativeforcing
oftheatmosphere.Science,2000,289:19221925.
[12] MaY.CausesandcountermeasuresofthedesertificationinHulunbeierGrassland.ScienceandTechnologyInformation,2007,34:653654.
[13] KangB,LiuSR,CaiDX,犲狋犪犾.SoilphysicalandchemicalcharacteristicsunderdifferentvegetationrestorationpatternsinChinasouthsub
tropicalarea.ChineseJournalofAppliedEcology,2010,21(10):24792486.
[14] WanQQ,PiaoQH,DingGD,犲狋犪犾.AnalysisonthereasonofsandydesertificationinHulunberSteppe.ResearchofSoilandWaterConser
vation,2007,14(4):263266.
[15] JinYX,XuB,YangXC,犲狋犪犾.Belowgroundbiomassandfeaturesofenvironmentalfactorsinthedegreeofgrasslanddesertification.Acta
PrataculturaeSinica,2013,22(5):4451.
[16] ZhaoHL,LiuRT,ZhouRL,犲狋犪犾.Propertiesandmechanismsofchangeofsoilmacrofaunacommunitiesinthedesertificationprocessof
Horqinsandygrassland.ActaPrataculturaeSinica,2013,22(3):7077.
221 草 业 学 报 第24卷
[17] LiuRT,YangXG,ChaiYQ,犲狋犪犾.Responseofgrounddwelingarthropodguildstoreseedingandcuttinginartificial犆犪狉犪犵犪狀犪犽狅狉狊犺犻狀狊犽犻犻
plantationsindesertsteppe.ActaPrataculturaeSinica,2013,22(3):7884.
[18] BaoSD.SoilandAgricultureChemistryAnalysis.3rded[M].Beijing:ChinaAgriculturePress,2008:3989.
[19] LiG,WangLJ,LiYJ,犲狋犪犾.EffectsofdifferentvegetationrestorationpatternsonthediversityofsoilnitrogenfixingmicrobesinHulun
beiersandyland,InnerMongoliaofNorthChina.ChineseJournalofAppliedEcology,2013,24(6):16391646.
[20] LiG,XiuW M,WangJ,犲狋犪犾.Communitystructureanddiversityofsoilammoniaoxidizingbacteriaunderdifferentvegetationrestoration
patternsinHulunbeiersandyland,InnerMongolia,China.JournalofAgroEnvironmentScience,2014,33(1):111120.
[21] WangQZ,XuQF,JiangPK,犲狋犪犾.DGGEanalysisofPCRof16SrDNAV3fragmentsofsoilbacteriacommunityinsoilundernaturalbroadleaffor
estinvadedbyphylostachypubescensinpubescensinTianmumountainnaturereserve.ActaPedologicaSinica,2009,46(4):662669.
[22] GauchHG.MultivariateAnalysisinCommunityEcology[M].UK:CambridgeUniversityPress,1982.
[23] SunJG,GaoJL,MaXT,犲狋犪犾.Molecularecologicaltechniquesforthestudyofdenitrifiersandprogressinrelatedareas.SoilandFertilizer
SciencesinChina,2007,(2):712.
[24] GuoLY,ShiF,YangLY.Advancesinfunctionalgenesandmolecularecologyindenitrifiers.MicrobiologyChina,2011,38(4):583590.
[25] YoshidaM,IshiS,OtsukaS,犲狋犪犾.Temporalshiftsindiversityandquantityof狀犻狉犛and狀犻狉犓inaricepaddyfieldsoil.SoilBiologyandBio
chemistry,2009,41:20442051.
[26] DjigalD,BaudoinE,PlilippotL,犲狋犪犾.Shiftsinsize,geneticstructureandactivityofthesoildenitrifiercommunitybynematodegrazing.Eu
ropeanJournalofSoilBiology,2010,46:112118.
[27] MarschnerP,YangCH,LiebereiR,犲狋犪犾.Soilandplantspecificeffectsonbacterialcommunitycompositionintherhizosphere.SoilBiology
andBiochemistry,2001,33(11):14371445.
[28] PengthamkeeratiP,MotavaliPP,KremerRJ.Soilmicrobialactivityandfunctionaldiversitychangedbycompaction,poultrylitterandcrop
pinginaclaypansoil.AppliedSoilEcology,2011,48(1):7180.
[29] WolsingM,PriemeA.Observationofhighseasonalvariationincommunitystructureofdenitrifyingbacteriainarablesoilreceivingartificial
fertilizerandcattlemanurebydeterminingTRFLPof狀犻狉genefragments.FEMSMicrobiologyandEcology,2004,48(2):261271.
[30] MergelA,SchmitzO,MalmannT,犲狋犪犾.Relativeabundanceofdenitrifyinganddinitrogenfixingbacteriainlayersofaforestsoil.FEMS
MicrobiologyandEcology,2001,36(1):3342.
[31] BardgettRD,MawdsleyJL,EdwardsPJ,犲狋犪犾.Plantspeciesandnitrogeneffectsonsoilbiologicalpropertiesoftemperateuplandgrass
lands.FunctionalEcology,1999,13(5):650660.
[32] ShanZH,DingLL,LongRJ,犲狋犪犾.Relationshipbetweensoilmicroorganism,abovegroundvegetation,andsoilenvironmentofdegraded
alpinemeadowintheheadwaterareaoftheYangtzeandYelowRiver,QinghaiTibetanplateau.ActaPrataculturaeSinica,2007,16(1):3440.
[33] SprentJI,ParsonsR.Nitrogenfixationinlegumeandnonlegumetrees.FieldCropsResearch,2000,65(23):183196.
[34] BanerasL,RuizRuedaO,LopezFloresR,犲狋犪犾.Theroleofplanttypeandsalinityintheselectionforthedenitrifyingcommunitystructure
intherhizosphereofwetlandvegetation.InternationalMicrobiology,2012,15(2):8999.
[35] BremerC,BrakerG,MatthiesD,犲狋犪犾.Impactofplantfunctionalgroup,plantspecies,andsamplingtimeonthecompositionof狀犻狉犓type
denitrifiercommunitiesinsoil.AppliedandEnvironmentalMicrobiology,2007,73(21):68766884.
[36] FalkS,HannigM,BrakerG,犲狋犪犾.狀犻狉犛containingdenitrifiercommunitiesinthewatercolumnandsedimentoftheBslticSea.Biogeoscienc
esDiscussions,2006,3:697727.
[37] YinC,FanFL,LiZJ,犲狋犪犾.InfluencesoflongtermapplicationoforganicandinorganicfertilizersonthecompositionandabundanceofnirS
typedenitrifiersinBlackSoil.EnvironmentalScience,2012,33(11):39673975.
[38] WangHT,ZhengTL,YangXR.Molecularecologyresearchprogressforsoildenitrificationandresearchstatusforitsinfluencingfactors.
JournalofAgroEnvironmentScience,2013,32(10):19151924.
参考文献:
[1] 张于光,李迪强,王慧敏,等.青藏高原三江源地区高寒草甸土反硝化细菌多样性的初步研究.科学通报,2006,51(6):715723.
[12] 马琰.呼伦贝尔草原沙化综合治理技术模式研究.科技信息,2007,34:653654.
[13] 康冰,刘世荣,蔡道雄.南亚热带不同植被恢复模式下土壤理化性质.应用生态学报,2010,21(10):24792486.
[14] 万勤琴,朴起亨,丁国栋,等.呼伦贝尔沙地草场沙漠化成因分析.水土保持研究,2007,14(4):263266.
[15] 金云翔,徐斌,杨秀春,等.不同沙化程度草原地下生物量及其环境因素特征.草业学报,2013,22(5):4451.
[16] 赵哈林,刘任涛,周瑞莲,等.沙漠化对科尔沁沙质草地大型土壤动物群落的影响及其成因分析.草业学报,2013,22(3):7077.
[17] 刘任涛,杨新国,柴永青,等.荒漠草原区柠条林地地面节肢动物功能群对补播牧草和平茬措施的响应.草业学报,2013,22(3):7884.
[18] 鲍士旦.土壤农化分析(第三版)[M].北京:中国农业出版社,2008:3989.
[19] 李刚,王丽娟,李玉洁,等.呼伦贝尔沙地不同植被恢复模式对土壤固氮微生物多样性的影响.应用生态学报,2013,24(6):16391646.
[20] 李刚,修伟明,王杰,等.呼伦贝尔沙地不同植被恢复模式土壤氨氧化细菌群落结构及多样性.农业环境科学学报,2014,33(1):111120.
[21] 王奇赞,徐秋芳,姜培坤,等.天目山毛竹入侵阔叶林后土壤细菌群落16SrDNAV3区片段PCR的DGGE分析.土壤学报,2009,46(4):
662669.
[23] 孙建光,高俊莲,马晓彤,等.反硝化微生物分子生态学技术及相关研究进展.中国土壤与肥料,2007,(2):712.
[24] 郭丽芸,时飞,杨柳燕.反硝化菌功能基因及其分子生态学研究进展.微生物学通报,2011,38(4):583590.
[32] 尚占环,丁玲玲,龙瑞军,等.江河源区退化高寒草地土壤微生物与地上植被及土壤环境的关系.草业学报,2007,16(1):3440.
[37] 尹昌,范分良,李兆君,等.长期施用有机和无机肥对黑土狀犻狉犛型反硝化菌种群结构和封堵的影响.环境科学,2012,33(11):39673975.
[38] 王海涛,郑天凌,杨小茹.土壤反硝化的分子生态学研究进展及其影响因素.农业环境科学学报,2013,32(10):19151924.
321第1期 李刚 等:不同植被恢复模式下呼伦贝尔沙地土壤反硝化细菌狀犻狉犓基因组成结构和多样性研究