丁鈺珮,杜宇佳,高廣磊,3,*,張 英,曹紅雨,朱賓賓,楊思遠(yuǎn),張儆醒,邱 業(yè),劉惠林

1 北京林業(yè)大學(xué)水土保持學(xué)院,水土保持國家林業(yè)和草原局重點(diǎn)實(shí)驗(yàn)室,北京 100083

2 北京師范大學(xué)地理科學(xué)學(xué)部,北京 100875

3 寧夏鹽池毛烏素沙地生態(tài)系統(tǒng)國家定位觀測研究站,鹽池 751500

4 呼倫貝爾市林業(yè)科學(xué)研究所,呼倫貝爾 021000

土壤細(xì)菌參與凋落物分解和腐殖質(zhì)形成等生態(tài)過程,聯(lián)接著地上和地下生態(tài)系統(tǒng),與植物和土壤聯(lián)系密切,是陸地生態(tài)系統(tǒng)物質(zhì)能量循環(huán)的重要環(huán)節(jié),深刻地影響著陸地生態(tài)系統(tǒng)的穩(wěn)定、調(diào)節(jié)和修復(fù)[1- 2]。一方面,植被可以通過凋落物和根系分泌物,改善土壤水分、養(yǎng)分和酸堿度條件,直接或間接地影響土壤細(xì)菌群落結(jié)構(gòu)[3]。另一方面,土壤細(xì)菌能夠通過代謝活動參與土壤生態(tài)過程,為植物提供所需的營養(yǎng)元素,促進(jìn)植物的定植、生長、發(fā)育和演替[4-5]。此外,土壤細(xì)菌對環(huán)境變化還具有高度敏感性,其結(jié)構(gòu)和功能變化能夠有效指示土壤環(huán)境或植被的變化[6]。因此,土壤細(xì)菌是土壤微生物研究的重要和熱點(diǎn)領(lǐng)域,得到了國內(nèi)外學(xué)者的廣泛關(guān)注[2]。

在研究早期,受制于相對落后的實(shí)驗(yàn)條件和研究方法,土壤細(xì)菌研究進(jìn)展較為緩慢。近年來,受益于現(xiàn)代生物技術(shù)的飛速發(fā)展,國內(nèi)外學(xué)者圍繞土壤細(xì)菌開展了大量卓有成效的研究工作,在土壤細(xì)菌群落結(jié)構(gòu)、多樣性及其驅(qū)動因子等方面取得了豐碩的研究成果,極大地深化了我們對于土壤細(xì)菌的認(rèn)知[4,7]。但是,由于土壤細(xì)菌豐富的物種、遺傳和功能多樣性及其與環(huán)境要素強(qiáng)烈的相互作用,復(fù)雜時(shí)空尺度條件和氣候變化背景下土壤細(xì)菌群落結(jié)構(gòu)依然是相關(guān)研究領(lǐng)域的重要內(nèi)容。同時(shí),土壤細(xì)菌功能也得到越來越多的關(guān)注,相關(guān)研究也逐漸轉(zhuǎn)向土壤細(xì)菌功能預(yù)測與分析[8]。其中,土壤細(xì)菌功能與環(huán)境因子的互饋?zhàn)饔贸蔀橥寥兰?xì)菌研究的重要方向[9- 10]。

PICRUSt(Phylogenetic Investigation of Communities by Reconstruction of Unobserved States)是一種基于系統(tǒng)發(fā)育與功能聯(lián)系假設(shè),判斷和預(yù)測標(biāo)記基因功能組成的微生物功能預(yù)測方法[11]。該方法采用16S rRNA測序基因與Greengenes數(shù)據(jù)庫進(jìn)行對比,依據(jù)基因相似性預(yù)測細(xì)菌代謝功能譜。與Tax4Fun和FAPROTAX等分析方法相比,PICRUSt數(shù)據(jù)庫完善,方法便捷,功能先進(jìn),成本低廉,結(jié)果準(zhǔn)確[12],已經(jīng)成功應(yīng)用于林地、草地、沙地等陸地生態(tài)系統(tǒng)土壤微生物功能分析,為深入理解土壤細(xì)菌功能提供了重要途徑[12-13]。

沙地樟子松(Pinussylvestrisvar.mongolica)具有耐寒、耐旱、耐瘠薄等優(yōu)良特性,是我國“三北”地區(qū)重要的生態(tài)環(huán)境建設(shè)樹種[14]。呼倫貝爾沙地是我國四大沙地之一,既是沙地樟子松的原產(chǎn)地,也是沙地樟子松人工林的重要分布區(qū),但營建沙地樟子松人工林對于土壤細(xì)菌群落結(jié)構(gòu)和功能有何影響尚不明確。鑒于此,以呼倫貝爾沙地樟子松人工林為研究對象,采用16S rRNA高通量測序和PICRUSt功能預(yù)測相結(jié)合的研究方法,鑒定分析土壤細(xì)菌群落結(jié)構(gòu)和功能特征,研究揭示土壤理化性質(zhì)對于土壤細(xì)菌群落的影響,以期豐富完善土壤細(xì)菌群落結(jié)構(gòu)、功能與環(huán)境因子互饋?zhàn)饔孟嚓P(guān)研究,并為沙地樟子松人工林的經(jīng)營管理提供理論依據(jù)和科技支撐。

1 研究區(qū)概況

研究區(qū)位于內(nèi)蒙古自治區(qū)呼倫貝爾市海拉爾國家森林公園(49°06′—49°28′ N,119°28′—120°34′ E),海拔約600—800 m,地處呼倫貝爾沙地,屬于中溫帶半干旱大陸性草原氣候。年均降水量約為346.5 mm,多集中于6—8月,年均氣溫-1.2℃,年均無霜期121 d,年均日照時(shí)數(shù)2701.5 h。境內(nèi)河流主要有額爾古納河、海拉爾河和輝河等,干支流總長約960 km。土壤類型為風(fēng)沙土,土層較厚。研究區(qū)林下基本無灌木,多為草本植物,主要包括羊草(Leymuschinensis(Trin.) Tzvel.)、針茅(StipacapillataLinn.)、地榆(SanguisorbaofficinalisL.)、委陵菜(PotentillachinensisSer.)和狗尾草(Setariaviridis(L.) Beauv.)等。

2 材料與方法

2.1 土壤樣品采集

2018年7月,在海拉爾國家森林公園用樹木生長錐測量各林分樹木年齡后,選取無經(jīng)營管理措施、林下植被相似的25 a、34 a和43 a沙地樟子松人工林為實(shí)驗(yàn)組,并選取鄰近地段沙質(zhì)草地作為對照組(表1)。每個(gè)人工林樣地內(nèi)分別選取3株標(biāo)準(zhǔn)木,在樹冠投影處采集土壤樣品,草地則在樣地內(nèi)隨機(jī)選取3個(gè)取樣點(diǎn)采集土壤樣品,樣品采集深度為0—10和10—20 cm。在每個(gè)取樣點(diǎn)的同一土層隨機(jī)采集3個(gè)樣品并充分混和為1個(gè)樣品,4個(gè)樣地共獲得24份土壤樣品。期間,為防止樣品污染,采用75%酒精對取土鉆和鐵鏟滅菌消毒,以保證土壤微生物樣品采集在無菌條件下進(jìn)行。同時(shí),分別用環(huán)刀和鋁盒采集土壤樣品,用于測定土壤孔隙度和土壤含水量。所有土壤樣品置于4℃便攜式保溫箱保存,運(yùn)回實(shí)驗(yàn)室后轉(zhuǎn)入-80℃冰箱中保存。

表1 樣地基本概況

2.2 土壤理化性質(zhì)分析

將各樣地的新鮮土壤樣品自然風(fēng)干后過2 mm篩進(jìn)行土壤理化性質(zhì)的分析。采用烘干法測定土壤含水量(soil water content, SWC),環(huán)刀法測定土壤孔隙度(soil porosity, SP),重鉻酸鉀稀釋熱法測定土壤有機(jī)質(zhì)(soil organic matter, SOM),靛酚藍(lán)比色法測定全氮(total nitrogen, TN),鉬銻抗比色法測定全磷(total phosphorus, TP),火焰分光光度法測定速效鉀(available potassium, AK)[15]。

2.3 土壤細(xì)菌分離鑒定

使用PowerSoil DNA Isolation Kit試劑盒(MoBio Laboratories, Carlsbad, CA)提取土壤細(xì)菌樣品DNA,完成抽提后,利用1%瓊脂糖凝膠電泳檢測抽提的基因組DNA。對16S rRNA基因的V3—V4區(qū)域進(jìn)行PCR擴(kuò)增,引物序列為338F(5′-CCTACGGGAGGCAGCAG-3′)、806R(5′-ATTACCGCGGCTGCTGG-3′)。擴(kuò)增條件設(shè)置：95℃預(yù)變性5 min,95℃變性45 s,55℃退火50 s,72℃延伸45 s,28個(gè)循環(huán)。每個(gè)樣品3次重復(fù),將同一樣品的DNA混合,作為每份土壤樣品的總DNA。將 PCR產(chǎn)物混合后用2%瓊脂糖凝膠電泳檢測,使用AxyPrepDNA 凝膠回收試劑盒(AXYGEN 公司)切膠回收PCR產(chǎn)物,用Tris-HCl緩沖液洗脫后,使用2%瓊脂糖凝膠電泳檢測。

16S rRNA基因的PCR產(chǎn)物構(gòu)建上機(jī)文庫,使用Qubit進(jìn)行檢測定量,Nanodrop 粗檢文庫濃度,Agilent 2100 檢測文庫片段,qPCR 精確定量文庫濃度。在Illumina Miseq平臺對最后的文庫進(jìn)行上機(jī)測序。測序結(jié)果使用uclust、uparse、cdhit等方法對97%的相似水平下的OTU進(jìn)行聚類統(tǒng)計(jì)分析,得到門、綱、目、科、屬各水平群落的生物信息[16]。

2.4 數(shù)據(jù)處理與分析

采用SPSS 23.0計(jì)算平均值和標(biāo)準(zhǔn)差,進(jìn)行單因素方差分析(One-way ANOVA),顯著性水平0.05。采用Mothur軟件計(jì)算Chao 1豐富度指數(shù)、Shannon多樣性指數(shù)和Pielou均勻度指數(shù)。將由QIIME軟件得到的biom文件上傳到Galaxy網(wǎng)站進(jìn)行PICRUSt功能基因預(yù)測分析,得到土壤細(xì)菌功能基因組成。選擇土壤細(xì)菌相對豐度與土壤性質(zhì)為數(shù)據(jù)矩陣,采用Canoco for Windows 4.5基于線性模型進(jìn)行冗余分析(Redundancy Analysis, RDA),并進(jìn)行蒙特卡洛檢驗(yàn)(Monte-Carlo permutation test),檢驗(yàn)重復(fù)數(shù)為499。

3 結(jié)果與分析

3.1 土壤細(xì)菌群落結(jié)構(gòu)與多樣性

呼倫貝爾沙地樟子松人工林和沙質(zhì)草地土壤中共獲得土壤細(xì)菌35門92綱109目210科267屬,其中相對豐度≥1%土壤細(xì)菌門類均為12個(gè)(圖1)。變形菌門(Proteobacteria)、放線菌門(Actinobacteria)和酸桿菌門(Acidobacteria)是研究區(qū)的主要土壤細(xì)菌優(yōu)勢門,平均相對豐度分別為24.29%±3.39%、23.72%±4.10%和23.40%±2.55%。隨著林齡的增加,沙地樟子松人工林土壤優(yōu)勢細(xì)菌門相對豐度呈現(xiàn)波動變化趨勢,25 a人工林0—10 cm層次變形菌門相對豐度顯著高于43 a人工林(P<0.05)；各林地間放線菌門和酸桿菌門相對豐度不存在顯著差異(P>0.05)。沙質(zhì)草地變形菌門和酸桿菌門相對豐度與人工林存在顯著差異(P<0.05),變形菌門相對豐度比人工林低6.05%,酸桿菌門相對豐度比人工林高4.54%。各樣地不同土層間土壤細(xì)菌組成較為相似,相對豐度無明顯變化。

圖1 呼倫貝爾沙地樟子松人工林土壤細(xì)菌相對豐度熱圖(門水平)

沙地樟子松人工林土壤細(xì)菌群落的Chao 1、Shannon和Pielou指數(shù)不存在顯著差異(P>0.05)(表2)。沙地樟子松人工林土壤細(xì)菌群落的Chao 1和Shannon指數(shù)在43 a人工林中最高,在25 a人工林中最低；Pielou指數(shù)在各林地中基本一致,土壤細(xì)菌分布較為均勻。沙質(zhì)草地土壤細(xì)菌群落的Chao 1和Shannon指數(shù)與人工林不存在顯著差異(P>0.05),Pielou指數(shù)僅在0—10 cm土層與25 a和34 a人工林存在顯著差異(P<0.05)。此外,研究區(qū)土壤細(xì)菌群落的Chao 1、Shannon和Pielou指數(shù)均隨土層深度增加而降低,說明表層土壤細(xì)菌群落多樣性水平相對較高,且分布更為均勻。

表2 呼倫貝爾沙地樟子松人工林細(xì)菌群落的Alpha多樣性指數(shù)

3.2 土壤細(xì)菌群落結(jié)構(gòu)與土壤理化因子的相關(guān)性

呼倫貝爾沙地樟子松人工林土壤理化性質(zhì)波動性較強(qiáng)(表3)。隨著林齡的增加,表層土壤含水量和有機(jī)質(zhì)含量顯著降低(P<0.05),土壤孔隙度先降低后升高,全氮、全磷和速效鉀含量均表現(xiàn)為先增加再減少。沙質(zhì)草地土壤的全氮和速效鉀含量顯著高于人工林(P<0.05),土壤孔隙度和全磷含量與人工林不存在顯著差異(P>0.05)。

表3 呼倫貝爾沙地樟子松人工林土壤理化性質(zhì)

冗余分析結(jié)果顯示(圖2),第一排序軸特征值為42.4%,第二排序軸特征值為10.8%,前兩軸累積解釋量占特征值的總和50%以上,RDA前兩軸能很好地反映土壤細(xì)菌群落與土壤因子之間關(guān)系的大部分信息[17- 18]。樣本點(diǎn)分布相對離散,說明不同樣地土壤細(xì)菌群落結(jié)構(gòu)存在差異且受土壤理化性質(zhì)影響較大。經(jīng)蒙特卡洛檢驗(yàn),速效鉀、全磷和全氮是研究區(qū)土壤細(xì)菌群落的主要影響因子(P<0.05)。不同林齡沙地樟子松人工林土壤細(xì)菌群落結(jié)構(gòu)與土壤理化因子的相關(guān)性有差異,25 a和34 a人工林土壤細(xì)菌的主要影響因子為全氮(P<0.05),43 a人工林土壤細(xì)菌的主要影響因子為土壤含水量(P<0.05)。沙質(zhì)草地土壤細(xì)菌群落的主要影響因子為全磷(P<0.05)。此外,各樣地0—10 cm土層細(xì)菌群落與土壤含水量呈正相關(guān),10—20 cm土層情況相反。

3.3 土壤細(xì)菌功能預(yù)測

呼倫貝爾沙地樟子松人工林土壤細(xì)菌共得到5個(gè)生物代謝通路(圖3)。其中,環(huán)境信息處理、代謝和遺傳信息處理是人工林土壤細(xì)菌功能的重要組成,其相對豐度分別為32.62%±0.98%、26.43%±0.31%和8.13%±0.28%。沙質(zhì)草地一級代謝功能種類與人工林相同。

沙地樟子松人工林中,43 a人工林0—10 cm土層環(huán)境信息處理功能基因相對豐度顯著低于其他林地,代謝、遺傳信息處理和有機(jī)系統(tǒng)功能基因相對豐度顯著高于其他林地(P<0.05)。與沙質(zhì)草地相比,25 a和34 a人工林0—10 cm土層環(huán)境信息處理和細(xì)胞過程功能基因相對豐度顯著高于沙質(zhì)草地,代謝、遺傳信息處理和有機(jī)系統(tǒng)功能基因相對豐度顯著低于沙質(zhì)草地(P<0.05)。此外,研究區(qū)土壤細(xì)菌代謝、遺傳信息處理、有機(jī)系統(tǒng)和細(xì)胞過程功能基因相對豐度均表現(xiàn)為0—10 cm>10—20 cm,沙質(zhì)草地和43 a人工林土壤細(xì)菌環(huán)境信息處理功能基因相對豐度表現(xiàn)為0—10 cm<10—20 cm。

圖2 呼倫貝爾沙地土壤細(xì)菌群落與環(huán)境因子的冗余分析

圖3 呼倫貝爾沙地樟子松人工林土壤細(xì)菌功能預(yù)測(一級功能層)

沙地樟子松人工林和沙質(zhì)草地土壤細(xì)菌功能基因的5個(gè)一級功能層共有31個(gè)子功能。其中,共有10個(gè)二級功能層相對豐度>1%(圖4),且膜運(yùn)輸功能占主導(dǎo)地位,其相對豐度為32.60%±0.98%,其他功能相對豐度均低于10%。沙地樟子松人工林中,43 a人工林0—10 cm土壤細(xì)菌復(fù)制和修復(fù)、酶家族、核苷酸代謝、能量代謝和碳水化合物代謝的功能基因相對豐度顯著高于其他林地,膜運(yùn)輸功能基因相對豐度顯著低于其他林地(P<0.05)。25 a與34 a人工林所有二級功能相對豐度均不存在顯著差異(P>0.05)。與沙質(zhì)草地相比,25 a和34 a人工林0—10 cm土層膜運(yùn)輸功能基因相對豐度顯著高于沙質(zhì)草地,其他二級功能基因相對豐度均顯著低于沙質(zhì)草地(P<0.05)。此外,人工林和沙質(zhì)草地各功能基因在不同土層深度相對豐度基本一致,無明顯變化規(guī)律。

圖4 呼倫貝爾沙地樟子松人工林土壤細(xì)菌功能預(yù)測熱圖(二級功能層)

4 討論

4.1 呼倫貝爾沙地樟子松人工林土壤細(xì)菌群落結(jié)構(gòu)特征

呼倫貝爾沙地樟子松人工林土壤細(xì)菌優(yōu)勢門為變形菌門、放線菌門和酸桿菌門,與此前學(xué)者在呼倫貝爾地區(qū)的研究結(jié)果一致[12]。酸桿菌屬于典型的寡營養(yǎng)型細(xì)菌,而變形菌屬于富營養(yǎng)型細(xì)菌,適生于養(yǎng)分含量高的土壤環(huán)境[19]。因此,變形菌門和酸桿菌門相對豐度可以用作土壤質(zhì)量評價(jià)[20]。人工林變形菌門相對豐度顯著高于沙質(zhì)草地,酸桿菌門顯著低于沙質(zhì)草地,這是由于人工林和沙質(zhì)草地的環(huán)境條件明顯不同,人工林產(chǎn)生的大量枯落物以及根系分泌物能夠提高土壤養(yǎng)分含量,導(dǎo)致土壤細(xì)菌生活的微生境及養(yǎng)分需求有差異,進(jìn)而改變細(xì)菌類群比例[21]。但隨林齡的增加,研究區(qū)沙地樟子松人工林中變形菌門相對豐度逐漸降低,這表明土壤養(yǎng)分狀況下降,可能與植被根系生物量的衰退以及凋落物分解速度降低有關(guān)[22]。此外,研究區(qū)厚壁菌門(Firmicutes)相對豐度僅為1%左右,與陜北沙化地區(qū)和內(nèi)蒙荒漠草原存在較大差異[23- 24],說明不同的植被能夠通過凋落物輸入及根系分泌物等差異改變土壤環(huán)境,進(jìn)而影響土壤細(xì)菌群落組成。

土壤細(xì)菌多樣性對提高土地生產(chǎn)力、增強(qiáng)脅迫環(huán)境下群落穩(wěn)定性有著重要影響。研究區(qū)各林地土壤細(xì)菌群落多樣性無顯著差異,說明隨著沙地樟子松人工林林齡的增加,土壤細(xì)菌群落結(jié)構(gòu)較為穩(wěn)定。同時(shí),人工林土壤表層細(xì)菌多樣性水平較高,這與其他學(xué)者的研究結(jié)果相似[25]。這是由于密集分布在土壤表層的樟子松根系及其分泌物能為細(xì)菌提供豐富的能源物質(zhì)[26]。

4.2 呼倫貝爾沙地樟子松人工林土壤細(xì)菌群落的影響因素

土壤細(xì)菌群落結(jié)構(gòu)與土壤理化性質(zhì)緊密相關(guān)且相互影響。隨著人工林林齡的增加,全磷含量先增加再降低,土壤有機(jī)質(zhì)和全氮含量逐漸下降,且林地全氮含量低于沙質(zhì)草地,這與前人研究存在差異[27]。呼倫貝爾沙地位于寒溫帶大陸性半干旱地區(qū),冬季漫長,植物生長期短,林地土壤表面養(yǎng)分含量與凋落物分解速率緊密相關(guān)。當(dāng)沙地樟子松為20—30齡時(shí),樟子松生長旺盛,凋落物量增加,土壤養(yǎng)分積累量增加,土壤養(yǎng)分含量升高[28]；而后,隨著林木生長,樟子松生長變緩,植物凋落物和根系分泌物的養(yǎng)分歸還量小于植物的養(yǎng)分吸收量,沙地樟子松人工林土壤養(yǎng)分庫存在消耗的現(xiàn)象[22,29]。呼倫貝爾沙地樟子松人工林土壤細(xì)菌群落的主要影響因子為速效鉀、全磷和全氮,與前人研究結(jié)果相似[30- 31],說明土壤細(xì)菌群落結(jié)構(gòu)與氮、磷、鉀等營養(yǎng)元素密切相關(guān),土壤限制性資源對土壤細(xì)菌相對豐度有顯著影響。25 a和34 a人工林土壤細(xì)菌的主要影響因子為全氮,可能是由于該地土壤氮資源與林地初級生產(chǎn)力緊密相關(guān)[32],兩林地內(nèi)全氮含量較高,充足的能源物質(zhì)也有利于土壤細(xì)菌的存活和繁衍；43 a人工林土壤細(xì)菌的主要影響因子為土壤含水量,該林地的土壤含水量顯著低于其他樣地,較低的土壤水分含量會限制微生物的代謝功能以及呼吸底物的擴(kuò)散,導(dǎo)致土壤呼吸作用減弱,并且水分脅迫還會通過改變根際環(huán)境來影響土壤細(xì)菌區(qū)系[33- 34]。此外,各林地0—10 cm土層細(xì)菌群落相對豐度與含水量呈正相關(guān),10—20 cm土層中兩者呈負(fù)相關(guān)。10—20 cm土層土壤含水量較低,僅為表層的0.55—0.74倍,進(jìn)一步說明了當(dāng)土壤中水分含量較低時(shí),土壤含水量會成為限制土壤細(xì)菌群落分布的主要因素。

4.3 呼倫貝爾沙地樟子松土壤細(xì)菌功能特征

植被對土壤細(xì)菌群落結(jié)構(gòu)的影響會進(jìn)一步影響土壤細(xì)菌功能組成[35]。呼倫貝爾沙地樟子松人工林土壤細(xì)菌主要涉及5個(gè)一級功能層和31個(gè)二級功能層。其中,土壤細(xì)菌代謝功能活躍且極為重要,如碳水化合物代謝與固氮、溶磷等作用密切相關(guān),有益于植物的氮、磷循環(huán)[36- 37]；氨基酸代謝可以幫助細(xì)菌吸收氨基酸,有利于細(xì)菌的存活與繁衍[38],部分土壤細(xì)菌類群還可以通過代謝產(chǎn)生抗生素、抑菌蛋白及促生長激素等,抑制病原菌并促進(jìn)植物生長[39]。43 a人工林代謝功能基因相對豐度顯著高于其他林地,說明沙地樟子松人工林能夠提高土壤細(xì)菌代謝活動,從而促進(jìn)自身對養(yǎng)分的吸收和利用。同時(shí),研究區(qū)土壤細(xì)菌活躍的膜運(yùn)輸功能有利于維持細(xì)菌各項(xiàng)生命活動,是土壤細(xì)菌正常發(fā)揮其功能的基礎(chǔ)[40]。沙質(zhì)草地的部分功能基因相對豐度顯著高于25 a和34 a人工林,這些功能的活躍也許能使土壤細(xì)菌在較為嚴(yán)苛的環(huán)境中更好地吸收養(yǎng)分。此外,土壤細(xì)菌功能隨不同樣地的變化均表現(xiàn)在0—10 cm土層,說明土層深度對細(xì)菌功能基因組成的影響較大,表層土壤細(xì)菌功能對環(huán)境變化更為敏感。

5 結(jié)論

呼倫貝爾沙地樟子松人工林土壤細(xì)菌隸屬35門92綱109目210科267屬,主要優(yōu)勢門為變形菌門、放線菌門和酸桿菌門。人工林變形菌門相對豐度顯著高于沙質(zhì)草地,酸桿菌門相對豐度顯著低于沙質(zhì)草地。不同林齡人工林土壤細(xì)菌群落多樣性無顯著差異。研究區(qū)沙地樟子松人工林土壤細(xì)菌群落結(jié)構(gòu)的主要影響因子為速效鉀、全磷和全氮。沙地樟子松人工林土壤細(xì)菌群落主要涉及5個(gè)一級功能層和31個(gè)二級功能層,43 a人工林土壤細(xì)菌代謝功能活躍,有利于促進(jìn)植物對養(yǎng)分的吸收和利用。沙地樟子松人工林的種植改善了土壤細(xì)菌群落結(jié)構(gòu),促進(jìn)土壤細(xì)菌代謝功能。研究成果有助于深入理解沙地樟子松與土壤細(xì)菌的互饋關(guān)系,并可為沙地樟子松人工林可持續(xù)經(jīng)營管理提供理論依據(jù)。