周柳婷,李建鵑,趙艷琳,羅 揚(yáng),白 瑩,陳 軍,吳則焰,3,4,*,林文雄

1 福建農(nóng)林大學(xué)生命科學(xué)學(xué)院,福州 350002 2 福建農(nóng)林大學(xué)林學(xué)院,福州 350002 3 作物生態(tài)與分子生理學(xué)福建省高校重點(diǎn)實(shí)驗(yàn)室,福州 350002 4 福建省農(nóng)業(yè)生態(tài)過(guò)程與安全監(jiān)控重點(diǎn)實(shí)驗(yàn)室,福州 350002

木麻黃(CasuarinaequisetifoliaForst.)是木麻黃屬常綠喬木,20世紀(jì)50年代作為庭院綠化樹(shù)種引入我國(guó)[1- 2],因其良好的適應(yīng)性及抗風(fēng)、耐貧瘠、耐鹽堿等生理特點(diǎn),在維持海岸帶生態(tài)系統(tǒng)穩(wěn)定、緩解沿海沙地鹽堿危害、防潮汐侵蝕等方面都起到極其重要的作用[3]。目前濱海沙地的適生樹(shù)種較少,木麻黃作為重要的防護(hù)林樹(shù)種,長(zhǎng)期受到嚴(yán)重的連栽障礙問(wèn)題困擾[4],具體表現(xiàn)為木麻黃平均木生物量、林分生物量及林分凈生產(chǎn)力均呈現(xiàn)逐代下降趨勢(shì),植株生長(zhǎng)緩慢,病蟲(chóng)害嚴(yán)重,造成產(chǎn)量、品質(zhì)下降[5- 6]。如何緩解或消減木麻黃連栽障礙,成為國(guó)內(nèi)外同行研究的前沿?zé)狳c(diǎn)。

前人研究認(rèn)為,植物連栽障礙是地力衰退[7]、化感物質(zhì)自毒作用[8- 9]、土壤微生態(tài)環(huán)境失衡等諸多方面綜合作用的結(jié)果。在地力衰退研究方面,有學(xué)者指出長(zhǎng)期連栽條件下植物對(duì)土壤礦質(zhì)營(yíng)養(yǎng)元素的片面吸收,使連栽土壤養(yǎng)分不均衡,進(jìn)而導(dǎo)致林分生產(chǎn)力下降、森林防護(hù)效益降低。更有甚者試圖通過(guò)增施肥料的措施來(lái)緩解連栽障礙但收效甚微,反而加劇了環(huán)境污染[10]。在化感物質(zhì)自毒作用研究方面,李鍵[11]、鄧蘭桂[12]等人從木麻黃小枝提取物中分離鑒定出了黃酮類(lèi)、鞣花酸類(lèi)化感物質(zhì),認(rèn)為自毒作用是引起木麻黃連栽障礙的原因。隨著研究的深入,國(guó)內(nèi)外同行逐漸認(rèn)識(shí)到化感物質(zhì)釋放到土壤后勢(shì)必受到微生物的加工、分解、轉(zhuǎn)化等[10,13],其存在只是誘因而不是導(dǎo)致植物連栽障礙形成的直接因素[14],因此土壤微生態(tài)成為連栽障礙研究的關(guān)注焦點(diǎn)。深入研究“植物-土壤-微生物”三者在根際的互作過(guò)程,對(duì)于探索土壤微生物與植物生長(zhǎng)發(fā)育之間的相互調(diào)控關(guān)系、揭示植物連栽障礙機(jī)制至關(guān)重要[15]。例如：Zhao等[16]研究表明,桉樹(shù)多代連栽后根際土壤中真菌多樣性增加；Wu等[17]發(fā)現(xiàn)連栽杉木使根際土壤中真菌/細(xì)菌比例顯著增加,且微生物群落多樣性和代謝活性顯著降低,從而打破了土壤微生態(tài)平衡；葉功富等[18]研究發(fā)現(xiàn),連栽木麻黃林地中細(xì)菌、放線(xiàn)菌、固氮菌與纖維素分解菌數(shù)量分別比一代林地減少22.46%、12.25%、32.66%和29.61%,真菌數(shù)量增加22.39%,微生物總量減少21.07%。目前,雖然已有少量關(guān)于連栽木麻黃土壤微生物數(shù)量變化的研究,但由于傳統(tǒng)微生物培養(yǎng)方法的局限性,未能充分揭示連栽條件下土壤微生物群落結(jié)構(gòu)及多樣性的變化規(guī)律[19]。本研究利用空間代替時(shí)間的方法,采用高通量測(cè)序技術(shù)擴(kuò)增細(xì)菌的16S rRNA基因高變區(qū),并對(duì)物種進(jìn)行注釋及相對(duì)豐度分析,探索連栽木麻黃根際土壤細(xì)菌群落變化規(guī)律,為揭示木麻黃連栽障礙機(jī)制提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)樣地選址于福建惠安赤湖國(guó)有防護(hù)林場(chǎng)(118°55′E,24°35′N(xiāo)),該林場(chǎng)占地面積約433 hm2,地處中國(guó)南亞熱帶氣候區(qū)。該地區(qū)夏季多臺(tái)風(fēng)和暴雨,秋冬季盛行東北風(fēng)。年均溫19.8℃,極端低溫1℃,極端高溫35℃,年降雨量1029 mm,年蒸發(fā)量2000 mm[20]。該林場(chǎng)栽植有3個(gè)不同代數(shù)的木麻黃人工林,其中第一代木麻黃林(First rotation plantation,FCP)栽植于1987年,第二代(Second rotation plantation,SCP)栽植于2011年,第三代(Third rotation plantation,TCP)栽植于2014年。不同樣地的伴生樹(shù)種和林下植被基本一致,伴生樹(shù)種均為少量的臺(tái)灣欒樹(shù)(Koelreuteriaelegans)和潺槁木姜子(Litseaglutinosa(Lour.)C. B. Rob.),林下植被均為霍香薊(AgeratumconyzoidesL.)和鬼針草(BidenspilosaL.)。土壤質(zhì)地均為風(fēng)積黃沙土。

1.2 土壤樣品的采集

2018年7月,于福建惠安赤湖境內(nèi)的國(guó)有防護(hù)林場(chǎng)內(nèi)設(shè)置樣地：在FCP、SCP和TCP設(shè)立3個(gè)林地概況相近的20 m×20 m實(shí)驗(yàn)樣地,每個(gè)樣地設(shè)3個(gè)重復(fù)樣方,共9個(gè)樣方。采集木麻黃根際土壤方法參考王圳等[21]的“抖落法”：沿“S”形路徑在每個(gè)樣方內(nèi)選擇20株平均胸徑、樹(shù)高相近的木麻黃,去除其落葉層后逐層挖去上層覆土,剪下細(xì)根分枝并輕輕抖動(dòng),仍粘在細(xì)根上的即為根際土壤,用小毛刷收集至自封袋中保存?zhèn)溆?。將每個(gè)樣方內(nèi)取得的20份根際土壤混勻?yàn)?份土樣,共獲得9份土壤樣品。土樣通過(guò)2 mm篩選后其中一部分保存于-80℃冰箱,另一部分土樣4℃冰箱保存用于測(cè)定土壤中微生物多樣性。

1.3 根際土壤DNA提取及PCR擴(kuò)增

利用試劑盒BioFast Soil Genomic DNA Extraction Kit(BioFlux公司,中國(guó))提取木麻黃根際土壤微生物基因組總DNA。DNA濃度通過(guò)Nanodrop進(jìn)行測(cè)定,并用1.2%瓊脂糖凝膠電泳檢測(cè)DNA 的質(zhì)量。再將基因組總DNA用無(wú)菌水稀釋至1 ng/μL。以稀釋后的基因組DNA為模板,采用 Phusion? High-Fidelity PCR Master Mix with GC Buffer和高效高保真酶(New England Biolabs公司,美國(guó))對(duì)細(xì)菌16S rRNA基因的V4區(qū)域進(jìn)行擴(kuò)增。擴(kuò)增引物為515F(5′-GTGCCAGCMGCCGCGGTAA- 3′)和806R(5′-GGACTACHVGGGTWTCTAAT- 3′)[22]。采用50 μL PCR擴(kuò)增體系,反應(yīng)程序?yàn)椋?8℃預(yù)變性1 min；30個(gè)循環(huán)包括(98℃,10 sec；50℃,30 sec；72℃,30 sec)；72℃,5 min。擴(kuò)增后用2%瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物的質(zhì)量。將PCR產(chǎn)物送至北京諾禾至源科技有限公司進(jìn)行高通量測(cè)序。

1.4 數(shù)據(jù)分析

測(cè)序所得的原始數(shù)據(jù)經(jīng)Fast QC軟件進(jìn)行質(zhì)控后,利用Cutadapt(V1.9.1)軟件去除短序列(<200 bp)及低質(zhì)量序列(q<25)[23]。通過(guò)Mothur方法與SILVA數(shù)據(jù)庫(kù)對(duì)97%相似度的有效序列分配操作分類(lèi)單元(Operational Taxonomic Units,OTUs)[24]。依據(jù)物種分類(lèi)信息繪制物種分類(lèi)條形圖和物種豐度熱圖[25]。再通過(guò)QI-IME 軟件計(jì)算樣品的覆蓋度和多樣性指標(biāo)[26- 28]。運(yùn)用PICRUSt軟件對(duì)菌群的功能和代謝途徑進(jìn)行預(yù)測(cè)分析[29]。

2 結(jié)果與分析

2.1 連栽木麻黃根際土壤細(xì)菌群落OTU組成及結(jié)構(gòu)

通過(guò)IonS5TMXL測(cè)序平臺(tái),從3個(gè)林地概況相近的不同連栽代數(shù)木麻黃根際土壤中總共得到338560條有效序列,平均每個(gè)樣品的有效序列為37618條。以97%的序列相似性閾值將有效序列聚類(lèi)為17627個(gè)OTU。其中FCP均值為1905個(gè)OTU,SCP均值為2046個(gè)OTU,TCP均值為1924個(gè)OTU。平均而言,約98.6%的有效序列可聚類(lèi)為門(mén)級(jí),但僅有38.68%以上的有效序列可聚類(lèi)為屬級(jí)。稀釋曲線(xiàn)可直接反映所抽取的優(yōu)化序列深度是否合理,并間接反映樣品中的物種豐富度[30]。由圖1可知,抽取的序列條數(shù)達(dá)到25000條以上,曲線(xiàn)趨于平坦,表明不同處理土壤所測(cè)的序列庫(kù)容都能夠較好地反映細(xì)菌群落種類(lèi)數(shù)量,測(cè)序數(shù)據(jù)量基本合理。

Venn圖可以直觀(guān)地體現(xiàn)不同代數(shù)木麻黃根際土壤細(xì)菌群落OTU組成的差異性及重疊情況(圖2)。維恩分析結(jié)果表明,在OTU水平上,FCP中特異性細(xì)菌OTU占總OTU序列數(shù)的9.86%(341),SCP中特異性細(xì)菌OTU占14.68%(508),TCP中特異性細(xì)菌OTU占9.62%(333)。此外,FCP與SCP共有的OTU數(shù)量為292(8.44%),SCP與TCP共有的OTU數(shù)量為345(9.97%),FCP與TCP共有的OTU數(shù)量為177(5.12%)。FCP、SCP、TCP根際土壤中共有的細(xì)菌OTU數(shù)量為1464(42.31%)。

圖1 木麻黃根際土壤樣本微生物群落的稀釋曲線(xiàn)Fig.1 Rarefaction curves of microbial communities in rhizosphere soils samples of C. equisetifolia FCP：第一代First rotation plantation；SCP：第二代Second rotation plantation；TCP：第三代Third rotation plantation

圖2 基于OTU豐度的不同代數(shù)木麻黃根際土壤細(xì)菌群落維恩圖Fig.2 Venn diagram of bacterial communities in rhizosphere soil of different rotations of C. equisetifolia based on OTUs abundance

2.2 連栽木麻黃根際土壤細(xì)菌群落結(jié)構(gòu)

2.2.1連栽木麻黃根際土壤細(xì)菌群落多樣性分析

通過(guò)土壤樣品中細(xì)菌群落多樣性分析指數(shù),對(duì)連栽木麻黃根際土壤細(xì)菌群落物種豐富度和均勻度進(jìn)行評(píng)估(cutoff=27075)。由表1可知,隨著連栽代數(shù)增加,Chao1指數(shù)、ACE指數(shù)均逐代下降,而Observed species、Shannon指數(shù)則呈現(xiàn)先上升后降低的趨勢(shì),但差異均不顯著(P>0.05),Simpson指數(shù)幾乎沒(méi)有發(fā)生變化。Shannon指數(shù)分析顯示,木麻黃根際土壤SCP(8.874)中細(xì)菌多樣性相對(duì)較高,FCP(8.780)最低。從Chao1指數(shù)分析來(lái)看,FCP物種豐富度最高(2086.012),SCP次之(2005.244),TCP最低(2002.447)。

表1 連栽木麻黃根際土壤細(xì)菌群落豐富度和多樣性指數(shù)

每列不同字母表示差異達(dá)顯著水平(P≤0.05,n=3)

圖3 不同代數(shù)木麻黃根際土壤樣本間的Weighted Unifrac距離和Unweighted Unifrac距離Fig.3 Weighted and unweighted Unifrac distances between rhizosphere soils samples of different rotations of C. equisetifolia

通過(guò)Weighted Unifrac距離和Unweighted Unifrac距離來(lái)比較兩個(gè)樣品間物種多樣性的相異程度。FCP和SCP之間的Weighted Unifrac距離和Unweighted Unifrac距離分別為0.181和0.416,FCP和TCP之間分別為0.174和0.405,SCP和TCP之間分別為0.123和0.391。

圖4 連栽木麻黃根際土壤細(xì)菌群落PcoA分析Fig.4 PcoA analysis of bacterial communities from rhizosphere soils of continuous rotations of C. equisetifolia

2.2.2連栽木麻黃根際土壤細(xì)菌群落PCoA分析與UPGMA聚類(lèi)分析

基于OTU的木麻黃根際土壤細(xì)菌群落主坐標(biāo)分析(Principal Co-ordinates Analysis,PCoA)結(jié)果見(jiàn)圖4,主成分1(PC1)與主成分2(PC2)分別解釋變量方差的30.32%、22.93%,兩者累計(jì)貢獻(xiàn)率達(dá)53.25%。PC1將FCP中的細(xì)菌群落與SCP、TCP明顯區(qū)分開(kāi),PC2將SCP中的細(xì)菌群落與FCP、TCP明顯區(qū)分開(kāi)。UPGMA(Unweighted Pair-group Method with Arithmetic Mean)聚類(lèi)分析結(jié)果見(jiàn)圖5,SCP和TCP的細(xì)菌群落結(jié)構(gòu)聚集為一個(gè)群體,再與FCP聚集在一起形成系統(tǒng)發(fā)生樹(shù),說(shuō)明FCP與SCP、TCP二者存在較大差異。

2.2.3連栽木麻黃根際土壤細(xì)菌群落組成及結(jié)構(gòu)變化

圖5 連栽木麻黃根際土壤的UPGMA聚類(lèi)分析Fig.5 UPGMA analysis of bacterial communities from rhizosphere soil of continuous rotations of C. equisetifolia

在木麻黃根際土壤細(xì)菌群落中共檢測(cè)到10個(gè)門(mén),優(yōu)勢(shì)類(lèi)群依次是變形菌門(mén)(Proteobacteria)(30.27%—34.42%)、酸桿菌門(mén)(Acidobacteria)(28.10%—32.21%)、放線(xiàn)菌門(mén)(Actinobacteria)(14.86%—18.84%)、疣微菌門(mén)(Verrucomicrobia)(3.79%—8.09%)(圖6)。由表2可知,在屬水平上,熱酸菌屬(Acidothermus)、CandidatusSolibacter、慢生根瘤菌屬(Bradyrhizobium)、根微菌屬(Rhizomicrobium)、酸桿菌屬(Acidibacter)、變桿菌屬(Variibacter)、布氏桿菌屬(Bryobacter)是根際土壤細(xì)菌群落的優(yōu)勢(shì)屬(豐度>1%),其中FCP中最高豐度屬是熱酸菌屬(Acidothermus),而SCP、TCP根際土壤中CandidatusSolibacter相對(duì)豐度最高。隨著木麻黃連栽代數(shù)增加,根際土壤中熱酸菌屬(Acidothermus)、根微菌屬(Rhizomicrobium)、酸桿菌屬(Acidibacter)、布氏桿菌屬(Bryobacter)的相對(duì)豐度大體呈下降趨勢(shì),其中酸桿菌屬(Acidibacter)的相對(duì)豐度在FCP、SCP、TCP之間的下降趨勢(shì)達(dá)顯著水平(P<0.05)。與之相反,慢生根瘤菌屬(Bradyrhizobium)的相對(duì)豐度隨連栽代數(shù)的增加而增加。此外,從木麻黃根際土壤中前35個(gè)優(yōu)勢(shì)細(xì)菌在屬水平的群落熱圖分析也可以看出,木麻黃連栽后土壤細(xì)菌群落組成結(jié)構(gòu)發(fā)生了顯著變化(圖7)。

圖6 連栽木麻黃根際土壤優(yōu)勢(shì)細(xì)菌在門(mén)水平的群落相對(duì)豐度Fig.6 The relative abundance of predominant bacteriain rhizospheric soil samples of continuous rotations of C. equisetifolia at the phylum level

3 結(jié)論與討論

植物根系與根際微生物形成的植物-微生物群落綜合體,是“植物-土壤-微生物”三者進(jìn)行根際互作過(guò)程的重要場(chǎng)所[31]。根際微生物是土壤生態(tài)系統(tǒng)能量流動(dòng)、物質(zhì)循環(huán)和信息傳遞最主要的驅(qū)動(dòng)者,對(duì)維持土壤生產(chǎn)力至關(guān)重要[32]。本研究通過(guò)高通量測(cè)序技術(shù)對(duì)連栽木麻黃根際土壤細(xì)菌群落變化進(jìn)行探究。結(jié)果表明連栽后根際土壤細(xì)菌群落的組成與結(jié)構(gòu)發(fā)生了顯著變化。由多樣性分析可知,木麻黃根際土壤細(xì)菌群落的Simpson、Chao1、ACE指數(shù)隨著連栽代數(shù)的增加而下降,Observed species、Shannon指數(shù)呈現(xiàn)先增加后減少的趨勢(shì),變化趨勢(shì)不一致的原因可能是由于土壤理化性質(zhì)不完全一致、林齡差異、森林凋落物等因素導(dǎo)致。這與韓亞飛[33]研究楊樹(shù)土壤細(xì)菌群落多樣性所提出的觀(guān)點(diǎn)相似,但與李延茂等[34]通過(guò)DGGE分析法研究杉木根際土壤細(xì)菌群落的結(jié)果不一致。推測(cè)可能是由于技術(shù)方法不同、樹(shù)種差異、地域差別、土壤異質(zhì)性等因素造成研究結(jié)果不一致。

進(jìn)一步分析連栽木麻黃根際土壤細(xì)菌群落組成與結(jié)構(gòu)變化。由圖6可知,不同連栽代數(shù)木麻黃根際土壤細(xì)菌在門(mén)水平上的最優(yōu)勢(shì)類(lèi)群是變形菌門(mén)(Proteobacteria),分別占FCP、SCP、TCP根際土壤總細(xì)菌豐度的34.42%、30.27%、32.34%,證實(shí)了變形菌門(mén)是陸地土壤生態(tài)系統(tǒng)細(xì)菌群落的優(yōu)勢(shì)類(lèi)群,與Bazylinski等[35]的研究結(jié)果一致。慢生根瘤菌屬(Bradyrhizobium)是變形菌門(mén)α變形菌綱的一種根瘤菌[36],可以與寄主植物形成共生關(guān)系,將大氣中的游離氮固定成寄主生物可以利用的形式,例如氨(NH3)或銨(NH4+)[37]。木麻黃是重要的結(jié)瘤共生固氮樹(shù)種之一[38],可以富集根瘤菌以適應(yīng)貧瘠的濱海沙地環(huán)境。在本研究中,隨著連栽代數(shù)的增加,慢生根瘤菌屬(Bradyrhizobium)的相對(duì)豐度隨之增加約132.57%、244.81%,且FCP、SCP兩者與TCP的差異達(dá)顯著水平(P<0.05)。與朱琳等[39]研究大豆連作土壤細(xì)菌群落多樣性所得的結(jié)論一致。與之相反,根際土壤中的熱酸菌屬(Acidothermus)、根微菌屬(Rhizomicrobium)、酸桿菌屬(Acidibacter)、布氏桿菌(Bryobacter)的相對(duì)豐度呈下降趨勢(shì)。變形菌門(mén)下的根微菌屬(Rhizomicrobium)是在土壤氮循環(huán)、解磷過(guò)程以及降解有機(jī)物過(guò)程中發(fā)揮重要作用的有益微生物[40]。在本研究中,隨著連栽代數(shù)的增加,根微菌屬(Rhizomicrobium)的相對(duì)豐度隨之降低約311.49%、282.16%,且FCP與SCP、TCP兩者的差異達(dá)顯著水平(P<0.05)。酸桿菌門(mén)下的酸桿菌屬(Acidibacter)是具有分解蛋白質(zhì)和攝取環(huán)境周?chē)嵝晕镔|(zhì)的有益菌屬[41]。在本研究中,隨著連栽代數(shù)的增加,酸桿菌屬(Acidibacter)的相對(duì)豐度隨之降低約176.07%、284.54%,且FCP、SCP與TCP三者之間的差異達(dá)顯著水平(P<0.05)。由此可見(jiàn),木麻黃連栽后根際土壤微生物群落結(jié)構(gòu)發(fā)生了改變。這種現(xiàn)象可能是由于根系化感物質(zhì)所產(chǎn)生的間接生態(tài)效應(yīng)引起的,從而導(dǎo)致根際微生物發(fā)生惡化偏移。化感物質(zhì)可以直接或間接影響根際微生物的代謝和生長(zhǎng)發(fā)育,進(jìn)而對(duì)根際微生物的種類(lèi)、數(shù)量和分布產(chǎn)生影響[42]。如田給林[43]研究發(fā)現(xiàn),外源添加阿魏酸能顯著促進(jìn)草莓專(zhuān)化型尖孢鐮刀菌的生長(zhǎng)及孢子萌發(fā),從而間接侵染草莓植物,造成連作障礙的發(fā)生。劉金光等[44]研究發(fā)現(xiàn),低濃度的苯甲酸對(duì)曲霉菌屬和鐮刀菌屬有顯著促進(jìn)作用,造成根際土壤微生物群落結(jié)構(gòu)失衡。

本研究證實(shí)了木麻黃長(zhǎng)期單一化栽植后,根際土壤細(xì)菌群落結(jié)構(gòu)發(fā)生了顯著變化,導(dǎo)致土壤微生態(tài)失衡,間接造成木麻黃連栽障礙問(wèn)題。然而,連栽障礙形成及加重的原因不是單一孤立的,是眾多因素彼此影響、相互制約的結(jié)果。唯有建立在充分了解“土壤-微生物-根際”相互關(guān)系的基礎(chǔ)上,才能系統(tǒng)地認(rèn)識(shí)連栽木麻黃的根際生態(tài)效應(yīng)。因此,對(duì)于木麻黃根際土壤化感物質(zhì)如何介導(dǎo)土壤微生物區(qū)系定向演變？以及關(guān)鍵的特異微生物如何在木麻黃連栽障礙形成的根際生態(tài)學(xué)中扮演重要角色,仍需實(shí)驗(yàn)驗(yàn)證討論。

圖7 連栽木麻黃根際土壤優(yōu)勢(shì)細(xì)菌在屬水平的群落熱圖分析Fig.7 Heat map analysis of the main bacteria in rhizospheric soil of continuous rotations of C. equisetifolia at the genus level

表2 連栽木麻黃根際土壤在屬水平的細(xì)菌群落信息

Table 2 Information regarding the bacterial community at the genus level inrhizospheric soil ofC.equisetifoliaunder successive rotations

屬Genus相對(duì)豐度Relative abundance/%FCPSCPTCP熱酸菌屬Acidothermus5.2841a2.498b3.5629bCandidatus Solibacter3.7661a3.7932a3.9471a慢生根瘤菌屬Bradyrhizobium1.6103b2.1348b3.9421a根微菌屬Rhizomicrobium3.4669a1.1130b1.2287b酸桿菌屬Acidibacter3.6959a2.0991b1.2989c變桿菌屬Variibacter3.6171a2.2727b3.4595ab布氏桿菌屬Bryobacter2.5793a1.5328b1.6707b伯克氏菌屬Burkholderia0.8815b1.3777a1.7593aH161.2533a0.7842b0.7707bTerracidiphilus0.7436b0.9037b1.2213a放線(xiàn)菌屬Actinospica0.1404b0.4481b0.8593a克洛氏菌屬Crossiella0.2425c0.4087b0.9701a地桿菌屬Geobacter0.2622b0.5134ab0.8199a分枝桿菌屬M(fèi)ycobacterium0.4986b0.5368b0.8384a假羅布里斯屬Pseudolabrys0.3804ab0.5405a0.3053b鏈霉菌屬Streptomyces0.2278c0.474b0.6316a埃達(dá)布特屬Edaphobacter0.4666a0.0246b0.0320b乳桿菌屬Acidobacterium0.5392a0.0813b0.1305b苯基桿菌屬Phenylobacterium0.4555a0.1108b0.0874b酸雙酯屬Acidicaldus0.4925a0.3533b0.2979b梭菌屬Clostridium sensu stricto 10.3029b0.4112a0.3977abMizugakiibacter0.3817a0.0111b0.0037b

每列不同字母表示差異達(dá)顯著水平(P<0.05,n=3)