川西北不同沙化程度草地土壤細(xì)菌群落特征

楊秉珣, 劉 泉, 董廷旭

(四川綿陽師范學(xué)院 資源環(huán)境工程學(xué)院, 四川 綿陽621000)

微生物是土壤生態(tài)系統(tǒng)的重要組成部分，微生物參與碳、氮循環(huán)，有機(jī)化合物分解以及能量傳輸、運(yùn)送等多種過程的維護(hù)和調(diào)節(jié)，是土壤生態(tài)系統(tǒng)中最具有生命活力的重要組成部分[1-2]。微生物群落多樣性與土壤生態(tài)系統(tǒng)的結(jié)構(gòu)、功能密切相關(guān)，在維持土壤肥力和土壤生態(tài)平衡中發(fā)揮著重要作用，其結(jié)構(gòu)和功能對于周圍環(huán)境條件特別敏感，是土壤變化的重要指示之一[3]。細(xì)菌作為微生物中含量最多、豐富度最高的類群，通常占土壤微生物的70%～90%，具有最為豐富的遺傳多樣性，能有效促進(jìn)有機(jī)質(zhì)分解、營養(yǎng)物質(zhì)的釋放，參與碳、氮等物質(zhì)循環(huán)過程，維持了生態(tài)系統(tǒng)能量流動和物質(zhì)循環(huán)[4-6]。在陸生系統(tǒng)中，細(xì)菌利用生產(chǎn)者提供的有機(jī)物質(zhì)和代謝物作為營養(yǎng)來源，分解為簡單的無機(jī)物，同時這些無機(jī)物又是生產(chǎn)者的營養(yǎng)來源。土壤中細(xì)菌多樣性變化會影響土壤生態(tài)過程，例如養(yǎng)分循環(huán)等[7-8]。此外，土壤中細(xì)菌多樣性變化還會受自然和人為因素的影響，例如植被、土壤養(yǎng)分、土地利用方式等。高通量測序作為第二代測序方法，由于無需構(gòu)建克隆，耗時少，通量高，能夠準(zhǔn)確全面地反映土壤微生物群落分布特征等優(yōu)勢，已經(jīng)逐漸取代傳統(tǒng)的測序方法[9]。

自20世紀(jì)90年代以來，隨著我國實施的大規(guī)模生態(tài)恢復(fù)工程，西北地區(qū)土地沙化和荒漠化得到遏制，然而從目前的研究成果來看，我國土地沙化的重點(diǎn)工作仍集中在北方干旱、半干旱區(qū)域，忽略了對西北地區(qū)草地沙化的研究[10]。川西北高寒草原是全國五大牧區(qū)之一，地處青藏高原東緣，生態(tài)地理位置極其重要，是四川省沙化最嚴(yán)重的區(qū)域，過去20(1994—2013年)年間沙化總面積增加了28.1%，目前草地面積剩余820萬hm2。隨著全球氣候變化和人類掠奪式發(fā)展(濫墾和亂挖)，導(dǎo)致該區(qū)草原植被數(shù)量級覆蓋度急劇下降，生態(tài)環(huán)境也極度惡化，退化沙化草地面積也逐年擴(kuò)大，形成了不同沙化程度的草地類型，這已經(jīng)嚴(yán)重影響了該區(qū)經(jīng)濟(jì)的可持續(xù)發(fā)展和生態(tài)安全問題，因此對該區(qū)不同沙化程度草地的治理和研究已經(jīng)迫在眉睫[11-12]。截止目前，該區(qū)沙化草地的研究大部分集中在如何治理及治理措施的選取等方面，在沙化草地退化過程中土壤養(yǎng)分急劇降低，導(dǎo)致微生物功能降低。大量的研究也表明，草地沙化降低了土壤微生物數(shù)量和活性，不同沙化程度對土壤微生物群落的影響不同，微生物種類及多寡是土壤質(zhì)量的重要指標(biāo)，它們對環(huán)境極為敏感，又是恢復(fù)環(huán)境的先鋒，提高了生態(tài)系統(tǒng)的緩沖能力，微生物指標(biāo)已經(jīng)用來評價退化生態(tài)系統(tǒng)中生物群系與恢復(fù)功能之間的聯(lián)系并能為退化土壤恢復(fù)提供有用信息[11-12]。目前針對不同程度沙化草地土壤微生物群落特征多樣性的具體變化情況尚不明確。鑒于此，本文以典型川西北不同沙化程度草地為研究對象，運(yùn)用高通量測序技術(shù)，研究土壤細(xì)菌多樣性、物種組成和豐富度，并結(jié)合土壤理養(yǎng)分探討影響細(xì)菌群落結(jié)構(gòu)的環(huán)境因素，對發(fā)揮土壤潛在肥力、了解土壤健康狀況，實現(xiàn)該區(qū)植被的管理與可持續(xù)利用有著重要的意義。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)域位于紅原縣(四川省西北部)，地處阿壩藏族羌族自治州中部(31°51′—33°19′N，101°51′—103°23′E)，屬于川西北高原灌叢和若爾蓋高原植被地區(qū)中的植被小區(qū)。該區(qū)屬高原寒冷地區(qū)，沙漠化土地總面積為6 915.4 hm2，海拔3 210～4 857 m，冬季長春秋季短，無夏季，日照充足，寒冷干燥，無霜期40～60 d，晝夜溫差大，年平均氣溫1.1℃，年均降水量753.0 mm，近十年極端最高氣溫為26℃，極端最低氣溫為-36℃，日照達(dá)2 417.9 h，日照率為55%。植被以亞高山草甸為主，沼澤植被占有優(yōu)勢，受高寒氣候的影響，該區(qū)域以高山草甸土、亞高山草甸土和沼澤土為主，伴有寒漠土和風(fēng)沙土，近年來風(fēng)沙土的范圍不斷擴(kuò)大。本研究的取樣點(diǎn)分布于紅原至若爾蓋公路旁，該區(qū)域不同沙化草地均有所表現(xiàn)。

1.2 試驗設(shè)計

2017年8—10月，對研究區(qū)進(jìn)行初步的資料收集和調(diào)查，了解研究區(qū)沙化土地分布格局，采樣時間為植被生長最旺盛期，根據(jù)研究區(qū)現(xiàn)狀選擇不同沙化程度的草地(未沙化草地、輕度沙化草地、中度沙化草地、嚴(yán)重沙化草地)。GPS采集樣地經(jīng)緯度和高程信息，詳細(xì)記錄樣地基礎(chǔ)環(huán)境和物種等指標(biāo)(表1)。每種沙化草地分別選取3個500 m×500 m的樣地，每個樣地隨機(jī)選取1 m×1 m樣方15個，進(jìn)行植物群落樣方調(diào)查。在每個樣方按照五點(diǎn)取樣法收取500 g左右鮮土，鮮土過2 mm篩除去大的礫石和根系，取樣深度為0—10 cm混合土樣(除去表層的枯落物和較多的根系)，充分混勻后液氮保存(-80℃)，帶回實驗室內(nèi)用于測定微生物多樣性；另一部分自然風(fēng)干15 d后過80目篩后測定土壤養(yǎng)分含量。

表1 不同程度沙化樣地特征

注：表中用1代表上坡位，2代表中坡位，3代表下坡位。坡向是順時針旋轉(zhuǎn)的角度來表示(以朝東為起點(diǎn)，為0°)，每45°為一個等級，其變化范圍在0～8之間，數(shù)字越大，表示越干熱[19-20]。

1.3 土壤樣品的測定

1.3.1 土壤養(yǎng)分的測定 一部分土壤樣品經(jīng)自然風(fēng)干20 d后，去除植物根系等雜物過2 mm篩。土壤pH值采用電極電位法(2.5∶1水土比浸提液)；有機(jī)碳和全氮含量采用元素分析儀；速效磷采用NaHCO3浸提—鉬銻抗比色法；堿解氮采用NaOH-H3BO3法；全鉀采用乙酸銨浸提—火焰光度計法[13]。

1.3.2 土壤細(xì)菌測序

(1) 土壤DNA提取和測序。土壤DNA的提取使用上海美吉公司提供的DNA試劑盒(Omega Bio-tek,Inc，Doravilla，GA，USA)利用1%瓊脂糖凝膠電泳進(jìn)行基因組DNA的檢測抽提。參照電泳初步定量的結(jié)果，使用Promega 公司的QuantiFluorTM-ST 藍(lán)色熒光系統(tǒng)對PCR 產(chǎn)物進(jìn)行定量檢測，之后根據(jù)每個樣本的測序量要求進(jìn)行相應(yīng)比例的混合。利用磁珠純化后的PCR產(chǎn)物對目的片段(16 SrRNA V3—V4)進(jìn)行檢測。采用上海美吉生物公司提供的MiSeq測序平臺，利用細(xì)菌的特異性引物(338 F_806 R)進(jìn)行16 SrRNA基因擴(kuò)增，使用AXYGEN公司的AxyPrepDNA凝膠回收試劑盒切膠回收PCR產(chǎn)物，Tris-HCl洗脫；2%瓊脂糖電泳檢測。下機(jī)的原始數(shù)據(jù)去接頭、去污染、去低質(zhì)量，篩選得到可信數(shù)據(jù)，根據(jù)索引號對比至各個樣品。

(2) 序列的優(yōu)化及去雜。高通量測序過程中通常會出現(xiàn)一些點(diǎn)突變和高分子均聚物等測序錯誤，隨著測序長度的增加造成序列末端的質(zhì)量降低，為了得到更高質(zhì)量及更精準(zhǔn)的生物信息分析結(jié)果，則應(yīng)對有效序列進(jìn)行去雜和修剪得到優(yōu)化數(shù)據(jù)。使用軟件Trimmomatic，F(xiàn)LSAH對數(shù)據(jù)進(jìn)行去雜。

(3) 生物信息分析及數(shù)據(jù)處理。根據(jù)Barcode 序列拆分各樣品數(shù)據(jù)，然后截取Barcode和引物的序列使用FLASH對進(jìn)行拼接；高通量序列數(shù)據(jù)使用QⅡME進(jìn)行數(shù)據(jù)處理：利用UPARSE軟件對全部有效序列進(jìn)行聚類，默認(rèn)以99%序列相似度對序列聚類獲得OTUs(operational taxonomic units，操作分類單元)；從每個OTUs選取一個代表性序列，以Greengenes數(shù)據(jù)庫為參考，使用RDP Classifier對OUT進(jìn)行物種注釋，置信水平為80%；最后對樣品數(shù)據(jù)均一化處理，以樣品中數(shù)量最少的為標(biāo)準(zhǔn)，統(tǒng)一抽取有效序列來研究細(xì)菌群落結(jié)構(gòu)。通過計算物種豐富度(OTUs)和土壤細(xì)菌群落多樣性。

計算細(xì)菌豐度(Community richness)的指數(shù)如下[14]：Chao指數(shù)是用Chao1算法估計群落中含OTU數(shù)目的指數(shù)，Chao1在生態(tài)學(xué)中常用來估計物種總數(shù)，由Chao最早提出。Ace指數(shù)用來估計群落中含有OTU數(shù)目的指數(shù)，由Chao提出，是生態(tài)學(xué)中估計物種總數(shù)的常用指數(shù)之一，與Chao1的算法不同。

細(xì)菌多樣性(community diversity)的指數(shù)如下[15]：Simpson指數(shù)用來估算樣品中微生物的多樣性指數(shù)之一，由Simpson提出，在生態(tài)學(xué)中常用來定量描述一個區(qū)域的生物多樣性，Simpson指數(shù)值越大，說明群落多樣性越低。Shannon用來估算樣品中微生物的多樣性指數(shù)之一。它與Simpson多樣性指數(shù)均為常用來反映α多樣性的指數(shù)，Shannon值越大，說明群落多樣性越高。

測序深度指數(shù)如下[16]：Coverage是指各樣品文庫的覆蓋率，其數(shù)值越高，則樣本中序列沒有被測出的概率越低。該指數(shù)反映了本次測序結(jié)果是否代表樣本的真實情況。對OTU列表中獲得的分類信息與豐度進(jìn)行整理，在門和綱分類水平下對各樣品進(jìn)行物種豐度統(tǒng)計及冗余分析(RDA)，可得到樣品中群落組成結(jié)構(gòu)、相似性以及群落結(jié)構(gòu)與環(huán)境因子的關(guān)系。

2 結(jié)果與分析

2.1 不同沙化草地土壤養(yǎng)分特征

由表2可知，川西北不同沙化草地土壤養(yǎng)分具有明顯差異，依次表現(xiàn)為：隨著沙化程度的增加，土壤pH值逐漸增加，而土壤有機(jī)碳、全氮、全鉀、堿解氮和速效磷逐漸降低；其中土壤pH值變化范圍7.56～8.63，均值為8.03，沙化草地pH值顯著高于未沙化草地(p<0.05)；土壤有機(jī)碳變化范圍為9.04～13.2 g/kg，均值為11.25 g/kg，ND和LD差異不顯著(p>0.05)，MD和HD差異不顯著(p>0.05)；土壤全氮變化范圍為1.01～1.23 g/kg，均值為1.11 g/kg，其中沙化草地土壤全氮顯著高于未沙化草地(p<0.05)，MD和HD差異不顯著(p>0.05)；土壤全鉀變化范圍為9.87～13.56 g/kg，均值為11.76 g/kg，ND和LD差民不顯著(p>0.05)，HD顯著高于其他草地(p<0.05)；土壤堿解氮變化范圍18.74～29.85 mg/kg，均值為24.24 mg/kg，不同草地土壤堿解氮差異均顯著(p<0.05)；土壤速效磷變化范圍18.74～29.85 mg/kg，均值為29.19 mg/kg，不同草地土壤速效磷差異均顯著(p<0.05)。

表2 不同沙化草地土壤養(yǎng)分特征

注：不同小寫字母表示差異顯著(p<0.05)，下同。

2.2 不同沙化草地土壤細(xì)菌多樣性

如表3所示，所有樣品的平均覆蓋率為98%，表明該測序效果理想。在3%分類水平下，不同沙化草地土壤細(xì)菌Chao，Ace、Simpson指數(shù)、Shannon指數(shù)、OTU數(shù)量有所差異。OUT數(shù)目、Chao指數(shù)、Ace指數(shù)、Shannon指數(shù)均表現(xiàn)為：ND>LD>MD>HD，其中不同沙化草地土壤細(xì)菌覆蓋率和Simpson指數(shù)差異不顯著(p>0.05)；不同沙化草地土壤Chao、Ace、Simpson指數(shù)、Shannon指數(shù)、OTU均顯著高于未沙化草地。

表3 不同沙化草地土壤細(xì)菌序列統(tǒng)計及多樣性指數(shù)

2.3 不同沙化草地土壤細(xì)菌的群落組成特征

高通量測序結(jié)果顯示不同沙化程度草地土壤樣品中共檢測到細(xì)菌的32個門，65個綱，169個目。所有樣品一共獲得756 234條有效序列，其中最少序列為65 231條，最多序列為70 218條。這些序列的99.1%可以分類到不同的細(xì)菌門類，0.9%與數(shù)據(jù)庫比對分類到古菌門類。如圖1A所示，通過MiSeq高通量測序發(fā)現(xiàn)不同沙化草地土壤中檢測到的主要門有：變形菌門(Proteobacteria)(29.8%～35.2%)、放線菌門(Actinobacteria)(11.3%～19.8%)、酸桿菌門(Acidobacteria)(11.3%～19.8%)、綠彎菌門(Chloroflexi)(8.3%～12.6%)、芽單胞菌門(Gemmatimonadetes)(3.2%～7.5%)、厚壁菌門(Firmicutes)(3.2%～7.2%)、擬桿菌門(Bacteroidetes)(1.6%～4.1%)、硝化螺旋菌門(Nitrospirae)(1.2%～3.9%)共8個門。其中，前6個門在土壤中占主導(dǎo)地位，約占到所有細(xì)菌的89%～95%。此外還發(fā)現(xiàn)了相對豐度小于1%的其他41個稀有門類，它們占總序列5.4%。

如圖1B所示，通過MiSeq高通量測序發(fā)現(xiàn)不同沙化草地下土壤中檢測到的主要綱有：α-變形菌綱(α-Proteobacteria)(21.9%～26.9%)、放線菌綱(Actinobacteria)(18.3%～29.3%)、酸桿菌綱(Acidobacteria)(15.3%～20.5%)、β-變形菌綱(β-Proteobacteria)(12.1%～13.5%)、嗜熱油菌綱(Thermoleophilia)(6.3%～9.8%)、芽單胞菌綱(Gemmatimonadetes)(5.4%～8.3%)、桿菌綱(Bacilli)(2.3%～7.2%)和δ-變形菌綱(δ-Proteobacteria)(1.7%～4.5%)共8個。其中，前6個綱在土壤中占主導(dǎo)地位，約占到所有所有綱的89%～95%。

圖1 不同沙化草地土壤細(xì)菌門和綱分類水平下的相對豐度

2.4 土壤養(yǎng)分與細(xì)菌群落結(jié)構(gòu)的相關(guān)性

土壤可以提供細(xì)菌群落生長繁殖的微環(huán)境，而不同的植被類型通過改變土壤的微環(huán)境間接影響了土壤細(xì)菌群落結(jié)構(gòu)的組成。土壤理化性質(zhì)和細(xì)菌群落結(jié)構(gòu)的相關(guān)性采用Mantel檢驗方法(Mantel tests)分析，從表4可以看出土壤細(xì)菌群落結(jié)構(gòu)與土壤pH值、土壤含水量、土壤有機(jī)碳和總氮呈現(xiàn)顯著相關(guān)(p<0.05)。相關(guān)性分析表明，Proteobacteria，Acidobacteria，F(xiàn)irmicutes，Nitrospirae的相對豐富度與土壤pH值呈顯著負(fù)相關(guān)，Actinobacteria的相對豐富度與土壤pH值呈顯著正相關(guān)；Actinobacteria的相對豐富度與SOC呈顯著負(fù)相關(guān)，Proteobacteria，Acidobacteria，Chloroflexi，F(xiàn)irmicutes的相對豐富度與SOC呈顯著正相關(guān)；Actinobacteria的相對豐富度與TN呈顯著負(fù)相關(guān)，Proteobacteria，Acidobacteria，ChloroflexiFirmicutes的相對豐富度與TN呈顯著正相關(guān)；Chloroflexi的相對豐富度與TK呈顯著正相關(guān)；Chloroflexi，Gemmatimonadetes的相對豐富度與AN呈顯著正相關(guān)；Acidobacteria的相對豐富度與AP呈顯著正相關(guān)。

表4 土壤養(yǎng)分與土壤優(yōu)勢細(xì)菌相對豐度相關(guān)性

注：*，**分別表示在0.05，0.01水平上差異顯著。

2.5 環(huán)境因素對微生物群落的影響

為了探討土壤環(huán)境對微生物群落組成的影響，本文通過將土壤理化性質(zhì)分別與綱和屬分類水平下細(xì)菌的群落組成關(guān)系進(jìn)行冗余分析，利用R語言bioENV篩選出最能體現(xiàn)土壤細(xì)菌群落結(jié)構(gòu)變化的6個土壤因子，研究結(jié)果與Mantel tests方法一致，將經(jīng)過篩選的土壤理化因子與細(xì)菌群落結(jié)構(gòu)進(jìn)行冗余分析RDA)，RDA二維排序圖可以直觀地給出研究對象與環(huán)境變量之間的關(guān)系，排序軸與箭頭連線夾角表示環(huán)境因子與排序軸的相關(guān)性，夾角越小表明關(guān)系越密切，而箭頭連線表示環(huán)境因子與研究對象相關(guān)程度的大小，連線越短，相關(guān)性越小，反之越大。所有樣品因環(huán)境因素的不同而聚類或分離的情況，蒙特卡羅檢驗(Monte Carlo test)發(fā)現(xiàn)，特征根的F=2.76，軸1，軸2，軸3和軸4的特征根的相關(guān)系數(shù)分別為0.953，0.967，0.853，0.981，到軸4的累計合理解釋變量為87.36%，故RDA分析排序結(jié)果可信。環(huán)境變量間夾角的余弦值表示二者的相關(guān)關(guān)系，ND和HD土壤細(xì)菌組成差異較大，MD和HD位于二、三象限，ND和LD位于一、四象限。SOC和pH值的射線較長，表明其對細(xì)菌群落組成影響較大，而AP值的射線較短，說明其對細(xì)菌群落組成影響較?。辉谕寥纏H值，SOC，TN對土壤細(xì)菌群落結(jié)構(gòu)具有較大影響，同時不同沙化草地土壤細(xì)菌群落發(fā)生明顯分異，ND土壤細(xì)菌群落主要處于土壤pH值較高區(qū)域，而HD土壤細(xì)菌群落主要聚類于土壤含水量較高的區(qū)域，不同沙化草地土壤細(xì)菌群落結(jié)構(gòu)較為分散，其中Actinobacteria，Acidobacteria群落與pH值顯著負(fù)相關(guān)，Proteobacteria群落與SOC，TN呈顯著正相關(guān)。

圖2 基于門(A)和綱(B)水平下的土壤細(xì)菌分布與環(huán)境因子的關(guān)系

3 討論與結(jié)論

3.1 不同沙化草地土壤細(xì)菌多樣性及群落組成

川西北草地生態(tài)系統(tǒng)是全球變化的敏感區(qū)域，其微生物組成和功能影響生態(tài)系統(tǒng)轉(zhuǎn)化的方向和進(jìn)程[8]。本研究中，未沙化草地土壤養(yǎng)分含量高于沙化草地，這可能與土壤含水量較高有關(guān)，而土壤水分是調(diào)節(jié)生態(tài)系統(tǒng)土壤微生物物質(zhì)轉(zhuǎn)化的關(guān)鍵因子，沙化草地土壤含水量較低，破壞較為嚴(yán)重，受干擾程度較高，因此，其土壤養(yǎng)分含量較低，其中重度沙化草地土壤含水量最低。土壤細(xì)菌是土壤微生物的重要組成部分，絕大多數(shù)土壤優(yōu)勢細(xì)菌種類基本相同，主要包括10個左右的細(xì)菌類群。Chu等[17]在青藏高原西北部研究表層土壤細(xì)菌群落結(jié)構(gòu)發(fā)現(xiàn)放線菌門、α-變形菌門、酸桿菌門、綠彎菌門和芽單胞菌門5大門類為研究區(qū)的優(yōu)勢類群。厚壁菌門、γ-變形菌門、β-變形菌門、δ-變形菌門、擬桿菌門、浮霉菌門6大門類為平均相對豐度小于5%大于1%的細(xì)菌門類，與本研究基本一致。而Yuan等[18]在念青唐古拉南面沿海拔梯度取樣研究高山草原土壤細(xì)菌群落發(fā)現(xiàn)酸桿菌門、變形菌門、芽單胞菌門為研究范圍內(nèi)最主要的3個細(xì)菌門類。Zhang等[19]在青藏高原北麓河流域分析草甸、草原和荒漠草原土壤細(xì)菌群落結(jié)構(gòu)差異發(fā)現(xiàn)主要優(yōu)勢菌群為變形菌門、放線菌門、酸桿菌門、擬桿菌門4個門類，這可能與當(dāng)?shù)氐臍夂颉⑺?、土壤、植被條件不同有關(guān)，似乎沒有兩種土壤包含完全一致的微生物群落結(jié)構(gòu)。本研究中土壤細(xì)菌優(yōu)勢菌群，與其他研究相比厚壁菌門相對豐度較高，這可能與厚壁菌門更能適應(yīng)川西北沙化草地土壤環(huán)境和利用有限的土壤養(yǎng)分有關(guān)，草地沙化引起了放線菌門相對豐度的增加，放線菌最適生長環(huán)境土壤pH值偏堿性，因此重度沙化草地土壤環(huán)境更適宜放線菌的生長。

大量研究表明，草地沙化、破壞和干擾會改變土壤細(xì)菌的群落多樣性，通過植被恢復(fù)促進(jìn)土壤形成團(tuán)聚體，從而提高土壤的穩(wěn)定性，在植被恢復(fù)的生態(tài)環(huán)境效應(yīng)方面起著重要作用[20-21]。而土壤細(xì)菌的生長主要受土壤的理化性質(zhì)與環(huán)境因子的影響。不同土地利用方式因其管理措施的不同，對土壤理化性質(zhì)的影響存在很大差別，對土壤細(xì)菌群落結(jié)構(gòu)的影響也不同。土壤微生物多樣性指數(shù)表示生物群落中的物種多寡，數(shù)值愈大表示物種越豐富；優(yōu)勢度指數(shù)越大，生物群落內(nèi)的優(yōu)勢種越突出，生物種群豐富和多樣性指數(shù)高是生態(tài)環(huán)境健康穩(wěn)定的重要表現(xiàn)[22-23]。根據(jù)表3可知，未沙化草地土壤細(xì)菌多樣性指數(shù)均高于沙化草地，與前人的研究結(jié)果一致，說明草地沙化降低了土壤微生物多樣性。草地沙化過程中，降低了植被凋落物，從而降低了細(xì)菌對養(yǎng)分的吸收、利用和轉(zhuǎn)化，進(jìn)而影響了土壤細(xì)菌群落生長。

3.2 不同沙化草地土壤細(xì)菌群落組成與土壤養(yǎng)分的關(guān)系

研究區(qū)土壤細(xì)菌群落結(jié)構(gòu)和α多樣性與土壤pH值、土壤有機(jī)碳和全氮顯著相關(guān)。土壤pH值是影響土壤細(xì)菌群落分布的最主要影響因子。Rousk等[24]對具有pH值梯度的農(nóng)田土壤180 m距離(小尺度)內(nèi)土壤微生物群落結(jié)構(gòu)和多樣性采用高通量測序分析，結(jié)果表明土壤細(xì)菌群落結(jié)構(gòu)組成和多樣性與土壤pH值密切相關(guān)。Fierer等[25]利用末端限制性片段長度多態(tài)性分析了美洲大陸從北到南土壤細(xì)菌群落多樣性和生物地理分布模式，發(fā)現(xiàn)土壤細(xì)菌群落豐富度和多樣性與土壤pH值相關(guān)。Shen等[26]在長白山研究6種典型植被類型，采用高通量技術(shù)研究土壤細(xì)菌群落組成和多樣性隨海拔的分布規(guī)律，發(fā)現(xiàn)土壤pH值是驅(qū)動土壤細(xì)菌分布的影響因子，與本研究結(jié)論一致。同時SOC和TN也是土壤細(xì)菌群落結(jié)構(gòu)和多樣性的主要影響因子。而Chu等[17]在對青藏高原西北部大尺度研究表層和亞表層土壤細(xì)菌群落分布時發(fā)現(xiàn)總碳和碳氮比是青藏高原微生物生物地理分布的主要影響因子。植被類型可能通過凋落物和根際分泌物影響土壤微環(huán)境并間接地改變了土壤細(xì)菌的群落結(jié)構(gòu)和多樣性。對于未沙化草地，具有密實的根系，土壤有機(jī)碳含量較高；而重度沙化草地草本覆蓋度低，土壤有機(jī)碳和總氮含量低，土壤微環(huán)境的營養(yǎng)差異可能改變土壤細(xì)菌群落結(jié)構(gòu)和多樣性。當(dāng)土壤微生物有適宜的營養(yǎng)物質(zhì)、土壤pH值和良好的水熱條件，有利于土壤微生物的生長和繁殖[27-29]?？傮w來看，土壤細(xì)菌群落組成與土壤pH值和含水量密切相關(guān)，而當(dāng)土壤pH值相似時，土壤養(yǎng)分可能是細(xì)菌結(jié)構(gòu)的重要影響因子。

隨土壤肥力及環(huán)境狀況的不同，土壤微生物種群數(shù)量也會存在不同程度的差別。本研究中，不同沙化草地之間土壤養(yǎng)分含量不同，土壤優(yōu)勢細(xì)菌門、綱的相對豐富度也不同，說明土壤營養(yǎng)元素含量的變化導(dǎo)致了土壤細(xì)菌的組成及群落結(jié)構(gòu)的變化。細(xì)菌的優(yōu)勢菌門有Proteobacteria，Actinobacteria，Acidobacteria，Chloroflexi，Gemmatimonadetes，F(xiàn)irmicutes，雖然不同沙化草地之間土壤細(xì)菌群落物種組成相似，但土壤細(xì)菌門、綱、種的相對豐度不同，可能是因為植物類型不同以及其向土壤提供的營養(yǎng)物質(zhì)的形式與含量存在差異。Actinobacteria能夠降解復(fù)雜的木質(zhì)素與纖維素，為土壤提供養(yǎng)分，圖1可知，Proteobacteria是最主要的優(yōu)勢菌門，其結(jié)果表示Proteobacteria為堿性土壤中主要優(yōu)勢菌群，在本研究優(yōu)勢菌Proteobacteria門中，α-Proteobacteria是最主要的亞門，其次是β-Proteobacteria，然后是δ-Proteobacteria，與Zhang等[19]研究結(jié)果不一致，在其研究結(jié)果中，β-Proteobacteria是Proteobacteria門中豐富度最高的亞門，這可能是因為所選取的研究區(qū)域不一樣，其土壤結(jié)構(gòu)、含水率等環(huán)境因素差異較大，土壤細(xì)菌的群落組成也就不一樣，優(yōu)勢菌主要亞門也就不相同。RDA分析發(fā)現(xiàn)，α-Proteobacteria主要受pH值SOC的影響，Acidobacteria與SOC顯著相關(guān)，說明草地沙化明顯改善了土壤的理化性質(zhì)，改變了土壤細(xì)菌的群落組成。本研究表明SOC和TN是該研究區(qū)域影響土壤細(xì)菌群落組成與多樣性特征的主要理化性質(zhì)因素，草地沙化明顯改變了土壤細(xì)菌群落組成。