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

    秸稈還田與施肥對(duì)稻田土壤微生物生物量及固氮菌群落結(jié)構(gòu)的影響

    2013-12-08 12:45:27劉驍蒨涂仕華孫錫發(fā)辜運(yùn)富張先琴張小平
    生態(tài)學(xué)報(bào) 2013年17期
    關(guān)鍵詞:固氮菌固氮條帶

    劉驍蒨,涂仕華,孫錫發(fā),辜運(yùn)富,張先琴,張小平,*

    (1. 四川農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院微生物系, 成都 611130; 2. 四川省農(nóng)業(yè)科學(xué)院土壤肥料研究所,成都 610066)

    秸稈還田與施肥對(duì)稻田土壤微生物生物量及固氮菌群落結(jié)構(gòu)的影響

    劉驍蒨1,涂仕華2,孫錫發(fā)2,辜運(yùn)富1,張先琴1,張小平1,*

    (1. 四川農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院微生物系, 成都 611130; 2. 四川省農(nóng)業(yè)科學(xué)院土壤肥料研究所,成都 610066)

    利用氯仿熏蒸法和變性梯度凝膠電泳法(PCR-DGGE)研究了秸稈覆蓋還田與施肥對(duì)灰棕沖積水稻土0—10cm和10—20cm土層土壤微生物生物量碳、氮和固氮菌群落結(jié)構(gòu)的影響。結(jié)果表明:土壤微生物量碳、氮和固氮菌多樣性從0—10cm土層到10—20cm土層均呈現(xiàn)降低趨勢(shì)。無秸稈覆蓋處理(對(duì)照組)的土壤微生物生物量碳(SMB-C)和微生物生物量氮(SMB-N)量最小。在秸稈覆蓋還田處理中,低氮和無鉀處理的SMB-C和SMB-N都顯著低于全量氮磷鉀肥處理。雖然無磷處理的SMB-N低于全量氮磷鉀處理, 但差異不顯著。說明秸稈覆蓋還田配施充足氮磷鉀肥能顯著提高土壤微生物生物量碳、氮。由DGGE圖譜多樣性指數(shù)分析得知,配施充足氮磷鉀肥的處理土壤的固氮菌多樣性最豐富。UPGMA聚類分析顯示,10種不同處理的聚類圖也不同,對(duì)照(無秸稈)處理0—10cm和10—20cm的微生物不同于其它處理單獨(dú)聚在了一個(gè)群里。DGGE條帶測(cè)序得知,14個(gè)條帶的近緣種大部分為非培養(yǎng)細(xì)菌nifH基因片段,主要優(yōu)勢(shì)菌群其歸屬于變形菌門(Proteobacteria)的β-變形菌綱(Betaproteobacteria)。應(yīng)用PCR-DGGE技術(shù)可以解釋灰棕沖積水稻土秸稈覆蓋不同肥料用量固氮菌分子群落結(jié)構(gòu)特點(diǎn)。

    秸稈還田;灰棕沖積水稻土;微生物生物量;PCR-DGGE;固氮菌群落結(jié)構(gòu)

    我國(guó)的秸稈資源十分豐富,據(jù)統(tǒng)計(jì)全國(guó)每年共產(chǎn)生秸稈數(shù)量為6.2億t,約占世界的20%—30%[1],農(nóng)作物秸稈中蘊(yùn)含有大量的氮、磷、鉀和中、微量元素等礦質(zhì)養(yǎng)分。農(nóng)作生產(chǎn)中產(chǎn)生的大量秸稈被焚燒,造成了營(yíng)養(yǎng)物質(zhì)流失和嚴(yán)重環(huán)境污染。秸稈還田作為秸稈利用的一種重要且有效的方式,既可以避免資源浪費(fèi)和環(huán)境污染,同時(shí)也可以增加土壤有機(jī)質(zhì)和養(yǎng)分的含量,提高孔隙度,降低土壤容重,還可以提高土壤微生物生物量及增強(qiáng)各種土壤酶的活性[2],從而使作物更好更快地生長(zhǎng),提高作物產(chǎn)量。

    土壤微生物是土壤生態(tài)系統(tǒng)的重要組成部分,在有機(jī)質(zhì)分解、土壤養(yǎng)分循環(huán)和植物養(yǎng)分吸收利用過程中發(fā)揮著至關(guān)重要的作用[3]。生物多樣性越高的生態(tài)系統(tǒng)越穩(wěn)定的Elton假說,受到大多數(shù)學(xué)者的認(rèn)同。Tilman通過直接控制生物多樣性,對(duì)147個(gè)草地進(jìn)行了實(shí)驗(yàn)研究。結(jié)果表明,草地生態(tài)系統(tǒng)的生物多樣性影響著草地短期和長(zhǎng)期的生產(chǎn)率及發(fā)展[4]。一些研究者指出隨著有機(jī)質(zhì)的施入和化學(xué)肥料的減少,土壤微生物生物量和微生物活性都明顯增加[5]。秸稈覆蓋與土壤微生物之間作用是相互促進(jìn)相互影響的,一方面秸稈覆蓋增加了土壤微生物生物量,提高微生物的活性[6],另一方面土壤微生物生物量及活性的提高也會(huì)加快秸稈的腐解[7]。

    傳統(tǒng)平板培養(yǎng)分離方法可分離出的微生物種類只占土壤微生物種類總數(shù)的0.1%—1%[8]。應(yīng)用PCR-DGGE技術(shù)可以有效地避免在傳統(tǒng)菌落培養(yǎng)分離過程由于微生態(tài)生境無法復(fù)制而造成的大量微生物多樣性的丟失,能夠更加全面而可靠的反映土壤微生物的組成情況,故其在微生物群落多樣性研究中得到越來越廣泛的應(yīng)用[9]。固氮細(xì)菌是生態(tài)系統(tǒng)中一類重要的功能菌群,它在土壤氮素循環(huán)發(fā)揮著不可取代的作用[10]。目前有關(guān)秸稈還田條件下配合化肥使用對(duì)土壤微生物量和特殊生理類群微生物群落結(jié)果影響的研究還不多,關(guān)于秸稈還田下施用不同肥料對(duì)四川灰棕沖積水稻土微生物量及固氮菌群落結(jié)構(gòu)的影響的研究更無報(bào)道。本研究將應(yīng)用PCR-DGGE技術(shù)研究秸稈還田下不同施肥土壤固氮菌群落結(jié)構(gòu),并對(duì)土壤中的主要固氮菌群進(jìn)行測(cè)序分析,從微生物角度了解秸稈還田土壤生態(tài)狀況,為秸稈還田與合理施肥、水分管理、保護(hù)土壤微生物等方面的工作提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)點(diǎn)及實(shí)驗(yàn)處理

    試驗(yàn)土壤采集于四川省農(nóng)業(yè)科學(xué)院土壤肥料研究所于2002年建立在四川省廣漢市西高鎮(zhèn)的試驗(yàn)基地。該地為亞熱帶濕潤(rùn)氣候,年均氣溫16.3℃,積溫5980.4℃,降水量為890.8mm,日照時(shí)數(shù)為1229.2h,無霜期達(dá)281d。試驗(yàn)開始于2004年10月,作物輪作為小麥水稻,一年兩熟。土壤為灰棕沖積水稻土,試驗(yàn)前的土壤基本理化性質(zhì)見表1。試驗(yàn)采取完全隨機(jī)區(qū)組設(shè)計(jì)方案,5個(gè)處理:N150P60K90(CK,無覆蓋),秸稈覆蓋+N90P60K90(SNL),秸稈覆蓋+N150P60K90(SCK),秸稈覆蓋+N150P60K0(SNP)和秸稈覆蓋+N150P0K90(SNK),3次重復(fù),小區(qū)面積為3m×4m。試驗(yàn)小區(qū)四周為20cm高田埂,重復(fù)之間有排水溝。各處理肥料用量分別為,氮肥(N,尿素) 90、150 kg/hm2;磷肥(P2O5,磷酸一銨) 0、60 kg/ hm2;鉀肥(K2O,氯化鉀) 0、90 kg/hm2。每季作物收獲后,秸稈全部均勻覆蓋于原試驗(yàn)小區(qū),連續(xù)覆蓋免耕。小麥和水稻施肥時(shí)間和比例分別為,氮肥30%作底肥施用,70%作分蘗肥施用;磷肥全部作底肥施用;鉀肥50%作底肥,50%作分蘗肥。

    表1 供試土壤基本理化性質(zhì)

    1.2 取樣方法

    于2009年7月水稻生長(zhǎng)旺盛期,在田間采用“梅花型”布點(diǎn)取樣,取樣深度分別為0—10cm(在處理代碼后加T下標(biāo))和10—20cm(在處理代碼后加s下標(biāo))。將每一小區(qū)內(nèi)同一取樣點(diǎn)中相同深度土壤混勻,用無菌PET樹脂袋封裝帶回實(shí)驗(yàn)室。立即進(jìn)行土壤微生物數(shù)量分析和DNA的提取。另取一部分放于室內(nèi)風(fēng)干,研細(xì)過2 mm篩,用于土壤理化性質(zhì)的分析。其它放在-20℃的冰箱里保存以備后續(xù)分析。

    1.3 土壤微生物生物量的測(cè)定。

    土壤微生物生物量土壤微生物生物量碳(MBC)、微生物量氮(MBN)的測(cè)定采用氯仿熏蒸法[11]。

    1.4 土壤固氮菌群落結(jié)構(gòu)的PCR-DGGE分析和系統(tǒng)發(fā)育樹的構(gòu)建

    1.4.1主要儀器和試劑DGGE所用儀器the D-code TM Universal Mutation Detection System (Bio-Rad Co)。采取巢式一步PCR擴(kuò)增nifH基因,巢式一輪引物為FGPH19和PolR,巢式第二輪引物為AQER和PolF-GC。(引物均由Invitrogen英杰生命技術(shù)公司合成),具體引物序列參見文獻(xiàn)[12]。切膠純化試劑盒(TaKaRa,大連寶生物)純化。

    1.4.2土壤微生物總DNA提取 采用土壤總DNA提取試劑盒(北京天恩澤基因科技有限公司),稱取0.5g于20℃保存的土樣,按試劑盒提供的操作步驟進(jìn)行土壤微生物總DNA的提取,在0. 8%的瓊脂糖凝膠中對(duì)提取的DNA進(jìn)行電泳檢測(cè)。

    1.4.3固氮菌nif-H基因的PCR擴(kuò)增采用Nested-PCR程序。

    第一次PCR 所用引物為固氮菌nif-H基因的引物FGPH 19和Po1R,反應(yīng)體系為:PCR Master Mix (Tiangen Biotech, Bejing) 25L,每種引物0.5μL (20pmol/μL),0.5μL 10倍稀釋的土壤總DNA提取液,加ddH2O至反應(yīng)終體積50μL。為減少擴(kuò)增過程中的非特異性產(chǎn)物,采用Touch-down PCR程序,反應(yīng)程序:94℃預(yù)變性5min,94℃ 1min,65—55℃ 50s (每個(gè)循環(huán)溫度降低0.5℃,72℃ 1.30min,20個(gè)循環(huán),然后在其它條件不變的情況下,在55℃的退火溫度下繼續(xù)擴(kuò)增15個(gè)循環(huán),72℃ 7 min,最后于4℃恒定保存。

    第二次PCR 第一次PCR產(chǎn)物1∶10稀釋后作為二次PCR反應(yīng)的模板進(jìn)行,引物為AQER和PolF-GC,反應(yīng)的體系為:10×PCR buffer緩沖液5μL, MgCl2(25mmol/μL) 3μL, dNTP(2.5mmol/μL) 1μL, PolF-GC (20pmol/μL) 0.5μL, AQER(20pmol/μL) 0.5μL, Taq DNA聚合酶(5u/μL) (Tiangen Biotech, Bejing)0.5μL,模板DNA1μL,加ddH2O至終體積50μL。反應(yīng)程序是在第一輪PCR的基礎(chǔ)上,退火溫度變?yōu)樵谇?0個(gè)循環(huán)為63—53℃,接下來的后15個(gè)循環(huán)的退火溫度為53℃,其它步驟與第一輪PCR反應(yīng)相同。

    1.4.4 PCR產(chǎn)物的變性梯度凝膠電泳(DGGE)分析

    取第一次PCR產(chǎn)物15μL進(jìn)行DGGE分析,變性劑梯度范圍為20%—50%,自配的聚丙烯酰胺凝膠濃度為8% (100%的變性劑為尿素7 mol/L和40%的去離子甲酰胺)。在1×TAE緩沖液中,50V電壓進(jìn)膠30 min,然后在150V 60℃緩沖液中電泳6h。電泳后采用銀染法進(jìn)行染色。然后用數(shù)碼相機(jī)拍照。DGGE指紋圖譜分析借助于Bio-Rad公司的凝膠成像系統(tǒng)(Quantity One Bio-Rad, USA)分析樣品電泳條帶。

    在圖像處理過程中,在DGGE電泳圖譜上有一些肉眼可見、但被軟件忽略掉的細(xì)小條帶進(jìn)行手動(dòng)添加,條帶的密度由該軟件自動(dòng)算出。

    1.4.5 DGGE條帶的切割、克隆與測(cè)序

    從DGGE凝膠上仔細(xì)切下有用的DGGE條帶裝于PCR反應(yīng)管中,用50μL無菌水反復(fù)沖洗3次后,加人20μL無菌水在4℃冰箱過夜。次日12000g離心5 min,收集上清液。以此作為模板進(jìn)行下一步的PCR擴(kuò)增,用引物AQER和PolF,進(jìn)行nif-H基因PCR擴(kuò)增。PCR反應(yīng)程序同1.4.3中第一次PCR反應(yīng)。將上述PCR產(chǎn)物連接到pMD19-T載體(Takara, 大連)上,轉(zhuǎn)化到E.coliDH5α感受態(tài)細(xì)胞中,在含有氨芐青霉素的LB培養(yǎng)基上選擇具有氨芐青霉素抗性的白色轉(zhuǎn)化子。采用T載體通用引物M13F, M13R進(jìn)行菌落PCR,PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測(cè)是否為陽(yáng)性克隆。測(cè)序工作由北京華大基因科技有限公司完成。

    1.4.6 固氮基因nif-H序列分析

    所得序列提交到GenBank數(shù)據(jù)庫(kù)中,并用BLAST程序進(jìn)行同源性比對(duì),獲得與本實(shí)驗(yàn)所得序列相似性最高的序列。將這些相似性較高序列與本實(shí)驗(yàn)獲得的固氮基因nif-H序列一起,使用Clustal X1.83工具進(jìn)行多序列比對(duì)。然后使用MEGA軟件中的Kimura雙參數(shù)修正模型進(jìn)行系統(tǒng)發(fā)育樹的構(gòu)建。

    1.5 數(shù)據(jù)處理

    用Quantity One分析軟件分析得到相似性指數(shù);以Shannon-Wiener指數(shù)反映固氮菌多樣性,香濃指數(shù)計(jì)算方法采用Eichner等[13]和Hedrick等[14]的方法;數(shù)據(jù)處理用Microsoft Excel 2000進(jìn)行,統(tǒng)計(jì)與顯著性檢測(cè)利用SPSS11.0軟件進(jìn)行。

    2 結(jié)果與分析

    2.1 秸稈覆蓋與不同施肥量對(duì)土壤微生物生物量碳、氮的影響

    農(nóng)田生態(tài)系統(tǒng)中土壤微生物生物量的變化受到多種因素的影響如環(huán)境條件、施肥管理等,同時(shí)它們能夠靈敏地反映環(huán)境因子的變化,故常被用于土壤質(zhì)量生物學(xué)的評(píng)價(jià)指標(biāo)之一[15]。本研究中由于覆蓋秸稈和不同的施肥處理,導(dǎo)致土壤微生物利用的碳源數(shù)量、方式發(fā)生變化,從而使不同處理的土壤微生物生物量碳、氮表現(xiàn)出差異。供試土壤微生物生物量碳的變化范圍介于156.97—229.37mg/kg,微生物生物量氮的變化則介于21.87—35.85mg/kg,最低值均是沒有覆蓋秸稈的對(duì)照組的下層(10—20cm)土樣。SCKT與CKT處理的微生物生物量碳、氮差異顯著說明覆蓋秸稈配施定量的氮磷鉀肥能夠使土壤微生物量碳、氮有明顯的提高。處理SCKT的微生物生物量碳、氮明顯高于SNLT,說明秸稈覆蓋還田后只有當(dāng)?shù)爻渥銜r(shí)微生物才更容易繁殖。SNPT處理的微生物量碳、氮明顯低于SCKT,表明秸稈覆蓋還田配施相同的氮肥磷肥,但不施鉀肥不利于微生物生長(zhǎng)。處理SNKT微生物量碳顯著低于處理SCKT,而其微生物量氮雖也較低但卻沒有顯著性差異,說明磷肥對(duì)于土壤微生物生物量碳、氮的影響不相同。同時(shí)說明磷對(duì)于微生物的代謝具有重要的作用[16]。不同處理的下層土壤之間的差異與上層相似,且同一處理的上下層土壤的微生物生物量碳有顯著性差異,而微生物生物量氮均沒有顯著性差異。

    土壤微生物生物量碳氮比(SMBC/SMBN)可反映微生物群落結(jié)構(gòu)信息,其顯著的變化表明微生物群落結(jié)構(gòu)變化直接影響微生物量的高低[17]。一般情況下細(xì)菌的碳氮比在5∶1左右。Petra Marschnera等[18]認(rèn)為土壤有機(jī)碳的含量和碳氮比能夠顯著影響土壤細(xì)菌群落結(jié)構(gòu)的變化。不同處理間土壤微生物群落結(jié)構(gòu)的相似性或多樣性變化主要是由于秸稈還田下不同施肥土壤的微域生境的改變所致,這種改變影響了土壤生境對(duì)多種微生物的適宜性。從表2看出,微生物生物量碳氮比的范圍5.52—7.26之間。同一層土樣不同處理的微生物碳氮比并沒有差異,同一處理的上層土壤與下層之間也沒有顯著性差異。說明秸稈還田為微生物提供了充足的碳源,配施適量氮磷鉀肥增加了土壤養(yǎng)分,但是對(duì)土壤微生物生物量碳、氮的影響卻不相同,這可能是因?yàn)閺?fù)雜的環(huán)境因素對(duì)微生物有著多方面的影響。只有微生物碳氮比合適,且具有較高有效利用性,才能促進(jìn)微生物繁殖,從而增加土壤微生物生物量。

    2.2 DGGE圖譜分析

    固氮菌被認(rèn)為是土壤氮素循環(huán)中固氮環(huán)節(jié)不可或缺的重要成員。對(duì)固氮菌群落結(jié)構(gòu)和與固氮相關(guān)的種群動(dòng)態(tài)研究對(duì)于人們了解不同環(huán)境下的生態(tài)系統(tǒng)過程十分重要。許多研究試圖通過不同的分離技術(shù)去研究土壤中具有固氮功能的固氮菌,但是對(duì)大多數(shù)生態(tài)系統(tǒng)來說,細(xì)菌中不可培養(yǎng)部分占了總數(shù)的絕大多數(shù)。通過靶定固氮酶編碼基因nifH的PCR擴(kuò)增和變性梯度凝膠電泳(DGGE)分析,為深入了解固氮菌種群提供了很好的途徑[22]。

    表2 秸稈還田下不同施肥處理的土壤微生物生物量碳氮

    *CK為無秸稈覆蓋處理(N10P4K6), SNL為秸稈覆蓋+N6P4K6(低量氮肥),SCK為秸稈覆蓋+N10P4K6,SNP為秸稈覆蓋+N10P4K0,SNK為秸稈覆蓋+N10P0K6,處理代碼的下標(biāo)T表示0—10cm土層,S表示10—20cm土層;表中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差;表中同列相同小寫字母表示P≥0.05水平無顯著性差異

    圖1 土壤固氮菌nifH基因的PCR-DGGE圖譜Fig.1 PCR-DGGE profiles of soil microbial nifH genes條帶近緣菌株依次為Closest relatives in turn:1. Azospirillum brasilense strain AWB12 (GQ161239.1,90%); 2. Uncultured soil bacterium clone T01(EU331520.1, 92%); 3. Uncultured bacterium clone KW-B4(HQ335637.1,96%); 4. Uncultured bacterium clone BNnif4 (HQ190124.1, 96%); 5. Uncultured bacterium clone U1C_12(DQ520301.1,97%); 6. Pseudacidovorax sp ptl-2(GQ249664.1, 94%); 7. Azospira oryzae strain(U97115.2,90%); 8. Uncultured bacterium clone(DQ140703.1,88%); 9. Uncultured bacterium clone pPS204 (FJ381614.1,95%); 10. Uncultured Burkholderia sp isolate DGGE(HM565855.1,94%); 11. Uncultured bacterium isolate DGGE (GU097344,95%); 12. Uncultured bacterium clone BNnif22 (HQ190133.1,95%); 13. Uncultured bacterium clone BNnif17 (HQ190131.1, 98%); 14. Ideonella dechloratans (EU542578.1,93%)

    土壤中固氮酶編碼基因nifH的擴(kuò)增片段通過DGGE被分成了許多條帶,不同的條帶代表不同的固氮細(xì)菌及基因片段。結(jié)果如圖1,從圖中我們可以初步看出:供試土壤在DGGE圖譜中電泳條帶數(shù)目、強(qiáng)度和遷移率均存在一定程度的差異,充分顯示了固氮細(xì)菌的多樣性。每個(gè)泳道中條帶的亮度反映出固氮菌相對(duì)生物量的多少,電泳條帶的多少可以直觀地反映樣品中固氮菌群落的遺傳多樣性。從指紋圖譜可以看出每個(gè)樣品均可分離到10到18條不等的條帶。條帶數(shù)目最少的是CKS處理,條帶數(shù)目最多的是SCKT處理,表明配施足夠的氮磷鉀肥+秸稈覆蓋能夠維持土樣營(yíng)養(yǎng)物質(zhì)的充足,防止植物與微生物爭(zhēng)氮,更有利于固氮細(xì)菌的生長(zhǎng)繁殖。圖譜中不同土壤間有一些共同的條帶如條帶1、4,說明供試土壤間存在共同的固氮菌類群,但是這些共有條帶的亮度不相同,說明土壤固氮菌的數(shù)量有一定的改變。

    2.3 土壤固氮菌多樣性指數(shù)

    根據(jù)電泳圖譜中每個(gè)條帶的信息,對(duì)不同處理土壤中固氮菌多樣性指數(shù)(H)、豐度(S)和均勻度(EH)等指標(biāo)進(jìn)行綜合分析,結(jié)果如表3所示。結(jié)果表明,多樣性指數(shù)和豐富度最高的處理是SCKT,其值分別為2.86和18.00;而最低的是CKS,這兩個(gè)值僅為2.38和10.33。統(tǒng)計(jì)學(xué)分析后得知,處理SCKT和CKT多樣性指數(shù)及豐富度均有顯著性差異,說明覆蓋秸稈后配施氮磷鉀肥能夠明顯的提高土壤固氮菌的多樣性。處理SCKT的多樣性指數(shù)和豐富度明顯大于SNLS說明覆蓋秸稈配施相同量的磷鉀肥,氮肥施用高的土壤固氮菌多樣性指數(shù)更大。土壤生境碳氮比需要維持在一定的水平,由于配施秸稈后土壤碳素增加,只有當(dāng)配施足夠的氮才能防止微生物與植物爭(zhēng)氮,更有利于微生物的生長(zhǎng)繁殖[2]。SNPT處理的多樣性指數(shù)和豐富度明顯低于SCKT處理,表明秸稈覆蓋后只配施氮磷肥、缺少鉀肥不利于固氮菌的生長(zhǎng)。SNKT處理的多樣性指數(shù)和豐富度低于SCKT處理但沒有明顯差異,表明磷肥對(duì)于固氮菌多樣性指數(shù)有一定影響但是并沒有統(tǒng)計(jì)學(xué)上的意義。從表3可以看出,同一個(gè)處理的上下層土壤多樣性指數(shù)及豐富度均沒有顯著差異。同一處理的上層或下層土壤固氮菌的均勻度除了處理SCKT和SCKS有差異外,其余上下層處理均沒有顯著差異。不同處理的上層土壤固氮菌的均勻度中處理SCKT只與處理SNPT有顯著性差異,不同處理的下層土壤中處理SCKs與處理SNPs也有顯著性差異,說明鉀肥對(duì)于固氮細(xì)菌均勻度也有一定的影響。

    表3 秸稈還田下不同施肥處理土壤固氮基因nifH多樣性指數(shù)

    2.4 土壤固氮菌群落相似性分析

    進(jìn)一步利用Quantity One軟件(Bio-Rad)對(duì)DGGE圖譜的條帶數(shù)和亮度進(jìn)行數(shù)字化處理,并依此為據(jù)進(jìn)行聚類分析,得到各處理間的聚類分析結(jié)果(圖2) 。10個(gè)處理聚類為三大族群, CKT、CKS為一種族群;SNLT、SNLS、SNPT、SNPS為一個(gè)族群;SCKT、SCKS、SNKT、SNKS為第三族群。通過DGGE聚類圖可以看出,相同的處理的上下層土壤聚到了一起,不同處理聚在了不同的族群里,說明秸稈覆蓋和不同的施肥方式都導(dǎo)致土壤理化性質(zhì)差異,從而導(dǎo)致不同處理間的土壤固氮菌群落結(jié)構(gòu)的影響不一致。從圖中可以看出,只施肥而不覆蓋秸稈的處理單獨(dú)聚成了一族,說明覆蓋秸稈使土壤固氮菌的群落結(jié)構(gòu)發(fā)生了明顯變化。

    圖2 DGGE圖譜聚類分析Fig.2 Cluster analysis of DGGE banding patterns

    2.5 DGGE條帶測(cè)序分析

    根據(jù)土壤固氮菌DGGE圖譜中出現(xiàn)的優(yōu)勢(shì)和特殊共有條帶進(jìn)行測(cè)序,將所得序列輸入NCBI網(wǎng)站將同源性較高的序列下載,通過Clustal X1.83, Mega Soft等軟件構(gòu)建系統(tǒng)發(fā)育樹。結(jié)果如圖3。從系統(tǒng)樹可見該固氮菌群有豐富的菌種組成多樣性,14個(gè)條帶的近緣種大部分為非培養(yǎng)細(xì)菌nifH基因片段,主要優(yōu)勢(shì)菌群其歸屬于變形菌門(Proteobacteria)的β-變形菌綱(Betaproteobacteria)非培養(yǎng)伯克氏菌屬、假食酸菌屬(Pseudacidovorax)、Azospira屬、Ideonella屬。14個(gè)條帶分屬于4個(gè)主要的簇,其中第一簇中的band 10、band 11、band 12、band 13與銅礦荒地的非培養(yǎng)伯克氏菌屬固氮細(xì)菌HM565855.1 十分相似,band 3與荷蘭Kollumerwaard土壤非培養(yǎng)細(xì)菌HQ335637.1十分相似、band 14與Ideonelladechloratans菌株EU542578.1相似性為94%,第二簇中band 2、band 5、band 6 與假食酸菌屬(Pseudacidovorax)菌株聚在一起,band 1、band 4、band 7與Azospira屬菌株聚為第三簇,band 8、band 9與其它序列差異較大單獨(dú)成簇,其為未確定歸屬的固氮細(xì)菌。

    圖3 DGGE切膠條帶的系統(tǒng)發(fā)育樹,序列的GenBank登錄號(hào)列于括號(hào)中,尺長(zhǎng)度為5%的核苷酸置換率Fig.3 Phylogenetic tree of sequences from DGGE bands. Numbers in parentheses represent the sequences accession number in GenBank. The scale bar: 5% sequence divergence

    3 討論

    3.1秸稈還田能有效改善土壤物理性結(jié)構(gòu),改善各類土壤細(xì)菌營(yíng)養(yǎng)條件與提高細(xì)菌活性,有機(jī)肥與秸稈還田措施均可以提高不同功能土壤微生物類群的活性,有效的提高了氮肥的利用效率,使農(nóng)作物產(chǎn)量得到增加[19]。微生物生物量碳、氮被認(rèn)為是土壤活性養(yǎng)分的儲(chǔ)存庫(kù),是植物生長(zhǎng)可利用養(yǎng)分的重要來源[11]。土壤微生物生物量碳、氮是土壤微生物軀體中所固定的碳素和氮素,其含量的多少?zèng)Q定于土壤中微生物的數(shù)量,是土壤微生物對(duì)碳素、氮素同化與固持作用這種動(dòng)態(tài)過程的綜合反映。Nie Jun[20]等研究表明秸稈覆蓋提高了微生物生物量碳和微生物生物量氮,與只施化肥田塊相比,分別提高了12.66%和15.07%;與對(duì)照相比,分別提高了7.76%和31.42%。

    本試驗(yàn)表明秸稈覆蓋還田配施適量氮磷鉀肥土壤微生物生物量碳氮的含量明顯高于對(duì)照處理,說明秸稈還田為微生物生長(zhǎng)和繁殖提供了大量的碳源和能源,同時(shí)施用的氮肥又為微生物的生長(zhǎng)提供了較多的氮源,刺激了微生物的生長(zhǎng)。在秸稈覆蓋的處理中SCKT、SCKS分別高于SNLT、SNLS,即施氮肥較高的處理具有更高的微生物生物量碳、氮值。秸稈還田后,隨著施氮量的增加,提供給微生物的氮素增多,促進(jìn)了微生物的繁殖,加快了秸稈腐解。處理SNPT、SNKT的微生物生物量碳氮低于SCKT,處理SNPS、SNKS的微生物生物量碳氮低于SCKS。說明秸稈覆蓋后需同時(shí)配施適量氮磷鉀才能更好地為微生物的生長(zhǎng)提供更全面的能源物質(zhì)。秸稈覆蓋下不同施肥配比對(duì)表層土壤影響更大,并隨著土層加深而減少。其原因是表層土壤能較好地與外界進(jìn)行物質(zhì)與能量交換,進(jìn)入土壤的有機(jī)物質(zhì)最先聚集于表層,養(yǎng)分條件好,土壤微生物活動(dòng)旺盛,有利于微生物的繁殖。本試驗(yàn)同一處理的上下層土壤的微生物生物量碳有顯著性差異,而微生物量氮均沒有顯著性差異。說明本研究中20cm表土層內(nèi)微生物生物量碳氮均受到秸稈覆蓋及施肥的影響,但是影響的程度和方式卻不相同。

    3.2PCR-DGGE技術(shù)是從微生物基因多樣性的角度研究土壤微生物群落結(jié)構(gòu)的方法,近年來己經(jīng)成為國(guó)內(nèi)外研究的熱點(diǎn)[21]。與純培養(yǎng)相比,PCR克隆、RFLP分析等非培養(yǎng)方法被用于研究不同環(huán)境下nifH基因的組成和多樣性,更有利于對(duì)固氮生物進(jìn)行完整和深入的探索和研究。采用PCR-DGGE方法對(duì)秸稈覆蓋還田下不同施肥處理的土壤微生物nifH基因多樣性及群落結(jié)構(gòu)組成進(jìn)行研究,從指紋圖譜我們可以看出不同處理土壤固氮微生物群落結(jié)構(gòu)有一定的差異。Juraeva等研究發(fā)現(xiàn),土壤總氮含量能夠影響固氮微生物的分布[22]。Keeling等研究表明,土壤碳氮比例影響著土壤固氮微生物的活性[23]。因此,土壤理化性質(zhì)的差異可以影響固氮微生物的群落結(jié)構(gòu)。對(duì)DGGE圖譜進(jìn)行聚類分析和多樣性統(tǒng)計(jì)表明同層間SCKT、SCKS處理具有最高的多樣性指數(shù)及豐富度,而僅CKT、CKS對(duì)照組這兩個(gè)值最小。結(jié)果說明秸稈覆蓋還田配合不同施肥措施改變了土壤的營(yíng)養(yǎng)狀況,改善土壤的生物學(xué)特征,導(dǎo)致土壤固氮微生物多樣性發(fā)生變化。Baumann等研究也指出植物殘留物能夠影響微生物的呼吸和微生物群落結(jié)構(gòu)[24]。從樣品DGGE優(yōu)勢(shì)條帶測(cè)序結(jié)果可以看出,這些條帶中最多的是非培養(yǎng)固氮細(xì)菌,其余對(duì)應(yīng)的同源性最高的微生物分別屬于變形菌門(Proteobacteria),β-變形菌綱(Betaproteobacteria)的非培養(yǎng)伯克氏菌屬、假食酸菌屬(Pseudacidovorax)、Azospira屬以及Ideonella屬。說明本實(shí)驗(yàn)土壤中固氮菌遺傳多樣性較為豐富。Bodelier等也證明了像水稻這樣的淹水作物的根際上壤中氧量較低,特別適合于進(jìn)行生物固氮[25]。本研究運(yùn)用PCR-DGGE方法得到的這些固氮細(xì)菌還未被培養(yǎng)或不可培養(yǎng),故很難鑒定它們?cè)诟魍寥乐械纳鷳B(tài)功能,要具體了解這些優(yōu)勢(shì)菌群的結(jié)構(gòu)與功能,還需要做更全面細(xì)致的研究。

    [1] Jiang Y H, Yu Z R, Ma Y L. The effect of stubble return on Agro-ecological system and crop growth. Chinese Journal of Soil Science, 2001, 32(5): 209- 213.

    [2] Hong C L, Wei Y Z, Huang J F, Wang R Y, Yang X E. Effects of total crop straw return on soil fertility and field ecological environment. Journal of Zhejiang University: Agriculture and Life Sciences, 2003, 29(6): 627- 633.

    [3] Steenwerth K L, Jackson L E, Calderón F J, Stromberg M R, Scow K M. Soil microbial community composition and land use history in cultivated and grassland ecosystems of coastal California. Soil Biology and Biochemistry, 2002, 34(11): 1599- 1611.

    [4] David Tilman D, Johannes Knops J, David Wedin D, Reich P, Ritchie M, Siemann E. The influence of functional diversity and composition on ecosystem processes. Science, 1997, 277(5330): 1300- 1302

    [5] Clough Y, Kruess A, Tscharntke T. Organic versus conventional arable farming systems: functional grouping helps understand staphylinid response. Agriculture, Ecosystems and Environment, 2007, 118(1/4): 285- 290.

    [6] Yu C Y, Liu H, Xing Y D, Manukovskyd N S, Kovalevd V S, Gurevichd Y L. Bioconversion of rice straw into a soil-like substrate. Acta Astronautica 2008, 63(7/10): 1037- 1042.

    [7] Bowen R M, Harper S H T. Decomposition of related compounds by fungi isolated from straw wheat straw and in arable soi1. Soil Biology Biochemistry, 1990, 22(3): 393- 399.

    [8] Chen H, Tang X S, Lin J, Zhang B S, Ren D M. Community constitute and phylogenetic analysis on soil uncultured microorganism. Acta Microbiologica Sinica, 2002, 42(4): 478- 483.

    [9] Kisand V, Wikner J. Limited resolution of 16S rDNA DGGE caused by melting properties and closely related DNA sequences. Journal of Microbiological Methods, 2003, 54(2): 83- 191.

    [10] Dixon R, Kahn D. Genetic regulation of biological nitrogen fixation. Nature Reviews Microbiology, 2004, 2(8): 621- 631.

    [11] Wu J S, Lin Q M, Huang Q Y, Xiao H A. The Method of Measuring Soil Microbial Biomass and Its Application. Beijing: China Meteorological Press, 2006.

    [12] Dlallo M D, Wlllems A, Vloemans N, Cousin S, Vandekerckhove T T, de Lajudie P, Neyra M, Vyverman W, Gillis M, Van der Gucht K. Polymerase chain reaction denaturing gradient gel electrophoresis analysis of the N2-fixing bacterial diversity in soil under Acacia tortilis ssp. raddiana and Balanites aegyptiaca in the dryland part of Senegal. Environmental Microbiology, 2004, 6(4): 400- 415.

    [13] Eichner C A, Erb R W, Timmis K N, Wagner-D?bler J. Thermal gradient gel electrophoresis analysis of bioprotection from pollutant shocks in the activated sludge microbial community. Applied and Environmental Microbiology, 1999, 65(1): 102- 109.

    [14] Hedrick D B, Peacock A, Stephen J R, Macnaughton S J, Brüggemann J, White D C. Measuring soil microbial community diversity using polar lipid fatty acid and denaturing gradient gel electrophoresis data. Journal of Microbiological Methods, 2000, 41(3): 235- 248.

    [15] Lin X G, Hu J L. Scientific connotation and ecological service function of soil microbial diversity. Acta Pedologica Sinica, 2008, 45(5): 892- 900.

    [16] Zheng S X, Hu J L, Chen K, Yao J, Yu Z N, Lin X G. Soil microbial activity measured by microcalorimetry in response to long-term fertilization regimes and available phosphorous on heat evolution. Soil Biology and Biochemistry, 2009, 41(10): 2094- 2099.

    [17] Lovell R D, Jarvis S C, Bardgett R D. Soil microbial biomass and activity in long-term grassland: effects of management changes. Soil Biology and Biochemistry, 1995, 27(7): 969- 975.

    [18] Marschnera P, Kandeler E, Marschner B. Structure and function of the soil microbial community in a long-term fertilizer experiment. Soil Biology and Biochemistry, 2003, 35(3): 453- 461.

    [19] Ye W P, Xie X L, Wang K R, Li Z G. Effects of rice straw manuring in different periods on growth and yield of rice. Chinese Journal of Rice Science, 2008, 22(1): 65- 70.

    [20] Nie J, Zhou J M, Wang H Y, Chen X Q, Du C W. Effect of long-term rice straw return on soil Glomalin, Carbon and Nitrogen. Pedosphere, 2007, 17(3): 295- 302.

    [21] Petra M, Pauline F G, Zed R. Microbial community composition and functioning in the rhizosphere of threeBanksiaspecies in native woodland in Western Australia. Applied Soil Ecology, 2005, 28(3): 191- 201.

    [22] Juraeva D, George E, Davranov K, Ruppel S. Detection and quantification of thenifHgene in shoot and root of cucumber plants Canadian Journal of Microbiology, 2006, 52(8): 731- 739.

    [23] Keeling A A, Cook J A, Wilcox A. Effects of carbohydrate application on diazotroph populations and nitrogen availability in grass swards established in garden waste compost. Bioresource Technology, 1998, 66(2): 3814- 3822.

    [24] Baumann K, Marschner P, Smernik R J, Baldock J A. Residue chemistry and microbial community structure during decomposition of eucalypt, wheat and vetch residues. Soil Biology and Biochemistry, 2009, 41(9): 1966- 1975.

    [25] Bodelier P L E. Interaction between oxygen-releasing roots and microbial processes in flooded soils and sediments//de Kroon H, Visser E J W, eds.Root Ecology. Ecology Studies. New York: Springer, 2003, 168: 331- 362.

    參考文獻(xiàn):

    [1] 江永紅, 宇振榮, 馬永良. 秸稈還田對(duì)農(nóng)田生態(tài)系統(tǒng)及作物生長(zhǎng)的影響. 土壤通報(bào), 2001, 32(5): 209- 213.

    [2] 洪春來, 魏幼璋, 黃錦法, 王潤(rùn)屹, 楊肖娥. 秸稈全量直接還田對(duì)土壤肥力及農(nóng)田生態(tài)環(huán)境的影響研究. 浙江大學(xué)學(xué)報(bào): 農(nóng)業(yè)與生命科學(xué)版, 2003, 29(6): 627- 633.

    [8] 陳灝, 唐小樹, 林潔, 張伯生, 任大明. 不經(jīng)培養(yǎng)的農(nóng)田土壤微生物種群構(gòu)成及系統(tǒng)分類的初步研究. 微生物學(xué)報(bào), 2002, 42(4): 478- 483.

    [11] 吳金水, 林啟美, 黃巧云, 肖和艾. 土壤微生物量測(cè)定方法及其應(yīng)用. 北京: 氣象出版社, 2006: 54- 61.

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

    [19] 葉文培, 謝小立, 王凱榮, 李志國(guó). 不同時(shí)期秸稈還田對(duì)水稻生長(zhǎng)發(fā)育及產(chǎn)量的影響. 中國(guó)水稻科學(xué), 2008, 22(1): 65- 70.

    Effectofdifferentfertilizercombinationsandstrawreturnonmicrobialbiomassandnitrogen-fixingbacteriacommunityinapaddysoil

    LIU Xiaoqian1, TU Shihua2, SUN Xifa2, GU Yunfu1, ZHANG Xianqin1, ZHANG Xiaoping1,*

    1CollegeofResourcesandEnvironment,SichuanAgriculturalUniversity,Chengdu611130,China2SoilandFertilizerInstitute,SichuanAcademyofAgriculturalSciences,Chengdu610066,China

    Chloroform fumigation and denaturing gradient gel electrophoresis(DGGE) were employed to study the influence of different fertilizer treatments with crop residue return on microbial biomass and nitrogen-fixing bacteria community in a gray brown alluvial paddy soil at soil depth of 0—10 cm and 10—20cm. Results show that soil microbial biomass carbon (SMB-C), microbial biomass nitrogen (SMB-N) and diversity of nitrogen-fixing bacteria decreased with soil depths. The control treatment with straw removal had the lowest amounts of SMB-C and SMB-N. Among the treatments receiving straw return, the quantity of SMB-C and SMB-N in the treatments with low nitrogen rate and omitting was significantly lower than the full rate NPK treatment. Omission of phosphate on the basis of NPK nonsignificantly reduced the SMB-N compared to the NPK treatment. It shows that straw mulching with addition of adequate NPK fertilizers can significantly increase the SMB-C and SMB-N. The diversity indices of DGGE profiles indicate that the treatment with adequate NPK had the most abundant nitrogen-fixing bacterial diversity. According to UPGMA cluster analysis, dendrogram of the 10 treatments were different, and the micro-organisms at the 0—10cm and 10—20cm in the control treatment differed from all the other treatments and stood alone in a separate group. DGGE bands sequencing further revealed that the 14 closely related species were largely uncultured bacterianifHgene fragment with the predominant community falling in the class of Proteobacteria in Betaproteobacteria. Application of PCR-DGGE can well interpret the community structure characteristics of the nitrogen-fixing bacteria molecules in the gray-brown alluvial paddy soil with straw return and different fertilizer treatments.

    straw return; gray-brown alluvial paddy soil; microbial biomass; PCR-DGGE; nitrogen-fixing bacteria community

    農(nóng)業(yè)部公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)資助(201003016)

    2012- 05- 29;

    2012- 10- 26

    *通訊作者Corresponding author.E-mail: luwenru@mail.caas.net.cn

    10.5846/stxb201205290788

    劉驍蒨,涂仕華,孫錫發(fā),辜運(yùn)富,張先琴,張小平.秸稈還田與施肥對(duì)稻田土壤微生物生物量及固氮菌群落結(jié)構(gòu)的影響.生態(tài)學(xué)報(bào),2013,33(17):5210- 5218.

    Liu X Q, Tu S H, Sun X F, Gu Y F, Zhang X Q, Zhang X P.Effect of different fertilizer combinations and straw return on microbial biomass and nitrogen-fixing bacteria community in a paddy soil.Acta Ecologica Sinica,2013,33(17):5210- 5218.

    猜你喜歡
    固氮菌固氮條帶
    解磷菌、解鉀菌和固氮菌的分離篩選與鑒定
    土壤中的天然化肥廠
    施用固氮菌肥的注意事項(xiàng)
    基于條帶模式GEOSAR-TOPS模式UAVSAR的雙基成像算法
    杉木與固氮樹種混交對(duì)土壤有機(jī)質(zhì)及氮含量的影響
    土壤中的天然化肥廠
    ——固氮微生物
    基于 Savitzky-Golay 加權(quán)擬合的紅外圖像非均勻性條帶校正方法
    一種基于MATLAB的聲吶條帶圖像自動(dòng)拼接算法
    海岸工程(2014年4期)2014-02-27 12:51:28
    一株具有固氮功能的煙草根際微生物的鑒定及其初步效應(yīng)
    一株具有ACC脫氨酶活性固氮菌的篩選與鑒定
    久久午夜综合久久蜜桃| 免费看美女性在线毛片视频| 午夜影院日韩av| 人妻夜夜爽99麻豆av| 99国产精品一区二区三区| 国产在线精品亚洲第一网站| 久久伊人香网站| 成熟少妇高潮喷水视频| 久久久久国产一级毛片高清牌| 亚洲五月婷婷丁香| 日韩高清综合在线| 99国产精品一区二区蜜桃av| 99在线人妻在线中文字幕| 亚洲国产看品久久| 国产高清视频在线播放一区| 午夜激情av网站| 黑人操中国人逼视频| 精品国产乱子伦一区二区三区| 亚洲第一欧美日韩一区二区三区| 成人av在线播放网站| 亚洲一码二码三码区别大吗| 国产三级黄色录像| 男女做爰动态图高潮gif福利片| 99久久久亚洲精品蜜臀av| 国产精品亚洲美女久久久| 欧美大码av| 50天的宝宝边吃奶边哭怎么回事| 日日干狠狠操夜夜爽| 1024香蕉在线观看| 亚洲乱码一区二区免费版| 欧美三级亚洲精品| 欧美日韩国产亚洲二区| 精品第一国产精品| 国产视频内射| 高清毛片免费观看视频网站| 久久人妻福利社区极品人妻图片| 欧美午夜高清在线| 巨乳人妻的诱惑在线观看| 国产一级毛片七仙女欲春2| 欧美乱色亚洲激情| 99热这里只有是精品50| 中出人妻视频一区二区| 欧美日韩福利视频一区二区| 日本 欧美在线| 嫩草影院精品99| 黄片小视频在线播放| 国产av在哪里看| 精品国产超薄肉色丝袜足j| 日韩欧美在线乱码| 精品一区二区三区av网在线观看| 色综合欧美亚洲国产小说| 最新在线观看一区二区三区| 最近在线观看免费完整版| 亚洲国产欧美一区二区综合| 久久久久久国产a免费观看| 美女黄网站色视频| 丰满的人妻完整版| 色噜噜av男人的天堂激情| 两个人的视频大全免费| 国产黄片美女视频| 国产精品国产高清国产av| 国产熟女xx| 国产1区2区3区精品| 国产av一区在线观看免费| 成人国语在线视频| 在线观看日韩欧美| 香蕉国产在线看| cao死你这个sao货| 国产亚洲精品久久久久5区| a级毛片a级免费在线| 免费电影在线观看免费观看| 啦啦啦免费观看视频1| av天堂在线播放| 欧美一区二区精品小视频在线| 在线a可以看的网站| 曰老女人黄片| 欧美绝顶高潮抽搐喷水| 成人欧美大片| 淫妇啪啪啪对白视频| 久久热在线av| 国产精品永久免费网站| 人妻久久中文字幕网| 三级毛片av免费| 亚洲成人久久爱视频| 两性夫妻黄色片| 91国产中文字幕| 高清毛片免费观看视频网站| 岛国在线免费视频观看| a级毛片a级免费在线| 国产午夜精品久久久久久| 天堂动漫精品| 露出奶头的视频| 日本五十路高清| 亚洲精品在线观看二区| 久久久久亚洲av毛片大全| 国产三级中文精品| 亚洲美女黄片视频| 男女午夜视频在线观看| 色在线成人网| 日本五十路高清| 欧美又色又爽又黄视频| 精品电影一区二区在线| 国产99白浆流出| 嫩草影视91久久| 91字幕亚洲| 久久婷婷人人爽人人干人人爱| 巨乳人妻的诱惑在线观看| 男女下面进入的视频免费午夜| 他把我摸到了高潮在线观看| 在线永久观看黄色视频| 亚洲精品色激情综合| 久久久水蜜桃国产精品网| 亚洲精品粉嫩美女一区| 九色成人免费人妻av| 麻豆成人午夜福利视频| 亚洲一码二码三码区别大吗| 日本a在线网址| 亚洲精品中文字幕一二三四区| 国产亚洲精品久久久久5区| 天天添夜夜摸| 国产精品日韩av在线免费观看| 亚洲熟妇熟女久久| 免费高清视频大片| 久久精品影院6| 天堂√8在线中文| 亚洲国产欧美人成| 亚洲狠狠婷婷综合久久图片| 久久久国产精品麻豆| 免费看美女性在线毛片视频| 最近最新中文字幕大全免费视频| 黄频高清免费视频| 日韩av在线大香蕉| 国产高清有码在线观看视频 | 露出奶头的视频| or卡值多少钱| 精品国产乱码久久久久久男人| 久久久久久久精品吃奶| 18美女黄网站色大片免费观看| 亚洲五月婷婷丁香| 亚洲色图 男人天堂 中文字幕| 美女 人体艺术 gogo| 欧美3d第一页| 美女扒开内裤让男人捅视频| svipshipincom国产片| 精品第一国产精品| 成人欧美大片| 久久精品国产综合久久久| 婷婷亚洲欧美| 亚洲欧美精品综合一区二区三区| 久久国产精品影院| 国产一级毛片七仙女欲春2| 悠悠久久av| 午夜影院日韩av| 国产精品九九99| 国内毛片毛片毛片毛片毛片| 中亚洲国语对白在线视频| 中文字幕精品亚洲无线码一区| 午夜成年电影在线免费观看| 高清在线国产一区| 免费在线观看日本一区| 真人做人爱边吃奶动态| 亚洲欧美日韩东京热| 麻豆久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆 | 欧美性猛交╳xxx乱大交人| 久9热在线精品视频| 特大巨黑吊av在线直播| 欧美国产日韩亚洲一区| 制服人妻中文乱码| 国产伦一二天堂av在线观看| av免费在线观看网站| 精品人妻1区二区| 天堂动漫精品| 两性夫妻黄色片| 男女那种视频在线观看| 国产精品免费视频内射| 欧美在线黄色| 男女之事视频高清在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲精品av麻豆狂野| ponron亚洲| 亚洲av成人av| 国模一区二区三区四区视频 | 国产97色在线日韩免费| 成人欧美大片| 久久久久久免费高清国产稀缺| 久久这里只有精品中国| 国产黄色小视频在线观看| 免费在线观看成人毛片| 中文字幕久久专区| 久久久久久免费高清国产稀缺| 两个人免费观看高清视频| 我要搜黄色片| 亚洲人成伊人成综合网2020| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽| 久久久久久久久久黄片| 一个人观看的视频www高清免费观看 | 国产在线精品亚洲第一网站| 在线播放国产精品三级| 日韩中文字幕欧美一区二区| 国产av不卡久久| 一进一出抽搐动态| 99精品在免费线老司机午夜| 午夜视频精品福利| 床上黄色一级片| 亚洲成人久久爱视频| 丁香六月欧美| 日韩免费av在线播放| 国产亚洲精品久久久久5区| 欧美黄色淫秽网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区三区四区久久| 成年人黄色毛片网站| 午夜精品一区二区三区免费看| 国产一区二区三区视频了| 波多野结衣高清作品| 国产真实乱freesex| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 亚洲av第一区精品v没综合| 精品欧美一区二区三区在线| 亚洲一区中文字幕在线| 亚洲欧美日韩无卡精品| 露出奶头的视频| 欧美另类亚洲清纯唯美| 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 午夜精品在线福利| 国产av麻豆久久久久久久| 18禁观看日本| 国产高清有码在线观看视频 | 国产精品一区二区精品视频观看| 欧美黑人精品巨大| 日韩高清综合在线| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 88av欧美| 毛片女人毛片| 丁香六月欧美| 神马国产精品三级电影在线观看 | 一级毛片高清免费大全| 999久久久国产精品视频| 亚洲无线在线观看| 国产麻豆成人av免费视频| 窝窝影院91人妻| 天天躁夜夜躁狠狠躁躁| 最近视频中文字幕2019在线8| 久久久国产成人精品二区| 老司机在亚洲福利影院| 久久天堂一区二区三区四区| 青草久久国产| 看片在线看免费视频| 欧美又色又爽又黄视频| 免费观看人在逋| 午夜成年电影在线免费观看| 亚洲熟女毛片儿| 波多野结衣高清作品| 在线观看日韩欧美| 欧美日本视频| 国产一区在线观看成人免费| 中亚洲国语对白在线视频| 久久久久久久久免费视频了| ponron亚洲| 欧美性猛交╳xxx乱大交人| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 日韩欧美三级三区| 精品熟女少妇八av免费久了| 亚洲中文日韩欧美视频| 亚洲乱码一区二区免费版| 亚洲一码二码三码区别大吗| 成人18禁高潮啪啪吃奶动态图| 亚洲一区二区三区色噜噜| 国产一区二区激情短视频| 老司机午夜十八禁免费视频| 国产熟女xx| 国产不卡一卡二| 国模一区二区三区四区视频 | 久久久久免费精品人妻一区二区| 亚洲国产精品999在线| 一级片免费观看大全| 久久久久久久久免费视频了| 亚洲专区中文字幕在线| 亚洲人成77777在线视频| 婷婷丁香在线五月| 熟女电影av网| 老汉色av国产亚洲站长工具| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 啦啦啦免费观看视频1| 又大又爽又粗| 黄色毛片三级朝国网站| 国产野战对白在线观看| 男女视频在线观看网站免费 | 99精品欧美一区二区三区四区| 日韩有码中文字幕| 中文资源天堂在线| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 精华霜和精华液先用哪个| 欧美乱妇无乱码| 精品福利观看| av超薄肉色丝袜交足视频| 国产伦一二天堂av在线观看| 国产三级黄色录像| 久久久久国产一级毛片高清牌| 久久国产精品影院| 午夜免费观看网址| 18禁黄网站禁片午夜丰满| 亚洲精品在线美女| 久久九九热精品免费| 午夜激情av网站| 亚洲av电影在线进入| 高潮久久久久久久久久久不卡| 亚洲激情在线av| 亚洲国产中文字幕在线视频| www日本黄色视频网| 欧美成狂野欧美在线观看| 国产激情偷乱视频一区二区| 精华霜和精华液先用哪个| 国模一区二区三区四区视频 | 欧美绝顶高潮抽搐喷水| 久久中文看片网| 精品一区二区三区四区五区乱码| 曰老女人黄片| 国产av不卡久久| 成年女人毛片免费观看观看9| 精品少妇一区二区三区视频日本电影| 国产av不卡久久| 精品久久蜜臀av无| 日本黄大片高清| 丝袜美腿诱惑在线| 久久国产精品人妻蜜桃| 伊人久久大香线蕉亚洲五| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清| 欧美av亚洲av综合av国产av| 俺也久久电影网| 亚洲国产日韩欧美精品在线观看 | 国产免费男女视频| 男女床上黄色一级片免费看| 一卡2卡三卡四卡精品乱码亚洲| 好男人电影高清在线观看| 免费人成视频x8x8入口观看| 成人欧美大片| 亚洲精品久久成人aⅴ小说| 国产97色在线日韩免费| 久久精品91蜜桃| а√天堂www在线а√下载| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av在线| 中国美女看黄片| 色综合婷婷激情| 1024香蕉在线观看| 国产日本99.免费观看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品美女久久av网站| 欧美中文日本在线观看视频| 在线观看www视频免费| 国内毛片毛片毛片毛片毛片| 久久热在线av| 一级a爱片免费观看的视频| 嫁个100分男人电影在线观看| 亚洲av片天天在线观看| 免费av毛片视频| av片东京热男人的天堂| 日日干狠狠操夜夜爽| 久久午夜综合久久蜜桃| 好男人电影高清在线观看| 欧美在线一区亚洲| 中文字幕高清在线视频| 日韩 欧美 亚洲 中文字幕| 99久久国产精品久久久| 精品福利观看| 露出奶头的视频| 久久久久精品国产欧美久久久| 18禁观看日本| 中亚洲国语对白在线视频| 一级毛片高清免费大全| 啦啦啦免费观看视频1| 制服诱惑二区| 精品福利观看| 99国产综合亚洲精品| 欧美一级a爱片免费观看看 | 色综合欧美亚洲国产小说| 国产爱豆传媒在线观看 | 黑人巨大精品欧美一区二区mp4| 亚洲九九香蕉| 欧美最黄视频在线播放免费| 国产高清激情床上av| 久久这里只有精品19| 亚洲成人免费电影在线观看| 男女之事视频高清在线观看| 黑人欧美特级aaaaaa片| 搡老岳熟女国产| av福利片在线| 曰老女人黄片| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 午夜福利欧美成人| 一本一本综合久久| 亚洲欧美日韩东京热| a级毛片a级免费在线| 国内精品久久久久精免费| 国产主播在线观看一区二区| 精品一区二区三区四区五区乱码| 在线a可以看的网站| 99国产综合亚洲精品| 欧美久久黑人一区二区| 丁香欧美五月| 亚洲人成伊人成综合网2020| 动漫黄色视频在线观看| 欧美在线一区亚洲| 欧美日韩国产亚洲二区| 亚洲熟妇中文字幕五十中出| 女人被狂操c到高潮| 日本免费一区二区三区高清不卡| 在线观看午夜福利视频| 日韩精品中文字幕看吧| a级毛片a级免费在线| 亚洲国产精品久久男人天堂| 亚洲欧美日韩无卡精品| 国产精品美女特级片免费视频播放器 | 日韩欧美国产在线观看| 国产精品久久久av美女十八| 国产成人aa在线观看| 妹子高潮喷水视频| 男女视频在线观看网站免费 | 国产亚洲精品久久久久5区| 一级毛片精品| 亚洲av美国av| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 无人区码免费观看不卡| 久久精品91蜜桃| a级毛片在线看网站| 欧美性长视频在线观看| 久久久精品欧美日韩精品| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 校园春色视频在线观看| 亚洲欧美日韩无卡精品| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 一级片免费观看大全| 欧美绝顶高潮抽搐喷水| 一个人免费在线观看的高清视频| 国内精品久久久久精免费| 久久99热这里只有精品18| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清无吗| 午夜精品久久久久久毛片777| 亚洲人成网站在线播放欧美日韩| 亚洲aⅴ乱码一区二区在线播放 | 精品国产乱码久久久久久男人| 国产av麻豆久久久久久久| 美女高潮喷水抽搐中文字幕| 亚洲第一电影网av| 国产成人欧美在线观看| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久| 天天添夜夜摸| 亚洲av片天天在线观看| videosex国产| 国产激情偷乱视频一区二区| 亚洲精品美女久久久久99蜜臀| 制服人妻中文乱码| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 久久国产乱子伦精品免费另类| 国内精品久久久久久久电影| 久久精品人妻少妇| 91国产中文字幕| 熟女电影av网| 日韩高清综合在线| 一二三四在线观看免费中文在| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品中文字幕一二三四区| 99久久99久久久精品蜜桃| 久久久久久久午夜电影| 亚洲 欧美一区二区三区| 亚洲18禁久久av| 男女那种视频在线观看| 宅男免费午夜| 午夜久久久久精精品| 久久天堂一区二区三区四区| svipshipincom国产片| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲人成网站在线播放欧美日韩| 国产精品久久视频播放| 亚洲欧美精品综合一区二区三区| 亚洲五月天丁香| 亚洲专区国产一区二区| 免费看十八禁软件| 麻豆成人av在线观看| 午夜福利在线观看吧| 亚洲欧美精品综合一区二区三区| 一个人观看的视频www高清免费观看 | 久久久国产成人免费| 老司机靠b影院| 久9热在线精品视频| 我的老师免费观看完整版| 亚洲色图 男人天堂 中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 国产精品日韩av在线免费观看| 日日爽夜夜爽网站| 国产一区二区在线av高清观看| 可以在线观看毛片的网站| 麻豆国产97在线/欧美 | 亚洲成人精品中文字幕电影| 亚洲一区高清亚洲精品| 久久久国产欧美日韩av| 狂野欧美激情性xxxx| 18禁黄网站禁片免费观看直播| 可以免费在线观看a视频的电影网站| 色在线成人网| 美女扒开内裤让男人捅视频| 99re在线观看精品视频| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 欧美中文综合在线视频| 久久这里只有精品中国| svipshipincom国产片| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站 | 精品第一国产精品| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 欧美一区二区国产精品久久精品 | 精品人妻1区二区| 欧美黄色淫秽网站| 神马国产精品三级电影在线观看 | 婷婷亚洲欧美| 99国产精品一区二区蜜桃av| 人人妻人人看人人澡| 1024香蕉在线观看| 50天的宝宝边吃奶边哭怎么回事| 长腿黑丝高跟| 成人精品一区二区免费| 听说在线观看完整版免费高清| 国产成+人综合+亚洲专区| 久久精品国产亚洲av高清一级| 午夜激情av网站| 午夜福利视频1000在线观看| 国产99白浆流出| 免费在线观看黄色视频的| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 亚洲国产精品sss在线观看| 色综合欧美亚洲国产小说| 日本 av在线| 制服人妻中文乱码| 久久精品影院6| 国产亚洲精品综合一区在线观看 | 熟妇人妻久久中文字幕3abv| 国产成人精品久久二区二区91| 亚洲欧美精品综合一区二区三区| 男插女下体视频免费在线播放| 亚洲成人久久性| 99热6这里只有精品| 亚洲美女黄片视频| 欧美成人免费av一区二区三区| 日韩精品免费视频一区二区三区| 蜜桃久久精品国产亚洲av| 宅男免费午夜| 日日摸夜夜添夜夜添小说| av在线天堂中文字幕| 757午夜福利合集在线观看| 亚洲五月天丁香| 亚洲中文字幕一区二区三区有码在线看 | 国内毛片毛片毛片毛片毛片| 久久久久国产精品人妻aⅴ院| 午夜福利在线观看吧| 身体一侧抽搐| aaaaa片日本免费| 亚洲国产看品久久| 中文字幕最新亚洲高清| 九色成人免费人妻av| 国产伦一二天堂av在线观看| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 好男人在线观看高清免费视频| 久久欧美精品欧美久久欧美| 成年女人毛片免费观看观看9| 婷婷六月久久综合丁香| 亚洲av成人精品一区久久| 久久久精品欧美日韩精品| 岛国在线观看网站| 一本久久中文字幕| 男人舔女人的私密视频| 免费看日本二区| 91麻豆精品激情在线观看国产| 窝窝影院91人妻| 国产伦人伦偷精品视频| 女同久久另类99精品国产91| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品一区av在线观看| 蜜桃久久精品国产亚洲av| 中文字幕人妻丝袜一区二区| 欧美高清成人免费视频www| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽| 亚洲人与动物交配视频| 巨乳人妻的诱惑在线观看| 制服人妻中文乱码| 琪琪午夜伦伦电影理论片6080| 不卡av一区二区三区| 嫁个100分男人电影在线观看| 91字幕亚洲| 国产成人啪精品午夜网站| АⅤ资源中文在线天堂| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看 |