祁連山中段退化高寒草地土壤細菌群落分布特征

李海云，姚拓，馬亞春，張慧榮，路曉雯，楊曉蕾，夏東慧，張建貴，高亞敏

(甘肅農業大學草業學院，甘肅 蘭州 730070)

草地在陸地生態系統中起著重要作用，尤其是在營養物質和能量循環過程扮演著重要角色[1]。近年來由于對草地資源的過度放牧和濫耕亂伐，使草地植被破壞嚴重，水土流失加劇，草原鼠蟲害嚴重[2-3]，導致草地退化、沙化、鹽堿化面積日益擴大，近年來草地退化等生態問題已受到人們的普遍關注和高度重視[4]。放牧活動是影響高寒草地退化的主要驅動因素之一，由于家畜的采食、踐踏和排泄糞便等改變了地表植被覆蓋狀況，造成地表特征的改變，進而影響土壤結構的破壞、養分的損耗、土壤微生物區系和數量的改變，最終導致草地退化[5]。隨著研究和實踐的深入，許多研究者認識到草地退化不僅是地表植被和土壤理化特性的退化，更重要的是草地土壤微生物群落結構和多樣性的變化。土壤微生物作為土壤生態系統的重要組成部分，在土壤有機質分解和營養物質循環等生態系統功能過程中起主導作用[6]。由于土壤微生物對生存的微環境十分敏感，能對土壤環境變化做出快速反應，進而影響植物群落物種多樣性和土壤結構的形成。因此，土壤微生物群落結構和多樣性的變化可作為衡量草地生態系統的健康變化和草地退化或恢復程度的重要指標[7]。目前，關于草地退化對土壤微生物的研究，主要集中在可培養微生物、磷脂脂肪酸(phospholipid fatty acid，PLFA)、PCR-DGGE、16S rRNA基因文庫構建技術等方面[8-11]。這些傳統的研究方法所能反映出的土壤微生物信息十分有限，在很大程度上低估了土壤微生物的物種組成并高估其豐度[12]。隨著高通量測序技術的興起，其具有高準確性、高靈敏度等優勢，能夠較為全面和準確地反映土壤微生物群落結構，被廣泛應用于各研究領域中[13]，彌補了前期土壤微生物傳統研究方法的不足。因此，為探究祁連山中段不同退化草地土壤細菌群落結構變化規律，采用高通量測序技術分析退化高寒草地土壤微生物群落變化趨勢，為高寒草地退化與恢復治理提供理論依據和數據參考價值。

1 材料與方法

1.1 研究區概況

研究區位于祁連山中部地段(37°53′ N、101°45′ E)，隸屬于甘肅省張掖市肅南縣皇城鎮，其南部為冷龍嶺山地，北部為蓋掌達坂山地，東連天祝藏族自治縣，西接山丹軍馬場。海拔2600～3500 m，屬高原大陸性氣候，全年平均氣溫0.6～3.8 ℃，無霜期不足140 d，絕對無霜期僅45～60 d，年降水量361.6 mm，氣候寒冷多變，四季不分明[14]。

1.2 樣地設置及土樣采集

依據草地退化程度劃分相關劃價標準[15]，在2017年8月對肅南縣皇城草地植被進行調查，選擇3種退化草地類型，分別為：輕度退化草地(lightly degraded grassland，LDG)、中度退化草地(moderately degraded grassland，MDG)和重度退化草地(severely degraded grassland，SDG)，具體樣地信息見表1。在每種樣地中，設置3個大樣方(50 m×50 m)，其內隨機設置3個小樣方(50 cm×50 cm)，記錄小樣方中植物種類、物種數、蓋度和高度。地上部生物量采用齊地面刈割，將鮮樣在65 ℃下烘干24 h后稱重。用土鉆(內徑53 mm)采集0～20 cm土壤樣品，土樣混勻后裝入已滅菌袋中迅速帶回，在-80 ℃冰箱中保存以備用。土樣分為兩份，一份用于土壤養分及酶活測定，另一份用于土壤微生物總DNA提取。

表1 樣地基本信息Table 1 Basic information of the soil samples

1.3 土壤理化性質及酶活性測定

采用鮑士旦[16]的方法測定土壤理化性質pH、含水量、電導率、有機質、全量養分和速效養分含量；采用關松蔭[17]的方法測定土壤脫氫酶、氧化還原酶、脲酶、中性磷酸酶、纖維素酶和蔗糖酶活性。

1.4 土壤總DNA提取及細菌16S rRNA基因擴增

采用OMEGA土壤DNA提取試劑盒(上海翊圣生物科技有限公司，上海)對土壤DNA基因組進行提取，每個樣品3次重復并進行混合，以降低DNA提取過程中造成的誤差。以提取的DNA為模板(1 ng·μL-1)，對細菌16S rRNA基因V4～V5區采用帶Barcode的特異性引物[18](515F和907R)進行PCR擴增，每個樣品3次重復。PCR擴增條件參照趙帆等[19]的方法進行擴增。PCR產物根據產物濃度進行等量混樣，充分混勻后使用2%的瓊脂糖凝膠電泳檢測PCR產物并切膠回收，產物送至天津諾禾致源生物信息科技有限公司，采用Illumina HiSeq PE250高通量測序技術平臺進行測序分析。

1.5 生物信息學分析

根據Barcode序列和PCR擴增引物序列從測定數據中拆分出各樣品數據，截去Barcode和引物序列后對每個樣品的reads進行拼接[20]后得到原始測序數據；原始測序數據經質檢和嵌合體去除得到有效數據(effective tags)。利用Uparse軟件對樣品的有效序列進行OTUs聚類(相似度97%以上)，用Mothur方法[21]與SILVA軟件的SSUrRNA數據庫[22]進行物種注釋(閾值：0.8～1.0)；采用PyNAST軟件[23]與GreenGene數據庫中數據信息進行多序列比對，最后對測序數據進行標準化處理，采用Qiime 1.7.0軟件計算Alpha多樣性指數和Beta多樣性。采用R 2.15.3軟件的Vegan軟件包繪制稀釋曲線，并進行無度量多維標定法分析(non-metric multi-dimensional scaling，NMDS)。

1.6 數據分析

所有數據均采用SPSS 21.0軟件進行處理分析，One-Way ANOVA和Duncan氏新復極差法分析差異顯著性。采用CANOCO 4.5軟件對細菌優勢群落進行去趨勢對應分析(DCA)，根據第一排序軸的梯度范圍(lengths of gradient)數值，小于3.0，采用冗余分析(RDA)；為3.0～4.0，采用RDA分析和典型相關分析(CCA)均可；大于4.0，采用CCA分析，對細菌群落與土壤環境指標間的相互關系進行分析。

2 結果與分析

2.1 植被特征

3種高寒草地植被特征間差異顯著(P<0.05)。隨著退化程度的加重，植被蓋度、高度、地上生物量均明顯降低，依次為輕度退化>中度退化>重度退化(表2)；群落物種多樣性和均勻度指數的變化為輕度退化>中度退化>重度退化。

表2 各樣地植被特征Table 2 The basic vegetation characteristics of different soil sampling sites

注：同行不同字母表示差異顯著(P<0.05)，下同。

Note: Different letters in the same line mean significant differences at 0.05 level, the same below.

2.2 土壤理化性質

從表3可以看出，各樣地土壤理化性質及酶活性間均差異顯著(P<0.05)。隨著草地退化程度的加重，土壤含水量、有機碳、有機質、全氮、全磷、全鉀、速效磷、速效氮、速效鉀均呈逐漸減小趨勢；土壤電導率呈逐漸增大趨勢；土壤pH、脲酶、中性磷酸酶、脫氫酶、蔗糖酶和氧化還原酶均呈先增大后減小趨勢；土壤纖維素酶呈先減小后增大趨勢。

2.3 土壤細菌群落豐度與Alpha多樣性分析

通過Illumina HiSeq PE250高通量測序，共得到有效序列257971條，質控過濾和去除嵌合體后得到219017條優質序列，聚類共得到2004個OTUs。各樣品文庫的覆蓋度均在99.6%以上，并結合樣品稀釋曲線均趨于平緩，說明本研究測序數據合理，能夠準確反映出土壤細菌群落的真實信息(圖1)。3種退化草地土壤細菌群落豐富度指數(Chao1和ACE指數)依次為LDG>MDG>SDG(表4)；細菌群落多樣性指數(Simpson和Shannon-Wiener指數)依次為LDG>SDG>MDG。如圖2所示，所有樣品中共有OTUs數目為1317個，其中LDG、MDG和SDG中所特有的OTUs數目分別為112、100和47個。

2.4 不同退化草地土壤細菌群落分布特征分析

2.4.1門水平上的群落組成 通過測序發現，3種退化草地土壤細菌中相對豐度>1%的菌門為：放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)、酸桿菌門(Acidobacteria)、變形菌門(Proteobacteria)、芽單胞菌門(Gemmatimonadetes)、綠彎菌門(Chloroflexi)、疣微菌門(Verrucomicrobia)、浮霉菌門(Planctomycetes)、擬桿菌門(Bacteroidetes)和硝化螺旋菌門(Nitrospirae)，豐度低于1%的類群占1.57%～3.07%。3種退化草地土壤細菌在門分類水平上，相對豐度存在一定的變化趨勢(圖3)。隨著草地退化程度加劇，放線菌門、酸桿菌門、芽單胞菌門、綠彎菌門、浮霉菌門、擬桿菌門相對豐度呈先減小后增大趨勢，厚壁菌門、疣微菌門和硝化螺旋菌門相對豐度呈先增大后減小趨勢，變形菌門相對豐度呈逐漸減小趨勢；功能菌群Proteobacteria/Acidobacteria呈先增大后減小的趨勢。

表3 各樣地土壤基本理化性質Table 3 The basic soil physiochemical indexes of different soil sampling sites

表4 樣品序列數統計、豐富度與多樣性指數Table 4 Sample sequence numbers statistics, richness and diversity index

2.4.2屬水平上的群落組成 在屬分類水平上，相對豐度>1%的菌屬為：乳球菌屬(Lactococcus)、土壤紅桿菌屬(Solirubrobacter)、未分類RB41屬、Gaiella屬、類諾卡氏屬(Nocardioides)、假單胞菌屬(Pseudomonas)、鞘氨醇單胞菌屬(Sphingomonas)、未分類H16屬、芽球菌屬(Blastococcus)、Pseudarthrobacter、鏈球菌屬(Streptococcus)和芽胞桿菌屬(Bacillus)，相對豐度低于1%的類群占48.94%～73.94%。3種退化草地土壤在屬分類水平上，相對豐度存在一定變化趨勢(圖4)。隨著退化程度加劇，乳球菌屬、假單胞菌屬、鞘氨醇單胞菌屬、未分類H16屬、Pseudarthrobacter和鏈球菌屬相對豐度呈先增大后減小趨勢，土壤紅桿菌屬、未分類RB41屬、Gaiella、類諾卡氏屬和芽球菌屬相對豐度呈先減小后增大趨勢，芽胞桿菌屬相對豐度呈逐漸減小趨勢。

2.5 Beta多樣性分析

LDG與MDG樣地間加權Unifrac距離相異系數分別為0.247和0.245(圖5)；LDG與SDG樣地間加權Unifrac距離相異系數分別為0.184和0.216；MDG和SDG樣地間加權Unifrac距離相異系數分別為0.294和0.258。結果表明，3種退化草地土壤間物種多樣性差異為：LDG>SDG>MDG。為了研究不同樣品間的相似性，對樣品進行了UPGMA聚類分析。在門分類水平上，通過加權Unifrac距離的UPGMA聚類分析，LDG與SDG樣地間細菌群落組成及豐度相似性較高，與MDG樣地間細菌群落組成差異較大(圖6)。

圖1 樣品稀釋曲線Fig.1 Rarefaction curves for samples

圖2 樣品韋恩圖Fig.2 Venn diagrams of samples

圖3 門分類水平下的細菌群落相對豐度Fig.3 Relative abundance of bacterial community at phylum level

圖4 屬分類水平下的細菌群落相對豐度Fig.4 Relative abundances of bacterial community at genus level

2.6 細菌優勢類群與土壤環境因子間RDA分析

通過RDA分析表明：第一、二排序軸累計解釋率分別為68.1%和11.9%(圖7)。其中，厚壁菌門與氧化還原酶呈顯著正相關(P<0.05)；放線菌門與中性磷酸酶、蔗糖酶、氧化還原酶、含水量呈極顯著負相關(P<0.01)，與電導率呈極顯著正相關(P<0.01)，與速效磷、速效鉀呈顯著負相關(P<0.05)；變形菌門與脲酶、磷酸酶、脫氫酶、蔗糖酶、含水量、全N、全P、全K、速效N、速效P、速效K、有機質、有機碳呈極顯著正相關(P<0.01)，與pH、電導率和纖維素酶呈極顯著負相關(P<0.01)，與氧化還原酶呈顯著正相關(P<0.05)；酸桿菌門與脲酶、脫氫酶、全N、全P、全K、速效N、有機質、有機碳呈極顯著正相關(P<0.01)，與pH、纖維素酶呈極顯著負相關(P<0.01)，與蔗糖酶、速效P、速效K呈顯著正相關(P<0.05)；芽單胞菌門與pH呈極顯著負相關(P<0.01)；浮霉菌門與全P、速效N呈極顯著正相關(P<0.01)，與pH呈極顯著負相關(P<0.01)，與脲酶、全N、有機質、有機碳呈顯著正相關(P<0.05)；疣微菌門與pH呈顯著正相關(P<0.05)；硝化螺旋菌門與脲酶、全N、全P、速效N、有機質、有機碳呈極顯著負相關(P<0.01)，與pH呈極顯著正相關(P<0.01)；擬桿菌門與pH呈顯著負相關(P<0.05)。綜上，土壤酶活性和土壤理化性質均對細菌群落的分布均有影響，其中，土壤pH是土壤細菌群落分布的主要驅動因子。

圖5 各樣地OTUs間的Heatmap分析Fig.5 Heatmap analysis between OTUs of sample plots level

3 討論

草地植被退化和土壤理化性質的改變，將直接影響土壤微生物的生物量和群落結構多樣性[24-27]。戴雅婷等[13]采用PCR-DGGE技術對內蒙古退化荒漠草原土壤細菌群落結構特征進行研究發現，隨退化程度的加劇，細菌多樣性指數由高到低依次為：輕度>重度>中度，本研究所得結果與其是一致的。本研究中相對豐度較高的細菌菌門為厚壁菌門、放線菌門、變形菌門和酸桿菌門，真菌菌門為子囊菌門和擔子菌門，這與多數有關高寒草地土壤微生物群落組成所得研究結果是基本一致的。Li等[28]對青藏高原地區高寒草地土壤微生物群落結構組成進行研究發現變形菌門，放線菌門和酸桿菌門是該地區的優勢種群。Zhou等[29]對青藏高原地區不同退化高寒草地土壤微生物群落結構組成進行研究發現放線菌門、變形菌門、酸桿菌門和綠彎菌門是該地區的優勢種群。這表明不同高寒地區的退化草地土壤類型間微生物種群結構分布雖存在明顯差異，但優勢菌群基本相同。隨著退化程度加劇，厚壁菌門整體上呈先升高后降低的趨勢。這是因為厚壁菌門適于生長在可利用養分低的環境中，而高寒地區的低溫環境又限制了有機物的可利用性，不利于它們的生長；另外，厚壁菌門在惡劣環境條件時進入休眠狀態，尤其是能形成芽孢和孢子。隨著退化程度加劇，放線菌門、酸桿菌門和變形菌門的豐度呈先減小后增大趨勢，是由于高寒環境中土壤整體營養狀況較低，且土壤pH屬堿性導致的[30]。許多研究表明，在多種土壤生態系統中pH通常與細菌群落結構有顯著的相關性，在近中性時多樣性最高，pH一旦偏離中性，微生物群落會受到環境脅迫壓力，多樣性也因環境脅迫的選擇而降低[31]。變形菌門是土壤中最主要的細菌類群，該類群的代謝活動是土壤中最主要的細菌活動。酸桿菌門一般存在于營養比較匱乏的環境，放線菌門能夠降解纖維素和幾丁質，是土壤養分供給的主要來源，并且產生的孢子能夠抵抗外界不利的環境條件，使得它成為高寒草地土壤的優勢菌群[32]。地上植被通過影響土壤中的有機碳、有機氮、土壤水分、溫度、通氣性和pH值等來影響土壤微生物多樣性。植物群落多樣性越豐富，凋落物和根系分泌物組成就越豐富，土壤微生物多樣性也就越高[33]。植物多樣性的增加往往伴隨著植物生物量的增加，增加了凋落物和根系等有機物質的吸積，從而有利于土壤OC和N的積累[34]。本研究樣地地處高寒環境，不同退化程度草地中土壤微生物群落結構的變化也相應地受到地上植被的影響，植被特征和細菌群落結構多樣性顯著相關(P<0.05)，這說明地上植被群落結構的發展變化對細菌多樣性產生影響。本研究發現，隨著草地退化演替的發展，土壤細菌群落結構組成和多樣性差異較大，這與植被和土壤理化性質是顯著相關的。通過植被、土壤理化性質與細菌群落結構組成間的RDA分析發現，影響祁連山高寒草地土壤細菌群落結構組成的重要驅動因子是pH。大量研究也報道了pH對微生物的分布規律存在著普遍的影響[30]。本次研究也發現，pH對土壤微生物群落有顯著性影響，但微生物種群對pH的響應是不同的，如土壤中酸桿菌門的豐度與pH值呈負相關，放線菌門的豐度與pH值呈正相關[30]。

圖6 門水平上物種組成的UPGMA聚類樹圖Fig.6 UPGMA clustering tree diagram of species composition at the gate

圖7 細菌群落與土壤環境因子的冗余分析Fig.7 Redundancy analysis for bacterial community and soil environmental factors

4 結論

1) 植被蓋度、高度、地上生物量和植物多樣性指數隨退化程度加劇均明顯降低，依次為輕度退化>中度退化>重度退化；各樣地土壤理化性質及酶活性間均差異顯著(P<0.05)。

2) 通過高通量測序，共得到有效序列257971條，優質序列219017條，OTUs 2004個；土壤細菌群落豐富度指數表現為LDG>MDG>SDG；多樣性指數表現為LDG>SDG>MDG；Beta多樣性為LDG>SDG>MDG；3種退化草地土壤細菌群落的變化主要在于微生物量的變化，對細菌群落結構的影響并不明顯。

3) RDA分析表明：土壤酶活性和土壤理化性質均對細菌群落的分布具有影響。其中土壤pH是影響土壤細菌群落分布的主要驅動因子。