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

    內(nèi)蒙古河套灌區(qū)不同鹽堿程度的土壤細(xì)菌群落多樣性

    2016-06-07 06:21:45楊銘德溫慧洋內(nèi)蒙古師范大學(xué)化學(xué)與環(huán)境科學(xué)學(xué)院內(nèi)蒙古呼和浩特010022
    中國(guó)環(huán)境科學(xué) 2016年1期
    關(guān)鍵詞:主成分分析細(xì)菌

    李 新,焦 燕,代 鋼,楊銘德,溫慧洋 (內(nèi)蒙古師范大學(xué)化學(xué)與環(huán)境科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010022)

    ?

    內(nèi)蒙古河套灌區(qū)不同鹽堿程度的土壤細(xì)菌群落多樣性

    李 新,焦 燕*,代 鋼,楊銘德,溫慧洋 (內(nèi)蒙古師范大學(xué)化學(xué)與環(huán)境科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010022)

    摘要:采用變性梯度凝膠電泳(DGGE)技術(shù)對(duì)內(nèi)蒙古河套灌區(qū)三種不同鹽堿程度土壤(鹽土,強(qiáng)度鹽化土,輕度鹽化土)不同深度(0~20cm和20~30cm)土壤細(xì)菌16S rDNA V3~V6 可變區(qū)擴(kuò)增片段進(jìn)行分析,并對(duì)土壤理化性質(zhì)進(jìn)行了測(cè)定.結(jié)果表明:細(xì)菌群落多樣性隨土壤鹽堿化程度的加深而減少(輕度鹽化土>強(qiáng)度鹽化土>鹽土),隨土壤深度的增加而降低(細(xì)菌群落多樣性0~20cm土層大于20~30cm土層).細(xì)菌Shannon-Wiener指數(shù)在輕度鹽化土中最大為3.36,在強(qiáng)度鹽化土和鹽土分別為3.05和2.49.不同鹽堿程度土壤以細(xì)菌相似系數(shù)聚類(lèi),分為0~20cm層與20~30cm層兩大族群,土壤細(xì)菌群落Shannon-Wiener指數(shù)在0~20cm層中(鹽土為3.04,強(qiáng)度鹽化土為3.29,輕度鹽化土為3.36)均大于在20~30cm層(鹽土為2.49,強(qiáng)度鹽化土為3.05,輕度鹽化土為3.14).相關(guān)性分析和典范對(duì)應(yīng)分析表明,土壤w(EC)、pH值、w(SOC)、w(TP)是土壤細(xì)菌群落結(jié)構(gòu)多樣性的顯著影響因素,不同鹽堿程度土壤中細(xì)菌群落的Shannon-Wiener指數(shù)與土壤w(EC)(r=-0.542,P<0.05)、pH(r=-0.526,P<0.05)呈顯著負(fù)相關(guān),與土壤w(SOC)(r=0.700,P<0.01)和w(TP)(r=0.805,P<0.01)呈極顯著正相關(guān).w(EC)和pH對(duì)鹽堿土壤細(xì)菌群落結(jié)構(gòu)的影響力最大.回收DGGE圖譜中20個(gè)優(yōu)勢(shì)條帶進(jìn)行測(cè)序分析,結(jié)果顯示,變形菌綱(α-變形菌綱、β-變形菌綱、γ-變形菌綱和δ-變形菌綱)是鹽堿土壤的主要類(lèi)群.

    關(guān)鍵詞:河套灌區(qū);鹽堿土壤;細(xì)菌;變性梯度凝膠電泳(DGGE);主成分分析

    * 責(zé)任作者, 教授, jiaoyan@imnu.edu.cn

    土壤微生物是土壤有機(jī)組分和生態(tài)系統(tǒng)中最活躍的部分,在促進(jìn)土壤質(zhì)量和植物健康方面發(fā)揮著重要的作用,被認(rèn)為是最敏感的土壤質(zhì)量生物學(xué)指標(biāo)[1].土壤微生物多樣性指土壤生態(tài)系統(tǒng)中所有的微生物種類(lèi)、它們擁有的基因以及這些微生物與環(huán)境之間相互作用的多樣化程度.一般認(rèn)為,土壤微生物多樣性存在于基因、物種、種群以及群落等4個(gè)層面,是土壤生態(tài)系統(tǒng)的一個(gè)基本生命特征,也是時(shí)間和空間的函數(shù)[2-3].國(guó)內(nèi)外對(duì)農(nóng)田、森林、草原的土壤微生物研究較多[4?5],但對(duì)鹽堿土壤的研究多局限于理化特性的測(cè)定分析,國(guó)內(nèi)針對(duì)鹽堿土壤的微生物研究也僅限于天津?yàn)I海區(qū)域,江蘇濱海鹽堿地,松嫩平原,新疆、甘肅[6]等地,康貽軍等[7]研究灘涂鹽堿土壤微生物生態(tài)特征,曹?chē)?guó)棟等[8]對(duì)扇緣帶鹽堿土上生長(zhǎng)的三種植被類(lèi)型下的土壤微生物類(lèi)群數(shù)量進(jìn)行了研究.張廣志等[9]在鹽堿土壤鑒定出假單胞菌屬和黃桿菌屬.牛世全等[10]對(duì)河西走廊鹽堿土細(xì)菌種群多樣性的16S rDNA研究表明,γ-變形菌為優(yōu)勢(shì)菌群.石偉等[11]對(duì)極端鹽堿土壤細(xì)菌的分離篩選及抗鹽特性研究也發(fā)現(xiàn)多數(shù)菌株屬于γ-變形菌門(mén).然而,針對(duì)內(nèi)蒙古河套灌區(qū)鹽堿土壤中的微生物研究數(shù)量、種群結(jié)構(gòu)、優(yōu)勢(shì)菌系以及鹽害與土壤微生物活動(dòng)之間的生態(tài)關(guān)系等相關(guān)報(bào)道較少.

    土壤鹽堿化已嚴(yán)重制約著內(nèi)蒙古河套灌區(qū)的農(nóng)業(yè)生產(chǎn),同時(shí)對(duì)本地區(qū)的糧食安全構(gòu)成了嚴(yán)重威脅[12-13].河套灌區(qū)鹽堿地面積約4.3×105hm2,巴彥淖爾市占2.4×105hm2,其土壤鹽漬化面積與比例如下:輕度鹽化面積1.3×105hm2,占52.83%;中度鹽化面積7.6×104hm2,占31.94%;重度鹽化面積3.6×104hm2,占15.23%[14].

    土壤中最活躍生物因子是土壤細(xì)菌,它既是土壤微生物的重要組成部分,占土壤微生物總數(shù)的70%~90%,也是土壤物質(zhì)流和能量流的主要推動(dòng)者,可以直接反映土壤肥力,已被公認(rèn)為土壤生態(tài)系統(tǒng)變化的預(yù)警及敏感指標(biāo)[15-19],在土壤養(yǎng)分循環(huán)中起著至關(guān)重要的作用.由于從環(huán)境樣品中分離和培養(yǎng)細(xì)菌困難,不經(jīng)分離培養(yǎng)步驟,直接從土壤中提取總DNA并分析其基因信息的分子生物學(xué)方法已被發(fā)展用來(lái)描述和鑒定微生物群落.近年來(lái)基于DNA方法的群落分析如PCR擴(kuò)增技術(shù)、克隆文庫(kù)法、熒光原位雜交法、限制性酶切片段長(zhǎng)度多態(tài)性法,變性和溫度梯度凝膠電泳法[20],高通量測(cè)序技術(shù)[21]等得到了迅速的發(fā)展其中變性梯度凝膠電泳(DGGE)技術(shù)由于有結(jié)果精確、可靠性強(qiáng)、重現(xiàn)率高等優(yōu)點(diǎn),已廣泛應(yīng)用于不同生態(tài)環(huán)境中微生物群落組成、多樣性及種群動(dòng)態(tài)變化的研究[22].如謝學(xué)軍等[23]借助DGGE技術(shù)研究了重金屬污染物對(duì)土壤微生物多樣性的影響;Dong等[24]利用DGGE分析方法研究了濕地土壤微生物多樣性;Smalla等[25]采用DGGE技術(shù)對(duì)不同植被根際微生物進(jìn)行了微生物多樣性分析.

    土壤耕作層一般疏松多孔,通透性良好,根系集中,而土壤犁底層孔隙度小,通氣性差,透水性不良.二者的差異勢(shì)必會(huì)對(duì)微生物群落結(jié)構(gòu)等產(chǎn)生影響.該文采用DGGE技術(shù),對(duì)內(nèi)蒙古河套灌區(qū)不同鹽堿程度不同深度土壤(0~20cm耕作層和20~30cm犁底層)中細(xì)菌的群落結(jié)構(gòu)和多樣性差異進(jìn)行分析,探討土壤理化特性和鹽堿化土壤微生物的相互作用機(jī)制,對(duì)于鹽堿地的綜合利用與開(kāi)發(fā)具有重要的實(shí)踐指導(dǎo)意義.

    1 材料與方法

    1.1 研究區(qū)概況

    供試土壤采于內(nèi)蒙古巴彥淖爾市烏拉特前旗,地處黃河北岸,河套平原東端,該地屬溫帶大陸性氣候,冬寒而長(zhǎng),夏熱而短,干旱少雨,春季風(fēng)沙較大,極端最高氣溫為39.7℃ ,最低氣溫-30.7℃ ,年平均氣溫7.7℃ .年平均日照3212.5h,無(wú)霜期167d.降水集中于7~9月,年平均降水量213.5mm,最大降水量為8月,極端日降水量達(dá)109.6mm.

    1.2 樣品采集

    為避免地形等因素影響,研究樣區(qū)選擇平坦地形并按照鄰近原則布置.試驗(yàn)于2014年5月未種植作物前(前茬季節(jié)作物葵花)選取3種不同鹽堿程度的鹽堿土壤,樣地S1、S2和S3(10m×10m)采用“S”型布點(diǎn)法采集樣方在距離地表0~20cm和20~30cm用土鉆分層取樣,每個(gè)樣方設(shè)10個(gè)樣點(diǎn)(1m×1m).各采樣點(diǎn)重復(fù)取樣3次,并將3次所得土壤樣品充分混勻,將可見(jiàn)的植物殘?bào)w(如根、莖和葉)和土壤動(dòng)物去除,裝于無(wú)菌聚乙烯自封袋中.一份風(fēng)干磨碎過(guò)2mm篩用于測(cè)定土壤理化特性,另一份迅速運(yùn)回實(shí)驗(yàn)室,然后將部分土樣分裝于若干無(wú)菌離心管中-80℃保存.試驗(yàn)土壤樣地基本情況見(jiàn)表1.

    表1 不同土壤鹽分含量(%)Table 1 Salt content in different types of soil in Hetao irrigation area (%)

    1.3 土壤鹽堿程度分析

    內(nèi)蒙古河套灌區(qū)地處干旱、半干旱區(qū),該研究三種土壤的含鹽總量依次為S1(1.69%)>S2(0.83%)> S3(0.12%),S1,S2和S3都是SO42-和Cl-含量最多,見(jiàn)表1,結(jié)合土壤鹽化分級(jí)標(biāo)準(zhǔn)(表2)可知,S1土壤為鹽土,S2土壤為強(qiáng)度鹽化土壤,S3為輕度鹽化土壤.

    表2 土壤鹽化分級(jí)指標(biāo)[26]Table 2 Soil salinization classification index[26]

    1.4 樣品分析

    1.4.1 樣品理化性質(zhì)測(cè)定 pH值以1:2.5土水比用復(fù)合電極測(cè)定;電導(dǎo)率(EC)以1:5土水比用復(fù)合電極測(cè)定;土壤容重的測(cè)定用環(huán)刀法;土壤質(zhì)地用比重計(jì)速測(cè)法;土壤有機(jī)碳(SOC)用重鉻酸鉀容量法-外加熱法;土壤全磷(TP)用HClO4-H2SO4測(cè)定,土壤全氮(TN)用凱氏定氮法.土壤鹽分:碳酸根、重碳酸根離子測(cè)定用電位滴定法;氯離子測(cè)定用硝酸銀滴定法;硫酸根離子測(cè)定用EDTA容量法;鈣、鎂離子測(cè)定用EDTA容量法;鉀、鈉離子測(cè)定用火焰光度法.理化性質(zhì)具體測(cè)定方法按照《土壤農(nóng)化分析(第三版)》進(jìn)行[27].1.4.2 土壤總DNA的提取及細(xì)菌16S rDNA片段的PCR擴(kuò)增 分別取0.5g土樣樣品,采用土壤DNA提取試劑盒(E.Z.N.A?.Soil DNA kit, Omega, USA)提取土壤總DNA.以樣品基因組DNA為模板,采用細(xì)菌通用正向引物338F在5‘末端加40個(gè)堿基的GC夾: 5c-CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG G ACTCCTACG GGA GGC AGC AG-3‘;和反向引物是518R:5‘-ATT ACC GCG GCT GCT GG –3’,40個(gè)GC夾加到338F前面,以保證DGGE實(shí)驗(yàn)的穩(wěn)定和片斷的分離,這一對(duì)引物能夠擴(kuò)增絕大多數(shù)的細(xì)菌.

    PCR擴(kuò)增體系(50μL):10×PCR buffer 5μL; dNTP(2.5mmol/L)3.2μL;rTaq(5U/μL)0.4μL;GC-3 38F(20μmol/L)1μL;518R(20μmol/L)1μL;模板DNA 50ng;補(bǔ)ddH2O至50μL.PCR擴(kuò)增程序:94℃預(yù)變性5min;94℃變性1min,55℃復(fù)性45s,72℃延伸1min,30個(gè)循環(huán);最終72℃延伸10min.PCR產(chǎn)物采用OMEGA公司DNA Gel Extraction Kit純化回收.PCR儀為Biometra公司生產(chǎn)的T-gradient,凝膠成像儀為Bio-Rad公司的Gel-Doc2000凝膠成像系統(tǒng).

    1.4.3 PCR產(chǎn)物的變性梯度凝膠電泳(DGGE)分析 取10μL PCR的產(chǎn)物進(jìn)行變性梯度凝膠電泳(DGGE)分析.采用變性梯度為35%~55%、濃度為8%的聚丙烯酰胺凝膠(化學(xué)變性劑為100%尿素7mol/L和40%(V/V)的丙烯酰胺)在1×TAE緩沖液中150V 60℃下電泳5h.變性梯度凝膠電泳(DGGE)完畢后、采用銀染法染色、步驟如下:首先固定液(乙醇50mL、冰醋酸2.5mL、定容500mL)固定15min;Milli-Q純水清洗、20s和2min各一次;其次銀染液(硝酸銀1g、37%甲醛0.75mL、定容500mL)染色15min;Milli-Q純水清洗、20s和2min各一次;再次顯色液(氫氧化鈉7.5g、37%甲醛2.5mL、定容500mL)顯色5-7min;最后用終止液(乙醇50mL、冰醋酸2.5mL、定容500mL)終止反應(yīng).

    1.4.4 DGGE圖譜中優(yōu)勢(shì)條帶的回收與測(cè)序用滅菌的手術(shù)刀切下待回收DGGE條帶,采用OMEGA公司Poly-Gel DNA Extraction Kit回收目的條帶.以2 μL回收產(chǎn)物為模板,338F/518R為引物進(jìn)行PCR擴(kuò)增.PCR擴(kuò)增體系(50 μL) 為:10×PCR buffer 5μL;dNTP(2.5mmol/L)3.2μL; rTaq(5U/μL)0.4μL;338F(20mmol/L)1μL;518R(20 mmol/L)1μL;模板DNA 1μL;補(bǔ)ddH2O至50μL.PCR擴(kuò)增程序?yàn)?94℃預(yù)變性4min;94℃變性30s,55℃復(fù)性30s,72℃延伸30s,30個(gè)循環(huán);最后,72℃延伸10min.將重新擴(kuò)增的DNA片段切膠回收、純化后,連接到Pmd18-T載體上,并轉(zhuǎn)化至DH5α感受態(tài)細(xì)胞中,篩選陽(yáng)性克隆,進(jìn)行序列測(cè)定.

    1.5 數(shù)據(jù)分析方法

    采用SPSS 22.0統(tǒng)計(jì)軟件進(jìn)行方差分析.細(xì)菌多樣性指數(shù)是研究群落物種數(shù)和個(gè)體數(shù)以及均勻度的綜合指標(biāo).根據(jù)電泳圖譜中樣品條帶數(shù)目及每個(gè)條帶的強(qiáng)度(灰度),對(duì)各樣品中細(xì)菌多樣性指數(shù)(Shannon-Wiener,H)、均勻度(Evenness, E)和豐富度(Richness,S)等指標(biāo)進(jìn)行分析.DGGE圖譜采用Quantity one軟件對(duì)每個(gè)樣品的電泳條帶數(shù)目、條帶密度進(jìn)行數(shù)字化分析,多樣性指數(shù)(H)、豐富度指數(shù)(S)和均勻度指數(shù)(E)等指標(biāo)被用來(lái)比較不同樣品的多樣性情況.其算法如下所示:

    式中:pi為樣品中單一條帶的強(qiáng)度在該樣品所有條帶總強(qiáng)度中所占的比率;N為DGGE圖譜單一泳道上所有條帶的豐富度;Ni為第i條帶的豐富度;S是某樣品中所有條帶數(shù)目總和.測(cè)序結(jié)果采用DNAstar和Cluster軟件進(jìn)行序列分析,下載最相似的菌株序列作為系統(tǒng)發(fā)育樹(shù)的參考序列.然后采用MEGA軟件,Neighbor-joining法構(gòu)建系統(tǒng)發(fā)育樹(shù),自展數(shù)(bootstrap)為1000.利用Canoco for Windows 4.5軟件進(jìn)行細(xì)菌群落和理化因子的典范對(duì)應(yīng)(canonical correspondence analysis CCA)分析.

    2 結(jié)果與分析

    2.1 土壤細(xì)菌群落多樣性分析

    2.1.1 土壤細(xì)菌群落DGGE圖譜 經(jīng)凝膠電泳檢測(cè),選取較高質(zhì)量土壤樣品細(xì)菌基因組DNA.各樣品的PCR擴(kuò)增產(chǎn)物經(jīng)純化后進(jìn)行細(xì)菌的DGGE電泳見(jiàn)圖1,利用Quantity One軟件對(duì)DGGE膠圖進(jìn)行條帶識(shí)別和圖譜分析,結(jié)果顯示,土壤樣品微生物的16S rDNA片段通過(guò)DGGE被分離出不同位置的若干條帶.根據(jù)膠圖上條帶遷移的位置和數(shù)目可以看出細(xì)菌的種類(lèi)豐富程度.樣品間的主要條帶不同,表明細(xì)菌優(yōu)勢(shì)種存在差別.不同泳道同一位置條帶的明暗程度(光密度值)則反映該細(xì)菌類(lèi)群在不同環(huán)境樣品中的相對(duì)豐度[28].不同土壤樣品在細(xì)菌種類(lèi)組成上差別明顯.

    鹽土、強(qiáng)度鹽化土和輕度鹽化土的表層(0~20cm)土壤中,對(duì)應(yīng)的泳道有許多相同的條帶,表明這些相同的條帶所代表的細(xì)菌種群為不同程度鹽堿地群落的表層土壤所共有,但是這些共有條帶的亮度不盡相同,表示它們?cè)诟髯喳}堿地群落表層土壤中的豐富度存在差異.此外,強(qiáng)度鹽化土和輕度鹽化土細(xì)菌群落的表層土壤有一些條帶,在鹽土中卻未見(jiàn),說(shuō)明強(qiáng)度鹽化土和輕度鹽化土表層土壤細(xì)菌群落與鹽土細(xì)菌群落有顯著差異,可能是鹽脅迫抑制了這些細(xì)菌的生長(zhǎng).相反,鹽土表層土壤中又有個(gè)別幾種細(xì)菌種群是其特有的,這可能與鹽土土壤中嗜鹽或耐鹽微生物有關(guān).在20~30cm層土壤中,不同鹽堿程度土壤類(lèi)型間的泳道條帶變異較大,輕度鹽化土細(xì)菌種群較鹽土和強(qiáng)度鹽化土少.同一鹽堿程度土壤在不同土壤深度的泳道條帶比較,可以看出各自的優(yōu)勢(shì)條帶明顯,說(shuō)明優(yōu)勢(shì)種群發(fā)生了變化.

    圖1 不同鹽堿程度不同土壤深度細(xì)菌群落DGGE圖譜Fig.1 DGGE fingerprint of bacteria from different soil layers in different saline-alkaline soil

    圖2 不同鹽堿程度不同深度土壤細(xì)菌群落聚類(lèi)分析Fig.2 The cluster analysis of bacterial community from different soil layers in different saline-alkaline soil數(shù)值代表各泳道與1泳道相比的相似性

    2.1.2 不同鹽堿程度不同深度土壤細(xì)菌群落相似性分析 依據(jù)樣品間的相似系數(shù)構(gòu)建的聚類(lèi)圖(UPGMA),如圖2所示.15個(gè)土壤樣品中,細(xì)菌群落組成總體分為兩大族群,3種類(lèi)型鹽堿土0~20cm層土壤細(xì)菌聚為一族群,20~30cm層土壤細(xì)菌聚為另一族群.說(shuō)明0~20cm層與20~30cm層土壤細(xì)菌種群的構(gòu)成變化較大.其中,0~20cm層土壤細(xì)菌又可分成兩大族群,鹽土0~20cm層的3個(gè)重復(fù)土壤細(xì)菌聚為一族群,強(qiáng)度鹽化土和輕度鹽化土土壤細(xì)菌聚為另一族群.這說(shuō)明鹽土的土壤細(xì)菌種群不同于強(qiáng)度鹽化土和輕度鹽化土的細(xì)菌種群.強(qiáng)度鹽化土中5#和6#的相似度為0.66,輕度鹽化土中8#和9#的相似度達(dá)0.61.一般認(rèn)為相似值高于0.60的兩個(gè)群體具有較好的相似性[27].強(qiáng)度鹽化土和輕度鹽化土0~20cm層細(xì)菌,他們的相似系數(shù)僅為0.43,說(shuō)明強(qiáng)度鹽化土與輕度鹽化土細(xì)菌種群結(jié)構(gòu)存在很大差異.

    2.1.3 土壤細(xì)菌群落多樣性指數(shù)分析 不同鹽堿程度不同深度土壤細(xì)菌群落多樣性指數(shù)如表3所示.不論在0~20cm層還是20~30cm層土壤中,輕度鹽化土和強(qiáng)度鹽化土的細(xì)菌Shannon-Wiener指數(shù)和Richness均為最大,并且與鹽土有顯著性差異.說(shuō)明輕度鹽化土和強(qiáng)度鹽化土不論在0~20cm層還是20~30cm層土壤中細(xì)菌群落都最為豐富,但其中各菌種分布的Evenness差異不顯著;隨著土壤深度的增加,Evenness不變,Shannon-Wiener指數(shù)和Richness指數(shù)均呈顯著降低.這表明從土壤0~20cm層至30cm的深度范圍內(nèi),隨著土層深度的增加,土壤細(xì)菌的Shannon-Wiener指數(shù)越小,Richness越低.

    表3 不同鹽堿程度不同深度土壤細(xì)菌群落多樣性指數(shù)分析Table 3 Analysis of bacterial genetic diversity in different soil layers in different saline-alkaline soil

    2.1.4 土壤細(xì)菌群落系統(tǒng)組成分析 DGGE凝膠條帶回收后,以338F/518R為引物進(jìn)行PCR擴(kuò)增,獲得約200bp的DNA片段.PCR產(chǎn)物純化后連接到pMD18-T載體上,轉(zhuǎn)化至DH5α感受態(tài)細(xì)胞中,篩選陽(yáng)性克隆測(cè)序.測(cè)序結(jié)果與GenBank中的序列進(jìn)行比對(duì),得到條帶所代表的細(xì)菌類(lèi)型.每個(gè)回收條帶選取3個(gè)克隆進(jìn)行了序列測(cè)定,如表4.根據(jù)這些DGGE帶譜所代表的微生物即可判定不同鹽堿程度土壤不同土層細(xì)菌群落的組成狀況.

    條帶Band4在輕度鹽化土0~20cm條帶多,屬于固氮弓菌屬(Azoarcus)為輕度鹽化土0~20cm層特有的優(yōu)勢(shì)菌.條帶Band11和Band36在鹽土0~20cm條帶多,屬于黃桿菌屬(Flavobacterium)為鹽土0~20cm層特有的優(yōu)勢(shì)菌.條帶Band23屬于纖維弧菌屬(Cellvibrio),條帶Band27屬于食堿菌屬(Alcanivorax),在3種不同程度鹽堿土0~20cm層都有出現(xiàn).條帶Band3屬于Rhodoligotrophos appendicifer為鹽化土20~30cm特有的優(yōu)勢(shì)菌.一般情況下認(rèn)為,兩種細(xì)菌16S rDNA 序列同源性小于98%時(shí)屬于不同種,而同源性小于93%~95%時(shí),則屬于不同的屬[28].所以條帶Band8、18、28、37、38和43屬于新菌的可能性比較大.

    將獲得的20條序列通過(guò)BLAST比對(duì)分析后,由表5可知,大致可分為5個(gè)大的細(xì)菌系統(tǒng)發(fā)育類(lèi)群,即變形菌門(mén)(Proteobacterium)、藍(lán)藻門(mén)(Cyanophyta)、綠彎門(mén)(Chloroflexi)、擬桿菌門(mén)(Bacteroidetes)和放線菌門(mén)(Actinobacteria)系統(tǒng).其中,變形菌門(mén)的12個(gè)條帶中有7個(gè)(Band1、12、15、19、23、24和27)在3種不同鹽堿度土壤0~20cm和20~30cm都出現(xiàn),屬于γ-變形菌綱(γ-proteobacteria),它是變形菌門(mén)最大的亞群,也是鹽堿土壤中的優(yōu)勢(shì)類(lèi)群;3個(gè)(Band3、6、34)屬于α-變形菌綱(α-Proteobacteria),1個(gè)(Band4)屬于β-變形菌綱(β-proteobacteria),1個(gè)(Band39)屬于δ-變形菌綱(δ–Deltaproteobacteria).γ-變形菌綱包括假單胞菌屬、纖維弧菌屬和食堿菌屬, α-變形菌綱包括Rhodoligotrophos appendicifer屬、新鞘脂菌屬和鞘脂單胞屬,β-變形菌綱包括固氮弓菌屬,δ-變形菌綱包括脫硫弧菌屬.藍(lán)藻門(mén)包括Crocosphaera watsonii屬,綠彎門(mén)包括Roseiflexus castenholzii,擬桿菌門(mén)包括黃桿菌屬和海藻菌屬,放線菌門(mén)包括Planobispora siamensis屬、庫(kù)茨勒菌屬和血桿菌屬.變形菌門(mén)(α-變形菌綱、β-變形菌綱、γ-變形菌綱和δ-變形菌綱)是鹽堿土壤的主要類(lèi)群.DGGE測(cè)定鹽堿土壤中的菌株大部分都是屬于耐鹽堿的細(xì)菌菌株.

    表4 測(cè)序克隆序列與其GenBank最相似序列的比對(duì)結(jié)果Table 4 Comparison between determined clone sequence and its most similar sequence in GenBank

    2.2 土壤理化特性對(duì)細(xì)菌群落結(jié)構(gòu)影響

    2.2.1 不同鹽堿程度不同深度土壤理化特性 由表5可知, 3種不同鹽堿度不同深度的土壤, w(EC)、pH值、w(TN)、w(TP)、w(SOC)、土壤水分和土壤容重的含量均存在差異.同一土壤深度,3種不同鹽堿程度土壤的pH和w(EC)表現(xiàn)為鹽土最大,并且與輕度鹽化土和強(qiáng)度鹽化土有顯著差異.土壤容重、w(SOC)、w(TP)和w(TN)表現(xiàn)為:鹽土<強(qiáng)度鹽化土;同一鹽堿程度土壤,不同土壤深度,土壤水分、w(EC)、土壤w(SOC)、w(TP)和w(TN)表現(xiàn)為:0~20cm層含量高于20~30cm層含量,與土壤容重變化趨勢(shì)相反.pH值在0~20cm層和20~30cm層土壤中則無(wú)明顯規(guī)律.

    表5 不同鹽堿程度不同深度土壤理化參數(shù)Table 5 The physical and chemical properties of different soil layers in different saline-alkaline soil

    2.2.2 相關(guān)性分析 表6為0~20cm層與20~ 30cm層土壤細(xì)菌群落基因型多樣性Shannon-Wiener 指數(shù),Evenness和Richness指數(shù)與表5七種理化參數(shù)的相關(guān)性分析結(jié)果.由表6可見(jiàn),3種不同鹽堿程度土壤細(xì)菌群落的Shannon-Wiener指數(shù)與Richness都與土壤w(EC)、pH值呈顯著負(fù)相關(guān),與土壤w(SOC)、w(TN)和 w(TP)呈極顯著正相關(guān).Evenness指數(shù)與土壤土壤水分和w(TN)呈顯著正相關(guān).w(EC)、pH值越小,土壤w(SOC)、w(TN)和w(TP)越大,土壤細(xì)菌群落的Shannon-Wiener和Richness指數(shù)越大.

    表6 不同鹽堿程度不同深度土壤細(xì)菌群落指數(shù)與理化參數(shù)的相關(guān)性分析Table 6 Pearson correlation analysis between different soil layers in different saline-alkaline soil bacterial community indices and physicochemical parameters

    圖3 不同鹽堿程度不同深度土壤細(xì)菌群落-環(huán)境CCA排序Fig.3 CCA biplot of samples and environmental variables1~3:鹽土0~20cm深度土壤3個(gè)平行樣,4~6:強(qiáng)度鹽化土0~20cm深度土壤3個(gè)平行樣,7~9:輕度鹽化土0~20cm深度土壤3個(gè)平行樣;10、11:鹽土20~30cm深度土壤2個(gè)平行樣,12、13:強(qiáng)度鹽化土20~30cm深度土壤2個(gè)平行樣,14、15:輕度鹽化土20~30cm深度土壤2個(gè)平行樣

    2.2.3 土壤細(xì)菌和理化性質(zhì)的CCA分析 將DGGE圖譜中各條帶的灰度值作為物種數(shù)據(jù),利用CANOCO進(jìn)行去趨勢(shì)對(duì)應(yīng)分析(DCA),所得結(jié)果中第1排序軸的梯度長(zhǎng)度值為2.012,當(dāng)該值介于1.5~3.0時(shí),選擇線性模型做冗余分析(RDA)或者選擇單峰模型做典范對(duì)應(yīng)(CCA)分析均合適.根據(jù)Legendre[31]分析理論,由于典范對(duì)應(yīng)分析(CCA)每一步計(jì)算結(jié)果都可以和環(huán)境因子進(jìn)行回歸,CCA排序能夠很好地反應(yīng)物種在特征因子梯度上的分布,因此選擇CCA更合適.將物種數(shù)據(jù)與0~20cm,20~30cm土壤理化參數(shù)w(EC)、pH值、w(TP)、w(TN)、w(SOC)、土壤水分w(W)和土壤容重做CCA,第1排序軸和第2排序軸解釋了土壤細(xì)菌群落的21.3%和15.5%,總排序軸即所選的土壤環(huán)境理化因子對(duì)物種的總解釋量為81.7%,損失的信息只有18.3%.這說(shuō)明土壤w(EC)、pH值、w(TP)、w(SOC)、土壤水分和土壤容重在研究鹽堿土土壤中的重要性.然而在鹽堿土壤中環(huán)境中可能還有其他理化因子對(duì)細(xì)菌群落結(jié)構(gòu)產(chǎn)生影響.

    圖3為不同鹽堿程度不同深度土壤細(xì)菌群落土壤-環(huán)境CCA排序圖.第1排序軸與w(EC)、pH值呈負(fù)相關(guān),相關(guān)系數(shù)分別為-0.7695、-0.7038,說(shuō)明橫坐標(biāo)從左到右w(EC)、pH值逐漸降低;第1排序軸與土壤容重呈正相關(guān),相關(guān)系數(shù)為0.7637,說(shuō)明從左到右土壤容重逐漸增大.第2排序軸與w(SOC)、w(TN)和w(TP)呈正相關(guān),相關(guān)系數(shù)分別為0.8871、0.8304和0.9659,說(shuō)明縱坐標(biāo)從下至上w(SOC)、w(TN)和w(TP)逐漸增大.w(EC)、pH值和土壤容重是影響鹽堿土壤中細(xì)菌群落結(jié)構(gòu)的主要影響因子.由圖3可見(jiàn),在0~20cm和20~30cm土層3種不同鹽堿程度土壤的分異明顯,而且各土壤平行樣一致性較好,均聚集在一起.w(EC)和pH值環(huán)境變量的長(zhǎng)度都長(zhǎng)于其他理化因子,也說(shuō)明其在鹽堿土壤中對(duì)細(xì)菌群落的影響力最大.樣方垂直投影于射線上,沿著變量箭頭方向環(huán)境變量值增大.輕度鹽化土和強(qiáng)度鹽化土0~20cm的細(xì)菌群落適宜在高w(SOC)、w(TN)和w(TP)中生存,鹽土0~20cm的細(xì)菌群落適宜在高w(EC)、pH值生存.在20~30cm土壤中,輕度鹽化土中的細(xì)菌群落在土壤容重較大時(shí)適宜生存.

    CCA分析結(jié)果表明,Monte Carlo測(cè)試的所有成分軸P <0.01,w(EC)、pH和土壤容重對(duì)鹽堿土壤的細(xì)菌群落存在顯著影響,并且w(EC)、pH值和土壤容重對(duì)鹽堿土壤細(xì)菌群落的影響力最大.

    3 討論

    在自然生態(tài)環(huán)境下的許多微生物不適宜在實(shí)驗(yàn)室培養(yǎng),因此用傳統(tǒng)的方法如平板計(jì)數(shù)法得到的微生物多樣性結(jié)果是非常片面的[32];碳素利用法(Biolog系統(tǒng))僅能鑒定快速生長(zhǎng)的部分微生物,只反應(yīng)了潛在的而不是原位的代謝多樣性[33];而以PCR-DGGE技術(shù)為代表的分子生態(tài)學(xué)方法,因?yàn)闊o(wú)需培養(yǎng)且檢驗(yàn)時(shí)間較短,相對(duì)較為先進(jìn)和可靠.PCR-DGGE能很好的揭示不同程度鹽堿土壤中相關(guān)細(xì)菌較高的種群多樣性、遺傳多樣性和生態(tài)多樣性.

    3.1 不同鹽堿程度土壤對(duì)細(xì)菌群落結(jié)構(gòu)的影響

    利用PCR-DGGE的研究方法得出,河套地區(qū)鹽堿土壤,在0~20cm層和20~30cm層土壤中,鹽土和強(qiáng)度鹽化土、輕度鹽化土細(xì)菌群落種類(lèi)之間差異顯著.隨著鹽堿化程度的加深,土壤細(xì)菌群落的Shannon-Wiener指數(shù)和豐富度指數(shù)都顯著性降低.這與孫佳杰等[34]研究一致,土壤鹽化程度越高,土壤微生物數(shù)量越少,土壤微生物中細(xì)菌的數(shù)量從大到小依次為輕度鹽化土、中度鹽化土、重度鹽化土、鹽土.本研究表明,土壤中細(xì)菌群落的Shannon-Wiener指數(shù)與Richness隨土層深度的增加而顯著降低.潘雪蓮等[35]研究表明, 在0~30cm土壤深度范圍內(nèi),黃土高原微生物多樣性指數(shù)隨土層深度的增加而減少;張社奇等[36]在研究黃土高原刺槐林地土壤微生物垂直分布時(shí)發(fā)現(xiàn)土壤細(xì)菌群落對(duì)深度的變化極為顯著.

    3.1.1 土壤鹽度和堿度對(duì)土壤對(duì)細(xì)菌群落結(jié)構(gòu)的影響 土壤中鹽堿程度的加深主要表現(xiàn)為土壤鹽分的增加和pH值升高.相關(guān)性分析表明,3種不同鹽堿程度的土壤中細(xì)菌群落的Shannon-Wiener指數(shù)與Richness與土壤w(EC)、pH值呈顯著負(fù)相關(guān).土壤鹽堿程度加深意味著土壤w(EC)和pH值越高,土壤細(xì)菌群落的多樣性和豐富度越低.土壤水溶性離子總量(簡(jiǎn)稱土壤全鹽量)是衡量土壤鹽害程度大小的重要指標(biāo),土壤水溶性離子又是強(qiáng)電解質(zhì),其導(dǎo)電能力可用電導(dǎo)率w(EC)表示,并且二者具有正相關(guān)性,所以土壤w(EC)越高,鹽分越高,鹽害程度越嚴(yán)重[37].土壤微生物對(duì)w(EC)變化十分敏感,當(dāng)土壤鹽分升高時(shí),鹽分脅迫直接改變了微生物的生存環(huán)境,造成土壤微生物滲透脅迫,抑制和降低土壤活性微生物種群數(shù)量[38].李新等在利用PLFA方法研究鹽堿土壤微生物時(shí)發(fā)現(xiàn),土壤鹽害程度越高,堿性越強(qiáng)土壤主要微生物多樣性越單一,反之則越豐富[39].康貽軍等[7],Yu等[40]指出,細(xì)菌數(shù)量與土壤全鹽含量呈顯著負(fù)相關(guān),土壤鹽害程度越高,微生物數(shù)量越少,土壤細(xì)菌和真菌的數(shù)量分布從大到小為輕度鹽化土、中度鹽化土、重度鹽化土和鹽土.許多研究證明,在多種生態(tài)系統(tǒng)中pH通常與細(xì)菌群落結(jié)構(gòu)有很好的相關(guān)性[41?42].

    3.1.2 土壤有機(jī)碳、全磷等理化性質(zhì)對(duì)土壤細(xì)菌群落結(jié)構(gòu)的影響 3種不同鹽堿程度的土壤中細(xì)菌群落的Shannon-Wiener指數(shù)與Richness都與土壤w(SOC)和w(TP)呈極顯著正相關(guān),與土壤容重呈顯著負(fù)相關(guān).即土壤鹽堿程度越低,土壤w(SOC)和w(TP)越高,土壤細(xì)菌群落的多樣性和豐富度也越高.這與許多學(xué)者研究一致,隨著鹽堿程度的增加,土壤養(yǎng)分含量逐漸下降,N、P等含量明顯降低[43-46].水分條件是影響土壤中微生物生存和活性的一個(gè)重要因素,適當(dāng)?shù)耐寥浪謼l件可以增加土壤微生物量和微生物活性[47].Griffiths等[48]研究指出,水分能明顯調(diào)節(jié)草地微生物多樣性,并會(huì)增加微生物本身對(duì)水的抗受性.而該研究表明,土壤水分w(W)與細(xì)菌Evenness指數(shù)呈顯著正相關(guān),和Shannon-Wiener指數(shù)并不存在相關(guān)關(guān)系,可能與微生物在河套灌區(qū)這種極端干旱和鹽堿脅迫的條件下,已經(jīng)適應(yīng)了缺水的環(huán)境,土壤中的其他營(yíng)養(yǎng)成分成為限制微生物活動(dòng)的主要因子.

    3.2 主要土壤細(xì)菌類(lèi)群的分析

    DGGE條帶經(jīng)序列測(cè)定后發(fā)現(xiàn),不同鹽堿程度土壤不同深度的優(yōu)勢(shì)菌群不同.變形菌門(mén)細(xì)菌(α-變形菌綱、β-變形菌綱、γ-變形菌綱和δ-變形菌綱)、擬桿菌門(mén)(黃桿菌屬)和綠彎門(mén)(玫瑰彎菌屬)在不同鹽堿度土壤中均有分布.張廣志[9]在鹽堿土壤中分離鑒定耐鹽細(xì)菌時(shí)發(fā)現(xiàn)了黃桿菌屬和假單胞菌屬,并且在鹽濃度高達(dá)25%時(shí),也能生存.有研究發(fā)現(xiàn),假單胞菌屬細(xì)菌是一種耐鹽堿菌群,大多數(shù)種不需要有機(jī)生長(zhǎng)因子,同時(shí)還可降解環(huán)境中廣泛存在的多種有機(jī)分子,在自然界和污染環(huán)境的礦化修復(fù)作用過(guò)程中有著重要地位.硫弧菌屬最適生長(zhǎng)溫度在25~30℃,是硫酸鹽還原細(xì)菌,主要以硫酸鹽或其他的氧化態(tài)硫化物作為呼吸鏈的末端電子受體,以氫氣為電子受體,還原產(chǎn)物硫化氫.玫瑰彎菌屬(Rosebacter Shiba)是一種在堿性熱溫泉中發(fā)現(xiàn)的細(xì)絲狀不產(chǎn)氧光合細(xì)菌[50].食堿菌屬只有在相對(duì)高堿性和高鹽量的鹽土和強(qiáng)度鹽化土0~20cm層分布.劉艷等[51]研究深海熱液區(qū)微生物對(duì)深海環(huán)境響應(yīng)機(jī)制時(shí)也發(fā)現(xiàn)了食堿菌屬,并證明了食堿菌屬具有適應(yīng)高溫和大量硫的代謝特征.鞘脂單胞菌屬大部分是鞘氨醇單胞菌,存在于輕度鹽化土0~30cm,鹽土和強(qiáng)度鹽化土的20~30cm層.它們能廣泛的生存在各種水體、大氣、土壤甚至地球上的極端環(huán)境中.它們還能降解比較復(fù)雜的大分子有機(jī)物,也能產(chǎn)生有價(jià)值的高分子,還能利用簡(jiǎn)單的小分子在極端貧瘠的環(huán)境中存活[52].

    綜上可知,不同鹽堿土壤中的主要細(xì)菌都是耐鹽堿的細(xì)菌菌株.變形菌綱(α-變形菌綱、β-變形菌綱、γ-變形菌和δ-變形菌綱)是鹽堿土壤的主要類(lèi)群.在所有測(cè)序條帶中,變形菌門(mén)占分析樣品總量的55%,其余為擬桿菌門(mén)(Bacteroidetes)、綠彎菌門(mén)(Chloroflexi)、放線菌門(mén)(Actinobacteria)和藍(lán)藻門(mén).牛世全等[8]對(duì)河西走廊鹽堿土細(xì)菌種群多樣性的16S rDNA研究發(fā)現(xiàn),原生鹽堿土中存在9大細(xì)菌類(lèi)群,γ-變形菌為優(yōu)勢(shì)菌群.石偉[11]對(duì)極端鹽堿土壤細(xì)菌的分離篩選及抗鹽特性研究顯示多數(shù)菌株屬于γ-變形菌綱(γproteobacteria).該研究表明內(nèi)蒙古河套灌區(qū)鹽堿土壤細(xì)菌具有豐富物種、遺傳和發(fā)育系統(tǒng)多樣性.16S rDNA序列分析技術(shù)能夠在屬以上水平很好的反應(yīng)細(xì)菌的分類(lèi).利用現(xiàn)代分子生物學(xué)技術(shù)分析細(xì)菌結(jié)構(gòu)的基礎(chǔ)上,進(jìn)一步對(duì)土壤微生物區(qū)系及其相互協(xié)同拮抗的代謝機(jī)制、以及與鹽堿環(huán)境之間的關(guān)系將成為今后研究的重點(diǎn).

    4 結(jié)論

    4.1 內(nèi)蒙古河套灌區(qū)鹽土、強(qiáng)度鹽化土和輕度鹽化土分別在0~20cm層和20~30cm層土壤細(xì)菌群落差異顯著,土壤鹽堿化程度越低,土壤細(xì)菌群落多樣性和豐富度越高,表現(xiàn)為輕度鹽化土>強(qiáng)度鹽化土>鹽土;0~20cm層與20~30cm層的土壤細(xì)菌群落也存在顯著性差異,隨著土層的加深,土壤細(xì)菌群落多樣性和豐富度減小.

    4.2 CCA分析表明, w(EC)、pH和土壤容重對(duì)鹽堿土壤細(xì)菌群落結(jié)構(gòu)的影響力最大.相關(guān)分析表明,土壤pH值、w(EC)越低,土壤w(SOC)、w(TN) 和w(TP)越高,土壤細(xì)菌群落的多樣性和豐富度也越高.

    4.3 變形菌門(mén)(α-變形菌綱、β-變形菌綱、γ-變形菌綱和δ-變形菌綱)是鹽堿土壤的主要類(lèi)群.鹽堿土壤中的菌株大部分都是屬于耐鹽堿的細(xì)菌菌株.

    參考文獻(xiàn):

    [1] Niemi R M, Heiskanen I, Wallenius K, et al.Extraction and purification of DNA in rhizosphere soil samples for PCR-DGGE analysis of bacterial consortia [J].Journal of Microbiological Methods, 2001,45(3):155-165.

    [2] 周德慶.微生物學(xué)教程 [M].北京:高等教育出版社, 1993.

    [3] 黃昌勇.土壤學(xué) [M].北京:中國(guó)農(nóng)業(yè)出版社, 2000.

    [4] 郭逍宇,董 志,宮輝力.再生水灌溉對(duì)草坪土壤微生物群落的影響 [J].中國(guó)環(huán)境科學(xué), 2006,26(4):482-485.

    [5] 侯憲文,吳建軍,徐建明.鉛-芐嘧磺隆對(duì)土壤微生物活性與群落結(jié)構(gòu)的影響 [J].中國(guó)環(huán)境科學(xué), 2007,27(6):738-742.

    [6] 李海云,牛世全,孔維寶,等.河西走廊石羊河下游地區(qū)鹽堿土中放線菌多樣性—以民勤縣為例 [J].中國(guó)環(huán)境科學(xué), 2015, 35(6):1805-1813.

    [7] 康貽軍,胡 建,董必慧.灘涂鹽堿土壤微生物生態(tài)特征的研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2007,26(S1):181-183.

    [8] 曹?chē)?guó)棟,陳接華,王紹明,等.不同鹽生植被類(lèi)型下土壤微生物特性研究 [J].新疆農(nóng)業(yè)科學(xué), 2012,49(3):523-530.

    [9] 張廣志,周紅姿,楊合同,等.鹽堿土壤中耐鹽細(xì)菌的分離與鑒定[J].山東農(nóng)業(yè)科學(xué), 2009,9:49-50.

    [10] 牛世全,景彩虹,廖世齊,等.河西走廊鹽堿土細(xì)菌種群結(jié)構(gòu)多樣性的研究 [J].西北師范大學(xué)學(xué)報(bào), 2013(2):90-95.

    [11] 石 偉,極端鹽堿土壤細(xì)菌的分離篩選及抗鹽特性研究 [D].哈爾濱:東北林業(yè)大學(xué), 2011.

    [12] 李鳳霞,郭永忠,許 興,等.鹽堿地土壤微生物生態(tài)特征研究進(jìn)展 [J].安徽農(nóng)業(yè)科學(xué), 2011,39(23):14065-14067.

    [13] 劉 鑫,魏占民,王長(zhǎng)生,等.基于ArcGIS的河套灌區(qū)土壤鹽堿化空間分析 [J].人民黃河, 2011,33(12):88-91.

    [14] 侯玉明,王 剛,王二英,等.河套灌區(qū)鹽堿土成因、類(lèi)型及有效的治理改良措施 [J].現(xiàn)代農(nóng)業(yè), 2011,(1):92-93.

    [15] Constancias F, Prevost-Boure NC, Terrat S, et al.Microscale evidence for a high decrease of soil bacterial density and diversity by cropping [J].Agronomy for Sustainable Development, 2014, 34(4):831-840.

    [16] Sugden A M.Ecology: diversity and ecosystem resilience [J].Science, 2000,290(5490):233-235.

    [17] 張 薇,魏海雷,高洪文,等.土壤微生物多樣性及其環(huán)境影響因子研究進(jìn)展 [J].生態(tài)學(xué)雜志, 2005,24(1):48-52.

    [18] Mann L, Tolbert V.Soil sustainability in renewable biomass plantings [J].Ambio, 2000,29(8):492-498.

    [19] Walter K D, Margaret K B, Courtney S C, et al.Biological properties of soil and subsurface sediments under abandoned pasture and cropland [J].Soil Biology Biochemistry, 1997,28(7): 837-846.

    [20] 馬悅欣, Holmstrom C, Webb J,等.變性梯度凝膠電泳(DGGE)在微生物生態(tài)學(xué)中的應(yīng)用 [J].生態(tài)學(xué)報(bào), 2003,23(8):1561-1569.

    [21] 秦 楠,栗東芳,楊瑞馥.高通量測(cè)序技術(shù)及其在微生物學(xué)研究中的應(yīng)用 [J].微生物學(xué)報(bào), 2011(4):445-457.

    [22] 李 丹,王秋玉.變性梯度凝膠電泳及其在土壤微生物生態(tài)學(xué)中的應(yīng)用 [J].中國(guó)農(nóng)學(xué)通報(bào), 2011,27(3):6-9.

    [23] 段學(xué)軍,黃春曉.重金屬鎘對(duì)水田土壤微生物基因多樣性的影響[J].應(yīng)用與環(huán)境生物學(xué)報(bào), 2008,14(4):510-513.

    [24] Dong X L, Reddy G B.Soil bacterial communities in constructed wetlands treated with swine wastewater using PCR-DGGE technique [J].Bioresource Technology, 2010,101(4):1175-1182.

    [25] Smalla K, Wieland G, Buchner A, et al.Bulk and rhizosphere,soil bacterial communities studied by denaturing gradient gel electrophoresis: plant-dependent enrichment and seasonal shifts revealed [J].Applied and Environmental Microbiology, 2001,67: 4742–4751.

    [26] 王遵親.中國(guó)鹽漬土 [M].北京:科學(xué)出版社, 1993.

    [27] 鮑士旦.土壤農(nóng)化分析 [M].3版.北京:中國(guó)農(nóng)業(yè)出版社, 2000.

    [28] Lyautey E, Lacoste B, Ten-Hage L, et al.Analysis of bacterial diversity in river biofilms using 16S r DNA PCR-DGGE: methodological settings and fingerprints interpretation [J].Water Research, 2005,39(2/3):380-388.

    [29] 陳法霖,張 凱,鄭 華,等.PCR-DGGE技術(shù)解析針葉和闊葉凋落物混合分解對(duì)土壤微生物群落結(jié)構(gòu)的影響 [J].應(yīng)用生物學(xué)報(bào), 2011,17(2):145-150.

    [30] Moschetti G, Aponte M, Blaiotta G, et al.Characterization of halophilic archaea isolated from different hypersaline ecosystems [J].Annals of microbiology, 2006,56(2):119-127.

    [31] Lehendre P, Gallagher E D.Ecologically meaningful transformations for ordination of species data [J].Oecologia, 2001,129(2):271-280.

    [32] 葉姜瑜,羅固源.微生物可培養(yǎng)性低的生態(tài)學(xué)釋因與對(duì)策 [J].微生物學(xué)報(bào), 2005,45(3):478-482.

    [33] Garland J L, Mills A L.Classification and characterization of heterotrophic microbial communities on the basis of patterns of community-level sole-carbon-source utilization [J].Applied and environmental microbiology, 1991,57(8):2351-2359.

    [34] 孫佳杰,尹建道,謝玉紅,等.天津?yàn)I海鹽堿土土壤微生物生態(tài)特征研究 [J].南京林大學(xué)學(xué)報(bào):自然科學(xué)版, 2010,34(3):57-61.

    [35] 潘雪蓮,黃 晟,方 昊,等.黃土高原土壤中細(xì)菌群落結(jié)構(gòu)多樣性的PCR-DGGE分析 [J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2009,25(3): 39-43.

    [36] 張社奇,王國(guó)棟,田 鵬,等.黃土高原刺槐林地土壤微生物的分布特征 [J].水土保持學(xué)報(bào), 2004,18(6):128?131.

    [37] 尹建道,孫佳杰,郝志強(qiáng),等.天津?yàn)I海地區(qū)土壤含鹽量與電導(dǎo)率的關(guān)系 [J].安徽農(nóng)業(yè)科學(xué), 2010,38(30):16882-16883.

    [38] Garcia-Gil J C, Plaza C, Soler-Rovira P, et al.Long-term effects of municipal solid waste compost application on soil enzyme activities and microbial biomass [J].Soil Biology and Biochemistry, 2000,(32):1907-1913.

    [39] 李 新,焦 燕,楊銘德.用磷脂脂肪酸(PLFA)譜圖技術(shù)分析內(nèi)蒙古河套灌區(qū)不同鹽堿程度土壤微生物群落多樣性 [J].生態(tài)科學(xué), 2014,33(3):488-494.

    [40] Yu S, He Z L, Huang C Y.Advances in their search of soil microorganisms and their mediated processes under heavy metal stress [J].The Journal of Applied Ecoloyg, 2003,14:618-622.

    [41] Christopher D C.Impact of cattle grazing and inorganic fertiliser additions to managed grasslands on the microbial communitycomposition of soils [J].Applied Soil Ecology, 2006,31(1/2):73?82.

    [42] Singh B K, Nunan N, Ridgway K P, et al.Relationship between assemblages of my corrhizal fungi and bacteria on grass roots [J].Environmental Microbiology, 2008,10(2):534?541.

    [43] Vega-Jarquin C, Garcia-Mendoza M, Jablonowski N, et al.Rapid immobilization of applied nitrogen in saline–alkaline soils [J].Plant and Soil, 2003,256:379–388.

    [44] Fernandez-Luqueno F, Marsch R, Espinosa-Victoria D, et al.Remediation of PAHs in a saline–alkaline soil amended with waste water sludge and the effect on dynamics of C and N [J] .Science of the Total Environment, 2008,402:18-28.

    [45] 張 杰,余 潮,王自海,等.不同植被群落表層土壤中細(xì)菌群落多樣性 [J].環(huán)境科學(xué)研究, 2013,26(8):866-872.

    [46] Araujo A S F, Borges C D, Tsai S M, et al.Soil bacterial diversity in degraded and restored lands of Northeast Brazil [J].Antonie van Leeuwenhoek, 2014,106(5):891-899.

    [47] Killham K.Soil ecology [M].Cambridge: Cambridge University Press, 1994.

    [48] Griffiths R I, Whitelely A S, O'donnell A G, et al.Physiological and community responses of established grassland bacterial populations to water Stress [J].Applied and Environmental Microbiology, 2003,69(12):6961-6968.

    [49] 邵明安,王全九,黃明斌.土壤物理學(xué) [M].北京:高等教育出版社, 2006.

    [50] Hanada S, Takaichi S, Matsuura K, et al.Roseiflexus castenholziigen.nov.sp.nov., a thermophilic, filamentous, photosynthetic bacterium that lacks chlorosome [J].International Journal of Systematic and Evolutionary Microbiology, 2002,52: 187-193.

    [51] 劉 艷.深海熱液區(qū)微生物的篩選鑒定及對(duì)深海環(huán)境的響應(yīng)機(jī)制 [D].山東:山東師范大學(xué), 2009.

    [52] 胡 杰,何曉紅,李大平,等.鞘氨醇單胞菌研究進(jìn)展 [J].應(yīng)用與環(huán)境生物學(xué)報(bào), 2007,13(12):431- 437.

    Soil bacterial community diversity under different degrees of saline-alkaline in the Hetao Area of Inner Mongolia.

    LI Xin, JIAO Yan*, DAI Gang, YANG Ming-de, WEN Hui-yang (Chemistry and Environment Science College, Inner Mongolia Normal University, Huhhot 010022, China).China Environmental Science, 2016,36(1):249~260

    Abstract:Denaturing gradient gel electrophoresis (DGGE) based on analyzing V3~V6 variable regions of bacterial 16S rDNA, were used to examine the microbial community structure and diversity in different saline-alkaline soils (saline soil, strongly salinized soil and slight salinized soil) and soil depths (0~20cm, 20~30cm) in the Hetao Area of Inner Mongolia.Their physical and chemical properties were measured, respectively.Experimental results showed that bacterial community structure and diversity decreased with the different saline-alkaline soil (slight salinized soil> strongly salinized soil> saline soil)and declined with soil depths (0~20cm>20~30cm).Shannon-Wiener index of bacterial community in slight salinized soil was the highest value (3.36), while it was only 3.05 and 2.49 in strongly salinized soil and saline soil.The cluster analysis of DGGE bands showed that the bacterial community structure and diversity formed two distinct clusters 0~20cm and 20~30cm.Shannon-Wiener index of bacterial community in 0~20cm layer (saline soil 3.04, strongly salinized soil 3.29, slight salinized soil 3.36) were higher than 20~30cm layer (saline soil2.49, strongly salinized soil 3.05, slight salinized soil 3.14).Analyses of Pearson correlation and CCA revealed that variations of bacterial community structure were mainly affected by contents of w(EC), pH, w(SOC), w(TP).The soil bacteria diversification index has a negative correlation with the w(EC)(r=-0.542,P<0.05)、pH(r=-0.526,P<0.05), and bacterial community diversity exhibited very significant positive correlation with w(SOC)(r=0.700,P<0.01) and w(TP)(r=0.805,P<0.01).w(EC) and pH were the the greatest influence in saline-alkaline soils.Twenty bands were excised from the DGGE gel and re-amplified for 16S rDNA sequencing.Based on the sequencing results, eleven bands can be identified as related to Proteobacteria.These results provide evidence that Proteobacteria are the domain bacterial communities in different saline-alkaline soilsbook=250,ebook=253of Hetao Area of Inner Mongolia.Saline-alkali soils were mostly bacterial strains belonging to tolerant salty.

    Key words:Hetao Area;saline-alkaline soils;bacterial;DGGE (denaturing gradient gel electrophoresis); principal component analysis

    中圖分類(lèi)號(hào):X172

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-6923(2016)01-0249-11

    收稿日期:2015-05-18

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(41165010,41375144,41565009);2013內(nèi)蒙古高等學(xué)?!扒嗄昕萍加⒉拧敝С钟?jì)劃資助(NJYT-13-B06);內(nèi)蒙古師范大學(xué)研究生科研創(chuàng)新基金項(xiàng)目(CXJJS13044)

    作者簡(jiǎn)介:李 新(1987-),女,內(nèi)蒙古烏蘭浩特人,在讀碩士研究生,主要從事土壤微生物多樣性研究.

    猜你喜歡
    主成分分析細(xì)菌
    偉大而隱秘的細(xì)菌
    細(xì)菌大作戰(zhàn)
    細(xì)菌怪來(lái)了
    細(xì)菌大作戰(zhàn)
    基于NAR模型的上海市房產(chǎn)稅規(guī)模預(yù)測(cè)
    主成分分析法在大學(xué)英語(yǔ)寫(xiě)作評(píng)價(jià)中的應(yīng)用
    江蘇省客源市場(chǎng)影響因素研究
    SPSS在環(huán)境地球化學(xué)中的應(yīng)用
    考試周刊(2016年84期)2016-11-11 23:57:34
    長(zhǎng)沙建設(shè)國(guó)家中心城市的瓶頸及其解決路徑
    服務(wù)貿(mào)易結(jié)構(gòu)優(yōu)化路徑研究
    热99国产精品久久久久久7| 成人影院久久| 97在线视频观看| 国产精品一国产av| 免费黄频网站在线观看国产| 波多野结衣av一区二区av| 成年女人在线观看亚洲视频| 久久精品熟女亚洲av麻豆精品| 少妇的逼水好多| 亚洲男人天堂网一区| 捣出白浆h1v1| 99精国产麻豆久久婷婷| 满18在线观看网站| 久久99精品国语久久久| 日韩伦理黄色片| 亚洲国产av影院在线观看| 久久精品国产自在天天线| 免费观看a级毛片全部| 欧美人与性动交α欧美精品济南到 | 免费观看在线日韩| 欧美xxⅹ黑人| 国产成人精品久久久久久| 亚洲一区二区三区欧美精品| 最新的欧美精品一区二区| 一区在线观看完整版| 国产一区二区三区av在线| 91aial.com中文字幕在线观看| 精品一区二区免费观看| 精品国产乱码久久久久久小说| 午夜福利网站1000一区二区三区| 考比视频在线观看| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区四区第35| 国产野战对白在线观看| 十八禁高潮呻吟视频| 国产在线一区二区三区精| 国产精品久久久av美女十八| 国产乱人偷精品视频| 精品亚洲乱码少妇综合久久| 久久精品国产亚洲av天美| 乱人伦中国视频| 男人舔女人的私密视频| 黑人巨大精品欧美一区二区蜜桃| 成年av动漫网址| 国产深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 欧美在线黄色| 国产精品香港三级国产av潘金莲 | 日本猛色少妇xxxxx猛交久久| 在线观看三级黄色| 久久影院123| 黑人欧美特级aaaaaa片| 成人漫画全彩无遮挡| 国产成人av激情在线播放| 久久午夜福利片| 一本久久精品| 美国免费a级毛片| 最近手机中文字幕大全| 青春草国产在线视频| 自线自在国产av| 黑人巨大精品欧美一区二区蜜桃| 久久久国产欧美日韩av| 高清欧美精品videossex| 亚洲成人手机| 一本大道久久a久久精品| 在线 av 中文字幕| 久久人人爽人人片av| 精品少妇一区二区三区视频日本电影 | 寂寞人妻少妇视频99o| 国产免费福利视频在线观看| 91午夜精品亚洲一区二区三区| 国产 一区精品| 国产欧美日韩综合在线一区二区| 桃花免费在线播放| 伦理电影免费视频| 天美传媒精品一区二区| 亚洲国产欧美网| 欧美人与性动交α欧美软件| 一边摸一边做爽爽视频免费| 婷婷成人精品国产| 久久精品国产亚洲av天美| 日韩制服丝袜自拍偷拍| 视频区图区小说| 男女边吃奶边做爰视频| 免费观看性生交大片5| 中文乱码字字幕精品一区二区三区| 不卡av一区二区三区| 在线免费观看不下载黄p国产| 自线自在国产av| 天天影视国产精品| 国产一区有黄有色的免费视频| 伊人久久国产一区二区| 成人二区视频| 日韩制服丝袜自拍偷拍| 伊人久久国产一区二区| 国产亚洲精品第一综合不卡| 日韩成人av中文字幕在线观看| 亚洲av福利一区| 高清av免费在线| 国产免费视频播放在线视频| 中文字幕最新亚洲高清| 精品少妇久久久久久888优播| 日日摸夜夜添夜夜爱| 在现免费观看毛片| 久久久国产欧美日韩av| 久久这里有精品视频免费| 精品酒店卫生间| 亚洲精品日韩在线中文字幕| 亚洲精品日韩在线中文字幕| 国产 一区精品| 自拍欧美九色日韩亚洲蝌蚪91| 97在线视频观看| 如日韩欧美国产精品一区二区三区| 国产成人午夜福利电影在线观看| 久久狼人影院| 巨乳人妻的诱惑在线观看| 日韩欧美一区视频在线观看| 蜜桃在线观看..| 亚洲,欧美精品.| 久久av网站| 国产白丝娇喘喷水9色精品| 99精国产麻豆久久婷婷| 熟女电影av网| 男女免费视频国产| 只有这里有精品99| 精品第一国产精品| 99九九在线精品视频| 欧美日韩精品网址| 一级爰片在线观看| 人妻一区二区av| 一区二区三区精品91| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 伦理电影大哥的女人| 美女主播在线视频| 精品少妇内射三级| 午夜影院在线不卡| 电影成人av| 亚洲国产精品999| 亚洲精品日韩在线中文字幕| 黑人猛操日本美女一级片| 在线 av 中文字幕| 久久亚洲国产成人精品v| 久久久亚洲精品成人影院| 免费av中文字幕在线| 三级国产精品片| 人人妻人人爽人人添夜夜欢视频| 久久精品国产自在天天线| 最近手机中文字幕大全| 国产成人精品一,二区| 香蕉丝袜av| √禁漫天堂资源中文www| 国产 精品1| 少妇被粗大的猛进出69影院| 欧美日韩视频精品一区| 在线观看免费高清a一片| 黄色 视频免费看| 在线精品无人区一区二区三| 日韩不卡一区二区三区视频在线| 精品久久久精品久久久| 波多野结衣av一区二区av| 日韩av在线免费看完整版不卡| 欧美中文综合在线视频| 国产 精品1| 午夜福利一区二区在线看| 一级片免费观看大全| 2018国产大陆天天弄谢| 一边摸一边做爽爽视频免费| 人妻一区二区av| 国产 一区精品| av国产精品久久久久影院| 在线亚洲精品国产二区图片欧美| 久久国产精品大桥未久av| 精品一区二区三区四区五区乱码 | 国产麻豆69| 极品人妻少妇av视频| 有码 亚洲区| 国产精品一区二区在线不卡| 啦啦啦视频在线资源免费观看| 欧美人与性动交α欧美软件| 在线观看免费高清a一片| 不卡av一区二区三区| 美女视频免费永久观看网站| 国产又爽黄色视频| 不卡视频在线观看欧美| 叶爱在线成人免费视频播放| 午夜福利网站1000一区二区三区| 99香蕉大伊视频| 国产精品三级大全| 男男h啪啪无遮挡| 国产片特级美女逼逼视频| 婷婷成人精品国产| 尾随美女入室| 免费观看av网站的网址| 欧美成人精品欧美一级黄| 一个人免费看片子| 欧美日本中文国产一区发布| 午夜日本视频在线| 欧美bdsm另类| 中文精品一卡2卡3卡4更新| 精品一区在线观看国产| 七月丁香在线播放| 最新的欧美精品一区二区| 国产精品国产三级专区第一集| 国产极品粉嫩免费观看在线| 七月丁香在线播放| 激情五月婷婷亚洲| 午夜免费鲁丝| 2021少妇久久久久久久久久久| 日韩中文字幕欧美一区二区 | 国产精品二区激情视频| 免费观看在线日韩| 一区在线观看完整版| 中文精品一卡2卡3卡4更新| 天堂8中文在线网| 我要看黄色一级片免费的| 女的被弄到高潮叫床怎么办| 久久这里有精品视频免费| 国产有黄有色有爽视频| 亚洲美女黄色视频免费看| 国产成人精品婷婷| 亚洲成色77777| 欧美xxⅹ黑人| 久久久久视频综合| 大片免费播放器 马上看| 最近2019中文字幕mv第一页| 如何舔出高潮| 欧美在线黄色| 国产成人精品在线电影| 亚洲人成电影观看| 精品午夜福利在线看| 校园人妻丝袜中文字幕| 黄色视频在线播放观看不卡| 极品少妇高潮喷水抽搐| 人人妻人人爽人人添夜夜欢视频| 男女免费视频国产| 亚洲精品国产色婷婷电影| 黄色 视频免费看| 人人妻人人澡人人爽人人夜夜| 美女国产视频在线观看| 日韩 亚洲 欧美在线| 亚洲精品美女久久av网站| 国产精品秋霞免费鲁丝片| 五月天丁香电影| 久久精品aⅴ一区二区三区四区 | 欧美日韩亚洲国产一区二区在线观看 | 男男h啪啪无遮挡| 亚洲第一av免费看| 亚洲av中文av极速乱| 国产免费视频播放在线视频| 午夜福利乱码中文字幕| 最黄视频免费看| 亚洲精品国产av成人精品| 亚洲人成网站在线观看播放| 亚洲欧美成人综合另类久久久| 我要看黄色一级片免费的| 国产精品麻豆人妻色哟哟久久| 伦精品一区二区三区| 国产精品秋霞免费鲁丝片| 在现免费观看毛片| www.av在线官网国产| 十分钟在线观看高清视频www| 国产精品99久久99久久久不卡 | kizo精华| 男人操女人黄网站| 一本久久精品| 国产精品国产av在线观看| 哪个播放器可以免费观看大片| 高清欧美精品videossex| 亚洲视频免费观看视频| 日韩大片免费观看网站| 婷婷色av中文字幕| 久久99一区二区三区| 中国三级夫妇交换| 国产成人aa在线观看| 久久久久国产网址| 一级毛片电影观看| 老熟女久久久| 天天躁夜夜躁狠狠久久av| 久久青草综合色| 日本欧美视频一区| 欧美bdsm另类| 国产黄色视频一区二区在线观看| 岛国毛片在线播放| 男女边吃奶边做爰视频| 久久人人97超碰香蕉20202| 寂寞人妻少妇视频99o| 久久人人爽av亚洲精品天堂| 十八禁高潮呻吟视频| 亚洲av国产av综合av卡| 色播在线永久视频| 国产成人一区二区在线| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 黄色配什么色好看| 国产精品免费大片| 亚洲国产精品国产精品| 91国产中文字幕| 精品国产国语对白av| 少妇人妻精品综合一区二区| 丝袜喷水一区| 国产午夜精品一二区理论片| 在线天堂中文资源库| 欧美日韩成人在线一区二区| 亚洲国产欧美在线一区| 丝袜美足系列| tube8黄色片| 国产男女超爽视频在线观看| 少妇人妻久久综合中文| 卡戴珊不雅视频在线播放| 免费看av在线观看网站| 精品人妻在线不人妻| 久久99精品国语久久久| 免费观看性生交大片5| 97在线人人人人妻| 男女免费视频国产| 97在线视频观看| 国产黄色视频一区二区在线观看| 少妇人妻精品综合一区二区| 九九爱精品视频在线观看| 老汉色av国产亚洲站长工具| 91aial.com中文字幕在线观看| 青青草视频在线视频观看| 亚洲三区欧美一区| 国产av一区二区精品久久| 中文天堂在线官网| 2018国产大陆天天弄谢| 97人妻天天添夜夜摸| 午夜免费观看性视频| 一级,二级,三级黄色视频| 国产精品三级大全| 亚洲婷婷狠狠爱综合网| 丝袜在线中文字幕| 亚洲国产日韩一区二区| 男人操女人黄网站| 亚洲欧洲精品一区二区精品久久久 | 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 久久精品国产亚洲av高清一级| av免费在线看不卡| av天堂久久9| 久久午夜福利片| 妹子高潮喷水视频| 欧美激情极品国产一区二区三区| 成人黄色视频免费在线看| 久久人人爽av亚洲精品天堂| 老汉色∧v一级毛片| 免费观看在线日韩| 成人国产av品久久久| 丰满乱子伦码专区| 最近2019中文字幕mv第一页| 久热这里只有精品99| 这个男人来自地球电影免费观看 | 一边亲一边摸免费视频| 不卡av一区二区三区| 久久精品国产a三级三级三级| 亚洲精品美女久久av网站| 1024视频免费在线观看| 欧美日韩精品成人综合77777| 好男人视频免费观看在线| av在线老鸭窝| 欧美xxⅹ黑人| 高清av免费在线| 免费黄色在线免费观看| 18禁观看日本| 久久久久国产网址| freevideosex欧美| 9色porny在线观看| 欧美成人午夜精品| 亚洲久久久国产精品| 99热国产这里只有精品6| 欧美日韩国产mv在线观看视频| 日本91视频免费播放| 国产福利在线免费观看视频| 一本色道久久久久久精品综合| a级毛片黄视频| 性色av一级| 18+在线观看网站| 欧美人与性动交α欧美精品济南到 | 国产精品一国产av| av在线播放精品| 韩国精品一区二区三区| 精品国产露脸久久av麻豆| av网站在线播放免费| 一边摸一边做爽爽视频免费| 精品人妻熟女毛片av久久网站| 国产国语露脸激情在线看| 人妻人人澡人人爽人人| 亚洲精品久久成人aⅴ小说| 丝袜人妻中文字幕| 91在线精品国自产拍蜜月| 精品亚洲乱码少妇综合久久| 亚洲国产精品一区二区三区在线| 亚洲欧洲精品一区二区精品久久久 | 久久精品国产亚洲av高清一级| 夫妻午夜视频| 欧美亚洲日本最大视频资源| 免费高清在线观看日韩| 国产精品免费大片| 美女视频免费永久观看网站| 99精国产麻豆久久婷婷| 男女边吃奶边做爰视频| 国产精品久久久久久精品古装| 亚洲成人av在线免费| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| videosex国产| 大片电影免费在线观看免费| 丰满饥渴人妻一区二区三| 伦精品一区二区三区| 美女国产视频在线观看| 欧美日韩av久久| 欧美日本中文国产一区发布| 桃花免费在线播放| 一二三四在线观看免费中文在| 亚洲在久久综合| 国产精品秋霞免费鲁丝片| av女优亚洲男人天堂| 免费看不卡的av| 国产欧美日韩综合在线一区二区| 国产精品久久久久成人av| 国产又色又爽无遮挡免| 国产欧美日韩综合在线一区二区| 超碰成人久久| 国产又色又爽无遮挡免| 人人妻人人澡人人看| 国产老妇伦熟女老妇高清| 欧美+日韩+精品| 日日撸夜夜添| 精品少妇一区二区三区视频日本电影 | 亚洲人成网站在线观看播放| 免费在线观看黄色视频的| av不卡在线播放| 欧美黄色片欧美黄色片| 国产男人的电影天堂91| 精品国产超薄肉色丝袜足j| 精品国产一区二区三区久久久樱花| 午夜久久久在线观看| 国产亚洲一区二区精品| 国产精品av久久久久免费| 夫妻午夜视频| 久久人人97超碰香蕉20202| 男女下面插进去视频免费观看| 999精品在线视频| 国产av码专区亚洲av| 下体分泌物呈黄色| 1024香蕉在线观看| 亚洲国产av影院在线观看| 欧美精品亚洲一区二区| 色94色欧美一区二区| 欧美国产精品一级二级三级| 天天躁夜夜躁狠狠躁躁| 母亲3免费完整高清在线观看 | 建设人人有责人人尽责人人享有的| 日韩熟女老妇一区二区性免费视频| 激情五月婷婷亚洲| 少妇的逼水好多| 国产成人免费观看mmmm| 精品99又大又爽又粗少妇毛片| 国产精品一国产av| 国产精品一二三区在线看| 老鸭窝网址在线观看| 免费日韩欧美在线观看| 国产精品av久久久久免费| 卡戴珊不雅视频在线播放| 亚洲精品美女久久久久99蜜臀 | 叶爱在线成人免费视频播放| 免费观看无遮挡的男女| 深夜精品福利| 国产人伦9x9x在线观看 | 国产女主播在线喷水免费视频网站| 日本爱情动作片www.在线观看| 午夜福利乱码中文字幕| 久久久精品国产亚洲av高清涩受| 久久久久国产精品人妻一区二区| 国产白丝娇喘喷水9色精品| 丰满乱子伦码专区| 老女人水多毛片| 欧美日韩成人在线一区二区| 最近手机中文字幕大全| 99国产精品免费福利视频| 亚洲av福利一区| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91 | 九草在线视频观看| 一级毛片 在线播放| 丝袜脚勾引网站| 人妻一区二区av| 国产精品人妻久久久影院| 亚洲精品成人av观看孕妇| 国产免费现黄频在线看| 国产免费视频播放在线视频| 美女中出高潮动态图| 免费播放大片免费观看视频在线观看| 免费黄色在线免费观看| 老司机影院毛片| 国产综合精华液| 美女中出高潮动态图| 免费在线观看黄色视频的| 亚洲国产精品成人久久小说| 国产成人精品福利久久| 亚洲精品国产色婷婷电影| 久久久a久久爽久久v久久| 亚洲一区二区三区欧美精品| 大片免费播放器 马上看| 欧美激情 高清一区二区三区| 七月丁香在线播放| 韩国高清视频一区二区三区| 在线看a的网站| 久久99热这里只频精品6学生| 天天影视国产精品| 日韩不卡一区二区三区视频在线| 一级黄片播放器| 午夜福利在线观看免费完整高清在| 黄片播放在线免费| 成年动漫av网址| 少妇人妻久久综合中文| 亚洲精品第二区| 久久人妻熟女aⅴ| 丰满迷人的少妇在线观看| 欧美在线黄色| 日韩不卡一区二区三区视频在线| 国产黄色视频一区二区在线观看| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看 | 日韩中文字幕欧美一区二区 | 18禁裸乳无遮挡动漫免费视频| 老汉色∧v一级毛片| 日韩av不卡免费在线播放| 亚洲成人手机| 久久久久精品人妻al黑| www.熟女人妻精品国产| 久久久a久久爽久久v久久| 久久精品亚洲av国产电影网| 亚洲一区中文字幕在线| 美女高潮到喷水免费观看| 黄片播放在线免费| 久久久精品免费免费高清| 可以免费在线观看a视频的电影网站 | 久久99蜜桃精品久久| 夜夜骑夜夜射夜夜干| 亚洲一区中文字幕在线| 蜜桃在线观看..| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美一区二区三区国产| 边亲边吃奶的免费视频| 女性被躁到高潮视频| 亚洲经典国产精华液单| 亚洲精品国产色婷婷电影| 国产免费又黄又爽又色| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 成年人免费黄色播放视频| 国产av码专区亚洲av| 亚洲色图综合在线观看| kizo精华| 国产亚洲午夜精品一区二区久久| 少妇被粗大的猛进出69影院| 国产亚洲午夜精品一区二区久久| 午夜免费观看性视频| 大码成人一级视频| 三上悠亚av全集在线观看| a级毛片在线看网站| 日本欧美国产在线视频| 国产av精品麻豆| 国产成人a∨麻豆精品| 精品少妇黑人巨大在线播放| 亚洲欧洲精品一区二区精品久久久 | 18禁国产床啪视频网站| www.熟女人妻精品国产| h视频一区二区三区| 精品午夜福利在线看| 高清在线视频一区二区三区| 日韩伦理黄色片| 晚上一个人看的免费电影| 香蕉国产在线看| 麻豆精品久久久久久蜜桃| 亚洲av日韩在线播放| 狂野欧美激情性bbbbbb| 精品人妻偷拍中文字幕| a 毛片基地| 永久网站在线| 老司机影院成人| 成人亚洲欧美一区二区av| 亚洲精品av麻豆狂野| 成年人免费黄色播放视频| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 人人妻人人添人人爽欧美一区卜| 欧美黄色片欧美黄色片| 亚洲国产精品一区三区| 亚洲男人天堂网一区| 成人漫画全彩无遮挡| 宅男免费午夜| 国产人伦9x9x在线观看 | 精品国产乱码久久久久久小说| 国产无遮挡羞羞视频在线观看| 欧美日韩成人在线一区二区| 国产精品国产三级国产专区5o| 国产精品国产三级专区第一集| 99热全是精品| 日产精品乱码卡一卡2卡三| 精品少妇一区二区三区视频日本电影 | 99久久中文字幕三级久久日本| 韩国精品一区二区三区| 国产熟女午夜一区二区三区| 日韩大片免费观看网站| 伊人久久大香线蕉亚洲五| 五月伊人婷婷丁香| 天天躁夜夜躁狠狠躁躁| 久久这里只有精品19| 国产野战对白在线观看| 欧美日韩精品成人综合77777| 韩国高清视频一区二区三区| 久久99精品国语久久久|