珠穆朗瑪峰不同海拔梯度上土壤細菌和真菌群落變化特征

張丹丹,張麗梅,沈菊培,旺 姆

1 中國科學院生態環境研究中心 城市與區域國家重點實驗室,北京 100085 2 中國科學院大學,北京 100049 3 西藏農牧學院,林芝 860000

土壤微生物數量和種類繁多,每克土壤中微生物的數量可達幾十億,上萬個物種[1]。微生物是土壤最具有生命力的組分,驅動著土壤中幾乎所有已知的生物反應過程,在有機物質分解、養分循環、污染物轉化與分解、溫室氣體產生和消減等過程中起著重要作用,對土壤肥力形成、土壤生態系統的結構和功能維持乃至全球變化都具有重要影響[1- 3]。鑒于土壤微生物的重要功能及其對外界環境變化的敏感響應,土壤微生物常被用作評價環境變化的重要指標[3- 4]。

全球氣候變化早已成為不爭的事實,據IPCC(2007)報告預測到21世紀末,全球地表溫度可能升高1.1—6.4℃[5]。高寒生態系統及干旱半干旱生態系統作為全球變化的敏感區域,主要分布在北極,北部森林,青藏高原和高山苔原等地,其氣候變暖速度高于地球平均水平;此外,氣候變化通過改變降水格局而加速高寒生態系統退化,影響土壤微生物的群落組成和生物多樣性[6]。因此,高寒生態系統對氣候變化的響應與反饋不僅影響著當地生態系統,也將對全球生態系統及人類生活產生巨大影響。近年來,人們對高寒生態系統土壤微生物展開了大量的研究。Bryant等對美國科羅拉多州落基山脈(海拔2400—3600m)細菌多樣性的研究發現細菌類群豐富度、系統發育多樣性均隨海拔高度上升而降低[7];Margesin等發現高山(海拔2300—2530m)和亞高山(海拔1500—1900m)土壤微生物活性隨海拔增加而降低,且微生物群落組成隨海拔增加變化,一些嗜冷細菌、真菌和革蘭氏陰性細菌主要在高海拔地區出現[8]。Shen等對我國長白山地區的研究也發現土壤細菌群落隨海拔明顯變化,細菌群落組成及優勢門的相對豐度與土壤pH顯著相關[9]。可見海拔相關因素在很大程度上影響著微生物的豐度和群落結構。

青藏高原是世界上海拔最高的高原,地勢險峻多變,地形復雜,素有“世界屋脊”和“地球第三極”之稱,是北半球氣候變化的啟張器和調節器,不僅直接驅動中國東部和西南部氣候的變化,而且對北半球具有巨大的影響,甚至對全球的氣候變化,也具有明顯的敏感性、超前性和調節性[10]。青藏高原板塊構造隆升形成了其獨特的自然地理環境：高寒,輻射強,日照時間長,缺氧,養分濃度低,氣溫隨高度和緯度的升高而降低。因此,表層土壤微生物群落在對抗上述極端環境過程中可能會形成特有的生存適應能力和群落組成差異。此外,由于青藏高原對全球氣候的敏感性,全球變暖和CO2濃度升高可能會對地下微生物產生復雜的影響[11]。如蘆曉飛利用克隆測序技術對西藏米拉山高寒草甸土壤(海拔3600m)細菌和古菌群落多樣性發現該區域存在著豐富的特殊微生物資源,但是由于受到氣候變化和人類活動影響,退化草甸土壤的微生物群落結構發生改變[12]。彭岳林等對藏北高原不同海拔(海拔4584—4956m)草地植物叢枝菌根(AM)真菌種群多樣性的調查研究發現Glomus屬真菌為最優勢屬,且草地類型、海拔高度和土壤pH值是影響AM真菌生態分布的重要因子[13]。李超男等對貢嘎山海拔梯度上(海拔1800—4100m)土壤甲烷氧化菌的研究發現土壤理化性質和氣候變化綜合作用是造成甲烷氧化菌群落結構和多樣性變化的主要因素[14]。

以上研究雖然都是針對青藏高原地區微生物群落開展的,但是主要集中在低海拔地區,尤其在珠峰永久雪線以上的研究非常少。鑒于青藏高原在全球變化中的特殊地位和作用,深入研究不同海拔梯度下細菌和真菌物種多樣性及群落組成沿海拔高度的變化及其規律,對于理解全球氣候變化對青藏高原生態系統的影響和及其響應機制具有重要意義。此外,高寒生態系統因其獨特的自然氣候條件,賦予了一些土壤微生物的獨特優勢,也是微生物資源的潛在寶庫,研究青藏高原珠穆朗瑪峰高原土壤微生物豐度和群落組成為合理開發利用高寒生態系統微生物資源提供參考。以往研究中,發現珠穆朗瑪峰高山土壤中(海拔4000—6550m)氨氧化細菌和古菌的豐度和群落組成隨海拔梯度發生明顯變化,尤其在雪線附近發生顯著轉變[15],但對于土壤中細菌和真菌群落整體組成情況不清楚。本研究利用活菌計數法(CFU)、實時定量PCR(real-time PCR)技術,并結合磷脂脂肪酸PLFA分析法,變性梯度凝膠電泳技術(DGGE)和克隆測序方法研究了青藏高原珠穆朗瑪峰坡底和永久雪線以上高海拔(4000—6550m)地區土壤細菌和真菌群落的組成特征,以期為理解全球變化對青藏高原生態系統的影響以及發掘高寒地區土壤微生物資源提供基礎信息。

1 材料和方法

1.1 土壤樣品采集

本研究土壤樣品為中國科學院第四次珠穆朗瑪峰遠征期間(2005年4月至6月)采集自珠峰北部斜坡。除3個樣品(M1—M3)取自西藏日喀則地區的農田土壤外(28°82′N,29°28′E;海拔4000m左右,年均溫6.3℃),其余樣品(M4—M12)均采自珠峰北坡不同海拔高度(28°01′570″—28°08′385″N,86°51′533″—86°56′686″E)的裸地,海拔高度分別為5350,5400,5700,5850,6000,6150,6300,6450,6550m。每個采樣點均取5個重復次樣點,采樣時用鏟子移除大石頭和積雪層后采集0—10cm深度的土壤或細砂礫,混合后于移動冰箱運回實驗室,部分樣品保存于4℃用于活菌計數和土壤理化性質分析,其余樣品保存到-20℃冰箱用于DNA提取和分析。

1.2 磷脂脂肪酸(PLFA)分析

磷脂脂肪酸的提取方法和分析依照Bardgett等人描述的方法[16]。操作步驟如下：稱取4g土壤樣品于特氟龍管內,用氯仿-甲醇-磷酸緩沖液(1∶2∶0.8)提取總脂類,通過離心分離得到的磷脂脂肪酸經50℃堿性甲酯化后利用氣體色譜Agilent 6820測定土壤微生物磷脂脂肪酸的含量,以標準定性C11到C20混合細菌酸甲酯為外源標準(Bacterial Acid Methyl Esters CP Mix, Matreya Inc.)。每個樣品進行3次重復實驗,脂肪酸的命名方法如Frostegard等人[17]描述。細菌生物量通過3-OH 12:0,14:0, 2-OH 14:0,i- 15:0,和a- 15:0,15:0,i- 16:0,16:1ω9c,16:0,i- 17:0,17:0cy和10Me18:0的脂肪酸進行估算[18- 21],真菌生物量用18:2ω9,12,18:1ω9c,18:1ω9t進行估算[17- 20]。此外,革蘭氏陽性細菌通過i- 15:0、a- 15:0、i- 16:0,i- 17:0和10Me18:0估算,革蘭氏陰性細菌利用3-OH 12:0,14:0, 2-OH 14:0,15:0,16:1ω9c,16:0和17:0cy計算[22- 24]。

1.3 活菌計數

樣品活菌計數采用稀釋平板法進行,為盡可能模仿原生境寡營養條件培養,細菌計數以10倍稀釋的PGY培養基(PGY：蛋白胨10g,酵母膏5g,葡萄糖1g,蒸餾水1L,pH=7.4),真菌計數以含50mg/L硫酸鏈霉素的10倍稀釋的馬鈴薯葡萄糖瓊脂(PDA：馬鈴薯200g,葡萄糖20g,瓊脂20g,蒸餾水1L,pH自然)培養基進行。操作過程為：稱取土壤或砂礫樣品10g,加入90mL無菌水后振蕩30min,取1mL土壤懸浮液轉移到裝有9mL無菌水試管中,充分混合,并進行系列稀釋。吸取不同稀釋程度的菌懸液各0.1mL至培養皿中,涂布均勻。在正式進行計數培養前,選擇了少量樣品進行溫度影響試驗,發現在4℃培養一個月后細菌和真菌的(Colony forming units,CFUs)數比22℃培養明顯少,故在進行正式計數培養時將平板置于22℃黑暗培養7d和14d后進行菌落計數,計算可培養微生物的數量。

1.4 DNA提取和PCR變性梯度凝膠電泳(DGGE)

土壤總DNA提取采用FastDNA SPIN Kit for Soil 試劑盒(Q·BIOgene, Inc., Carlsbad, CA)進行,所有操作依照產品說明書進行,鮮土稱重0.8—1.0g,破碎速度設定為5.5m/s,破碎時間20s。提取的DNA樣品在1%的瓊脂糖凝膠中進行電泳檢測后使用微量分光光度計Nanodrop? ND- 1000測定濃度和純度(NanoDrop Technologies, Wilmington, DE)。

以引物f954-gc和r1369[25]和引物NS1和Fung GC[26](表1)分別對細菌16S rRNA基因和真菌18S rRNA基因進行普通PCR擴增。PCR擴增在PTC- 200 thermal cycler (MJ-Research, Waltham, MA) 進行,擴增程序如表1。PCR反應混合物的體積為50μL,包括1 × PCR buffer, 2mmol/L MgCl2, 400nM正反向引物濃度, 250μM dNTP, 2.5U ExTaqDNA polymerase (TaKaRa Bio Inc. Otsu, Shiga, Japan),0.4mg/mL牛血清白蛋白(BSA)和DNA模板10ng。所獲得的PCR產物在6%聚丙烯酰胺膠中進行DGGE電泳,細菌所用的變性膠濃度范圍是40%—60%,真菌的變性膠濃度范圍為10%—30%。電泳在DCode Universal Detection System Instrument (Bio-Rad Laboratories, Hercules, CA, USA)進行。電泳條件為120v,60℃恒溫條件,7h。依照操作說明將凝膠在1:10000 SYBR Gold 染料 (Invitrogen Molecular Probes, Eugene, USA)染色30min后利用凝膠呈像儀(G:BOX HR,Syngene, Frederick, USA)進行呈像掃描,并用Quantity One (Bio-Rad Laboratories, Hercules, USA)軟件分析結果。

對DGGE圖譜中清晰條帶進行切膠回收,并將膠條轉移到50μL 10mmol/L Tris-HCl (pH 8.0)緩沖液中,反復凍融3次以上。取上清液2μL用無GC夾子原引物再次擴增。PCR產物以膠回收試劑盒(TaKaRa Bio Inc., Shiga, Japan)切膠純化后,連接到pGEM-T Easy Vector上(Promega, Madison, USA)試劑盒連接載體,轉化大腸桿菌 JM109感受態細胞(JM109)(TaKaRa, Japan),涂布到含有氨芐青霉素(Ampicillin)/IPTG/X-Gal的LB (Luria-Bertani)培養基上,37℃下培養16—18h。隨機選取若干白色克隆子,采用菌體直接擴增方式,用pGEM-T Easy Vector通用引物T7/SP6擴增外源插入片段。每個條帶挑選3個陽性克隆子測序,利用DNAMAN 6.0.3.48 (Lynnon Biosoft, USA)進行序列分析,并把所得正確長度的序列提交NCBI數據庫進行Blast序列對比,下載最相似菌種序列作為系統發育樹的參考序列。然后采用Clustal W和Mega 4.0軟件建立Neighbor-Joining系統發育樹,系統發育樹各分支置信度由自舉分析方法(Bootstrap)檢驗,重復1000次。

1.5 熒光定量PCR分析(real-time PCR)

利用實時熒光定量PCR對土壤中的細菌和真菌豐度進行定量分析,所用引物和探針具體信息詳見表1。細菌熒光定量PCR分析采用探針法：使用TaqMan探針TM1389(5′-CTTGTACACACCGCCCGTC- 3′)(Takara Bio Inc. Otsu, Shiga, Japan),定量引物分別為BACT1369F(5′-CGGTGAATACGTTCYCGG- 3′)/PROK1541R(5′-AAGGAGGTGATCCRGCCGCA- 3′)[27];真菌擴增使用引物ITS1f(5′-TCCGTAGGTGAACCTGCGG- 3′)/5.8s(5′-CGCTGCGTTCTTCATCG- 3′)[28]。實時熒光定量PCR反應體系為25μL,細菌16S rRNA基因的擴增使用Premix Ex TaqTM(TaKaRa),真菌ITS區擴增采用SYBR? Premix Ex TaqTM試劑盒(TaKaRa),反應條件和引物濃度參見表1。已提取的DNA經5或10倍稀釋后作為模版DNA添加到混合體系中,其終濃度為1—10ng,同時加入0.4mg/mL的牛血清白蛋白(BSA)用于減少雜質對PCR反應的抑制作用。每個樣品做3個平行。

1.6 數據處理

利用quantity one軟件對DGGE膠圖像中條帶的位置和亮度進行處理,依據數字化結果,計算多樣性指數。采用canonical correspondence analyses (CCAs)分析環境和土壤變量對微生物種間群落結構的影響,在CANOCO 4.5軟件中進行。CFU數據和基因拷貝數經對數轉換后進行ANOVA分析。本研究中的統計分析在SPSS 13.0軟件中進行,其中多組數據間的方差分析采用單因素方差分析(Duncan,P<0.05);相關分析采用Spearman相關分析。

細菌和真菌rRNA基因序列GenBank號為：(遞交中)

表1 本研究中用到的引物,探針和PCR條件

2 結果

2.1 土壤的基本理化性質

土壤的基本理化性質如我們之前的研究報道[15]。海拔高度在4000m左右的農田土壤(M1—M3)pH值為8.6—8.7,海拔高于5000m高度的樣品(M4—M12)pH值在9.0至9.1之間(圖1A)。

圖1 青藏高原珠穆朗瑪峰不同海拔梯度土壤樣品的理化性質及其細菌和真菌的豐度變化Fig.1 Physicochemical properties and abundance changes of bacteria and fungi of soil samples along different altitudes from Mount Everest, on the Tibetan Plateau

在低海拔土壤樣品中(M1—M3,4000m),有機質含量范圍在8.4—12.6g/kg干土之間,土壤全氮在0.51—0.66g/kg干土之間。而高海拔土壤樣品(M4—M12,≥5350m)土壤有機質和全氮明顯低于低海拔土壤樣品,其有機質含量為2.9—6.3g/kg干土,全氮為0.06—0.21g/kg干土(圖1)。

2.2 細菌/真菌活菌計數和16S rRNA基因豐度

利用10倍稀釋的PGY和PDA培養基對細菌和真菌進行的平板計數分析結果如圖1所示。樣品中可培養細菌數量隨海拔增加而降低,在低海拔農田土壤樣品(M1—M3)中達每克干土含(1.53—5.68)×107CFU/g,高海拔土壤/砂礫樣品(M4—M12,≥5350m)(每克干重含7.36×104—8.83×106CFU/g)高1到2個數量級。所有土壤中可培養真菌的數量均較細菌低2個數量級,但變化趨勢與細菌相似,在低海拔農田土壤中為每克干土(2.21—6.09)×105CFU/g之間,高海拔樣品(M4—M12,≥5350m)中為1.67×103—1.16×104CFU/g干土。相關分析結果表明,細菌和真菌的CFU均與海拔高度呈負相關(r1=-0.797,r2=-0.656,P<0.05),與全氮含量呈顯著正相關(r1=0.534,r2=0.565,P<0.05)。

利用實時熒光定量PCR技術對細菌的16S rRNA基因豐度和真菌ITS豐度進行的分析結果表明,細菌16S rRNA基因豐度隨海拔高度增加而降低,在低海拔農田土壤樣品(M1—M3)達4.13×109—1.09×1010拷貝數/g干土,高于高海拔樣品(M4—M12,≥5350m)(8.82×107—5.63×109拷貝數/g干土)1到2個數量級;所有土壤真菌ITS豐度較細菌低1個數量級,但其變化趨勢與細菌相類似,在低海拔土壤樣品中豐度范圍在6.27×108拷貝數/g干土至1.20×109拷貝數/g干土之間,在高海拔土壤樣品中豐度在4.63×106—2.77×108拷貝數/g干土之間變化(圖1)。與活菌計數結果類似,細菌的rRNA基因和真菌ITS的豐度與土壤海拔高度呈顯著負相關(r1=-0.481,r2=-0.543,P<0.05),與土壤有機質(r1=0.606,r2=0.636,P<0.05)和全氮(r1=0.667,r2=0.697,P<0.05)呈正相關(表2)。

表2 土壤基本屬性與多種微生物多樣性相關分析

2.3 細菌和真菌磷脂脂肪酸含量和組分

磷脂脂肪酸分析結果表明,不同海拔梯度樣品的磷脂脂肪酸總量(total PLFA)在20.86和177.58nmol/g之間,其中細菌來源的磷脂脂肪酸含量在8.87—108.42nmol/g之間,真菌來源的為7.82—58.16nmol/g,無論是總的PLFA含量,還是細菌和真菌的PLFA含量均隨著海拔增加而降低(圖1,圖2),與海拔呈顯著負相關(P<0.05),與有機質和全氮含量呈顯著正相關(P<0.05)(表2)。

在供試的12個土壤樣品中,共檢測到15種C12—C18的PLFAs組分,所有樣品中均檢測到飽和脂肪酸(3-OH 12:0,14:0, 2-OH 14:0,i- 15:0, a- 15:0,15:0,i- 16:0, 16:0,i- 17:0和10Me18:0)、單不飽和脂肪酸、多不飽和脂肪酸和環丙烷脂肪酸。隨著海拔高度增加,樣品中PLFA種類和數量逐漸降低,其中低海拔土壤樣品(M1—M3)含有15個單體PLFAs,M4和M5其次,含有12個單體PLFAs,而在海拔高度高于5700m的樣品中(M6—M12),僅含有7種PLFAs,但這7種PLFAs(2-OH 14:0, 16:1ω9c, 16:0, 18:0,18:2ω9,12, 18:1ω9c, 18:1ω9t)在所有海拔高度樣品中均有分布,其中2-OH 14:0, 16:1ω9c, 16:0代表革蘭氏陰性細菌的特征脂肪酸[22-23],18:2ω9,12, 18:1ω9c和18:1ω9t代表真菌的特征脂肪酸[22-23]。而代表革蘭氏陽性細菌的i- 15:0,i- 17:0和代表革蘭氏陰性細菌的17:0△在海拔大于5350m的樣品(M4—M5)中未檢測到,代表革蘭氏陰性細菌的3-OH 12:0,14:0,15:0和代表革蘭氏陽性細菌的i- 15:0和i- 16:0的組份在海拔≥5700m的樣品(M6—M12)中均未檢測到(圖2),表明這些微生物類群不能適應于高海拔環境。

圖2 青藏高原珠穆朗瑪峰不同海拔梯度樣品中的PLFA種類分布Fig.2 The content of PLFA composition in samples from Mount Everest, on the Tibetan PlateauM1—M3代表低海拔農田土壤樣品(4000m),M4—M12代表高海拔土壤樣品(≥5350m)

對磷脂脂肪酸PLFAs組成與環境變量進行的典范對應分析(Canonical correspondence analysis)結果顯示,具有相似海拔高度的樣品聚集在一起,其中,低海拔農田土壤樣品(M1—M3)和高海拔土壤樣品(M4—M12)沿y軸分開,而來自最高海拔的樣品(M10—M12)與海拔稍低的樣品(M4—M9)沿x軸被分開,x軸和y軸分別解釋了87.2%和7.4%的變異。此外,海拔高度、土壤有機質和全氮對PLFAs的組成影響顯著(P<0.05)(圖3)。

圖3 微生物磷脂脂肪酸群落組成與環境變量的典范對應分析(CCA) Fig.3 Canonical correspondence analysis of PLFA composition and environmental variablesM1—M3代表低海拔農田土壤樣品(4000m),M4—M12代表高海拔土壤樣品(≥5350m)

2.4 不同海拔梯度樣品細菌和真菌群落組成

對細菌16S rRNA基因和真菌18S rRNA基因的變性梯度凝膠電泳(DGGE)后,利用quantity one 軟件對DGGE膠圖像中條帶的位置和亮度進行提取后,進行UPGMA聚類分析結果顯示,細菌UPGMA聚類分布圖中低海拔農田土樣(M1—M3)聚成一個小支,與高海拔樣品明顯分開(圖4),表明低海拔農田土壤和高海拔樣品中細菌群落組成明顯不同。對于真菌來說,則只觀察到高海拔樣品(M10—M12)單獨聚集(圖5)。總的來說,低海拔農田土壤樣品(M1—M3)DGGE圖譜更相似,而在高海拔土壤樣品(M4—M12,≥5350m)中,細菌和真菌條帶數量均在M5樣品中最高。

選擇細菌DGGE圖譜中代表低海拔土壤樣品M1和代表高海拔土壤樣品M4、M5、M6和M9清晰條帶進行克隆測序并構建細菌系統發育樹(圖6),結果發現,在M1樣品中檢測到變形菌綱(Proteobacteria)、酸桿菌綱(Acidobacteria)和放線菌綱(Actinobacteria),M4樣品中檢測到變形菌綱(Proteobacteria),M5中檢測到變形菌綱(Proteobacteria),芽單胞菌綱(Gemmatimonadetes),藍藻菌綱(Cyanobacteria),異常球菌屬(Deinococcus-Thermus),厚壁菌門(Fimicutes),裝甲菌門(Armatimonadetes),擬桿菌綱(Bacteroidetes)和放線菌綱(Actinobacteria),M6中檢測到放線菌綱(Actinobacteria),M9中檢測到酸桿菌綱(Acidobacteria),其中5個樣品中處于同一位置的DGGE條帶測序結果一致。從圖6可以看到,所有土樣中均含有區域1(條帶24)和區域3(條帶10)的條帶,測序結果顯示這兩個區域的條帶代表變形菌綱(Proteobacteria),說明變形菌是青藏高原地區的優勢菌群;區域2條帶6和25代表芽單胞菌綱(Gemmatimonadetes)僅分布在高海拔樣品中;此外,區域4條帶3代表的放線菌綱(Actinobacteria)主要分布在低海拔樣品中,在高海拔樣品中較為少見。由于真菌DGGE圖譜中不同海拔高度土壤樣品清晰條帶較少,將所有清晰條帶克隆測序并構建了真菌系統發育樹(圖7),結果顯示真菌序列分布在子囊菌門(Ascomycetes)、擔子菌門(Basidiomycota)和其他真核生物(Cercozoa)三個簇中;結合DGGE圖譜可以看出區域1條帶4,5和30代表Cercozoa僅分布在高海拔樣品中,而區域2條帶17,20代表子囊菌門(Ascomycetes)在所有土樣中均有分布,說明子囊菌門(Ascomycetes)為青藏高原地區的優勢菌群(圖7)。

圖4 青藏高原珠穆朗瑪峰不同海拔梯度下土壤樣品中細菌DGGE圖譜及其UPGMA聚類分布Fig.4 Denaturing gradient gel electrophoresis (DGGE) profiles and the UPGMA dendrogram of bacteria in soils from Mount Everest, on the Tibetan Plateau

圖5 青藏高原珠穆朗瑪峰不同海拔梯度下土壤樣品中真菌DGGE圖譜及其UPGMA聚類分布Fig.5 Denaturing gradient gel electrophoresis (DGGE) profiles and the UPGMA dendrogram of fungi in soils from Mount Everest, on the Tibetan Plateau

圖6 細菌16S rRNA基因系統發育樹Fig.6 Phylogenetic tree of bacterial 16S rRNA gene in soils from Mount Everest, on the Tibetan Plateau

圖7 真菌ITS基因系統發育樹Fig.7 Phylogenetic relationships among fungal sequences in soils from Mount Everest, on the Tibetan Plateau

3 討論

3.1 細菌和真菌數量隨海拔增加的變化特征

本研究中,用10倍稀釋的細菌和真菌培養基進行的活菌計數結果表明低海拔土壤樣品細菌CFU數(1.53—5.68)×107CFU/g明顯高于高海拔樣品(7.36×104—8.83×106CFU/g)。類似地,真菌CFU分布情況與細菌一致,低海拔土壤樣品CFU(2.21—6.09×105CFU/g)較高海拔樣品高1—2個數量級(1.67×103—1.16×104CFU/g)。利用熒光定量PCR技術對細菌16S rRNA和真菌ITS基因豐度的分析結果與活菌和真菌計數結果一致,低海拔土壤樣品細菌和真菌豐度明顯高于高海拔地區土壤樣品1到2個數量級。相關分析結果表明無論細菌和真菌的CFU數量,還是16S rRNA和ITS基因拷貝數均與全氮和有機質顯著正相關,與海拔高度顯著負相關性。此外,細菌和真菌PLFA量隨海拔增加而降低,與土壤海拔呈負相關,與有機質和全氮呈顯著正相關。隨海拔高度增加,土壤養分、氧氣、溫度降低,輻射增加,這些都可能是導致微生物生長代謝減弱、生物量隨海拔增加而減少的決定性因素。結果與以往在青藏高原腹地海拔高度在3400—4300m的三江源自然保護區高寒草原的研究結果一致,說明即使在嚴酷的生存環境下,青藏高原仍存有豐富的微生物類群,且高海拔高寒等極端環境明顯抑制了微生物的活動從而導致其微生物量逐漸降低[31- 32]。此外,CFU和熒光定量PCR結果顯示青藏高原珠穆朗瑪峰土壤樣品中細菌豐度明顯高于真菌豐度一個數量級,且細菌和真菌變化趨勢呈正相關關系(r=0.909,P<0.05),這與任佐華等人對青藏高原高寒草地土壤(海拔3400—4200m)微生物活性和微生物量的研究結果相類似[31-33],導致這種現象的原因可能是細菌和真菌占據的生態位相似,在青藏高原極端生存環境下,細菌由于個體小,繁殖快在資源競爭中處于優勢[34- 36]。

3.2 不同海拔梯度上細菌和真菌的群落組成特征

磷脂脂肪酸(PLFA)常用來表征復雜環境中具有活性的微生物群落組成。在12個供試樣品中共檢測到15種PLFAs,包括代表革蘭氏陽性細菌的脂肪酸組份i- 15:0、a- 15:0、i- 16:0和i- 17:0,代表革蘭氏陰性細菌的組份3-OH 12:0,14:0, 2-OH 14:0,15:0,16:1ω9c,16:0,17:0△和10Me 18:0,和代表真菌的組份18:2ω9,12,18:1ω9c,18:1ω9t。無論是細菌還是真菌的脂肪酸組分均隨海拔增加而減少,主要表現為海拔高度大于5350m(雪線)時,表征革蘭氏陽性細菌的脂肪酸(a- 15:0和i- 17:0)未被檢測到,當海拔高度增加到5700m以上時,表征革蘭氏陽性細菌的所有脂肪酸(i- 15:0、a- 15:0、i- 16:0和i- 17:0)均未被檢測到。同時,表征革蘭氏陰性細菌的脂肪酸17:0△在海拔高度大于5350m時未被檢測到,當海拔增加到5400m時,3-OH12:0和14:0也未被檢測到,但2-OH 14:0,16:1ω9c和16:0脂肪酸在全部樣品中均有分布,說明青藏高原珠穆朗瑪峰土壤中細菌多樣性隨海拔增加而降低,且革蘭氏陽性細菌對海拔高度的變化較革蘭氏陰性細菌更為敏感,可能的原因是革蘭氏陽性細菌細胞壁結構簡單,只含有90%肽聚糖和10%磷壁酸,其對抗外界環境變化的緩沖能力較弱,而革蘭氏陰性細菌細胞壁較為復雜,不僅含有肽聚糖,其細胞壁外還有由脂多糖、磷脂和脂蛋白等若干種蛋白質組成的外膜能夠有效的阻止或減緩有害物質或惡劣環境對細胞的迫害,從而保護細胞的正常生長代謝[37]。此外,表征真菌的PLFAs(18:2ω9,12,18:1ω9c,18:1ω9t)在所有土壤樣品中均有分布,說明海拔高度增加未能影響真菌在高山土壤的分布,與蔡曉布等在西藏高原草地生態系統中觀察到的結果一致[38],這可能是由于真菌在生理上更耐受低溫,以及其通過產生孢子抵御不良環境等所致。該結果與DGGE圖譜的UPGMA聚類分析結果相一致,再次表明不同海拔高度能夠顯著影響細菌群落組成,但對真菌群落變化沒有顯著影響。

對微生物磷脂脂肪酸PLFA組成與環境變量的典范對應分析(Canonical correspondence analysis)結果表明海拔高度、土壤有機質和全氮是決定極端環境下高山土壤微生物群落的組成和多樣性的主要因子。與對微生物量的影響相似,隨著海拔上升,植被逐漸退化,土壤養分減少,不僅影響土壤微生物數量,也會對土壤中微生物類群造成極大的危害,導致土壤微生物多樣性降低[32]。與此同時,高山溫度會隨海拔升高而降低,Lipson等人針對不同季節和海拔梯度下細菌群落結構與溫度的關系研究中發現,低海拔地區微生物群落對溫度較高海拔地區更為敏感[8,39]。

對細菌DGGE主要條帶的克隆測序分析結果顯示,變形菌綱(Proteobacteria)在低海拔和高海拔樣品中均有分布,說明變形菌是青藏高原土壤的優勢菌群;芽單胞菌僅出現在了高海拔樣品中,這可能與芽單胞菌易受環境濕度的影響有關[40]。DeBruyn等通過對南極旱地土壤、高山干旱土壤和高山苔原等13個不同類型土壤中芽單胞菌豐度和群落組成的研究發現,芽單胞菌相對豐度與水分負相關,芽單胞菌適宜在較為干燥的環境中生存[41]。本研究中低海拔樣品主要分布在日喀則地區,其降水含量較高海拔樣品豐沛,相對較高的降水量會使低海拔樣品的芽單胞菌代謝受到抑制,而高海拔相對干燥的環境促進了芽單胞菌的生長繁殖,使芽單胞菌成為高海拔樣品中的優勢菌群。放線菌綱(Actinobacteria)主要分布在低海拔樣品中出現,較少在高海拔樣品中發現,說明放線菌隨著海拔高度的增加而減少,這與以往對西藏土壤放線菌研究結果相似,即放線菌數量隨海拔升高,凍結期增長顯著減少[42]。

對真菌DGGE主要條帶的克隆測序結果顯示子囊菌門(Ascomycetes)廣泛分布在低海拔和高海拔樣品中,在4000m海拔樣品中也檢測到部分擔子菌門(Basidiomycota),但子囊菌門(Ascomycetes)是主要優勢菌群。此外,利用真菌引物NS1/Fung GC還回收到部分與絲足蟲類(Cercozoa)原生生物如Cryothecomonas,Lecythium等高度相似的序列,但這些序列僅分布在高海拔冰雪覆蓋的樣品中,表明這些序列所代表的生物具有較強的耐寒能力,與以往的研究相類似。如Cryothecomonas最早在南極海冰和冰緣區被發現[43- 44],并相繼在冰芯,沉積物和水柱[45- 48]中被發現。此外,Thaler等利用特異性熒光原位雜交技術對北極和亞北極海域中的Cryothecomonas研究發現Cryothecomonas豐度與冰和融雪量顯著正相關,說明Cryothecomonas可以作為北極冰柱融化的敏感指標[48]。雖然本研究選擇DGGE圖譜測序的代表條帶可能不會完全反映土壤中微生物群落結構,但至少代表了樣品中優勢微生物類群的變化特征。

4 結論

本文綜合采用平板活菌計數,定量PCR,PLFA和DGGE技術探討了青藏高原珠穆朗瑪峰北坡從坡底牧草和農業土壤,到海拔6550m高度的高山土壤或風化不完全的砂礫樣品中微生物的生物量和群落組成特征。結果表明受海拔依賴的相關因子,如養分(有機質、總氮)、溫度等的影響,無論是活性細菌和真菌的CFU數,細菌16S rRNA和真菌ITS基因豐度,還是PLFA所代表的微生物量和物種多樣性均隨海拔升高明顯降低。此外,PLFA中代表革蘭氏陽性細菌的組分在高于5350m以上的樣品中缺失,而代表革蘭氏陰性細菌的部分組分和真菌組分在所有樣品中都有分布,表明革蘭氏陽性細菌比革蘭氏陰性細菌和真菌對海拔梯度變化更為敏感,該結果與不同樣品中細菌DGGE指紋圖譜隨海拔高度聚集成簇,而真菌DGGE指紋圖譜隨海拔高度變化無明顯聚集的結果相一致,表明海拔高度顯著影響細菌群落組成,但對真菌群落變化沒有顯著影響。此外,對細菌和真菌DGGE主要條帶的回收測序分析結果顯示變形菌綱是珠峰不同海拔高度土壤/砂礫樣品中的優勢菌群,芽單胞菌是高海拔樣品中的優勢菌群,而放線菌主要分布在低海拔樣品中;真菌以子囊菌門為優勢菌群。此外,利用真菌的引物還在高海拔冰雪覆蓋樣品回收到絲足蟲類原生生物,代表了冰雪覆蓋樣品中所特有真核生物類群。這些結果為認識全球變化對高寒區生態系統的影響以及發掘高寒地區土壤微生物資源提供了重要信息。

參考文獻(References):

[1] 傅伯杰. 中國學科發展戰略: 土壤生物學. 北京: 科學出版社, 2016.

[2] 黃昌勇. 土壤學: 面向21世紀課程教材. 北京: 中國農業出版社, 2000.

[3] Nannipieri P, Ascher J, Ceccherini M T, Landi L, Pietramellara G, Renella G. Microbial diversity and soil functions. European Journal of Soil Science, 2003, 54(4): 655- 670.

[4] Romaniuk R, Giuffré L, Costantini A, Nannipieri P. Assessment of soil microbial diversity measurements as indicators of soil functioning in organic and conventional horticulture systems. Ecological Indicators, 2011, 11(5): 1345- 1353.

[5] Djukic I, Zehetner F, Mentler A, Gerzabek M H. Microbial community composition and activity in different Alpine vegetation zones. Soil Biology and Biochemistry, 2010, 42(2): 155- 161.

[6] Walker D A, Jia G J, Epstein H E, Raynolds M K, Chapin III F S, Copass C, Hinzman L D, Knudson J A, Maier H A, Michaelson G J, Nelson F, Ping C L, Romanovsky V E, Shiklomanov N. Vegetation-soil-thaw-depth relationships along a low-arctic bioclimate gradient, Alaska: synthesis of information from the ATLAS studies. Permafrost and Periglacial Processes, 2003, 14(2): 103- 123.

[7] Bryant J A, Lamanna C, Morlon H, Kerkhoff A J, Enquist B J, Green J L. Microbes on mountainsides: contrasting elevational patterns of bacterial and plant diversity. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(S1): 11505- 11511.

[8] Margesin R, Jud M, Tscherko D, Schinner F. Microbial communities and activities in alpine and subalpine soils. FEMS Microbiology Ecology, 2009, 67(2): 208- 218.

[9] Shen C C, Xiong J B, Zhang H Y, Feng Y Z, Lin X G, Li X Y, Liang W J, Chu H Y. Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biology and Biochemistry, 2013, 57: 204- 211.

[10] Wu S H, Yin Y H, Zheng D, Yang Q Y. Climatic trends over the Tibetan Plateau during 1971- 2000. Journal of Geographical Sciences, 2007, 17(2): 141- 151.

[11] Xu Z H, Chen C R, He J Z, Liu J X. Trends and challenges in soil research 2009: linking global climate change to local long-term forest productivity. Journal of Soils and Sediments, 2009, 9(2): 83- 88.

[12] 蘆曉飛. 西藏米拉山高寒草甸土壤微生物多樣性研究[D]. 北京: 中國農業科學院, 2009.

[13] 蔡曉布, 彭岳林. 西藏高原不同海拔區域叢枝菌根真菌群落的變化. 應用生態學報, 2015, 26(9): 2803- 2810.

[14] 李超男, 李家寶, 李香真. 貢嘎山海拔梯度上不同植被類型土壤甲烷氧化菌群落結構及多樣性. 應用生態學報, 2017, 28(3): 805- 814.

[15] Zhang L M, Wang M, Prosser J I, Zheng Y M, He J Z. Altitude ammonia-oxidizing bacteria and archaea in soils of Mount Everest. FEMS Microbiology Ecology, 2009, 70(2): 208- 217.

[16] Bardgett R D, Hobbs P J, Frosteg?rd ?. Changes in soil fungal: bacterial biomass ratios following reductions in the intensity of management of an upland grassland. Biology and Fertility of Soils, 1996, 22(3): 261- 264.

[17] Frosteg?rd ?, B??th E, Tunlio A. Shifts in the structure of soil microbial communities in limed forests as revealed by phospholipid fatty acid analysis. Soil Biology and Biochemistry, 1993, 25(6): 723- 730.

[18] Kemp P F, Cole J J, Sherr B F, Sherr E B. Handbook of Methods in Aquatic Microbial Ecology. Florida: Lewis Publishers, 1993: 271- 284.

[19] 顏慧, 蔡祖聰, 鐘文輝. 磷脂脂肪酸分析方法及其在土壤微生物多樣性研究中的應用. 土壤學報, 2006, 43(5): 851- 859.

[20] Tunlid A, White D C. Biochemical analysis of biomass, community structure, nutritional status, and metabolic activity of microbial communities in soil // Stotzky G, Bollag J M, eds. Soil Biochemistry. New York, Basel, Hong Kong: Marcel Dekker Inc., 1992, 7: 229- 262.

[21] Harwood J L, Russell N J. Lipids in Plants and Microbes. Netherlands: Springer, 1984.

[22] Stahl P D, Klug M J. Characterization and differentiation of filamentous fungi based on fatty acid composition. Applied and Environmental Microbiology, 1996, 62(11): 4136- 4146.

[23] Klamer M, B??th E. Estimation of conversion factors for fungal biomass determination in compost using ergosterol and PLFA 18: 2ω6, 9. Soil Biology and Biochemistry, 2004, 36(1): 57- 65.

[24] Schelble R T, McDonald G D, Hall J A, Nealson K H. Community structure comparison using FAME analysis of desert varnish and soil, Mojave Desert, California. Geomicrobiology Journal, 2005, 22(7/8): 353- 360.

[25] Yu Z T, Morrison M. Comparisons of different hypervariable regions ofrrsgenes for use in fingerprinting of microbial communities by PCR-Denaturing gradient gel electrophoresis. Applied and Environmental Microbiology, 2004, 70(8): 4800- 4806.

[26] May L A, Smiley B, Schmidt M G. Comparative denaturing gradient gel electrophoresis analysis of fungal communities associated with whole plant corn silage. Canadian Journal of Microbiology, 2001, 47(9): 829- 841.

[27] Suzuki M T, Taylor L T, Delong E F. Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5′-nuclease assays. Applied and Environmental Microbiology, 2000, 66(11): 4605- 4614.

[28] Fierer N, Jackson J A, Vilgalys R, Jackson R B. Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays. Applied and Environmental Microbiology, 2005, 71(7): 4117- 4120.

[29] Hu Y J, Xiang D, Veresoglou S D, Chen F L, Chen Y L, Hao Z P, Zhang X, Chen B D. Soil organic carbon and soil structure are driving microbial abundance and community composition across the arid and semi-arid grasslands in northern China. Soil Biology and Biochemistry, 2014, 77(1): 51- 57.

[30] 吳愉萍. 基于磷脂脂肪酸(PLFA)分析技術的土壤微生物群落結構多樣性的研究[D]. 杭州: 浙江大學, 2009.

[31] 任佐華, 張于光, 李迪強, 肖啟明, 蔡重陽. 三江源地區高寒草原土壤微生物活性和微生物量. 生態學報, 2011, 31(11): 3232- 3238.

[32] 林超峰, 陳占全, 薛泉宏, 來航線, 陳來生, 張登山. 青海三江源區植被退化對土壤養分和微生物區系的影響. 應用與環境生物學報, 2007, 13(6): 788- 793.

[33] 曾智科. 三江源區高寒草甸土壤微生物季節動態及對草地退化的響應[D]. 西寧: 青海師范大學, 2009.

[34] Frey S D, Knorr M, Parrent J L, Simpson R T. Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests. Forest Ecology and Management, 2004, 196(1): 159- 171.

[35] Wallenstein M D, McNulty S, Fernandez I J, Boggs J, Schlesinger W H. Nitrogen fertilization decreases forest soil fungal and bacterial biomass in three long-term experiments. Forest Ecology and Management, 2006, 222(1/3): 459- 468.

[36] Janssens I A, Dieleman W, Luyssaert S, Subke J A, Reichstein M, Ceulemans R, Ciais P, Dolman A J, Grace J, Matteucc G, Papale D, Piao S L, Schulze E D, Tang J, Law B E. Reduction of forest soil respiration in response to nitrogen deposition. Nature Geoscience, 2010, 3(5): 315- 322.

[37] 沈萍, 陳向東. 微生物學(第八版). 北京: 高等教育出版社, 2016.

[38] 蔡曉布, 彭岳林, 馮固, 錢成. 西藏高原草地植物AM真菌多樣性及其環境影響因子研究. 土壤學報, 2005, 42(4): 642- 651.

[39] Lipson D A. Relationships between temperature responses and bacterial community structure along seasonal and altitudinal gradients. FEMS Microbiology Ecology, 2007, 59(2): 418- 427.

[40] Fawaz M N. Revealing the ecological role ofGemmatimonadetesthrough cultivation and molecular analysis of agricultural soils[D]. Knoxville: University of Tennessee, 2013.

[41] DeBruyn J M, Nixon L T, Fawaz M N, Johnson A M, Radosevich M. Global biogeography and quantitative seasonal dynamics ofGemmatimonadetesin soil. Applied and Environmental Microbiology, 2011, 77(17): 6295- 6300.

[42] 薛泉宏, 譚志遠, 朱銘莪, 張曉琳. 西藏土壤放線性菌初步研究. 西北農業大學學報, 1999, 27(1): 28- 32.

[43] Garrison D L, Buck K R. The biota of Antarctic pack ice in the Weddell Sea and Antarctic Peninsula regions. Polar Biology, 1989, 10(3): 211- 219.

[44] Garrison D L, Buck K R. Protozooplankton in the Weddell Sea, Antarctica: abundance and distribution in the ice-edge zone. Polar Biology, 1989, 9(6): 341- 351.

[45] Sime-Ngando T, Juniper S K, Demers S. Ice-brine and planktonic microheterotrophs from Saroma-ko Lagoon, Hokkaido (Japan): quantitative importance and trophodynamics. Journal of Marine Systems, 1997, 11(1/2): 149- 161.

[46] Thomsen H A, Buck K R, Bolt P A, Garrison D L. Fine structure and biology ofCryothecomonasgen. nov. (Protista incertae sedis) from the ice biota. Canadian Journal of Zoology, 1991, 69(4): 1048- 1070.

[47] Tillmann U, Hesse K J, Tillmann A. Large-scale parasitic infection of diatoms in the Northfrisian Wadden Sea. Journal of Sea Research, 1999, 42(3): 255- 261.

[48] Thaler M, Lovejoy C. Distribution and diversity of a protist predatorCryothecomonas(Cercozoa) in Arctic marine waters. Journal of Eukaryotic Microbiology, 2012, 59(4): 291- 299.