• <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)化路徑研究
    你懂的网址亚洲精品在线观看| 国产精品国产三级国产av玫瑰| 黄片播放在线免费| 韩国高清视频一区二区三区| 美女中出高潮动态图| 免费av不卡在线播放| 国产极品天堂在线| av一本久久久久| 亚洲精品日本国产第一区| 亚洲精品一区蜜桃| 国产av精品麻豆| 亚洲精品日韩av片在线观看| 日本av手机在线免费观看| 中文乱码字字幕精品一区二区三区| 精品人妻一区二区三区麻豆| 国产有黄有色有爽视频| 99久久中文字幕三级久久日本| 黄色怎么调成土黄色| 日韩一区二区三区影片| 少妇的逼水好多| 亚洲国产精品一区二区三区在线| 久久人人爽人人爽人人片va| 最后的刺客免费高清国语| 18禁观看日本| 亚洲成色77777| 精品国产乱码久久久久久小说| 综合色丁香网| 中文字幕亚洲精品专区| 国产白丝娇喘喷水9色精品| 国产熟女午夜一区二区三区 | 简卡轻食公司| 久久久久久久久久人人人人人人| 国产成人精品在线电影| av.在线天堂| 99国产精品免费福利视频| 99热6这里只有精品| 国产亚洲一区二区精品| 国产日韩欧美视频二区| 亚洲精品456在线播放app| 91精品国产九色| 成人无遮挡网站| 日本欧美国产在线视频| 男人爽女人下面视频在线观看| 大香蕉97超碰在线| 国产精品免费大片| 精品一区二区免费观看| 久久精品国产亚洲网站| 色婷婷久久久亚洲欧美| 老司机亚洲免费影院| 秋霞在线观看毛片| 久久人人爽人人片av| 国产亚洲精品久久久com| 简卡轻食公司| a级毛色黄片| 亚洲国产毛片av蜜桃av| 免费观看av网站的网址| 一级片'在线观看视频| 国产亚洲精品久久久com| 又黄又爽又刺激的免费视频.| 美女脱内裤让男人舔精品视频| 亚洲怡红院男人天堂| 啦啦啦中文免费视频观看日本| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 国产成人91sexporn| 亚洲高清免费不卡视频| 丝袜喷水一区| 久久人妻熟女aⅴ| 五月玫瑰六月丁香| 一本一本综合久久| 久久婷婷青草| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区在线不卡| 婷婷色av中文字幕| 日本黄色日本黄色录像| 少妇被粗大猛烈的视频| 欧美日韩国产mv在线观看视频| 亚洲国产精品国产精品| 亚洲国产精品一区二区三区在线| 在线观看国产h片| 亚洲五月色婷婷综合| 另类亚洲欧美激情| 男人操女人黄网站| 日日摸夜夜添夜夜爱| 国产爽快片一区二区三区| 黄色一级大片看看| 国产一区亚洲一区在线观看| 亚洲精品成人av观看孕妇| 蜜桃在线观看..| 亚洲精品,欧美精品| 少妇被粗大猛烈的视频| 性色avwww在线观看| videos熟女内射| 免费观看av网站的网址| 欧美国产精品一级二级三级| 国产精品久久久久久精品古装| 久久精品国产亚洲av涩爱| 久久精品久久久久久噜噜老黄| 美女内射精品一级片tv| 黄色欧美视频在线观看| 免费观看av网站的网址| 制服人妻中文乱码| 亚洲av综合色区一区| 欧美日本中文国产一区发布| 女人精品久久久久毛片| 老司机影院成人| 国产免费又黄又爽又色| 一区二区三区精品91| 视频中文字幕在线观看| 久久久久久久国产电影| 插逼视频在线观看| 国产成人精品婷婷| 欧美3d第一页| 99国产综合亚洲精品| 免费不卡的大黄色大毛片视频在线观看| 我要看黄色一级片免费的| 国产免费又黄又爽又色| 日本av免费视频播放| 国产高清有码在线观看视频| 亚洲欧美色中文字幕在线| 天天躁夜夜躁狠狠久久av| 亚洲欧洲精品一区二区精品久久久 | 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 乱人伦中国视频| 精品国产露脸久久av麻豆| 亚洲国产欧美在线一区| 国产亚洲一区二区精品| 国产精品久久久久久久电影| 日韩强制内射视频| 亚洲图色成人| 欧美日韩一区二区视频在线观看视频在线| 亚洲av电影在线观看一区二区三区| 妹子高潮喷水视频| 国产高清不卡午夜福利| 精品午夜福利在线看| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 街头女战士在线观看网站| av国产精品久久久久影院| 久久精品国产鲁丝片午夜精品| 一级毛片黄色毛片免费观看视频| 国产高清国产精品国产三级| av黄色大香蕉| 欧美精品一区二区大全| 国精品久久久久久国模美| 精品人妻一区二区三区麻豆| 欧美日韩成人在线一区二区| 午夜福利视频精品| 午夜免费观看性视频| 成人手机av| 好男人视频免费观看在线| 青春草视频在线免费观看| 婷婷色麻豆天堂久久| 日本av手机在线免费观看| 亚洲少妇的诱惑av| 蜜臀久久99精品久久宅男| 一区二区三区四区激情视频| 在线观看一区二区三区激情| 男人操女人黄网站| 亚洲国产精品成人久久小说| 少妇高潮的动态图| 多毛熟女@视频| 亚洲精品456在线播放app| 最近中文字幕2019免费版| 亚洲怡红院男人天堂| 久久97久久精品| 日本黄色片子视频| av专区在线播放| 制服丝袜香蕉在线| 日韩中字成人| 亚洲欧洲日产国产| 九色成人免费人妻av| 国产成人a∨麻豆精品| 国产精品久久久久久精品古装| 日韩欧美一区视频在线观看| 国产欧美亚洲国产| 国产成人精品福利久久| 伊人亚洲综合成人网| 久久久亚洲精品成人影院| 日韩伦理黄色片| 日韩av免费高清视频| 青春草亚洲视频在线观看| 亚洲av免费高清在线观看| 乱人伦中国视频| 欧美老熟妇乱子伦牲交| 在线天堂最新版资源| 26uuu在线亚洲综合色| xxx大片免费视频| 有码 亚洲区| 国产乱来视频区| 女人久久www免费人成看片| 久久久国产一区二区| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 久久久久人妻精品一区果冻| 成人综合一区亚洲| 亚洲欧美成人综合另类久久久| 国产片特级美女逼逼视频| 精品少妇黑人巨大在线播放| 亚洲av福利一区| 丝袜在线中文字幕| 午夜福利在线观看免费完整高清在| 搡老乐熟女国产| 久久精品国产亚洲av天美| 纵有疾风起免费观看全集完整版| 国模一区二区三区四区视频| 91精品三级在线观看| 精品久久蜜臀av无| av线在线观看网站| 91国产中文字幕| 99久久精品国产国产毛片| 久久久久久久久久成人| 午夜久久久在线观看| 国产欧美日韩一区二区三区在线 | 免费黄色在线免费观看| 国产有黄有色有爽视频| 中文字幕av电影在线播放| 国产淫语在线视频| 久久久久久伊人网av| 国产乱来视频区| 人体艺术视频欧美日本| 在线观看国产h片| 日本91视频免费播放| 免费高清在线观看视频在线观看| 插阴视频在线观看视频| 99国产精品免费福利视频| 亚洲国产欧美在线一区| av在线app专区| 欧美日韩成人在线一区二区| 肉色欧美久久久久久久蜜桃| 综合色丁香网| 免费看不卡的av| 中国三级夫妇交换| 久久精品国产亚洲av天美| 亚洲精品久久久久久婷婷小说| 免费观看无遮挡的男女| 尾随美女入室| 国产精品欧美亚洲77777| 久久精品国产a三级三级三级| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 久久av网站| 美女国产视频在线观看| 99热这里只有精品一区| 久久99蜜桃精品久久| 老司机影院成人| 性高湖久久久久久久久免费观看| 激情五月婷婷亚洲| 欧美97在线视频| 视频区图区小说| 97精品久久久久久久久久精品| 免费高清在线观看日韩| 91精品国产国语对白视频| 又黄又爽又刺激的免费视频.| 免费观看a级毛片全部| 亚洲av免费高清在线观看| 日韩欧美一区视频在线观看| 亚洲一区二区三区欧美精品| 久久午夜综合久久蜜桃| 成年av动漫网址| 久久国产精品大桥未久av| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区| 午夜免费观看性视频| 日韩一本色道免费dvd| 成人国产av品久久久| 亚洲不卡免费看| 一级爰片在线观看| a级毛片在线看网站| a级毛片在线看网站| 91久久精品国产一区二区三区| 国模一区二区三区四区视频| 少妇 在线观看| 国产精品国产av在线观看| 天天操日日干夜夜撸| 久久久久精品久久久久真实原创| 亚洲一级一片aⅴ在线观看| 久久久国产一区二区| 亚洲av免费高清在线观看| 99热这里只有精品一区| 精品亚洲乱码少妇综合久久| 波野结衣二区三区在线| 有码 亚洲区| 美女视频免费永久观看网站| 精品酒店卫生间| 亚州av有码| 欧美成人午夜免费资源| 国产极品粉嫩免费观看在线 | 亚洲av日韩在线播放| 18禁在线无遮挡免费观看视频| 人人妻人人爽人人添夜夜欢视频| 成人影院久久| 国产av一区二区精品久久| 99热网站在线观看| 亚洲人成网站在线播| 91aial.com中文字幕在线观看| 国产精品嫩草影院av在线观看| 久久久久久久久久久免费av| 国产精品免费大片| 日本欧美视频一区| 久久久久久伊人网av| 亚洲精品国产av蜜桃| 2022亚洲国产成人精品| 国产精品99久久99久久久不卡 | 久久久久国产精品人妻一区二区| 日本黄大片高清| 国产av国产精品国产| 国产成人午夜福利电影在线观看| 一级毛片电影观看| 精品视频人人做人人爽| 国产男女内射视频| 亚洲三级黄色毛片| 精品午夜福利在线看| 精品人妻一区二区三区麻豆| 久久久久久久久大av| 熟妇人妻不卡中文字幕| 亚洲美女视频黄频| 九九爱精品视频在线观看| 亚洲美女黄色视频免费看| 国产精品一区www在线观看| 夜夜看夜夜爽夜夜摸| 99久久综合免费| 免费看光身美女| 中文字幕最新亚洲高清| 热99久久久久精品小说推荐| 夫妻午夜视频| 另类亚洲欧美激情| 丰满迷人的少妇在线观看| 国产精品久久久人人做人人爽| 大片电影免费在线观看免费| 国产xxxxx性猛交| 国产成人免费观看mmmm| 免费黄频网站在线观看国产| 老司机午夜福利在线观看视频 | 国产极品粉嫩免费观看在线| 日韩制服丝袜自拍偷拍| 亚洲国产av影院在线观看| 亚洲专区中文字幕在线| 久久久精品国产亚洲av高清涩受| 亚洲全国av大片| av福利片在线| 丝瓜视频免费看黄片| 国产精品美女特级片免费视频播放器 | 啪啪无遮挡十八禁网站| av免费在线观看网站| 精品福利永久在线观看| 在线观看66精品国产| 久久亚洲真实| 亚洲五月色婷婷综合| 最近最新免费中文字幕在线| 一边摸一边做爽爽视频免费| 最新的欧美精品一区二区| 精品福利观看| 亚洲精华国产精华精| 黄色视频在线播放观看不卡| 一级片免费观看大全| 女性被躁到高潮视频| 亚洲av日韩在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 每晚都被弄得嗷嗷叫到高潮| 久久午夜亚洲精品久久| 精品国产亚洲在线| 波多野结衣一区麻豆| 久久婷婷成人综合色麻豆| 大香蕉久久成人网| 国产一区二区三区视频了| 日本av免费视频播放| 热re99久久精品国产66热6| 久久性视频一级片| 老司机靠b影院| 成人黄色视频免费在线看| 999精品在线视频| 欧美成狂野欧美在线观看| 超碰97精品在线观看| 免费人妻精品一区二区三区视频| 美女视频免费永久观看网站| 视频区图区小说| 99久久人妻综合| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 少妇的丰满在线观看| 在线观看一区二区三区激情| 国产亚洲午夜精品一区二区久久| 考比视频在线观看| 久久人人爽av亚洲精品天堂| 黑人操中国人逼视频| 老司机影院毛片| 搡老熟女国产l中国老女人| 日韩视频一区二区在线观看| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频| 国产精品99久久99久久久不卡| 1024视频免费在线观看| 日本黄色视频三级网站网址 | av国产精品久久久久影院| 19禁男女啪啪无遮挡网站| 国产精品电影一区二区三区 | 性少妇av在线| 波多野结衣av一区二区av| 99热网站在线观看| 午夜福利欧美成人| av天堂在线播放| 国产区一区二久久| 男女午夜视频在线观看| 午夜老司机福利片| 丁香欧美五月| 多毛熟女@视频| 国产91精品成人一区二区三区 | 五月开心婷婷网| 99久久99久久久精品蜜桃| 俄罗斯特黄特色一大片| 一边摸一边抽搐一进一出视频| 交换朋友夫妻互换小说| 亚洲熟女毛片儿| 波多野结衣一区麻豆| 999久久久国产精品视频| 国产国语露脸激情在线看| av线在线观看网站| 精品午夜福利视频在线观看一区 | 亚洲成国产人片在线观看| 亚洲国产看品久久| 亚洲欧洲精品一区二区精品久久久| 人人妻人人爽人人添夜夜欢视频| 美女福利国产在线| 国产亚洲欧美在线一区二区| 亚洲情色 制服丝袜| tube8黄色片| a在线观看视频网站| 在线观看66精品国产| 久久亚洲真实| 欧美乱码精品一区二区三区| 十八禁网站免费在线| 最新美女视频免费是黄的| 精品人妻在线不人妻| 国产区一区二久久| 天堂俺去俺来也www色官网| 香蕉丝袜av| 成人18禁在线播放| 丰满饥渴人妻一区二区三| 肉色欧美久久久久久久蜜桃| 久久 成人 亚洲| 999精品在线视频| 久久免费观看电影| 我的亚洲天堂| 人妻 亚洲 视频| 久久久精品国产亚洲av高清涩受| 亚洲中文字幕日韩| 人成视频在线观看免费观看| 久久精品国产综合久久久| 久久亚洲精品不卡| 久久人妻福利社区极品人妻图片| 午夜激情av网站| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 首页视频小说图片口味搜索| 一区福利在线观看| 青青草视频在线视频观看| 久久精品人人爽人人爽视色| 五月天丁香电影| 高清毛片免费观看视频网站 | 欧美大码av| 国产一区二区激情短视频| 一本大道久久a久久精品| 夜夜爽天天搞| 国产精品久久久久久精品古装| 99热国产这里只有精品6| 涩涩av久久男人的天堂| 男女高潮啪啪啪动态图| 五月开心婷婷网| 99国产精品一区二区三区| 亚洲av国产av综合av卡| 好男人电影高清在线观看| 精品久久蜜臀av无| 久久av网站| 色综合婷婷激情| 大型av网站在线播放| 亚洲成人免费电影在线观看| 亚洲人成伊人成综合网2020| 精品人妻在线不人妻| 精品国内亚洲2022精品成人 | 亚洲精品自拍成人| 久久久久久久大尺度免费视频| 首页视频小说图片口味搜索| 亚洲精品一二三| 成人精品一区二区免费| 亚洲av欧美aⅴ国产| 丰满饥渴人妻一区二区三| 国产国语露脸激情在线看| 久久国产精品男人的天堂亚洲| 久久国产精品人妻蜜桃| 黄色视频,在线免费观看| 水蜜桃什么品种好| 精品国产乱码久久久久久小说| 国产精品久久久av美女十八| 制服人妻中文乱码| 国产精品自产拍在线观看55亚洲 | 国产精品国产av在线观看| 啦啦啦视频在线资源免费观看| 日本a在线网址| 激情视频va一区二区三区| 亚洲国产av影院在线观看| av欧美777| 亚洲精品久久午夜乱码| 国产av精品麻豆| 久久精品亚洲熟妇少妇任你| 欧美乱妇无乱码| 色在线成人网| 国产在线精品亚洲第一网站| 亚洲精品粉嫩美女一区| 欧美国产精品va在线观看不卡| 国产精品成人在线| 丝袜喷水一区| 91老司机精品| 91av网站免费观看| 动漫黄色视频在线观看| 国产97色在线日韩免费| a级片在线免费高清观看视频| 亚洲成人手机| 免费看十八禁软件| 日本五十路高清| 欧美精品av麻豆av| 大陆偷拍与自拍| 国产精品亚洲一级av第二区| 操美女的视频在线观看| 日韩欧美三级三区| 十八禁高潮呻吟视频| 老司机福利观看| 亚洲精品美女久久av网站| 高清欧美精品videossex| 在线观看免费日韩欧美大片| 波多野结衣av一区二区av| 亚洲情色 制服丝袜| 日韩欧美免费精品| 久久久久久免费高清国产稀缺| 欧美黄色淫秽网站| 亚洲成人国产一区在线观看| av电影中文网址| 青草久久国产| 一个人免费看片子| 黑人欧美特级aaaaaa片| 黄频高清免费视频| 国产一区二区三区视频了| 精品久久蜜臀av无| 亚洲熟妇熟女久久| 精品福利观看| 十分钟在线观看高清视频www| 亚洲,欧美精品.| 一本大道久久a久久精品| 国产福利在线免费观看视频| 精品国产亚洲在线| av又黄又爽大尺度在线免费看| 十八禁网站免费在线| www日本在线高清视频| 欧美日韩视频精品一区| av在线播放免费不卡| 精品少妇久久久久久888优播| 满18在线观看网站| 亚洲专区中文字幕在线| 伦理电影免费视频| 精品国产乱码久久久久久男人| 亚洲三区欧美一区| 精品欧美一区二区三区在线| 国产在线免费精品| 麻豆成人av在线观看| 亚洲免费av在线视频| 香蕉丝袜av| 999久久久精品免费观看国产| 18禁裸乳无遮挡动漫免费视频| www.999成人在线观看| 一边摸一边抽搐一进一出视频| 国产国语露脸激情在线看| 国产欧美日韩一区二区三| 久久精品亚洲精品国产色婷小说| 欧美成狂野欧美在线观看| 久久狼人影院| 亚洲av片天天在线观看| 久久精品亚洲av国产电影网| 欧美在线黄色| 热re99久久国产66热| 亚洲人成电影观看| 天天影视国产精品| 国产1区2区3区精品| 这个男人来自地球电影免费观看| 亚洲av成人不卡在线观看播放网| 国产aⅴ精品一区二区三区波| 久久人妻熟女aⅴ| 老熟妇仑乱视频hdxx| av一本久久久久| 少妇裸体淫交视频免费看高清 | 久久 成人 亚洲| 亚洲av国产av综合av卡| 国产成人系列免费观看| 日日摸夜夜添夜夜添小说| 18在线观看网站| 新久久久久国产一级毛片| 淫妇啪啪啪对白视频| 夜夜爽天天搞| videos熟女内射| 国产99久久九九免费精品| 1024视频免费在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久午夜乱码| 久久精品91无色码中文字幕| 一个人免费看片子| 丁香六月欧美| 亚洲精华国产精华精| 日韩成人在线观看一区二区三区| 日日爽夜夜爽网站| 视频在线观看一区二区三区| 国产成人系列免费观看| 美女福利国产在线| 青草久久国产| tube8黄色片| 成人精品一区二区免费|