• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    呼倫貝爾沙區(qū)土壤細(xì)菌群落結(jié)構(gòu)與功能預(yù)測(cè)

    2019-11-28 10:13:24杜宇佳高廣磊陳麗華丁國(guó)棟曹紅雨
    中國(guó)環(huán)境科學(xué) 2019年11期
    關(guān)鍵詞:樟子松天然林沙地

    杜宇佳,高廣磊*,陳麗華,丁國(guó)棟,張 英,曹紅雨

    呼倫貝爾沙區(qū)土壤細(xì)菌群落結(jié)構(gòu)與功能預(yù)測(cè)

    杜宇佳1,2,高廣磊1,2*,陳麗華1,丁國(guó)棟1,2,張 英1,曹紅雨1,2

    (1.北京林業(yè)大學(xué)水土保持學(xué)院,水土保持國(guó)家林業(yè)與草原局重點(diǎn)實(shí)驗(yàn)室,北京 100083;2.寧夏鹽池毛烏素沙地生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站,寧夏 鹽池 751500)

    以呼倫貝爾沙區(qū)裸沙地、草地、沙地樟子松(var.)人工林和沙地樟子松天然林四種生境土壤為研究對(duì)象,采用野外調(diào)查、16S rRNA基因高通量測(cè)序和PICRUSt功能預(yù)測(cè)相結(jié)合的研究方法比較分析不同生境土壤細(xì)菌群落結(jié)構(gòu)和潛在功能組成特征.結(jié)果顯示:呼倫貝爾沙區(qū)沙地樟子松天然林土壤細(xì)菌多樣性最高,人工林土壤細(xì)菌多樣性最低,Shannon指數(shù)分別為(8.623±0.193)和(7.432±0.028),不同生境土壤細(xì)菌alpha和beta多樣性存在顯著差異.草地、沙地樟子松人工林和天然林土壤中變形菌門(Proteobacteria)相對(duì)豐度最高,均值分別為29.83%±1.14%、34.73%±1.99%、31.95%±0.21%,裸沙地土壤放線菌門(Actinobacteria)相對(duì)豐度最高,均值為26.13%±0.43%.不同生境土壤細(xì)菌主要優(yōu)勢(shì)屬為慢生根瘤菌屬()、RB41,其相對(duì)豐度在四種生境中的均值分別為5.29%±2.24%、4.22%±1.23%.PICRUSt功能預(yù)測(cè)共得到6個(gè)一級(jí)功能層,40個(gè)二級(jí)功能層,土壤細(xì)菌功能較為豐富,土壤細(xì)菌群落在環(huán)境信息處理、代謝、遺傳信息處理和有機(jī)系統(tǒng)方面功能活躍.沙地樟子松天然林核苷酸代謝、酶家族、氨基酸代謝、碳水化合物代謝功能基因較為豐富,保證了沙地樟子松天然林土壤細(xì)菌的存活,使其具有較高的多樣性.呼倫貝爾沙區(qū)不同生境土壤細(xì)菌功能基因豐度波動(dòng),反映了四種生境的土壤細(xì)菌群落組成及多樣性的變化,指示了不同生境功能基因?qū)ν寥兰?xì)菌群落的影響規(guī)律,可為預(yù)測(cè)和理解沙區(qū)土壤細(xì)菌代謝潛力和功能提供參考借鑒.

    呼倫貝爾沙區(qū);生境;細(xì)菌群落結(jié)構(gòu);16S rRNA;功能預(yù)測(cè)

    土壤細(xì)菌是土壤微生物中種類最多、數(shù)量最大、功能最豐富的類群[1],是陸地生態(tài)系統(tǒng)的重要組成部分.土壤細(xì)菌驅(qū)動(dòng)地球生物化學(xué)循環(huán),是生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)的關(guān)鍵環(huán)節(jié),對(duì)于穩(wěn)定、調(diào)節(jié)和修復(fù)陸地生態(tài)系統(tǒng)發(fā)揮著重要作用[2].土壤細(xì)菌的空間分布格局存在顯著的生境依賴性特征[3],受土壤性質(zhì)、植被、氣候等多種環(huán)境因素的綜合影響.在早期研究中,相關(guān)學(xué)者主要關(guān)注土壤細(xì)菌的群落結(jié)構(gòu)與多樣性[4-5],而后聚焦于土壤細(xì)菌群落結(jié)構(gòu)與分布的影響因素與驅(qū)動(dòng)過(guò)程[6].近年來(lái),隨著分子生物學(xué)技術(shù)與方法的不斷發(fā)展,土壤細(xì)菌相關(guān)研究開(kāi)始由結(jié)構(gòu)向功能轉(zhuǎn)變[7-8].各國(guó)學(xué)者共同致力于土壤細(xì)菌功能預(yù)測(cè)與分析,并試圖通過(guò)土壤細(xì)菌功能研究,揭示土壤細(xì)菌在陸地生態(tài)系統(tǒng)過(guò)程中發(fā)揮的重要作用[9-11].

    目前,生物信息技術(shù)發(fā)展迅速.PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), Tax4Fun和FAPROTAX (Functional Annotation of Prokaryotic Taxa)等方法基于系統(tǒng)發(fā)育與功能充分聯(lián)系的假設(shè),根據(jù)16S rRNA測(cè)序基因與參考基因組數(shù)據(jù)庫(kù)的相似性,進(jìn)而判斷和預(yù)測(cè)標(biāo)記基因的功能組成,為我們深入理解微生物功能提供了新途徑[12-13].其中,FAPROTAX依據(jù)物種名進(jìn)行功能預(yù)測(cè),但僅涉及元素循環(huán)相關(guān)功能,無(wú)法對(duì)整體功能進(jìn)行預(yù)測(cè),算法存在明顯的局限性;Tax4Fun基于SILVA數(shù)據(jù)庫(kù)進(jìn)行基因功能預(yù)測(cè),雖然該數(shù)據(jù)庫(kù)擁有細(xì)菌、真菌、古菌rRNA基因序列,但數(shù)據(jù)庫(kù)剛剛建立,尚不完善,應(yīng)用仍處于起步階段;PICRUSt的比對(duì)數(shù)據(jù)庫(kù)為Greengenes,該數(shù)據(jù)庫(kù)涵蓋細(xì)菌、古菌基因序列,數(shù)據(jù)庫(kù)相對(duì)完善,具有預(yù)測(cè)準(zhǔn)確,功能先進(jìn),方便快捷,成本低廉的鮮明特點(diǎn).目前,PICRUSt算法已在土壤、海洋和湖泊細(xì)菌功能分析得到應(yīng)用[14-15],為土壤微生物群落結(jié)構(gòu)與功能研究提供了新途徑[16-17].

    呼倫貝爾沙區(qū)是我國(guó)防沙治沙的重點(diǎn)地區(qū)之一,也是構(gòu)建“兩屏三帶”國(guó)家生態(tài)格局“北方防沙帶”的核心地區(qū).相關(guān)學(xué)者已對(duì)呼倫貝爾沙區(qū)土壤微生物數(shù)量、多樣性、群落結(jié)構(gòu)進(jìn)行了大量的研究探索,取得了豐碩的研究成果[18-19].隨著研究的不斷深入,相關(guān)研究重點(diǎn)已由土壤微生物群落結(jié)構(gòu)向功能轉(zhuǎn)移.鑒于此,本研究以呼倫貝爾沙區(qū)裸沙地、草地、樟子松天然林、樟子松人工林四種典型生境為研究對(duì)象,采用PICRUSt方法分析土壤細(xì)菌群落結(jié)構(gòu)與功能特征,以期研究揭示不同生境土壤微生物群落結(jié)構(gòu)與功能的相互關(guān)系,為預(yù)測(cè)和理解沙區(qū)土壤細(xì)菌代謝潛力和功能提供參考借鑒.

    1 材料與方法

    1.1 研究區(qū)概況

    呼倫貝爾沙區(qū)位于內(nèi)蒙古自治區(qū)東北部(115°31′~126°04′E,47°05′~53°20′N),屬半濕潤(rùn)大陸性季風(fēng)氣候區(qū).年均氣溫-0.3℃,年均降水量280~ 355mm,年均蒸發(fā)量1400~1900mm,無(wú)霜期121d.地帶性土壤主要由黑鈣土和栗鈣土組成,風(fēng)沙土則多分布于沙帶沙丘及沙質(zhì)草原沙區(qū),境內(nèi)包括海拉爾河流域和伊敏河流域.研究區(qū)喬木植被主要有沙地樟子松(var)、白樺()等,灌木植被主要有小葉錦雞兒()、黃柳()和沙蒿()等,草本植被主要有羊草()、針茅()和沙蓬()等.

    1.2 土壤樣品采集

    2017年7月,在裸沙地、草地、沙地樟子松人工林和沙地樟子松天然林4種典型生境地勢(shì)平坦地段分別布設(shè)3個(gè)20×20m的樣地(表1),開(kāi)展土壤樣品采集工作,樣品采集深度為0~10和10~20cm.其中,沙地樟子松天然林和人工林樣地內(nèi)分別選取3株長(zhǎng)勢(shì)相當(dāng)?shù)臉?biāo)準(zhǔn)木,在樹(shù)冠投影處采集土壤樣品,草地和裸沙地則在樣地內(nèi)隨機(jī)選取3個(gè)采樣點(diǎn)采集土壤樣品.新鮮土壤樣品去除石子、根系和碎草等雜質(zhì)后裝入無(wú)菌袋,置于4 ℃便攜式保溫箱中冷凍保存,運(yùn)回實(shí)驗(yàn)室后采用-80 ℃的冰箱保存.

    表1 樣地基本概況

    注:HLB:bare sandy land,裸沙地;HLG:grass land,草地;HLN:natural forest,樟子松天然林;HLP:plantation,樟子松人工林.

    1.3 土壤樣品處理及分析

    1.3.1 土壤樣品處理及理化性質(zhì)分析 將土壤樣品風(fēng)干后,過(guò)2mm篩后進(jìn)行樣品理化性質(zhì)的分析.土壤含水量采用烘干測(cè)定法;土壤pH值采用電位法使用便攜式pH計(jì)測(cè)定;土壤有機(jī)質(zhì)采用稀釋熱法測(cè)定;土壤全氮含量、速效氮含量采用靛酚藍(lán)比色法測(cè)定;土壤全磷含量、速效磷含量采用鉬銻抗比色法測(cè)定[20].

    1.3.2 土壤細(xì)菌的分離鑒定 取5g采集的土壤樣品于加入液氮的研缽中研磨,使用Power Soil DNA Kit(MoBio,USA)試劑盒提取土壤微生物基因組DNA,置-20 ℃的冰箱中保存.對(duì)16S rRNA基因的V3-V4區(qū)域進(jìn)行擴(kuò)增,引物序列為338F(ACTCCT- ACGGGAGGCAGCAG)、806R(GGACTACHVGG- GTWTCTAAT).PCR采用25mL反應(yīng)體系:1mL× Forward Primer(5mmol/L),1mL×Reverse Primer (5mmol/L),3mL BSA(2ng/mL),12.5mL 2×Taq PCR MasterMix,30ng DNA模板,最后用7.5mL ddH2O補(bǔ)足至總體積25mL.PCR擴(kuò)增程序?yàn)?95℃預(yù)變性300s; 95℃變性45s,55℃退火50s,72℃延伸45s,持續(xù)28個(gè)循環(huán)周期;最后72℃延伸10min.將樣本PCR產(chǎn)物混合后用2%瓊脂糖凝膠電泳檢測(cè),使用AxyPrepDNA凝膠回收試劑盒(AXYGEN公司)切膠回收PCR產(chǎn)物,Tris-HCl洗脫;2%瓊脂糖電泳檢測(cè).合格的PCR之后進(jìn)行MiSeq文庫(kù)構(gòu)建并測(cè)序.本研究涉及的所有測(cè)序數(shù)據(jù)可在NCBI Sequence Read Archive下載(PRJNA551926).

    1.3.3 測(cè)序數(shù)據(jù)處理 Miseq測(cè)序得到的序列利用Trimmomatic對(duì)雙端序列數(shù)據(jù)進(jìn)行過(guò)濾處理,過(guò)濾read尾部質(zhì)量值20以下的堿基,設(shè)置50bp的窗口,如果窗口內(nèi)的平均質(zhì)量值低于20,從窗口開(kāi)始截去后端堿基,過(guò)濾質(zhì)控后50bp以下的read;利用Pear對(duì)含有barcode的數(shù)據(jù)進(jìn)行質(zhì)控過(guò)濾(根據(jù)-value和-value校正,-value值0.0001);然后根據(jù)PE測(cè)序的overlap關(guān)系,利用FLASH將成對(duì)的序列拼接(merge)成一條序列,拼接序列的overlap區(qū)允許的最大錯(cuò)配比率為0.1.根據(jù)序列首尾兩端的barcode和引物區(qū)分樣品,并調(diào)整序列方向,barcode允許的錯(cuò)配數(shù)為0;利用usearch軟件去除嵌合體.下機(jī)數(shù)據(jù)在去除barcode和primer并拼接后得到raw tags,raw tags經(jīng)進(jìn)一步去除嵌合體、短序列后得到優(yōu)質(zhì)序列clean tags,共得到優(yōu)質(zhì)序列1472261條.其中,裸沙地、草地、沙地樟子松天然林和人工林優(yōu)質(zhì)序列分別有303383、353675、422576、392627條.為最小化樣本變化對(duì)測(cè)序效率的影響[21],根據(jù)優(yōu)質(zhì)序列進(jìn)一步去除singlton序列,得樣本最低序列數(shù)99173條,按此樣本進(jìn)行隨機(jī)抽平.利用QIIME v.1.8軟件對(duì)優(yōu)質(zhì)序列進(jìn)行質(zhì)控,通常對(duì)97%的相似水平的物種分類單元(OTU)進(jìn)行信息統(tǒng)計(jì)分析[22].

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

    將QIIME軟件得到的biom文件上傳到Galaxy網(wǎng)站進(jìn)行PICRUSt功能基因預(yù)測(cè)分析,得到土壤細(xì)菌功能基因組成.采用mothur軟件計(jì)算細(xì)菌群落Chao1豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)、Simpson優(yōu)勢(shì)度指數(shù).采用SPSS 24.0進(jìn)行單因素方差分析(ANOVA)、最小顯著差異法(LSD)比較、Pearson相關(guān)性分析.weighted Nearest Sequenced Taxon Index(weighted NSTI)是將樣本中的每個(gè)OTU與已知數(shù)據(jù)庫(kù)中參考細(xì)菌基因組分開(kāi)的平均分支長(zhǎng)度,通過(guò)樣本中OTU的豐度加權(quán)得到,使用predict_metagenomes.py命令-a選項(xiàng)計(jì)算NSTI.采用Excel 2010繪制相對(duì)豐度圖.采用Origin 2018繪制土壤細(xì)菌群落組成與功能熱圖.采用Canoco for Windows進(jìn)行土壤細(xì)菌PCA與RDA分析.

    2 結(jié)果與分析

    2.1 不同生境土壤細(xì)菌多樣性特征

    圖1 不同生境土壤細(xì)菌群落多樣性指數(shù)

    不同大小寫字母分別表示0~10cm、10~20cm土層不同生境土壤細(xì)菌多樣性差異顯著(<0.05).橫坐標(biāo)中a為0~10cm土層深度,b為10~20cm 土層深度,下同

    呼倫貝爾沙區(qū)不同生境土壤細(xì)菌的Chao1豐富度指數(shù)、Shannon多樣性指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)組內(nèi)和組間差異顯著(<0.05)(圖1).草地和沙地樟子松天然林的土壤細(xì)菌多樣性指數(shù)較高,分別為8.60±0.06和8.62±0.19,其次是裸沙地為8.52±0.26,沙地樟子松人工林土壤細(xì)菌多樣性指數(shù)最低為7.43±0.03.草地和沙地樟子松天然林土壤細(xì)菌群落豐富度指數(shù)無(wú)顯著性差異(0.05),且明顯高于裸沙地與沙地樟子松人工林.裸沙地、草地、沙地樟子松天然林土壤細(xì)菌多樣性指數(shù)和優(yōu)勢(shì)度指數(shù)均無(wú)顯著性差異(0.05),均高于沙地樟子松人工林.由四種生境土壤細(xì)菌分布的垂直空間分析,裸沙地、草地、沙地樟子松人工林,土壤細(xì)菌群落豐富度指數(shù)和多樣性指數(shù)均隨著土層深度的增加而增加,而沙地樟子松天然林變化規(guī)律與之相反.

    2.2 不同生境土壤細(xì)菌群落結(jié)構(gòu)特征

    不同生境土壤細(xì)菌類群主要優(yōu)勢(shì)細(xì)菌門(圖2)是變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria),其在四種生境中的均值分別為29.76%±4.43%、22.11%±4.46%、18.71%±4.62%,主要優(yōu)勢(shì)細(xì)菌屬是慢生根瘤菌屬()、RB41,其相對(duì)豐度的均值分別為5.29%±2.24%、4.22%±1.23%.草地、沙地樟子松人工林、沙地樟子松天然林中變形菌門的相對(duì)豐度最高,其均值分別為29.83%±1.14%、34.73%±1.99%和31.95%±0.21%,裸沙地中放線菌門的相對(duì)豐度最高,為26.13%±0.43%.酸桿菌門在沙地樟子松天然林中相對(duì)豐度最高,為25.80%±4.33%.變形菌門的相對(duì)豐度在不同土層分布無(wú)差異,放線菌門和酸桿菌門在不同土層深度其相對(duì)豐度表現(xiàn)出波動(dòng).

    PC1軸的可信度是27.50%,PC2軸的可信度為17.40%(圖3).同一生境樣點(diǎn)相對(duì)聚集,不同生境樣點(diǎn)相對(duì)分開(kāi),表明四種生境細(xì)菌群落組成存在明顯差異.草地、沙地樟子松人工林和天然林3種生境內(nèi)樣點(diǎn)相對(duì)聚集,表明這3種生境不同土層細(xì)菌群落組成較為接近.裸沙地內(nèi)樣點(diǎn)距離其他3種生境最遠(yuǎn),由此推斷,裸沙地中土壤細(xì)菌群落結(jié)構(gòu)與其他3種生境相比存在較大差異.裸沙地0~10和10~20cm兩個(gè)土層中,各有一個(gè)點(diǎn)與樣地內(nèi)其他樣點(diǎn)距離較遠(yuǎn),表明裸沙地0~10和10~20cm土層細(xì)菌群落組成存在變異性

    圖2 不同生境土壤細(xì)菌群落組成(門水平)

    圖中僅列出相對(duì)豐度大于1%的細(xì)菌門,小于此值則計(jì)為other

    圖3 不同生境土壤細(xì)菌群落組成主成分分析

    2.3 不同生境土壤理化性質(zhì)與細(xì)菌群落的RDA分析及與細(xì)菌多樣性指數(shù)的相關(guān)性

    土壤理化性質(zhì)見(jiàn)表2.由圖4可見(jiàn),第一排序軸解釋了不同生境土壤理化性質(zhì)與細(xì)菌群落關(guān)系的63.7%,第二排序軸解釋了18.6%的變異,2個(gè)排序軸累計(jì)貢獻(xiàn)率為82.3%.不同生境樣本點(diǎn)分布相對(duì)離散且距離原點(diǎn)較遠(yuǎn),說(shuō)明不同生境細(xì)菌群落結(jié)構(gòu)特征存在差異且受土壤理化性質(zhì)影響較大.裸沙地土壤細(xì)菌群落結(jié)構(gòu)與土壤含水量有相關(guān)性(<0.05);沙地樟子松人工林土壤細(xì)菌群落結(jié)構(gòu)與土壤有機(jī)質(zhì)含量、土壤速效磷含量有顯著的相關(guān)性(<0.05);沙地樟子松天然林土壤細(xì)菌群落結(jié)構(gòu)主要受到土壤pH值的影響.土壤理化性質(zhì)與土壤細(xì)菌群落多樣性指數(shù)有相關(guān)關(guān)系(表3),其中,土壤細(xì)菌群落豐富度指數(shù)與土壤速效氮含量達(dá)到顯著相關(guān)(<0.05),與土壤速效磷含量達(dá)到極顯著相關(guān)(<0.01),土壤細(xì)菌群落Shannon指數(shù)與土壤速效磷含量達(dá)到顯著負(fù)相關(guān)(<0.05),表現(xiàn)為沙地樟子松人工林土壤細(xì)菌多樣性在4種生境中最低且土壤速效磷含量最高.

    表2 不同生境土壤理化性質(zhì)

    圖4 不同生境土壤細(xì)菌群落結(jié)構(gòu)與環(huán)境因子的冗余分析(RDA)

    SWC:土壤含水量,pH:土壤pH值,SOM:土壤有機(jī)質(zhì)含量,TN:土壤全氮含量,TP:土壤全磷含量,AN:土壤速效氮含量,AP:土壤速效磷含量

    表3 土壤理化性質(zhì)與土壤細(xì)菌多樣性指數(shù)之間的相關(guān)關(guān)系

    注:**:相關(guān)性顯著(<0.01),*:相關(guān)性顯著(<0.05).

    2.4 不同生境土壤細(xì)菌PICRUSt基因預(yù)測(cè)

    2.4.1 功能基因家族組成 Weighted NSTI是用來(lái)評(píng)估給定樣品中微生物與測(cè)序基因組相關(guān)的程度.裸沙地、草地、沙地樟子松天然林和人工林NSTI指數(shù)分別為0.202~0.223、0.184~0.187、0.172~0.182、0.171~0.172(表4).NSTI指數(shù)的數(shù)值越小,說(shuō)明樣品功能預(yù)測(cè)與已知數(shù)據(jù)庫(kù)的匹配度越高.

    表4 不同生境土壤細(xì)菌NSTI指數(shù)

    呼倫貝爾沙區(qū)不同生境中獲得的一級(jí)生物代謝功能通路分析相對(duì)豐度大于1%的功能基因包括4種(圖5):環(huán)境信息處理(Environmental Information Processing)、代謝(Metabolism)、遺傳信息處理(Genetic Information Processing)、有機(jī)系統(tǒng)(Organismal Systems).每種功能基因在不同生境中相對(duì)豐度基本一致,除沙地樟子松天然林功能基因相對(duì)豐度隨土層深度的變化而改變,其他3種生境在不同土層功能基因的相對(duì)豐度沒(méi)有改變.且呼倫貝爾沙區(qū)不同生境中樟子松天然林0~10cm土層功能基因與其他3種生境0~10cm土層功能基因存在顯著差異(<0.05).沙地樟子松天然林中代謝、遺傳信息處理、有機(jī)系統(tǒng)功能基因的相對(duì)豐度顯著高于其他生境,環(huán)境信息處理功能基因的相對(duì)豐度顯著低于其他生境.

    圖5 不同生境土壤細(xì)菌功能預(yù)測(cè)(一級(jí)功能層)

    以代謝通路相對(duì)豐度大于1%作圖

    2.4.2 功能基因家族差異分析 6個(gè)一級(jí)功能層預(yù)測(cè)得到40個(gè)二級(jí)功能層中,其中相對(duì)豐度在0.1%以上的二級(jí)功能層有13個(gè)(圖6).4種生境的功能基因在二級(jí)功能層存在差異.其中,0~10cm土層中沙地樟子松天然林膜運(yùn)輸(Membrane Transport)的功能基因豐度顯著低于其他生境,核苷酸代謝(Nucleotide Metabolism)、酶家族(Enzyme Families)、氨基酸代謝(Amino Acid Metabolism)、碳水化合物代謝(Carbohydrate Metabolism)的功能基因豐度顯著高于其他生境.10~20cm土層中裸沙地中核苷酸代謝與酶家族的功能基因豐度顯著高于其他生境,草地和沙地樟子松天然林在該土層中氨基酸代謝的功能基因豐度偏低.

    PC1軸的可信度是72.1%,PC2軸的可信度是19.7%(圖7).草地與沙地樟子松人工林兩種生境樣點(diǎn)相對(duì)聚集,表明兩種生境內(nèi)土壤細(xì)菌功能組成情況較為相似.裸沙地與沙地樟子松天然林內(nèi)樣點(diǎn)與另兩種生境距離較遠(yuǎn),由此推斷,這兩種生境與草地和沙地樟子松人工林土壤細(xì)菌功能組成存在較大差異.沙地樟子松天然林10~20cm土層中,有一個(gè)點(diǎn)與樣地內(nèi)其他樣點(diǎn)距離較遠(yuǎn),表明該生境10~20cm土層細(xì)菌功能組成情況存在變異性.同一生境距離之間的歐幾里得距離較同一土層之間的歐幾里得距離較近,因此,不同生境對(duì)細(xì)菌功能組成情況影響大于不同土層深度對(duì)細(xì)菌功能組成情況的影響.

    圖6 不同生境土壤細(xì)菌功能預(yù)測(cè)熱圖(二級(jí)功能層)

    僅展示代謝通路相對(duì)豐度大于0.1%功能

    圖7 不同生境土壤細(xì)菌功能組成主成分分析

    3 討論

    3.1 生境對(duì)土壤細(xì)菌多樣性影響

    呼倫貝爾沙區(qū)不同生境土壤細(xì)菌多樣性存在差異,這與前人研究結(jié)果一致.主要是因?yàn)椴煌惩寥乐兴牡蚵湮飻?shù)量、有機(jī)質(zhì)含量、根系分泌物均存在差異,這使得每種生境有其特定的微環(huán)境,從而影響細(xì)菌多樣性[23-24],土壤細(xì)菌多樣性與土壤理化性質(zhì)有關(guān),沙地樟子松人工林土壤含水量最低,pH值酸性最強(qiáng),全磷、速效磷含量顯著降低,從而使得該生境土壤細(xì)菌多樣性最低[25].不同生境土壤細(xì)菌群落多樣性隨土層深度的變化而改變,這與前人研究結(jié)果一致,裸沙地、草地、沙地樟子松人工林3種生境中10~20cm土層土壤含水量高于0~10cm土層土壤含水量,該生境土壤細(xì)菌多樣性隨土層深度增加而增加.但在沙地樟子松天然林中卻出現(xiàn)了相反的趨勢(shì),由于天然林表層土壤有良好的營(yíng)養(yǎng)和通氣條件,土壤結(jié)構(gòu)疏松,細(xì)菌群落結(jié)構(gòu)與土壤孔隙顯著相關(guān),天然林土壤的孔隙度不同于其他3種生境,細(xì)菌在表層土壤生長(zhǎng)迅速,從而提高了0~10cm土層的細(xì)菌多樣性[26-27].磷元素是植物生長(zhǎng)過(guò)程中所需要的最重要元素之一,植物對(duì)土壤中相對(duì)較高含量的磷元素有相對(duì)較低的吸收率[28],土壤速效磷含量與變形菌門、酸桿菌門存在相關(guān)性[29-30],土壤中存在的這些微生物能夠幫助植物吸收磷元素[31].

    3.2 生境對(duì)土壤細(xì)菌群落結(jié)構(gòu)的影響

    呼倫貝爾沙區(qū)土壤細(xì)菌優(yōu)勢(shì)門與陜北沙化區(qū)、黃土高原區(qū)、內(nèi)蒙荒漠草原的細(xì)菌優(yōu)勢(shì)門類別一致[32-34].不同生境中土壤優(yōu)勢(shì)細(xì)菌門的相對(duì)豐度不同是由各生境環(huán)境狀況的不同引起的,不同生境會(huì)影響土壤氮、磷含量,導(dǎo)致細(xì)菌群落結(jié)構(gòu)發(fā)生變化[35].由于變形菌門是營(yíng)養(yǎng)性的,對(duì)碳有很高的利用率[36],草地、沙地樟子松人工林和天然林中相對(duì)較高的土壤有機(jī)質(zhì)含量保證了沙地地區(qū)細(xì)菌的生長(zhǎng).由于放線菌門除對(duì)干旱有很強(qiáng)的耐受性外,還可以產(chǎn)生刺激促進(jìn)植物生長(zhǎng)的物質(zhì)[37].酸桿菌門在沙地樟子松天然林中相對(duì)豐度最高,這與前人研究不一致[38-39].沙地樟子松天然林與土壤pH值有顯著相關(guān)性,且該生境土壤pH值在4種生境中酸性最強(qiáng),多數(shù)研究認(rèn)為酸桿菌是一種土壤寡營(yíng)養(yǎng)菌,其生長(zhǎng)與土壤酸堿性密切相關(guān)[40-41].沙地樟子松天然林中土壤細(xì)菌群落豐富度和多樣性指數(shù)最高,且表層土有較高的營(yíng)養(yǎng)性,酸桿菌門起到降解植物纖維素的作用[42].不同生境土壤理化性質(zhì)與土壤細(xì)菌群落結(jié)構(gòu)有顯著的相關(guān)性[43-44],其中,裸沙地細(xì)菌群落結(jié)構(gòu)與土壤含水量有極強(qiáng)的相關(guān)性,土壤含水量影響土壤微生物呼吸,細(xì)菌群落結(jié)構(gòu)對(duì)土壤含水量變化敏感[45-46].

    3.3 生境對(duì)土壤細(xì)菌功能基因的影響

    在呼倫貝爾沙區(qū)四種典型生境的一級(jí)功能層中,因?yàn)榧?xì)菌群落的功能特征會(huì)影響細(xì)菌群落結(jié)構(gòu)與多樣性,沙地樟子松天然林中較高豐度的的代謝、遺傳信息處理、有機(jī)系統(tǒng)3類功能基因使土壤細(xì)菌代謝旺盛,生長(zhǎng)力好,從而提高了細(xì)菌群落結(jié)構(gòu)的多樣性[47].二級(jí)功能層的功能基因存在差異性,表現(xiàn)為樟子松天然林0~10cm土層中核苷酸代謝、酶家族、氨基酸代謝、膜運(yùn)輸?shù)墓δ芑蜇S度顯著高于其他樣本.功能基因代謝主要作用是從土壤中吸收營(yíng)養(yǎng)物質(zhì),通過(guò)攝取氨基酸、能量、碳水化合物等來(lái)保證細(xì)菌的存活,膜運(yùn)輸?shù)墓δ芑蚩梢匀芙忤F和小分子,保證細(xì)菌快速成活[48].樟子松天然林通過(guò)這些功能基因從土壤中攝取更多的核苷酸、氨基酸,提高了細(xì)菌多樣性.功能基因的豐度通過(guò)影響微生物過(guò)程進(jìn)而影響生態(tài)過(guò)程的轉(zhuǎn)化[49].

    4 結(jié)論

    4.1 呼倫貝爾沙區(qū)不同生境土壤細(xì)菌多樣性指數(shù)存在顯著性差異,且沙地樟子松天然林土壤細(xì)菌多樣性最高,沙地樟子松人工林中土壤細(xì)菌多樣性最低.其Shannon指數(shù)分別為8.623±0.193、7.432±0.028.

    4.2 呼倫貝爾沙區(qū)不同生境土壤細(xì)菌群落組成及優(yōu)勢(shì)細(xì)菌門的豐度存在差異.草地、沙地樟子松人工林、沙地樟子松天然林3種生境中變形菌門相對(duì)豐度最高,其均值分別為29.83%±1.14%、34.73%±1.99%、31.95%±0.21%,裸沙地的主要優(yōu)勢(shì)細(xì)菌門為放線菌門,其相對(duì)豐度的均值為26.13%±0.43%.

    4.3 呼倫貝爾沙區(qū)不同生境土壤細(xì)菌功能較為豐富.四種生境中土壤細(xì)菌群落在環(huán)境信息處理、代謝、遺傳信息處理、有機(jī)系統(tǒng)方面表現(xiàn)活躍.

    [1] 厲桂香,馬克明.北京東靈山樹(shù)線處土壤細(xì)菌的PICRUSt基因預(yù)測(cè)分析[J]. 生態(tài)學(xué)報(bào), 2018,38(6):2180-2186.Li G X, Ma K M. PICRUSt-based predicted metagenomic analysis of treeline soil bacteria on Mount Dongling, Beijing [J]. Acta Ecologica Sinica, 2018,38(6):2180-2186.

    [2] 李 婷,張 威,劉光琇,等.荒漠土壤微生物群落結(jié)構(gòu)特征研究進(jìn)展[J]. 中國(guó)沙漠, 2018,38(2):329-338. Li T, Zhang W, Liu G Y, et al. Advances in the study of microbial ecology in desert soil [J]. Desert of China, 2018,38(2):329-338.

    [3] 黃 藝,黃木柯,柴立偉,等.干旱半干旱區(qū)土壤微生物空間分布格局的成因[J]. 生態(tài)環(huán)境學(xué)報(bào), 2018,27(1):191-198. Huang Y, Huang M K, Chai L W, et al. Drivers of the spatial patterns of soil microbial communities in arid and semi-arid regions [J]. Ecology and Environmental Sciences, 2018,27(1):191-198.

    [4] 李 新,焦 燕,代 鋼,等.內(nèi)蒙古河套灌區(qū)不同鹽堿程度的土壤細(xì)菌群落多樣性[J]. 中國(guó)環(huán)境科學(xué), 2016,36(1):249-260. Li X, Jiao Y, Dai G, et al. Soil bacterial community diversity under different degrees of saline-alkaline in the Hetao Area of Inner Mongolia [J]. China Environmental Science, 2016,36(1):249-260.

    [5] Suleiman A A, Tojo Soler C M, Hoogenboom G. Determining FAO- 56 crop coefficients for peanut under different water stress levels [J]. Irrigation Science, 2013,31(2):169-178.

    [6] 程 亮,王 信,郭青云.青藏高原不同生境土壤細(xì)菌群落結(jié)構(gòu)特征及其與環(huán)境的關(guān)系[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2019,37(1):18-26. Cheng L, Wang X, Guo Q Y. Relationship of soil bacterial community composition from different ecosystems on Qinghai-Tibet Plateau with environment factors [J]. Agricultural Research in the Arid Areas, 2019,37(1):18-26.

    [7] Peter H, Sommaruga R. Shift in diversity and function of lake bacteria communities upon glacier retreat [J]. The ISME Journal, 2016,10(7): 1545-1554.

    [8] Mohammad B, Falk H, Forslund S K, et al. Structure and function of the global topsoil microbiome [J]. Nature, 2018,560:233-237.

    [9] Wei H, Peng C H, Yang B, et al. Contrasting soil bacterial community, diversity, and function in two forests in China [J]. Frontiers in Microbiology, 2018,9:1693.

    [10] Maestre F T, Delgado-Baquerizo M, Jeffries T C, et al. Increasing aridity reduces soil microbial diversity and abundance in global drylands [J]. Proceedings of the National Academy of Sciences, 2015,112(51):15684-15689.

    [11] Langille M G, Zaneveld J, Caporaso J G, et al. Predictive functional profiling of microbial communities using 16S rRNA 15marker gene sequences [J]. Nature Biotechnology, 2013,31(9):814-821.

    [12] Asshauer K P, Wemheuer B, Daniel R, et al. Tax4 Fun: Predicting functional profiles from metagenomic 16S rRNA data [J]. Bioinformatics, 2015,31(17):2882-2884.

    [13] Louca S, Jacques S M S, Pires A P F, et al. High taxonomic variability despite stable functional structure across microbial communities [J]. Nature Ecology & Evolution, 2016,1(1):0015.

    [14] 董志穎,洪 慢,胡晗靜,等.過(guò)量氮輸入對(duì)寡營(yíng)養(yǎng)海水細(xì)菌群落代謝潛力的影響[J]. 環(huán)境科學(xué)學(xué)報(bào), 2018,38(2):457-466. Dong Z Y, Hong M, Hu H J, et al. Effect of excess nitrogen loading on the metabolic potential of the bacteria community in oligotrophic coastal water [J]. Acta Scientiae Circumstantiae, 2018,38(2):457-466.

    [15] 孫 峰,田 偉,張 菲,等.丹江口庫(kù)區(qū)庫(kù)濱帶植被土壤細(xì)菌群落多樣性及PICRUSt功能預(yù)測(cè)分析[J]. 環(huán)境科學(xué), 2019,40(1):421-429. Sun F, Tian W, Zhang F, et al. Composition and predictive functional analysis of rhizosphere bacterial communities in riparian buffer strips in the Danjiangkou reservoir, China [J]. Environmental Science, 2019,40(1):421-429.

    [16] Langille M G I, Zaneveld J, Caporaso J G, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences [J]. Nature Biotechnology, 2013,31(9):814-821.

    [17] Chu H, Sun H, Tripathi B M, et al. Bacterial community dissimilarity between the surface and subsurface soils equals horizontal differences over several kilometers in the western Tibetan Plateau [J]. Environmental Microbiology, 2016,18(5):1523-1533.

    [18] Wendu R L, Li G, Zhang J N, et al. The study of soil microbial biomass and soil enzyme avtivity on different grassland in Hulunbeier, Inner Mongolia [J]. Acta Prataculturae Sinica, 2010,19(5):94-102.

    [19] Zhang H F, Song X L, Wang C L, et al. The effects of different vegetation restoration patterns on soil bacterial diversity for sandy land in Hulunbeier [J]. Acta Ecologica Sinica, 2013,33(4):211-216.

    [20] 查同剛.土壤理化分析[M]. 北京:中國(guó)林業(yè)出版社, 2017:18-50. Zha T G. Soil physicochemical analysis [M]. Beijing: China Forestry Publishing House, 2017:18-50.

    [21] Rintala A, Pietil? S, Munukka E, et al. Gut microbiota analysis results are highly dependent on the 16S rRNA gene target region, whereas the impact of DNA extraction is minor [J]. Biomol Tech, 2017,28(1): 19-30.

    [22] Kuczynski J, Stombaugh J, Walters W A, et al. Using QIIME to analyze 16S rRNA gene sequences from microbial communities [J]. Current Protocols in Bioinformatics, 2011,Chapter 10(7):Unit 10.7.

    [23] 武志華,夏冬雙,王雪寒,等.利用PCR-DGGE技術(shù)分析內(nèi)蒙古西部地區(qū)土壤細(xì)菌的多樣性[J]. 生態(tài)學(xué)報(bào), 2019,39(7):2545-2557. Wu Z H, Xia D S, Wang X H, et al. Analysis of bacterial diversity of soil in western Inner Mongolia using PCR-DGGE [J]. Acta Ecologica Sinica, 2019,33(7):2545-2557.

    [24] Martirosyan V, Unc A, Miller G, et al. Desert perennial shrubs shape the microbial-community miscellany in laimosphere and phyllosphere Space [J]. Microbial Ecology, 2016,72(3):659-668.

    [25] 劉 洋,曾全超,黃懿梅.基于454高通量測(cè)序的黃土高原不同喬木林土壤細(xì)菌群落特征[J]. 中國(guó)環(huán)境科學(xué), 2016,36(11):3487-3494. Liu Y, Zeng Q C, Huang Y M. Soil microbial communities by 454 prosequencing under different arbor forests on the Loess Plateau [J]. China Environmental Science, 2016,36(11):3487-3494.

    [26] 張博文,楊彥明,李金龍,等.連續(xù)深松對(duì)黑土水熱酶特性及細(xì)菌群落的影響[J]. 生態(tài)學(xué)雜志, 2018,37(11):3323-3332. Zhang B W, Yang Y M, Li J L, et al. Effects of continuous subsoiling on temperature, water content, enzyme activity and bacterial community in black soil [J]. Chinese Journal of Ecology, 2018,37(11): 3323-3332.

    [27] 韓世忠,高 人,李愛(ài)萍,等.中亞熱帶地區(qū)兩種森林植被類型土壤微生物群落結(jié)構(gòu)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015,26(7):2151-2158. Han S Z, Gao R, Li A P, et al. Soil microbial community structure of two types of forests in the mid-subtropics of China [J]. Chinese Journal of Applied Ecology, 2015,26(7):2151-2158.

    [28] Mackay J E, Macdonald L M, Smernik R J, et al. Organic amendments as phosphorus fertilisers: chemical analyses, biological processes and plant puptake [J]. Soil Biology and Biochemistry, 2017,107:50-59.

    [29] 納小凡,鄭國(guó)琦,彭 勵(lì),等.不同種植年限寧夏枸杞根際微生物多樣性變化[J]. 土壤學(xué)報(bào), 2016,53(1):241-252. Na X F, Zhen G Q, Peng L, et al. Microbial biodiversity in rhizosphere ofL. relative to cultivation history [J]. Acta Pedologica Sinica, 2016,53(1):241-252.

    [30] 張凱煜,谷 潔,王小娟,等.微生物有機(jī)肥對(duì)櫻桃園土壤細(xì)菌群落的影響[J]. 中國(guó)環(huán)境科學(xué), 2019,39(3):1245-1252. Zhang K Y, Gu J, Wang X J, et al. Effects of bio-organic fertilizer on the soil bacterial community in a cherry orchard [J]. China Environmental Science, 2019,39(3):1245-1252.

    [31] Saleemi M, Kiani M Z, Sultan T, et al. Integrated efect of plant growth-promoting rhizobacteria and phosphate-solubilizing microorganisms on growth of wheat (L.) under rainfed condition [J]. Agriculture and Food Security, 2017,6(1):46.

    [32] 羅 旦,陳吉祥,程 琳,等.陜北沙區(qū)3種主要植物根際土壤細(xì)菌多樣性與土壤理化性質(zhì)相關(guān)性分析[J]. 干旱區(qū)資源與環(huán)境, 2019, 33(3):151-157. Luo D, Chen J X, Cheng L, et al. Analysis of bacteria diversity in the rhizosphere soil of three main plants and its correlation with the soil physical and chemical properties in the desertification area of northern Shanxi [J]. Journal of Arid Land Resources and Environment, 2019,33(3):151-157.

    [33] 劉 洋,黃懿梅,曾全超.黃土高原不同植被類型下土壤細(xì)菌群落特征研究[J]. 環(huán)境科學(xué), 2016,37(10):3931-3938. Liu Y, Huang Y M, Zeng Q C. Soil bacterial communities under different vegetation types in the Loess Plateau [J]. Environmental Science, 2016,37(10):3931-3938.

    [34] 高雪峰,韓國(guó)棟,張國(guó)剛.短花針茅荒漠草原土壤微生物群落組成及結(jié)構(gòu)[J]. 生態(tài)學(xué)報(bào), 2017,37(15):5129-5136. Gao X F, Han G D, Zhang G G. Soil microbial community structure and composition of Stipa Breiflora on the desert steppe [J]. Acta Ecologica Sinica, 2017,37(15):5129-5136.

    [35] 劉 睿,周孝德.水沙環(huán)境變化對(duì)季節(jié)性多沙河流沉積物菌群特征的影響[J]. 中國(guó)環(huán)境科學(xué), 2017,37(11):4342-4352. Liu R, Zhou X D. Influence of sediment characteristics and overlying water quality on sediment bacterial communities in a seasonal sandy river [J]. China Environmental Science, 2017,37(11):4342-4352.

    [36] Lin Y T, Jangid K, Whitman W B, et al. Soil bacterial communities in native and regenerated perhumid montane forests [J]. Applied Soil Ecology, 2011,47(2):111-118.

    [37] Stevenson A, Hallsworth J E. Water and temperature relations of soil Actinobacteria [J]. Environmental Microbiology Reports, 2014,6(6): 744-755.

    [38] 楊安娜,陸云峰,張俊紅,等.杉木人工林土壤養(yǎng)分及酸桿菌群落結(jié)構(gòu)變化[J]. 林業(yè)科學(xué), 2019,55(1):119-127. Yang A N, Lu Y F, Zhang J H, et al. Changes in soil nutrients and Acidobacteria Community Structure in Cunninghamia lanceolata Plantations [J]. Scientia Silvae Sinicae, 2019,55(1):119-127.

    [39] Bachar A, Soares M I M, Gillor O. The effect of resource islands on abundance and diversity of bacteria in arid soils [J]. Microbial Ecology, 2012,63(3):694-700.

    [40] Nie M, Meng H, Li K, et al. Comparison of bacterial and fungal communities between natural and planted pine forests in subtropical China [J]. Current Microbiology, 2012,64(1):34-42.

    [41] 王光華,劉俊杰,于鎮(zhèn)華,等.土壤酸桿菌門細(xì)菌生態(tài)學(xué)研究進(jìn)展[J]. 生物技術(shù)通報(bào), 2016,32(2):14-20. Wang G H, Liu J J, Yu Z H, et al. Research progress of Acidobacteria ecology in soils [J]. Biotechnology Bulletin, 2016,32(2):14-20.

    [42] Xiong J B, Liu Y Q, Lin X G, et al. Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau [J]. Environmental Microbiology, 2012,14(9):2457-2466.

    [43] Liu G Z, Jin M, Cai C T, et al. Soil microbial community structure and physicochemical properties in amomum tsaoko -based agroforestry systems in the Gaoligong Mountains, Southwest China [J]. Sustainability, 2019,11(2):1-14.

    [44] 劉旻霞,李 瑞,張 燦,等.蘭州市南山季節(jié)性土壤微生物特征及影響因素[J]. 中國(guó)環(huán)境科學(xué), 2018,38(7):2722-2730. Liu M X, Li R, Zhang C, et al. Seasonal characteristics and influencing factors of soil microbial in Nanshan, Lanzhou [J]. China Environmental Science, 2018,38(7):2722-2730.

    [45] Beth F T, Brockett, Cindy E, et al. Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada [J]. Soil Biology and Biochemistry, 2012,44(1):9-20.

    [46] 白玉濤,周 玉,趙 吉.內(nèi)蒙古高原干涸湖泊反硝化及甲烷氧化細(xì)菌的群落分析[J]. 中國(guó)環(huán)境科學(xué), 2012,32(7):1293-1301. Bai Y T, Zhou Y, Zhao J. community structures of denitrifying bacteria and methanotrophs in wetland soils of dry-up lake in the Inner Mongolia Plateau [J]. China Environmental Science, 2012,32(7): 1293-1301.

    [47] Kanokratana P, Uengwetwanit T, Rattanachomsri U, et al. Insights into the phylogeny and metabolic potential of a primary tropical peat swamp forest microbial community by metagenomic analysis [J]. Microbial Ecology, 2011,61(3):518-528.

    [48] Wu Z X, Hao Z P, Sun Y Q, et al. Comparison on the structure and function of the rhizosphere microbial community between healthy and root-rot Panax notoginseng [J]. Applied Soil Ecology, 2016,107:99- 107.

    [49] Fabina S, Paula Jorge L M, Rodrigues Jizhong Z, et al. Land use change alters functional gene diversity, composition and abundance in Amazon forest soil microbial communities [J]. Molecular Ecology, 2014,23(12):2988-2999.

    Soil bacteria community structure and function prediction in the Hulun Buir Sandy Area.

    DU Yu-jia1,2,GAO Guang-lei1,2*, CHEN Li-hua1, DING Guo-dong1,2, ZHANG Ying1, CAO Hong-yu1,2

    (1.Key Laboratory of State Forestry and Grassland Administration on Soil and Water Conservation, School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China;2.Yanchi Ecology Research Station of the Mu Us Desert, Yanchi 751500, China)., 2019,39(11):4840~4848

    Soil samples were collected from four habitats in the Hulun Buir sandy area including bare sandland, grassland,varplantation andnatural forest. Subsequently, 16S rRNA genes high-throughput sequencing and PICRUSt-based functional prediction were performed to detect soil bacterial community structure and potential functional component. The results indicated that in the Hulun Buir sandy area, the natural forest and plantation had the highest and lowest soil bacterial Shannon Index (8.623±0.193 and 7.432±0.028), respectively. There were significant differences in alpha and beta diversity of soil bacteria in different habitats. The relative abundance of Proteobacteria was the highest in the grassland (29.83%±1.14%),plantation (34.73%±1.99%) and natural forest (31.95%±0.21%).The relative abundance of Actinobacteriawas the highest in the bare sand (26.13%±0.43%). The dominant soil bacteria genera with highest relative abundance were, RB41with the mean values of 5.29%±2.24% and 4.22%±1.23%, respectively. Soil bacterial functions were classified into 6 and 40 functional categories at hierarchy level 1and 2,which implied the abundant soil bacteria functions. Soil bacteria was active in environmental information processing, metabolism, genetic information processing and organic systems. Further, the functional genes of soil bacterial from the natural forest was abundant in nucleotide metabolism, enzyme family, amino acid metabolism and carbohydrate metabolism, which ensured soil bacteria survival with higher diversity. Conclusively, the functional genes fluctuation of soil bacteria associated with different habitats in the Hulun Buir sandy area reflected the changes of soil bacterial community structure and diversity, and indicated the effects of functional genes on soil bacterial community. Our study will provide a firm basis for better prediction and understanding of soil bacteirial metabolic potential and functions in sandy area.

    Hulun Buir sandy area;habitat;community structure;16S rRNA;functional prediction

    X172

    A

    1000-6923(2019)11-4840-09

    杜宇佳(1996-),女,山西忻州人,北京林業(yè)大學(xué)碩士研究生,主要從事荒漠生態(tài)學(xué)研究.

    2019-04-25

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2018YFC0507101);國(guó)家自然科學(xué)基金項(xiàng)目(31600583);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2017PT03, 2015ZCQ-SB-02)

    * 責(zé)任作者, 副教授, gaoguanglei@bjfu.edu.cn

    猜你喜歡
    樟子松天然林沙地
    能鉆過(guò)柔軟沙地的蛇形機(jī)器人
    軍事文摘(2021年18期)2021-12-02 01:27:58
    山西省天然林保護(hù)工程二期建設(shè)成效與經(jīng)驗(yàn)探討
    呼倫貝爾沙地實(shí)現(xiàn)良性逆轉(zhuǎn)
    沙地迷宮
    天津:全面停止天然林商業(yè)性采伐
    天保工程:把天然林都保護(hù)起來(lái)
    天然林保護(hù)任重而道遠(yuǎn)
    風(fēng)滾草
    塞罕壩樟子松幼林撫育與管理
    初探北方樟子松栽培關(guān)鍵技術(shù)
    久久久久久伊人网av| 中文字幕免费在线视频6| 久久婷婷青草| 成人亚洲精品一区在线观看| 国产一区二区在线观看日韩| 综合色丁香网| 国产亚洲一区二区精品| 高清视频免费观看一区二区| 欧美日韩视频精品一区| 亚洲av电影在线观看一区二区三区| 成人无遮挡网站| 蜜桃国产av成人99| av有码第一页| 亚洲美女黄色视频免费看| 精品久久久久久久久av| 简卡轻食公司| 一级毛片aaaaaa免费看小| 亚洲第一区二区三区不卡| 亚洲国产av影院在线观看| 国产黄片视频在线免费观看| 亚洲精品中文字幕在线视频| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| √禁漫天堂资源中文www| 2021少妇久久久久久久久久久| 自线自在国产av| 精品少妇黑人巨大在线播放| 欧美日韩精品成人综合77777| 青春草亚洲视频在线观看| 九色亚洲精品在线播放| 热99久久久久精品小说推荐| 国产精品偷伦视频观看了| 午夜激情久久久久久久| 最近2019中文字幕mv第一页| 熟女av电影| 亚洲av中文av极速乱| a级毛片免费高清观看在线播放| 中文乱码字字幕精品一区二区三区| 在线观看人妻少妇| 精品熟女少妇av免费看| 日韩av免费高清视频| 在现免费观看毛片| 久久精品熟女亚洲av麻豆精品| 国产成人精品久久久久久| 99久久中文字幕三级久久日本| 大片免费播放器 马上看| 在线亚洲精品国产二区图片欧美 | 久久久久久久精品精品| 久久久久久久国产电影| 色94色欧美一区二区| 最新中文字幕久久久久| 插阴视频在线观看视频| 男女无遮挡免费网站观看| 夜夜看夜夜爽夜夜摸| 七月丁香在线播放| 久久久久精品性色| 亚洲精品一二三| 国产黄色免费在线视频| 亚洲国产精品一区三区| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 高清av免费在线| 男人爽女人下面视频在线观看| 成人综合一区亚洲| 免费高清在线观看视频在线观看| 啦啦啦中文免费视频观看日本| 亚洲成人一二三区av| 免费少妇av软件| 日本黄大片高清| 国产一区亚洲一区在线观看| 免费黄色在线免费观看| 亚洲精品美女久久av网站| 亚洲精品日韩av片在线观看| av国产久精品久网站免费入址| 日本91视频免费播放| 国产男女超爽视频在线观看| 2018国产大陆天天弄谢| 欧美3d第一页| 国产精品一区二区三区四区免费观看| 日日撸夜夜添| 国产精品99久久99久久久不卡 | 色视频在线一区二区三区| 性色av一级| 色5月婷婷丁香| 亚洲美女黄色视频免费看| 久久99热这里只频精品6学生| 欧美亚洲日本最大视频资源| 人人妻人人澡人人爽人人夜夜| www.色视频.com| 最后的刺客免费高清国语| 亚洲少妇的诱惑av| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 91精品国产九色| 十八禁网站网址无遮挡| 一本—道久久a久久精品蜜桃钙片| 免费观看的影片在线观看| 国产免费又黄又爽又色| 只有这里有精品99| √禁漫天堂资源中文www| 狂野欧美激情性bbbbbb| 精品人妻熟女毛片av久久网站| 69精品国产乱码久久久| 亚洲精品国产av成人精品| 在线观看人妻少妇| 色94色欧美一区二区| 久久免费观看电影| 精品人妻熟女毛片av久久网站| 赤兔流量卡办理| 欧美激情极品国产一区二区三区 | a级毛片在线看网站| 精品久久国产蜜桃| 狂野欧美激情性bbbbbb| av视频免费观看在线观看| 亚洲av欧美aⅴ国产| 亚洲国产欧美在线一区| 国产精品嫩草影院av在线观看| 久久久久久久亚洲中文字幕| 狠狠婷婷综合久久久久久88av| 一本—道久久a久久精品蜜桃钙片| 高清欧美精品videossex| 亚洲国产毛片av蜜桃av| 免费观看av网站的网址| 爱豆传媒免费全集在线观看| 国产老妇伦熟女老妇高清| 男人添女人高潮全过程视频| 成人毛片60女人毛片免费| 色婷婷av一区二区三区视频| xxxhd国产人妻xxx| 国内精品宾馆在线| 亚洲人成网站在线播| 尾随美女入室| 免费av不卡在线播放| 青春草视频在线免费观看| 欧美变态另类bdsm刘玥| 日本-黄色视频高清免费观看| 又大又黄又爽视频免费| 青春草亚洲视频在线观看| 韩国av在线不卡| 国产成人午夜福利电影在线观看| 制服丝袜香蕉在线| av又黄又爽大尺度在线免费看| 卡戴珊不雅视频在线播放| 热re99久久精品国产66热6| 国产成人午夜福利电影在线观看| 我的老师免费观看完整版| av福利片在线| 中国国产av一级| 水蜜桃什么品种好| 18+在线观看网站| 国产av码专区亚洲av| 最新中文字幕久久久久| 国产精品久久久久久av不卡| 一级黄片播放器| 亚洲少妇的诱惑av| 少妇的逼水好多| 熟女电影av网| 91aial.com中文字幕在线观看| 国国产精品蜜臀av免费| 午夜激情av网站| av在线老鸭窝| 我的老师免费观看完整版| 最新中文字幕久久久久| 国产高清国产精品国产三级| av卡一久久| 特大巨黑吊av在线直播| 蜜桃在线观看..| 欧美精品人与动牲交sv欧美| 色吧在线观看| 午夜福利视频精品| 亚洲精品中文字幕在线视频| 午夜视频国产福利| 22中文网久久字幕| 亚州av有码| 一本久久精品| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| 最新的欧美精品一区二区| 亚洲欧洲日产国产| 五月天丁香电影| 成人综合一区亚洲| 少妇被粗大猛烈的视频| 欧美日韩av久久| 亚洲国产日韩一区二区| 午夜久久久在线观看| 久热久热在线精品观看| a级毛片在线看网站| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 日本与韩国留学比较| 一级毛片电影观看| 少妇的逼好多水| 日本色播在线视频| 777米奇影视久久| 午夜免费观看性视频| 永久网站在线| 少妇人妻 视频| 黄色怎么调成土黄色| 一级黄片播放器| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 久久国产精品男人的天堂亚洲 | 色5月婷婷丁香| av在线观看视频网站免费| 欧美少妇被猛烈插入视频| 五月天丁香电影| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 国产精品人妻久久久影院| 国产精品久久久久久av不卡| 成年av动漫网址| 国产成人精品在线电影| 国产成人av激情在线播放 | 少妇被粗大猛烈的视频| 欧美精品亚洲一区二区| 欧美三级亚洲精品| 丝袜美足系列| 日本黄色日本黄色录像| 丰满乱子伦码专区| 黑人巨大精品欧美一区二区蜜桃 | 夫妻午夜视频| 国产亚洲精品第一综合不卡 | 午夜免费鲁丝| 国产精品国产av在线观看| 亚洲国产精品专区欧美| 超碰97精品在线观看| 男女免费视频国产| 王馨瑶露胸无遮挡在线观看| 一个人看视频在线观看www免费| 各种免费的搞黄视频| 超碰97精品在线观看| 寂寞人妻少妇视频99o| 天堂8中文在线网| 国产一区有黄有色的免费视频| 国产精品成人在线| 午夜激情av网站| 免费av不卡在线播放| 久久久久人妻精品一区果冻| 青春草视频在线免费观看| 看非洲黑人一级黄片| 狂野欧美激情性bbbbbb| 久久99一区二区三区| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 欧美国产精品一级二级三级| 午夜影院在线不卡| 免费人妻精品一区二区三区视频| 一区二区三区四区激情视频| 极品少妇高潮喷水抽搐| 亚洲精品aⅴ在线观看| 2018国产大陆天天弄谢| 久久精品久久久久久久性| 国产成人一区二区在线| 国产国语露脸激情在线看| 久久韩国三级中文字幕| 日本欧美视频一区| 蜜臀久久99精品久久宅男| 亚洲av电影在线观看一区二区三区| 在线免费观看不下载黄p国产| 蜜桃在线观看..| 大香蕉久久成人网| 在线观看美女被高潮喷水网站| 蜜桃国产av成人99| 人体艺术视频欧美日本| 18禁裸乳无遮挡动漫免费视频| 2021少妇久久久久久久久久久| 国产成人精品一,二区| 在线看a的网站| 777米奇影视久久| 国产片特级美女逼逼视频| 中文字幕久久专区| 日本与韩国留学比较| 天堂8中文在线网| 如何舔出高潮| 精品熟女少妇av免费看| 丰满乱子伦码专区| 国产欧美日韩一区二区三区在线 | 一区二区三区四区激情视频| 亚洲av在线观看美女高潮| .国产精品久久| 国产成人免费观看mmmm| 国内精品宾馆在线| 老司机影院毛片| 观看美女的网站| 精品人妻在线不人妻| 成人国语在线视频| 亚洲少妇的诱惑av| 亚洲怡红院男人天堂| 国产乱来视频区| 天堂中文最新版在线下载| 久久韩国三级中文字幕| 国产在线免费精品| 中国国产av一级| 久久狼人影院| av又黄又爽大尺度在线免费看| 蜜桃久久精品国产亚洲av| 一区二区三区精品91| 国产片特级美女逼逼视频| 大陆偷拍与自拍| 久久97久久精品| 国产伦理片在线播放av一区| 制服人妻中文乱码| 久久国产精品男人的天堂亚洲 | 99九九在线精品视频| 欧美人与善性xxx| 久久久精品94久久精品| 精品人妻在线不人妻| 爱豆传媒免费全集在线观看| 午夜福利视频在线观看免费| 亚洲国产av新网站| 99热这里只有精品一区| 国产一区二区在线观看日韩| 韩国高清视频一区二区三区| 大香蕉久久成人网| 久久99热6这里只有精品| 边亲边吃奶的免费视频| 亚洲精品日本国产第一区| 亚洲成人av在线免费| 久久99热这里只频精品6学生| 亚洲精品一二三| 精品久久久久久久久av| 在线观看免费视频网站a站| 成人国语在线视频| 国产高清三级在线| 99久久中文字幕三级久久日本| 熟妇人妻不卡中文字幕| 亚洲欧美成人综合另类久久久| 99热这里只有是精品在线观看| av网站免费在线观看视频| 日本欧美国产在线视频| 99国产综合亚洲精品| av在线播放精品| 欧美激情 高清一区二区三区| 国产亚洲午夜精品一区二区久久| 午夜免费观看性视频| 777米奇影视久久| 丁香六月天网| 中国三级夫妇交换| 亚洲av成人精品一区久久| 日日撸夜夜添| 免费黄频网站在线观看国产| 亚州av有码| 在线观看人妻少妇| 各种免费的搞黄视频| 97在线人人人人妻| 色吧在线观看| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 欧美激情极品国产一区二区三区 | 国产一区亚洲一区在线观看| 国产 精品1| 人妻制服诱惑在线中文字幕| 18禁裸乳无遮挡动漫免费视频| 一区二区三区精品91| 99久久综合免费| 国产精品99久久久久久久久| 国产精品人妻久久久影院| 蜜桃国产av成人99| 三级国产精品欧美在线观看| 精品少妇久久久久久888优播| 精品久久久噜噜| av福利片在线| 丝袜脚勾引网站| 女人久久www免费人成看片| 日韩av免费高清视频| 亚洲美女搞黄在线观看| 国产精品 国内视频| 少妇猛男粗大的猛烈进出视频| 99久久精品一区二区三区| 久久久久精品久久久久真实原创| 欧美日韩av久久| .国产精品久久| 搡老乐熟女国产| videos熟女内射| 亚洲熟女精品中文字幕| 国产精品久久久久成人av| 久久av网站| 色哟哟·www| 大又大粗又爽又黄少妇毛片口| 在线观看三级黄色| 美女中出高潮动态图| 在线观看三级黄色| 午夜免费观看性视频| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 伦理电影免费视频| 男女无遮挡免费网站观看| 亚洲美女视频黄频| 国产亚洲一区二区精品| 在线免费观看不下载黄p国产| 伊人久久精品亚洲午夜| 一级a做视频免费观看| 熟妇人妻不卡中文字幕| 边亲边吃奶的免费视频| 一区二区三区乱码不卡18| 日日摸夜夜添夜夜添av毛片| 日韩三级伦理在线观看| 久久久久久久久久久丰满| 成人手机av| 午夜激情福利司机影院| 国产成人精品久久久久久| 王馨瑶露胸无遮挡在线观看| 尾随美女入室| 大码成人一级视频| 欧美日韩亚洲高清精品| 婷婷色麻豆天堂久久| 亚洲欧美清纯卡通| 人妻制服诱惑在线中文字幕| 一级毛片电影观看| 欧美97在线视频| 18禁裸乳无遮挡动漫免费视频| 午夜精品国产一区二区电影| 在线观看免费日韩欧美大片 | 欧美精品高潮呻吟av久久| 少妇丰满av| 简卡轻食公司| 人体艺术视频欧美日本| 一二三四中文在线观看免费高清| 考比视频在线观看| 国产亚洲欧美精品永久| 亚洲美女搞黄在线观看| 日本色播在线视频| 少妇熟女欧美另类| 久久久久久久久久久久大奶| 大香蕉97超碰在线| 成人亚洲精品一区在线观看| 观看美女的网站| 精品视频人人做人人爽| 国产精品国产三级国产av玫瑰| 国产一级毛片在线| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 国产男女超爽视频在线观看| 久久精品久久精品一区二区三区| 两个人免费观看高清视频| 在线观看三级黄色| 欧美成人精品欧美一级黄| av在线观看视频网站免费| 亚洲人与动物交配视频| 大话2 男鬼变身卡| 午夜老司机福利剧场| 女性生殖器流出的白浆| 久久综合国产亚洲精品| 成年女人在线观看亚洲视频| 亚洲欧洲国产日韩| 性色av一级| 永久免费av网站大全| 久久青草综合色| 成人综合一区亚洲| 爱豆传媒免费全集在线观看| 婷婷成人精品国产| 日韩av不卡免费在线播放| 夜夜骑夜夜射夜夜干| 久久久久国产网址| 国产精品女同一区二区软件| 日韩一区二区视频免费看| 亚洲精品色激情综合| 中文字幕亚洲精品专区| 欧美日韩视频精品一区| 少妇 在线观看| 日本-黄色视频高清免费观看| 嫩草影院入口| 人妻一区二区av| 国产精品免费大片| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 香蕉精品网在线| 免费人妻精品一区二区三区视频| 欧美老熟妇乱子伦牲交| 男女免费视频国产| 女性被躁到高潮视频| 亚洲欧美一区二区三区黑人 | 欧美+日韩+精品| 国产国拍精品亚洲av在线观看| 一级爰片在线观看| 女人精品久久久久毛片| 内地一区二区视频在线| 制服诱惑二区| 韩国高清视频一区二区三区| 亚洲精品亚洲一区二区| 午夜老司机福利剧场| 亚洲精品久久成人aⅴ小说 | 最新中文字幕久久久久| 一区二区三区乱码不卡18| 午夜视频国产福利| 少妇丰满av| a级毛片黄视频| 99国产精品免费福利视频| 搡女人真爽免费视频火全软件| 一本一本综合久久| 国产成人一区二区在线| 精品一区二区免费观看| a级毛片黄视频| 国产探花极品一区二区| av免费观看日本| 亚洲四区av| 久久精品国产亚洲av涩爱| 我的女老师完整版在线观看| 国产精品99久久99久久久不卡 | 久久久久久人妻| 久久热精品热| 国产亚洲欧美精品永久| 日产精品乱码卡一卡2卡三| 欧美xxⅹ黑人| videossex国产| 色94色欧美一区二区| 热re99久久国产66热| 亚洲av二区三区四区| 黄色毛片三级朝国网站| 亚洲久久久国产精品| 国产综合精华液| 国产成人freesex在线| 久久精品久久久久久噜噜老黄| 久久久国产一区二区| 亚洲三级黄色毛片| 有码 亚洲区| 91成人精品电影| 久久久久精品久久久久真实原创| 久久久久久久国产电影| 国产伦精品一区二区三区视频9| 美女视频免费永久观看网站| 国产精品成人在线| 亚洲精品456在线播放app| 免费播放大片免费观看视频在线观看| 伊人久久国产一区二区| 欧美bdsm另类| 日日啪夜夜爽| www.色视频.com| 亚洲精品一区蜜桃| 日韩不卡一区二区三区视频在线| 狠狠精品人妻久久久久久综合| 日本wwww免费看| 亚洲精品一区蜜桃| 免费久久久久久久精品成人欧美视频 | 久久久午夜欧美精品| 麻豆精品久久久久久蜜桃| 91国产中文字幕| 亚洲精品456在线播放app| 成人国产麻豆网| 新久久久久国产一级毛片| 久久久久精品久久久久真实原创| 午夜福利视频在线观看免费| a级毛片免费高清观看在线播放| 精品亚洲成a人片在线观看| 亚洲一级一片aⅴ在线观看| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 久热这里只有精品99| 中国国产av一级| 国产探花极品一区二区| 日本av手机在线免费观看| 国产成人精品福利久久| 精品卡一卡二卡四卡免费| 狠狠婷婷综合久久久久久88av| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 国产乱来视频区| 日韩制服骚丝袜av| 日本欧美视频一区| 亚洲精品色激情综合| av视频免费观看在线观看| 国产日韩一区二区三区精品不卡 | 免费高清在线观看视频在线观看| 在线亚洲精品国产二区图片欧美 | 99久久人妻综合| 国产成人精品婷婷| 日韩精品免费视频一区二区三区 | 哪个播放器可以免费观看大片| 久久久久久久国产电影| 极品人妻少妇av视频| 欧美97在线视频| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 亚州av有码| 欧美成人午夜免费资源| 免费观看a级毛片全部| 亚洲欧美日韩另类电影网站| 国产精品一区二区三区四区免费观看| 国产成人av激情在线播放 | 日韩,欧美,国产一区二区三区| 观看美女的网站| 考比视频在线观看| 国产一区亚洲一区在线观看| 黑人巨大精品欧美一区二区蜜桃 | 在线看a的网站| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| 男女边吃奶边做爰视频| 人妻 亚洲 视频| 亚洲国产av新网站| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 九色亚洲精品在线播放| 国产成人精品久久久久久| 午夜老司机福利剧场| 欧美日韩综合久久久久久| 插阴视频在线观看视频| 大陆偷拍与自拍| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 一级二级三级毛片免费看| 欧美一级a爱片免费观看看| 亚洲av.av天堂| 自线自在国产av| 欧美日韩精品成人综合77777| 99热全是精品| 男女啪啪激烈高潮av片| 五月天丁香电影| 国产精品久久久久久久电影| 精品卡一卡二卡四卡免费| 亚洲中文av在线| 亚洲精品久久午夜乱码| 国产毛片在线视频|