張 玥，姜愛霞，郭篤發(fā)

(山東師范大學人口·資源與環(huán)境學院，山東濟南 250014)

基于16S rDNA基因文庫的黃河三角洲鹽生植被土壤細菌群落多樣性研究

張 玥，姜愛霞，郭篤發(fā)*

(山東師范大學人口·資源與環(huán)境學院，山東濟南 250014)

[目的]研究黃河三角洲光板地及2種鹽生植被(獐茅和羅布麻)下土壤細菌群落結(jié)構(gòu)多樣性，探究其與植被演替的關(guān)系。[方法]采用細菌16S rDNA基因文庫方法，各文庫挑選200個陽性克隆子進行測定。[結(jié)果]光板地、獐茅和羅布麻3個文庫中分別得到121、150、155條有效序列。獐茅土壤細菌的Shannon指數(shù)和ACE指數(shù)最高。土壤中檢測到變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、酸桿菌門(Acidobacteria)等共8門。變形菌門在光板地、獐茅土壤和羅布麻土壤中的相對豐度分別為37.45%、51.09%、37.11%，擬桿菌門分別為33.02%、3.03%、4.47%，其他細菌相對豐度均較低。[結(jié)論]3種覆被類型下土壤細菌在種群組成與分布上差異明顯，變形菌為優(yōu)勢菌群。當鹽生植被處于不同演替階段時，土壤細菌群落結(jié)構(gòu)差別較大。

16S rDNA；黃河三角洲；細菌群落多樣性；鹽生植被

土壤微生物是生態(tài)系統(tǒng)的重要組成部分，其對土壤中幾乎所有過程都起著直接或間接的作用。在土壤微生物中，細菌的數(shù)量和種類最多[1]，其群落結(jié)構(gòu)多樣性的變化對土壤肥力、植被生長和污染物降解等具有重要影響，是評價自然及人為因素引起的土壤質(zhì)量變化的重要指示因子[2]。

土壤細菌起初被認為呈全球隨機性分布[3-4]，后來的研究將多種影響土壤細菌群落結(jié)構(gòu)多樣性的因素歸納為自然和人為兩大類[5]，自然因素中的植物群落演替與土壤細菌群落結(jié)構(gòu)的關(guān)系逐漸成為研究熱點。研究表明，微生物群落結(jié)構(gòu)變化與植被出現(xiàn)與否無關(guān)，表明植被對微生物的作用可能還與植被類型或植被演替階段有關(guān)[6]。盡管植被群落與土壤微生物關(guān)系的研究取得了較大進展，但在濱海鹽漬環(huán)境中，土壤細菌群落多樣性的變化與鹽生植被演替的關(guān)系尚不明晰。

筆者采用近年來在細菌多樣性領(lǐng)域應用廣泛的16S rDNA克隆文庫技術(shù)，研究黃河三角洲鹽生植被土壤細菌群落多樣性，探究其與植被演替的關(guān)系，為進一步研究黃河三角洲自然保護區(qū)濕地生態(tài)系統(tǒng)提供理論依據(jù)。

1 材料與方法

1.1 土樣采集在黃河三角洲研究區(qū)內(nèi)沿東西向路線(橫跨現(xiàn)代三角洲和近代三角洲)，依照鹽生植被正向演替順序，于2012年6月選擇光板地、一種強耐鹽植被獐茅、一種輕耐鹽植被羅布麻，每個植被群落取3個樣地，每個樣地采用對角線五點采樣法采集土壤樣品。

采集0～20 cm深的土壤樣品，并混合均勻，除去較大的根系等雜物，分成2份，每份約200 g。一份土樣封于滅菌袋中，保存在4 ℃條件下，用于微生物培養(yǎng)和土壤理化性質(zhì)分析；另一份土樣放入液氮罐帶回實驗室，保存在-80 ℃條件下，用于分子生物學研究。

1.2 DNA提取使用Omega DNA提取試劑盒對土壤樣品DNA進行提取，具體操作步驟見試劑盒操作說明書。

1.3 細菌16S rDNA 的PCR擴增PCR擴增所用引物為：27F(5′-AGAGTTTGATCCTGGCTCAG-3′)、1525R (5′-AGA-AAGGAGGTGATCCAGCC-3′)。PCR反應體系的各組分組成為：引物(10 pmol/L)各0.4 μl ，DNA模版1.0 ng，MgCl2(2.5 mmol/L) 0.3 μl，10×Buffer 2.0 μl，dNTP(2.5 mmol/L)0.4 μl，DNA聚合酶0.2 μl，Milli-Q Water補至20 μl。PCR擴增反應條件：94 ℃預變性5 min；94 ℃變性30 s,55 ℃退火1 min,72 ℃延伸2 min，30個循環(huán)；72 ℃延伸10 min。PCR產(chǎn)物采用熒光染料作為標記，用0.8%(W/V)的瓊脂糖凝膠在1×TAE的緩沖液中進行電泳測定。

1.4 擴增產(chǎn)物的克隆和測序使用試劑盒(Tiangen Biotech Beijing)回收16S rDNA片斷并進行純化。將純化產(chǎn)物連接至pMD 19-T Vector載體上，再轉(zhuǎn)入大腸桿菌(Escherichiacoli)感受態(tài)細胞中。在含有氨芐青霉素LB培養(yǎng)基平板上分別涂布，37 ℃下培養(yǎng)過夜。利用T7和SP6進行菌落PCR。PCR擴增后檢驗，挑選插入目的片斷的陽性克隆子，送山東農(nóng)業(yè)科學院生物信息技術(shù)服務中心測序。

1.5 序列分析及發(fā)育樹的構(gòu)建將測序所得序列在BLASTN中進行同源性比對，挑選與序列親緣關(guān)系最相似序列，使用CHECK-CHIMERA軟件分析檢查，去除嵌合及怪異序列。利用ClustalW2軟件將測定的序列和部分網(wǎng)上下載的16S rDNA相似序列進行多重分析比對，找出可操作分類單元(OTU)。使用Mega5.0軟件中的鄰接法建構(gòu)系統(tǒng)發(fā)育樹。

1.6 多樣性指數(shù)計算使用Mothur軟件進行土壤細菌α多樣性分析，計算α多樣性指數(shù)，包括豐富度指數(shù)(Richness)、覆蓋率(Coverage)、香農(nóng)指數(shù)(Shannon)、ACE指數(shù)、Chao1指數(shù)。

2 結(jié)果與分析

2.1 土壤理化性質(zhì)黃河三角洲土壤可溶性鹽分主要為NaCl，土壤電導率與鹽分含量能用線性進行模擬[7]，因此可用土壤電導率表征土壤含鹽量。由表1可知，3種鹽生植被下土壤1∶5土水比浸提液的電導率在0.04～1.95 dS/m。光板地、獐茅地為海積潮濕正常鹽成土，其中光板地鹽化程度比獐茅地高很多，羅布麻地為弱鹽淡色潮濕雛形土；根據(jù)國際土壤質(zhì)地分類制，光板地屬于粉砂壤土，獐茅地屬于粉砂黏壤土。有機質(zhì)、全氮、堿解氮在獐茅和羅布麻土壤中的含量均高于光板地，羅布麻地速效磷含量高于光板地和獐茅地，差異明顯。由此可知，土壤速效磷隨鹽生植被正向演替有不斷升高的趨勢。光板地的溫度高于其他覆被類型。

表1 不同覆被類型土壤理化性質(zhì)

2.2 土壤細菌測序深度及多樣性分析由圖1可知，OTU數(shù)目隨序列數(shù)增多而增加，但增幅逐漸變小，曲線趨于平緩，但仍未達到飽和，表明挑取的克隆數(shù)能基本反映土壤樣品中細菌群落的結(jié)構(gòu)組成。由表2可知，3種覆被類型下土壤細菌覆蓋率在52.67%～61.67%，表明黃河三角洲土壤細菌具有較高的多樣性。獐茅土壤Shannon指數(shù)最高，為8.52，說明獐茅土壤細菌具有較高的多樣性，光板地Shannon指數(shù)最低，為7.19。ACE指數(shù)和Chao1指數(shù)在3種土壤中差異明顯，獐茅土壤ACE指數(shù)最高，為31 554.68，獐茅土壤ACE指數(shù)和Chao1指數(shù)均最高，分別為3 155.68和18 806.41。相關(guān)性統(tǒng)計分析表明，土壤理化性質(zhì)與各多樣性指數(shù)無顯著相關(guān)關(guān)系。

表2 不同土壤細菌多樣性指數(shù)

2.3 土壤細菌16S rDNA序列分析通過Mothur分析，將相似性大于97%的序列歸為同一種OTU，得到有效序列和OTU分析結(jié)果：光板地有效序列121個，OTU 72個；獐茅有效序列150個，OTU 83個；羅布麻有效序列155個，OTU 88個。

對系統(tǒng)進化樹進行分析，發(fā)現(xiàn)細菌序列主要歸屬于8個門類：變形菌門(Proteobacteria)，擬桿菌門(Bacteroidetes)，異常球菌-棲熱菌門(Deinococcus-Thermus)，酸桿菌門(Acidobacteria)，綠彎菌門(Chloroflexi)，厚壁菌門(Firmicutes)，疣微菌門(Verrucomicrobia)和放線菌門(Actinobacteria)。8種土壤細菌門在3種土壤中的分布情況見圖2。由圖2可知，變形菌門在3種土壤中相對豐度為37.11%～51.09%，為優(yōu)勢菌群，經(jīng)相關(guān)性分析，δ-變形菌(δ-Proteobacteria)與土壤全氮相關(guān)系數(shù)為0.89 (P<0.05)，顯著相關(guān)；擬桿菌門相對豐度次于變形菌門，且在3種土壤中的分布差異明顯，為3.03%～33.02%，光板地中相對豐度最高，經(jīng)相關(guān)性分析，擬桿菌門與土壤電導率(含鹽量)相關(guān)系數(shù)為0.99 (P<0.01)，極顯著相關(guān)；酸桿菌門在3種土壤中相對豐度為0.21%～23.87%，在羅布麻土壤中相對豐度最高；放線菌門在3種土壤中相對豐度為2.32%～5.90%，在獐茅土壤中相對豐度最高；厚壁菌門在3種土壤中相對豐度為1.43%～4.22%，在羅布麻土壤中相對豐度最高，經(jīng)相關(guān)性分析，厚壁菌門與土壤含水率相關(guān)系數(shù)為0.94 (P<0.05)，顯著相關(guān)；綠彎菌門在3種土壤中相對豐度為0.81%～3.43%，在獐茅土壤中相對豐度最高；疣微菌門在3種土壤中相對豐度為0.37%～1.59%，在羅布麻土壤中相對豐度最高；異常球菌-棲熱菌門只在光板地中檢測到，相對豐度為0.23%。

2.3.1變形菌門。在獐茅土壤中相對豐度最高，為51.09%，有40個OTU。α-變形菌(α-Proteobacteria)含量最高，在獐茅土壤中相對豐度為14.67%，有17個OTU，其中12個屬于Rhodospirillales，3個屬于Rhizobiales，2個屬于Rhodospirillales；γ-變形菌(γ-Proteobacteria)相對豐度次之，為13.33%，有7個OTU，3個屬于綠膿桿菌(Pseudomonas)，1個屬于硫發(fā)菌目(Thiotrichales)，1個屬于海螺菌目(Oceanospirillales)，1個屬于交替單胞菌目(Alteromonadales)；β-變形菌(β-Proteobacteria)相對豐度為2.67%，有4個OTU，均屬于伯克氏菌目(Burkholderiales)；δ-變形菌(δ-Proteobacteria) 相對豐度最低，為2.00%，有3個OTU，1個OTU屬于互營桿菌目(Syntrophobacterales)，1個屬于脫硫菌目(Desulfobacterales)，1個OTU屬于粘球菌目(Myxococcales)。

在光板地中相對豐度次之，為37.45%，有16個OTU。δ-變形菌含量最高，在光板地中相對豐度為5.79%，有7個OTU，其中6個OTU屬于脫硫菌目，1個屬于粘球菌目；γ-變形菌相對豐度次之，為4.96%，有5個OTU；α-變形菌相對豐度為1.65%；β-變形菌相對豐度為1.65%，有1個OTU，均屬于伯克氏菌目叢毛單胞菌科叢毛單胞菌屬(Comamonas)。

在羅布麻土壤中相對豐度最低，為37.11%，有13個OTU。γ-變形菌在羅布麻土壤中相對豐度最高，為9.68%，有4個OTU；δ-變形菌相對豐度次之，為3.23%，有5個OTU，全部屬于粘球菌目，β-變形菌相對豐度最低，為1.94%，有3個OTU，均屬于伯克氏菌目。

2.3.2擬桿菌門。擬桿菌門在光板地中相對豐度最高，為32.02%。有17個克隆體為擬桿菌綱(Bacteroidetes)，10個克隆體為海洋滑行菌屬(Marinilabilia)，4個克隆體屬于鞘脂桿菌綱(Sphingobacteria)；在羅布麻土壤中相對豐度次之，為4.47%，有13個OTU，全部為鞘脂桿菌綱；在獐茅土壤中相對豐度最低，為3.03%，有3個OTU，全部為鞘脂桿菌綱。

2.3.3酸桿菌門。在羅布麻土壤中相對豐度最高，為16.65%，有30個OTU，其中包含Gp4 6個OTU，Gp9 3個OTU，Gp5 3個OTU，Gp7 2個OTU，Gp26、Gp17、Gp1、Gp3各1個OTU；在獐茅地土壤中相對豐度次之，為7.01%，有5個OTU，分別屬于Gp26、Gp21、Gp4、Gp10；在光板地中相對豐度最低，只有0.21%，有2個OTU，分別屬于Gp3、Gp10。經(jīng)相關(guān)性分析，Gp7、Gp15與土壤有效磷相關(guān)系數(shù)均為0.90 (P<0.05)，顯著相關(guān)。

2.4 土壤細菌群落結(jié)構(gòu)比較由圖3可知，光板地與屬于重鹽土壤的獐茅土壤細菌及屬于輕鹽土壤的羅布麻土壤細菌間共有的OTU序列為16個，分別占三者OTU總數(shù)的66.67%、53.33%和80.00%；光板地與屬于重鹽土壤的獐茅土壤細菌間共有的OTU序列為5個，分別占兩者OTU總數(shù)的20.83%和16.67%；光板地與屬于輕鹽土壤的羅布麻土壤細菌間共有的OTU序列為0個；屬于重鹽土壤的獐茅土壤與屬于輕鹽土壤的羅布麻土壤細菌間共有的OTU序列為3個，分別占兩者OTU總數(shù)的10.00%和15.00%。由此可知，當鹽生植被在不同演替階段時，土壤細菌群落結(jié)構(gòu)差別較大。

3 結(jié)論與討論

該研究利用細菌16S rDNA基因文庫克隆測序技術(shù)，對黃河三角洲不同演替階段的代表性鹽生植被：光板地、獐茅和羅布麻的土壤細菌群落結(jié)構(gòu)多樣性進行分析，進而明確其與植被演替的關(guān)系。變形菌門在3種植被土壤中為優(yōu)勢菌群，在屬于重鹽土壤的獐茅土壤中相對豐度最大，這與以往利用細菌16S rDNA克隆文庫方法對土壤樣品進行研究所得結(jié)論基本一致。Tait等[8]對美國馬薩諸塞州礁湖的研究、Schloss等[9]對農(nóng)田土壤樣品細菌的研究、Axelrood等[10]對英國哥倫比亞森林土壤的研究，均發(fā)現(xiàn)變形菌門為土壤中的優(yōu)勢菌群。擬桿菌門在3種土壤中相對豐度次于變形菌門，在光板地中相對豐度最高，其余2種土壤中相對豐度均低，且與土壤電導率(含鹽量)呈極顯著相關(guān)關(guān)系，表明其分布可能與土壤鹽分有關(guān)。

酸桿菌門是最近被劃分出的一門細菌，多為嗜酸菌，目前對其研究較少，但已有研究表明酸桿菌門廣泛存在于土壤中，且為優(yōu)勢菌種[10]。酸桿菌門主要分布于陸地、海洋沉積物和活性淤泥中，酸性土壤有利于其生長，該研究中，酸桿菌門在光板地、重鹽土壤(獐茅土壤)中的相對豐度明顯小于其在輕鹽土壤(羅布麻土壤)中。關(guān)于酸桿菌門中的Gp7、Gp15與土壤有效磷顯著相關(guān)的原因，有待進一步研究。

異常球菌-棲熱菌門包括一些能抵抗嚴酷環(huán)境的球狀細菌，有較強的耐熱、抗輻射、抗砷能力[11]。該研究只在光板地中檢測到異常球菌-棲熱菌門的出現(xiàn)，占有較大比重，且與土壤電導率(含鹽量)呈極顯著相關(guān)關(guān)系，這或許說明其也有高強度的耐鹽性。Saccharospirillum在輕鹽土壤(羅布麻土壤)中豐度大于其他類型土壤，今后可對其進行深入研究，探討將其作為黃河三角洲脫鹽演化指示性細菌的可能性。

許多細菌類群的分布并未發(fā)現(xiàn)有規(guī)律的演替變化，如彎菌門和疣微菌門，呈現(xiàn)一定的隨機性，它們的分布除受土壤含鹽量影響外，可能還與土壤植被類型以及人為影響等多種因素的共同作用有關(guān)[12]。

目前，對環(huán)境樣品16S rDNA克隆文庫的分析，通常隨機選取50～100個克隆子探討其細菌類群和多樣性[13]。由于克隆文庫較小，不能包含樣品所有細菌的種類，難以獲得整個細菌群落的所有信息。該研究每個克隆文庫都選取了200個克隆子，基本能夠科學反映所研究樣品的優(yōu)勢細菌類群和多樣性。

[1] 林先貴,胡君利.土壤微生物多樣性的科學內(nèi)涵及其生態(tài)服務功能[J].土壤學報,2008,45(5):892-900.

[2] 蔣婧,宋明華.植物與土壤微生物在調(diào)控生態(tài)系統(tǒng)養(yǎng)分循環(huán)中的作用[J].植物生態(tài)學報,2010,34(8):979-988.

[3] O’MALLEY M A.The nineteenth century roots of everything is everywhere[J].Nature Reviews Microbiology,2007,5(8):647-651.

[4] 劉銀銀,李峰,孫慶業(yè),等.濕地生態(tài)系統(tǒng)土壤微生物研究進展[J].應用與環(huán)境生物學報,2013,19(3):547-552.

[5] PERALTA A L,MATTHEWS J W,KENT A D.Microbial community structure and denitrification in a wetland mitigation bank[J].Appl Environ Microbiol,2010,76(13):4207-4215.

[6] ANDERSEN R,GRASSET L,THORMANN M N,et al.Changes in microbial community structure and function following Sphagnum peatland restoration[J].Soil Biology and Biochemistry,2010,42:291-301.

[7] 翁永玲,官鵬.黃河三角洲鹽漬土鹽分特征研究[J].南京大學學報：自然科學版,2006,42(6):602-610.

[8] TAIT E,CARMAN M,SIEVERT S M.Phylogenetic diversity of bacteria associated with ascidians in Eel Pond (Woods Hole,Massachusetts,USA)[J].Journal of Experimental Marine Biology and Ecology,2007,342(1):138-146.

[9] SCHLOSS P D,HANDELSMAN J.Status of the microbialcensus[J].Microbiology and Molecular Biology Reviews,2004,68(4):686-691.

[10] AXELROOD P E,CHOW M L,RADOMSKI C C,et al.Molecular characterization of bacterial diversity from British Columbia forest soils subjected to disturbance[J].Canadian Journal of Microbiology,2002,48(7):655-674.

[11] 張玉琴,李文均,郝滌非,等.異常球菌屬的分類及應用研究進展[J].微生物學通報,2006,33(6):133-137.

[12] 張瑞娟,李華,林勤保,等.土壤微生物群落表征中磷脂脂肪酸(PLFA)方法研究進展[J].山西農(nóng)業(yè)科學,2011,39(9):1020-1024.

[13] 陳彥闖,辛明秀.用于分析微生物種類組成的微生物生態(tài)學研究方法[J].微生物學雜志,2009,29(4):79-83.

Analysis of Saline Vegetation Diversity of Soil Bacterial Community in the Yellow River Delta Based on 16S rDNA Gene Library

ZHANG Yue, JANG Ai-xia, GUO Du-fa*

(College of Population, Resources and Environment, Shandong Normal University, Jinan, Shandong 250014)

[Objective]The aim of the study was to investigate the soil bacterial community structure of three kinds of saline vegetation (Bare Board, Angiospermae andA.venetum) in the Yellow River Delta, and reveal its relationship with saline vegetation succession. [Method]Phylogenetic trees were constructed by the use of bacterial 16S rDNA gene library, for each of the five libraries constructed, 200 positive clones were picked randomly and sequenced, and the data were statistically analyzed. [Result]A total of 121, 150 and 155 valid sequences were obtained from bare board, Angiospermae and A. venetum. The Angiospermae soil had the highest Shannon index and ACE index. Bacteria belonging to 8 phyla were identified, like Proteobacteria, Bacteroidetes, Acidobacteria and so on. Among them, Proteobacteria accounted for 37.45%, 51.09%, 37.11% form bare board, Angiospermae andA. venetum soil, Bacteroidetes accounted for 33.02%, 3.03%, 4.47%, while the other bacterial phyla each accounted less. [Conclusion]The soil bacterial community structure of the three soil bacteria differed significantly, and Proteobacteria was the advantage microflora. The community structure of soil bacteria had significantly differences when the saline vegetation succession at the different stage.

16S rDNA; The Yellow River Delta; Bacterial community diversity; Saline vegetation

山東省自然科學基金項目(2012ZRB019DW)。

張玥(1990-)，女，山東煙臺人，碩士研究生，研究方向：環(huán)境微生物。*通訊作者。

2015-03-25

S 182

A

0517-6611(2015)13-055-03