梭梭和檉柳土壤微生物多樣性初步分析①

楊瑞紅，趙成義，王新軍，馬亞麗

(1 新疆教育學(xué)院，烏魯木齊 830043；2 中國科學(xué)院新疆生態(tài)與地理研究所，烏魯木齊 830011；3 新疆農(nóng)業(yè)大學(xué)，烏魯木齊 830052)

梭梭和檉柳土壤微生物多樣性初步分析①

楊瑞紅1,2,3，趙成義2*，王新軍3，馬亞麗1

(1 新疆教育學(xué)院，烏魯木齊 830043；2 中國科學(xué)院新疆生態(tài)與地理研究所，烏魯木齊 830011；3 新疆農(nóng)業(yè)大學(xué)，烏魯木齊 830052)

直接從沙漠土壤中提取混合微生物DNA，利用Illuminamiseq測序平臺，對16S rDNA進(jìn)行測序和分析。結(jié)果表明：①荒漠植被土壤微生物數(shù)量很少，生物活性極弱，DNA提取難度大；②沙漠植被土壤細(xì)菌多樣性豐富，所有測得細(xì)菌分屬到 23個門和 316個屬，其中還存在一定數(shù)量的微生物新種，部分測得代表新屬和種的序列提交GenBank，獲得序列號 (KT984242~KT984249)；③不同沙漠土壤樣品微生物群落有相似性，但也有較明顯的差異，其差異與土壤理化因子等因素有關(guān)：如土壤含水量與酸桿菌和變形菌的分布相關(guān)；相對于空地，梭梭和檉柳群落土壤獨有的迷蹤菌門(Elusimicrobia)與固氮菌密切相關(guān)；梭梭的土壤環(huán)境pH明顯高于檉柳土壤環(huán)境這決定了它們不同的優(yōu)勢種；土壤微生物生物量碳的大小不能反映微生物量的種類多少，但可以反映微生物數(shù)量的多少。

梭梭；檉柳；DNA提??；土壤微生物多樣性；系統(tǒng)發(fā)育分析

自20世紀(jì)50年代以來，生物多樣性與生態(tài)系統(tǒng)穩(wěn)定性的關(guān)系一直是生態(tài)學(xué)中重點討論的理論問題之一。在荒漠生態(tài)系統(tǒng)中，關(guān)注重點多集中在地上生態(tài)系統(tǒng)；而對地下生態(tài)系統(tǒng)，尤其是對微生物多樣性與系統(tǒng)穩(wěn)定性關(guān)系的研究尚重視不夠。2000年，Nature在第406期推出以“Ecology Goes Underground”為題的封面文章[1]，標(biāo)志著地下生態(tài)系統(tǒng)開始成為生物多樣性研究的關(guān)注點。事實上，植被退化導(dǎo)致根系代謝產(chǎn)物變化，直接影響土壤微生物生態(tài)環(huán)境，而土壤中微生物又維持著土壤與植被的養(yǎng)分平衡，有著重要的交互作用[2–3]，這種作用的互饋關(guān)系對整個荒漠生態(tài)系統(tǒng)功能有重要意義[4]。所以，對荒漠微生物學(xué)的研究是認(rèn)識流沙固定過程和沙地生態(tài)系統(tǒng)發(fā)育穩(wěn)定性的重要內(nèi)容之一[5]。但是荒漠植被土壤環(huán)境的極端性和大多數(shù)土壤微生物難以培養(yǎng)與過程土壤生物學(xué)機制認(rèn)識的薄弱已成為認(rèn)識固沙過程中土壤微生物和土壤養(yǎng)分循環(huán)的瓶頸，這種認(rèn)識的薄弱不僅表現(xiàn)在對特定荒漠生態(tài)系統(tǒng)研究積累的不足，而且表現(xiàn)在理論和方法上仍處于探索階段。國內(nèi)外開展關(guān)于荒漠生態(tài)系統(tǒng)土壤微生物的研究集中在：土壤微生物量的變化、種類的變化、動態(tài)變化等過程[6–8]。隨著分子技術(shù)、電子顯微鏡和土壤樣品測定分析技術(shù)的發(fā)展，人們對土壤微生物組成和多樣性的認(rèn)識也更為深入，不斷開展了荒漠土壤微生物分子實驗的研究[9–10]。Chanal等[11]分析了非洲南突尼斯沙丘的微生物多樣性，盡管環(huán)境非常惡劣但是卻發(fā)現(xiàn)較高的微生物多樣性。Subramanya等[11]用實時定量PCR對比研究了印度荒漠中沙丘、干旱、半干旱等6種土樣的微生物多樣性，共得到1 240種微生物。從維諾格拉德斯基時代的純培養(yǎng)方法到現(xiàn)在免培養(yǎng)的分子生物學(xué)方法的高通量測序等技術(shù)的進(jìn)步為更加準(zhǔn)確地揭示土壤微生物多樣性組成提供了支撐。但目前，國內(nèi)外研究結(jié)論中關(guān)于荒漠生態(tài)系統(tǒng)土壤微生物特征有其地域的限制性[12]。

為此，本研究選擇沙漠優(yōu)勢種群梭梭、檉柳土壤為研究對象，直接提取總DNA進(jìn)行高通量測序，對細(xì)菌多樣性進(jìn)行分析，較全面了解其群落組成和多樣性特點，為退化生態(tài)系統(tǒng)的恢復(fù)和重建提供土壤微生物學(xué)方面的理論數(shù)據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本研究在古爾班通古特沙漠南緣進(jìn)行調(diào)查。該沙漠面積達(dá)4.88 萬km2，為中國第二大沙漠。沙丘類型呈固定、半固定、流動狀態(tài)，以梁窩狀和樹枝狀沙壟為主，長度0.1 ~ 10 km，沙丘大多南北走向，沙壟西坡緩而長，東坡短且陡。研究區(qū)遠(yuǎn)離海洋，為溫帶干旱荒漠氣候，穩(wěn)定積雪日數(shù)100 ~ 160天，最大積雪深度多在20 cm以上，全年降水量僅在70 ~ 150 mm之間，冬春兩季降水約占總量的30% ~ 45%。全年日照總數(shù)為2 600 ~ 3 100 h，≥10℃ 積溫3 267 ~ 3 661℃，潛在年蒸發(fā)量2 000 ~ 2 800 mm。年平均氣溫5 ~ 5.7℃，極端最低氣溫小于 –40℃，極端最高溫度 40℃以上。春夏季為風(fēng)季，最大風(fēng)速20 m/s，地下水13 m左右。該沙漠沙丘表面植被覆蓋度可達(dá) 15% ~ 50%，物種豐富度相對較高，灌木種類多，藜科植物在多樣性組成中具有明顯的優(yōu)勢地位。沙漠中小半灌木、半喬木群落廣泛分布，代表植物有梭梭(Haloxylon ammnodendron)、檉柳(Tamarix ramosissima)、心葉駝絨藜(Ceretocarpus ewersmanniana)、琵琶柴(Reaumuria soongorica)、羽毛三芒草(Aristida pennata)、沙拐棗(Calligonum mongolicum)、角果藜(Ceratocarpus arenarius)等[13]。

1.2 樣品的采集及處理

2014年 10月在古爾班通古特沙漠南緣(44°22￠20.012N，87°55￠06.112E)15年以上生的梭梭、檉柳及丘間空地采集土樣，根據(jù)植被類型、土層的不同，共選擇了3個具有代表性的樣地作為對照和處理(表 1)。隨機在樣方中四角及中部分別選取土壤樣，從地表向下分層0 ~ 20 cm和20 ~ 40 cm取土樣混勻分析。取土壤樣品密封帶回實驗室進(jìn)行土壤理化性質(zhì)和分子生物實驗分析。

表1 試驗樣地基本狀況Tabble 1 Basic conditions of tested plots

1.3 主要試劑和儀器

土壤樣品的DNA用OMEGA的土壤DNA試劑盒提取。Tag DNA聚合酶和pGEM _T Easy Vector 載體購自Promega公司；蛋白酶K購自Merck公司；T4 DNA連接酶購自TaKaRa公司；PCR引物由上海Genecore公司合成；其他試劑購自Sigma等公司。

1.4 土壤樣品理化因子測定

土壤含水量測定采用105℃ 烘干比重法；土樣過0.25 mm篩，采用重鉻酸鉀氧化外加熱法和凱氏法測定土壤全碳和氮；采用蒸餾水浸提土壤測定pH；用交流測量法測定電導(dǎo)率。有效磷含量測定采用碳酸氫鈉浸提，鉬銻抗比色法測定；速效鉀含量用乙酸銨浸提火焰光度法測定。采用氯仿熏蒸法測定土壤微生物生物量碳。提取樣品設(shè)置3個重復(fù)，取其平均值。

1.5 DNA的提取

參考國內(nèi)外相關(guān)報道[14]，結(jié)合沙漠土壤樣品本身的特殊性質(zhì)，我們建立了一套有效提取沙漠土壤樣品混合基因組DNA (也叫元基因組DNA) 的方法：取新鮮的室溫過篩后的土樣10 g，用50 ml 蒸餾水充分混合搖勻溶解，1 000 g 離心10 min，棄去沉淀，再用蒸餾水溶解，12 000 g 離心10 min棄去上清，用5 ml TE緩沖液懸浮沉淀。取上述處理的沙土樣品1 ml，加入DNA抽提液5 ml(Tris_HCl 100 mmol/L，EDTA 50 mmol/L，十二烷基肌氨酸2%)，用振蕩器充分混勻；加入50 mg/ml蛋白酶K 100 μl，液氮中速凍5 min， 然后65℃水浴溶解10 min，重復(fù)3次；再使用DNA試劑盒提取，最后以75% 乙醇洗滌2 ~ 3次，干燥，溶解于適量的TE緩沖液。

1.6 PCR擴增、克隆、測序與多樣性分析

根據(jù)細(xì)菌V3+V4區(qū)設(shè)計引物，F(xiàn)：ARACTYCT ACGGRAGGCWG，R：GACTACNVGGGTATCTAA TCC，進(jìn)行PCR擴增并對其產(chǎn)物進(jìn)行純化、定量和均一化形成測序文庫。對測得序列進(jìn)一步分析，通過OTU的豐度文件，挑選出各分類學(xué)水平上豐度最高的 OTU 的序列作為該分類學(xué)水平的代表序列，使用Mothur軟件進(jìn)行多序列比對產(chǎn)生距離矩陣，采用鄰接法(Neighbor-Joining)構(gòu)建系統(tǒng)進(jìn)化樹。通過構(gòu)建物種的進(jìn)化樹，可以更深一步了解特定環(huán)境下物種的進(jìn)化關(guān)系。進(jìn)化樹中枝長的長短表示進(jìn)化距離的差異。系統(tǒng)關(guān)系越近的物種，在進(jìn)化樹中距離越近。

2 結(jié)果與分析

2.1 土樣性質(zhì)分析

土壤類型為沙地，土壤基本情況見表2。從表中可以看出，土壤0 ~ 20、20 ~ 40 cm平均含水量和土壤平均有機碳含量、有機質(zhì)含量、有效氮含量、pH及平均電導(dǎo)率有明顯的差異。平均含水量0 ~ 20 cm土層低于20 ~ 40 cm土層；有機質(zhì)含量、有效磷、有效氮等0 ~ 20 cm土層明顯高于20 ~ 40 cm土層；梭梭土樣養(yǎng)分稍高于檉柳土樣，梭梭的土壤環(huán)境明顯比檉柳土壤環(huán)境偏堿性；地表0 ~ 20 cm土樣養(yǎng)分明顯高于20 ~ 40 cm土樣的養(yǎng)分；空地的養(yǎng)分明顯少于梭梭活株和檉柳的活株；土壤微生物生物量碳梭梭土樣最高，空間地的土樣最小。

表2 不同固沙植被土壤樣品養(yǎng)分平均含量情況Table 2 Soil physical and chemical properties under different sand-bingding vegetation

2.2 土樣微生物基因組DNA提取及16S rDNA的PCR

從0.5 ~ 10 g的土樣中提取DNA，發(fā)現(xiàn)提取DNA難度較大，不僅含量偏少，而且DNA質(zhì)量也不好。當(dāng)土壤樣品經(jīng)過冷凍處理后提取 DNA會出現(xiàn)降解(圖1A)，而常溫保存提取DNA效果較好(圖1B，1C)。對 DNA提取后經(jīng)過電泳圖分析有無主帶和降解情況，測得OD260/280和OD260/230來判斷所有樣品檢測結(jié)果合格為止。以總DNA為模板和設(shè)計的引物進(jìn)行一次PCR擴增(圖1D)。利用AMPure XP磁珠對PCR產(chǎn)物進(jìn)行純化。然后定量和均一化形成測序文庫。

圖1 固沙植被土壤細(xì)菌總DNA和16S rDNA的PCR擴增瓊脂糖凝膠電泳圖Fig. 1 Amplification of bacterial 16S rDNA and bacterial total DNA from sand-binding plant soil

2.3 測序深度評估

對構(gòu)建的文庫進(jìn)行高通量測序(IlluminaMiSeq測序平臺)得到原始圖像數(shù)據(jù)文件，經(jīng)堿基識別(Base Calling)分析轉(zhuǎn)化為原始測序序列(Sequenced Reads)，結(jié)果以 FASTQ(簡稱為 fq)文件格式存儲，其中包含測序序列(reads)的序列信息以及其對應(yīng)的測序質(zhì)量信息。對原始測序序列進(jìn)行過濾，去除帶接頭(adapter)的reads和低質(zhì)量的reads得到高質(zhì)量的reads，再利用 Mothur(version 1.34.4，http://www.mothur.org/)軟件[15]，依據(jù) Needleman-wunsch算法將原始的雙端Reads拼接成一條序列即Raw Tags，進(jìn)行一系列的去冗余處理和優(yōu)化后所有樣品共得到Clean Tags條13 0429條序列。為了便于進(jìn)行物種多樣性分析將Clean Tags聚類為OTU。利用Mothur軟件以平均鄰近聚類算法(Average neighbor clustering algorithm)在 0.03(或 97% 的相似度)水平下對 Clean Tags進(jìn)行OTU的聚類，并統(tǒng)計獲得 OTU 的個數(shù)。依據(jù)每個OTU至少被測序到所有序列的 0.005% 次的原則，對獲得的OTU進(jìn)行過濾[16]，5個樣品最終共產(chǎn)生1 703個OTU(表3)。采用對測序序列進(jìn)行隨機抽樣的方法，以抽到的序列數(shù)與它們所能代表OTU的數(shù)目構(gòu)建曲線(圖2)，即稀釋性曲線[17]，它可以用來比較測序數(shù)量不同的樣本物種的豐富度。稀釋性曲線圖中，當(dāng)曲線趨向平坦時，說明測序數(shù)量合理，更多的數(shù)據(jù)量對于發(fā)現(xiàn)新的OTU貢獻(xiàn)很小，反之則表明繼續(xù)測序還可能產(chǎn)生較多新的 OTU。因此，通過作稀釋性曲線，可以對測序深度進(jìn)行評估。從圖4可以看出克隆數(shù)幾乎可以較全面地反映微生物的多樣性。覆蓋度 (Coverage C) 64.5% 從理論上講能反映環(huán)境中超過一半的微生物。

表3 樣品OTU和序列測定統(tǒng)計Table 3 The statistics of OTU and sequencing

圖2 固沙植被土壤細(xì)菌16S rRNA基因克隆文庫Rarefaction曲線Fig. 2 Rarefaction curves of 16S rRNAgene in clone library of soil bacteria under different sand-binding vegetation

表4 固沙植被土壤細(xì)菌門水平豐度統(tǒng)計Table 4 Statistics of Phylum_community (%)

2.4 細(xì)菌多樣性分析

通過Mothur軟件，將OTU中包含的51% 以上的Clean Tags都注釋為同一物種，達(dá)不到51% 的標(biāo)準(zhǔn)，則在上一個物種分類等級再進(jìn)行OTU對應(yīng)物種分析，把序列相似性大于99% 的歸結(jié)為同一個細(xì)菌的序列。結(jié)合每個OTU內(nèi)所含的Clean Tags和物種組成信息，可得到各分類學(xué)水平上每個物種的相對豐度(表4為門水平上的豐度)。發(fā)現(xiàn)測得序列可代表23個門，74個綱，131個目，218個科，316個屬細(xì)菌。其中在屬水平上，有148個為已知的屬，178個為無類別(unclassified)的屬；科水平上，121個為已知的科，97個為無類別的科；目水平上，95個為已知的目，36個為無類別的目；綱水平上，67個為已知綱，7個為無類別的綱；門水平上，22個已知門，未知的門有25種OTU。得到主要細(xì)菌類群有(表4)：綠彎菌門(Chloroflexi)、藍(lán)藻門(Cyanobacteria)、變形菌門(Proteobacteria)、厚壁菌門(Firmicutes)、放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)、芽單胞菌門(Gemmatimonadetes)、酸桿菌門(Acidobacteria)、衣原體門(Chlamydiae)、綠細(xì)菌門(Chlorobi)、迷蹤菌門(Elusimicrobia)、纖維桿菌門(Fibrobacteres)、硝化螺旋菌門(Nitrospirae)、浮霉菌門(Planctomycetes)、螺旋體門(Spirochaetes)、棲熱菌門(Thermi)等。從表4可以看出，總的來說，放線菌門為優(yōu)勢菌群占所測的菌類 36.58%；其次是變形菌門約占所測的菌類25.05%，擬桿菌門約占總數(shù)的 16.56%，其余各門所占的比例較少。從土壤層和植被上看，不同土樣各門類也有所不同，在檉柳土壤中優(yōu)勢菌群依次為：擬桿菌門、放線菌門、變形菌門，在梭梭植被土壤中優(yōu)勢菌群依次為：放線菌門、變形菌門、綠彎菌門。在空地里土壤中優(yōu)勢菌群依次為：放線菌門、變形菌門、綠彎菌門，這與梭梭土樣比較接近。

將測序中代表主要細(xì)菌類群的 OTUs 的 16S rRNA基因序列與數(shù)據(jù)庫中具有最高相似性的菌株或克隆的序列經(jīng)比對分析后構(gòu)建系統(tǒng)進(jìn)化發(fā)育樹(圖3)。由圖3可以看出不同門細(xì)菌之間的系統(tǒng)發(fā)育關(guān)系。第一優(yōu)勢菌放線菌門(Actinobacteria)分布在 74個明確分類地位的屬中，除了優(yōu)勢菌鏈霉菌(Streptomyces)外，還發(fā)現(xiàn)了擬諾卡氏菌屬(Nocardiopsis)、擬無枝酸菌屬(Amycolatopsis)、疣孢菌屬(Verrucosispore)、喜冷桿菌屬(Cryobacterium)、異壁放線菌屬(Actinoalloteichus)、假諾卡氏菌屬(Pseudonocardia)、普氏菌屬(Prauserella)、鏈單孢菌屬(Streptomonospora)、貧養(yǎng)桿菌屬(Modestobacter)、放線產(chǎn)孢菌屬(Actinomycetospora)、芽生球菌屬(Blastococcus)等稀有放線菌屬，還有29種在屬水平上未歸類無明確分類地位的屬，4個在科水平上未歸類，2個在目水平上未歸類，1個在綱水平上未歸類。擬桿菌門(Bacteroidetes)分布在23個屬中，有黃桿菌屬(Flavobacterium)、擬桿菌屬(Bacteroides)、地桿菌屬(Pedobacter)、噬纖維菌目的嗜寒嗜堿的Rhodonellum屬等已確定分類地位的屬，還有20種在屬水平上未歸類。第二優(yōu)勢菌變形菌門(Proteobacteria)里除了優(yōu)勢菌α-變形細(xì)菌生絲微菌屬(Hyphomicrobium)外，還發(fā)現(xiàn)了不動桿菌(Acinetobacter)、蛭弧菌屬(Bdellovibrio)、脫硫弧菌屬(Desulfovibrio)、亞硝化球菌屬(Nitrosococcus)、假單胞菌屬(Pseudomonas)、腸桿菌屬(Enterobacter)、厭氧螺菌屬(Anaerobiospirillum)、代爾夫特菌屬(Delftia)等共110個屬，其中52個在屬水平上未歸類的潛在新屬；有63個科，其中有20個在科水平上未歸類；有35個目類，其中有5個在目水平上未歸類；3個明確分類的綱，有一個在綱水平上未歸類。厚壁菌門(Firmicutes)中包括15個明確分類地位的屬，如梭狀芽孢桿菌綱(Clostridia)的芽孢桿菌屬(Bacillus)、類芽孢桿菌屬(Paenibacillus)、桿菌屬(Lysinibacillus)、顫螺旋菌屬(Oscillospira)等，還有另外2個在屬水平上無歸類的屬。綠彎菌門(Chloroflexi)中有厭氧繩菌綱(Anaerolineae)、熱微菌綱(Thermomicrobia)、綠彎菌綱(Chloroflexi)、爬管菌目(Herpetosiphonales)等明確分類地位的綱和目，而大部分都是在科和屬水平上未歸類的(unclassified)、無明確分類地位的共23個。綠菌門(Chlorobi)、衣原體門(Chlamydiae)、藍(lán)藻門(Cyanobacteria)等包含 OTU數(shù)量極少就不再羅列。本研究分析獲得的菌類大部分都屬于好氣性細(xì)菌，極個別是厭氧性細(xì)菌，說明好氣細(xì)菌的數(shù)量多于厭氧細(xì)菌的數(shù)量，但也有可能是厭氧細(xì)菌的DNA相對好氣細(xì)菌DNA比較難獲得導(dǎo)致的。獲得OTU中有與多種來源各環(huán)境如土壤、沙漠、火山、湖泊、海洋等的非培養(yǎng)菌類和分離獲得菌株有極高的相似性。例如有與來源于鹽堿土和鹽湖純化分離的中度嗜鹽菌株(NR_044397) 相似性為 98%；有與從沙漠公園里分離得到的噬纖維細(xì)菌科的Pontibacter屬的Pontibacter solisp.新種同源性達(dá)到99%[18]。在未歸類(unclassified)的 OTU中其相似序列大部分為非培養(yǎng)菌類。這些序列尚難確定其分類地位，可能是代表新屬和種的序列，說明固沙植被土壤存在一定數(shù)量的潛在微生物新種。這表明沙漠土壤中蘊含著豐富的未培養(yǎng)微生物資源，可以為開發(fā)和利用極端環(huán)境微生物資源提供參考。

3 討論

土壤植被可改變土壤結(jié)構(gòu)與環(huán)境，不同植物與土壤微生物的互饋效應(yīng)有所不同。所以，固沙植被土壤與空地土壤微生物有明顯的差異。從表 4 可以看出空地沒有而梭梭和檉柳群落獨有的細(xì)菌門為衣原體門(Chlamydiae)和迷蹤菌門(Elusimicrobia)。文獻(xiàn)報道迷蹤菌門(Elusimicrobia)與固氮菌密切相關(guān)且大部分屬于厭氧細(xì)菌[19]，所以在表層和空地里沒有獲得該門的DNA序列。圖4縱坐標(biāo)代表相對豐度，橫坐標(biāo)代表梭梭土樣、檉柳土樣和空地土樣，用不同的顏色代表不同科和綱?？梢钥闯鏊笏?、檉柳和空地樣品的物種和豐度有明顯的差異。梭梭的土壤環(huán)境pH明顯高于檉柳土壤環(huán)境，導(dǎo)致其優(yōu)勢菌群有所不同。有研究[20]證明放線菌具有喜熱耐堿的特性，喜歡在偏堿性環(huán)境中生長，因而在堿性較強的梭梭土壤中放線菌數(shù)量占優(yōu)勢，而在接近中性的檉柳土壤中擬桿菌占優(yōu)勢。

圖3 物種系統(tǒng)進(jìn)化樹——OTU 在屬分類學(xué)水平上的系統(tǒng)進(jìn)化樹(圖中每個葉節(jié)點對應(yīng)一個 OTU)Fig. 3 Phylogenetic tree of bacterial 16S rRNA gene clone library from soil under sand-binding vegetation

不僅不同植被的土壤樣品有差異，即使是同種植被土樣不同土層其樣品也有明顯的差異。在門和科的水平上對樣品和樣品所含菌種類進(jìn)行聚類，對聚類后的各樣品中不同OTU代表的門和科所含序列的豐度作出heatmap(圖5)。用顏色變化來反映菌群的豐度信息，可以直觀地將菌群豐度值用定義的顏色深淺表示出來，顏色越深表示相對豐度越高，顏色越淺豐度越低。從圖5可以看出在天然荒漠不同灌叢沙土和不同土壤層的土樣土壤微生物的組成整體變化具有很大的相似性，都可以分為5大枝，有一定的共性。在不同土樣的群落結(jié)構(gòu)中，變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、厚壁菌門(Firmicutes)、放線菌門(Actinobacteria)的序列總和占全部序列的大部分，這些微生物在其他相關(guān)研究的土壤中也曾被報道為優(yōu)勢菌群[21–22]。這表明盡管取樣地點與研究目的不同，但是處于相同生境中微生物類群具有相似性。

圖4 樣品科和綱水平上的物種豐度Fig. 4 Abundances of sample Family and Class levels

圖5 樣品門和科水平物種Heatmap圖Fig. 5 Heatmap of sample Family and Class levels

但由于環(huán)境因子的差異，沙土覆蓋植被和土壤深淺的變化的過程中，各種微生物受到了不同程度的影響。每種細(xì)菌的豐度有明顯的差異。圖5中如：放線菌門 (Actinobacteria)的相對含量從6.72% (BL) 上升至 57%；未分類門的含量也變化明顯不同，其中，M05樣品中的未分類門含量最多。到了種水平物種分析，通過對5種樣品的OTU進(jìn)行了統(tǒng)計制得Van圖(圖6)，可以看出5個樣品即具有相同的物種又具有自己特有的物種， 5個樣品有215種相同的物種，每個樣品又具有自己特有的物種，其中，M01、M05、M03具有60種以上自己獨有的物種。

圖6 樣品OTU Van 圖Fig. 6 Sample OTU Van map

自然界中可以影響土壤微生物群落多樣性的因素很多，如植被類型、土壤結(jié)構(gòu)、化學(xué)組成、氣候變化等[23–25]。郝金娥[26]以退化草地土壤微生物為研究對象，得出微生物生物量碳與細(xì)菌和微生物總數(shù)量相關(guān)性極顯著。本研究表明土微生物生物量碳土壤表層(0 ~ 20 cm)高于下層(20 ~ 40 cm)，這與前人[27]研究結(jié)果一致，即微生物生物量碳最大值都出現(xiàn)在0 ~ 10 cm土層，土壤微生物具有一定的垂直分布規(guī)律，因為0 ~ 10 cm 土層植被根系發(fā)達(dá)、土壤透氣性較高，這一土層為微生物提供了生長空間與絕大部分的營養(yǎng)物質(zhì)，相比較深的土層，這一土層微生物大量繁殖，活性最高。然而隨著土壤層次的加深，植物根系分布減少、土壤透氣性降低，這都不利于微生物通過吸收植物營養(yǎng)或固定大氣中的 CO2來維持其生命活動所需。而從本研究結(jié)果可以看出土壤表層微生物總量較高，但微生物種類不一定最大，這個結(jié)論有待進(jìn)一步研究。有結(jié)果表明沙漠土壤植被水分情況的變化可導(dǎo)致土壤和沉積物中含氧量的變化，進(jìn)而影響微生物的群落組成[28]。也有研究表明水分含量對于決定濕地微生物群落結(jié)構(gòu)十分重要[29]，我們的研究結(jié)果與此一致。本研究測得表層土壤水分低于深層土壤，獲得的OTU數(shù)量也是表層土壤低于深層土壤，說明在土壤表層微生物的種類略多于深層(20 ~ 40 cm)土壤。在對濕地土壤細(xì)菌多樣性的研究中[30]，得出土壤含水量的增加減少了酸桿菌的分布，而增加了變形菌的分布，我們的研究結(jié)果中也可以看出相似的趨勢。

本研究測得序列經(jīng)過與已發(fā)表的序列比對發(fā)現(xiàn)除了與已經(jīng)有明確分類地位且可以培養(yǎng)的細(xì)菌具有很近的同源關(guān)系，還與一些未培養(yǎng)細(xì)菌和沒有分類地位的菌同源關(guān)系很近。分析得到豐富的細(xì)菌多樣性與通過傳統(tǒng)培養(yǎng)細(xì)菌法在沙漠土壤中分離得到種類數(shù)量較多的微生物具有相似性[31–32]，得到的菌類數(shù)量涵蓋了 Ajar在印度的北喜馬拉雅山脈西部寒冷的沙漠土壤里收集獲得的全部代表種類(4個門31個屬82個種的232個分離菌)[33]。通過分類分析獲得的種類大部分屬好氣性細(xì)菌，群落結(jié)構(gòu)具有普遍性，即與目前已經(jīng)獲得的土壤微生物群落結(jié)構(gòu)相似，但又有其獨特的菌類，比如棲熱菌門(Thermi)。在黃河三角洲土壤微生物群落結(jié)構(gòu)對鹽生植被演替的響應(yīng)研究時，邢平平[34]將異常球菌–棲熱菌門(Thermi)的出現(xiàn)作為土壤環(huán)境惡化的標(biāo)志類群，說明此門細(xì)菌是適應(yīng)土壤缺水和營養(yǎng)而產(chǎn)生的菌類，它普遍存在于沙漠土壤和貧瘠的土壤中。

值得注意的是，由于脅迫帶來的選擇壓力可以促進(jìn)初始群落迅速進(jìn)化。本研究中樣本間的差異除了已知的菌屬外，還有很多不明確的菌屬未知種類，說明本土樣品微生物群落有很多潛在的新的未知微生物種類，其功能特性還有待于進(jìn)一步開發(fā)。由于研究的局限，很多迅速進(jìn)化的菌群還未得到認(rèn)知，所以極端環(huán)境帶來的迅速進(jìn)化可能是微生物群落的結(jié)構(gòu)和功能對特定脅迫尤其是長期脅迫產(chǎn)生抵抗力和恢復(fù)力的重要機制[35]。這也可能是沙漠植被維持沙漠生態(tài)系統(tǒng)穩(wěn)定性的重要機制之一。

4 結(jié)論與展望

通過本研究得到以下結(jié)論：①荒漠植被土壤微生物數(shù)量較一般土壤少，生物活性極弱，DNA提取難度大；②沙漠植被土壤細(xì)菌多樣性豐富，所有測得細(xì)菌分屬到23個門和316個屬，還存在一定數(shù)量的潛在微生物新種；③土壤理化因子等因素導(dǎo)致樣品分析結(jié)果有較明顯的差異，如土壤含水量與酸桿菌和變形菌的分布相關(guān)；相對于空地，梭梭和檉柳群落土壤獨有的迷蹤菌門(Elusimicrobia)與固氮菌密切相關(guān)；梭梭的土壤環(huán)境pH明顯高于檉柳土壤環(huán)境決定了它們不同的優(yōu)勢種；土壤微生物生物量碳的大小不能反映微生物生物量的種類多少，但可以反映微生物數(shù)量的多少等。

本研究通過高通量測序針對固沙植物群落土壤微生物群落進(jìn)行了研究，其成果為植被恢復(fù)和荒漠化控制提供了寶貴的參考資料，但仍有一些問題需要進(jìn)一步探討，如對于固沙植物群落深層土壤微生物的研究和微生物功能數(shù)量與土壤養(yǎng)分和植被的關(guān)系，如何使用更多手段和方法用來研究土壤微生物群落對環(huán)境脅迫的抵抗力和恢復(fù)力等，需要再進(jìn)一步研究。而且，特定的長期脅迫的環(huán)境帶來的微生物迅速進(jìn)化下樣品含有大量未知菌群，這可能為認(rèn)識沙漠土壤生態(tài)系統(tǒng)穩(wěn)定性另辟蹊徑。

[1] Copley J. Ecology goes underground[J]. Nature, 2000, 406: 452–454

[2] Abbott L K, Tang C X, Reuter D. Soil-plant-microbe interactions frommicroscopy to field practice[J]. Plant Soil, 2011, 348: 1–5

[3] Stroud J L, Paton G I, Semple K T. Microbe-aliphatic hydrocarbon interactions in soil: Implications for biodegradation and bioremediation[J]. Appl. Microbiol., 2007, 102: 1 239–1 253

[4] Cavigelly M A, Robertson P. The functional significance of denitrifier community composition in a terrestrial ecosystem[J]. Ecology, 2000:1 402–1 414

[5] 陳華癸. 土壤微生物學(xué)[M]. 上海: 上?？茖W(xué)技術(shù)出版社, 1981: 253–272

[6] 趙麗莉, 李俠, 許冬梅. 鹽池縣草地沙漠化過程中土壤微生物的變化[J]. 西北農(nóng)業(yè)學(xué)報, 2013, 22(7): 187–192

[7] Dose K, Bieger-Dose A, Ernst B, et al. Survival of microorganisms under the extreme conditions of the Atacama Desert[J]. Orig. Life Evol. Biosph., 2001, 31: 287–303

[8] Vishnevetsky S, Steinberger Y. Bacterial and fungal dynamics and their contribution to microbial biomass in desert soil[J]. J. Arid. Environ., 1997, 37(7): 83–90

[9] Ziehr L. Microbial biodiversity along an arid watershed[D]. Lubbock: Texas Tech University, 1997

[10] Chanal A, Chapon V, Benzerara K, et al. The desert of Tataouine: An extreme environment that hosts a wide diversity of microorganisms and radiotolerant bacteria[J]. Environ. Microbiol., 2006, 8: 514–525

[11] Subramanya R, Yuki C, Donnabella C. et al. Pointing microbial diversity in soil, sand dune and rock substrates of the Thar Monsoon Desert, India[J]. Indian Journal of Microbiology, 2016, 56(1): 35–45

[12] Saul-Tcherkas V, Steingerber Y. Temporal and shrub adaptaion affact on soil microbial functional diversity in a desert system[J]. European Journal of Soil Science, 2009, 60: 871–882

[13] 錢亦兵, 張立運, 吳兆寧, 等. 新疆準(zhǔn)噶爾盆地邊緣部分地段生態(tài)環(huán)境特征[J]. 干旱區(qū)地理, 2003, 26(1): 30–36

[14] Ana-Maria T, Ioana M, Iulia C, et al. Comparison of total DNA extraction methods for microbial community form polluted soil original research[J]. Agriculture and Agricultural Science Procedia, 2015, 6: 616–622

[15] Schloss P D, Westcott S L, Ryabin T. Introducing mothur: Open-source, platform-independent, community-cupported software for describing and comparing microbial communities[J]. Appl. Environ. Microbiol., 2009, 5(23): 7 537–7 541

[16] Bokulich N A, Subramanian S, Faith J J, et al. Quality-filtering vastly improves diversity estimates from Illuminaamplicon sequencing[J]. Nature methods, 2013, 10(1): 57–59

[17] Amato K R, Yeoman C J, Kent A, et al. Habitat degradation impacts black howler monkey (Alouattapigra) gastrointestinal microbiomes[J]. The ISME Journal, 2013, 7(7): 1 344–1 353

[18] Dai J, Xu M B, Peng F, et al. Pontibacter soli sp. nov., isolated from the soil of a Populus rhizosphere in Xinjiang, China[J]. Antonie Van Leeuwenhoek, 2014, 105(1): 65–72

[19] Zheng H, Dietrich C, Radek R. Endomicrobium proavitum, the first isolate of Endomicrobia class. nov. (phylum Elusimicrobia)–an ultramicrobacterium with an unusual cell cycle that fixes nitrogen with a Group IV nitrogenase[J]. Environmental Microbiology, 2016, 18(1): 191–204

[20] 李鑫, 張會慧, 岳冰冰, 等. 桑樹–大豆間作對鹽堿土碳代謝微生物多樣性的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2012, 23(7): 1 825–1 831

[21] Hollister E B, Engledow A S, Hammett A J, et al. Shifts in microbial community structure along an ecological gradient of hypersaline soils and sediments[J]. ISME J., 2010, 4(6): 829–838

[22] Will C, ThurmerA, Wollherr A, et al. Horizon-specific bacterial community composition of German grassland soils, as revealed by pyrosequencing-based analysis of 16S rRNA genes[J]. Applied Environmental Microbiology, 2010, 76: 6 751–6 758

[23] Truu M J, Juhanson J T. Microbial biomass, activity and community composition in constructed wetlands[J]. Science of Total Environment, 2009, 407: 3 958–3 971

[24] ZhaoY J, Liu B, Zhang W G. An Effects of plant and influent C︰N︰P ratio on microbial diversity in pilot-scale constructed wetlands[J]. Ecological Engineering, 2010, 36: 441–449

[25] 鐘芳, 柴曉虹, 王國基, 等. 植被恢復(fù)方式對黃土丘陵區(qū)土壤理化性質(zhì)及微生物特性的影響[J]. 中國沙漠, 2014, 34(4): 1 064–1 072

[26] 郝金娥. 江河源區(qū)不同植被演替階段土壤微生物數(shù)量及微生物生物量碳的特性[D]. 西寧: 青海大學(xué), 2010

[27] 盧虎, 姚拓, 李建宏, 等. 高寒地區(qū)不同退化草地植被和土壤微生物特性及其相關(guān)性研究[J]. 草業(yè)學(xué)報, 2015, 24(5): 34–43

[28] 楊航宇, 劉艷梅, 王廷璞. 荒漠區(qū)生物土壤結(jié)皮對土壤酶活性的影響[J]. 土壤學(xué)報, 2015, 52(3): 654–664

[29] TangY S, Wang L, Jia J W, et al. Response of soil microbial community in Jiuduansha wetland to different successional stages and its implications for soil microbial respiration and carbon turnover[J]. Soil Biology and Biochemistry, 2011, 43: 638–646

[30] 隋心, 張榮濤, 鐘海秀, 等. 利用高通量測序?qū)θ皆∪~章濕地土壤細(xì)菌多樣性的研究[J]. 土壤, 2015, 47(5): 919–925

[31] 羅明, 單娜娜, 文啟凱, 等. 幾種固沙植物根際土壤微生物特性研究[J]. 應(yīng)用與環(huán)境生物學(xué)報, 2002, 8(6): 618–622

[32] Hanna A L, Youssef H H, Amer W M, et al. Diversity of bacteria nesting the plant cover of north Sinai deserts, Egypt[J]. Journal of Advanced Research, 2013, 4:13–26

[33] Yadav A N, Sachan S G, Verma P, et al. Prospecting cold deserts of north western Himalayas for microbial diversity and plant growth promoting attributes[J]. Journal of Bioscience and Bioengineering, 2015, 119(6): 683–693

[34] 邢平平. 黃河三角洲土壤微生物群落結(jié)構(gòu)對鹽生植被演替的響應(yīng)[D]. 濟南: 山東師范大學(xué), 2013

[35] 關(guān)統(tǒng)偉. 新疆羅布泊鹽湖放線菌多樣性及多相分類[D].雅安: 四川農(nóng)業(yè)大學(xué), 2012

Phylogenetic Diversity Preliminary Analysis ofHaloxylon ammnodendronandTamarix ramosissimaSoil Bacteria

YANG Ruihong1,2,3, ZHAO Chengyi2*, WANG Xinjun3, MA Yali1
(1Xinjiang Education Institute,Urumqi840043,China; 2Xingjiang Institute of Ecology and Geography,Chinese Academy of Sciences,Urumqi830011,China; 3Xijiang Agriculture University,Urumqi830052,China)

Microbiological mixed DNA (genomic DNA) was directly extracted from soil in the Gurbantunggut desert. Using Illuminamiseq sequencing platform, diversity of community structure was analyzed from the bacterial 16S rDNA. The research showed that: 1) Because of the low quantity of microbe in desert vegetation soil and the low activity, there were many difficulties in the DNA extraction; 2) In total, measured sequences belonged to 23 phyla and 316 genera, among them, some sequences were uncultured groups that had an uncertain affiliation, which had been submitted to the GenBank PEMBL PDDBJ databases under accession numbers KT984242–KT984249. 3)There were some discrepancies between groups. The difference was affected by factors such as soil physical and chemical factors. The water content influenced the amount ofProteobacteriaandAcidobacteria. The uniqueElusimicrobiawas the soil of plantHaloxylon ammnodendronandTamarix ramosissima,which was closely related to the nitrogen fixing bacteria. The soil microbial biomass carbon could reflect the quantity of microorganism but not the types of microorganism.

Haloxylon ammnodendron;Tamarix ramosissima; DNA extraction; Microbiology diversity; Phylogenetic analysis

Q938

10.13758/j.cnki.tr.2016.06.009

國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)項目(2013CB429905)和國家自然科學(xué)基金項目(41301205、2014211B014)資助。

* 通訊作者(zcy@ms.xjb.ac.cn)

楊瑞紅(1980—), 女, 新疆人, 副教授, 博士研究生, 主要研究方向為荒漠土壤生態(tài)與微生物學(xué)。E-mail: yangruihong99@163.com