• <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脫氨酶活性固氮菌的篩選與鑒定
    国产成人精品久久久久久| 色吧在线观看| 国产午夜精品久久久久久一区二区三区| 日韩一本色道免费dvd| av天堂中文字幕网| 寂寞人妻少妇视频99o| 久久精品国产自在天天线| 国产精品爽爽va在线观看网站| 观看美女的网站| av国产精品久久久久影院| 天美传媒精品一区二区| 啦啦啦啦在线视频资源| 亚洲美女搞黄在线观看| 中国美白少妇内射xxxbb| 中文字幕人妻熟人妻熟丝袜美| 欧美精品国产亚洲| 美女xxoo啪啪120秒动态图| 伊人亚洲综合成人网| 欧美人与善性xxx| 大香蕉久久网| 精品国产露脸久久av麻豆| 国产精品偷伦视频观看了| www.色视频.com| 亚洲av电影在线观看一区二区三区| 成人无遮挡网站| 青春草国产在线视频| 亚洲少妇的诱惑av| 久久久精品免费免费高清| 人人妻人人澡人人爽人人夜夜| 男女啪啪激烈高潮av片| 狠狠婷婷综合久久久久久88av| 超色免费av| 黑人高潮一二区| 久久人人爽av亚洲精品天堂| 国产亚洲精品久久久com| 国产一级毛片在线| 午夜影院在线不卡| 色94色欧美一区二区| 汤姆久久久久久久影院中文字幕| 黑人高潮一二区| 99久久综合免费| 免费看av在线观看网站| 国产一区有黄有色的免费视频| 国产成人精品福利久久| 一级,二级,三级黄色视频| 最黄视频免费看| 久久久久久久久大av| 亚洲欧洲精品一区二区精品久久久 | av免费在线看不卡| 大片电影免费在线观看免费| 黑人欧美特级aaaaaa片| 少妇猛男粗大的猛烈进出视频| 亚洲av免费高清在线观看| 成人二区视频| 精品亚洲成国产av| 国产精品人妻久久久久久| 91精品一卡2卡3卡4卡| 插逼视频在线观看| 国产熟女午夜一区二区三区 | 精品酒店卫生间| 免费观看a级毛片全部| 青青草视频在线视频观看| 韩国高清视频一区二区三区| 成人免费观看视频高清| 最近最新中文字幕免费大全7| av卡一久久| 亚洲精品乱码久久久v下载方式| 久久久国产精品麻豆| 最新中文字幕久久久久| 插阴视频在线观看视频| 久久精品熟女亚洲av麻豆精品| 满18在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 久久狼人影院| freevideosex欧美| av国产精品久久久久影院| 国产视频内射| 亚洲无线观看免费| 亚洲欧洲国产日韩| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 成人国产麻豆网| 热re99久久国产66热| av线在线观看网站| 国产成人精品无人区| 精品久久久久久久久av| 丰满饥渴人妻一区二区三| 爱豆传媒免费全集在线观看| 亚洲综合精品二区| 免费黄频网站在线观看国产| 91精品一卡2卡3卡4卡| 尾随美女入室| 国产男人的电影天堂91| 激情五月婷婷亚洲| 天美传媒精品一区二区| 免费黄色在线免费观看| 国产日韩欧美视频二区| 国产 精品1| 99热这里只有是精品在线观看| 久久精品国产鲁丝片午夜精品| 国产成人午夜福利电影在线观看| 五月开心婷婷网| 日韩人妻高清精品专区| 亚洲欧美成人综合另类久久久| 精品久久久久久久久av| 午夜福利视频精品| 一级黄片播放器| 新久久久久国产一级毛片| 午夜福利,免费看| 日产精品乱码卡一卡2卡三| 欧美性感艳星| av在线app专区| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 18禁观看日本| 国产成人精品无人区| 激情五月婷婷亚洲| 特大巨黑吊av在线直播| 男女免费视频国产| 亚洲精品久久午夜乱码| 丁香六月天网| 蜜臀久久99精品久久宅男| 午夜免费鲁丝| 国产一区有黄有色的免费视频| 久久久久久久久久久久大奶| 亚洲色图综合在线观看| 美女中出高潮动态图| 一级,二级,三级黄色视频| 18禁在线播放成人免费| 国产日韩欧美亚洲二区| 精品国产露脸久久av麻豆| 亚洲伊人久久精品综合| 久久精品夜色国产| 免费av不卡在线播放| 亚洲精品色激情综合| 黄色视频在线播放观看不卡| 黄色配什么色好看| 菩萨蛮人人尽说江南好唐韦庄| 99久国产av精品国产电影| 欧美精品高潮呻吟av久久| 免费观看在线日韩| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影| 男女免费视频国产| 啦啦啦中文免费视频观看日本| 高清毛片免费看| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 高清毛片免费看| 国产成人精品久久久久久| 欧美国产精品一级二级三级| 内地一区二区视频在线| 精品人妻偷拍中文字幕| 国产黄片视频在线免费观看| av在线老鸭窝| 亚洲色图 男人天堂 中文字幕 | 亚洲情色 制服丝袜| 中文字幕免费在线视频6| 欧美日韩一区二区视频在线观看视频在线| 丰满迷人的少妇在线观看| 丝袜在线中文字幕| 看十八女毛片水多多多| 免费av中文字幕在线| 亚洲精品中文字幕在线视频| 大陆偷拍与自拍| 国产亚洲午夜精品一区二区久久| 亚洲精品乱码久久久v下载方式| 三级国产精品欧美在线观看| 成人漫画全彩无遮挡| 亚洲精品乱码久久久v下载方式| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| av在线老鸭窝| 26uuu在线亚洲综合色| 一本久久精品| 亚洲,欧美,日韩| 七月丁香在线播放| 香蕉精品网在线| 最近中文字幕2019免费版| a级毛片免费高清观看在线播放| 国产欧美日韩一区二区三区在线 | 韩国av在线不卡| 国产免费现黄频在线看| 国产精品一区二区在线不卡| 国产在线一区二区三区精| 永久免费av网站大全| 色视频在线一区二区三区| 在线观看人妻少妇| 日韩中字成人| 欧美xxxx性猛交bbbb| 看十八女毛片水多多多| 亚洲内射少妇av| 91精品国产九色| 免费观看在线日韩| 综合色丁香网| 免费少妇av软件| 午夜av观看不卡| 春色校园在线视频观看| 国产欧美亚洲国产| 激情五月婷婷亚洲| 高清在线视频一区二区三区| 男人操女人黄网站| 国产成人aa在线观看| 亚洲精品美女久久av网站| 国产一区亚洲一区在线观看| 午夜精品国产一区二区电影| 美女大奶头黄色视频| 久久韩国三级中文字幕| 色5月婷婷丁香| 成人国产av品久久久| 国产免费福利视频在线观看| 啦啦啦视频在线资源免费观看| 欧美97在线视频| 精品久久久久久久久亚洲| 亚洲av免费高清在线观看| 不卡视频在线观看欧美| 成人国语在线视频| 中文字幕亚洲精品专区| 母亲3免费完整高清在线观看 | 3wmmmm亚洲av在线观看| 亚洲精品第二区| 免费高清在线观看日韩| 久久久久人妻精品一区果冻| av在线app专区| 欧美亚洲 丝袜 人妻 在线| 精品国产国语对白av| 另类亚洲欧美激情| 777米奇影视久久| 成人二区视频| 美女大奶头黄色视频| 亚洲成人手机| 国产一区二区三区综合在线观看 | 欧美日韩国产mv在线观看视频| 国产午夜精品一二区理论片| 国产欧美另类精品又又久久亚洲欧美| 高清毛片免费看| 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 国产伦理片在线播放av一区| 高清av免费在线| a级毛色黄片| 久久鲁丝午夜福利片| 精品久久久久久久久av| 精品亚洲乱码少妇综合久久| 男人操女人黄网站| 在线免费观看不下载黄p国产| 欧美日韩综合久久久久久| 人人澡人人妻人| 国产一区有黄有色的免费视频| 国产熟女午夜一区二区三区 | 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| 日本黄色日本黄色录像| 国产精品女同一区二区软件| 欧美变态另类bdsm刘玥| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 国产一区二区在线观看日韩| 婷婷成人精品国产| 亚洲精品视频女| 春色校园在线视频观看| .国产精品久久| 9色porny在线观看| 午夜视频国产福利| 高清av免费在线| 亚洲国产精品国产精品| 久久99精品国语久久久| 日韩亚洲欧美综合| 亚洲成色77777| 18禁在线播放成人免费| kizo精华| 国产片内射在线| 又黄又爽又刺激的免费视频.| 亚洲熟女精品中文字幕| 在线免费观看不下载黄p国产| 日本-黄色视频高清免费观看| 一二三四中文在线观看免费高清| 久久av网站| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 亚洲精品,欧美精品| 中文字幕av电影在线播放| 久久久久久久亚洲中文字幕| 少妇人妻久久综合中文| 中文字幕免费在线视频6| 欧美日韩成人在线一区二区| 秋霞在线观看毛片| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 欧美日韩一区二区视频在线观看视频在线| 成人无遮挡网站| 国产在线视频一区二区| 日韩成人av中文字幕在线观看| 午夜视频国产福利| 久久久精品94久久精品| 日韩电影二区| 自拍欧美九色日韩亚洲蝌蚪91| 乱人伦中国视频| 婷婷色综合大香蕉| 特大巨黑吊av在线直播| 国产精品无大码| 狂野欧美激情性bbbbbb| 久久热精品热| 国产精品一区www在线观看| 熟女电影av网| 国产综合精华液| 免费看不卡的av| 我的女老师完整版在线观看| 久久精品国产鲁丝片午夜精品| 3wmmmm亚洲av在线观看| 在线亚洲精品国产二区图片欧美 | 乱人伦中国视频| 日韩中字成人| 在线亚洲精品国产二区图片欧美 | 在线观看人妻少妇| 国产精品国产三级国产av玫瑰| 99久久人妻综合| 日本午夜av视频| 在线观看www视频免费| 天天躁夜夜躁狠狠久久av| 夜夜骑夜夜射夜夜干| 国产成人精品在线电影| 精品少妇内射三级| 国产精品一区二区在线观看99| 久久久久人妻精品一区果冻| 高清毛片免费看| 久久热精品热| 婷婷色综合大香蕉| 亚洲色图 男人天堂 中文字幕 | 久久久久精品久久久久真实原创| 丰满少妇做爰视频| 日韩视频在线欧美| 免费人妻精品一区二区三区视频| 91精品国产国语对白视频| 国产色婷婷99| 成人黄色视频免费在线看| 人妻人人澡人人爽人人| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 精品人妻熟女毛片av久久网站| 99久久中文字幕三级久久日本| 久久精品夜色国产| 波野结衣二区三区在线| 国国产精品蜜臀av免费| 熟妇人妻不卡中文字幕| 午夜福利网站1000一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 中文字幕人妻丝袜制服| 亚洲欧美清纯卡通| 天堂中文最新版在线下载| 午夜91福利影院| 国产免费现黄频在线看| 美女内射精品一级片tv| 免费av不卡在线播放| 91精品伊人久久大香线蕉| 人妻制服诱惑在线中文字幕| 日日爽夜夜爽网站| 国产免费现黄频在线看| 大片免费播放器 马上看| 久久久久国产网址| 一二三四中文在线观看免费高清| 女性生殖器流出的白浆| 99九九在线精品视频| 色94色欧美一区二区| 少妇猛男粗大的猛烈进出视频| 国产成人午夜福利电影在线观看| 亚洲熟女精品中文字幕| 亚洲一区二区三区欧美精品| 欧美性感艳星| 自线自在国产av| 高清黄色对白视频在线免费看| 国精品久久久久久国模美| 少妇的逼水好多| 最新中文字幕久久久久| 一本一本综合久久| 岛国毛片在线播放| 青春草亚洲视频在线观看| 在线亚洲精品国产二区图片欧美 | 亚洲精品乱久久久久久| 亚洲国产精品成人久久小说| 波野结衣二区三区在线| 有码 亚洲区| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| www.色视频.com| 欧美成人精品欧美一级黄| 日韩人妻高清精品专区| 我要看黄色一级片免费的| 欧美+日韩+精品| 国产亚洲欧美精品永久| 丝袜在线中文字幕| 午夜福利在线观看免费完整高清在| 黑人欧美特级aaaaaa片| 精品亚洲成a人片在线观看| 69精品国产乱码久久久| 韩国av在线不卡| 欧美精品人与动牲交sv欧美| 黄色怎么调成土黄色| 精品亚洲成国产av| av福利片在线| 欧美精品亚洲一区二区| 亚洲在久久综合| 97超视频在线观看视频| 国产 精品1| 99久久精品一区二区三区| 久久精品久久久久久噜噜老黄| 秋霞在线观看毛片| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美在线一区| 久热久热在线精品观看| 欧美日韩一区二区视频在线观看视频在线| 国产成人av激情在线播放 | 日本黄大片高清| 久久午夜福利片| 久久久久国产网址| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品色激情综合| 黑人巨大精品欧美一区二区蜜桃 | 少妇人妻久久综合中文| 99热这里只有是精品在线观看| 满18在线观看网站| 国产免费一区二区三区四区乱码| 亚洲在久久综合| 久久久久人妻精品一区果冻| 久久狼人影院| 极品人妻少妇av视频| 人妻少妇偷人精品九色| 女人久久www免费人成看片| 国产免费视频播放在线视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 女的被弄到高潮叫床怎么办| 熟女电影av网| 国产精品久久久久久av不卡| 亚洲性久久影院| 亚洲欧美清纯卡通| 亚洲欧美成人精品一区二区| 亚洲怡红院男人天堂| 久久99热6这里只有精品| 免费高清在线观看视频在线观看| 欧美精品亚洲一区二区| 永久网站在线| 一本大道久久a久久精品| 99视频精品全部免费 在线| 亚洲情色 制服丝袜| 国产精品人妻久久久久久| 极品少妇高潮喷水抽搐| .国产精品久久| 日本爱情动作片www.在线观看| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| 久久国产精品男人的天堂亚洲 | 黑人欧美特级aaaaaa片| 免费看不卡的av| 男女边吃奶边做爰视频| 久久精品熟女亚洲av麻豆精品| 亚洲综合色网址| 亚洲少妇的诱惑av| 精品久久久久久电影网| 在线 av 中文字幕| 成年美女黄网站色视频大全免费 | 亚洲国产精品一区二区三区在线| 热re99久久国产66热| av在线老鸭窝| 国产成人91sexporn| 精品人妻一区二区三区麻豆| 亚洲国产精品专区欧美| 日本黄色片子视频| 亚洲av中文av极速乱| 欧美xxxx性猛交bbbb| 久久免费观看电影| 国产精品久久久久成人av| 欧美变态另类bdsm刘玥| 亚洲综合色惰| 丝瓜视频免费看黄片| 观看av在线不卡| 免费黄色在线免费观看| 欧美xxⅹ黑人| 免费日韩欧美在线观看| 人体艺术视频欧美日本| 大片免费播放器 马上看| 99国产综合亚洲精品| 能在线免费看毛片的网站| 亚洲成人av在线免费| 卡戴珊不雅视频在线播放| 国产精品99久久久久久久久| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 亚洲综合色网址| 日韩视频在线欧美| 色哟哟·www| 亚洲丝袜综合中文字幕| 高清av免费在线| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频 | 女人精品久久久久毛片| 亚洲av在线观看美女高潮| 99热国产这里只有精品6| 免费黄色在线免费观看| 中文天堂在线官网| 国产老妇伦熟女老妇高清| 国产av码专区亚洲av| 18禁在线无遮挡免费观看视频| 成年人免费黄色播放视频| av天堂久久9| 日韩精品免费视频一区二区三区 | 69精品国产乱码久久久| 午夜激情福利司机影院| 人成视频在线观看免费观看| 老司机亚洲免费影院| 高清黄色对白视频在线免费看| 亚洲国产精品成人久久小说| 69精品国产乱码久久久| 精品一区二区三卡| 一区二区三区四区激情视频| 国产女主播在线喷水免费视频网站| 成人二区视频| 777米奇影视久久| 高清视频免费观看一区二区| 中文字幕精品免费在线观看视频 | 在线观看美女被高潮喷水网站| 免费高清在线观看日韩| 麻豆乱淫一区二区| 女人精品久久久久毛片| 亚洲无线观看免费| 久久国内精品自在自线图片| 91aial.com中文字幕在线观看| videos熟女内射| 中文字幕最新亚洲高清| 人人妻人人澡人人爽人人夜夜| 人妻少妇偷人精品九色| 成人国产麻豆网| 成人午夜精彩视频在线观看| 啦啦啦视频在线资源免费观看| 亚洲欧美精品自产自拍| 国产无遮挡羞羞视频在线观看| 97超视频在线观看视频| 精品少妇黑人巨大在线播放| 一本一本综合久久| 哪个播放器可以免费观看大片| 国产高清三级在线| 午夜免费鲁丝| 成人毛片60女人毛片免费| 欧美 日韩 精品 国产| 国产亚洲欧美精品永久| 五月开心婷婷网| 最近2019中文字幕mv第一页| 日本wwww免费看| av专区在线播放| 午夜福利视频精品| 国产黄色免费在线视频| 自线自在国产av| 国产视频首页在线观看| 最近中文字幕2019免费版| 久久99热6这里只有精品| 成人国语在线视频| 伊人久久精品亚洲午夜| 国产精品久久久久久精品电影小说| 精品久久久久久久久亚洲| av专区在线播放| 好男人视频免费观看在线| 午夜日本视频在线| 九九爱精品视频在线观看| 亚州av有码| 亚洲在久久综合| 亚洲欧美中文字幕日韩二区| 亚洲国产av影院在线观看| 91精品一卡2卡3卡4卡| 高清欧美精品videossex| 日本午夜av视频| 一区二区日韩欧美中文字幕 | 一本大道久久a久久精品| 欧美日韩视频精品一区| 日韩精品免费视频一区二区三区 | 九色亚洲精品在线播放| 日本vs欧美在线观看视频| 久久亚洲国产成人精品v| 大香蕉97超碰在线| 国产一区二区在线观看日韩| 少妇被粗大的猛进出69影院 | 国产精品国产三级国产专区5o| 爱豆传媒免费全集在线观看| 欧美激情极品国产一区二区三区 | 国产午夜精品一二区理论片| av电影中文网址| 日本午夜av视频| 少妇的逼水好多| 一级毛片我不卡| 日本av免费视频播放| 青春草视频在线免费观看| 亚洲成人手机| 国产午夜精品久久久久久一区二区三区| 免费少妇av软件| 亚洲成人手机| 欧美精品一区二区免费开放| 成人综合一区亚洲| 亚洲婷婷狠狠爱综合网| 大香蕉久久成人网| 免费少妇av软件| 亚洲成人手机| 欧美精品一区二区免费开放| 久久婷婷青草| 国产成人a∨麻豆精品| 超色免费av| 久久狼人影院| 欧美日韩国产mv在线观看视频| 午夜av观看不卡| 80岁老熟妇乱子伦牲交| 亚洲综合精品二区| 国产午夜精品久久久久久一区二区三区| 国产色婷婷99| 99九九在线精品视频| 国产熟女午夜一区二区三区 |