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

    牛糞堆肥中氨氧化細(xì)菌群落結(jié)構(gòu)及其與環(huán)境因子相關(guān)性研究

    2017-02-14 01:41:28孫雪微許修宏孟慶欣成利軍張文浩門夢(mèng)琪許本姝
    關(guān)鍵詞:銨態(tài)氮牛糞群落

    孫雪微,許修宏,孟慶欣,成利軍,張文浩,門夢(mèng)琪,許本姝,孫 瑜

    (東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,哈爾濱 150030)

    牛糞堆肥中氨氧化細(xì)菌群落結(jié)構(gòu)及其與環(huán)境因子相關(guān)性研究

    孫雪微,許修宏*,孟慶欣,成利軍,張文浩,門夢(mèng)琪,許本姝,孫 瑜

    (東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,哈爾濱 150030)

    采用PCR-DGGE(聚合酶鏈反應(yīng)-變性梯度凝膠電泳)和Real-time PCR技術(shù),分析牛糞堆肥期間氨氧化細(xì)菌群落結(jié)構(gòu)和數(shù)量的變化及其與理化因子的相關(guān)性。結(jié)果表明,在牛糞堆肥過程中氨氧化細(xì)菌數(shù)量變化范圍為1.62×106~3.80×107copies·g-1,在第14 d氨氧化細(xì)菌數(shù)量達(dá)到最高值;DGGE圖譜顯示,堆肥各個(gè)時(shí)期氨氧化細(xì)菌群落組成存在明顯的差異;Shannon指數(shù)隨著溫度的變化而變化,在第7 d達(dá)到最大值2.671 8;系統(tǒng)發(fā)育分析表明,氨氧化細(xì)菌均隸屬于β-變形菌綱不可培養(yǎng)的亞硝化螺菌屬Uncultured Nitrosospira sp.和不可培養(yǎng)的亞硝化單胞菌屬Uncultured Nitrosomonas sp.,其中與Uncultured Nitrosospira sp.相似性高的氨氧化細(xì)菌為優(yōu)勢(shì)菌,占60%左右;冗余分析(Redundancy analysis,RDA)顯示,溫度和銨態(tài)氮等對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)的影響較大。研究表明在牛糞堆肥期間氨氧化細(xì)菌群落結(jié)構(gòu)發(fā)生了明顯的變化,同時(shí)證實(shí)了氨氧化細(xì)菌群落結(jié)構(gòu)的變化受堆肥理化因子的影響。

    牛糞堆肥;PCR-DGGE;Real-time PCR;氨氧化細(xì)菌;冗余分析

    隨著我國(guó)畜禽養(yǎng)殖業(yè)的大規(guī)模發(fā)展,畜禽糞便排放量也在逐年遞增,給環(huán)境帶了巨大的壓力[1]。好氧堆肥是一種對(duì)牛糞等有機(jī)固體廢棄物進(jìn)行資源化、無害化和穩(wěn)定化的經(jīng)濟(jì)高效的處理方法[2]。畜禽糞便經(jīng)堆肥處理后轉(zhuǎn)變成有機(jī)肥料,實(shí)現(xiàn)了資源再利用,但在堆肥過程中,往往伴隨著嚴(yán)重的氮素?fù)p失,一般損失量為初始總氮量的17%~76%[3],既降低了堆肥產(chǎn)品的質(zhì)量,也會(huì)釋放臭氣和氮氧化物等氣體,對(duì)環(huán)境造成二次污染[4]。因此,控制堆肥過程氨氧化反應(yīng)造成的氮素?fù)p失,已成為提高堆肥質(zhì)量的重要途徑之一。

    硝化作用作為全球氮循環(huán)最重要的環(huán)節(jié)之一,以亞硝酸鹽作為載體將氨鹽氧化成硝酸鹽[5]。由氨氧化細(xì)菌催化的氨氧化反應(yīng)是硝化作用的限速步驟,在全球氮素生物化學(xué)循環(huán)中起著重要的作用[6]。氨氧化細(xì)菌屬于專性化能自養(yǎng)的功能性微生物[7],其驅(qū)動(dòng)的氨氧化反應(yīng)是通過氨單加氧酶(Ammonia monooxygenase,AMO)催化,編碼AMO基因來完成的[8]。隨著現(xiàn)代分子生物學(xué)技術(shù)的發(fā)展,對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)的研究也引起了國(guó)內(nèi)外研究者的關(guān)注,如Yamada等[9]在氨含量較高的牛糞堆肥中檢測(cè)到了氨氧化細(xì)菌群落的動(dòng)態(tài)變化,邱珊蓮等[10]在雞糞堆肥中發(fā)現(xiàn)氨氧化細(xì)菌占據(jù)優(yōu)勢(shì)地位。但在堆肥領(lǐng)域中,氨氧化細(xì)菌的理論知識(shí)仍未完善,無法全面地對(duì)生產(chǎn)實(shí)踐起到理論指導(dǎo)的作用,所以,將氨氧化細(xì)菌作為研究對(duì)象,研究其群落結(jié)構(gòu)的動(dòng)態(tài)演替及其與環(huán)境因子之間的相關(guān)性,對(duì)于揭示堆肥過程中氮素轉(zhuǎn)移與變化規(guī)律具有深遠(yuǎn)的意義。

    本研究采用變性梯度凝膠電泳(Denaturing gra dient gel electrophoresis,DGGE)和Real-time PCR技術(shù)檢測(cè)牛糞堆肥過程中氨氧化細(xì)菌群落結(jié)構(gòu)及數(shù)量的動(dòng)態(tài)變化,同時(shí)采用冗余分析(Redundancy analy sis,RDA)方法分析氨氧化細(xì)菌群落結(jié)構(gòu)與堆肥環(huán)境因子之間的相關(guān)性,來了解堆肥過程中氮素的轉(zhuǎn)移與轉(zhuǎn)化規(guī)律,為減少氮素?fù)p失提供理論依據(jù)。

    1 材料與方法

    1.1 堆肥及樣品采集

    堆肥材料由牛糞和水稻秸稈組成(表1)。牛糞源于哈爾濱幸福鄉(xiāng)奶牛養(yǎng)殖場(chǎng),水稻秸稈源于哈爾濱市香坊區(qū)東北農(nóng)業(yè)大學(xué)實(shí)驗(yàn)基地。將水稻秸稈截成5 cm左右,與牛糞按質(zhì)量1∶3.5的比例均勻混合,并調(diào)整水分含量為65%~70%,置于規(guī)格為50 cm×50 cm× 100 cm(長(zhǎng)×寬×高)的堆肥模擬容器中,每隔30 min通風(fēng)一次,采用上、中、下三點(diǎn)取樣方式,均勻混合,每次取樣后翻堆。堆肥過程共36 d,取樣時(shí)間為:原始樣品第1 d,升溫階段第2 d,高溫階段第3、4、7、11、14 d,腐熟階段第19、25、36 d,樣品編號(hào)分別為D1、D2、D3、D4、D7、D11、D14、D19、D25和D36。堆肥樣品分為兩份:一份在-4℃儲(chǔ)存,用于理化指標(biāo)的測(cè)定;另一份在-80℃儲(chǔ)存,用于DNA提取。

    表1 堆肥材料的主要成分Table 1 Properties of raw materials for composting

    1.2 理化指標(biāo)測(cè)定

    用LNI-TUT325數(shù)字電子溫度計(jì)測(cè)定堆體上、中、下三個(gè)位置的溫度。堆肥樣品的pH采用鮑士旦[11]的農(nóng)化分析方法進(jìn)行測(cè)定。水溶性銨態(tài)氮和硝態(tài)氮以2 mol·L-1氯化鉀溶液浸提,采用連續(xù)流動(dòng)分析儀(AA3,德國(guó))測(cè)定其含量。堆肥潛在硝化能力(PAO)的測(cè)定采用He等[12]的研究方法進(jìn)行測(cè)定。稱取10 g堆肥樣品,置于100 mL 25 mg N·L-1NH4Cl緩沖液中,25℃,150 r·min-1恒溫?fù)u床振蕩24 h,每隔4 h取培養(yǎng)液一次,用連續(xù)流動(dòng)分析儀測(cè)定NO2-N的含量。全碳和全氮分別采用重鉻酸鉀容量法和凱氏定氮法[11]測(cè)定。

    1.3 DNA提取和PCR擴(kuò)增

    取堆肥樣品3~5 g,采用Liu等[13]的方法提取堆肥樣品的總DNA,用1%的瓊脂糖凝膠電泳檢測(cè)后,用DNA純化試劑盒(Omega,美國(guó))純化回收DNA,再用1%的瓊脂糖凝膠電泳檢測(cè)結(jié)果,于-20℃的冰箱保存?zhèn)溆谩?/p>

    PCR擴(kuò)增采用氨氧化細(xì)菌的16S rRNA特異性引物進(jìn)行巢式PCR(Nested PCR)擴(kuò)增,引物序列和反應(yīng)程序見表2。第一輪,反應(yīng)體系為25 μL,包含400 μmol·L-1dNTP 2.5 μL,10×PCR Buffer(含MgCl2)2.5 μL,引物對(duì)各0.2 μL(上游引物為CTO189fAB與CTO189fC按照 2∶1的比例混合,下游引物為CTO654r)[14],Taq聚合酶0.3 μL,模板DNA 0.5 μL,添加無菌水調(diào)整終體積到25 μL。PCR擴(kuò)增產(chǎn)物大小為465 bp,將其稀釋100倍作為第二輪PCR反應(yīng)的DNA模板。第二輪PCR擴(kuò)增反應(yīng)條件同上,引物選用細(xì)菌特異性引物338f-GC和518r[15],PCR擴(kuò)增產(chǎn)物大小為234 bp。最后用1%的瓊脂糖凝膠電泳檢測(cè)PCR擴(kuò)增產(chǎn)物。

    表2 PCR引物及反應(yīng)條件Table 2 Primer and reaction condition of PCR amplification

    1.4 DGGE

    使用Bio-Rad DCodeTMsystem(Bio-Rad,美國(guó))對(duì)氨氧化細(xì)菌的16S rRNA PCR擴(kuò)增片段進(jìn)行DGGE分析,最適的變性梯度為30%~60%(100%的變性劑量中含有40%的去離子甲酰胺和7 mol·L-1尿素),聚丙烯酰胺膠的濃度為8%[丙烯酰胺∶雙丙烯酰胺= 37.5∶1(m·V-1)]。電泳運(yùn)行條件:在1×TAE電泳緩沖液中,恒溫60℃,恒壓120 V,電泳12 h。電泳完成后,采用硝酸銀銀染法染色,拍照后將回收的優(yōu)勢(shì)DGGE條帶置于滅菌的1.5 mL離心管中并編號(hào),加30 μL無菌水,-4℃保存過夜,采用引物338f和518r進(jìn)行PCR擴(kuò)增后,將擴(kuò)增產(chǎn)物送測(cè)序公司完成測(cè)序。

    1.5 Real-time PCR

    Real-time PCR反應(yīng)采用 SYBY Green法,在CFX-96TMReal-Time PCR儀上進(jìn)行,采用的氨氧化細(xì)菌引物對(duì)為CTO189fAB、CTO189fC和CTO654r,反應(yīng)體系為20 μL,包含10 μL SYBR Premix Ex TaqⅡ、1.6 μL稀釋10倍的DNA模板以及上下游引物各0.8 μL,其余以6.8 μL無菌水補(bǔ)齊。Real-time反應(yīng)程序?yàn)椋?5℃30 s;95℃5 s,60℃30 s共40個(gè)循環(huán)。RealtimePCR標(biāo)準(zhǔn)曲線通過特異性引物PCR擴(kuò)增DNA,擴(kuò)增產(chǎn)物純化后,連接到PMD-18T載體上,再轉(zhuǎn)化到大腸桿菌DH5α感受態(tài)細(xì)胞中,通過藍(lán)白斑試驗(yàn),利用引物M13F(5′-CGCCAGGGTTTTCCCAGTCACGAC-3′)和M13R(5′-AGCGGATAACAATTTCACACAGGA-3′)篩選陽(yáng)性克隆,提取質(zhì)粒后,采用A260/A280核酸蛋白測(cè)定儀測(cè)定其濃度,以8倍稀釋梯度獲取氨氧化細(xì)菌的標(biāo)準(zhǔn)曲線。標(biāo)準(zhǔn)曲線的擴(kuò)增效率和相關(guān)系數(shù)分別為99.9%和0.94。

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

    采用Quantity One軟件(Bio-Rad,美國(guó))對(duì)DGGE圖譜進(jìn)行分析,采用UPGAMA法進(jìn)行聚類分析。采用多樣性指數(shù)(H),均勻度指數(shù)(E)和豐富度(S)來分析氨氧化細(xì)菌16S rRNA基因的多樣性,計(jì)算公式如下:

    式中:H為Shannon多樣性指數(shù);E為均勻度指數(shù);S為DGGE膠中條帶數(shù)目,即豐富度;Pi為第i條帶灰度占該樣品總灰度的比率[16]。

    通過Blast程序?qū)⑺脙?yōu)勢(shì)條帶的序列與已知序列進(jìn)行比對(duì),采用Mega 6.0軟件構(gòu)建氨氧化細(xì)菌16S rRNA基因序列的系統(tǒng)發(fā)育樹;采用Canoco 4.5軟件對(duì)堆肥樣品中氨氧化細(xì)菌群落結(jié)構(gòu)與環(huán)境因子進(jìn)行冗余分析。

    2 結(jié)果與討論

    2.1 堆肥過程中理化指標(biāo)的變化

    理化指標(biāo)的測(cè)定結(jié)果見表3。隨著堆肥的發(fā)酵,微生物快速分解有機(jī)質(zhì),釋放大量的熱量,使堆體溫度急劇上升,在第7 d堆體溫度達(dá)到最高溫度,并且溫度超過55℃的堆肥時(shí)間持續(xù)了11 d,殺死大量的病原菌,達(dá)到了基本堆肥的衛(wèi)生標(biāo)準(zhǔn)[17],堆體溫度在第14 d開始持續(xù)下降,且堆肥后期溫度接近室溫。在堆肥過程中,pH變化范圍為8.85~8.14,在第4 d達(dá)到最大值,一直處于弱堿性的狀態(tài)。在堆肥前期,含氮有機(jī)物快速降解,銨態(tài)氮大量積累,從而使堆體的pH值不斷上升,而在堆肥后期,銨態(tài)氮減少,堆體pH值開始下降[18]。C/N是評(píng)價(jià)堆肥腐熟度的一種方法,在整個(gè)堆肥過程中C/N持續(xù)下降,在堆肥結(jié)束時(shí)下降到19.6,據(jù)此可以認(rèn)為堆肥已達(dá)到腐熟標(biāo)準(zhǔn)[19]。

    在堆肥過程中,銨態(tài)氮含量呈現(xiàn)先升后降的趨勢(shì),因?yàn)樵诙逊食跗?,微生物氨化作用及礦化作用加速了氮素的分解,所以堆肥前期銨態(tài)氮含量不斷增多;而堆肥后期,由于部分銨態(tài)氮轉(zhuǎn)化為硝態(tài)氮和氨氣揮發(fā)等原因,導(dǎo)致銨態(tài)氮含量下降[20]。硝態(tài)氮含量呈現(xiàn)波動(dòng)性上升的趨勢(shì),含量增加速度緩慢,可能是由于過多的銨態(tài)氮和過高的溫度,在一定程度上抑制了微生物的活性,使硝態(tài)氮生成速度減緩[21]。PAO的變化范圍是33.0~94.5 ng N·g-1·h-1,在第7 d達(dá)到最大值。

    表3 堆肥過程中理化性質(zhì)的變化Table 3 Changes in physic-chemical properties in composting

    2.2 樣品總DNA的提取和PCR的擴(kuò)增

    采用Liu等[13]的研究方法成功地提取了所用堆肥樣品的總DNA,其長(zhǎng)度約23 kb。將純化后的DNA作為進(jìn)行巢式PCR第一輪的模板,所有樣品均獲得了特異性擴(kuò)增產(chǎn)物,其大小為465 bp左右,將其稀釋100倍作為第二輪PCR的DNA模板,均得到了234 bp左右的特異性片段。

    2.3 氨氧化細(xì)菌豐度的分析

    采用Real-time PCR技術(shù)對(duì)堆肥樣品進(jìn)行氨氧化細(xì)菌16S rRNA基因的定量分析,獲得的氨氧化細(xì)菌基因豐度如圖1所示。在整個(gè)堆肥過程中,氨氧化細(xì)菌基因豐度較低,基因拷貝數(shù)的變化范圍為1.62× 106~3.8×107copies·g-1,在堆肥前期氨氧化細(xì)菌的基因拷貝數(shù)持續(xù)下降,第14 d達(dá)到最大值。氨氧化細(xì)菌基因豐度基本上都處于同一數(shù)量級(jí),波動(dòng)范圍并不是特別大。這與Lu等[22]的研究結(jié)果一致,在淡水水產(chǎn)養(yǎng)殖池塘中氨氧化細(xì)菌的基因豐度較低。隨著堆肥進(jìn)入發(fā)酵高溫期,氨氧化細(xì)菌基因拷貝數(shù)明顯降低。這可能與溫度影響有關(guān),溫度過高,抑制了多數(shù)微生物的生長(zhǎng)繁殖,僅一些耐熱微生物在快速生長(zhǎng),所以高溫期的氨氧化細(xì)菌相對(duì)較少[23],而在堆肥第14 d,堆體溫度達(dá)到55℃,隨后進(jìn)入降溫階段,有利于微生物快速生長(zhǎng),所以氨氧化細(xì)菌的基因拷貝數(shù)明顯增多。在本研究中高銨離子濃度可能是導(dǎo)致氨氧化細(xì)菌豐度低的另一個(gè)關(guān)鍵原因,隨著銨態(tài)氮的逐漸積累,在一定程度上抑制了氨氧化細(xì)菌的生長(zhǎng),在堆肥后期受銨態(tài)氮濃度的影響,氨氧化細(xì)菌的生長(zhǎng)減緩,導(dǎo)致基因豐度下降[24]。這與Wuchter等[25]的研究結(jié)果一致,在海水中氨氧化微生物的基因豐度與銨態(tài)氮的濃度顯著負(fù)相關(guān)。

    圖1 堆肥過程中氨氧化細(xì)菌的基因拷貝數(shù)Figure 1 amoA gene copy numbers per gram compost sample during composting

    2.4 氨氧化細(xì)菌16S rRNA基因群落多樣性的分析

    采用氨氧化細(xì)菌16S rRNA基因特異性引物對(duì)擴(kuò)增產(chǎn)物進(jìn)行DGGE,得到的氨氧化細(xì)菌DGGE圖譜如圖2所示。圖譜中各個(gè)泳道的條帶有較大的差異,說明堆肥各時(shí)期氨氧化細(xì)菌的群落變化較劇烈。利用Quantity One軟件進(jìn)行數(shù)字化處理,聚類分析(圖3)結(jié)果也顯示了氨氧化細(xì)菌在堆肥過程中群落結(jié)構(gòu)存在明顯的差異,泳道3和8、6和10、5和7、4和9彼此聚合,相似性系數(shù)分別達(dá)49%、31%、46%和59%。這說明各泳道群落組成相似性較高,而泳道1和2并未與其他泳道聚為一類,表明進(jìn)入高溫期后氨氧化細(xì)菌的菌群組成發(fā)生了劇烈的變化。

    圖2 牛糞堆肥中氨氧化細(xì)菌DGGE圖譜Figure 2 DGGE profile of ammonia oxidizing bacteria in cow manure composting

    圖3 氨氧化細(xì)菌DGGE圖譜聚類分析Figure 3 Cluster analysis of DGGE profile of ammonia oxidizing bacteria in cow manure composting

    通過Quantity One軟件計(jì)算出氨氧化細(xì)菌群落的多樣性指數(shù)、均勻度指數(shù)和豐富度指數(shù)(表4),結(jié)果顯示,均勻度指數(shù)變化不大,但多樣性指數(shù)和豐富度指數(shù)出現(xiàn)了明顯的改變。隨著堆肥的進(jìn)行,氨氧化細(xì)菌的多樣性指數(shù)呈現(xiàn)先升后降的趨勢(shì),在第7 d達(dá)到最大值2.671 8。這可能是因?yàn)樵诙逊是捌?,銨態(tài)氮等硝化作用底物含量充足,營(yíng)養(yǎng)物質(zhì)豐富,利于氨氧化細(xì)菌的生長(zhǎng),導(dǎo)致多樣性指數(shù)上升;而隨著堆肥的不斷發(fā)酵、溫度的升高和銨態(tài)氮濃度的下降,部分氨氧化細(xì)菌的生長(zhǎng)受到了影響而被淘汰,使多樣性指數(shù)逐漸減小[26]。

    2.5 系統(tǒng)發(fā)育分析

    氨氧化細(xì)菌DGGE圖譜優(yōu)勢(shì)條帶切膠測(cè)序比對(duì)結(jié)果見表5。序列比對(duì)結(jié)果表明,大多數(shù)序列與不可培養(yǎng)的亞硝化螺菌屬Uncultured Nitrosospira sp.和不可培養(yǎng)的亞硝化單胞菌屬Uncultured Nitrosomonas sp.這兩個(gè)屬的細(xì)菌相似度較高,在90%~100%之間,其中,與Uncultured Nitrosospira sp.相似度高的氨氧化細(xì)菌為優(yōu)勢(shì)菌,約占60%,而Uncultured Nitrosomonas sp.菌群約占40%。由所獲得優(yōu)勢(shì)菌序列構(gòu)建的氨氧化細(xì)菌系統(tǒng)發(fā)育樹見圖4。由圖可知,氨氧化細(xì)菌序列被劃分成Uncultured Nitrosospira sp.和Uncultured Nitrosomonas sp.兩大分支,其中Uncultured Nitrosospira sp.分支占優(yōu)勢(shì)。

    表4 堆肥過程中氨氧化細(xì)菌Shannon指數(shù)、均勻度指數(shù)和豐富度指數(shù)的變化Table 4 Changes in Shannon,evenness and richness index of ammonia oxidizing bacteria in composting

    表5 DGGE優(yōu)勢(shì)條帶序列比對(duì)結(jié)果Table 5 Results of BLAST analysis on the sequences of dominant DGGE bands

    通過氨氧化細(xì)菌16S rRNA基因系統(tǒng)發(fā)育樹分析可知,在牛糞堆肥過程中氨氧化細(xì)菌都隸屬于β-變形菌綱,多數(shù)序列與Uncultured Nitrosospira sp.和Uncultured Nitrosomonas sp.的相似度較高。這與虞泳等[27]研究結(jié)果一致,在農(nóng)業(yè)好氧堆肥中也發(fā)現(xiàn)了Nitrosomonas sp.和Nitrosospira sp.優(yōu)勢(shì)菌群。由此推測(cè),在堆肥過程中氨氧化細(xì)菌主要是Nitrosomonas sp.和Nitrosospirasp.菌。在本研究中,以Nitrosospira sp.屬的氨氧化細(xì)菌為優(yōu)勢(shì)菌群,可能是由于高濃度的銨態(tài)氮含量所引起的。Enwall等[28]的研究結(jié)果表明,在銨態(tài)氮濃度較高的條件下,Nitrosospira sp.菌生長(zhǎng)更為旺盛。

    圖4 牛糞堆肥中氨氧化細(xì)菌16S rRNA基因的系統(tǒng)發(fā)育樹Figure 4 Phylogenetic tree of ammonia oxidizing bacteria based on 16S rRNA gene in cow manure composting

    2.6 氨氧化細(xì)菌群落結(jié)構(gòu)與堆肥環(huán)境因子的冗余分析

    采用冗余分析方法分析堆肥過程中氨氧化細(xì)菌群落結(jié)構(gòu)變化與理化因子的相關(guān)性,結(jié)果見圖5。二維排序軸的前兩個(gè)排序軸分別解釋了種群與環(huán)境的積累變化率為79.2%和24.4%。在二維排序圖中,箭頭連線的長(zhǎng)短表示群落與該環(huán)境因子之間相關(guān)系數(shù)的大小,箭頭連線與排序軸的夾角表示該環(huán)境因子與排序軸相關(guān)性的大小,夾角越大,其相關(guān)性越小[29]。硝態(tài)氮和PAO箭頭連線較長(zhǎng),銨態(tài)氮、C/N、溫度和pH的箭頭連線較短,所以要據(jù)箭頭連線的長(zhǎng)度可知,環(huán)境因子對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)影響的相關(guān)性大小順序?yàn)镻AO>硝態(tài)氮>銨態(tài)氮>C/N>溫度>pH。從圖中也可以看到,隨著堆肥的進(jìn)行,氨氧化細(xì)菌的群落結(jié)構(gòu)變化顯著,在升溫階段(day2和day3),氨氧化細(xì)菌群落結(jié)構(gòu)較為相似,而在高溫階段與降溫階段,氨氧化細(xì)菌群落結(jié)構(gòu)劇烈變化。

    通過冗余分析可知,溫度、PAO、銨態(tài)氮、硝態(tài)氮和pH對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)的影響較大,基于手動(dòng)選擇的冗余分析表明溫度和銨態(tài)氮的影響較為顯著(P<0.05),這與Zhang等[30]研究結(jié)果相似,即環(huán)境因子顯著影響著氨氧化細(xì)菌的群落結(jié)構(gòu)變化。溫度能夠影響氨氧化細(xì)菌的生長(zhǎng)速度、酶反應(yīng)動(dòng)力學(xué)以及化合物的溶解度等[31],會(huì)對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)產(chǎn)生明顯的影響。在本研究的堆肥體系中,隨著溫度的升高,氨氧化細(xì)菌的群落結(jié)構(gòu)發(fā)生了明顯的變化,表明氨氧化細(xì)菌對(duì)溫度的適應(yīng)情況各不相同,僅部分耐高溫氨氧化細(xì)菌存活。這與Yamamoto等[32]的研究結(jié)果一致,即在堆肥高溫期氨氧化細(xì)菌是唯一的耐熱硝化細(xì)菌。

    圖5 氨氧化細(xì)菌群落結(jié)構(gòu)的冗余分析(RDA)Figure 5 Redundancy analysis(RDA)of ammonia oxidizing bacterial community structure during composting

    銨態(tài)氮作為氨氧化反應(yīng)的直接底物,其濃度與氨氧化細(xì)菌的數(shù)量與種類有直接關(guān)系[33]。宋亞娜等[34]的研究指出,提高銨濃度可以促進(jìn)間作土氨氧化細(xì)菌的生長(zhǎng),對(duì)氨氧化細(xì)菌的群落結(jié)構(gòu)和多樣性產(chǎn)生影響。在本研究中,堆肥前期銨態(tài)氮不斷積累,營(yíng)養(yǎng)物質(zhì)豐富,氨氧化細(xì)菌生長(zhǎng)旺盛,而在堆肥后期,由于銨態(tài)氮濃度逐漸下降,氨氧化細(xì)菌因底物濃度減小而衰退,表明銨態(tài)氮對(duì)氨氧化細(xì)菌的群落結(jié)構(gòu)的變化有十分顯著的影響。

    3 結(jié)論

    (1)堆肥樣品氨氧化細(xì)菌16S rRNA基因豐度的變化范圍是1.62×106~3.80×107copies·g-1,在第14 d氨氧化細(xì)菌的基因拷貝數(shù)達(dá)到最大值。

    (2)不同時(shí)期堆肥樣品氨氧化細(xì)菌多樣性指數(shù)存在明顯的差異,隨溫度達(dá)到最高值,氨氧化細(xì)菌的多樣性指數(shù)也達(dá)到最大。

    (3)通過序列比對(duì)和系統(tǒng)發(fā)育分析可知,氨氧化細(xì)菌均隸屬于β-變形菌綱,多數(shù)序列與Uncultured Nitrosospira sp.和Uncultured Nitrosomonas sp.屬的細(xì)菌相似度較高,與Uncultured Nitrosospira sp.屬相似度高的氨氧化細(xì)菌為優(yōu)勢(shì)菌。

    (4)冗余分析表明,堆肥環(huán)境因子對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)的影響顯著,其相關(guān)性大小為PAO>硝態(tài)氮>銨態(tài)氮>C/N>溫度>pH,其中溫度和銨態(tài)氮對(duì)氨氧化細(xì)菌群落結(jié)構(gòu)的影響更為顯著。

    [1]朱兆良,David N,孫 波.中國(guó)農(nóng)業(yè)面源污染控制對(duì)策[M].北京:中國(guó)環(huán)境科學(xué)出版社,2006.

    ZHU Zhao-liang,David N,SUN Bo.Agricultural non-point source pollution control measures in China[M].Beijing:China Environmental Science Press,2006.

    [2]Wei Y S,Fan Y B,Wang M J,et al.Composting and compost application in China[J].Resources Conservation&Recycling,2000,30(4):277-300.

    [3]Barrington S C D,Trigui M,Knight W.Effect of carbon source on compost nitrogen and carbon losses[J].Bioresource Technology,2002,83 (3):189-194.

    [4]殷 震,劉景華.動(dòng)物病毒學(xué)[M].二版.北京:科學(xué)出版社,1997.

    YIN Zhen,LIU Jing-hua.Animal virology[M].2nd Edition.Beijing:Science Press,1997.

    [5]K?nneke M,Bernhard A E,Torre J R D L,et al.Isolation of an autotrophic ammonia-oxidizing marine archaeon[J].Nature,2005,437 (7058):543-546.

    [6]Lullo L D,Addesse R,Comegna C,et al.Influence of effluent irrigation on community composition and function of ammonia-oxidizing bacteria in soil[J].Applied&Environmental Microbiology,2001,67(8):3426-3433.

    [7]Teske A,Alm E,Regan J M,et al.Evolutionary relationships among ammonia-and nitrite-oxidizing bacteria[J].Journal of Bacteriology,1994, 176(21):6623-6630.

    [8]Zeng J,Zhao D Y,Huang R,et al.Abundance and community composition of ammonia-oxidizing archaea and bacteria in two different zones of Lake Taihu[J].Canadian Journal of Microbiology,2012,58(8):1018-1026.

    [9]Yamada T,Araki S,Ikeda-Ohtsubo W,et al.Community structure and population dynamics of ammonia oxidizers in composting processes of ammonia-rich livestock waste[J].Systematic&Applied Microbiology, 2013,36(5):359-367.

    [10]邱珊蓮,張少平,翁伯琦,等.雞糞菌渣好氧堆肥過程中氨氧化古菌及氨氧化細(xì)菌群落的動(dòng)態(tài)變化[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2016,35 (3):573-583.

    QIU Shan-lian,ZHANG Shao-ping,WENG Bo-qi,et al.Dynamics of ammonia-oxidizing archaea and ammonia-oxidizing bacteria during composting of chicken manure and mushroom cultural waste[J].Journal of Agro-Environment Science,2016,35(3):573-583.

    [11]鮑士旦.土壤農(nóng)化分析[M].三版.北京:中國(guó)農(nóng)業(yè)出版社,2008:39-89.

    BAO Shi-dan.Soil and agriculture chemistry analysis[M].3rd Edition. Beijing:China Agriculture Press,2008:39-89.

    [12]He J Z,Shen J P,Zhang L M,et al.Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammoniaoxidizing archaea of a Chinese upland red soil under long-term fertilization practices[J].Environmental Microbiology,2007,9(9):2364-2374.

    [13]Liu J,Xu X H,Li H T,et al.Effect of microbiological inocula on chemical and physical properties and microbial community of cow manure compost[J].Biomass&Bioenergy,2011,35(8):3433-3439.

    [14]Innerebner G,Knapp B,Vasara T,et al.Traceability of ammonia-oxidizing bacteria in compost-treated soils[J].Soil Biology&Biochemistry,2006,38(5):1092-1100.

    [15]Li A J,Yang S F,Li X Y,et al.Microbial population dynamics during aerobic sludge granulation at different organic loading rates[J].Water Research,2008,42(13):3552-3560.

    [16]Yang Q,Zhang J Z,Wang S,et al.Bacterial community structure and diversity in a black soil as affected by long-term fertilization[J].Pedosphere,2008,18(5):582-592.

    [17]唐 珠.豬糞堆肥的理化特征及腐熟度評(píng)價(jià)研究[D].南京農(nóng)業(yè)大學(xué),2011.

    TANG Zhu.Characterization and identification of maturity of pig manure composts[D].Nanjing Agricultural University,2011.

    [18]張亞寧.堆肥腐熟度快速測(cè)定指標(biāo)和方法的建立[D].北京:中國(guó)農(nóng)業(yè)大學(xué),2004.

    ZHANG Ya-ning.Studies on the foundation of simple indexes and mensuration of compost maturity[D].Beijing:China Agriculture University,2004.

    [19]廖新俤,吳銀寶,汪植三,等.堆體大小對(duì)豬糞堆肥影響和袋裝堆肥的研究[J].農(nóng)業(yè)工程學(xué)報(bào),2003,19(4):287-290.

    LIAO Xin-di,WU Yin-bao,WANG Zhi-san,et al.Effects of pile size and packaging on swine manure composting[J].Transactions of the Chinese Society of Agricultural Engineering,2003,19(4):287-290.

    [20]徐大勇,黃為一.人工接種堆肥和自然堆肥微生物區(qū)系與分子多態(tài)性的變化[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2006,22(1):29-33.

    XU Da-yong,HUANG Wei-yi.Variation of microflora and molecular polymorphism in inoculated and natural composts[J].Journal of Ecology and Rural Environment,2006,22(1):29-33.

    [21]李海波,廉 靜,李 敏,等.反硝化菌株GW1的篩選及特性研究[J].河北科技大學(xué)學(xué)報(bào),2012,33(2):184-189.

    LI Hai-bo,LIAN Jing,LI Min,et al.Isolation and characteristics of denitrifying bacterium strain GW1[J].Journal of Hebei University of Science and Technology,2012,33(2):184-189.

    [22]Lu S M,Liao M J,Xie C X,et al.Seasonal dynamics of ammonia-oxidizing microorganisms in freshwater aquaculture ponds[J].Annals of Microbiology,2014,65(2):651-657.

    [23]諸葛健,李華鐘.微生物學(xué)[M].北京:科學(xué)出版社,2009.

    ZHU Ge-jian,LI Hua-zhong.Microbiology[M].Beijing:Science Press, 2009.

    [24]董蓮華,楊金水,袁紅莉.氨氧化細(xì)菌的分子生態(tài)學(xué)研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2008,19(6):1381-1388.

    DONG Lian-hua,YANG Jin-shui,YUAN Hong-li.Research advances in molecular ecology of ammonia oxidizing bacteria[J].Chinese Journal of Applied Ecology,2008,19(6):1381-1388.

    [25]Wuchter C,Abbas B,Coolen M J,et al.Archaeal nitrification in the ocean[J].Proceedings of the National Academy of Sciences of the United States of America,2006,103(33):12317-12322.

    [26]杜雪晴,廖新俤,吳銀寶,等.有機(jī)廢棄物好氧堆肥系統(tǒng)中氨氧化微生物的研究進(jìn)展[J].家畜生態(tài)學(xué)報(bào),2014(9):1-7.

    DU Xue-qing,LIAO Xin-di,WU Xin-bao,et al.Advances of ammonia oxidizing microorganisms in organic waste aerobic composting system [J].ActaEcologae Animalis Domastici,2014(9):1-7.

    [27]虞 泳,曾光明,陳耀寧,等.農(nóng)業(yè)廢物好氧堆肥中氨氧化細(xì)菌的群落結(jié)構(gòu)[J].環(huán)境科學(xué),2011,32(10):3067-3072.

    YU Yong,ZENG Guang-ming,CHEN Yao-ning,et al.Succession of ammonia-oxidizing bacteria community structure during agricultural waste composting[J].Chinese Journal of Environment Science,2011, 32(10):3067-3072.

    [28]Enwall K,Nyberg K,Bertilsson S,et al.Long-term impact of fertilization on activity and composition of bacterial communities and metabolic guilds in agricultural soil[J].Soil Biology&Biochemistry,2007,39(1):106-115.

    [29]Valentín-Vargas A,Root R A,Neilson J W,et al.Environmental factors influencing the structural dynamics of soil microbial communitiesduring assisted phytostabilization of acid-generating mine tailings:A mesocosm experiment[J].Science of the Total Environment,2014, 500/501:314-324.

    [30]Zhang J C,Zeng G M,Chen Y N,et al.Effects of physico-chemical parameters on the bacterial and fungal communities during agricultural waste composting[J].Bioresource Technology,2011,102(3):2950-2956.

    [31]曾光明,黃國(guó)和,袁興中,等.堆肥環(huán)境生物與控制[M].北京:科學(xué)出版社,2006.

    ZENG Guang-ming,HUANG Guo-he,YUAN Xing-zhong,et al.Compost environmental biology and control[M].Beijing:Science Press, 2006.

    [32]Yamamoto N,Oishi R,Suyama Y,et al.Ammonia-oxidizing bacteria rather than ammonia-oxidizing archaea were widely distributed in animal manure composts from field-scale facilities[J].Microbes and Environments,2012,27:519-524.

    [33]Wang S Y,Wang Y,Feng X J,et al.Quantitative analyses of ammoniaoxidizing archaea and bacteria in the sediments of four nitrogen-rich wetlands in China[J].Applied Microbiology&Biotechnology,2011,90 (2):779-787.

    [34]宋亞娜,Marschner P,張福鎖,等.小麥/蠶豆、玉米/蠶豆和小麥/玉米間作對(duì)要際細(xì)菌群落結(jié)構(gòu)的影響[J].生態(tài)學(xué)報(bào),2006,26(7):2268-2274.

    SONG Ya-na,Marschner P,ZHANG Fu-suo,et al.Effect of intercropping on bacterial community composition in rhizoshpere of wheat (Triticum aestivum L.),maize(Zea mays L.)and faba bean(Vicia fata L.)[J].ActaecologicaSinica,2006,26(7):2268-2274.

    Research of ammonia oxidizing bacterial community structure and its correlation with environmental factors in cow manure composting

    SUN Xue-wei,XU Xiu-hong*,MENG Qing-xin,CHENG Li-jun,ZHANG Wen-hao,MEN Meng-qi,XU Ben-shu,SUN Yu (College of Resources and Environment,Northeast Agricultural University,Harbin 150030,China)

    The dynamic succession of ammonia oxidizing bacterial community structure and quantity were investigated using PCR-DGGE (polymerase chain reaction and denaturing gradient gel electrophoresis)and real-time PCR,the correlations between ammonia oxidizing bacterial community structure and physic-chemical factors were evaluated using redundancy analysis(RDA)in cow manure composting. The results showed that during composting the range of ammonia oxidizing bacteria were 1.62×106~3.8×107copies·g-1,the amount of ammonia oxidizing bacteria reached the peak on 14 day;DGGE profile analysis showed that the community composition of ammonia oxidizing bacteria obviously varied in different periods of composting.Shannon index of ammonia oxidizing bacteria changed with composting temperature,reaching the maximum(2.671 8)on the 7th day.Phylogenetic analysis showed that ammonia oxidizing bacteria were clustered into two groups,Uncultured Nitrosomonas sp.and Uncultured Nitrosospira sp.,which belonged to the β-Proteobacteria.the Uncultured Nitrosospira sp.which accounted for about 60%of ammonia oxidizing bacteria were predominant.RDA showed that ammonia and temperature(P<0.05) effected significantly on the community structure of ammonia oxidizing bacteria.Significant changes of ammonia oxidizing bacterial community structure occurred in cow manure composting,and it is confirmed that the changes of ammonia oxidizing bacterial community structure were affected by environmental factors.

    cow manure composting;PCR-DGGE;Real-time PCR;ammonia oxidizing bacteria;redundancy analysis

    X713

    A

    1672-2043(2017)01-0189-09

    10.11654/jaes.2016-0960

    2016-07-25

    孫雪微(1991—),女,碩士研究生,研究方向?yàn)槲⑸飳W(xué)。E-mail:910236195@qq.com

    *通信作者:許修宏 E-mail:xuxiuhong@neau.edu.cn

    國(guó)家自然科學(xué)基金項(xiàng)目(31372351,31672469)

    Project supported:T卜e National Natural Science Foundation of C卜ina(31372351,31672469)

    孫雪微,許修宏,孟慶欣,等.牛糞堆肥中氨氧化細(xì)菌群落結(jié)構(gòu)及其與環(huán)境因子相關(guān)性研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2017,36(1):189-197.

    SUN Xue-wei,XU Xiu-hong,MENG Qing-xin,et al.Research of ammonia oxidizing bacterial community structure and its correlation with environmental factors in cow manure composting[J].Journal of Agro-Environment Science,2017,36(1):189-197.

    猜你喜歡
    銨態(tài)氮牛糞群落
    石牛糞金
    意林彩版(2022年2期)2022-05-03 10:23:56
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    丟失的牛糞
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    不同質(zhì)地土壤銨態(tài)氮吸附/解吸特征
    不同鹽堿化土壤對(duì)NH+4吸附特性研究
    屎殼郎大戰(zhàn)牛糞
    野保糗事之撿牛糞
    有機(jī)質(zhì)對(duì)城市污染河道沉積物銨態(tài)氮吸附-解吸的影響*
    春季和夏季巢湖浮游生物群落組成及其動(dòng)態(tài)分析
    久久久国产欧美日韩av| 女人高潮潮喷娇喘18禁视频| 亚洲欧美日韩无卡精品| 日韩欧美一区视频在线观看| 一区福利在线观看| 久久亚洲精品不卡| 精品国产美女av久久久久小说| 亚洲片人在线观看| 精品乱码久久久久久99久播| 日本精品一区二区三区蜜桃| 久久香蕉国产精品| 一个人免费在线观看的高清视频| 黄色视频,在线免费观看| 交换朋友夫妻互换小说| 国产精华一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 国产亚洲精品一区二区www| 午夜精品国产一区二区电影| 色综合站精品国产| 别揉我奶头~嗯~啊~动态视频| 国产片内射在线| 9热在线视频观看99| 成人三级黄色视频| 天堂中文最新版在线下载| 日本wwww免费看| 亚洲色图综合在线观看| 好男人电影高清在线观看| 国产欧美日韩精品亚洲av| 中亚洲国语对白在线视频| 精品欧美一区二区三区在线| 亚洲一区中文字幕在线| 国产又爽黄色视频| 亚洲黑人精品在线| 国产精品久久电影中文字幕| 亚洲国产精品合色在线| 久久精品aⅴ一区二区三区四区| 男女之事视频高清在线观看| 久久久国产一区二区| 亚洲自偷自拍图片 自拍| 新久久久久国产一级毛片| 视频区欧美日本亚洲| 国产熟女午夜一区二区三区| 麻豆一二三区av精品| 黄色 视频免费看| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 国产精品一区二区在线不卡| 久久久久久免费高清国产稀缺| 国产精品国产av在线观看| 一个人观看的视频www高清免费观看 | 午夜精品国产一区二区电影| 久久久精品欧美日韩精品| 久99久视频精品免费| 精品国产乱子伦一区二区三区| 久久精品91蜜桃| 久久天堂一区二区三区四区| 大香蕉久久成人网| 两个人看的免费小视频| 999久久久精品免费观看国产| 国产一卡二卡三卡精品| 久久精品亚洲av国产电影网| 99精品在免费线老司机午夜| 欧美日韩亚洲高清精品| 亚洲精品美女久久久久99蜜臀| 国产精品成人在线| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 女人精品久久久久毛片| 岛国在线观看网站| 国产又色又爽无遮挡免费看| 18美女黄网站色大片免费观看| 多毛熟女@视频| 亚洲av五月六月丁香网| 欧美黄色片欧美黄色片| 欧美日韩福利视频一区二区| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 男女之事视频高清在线观看| 国产精品免费一区二区三区在线| 欧美精品一区二区免费开放| 免费观看精品视频网站| 国产精品影院久久| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区 | 国产成人影院久久av| 最近最新中文字幕大全电影3 | 丰满饥渴人妻一区二区三| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 法律面前人人平等表现在哪些方面| 精品福利观看| 丰满的人妻完整版| 18禁裸乳无遮挡免费网站照片 | 成年人黄色毛片网站| 黄片播放在线免费| 神马国产精品三级电影在线观看 | 亚洲专区国产一区二区| 91精品三级在线观看| 精品熟女少妇八av免费久了| 欧美乱色亚洲激情| 日韩欧美一区视频在线观看| 日韩有码中文字幕| www日本在线高清视频| 精品一区二区三卡| 免费看十八禁软件| 一级毛片精品| 国产成人系列免费观看| 久久精品成人免费网站| 久久国产亚洲av麻豆专区| 99国产极品粉嫩在线观看| 黑丝袜美女国产一区| 久久久国产成人精品二区 | 欧美乱妇无乱码| 我的亚洲天堂| 99国产极品粉嫩在线观看| 欧美日韩乱码在线| 久久国产精品人妻蜜桃| 久久久久亚洲av毛片大全| 亚洲狠狠婷婷综合久久图片| 高清在线国产一区| 国产精品影院久久| 一边摸一边抽搐一进一小说| 久久人人爽av亚洲精品天堂| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 午夜免费鲁丝| 精品久久久久久,| 啪啪无遮挡十八禁网站| 看免费av毛片| 久久人人精品亚洲av| 99国产精品一区二区蜜桃av| 亚洲avbb在线观看| 成人特级黄色片久久久久久久| 欧美日本中文国产一区发布| 99久久综合精品五月天人人| 热re99久久精品国产66热6| 在线天堂中文资源库| 国产成人一区二区三区免费视频网站| 国产乱人伦免费视频| 久久中文看片网| 久久久国产成人免费| 国产亚洲精品综合一区在线观看 | 国产精品一区二区在线不卡| av天堂在线播放| 久久99一区二区三区| 国产精品永久免费网站| 久久久久亚洲av毛片大全| 亚洲国产精品一区二区三区在线| 人人妻人人添人人爽欧美一区卜| 国产野战对白在线观看| 亚洲五月天丁香| 中文字幕最新亚洲高清| 最好的美女福利视频网| 亚洲国产欧美一区二区综合| xxxhd国产人妻xxx| 色在线成人网| 91国产中文字幕| 免费久久久久久久精品成人欧美视频| 国产激情欧美一区二区| 亚洲精品一卡2卡三卡4卡5卡| 大型黄色视频在线免费观看| 国产一卡二卡三卡精品| 久久亚洲真实| 久久久久国产一级毛片高清牌| 我的亚洲天堂| 亚洲欧美激情在线| 欧美午夜高清在线| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 日韩大尺度精品在线看网址 | svipshipincom国产片| 咕卡用的链子| 国产精品综合久久久久久久免费 | 久久人人97超碰香蕉20202| 免费在线观看完整版高清| 亚洲aⅴ乱码一区二区在线播放 | 我的亚洲天堂| 亚洲精品成人av观看孕妇| 日本欧美视频一区| 亚洲九九香蕉| 99久久综合精品五月天人人| 亚洲五月色婷婷综合| 免费高清视频大片| 久久久久久亚洲精品国产蜜桃av| 美女国产高潮福利片在线看| 免费在线观看完整版高清| 久久这里只有精品19| 成年版毛片免费区| 男女床上黄色一级片免费看| 国产精品秋霞免费鲁丝片| a级毛片黄视频| 中文字幕人妻丝袜一区二区| 老司机午夜十八禁免费视频| 日韩精品中文字幕看吧| 深夜精品福利| 日本wwww免费看| 亚洲熟女毛片儿| 91精品三级在线观看| 性少妇av在线| 亚洲专区国产一区二区| 亚洲精品国产精品久久久不卡| 99久久综合精品五月天人人| 动漫黄色视频在线观看| 国产成人精品在线电影| 精品久久久久久成人av| 午夜视频精品福利| 亚洲精品中文字幕一二三四区| 成人18禁高潮啪啪吃奶动态图| 国产精品影院久久| 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 韩国av一区二区三区四区| 国产一区在线观看成人免费| 日韩欧美免费精品| 操美女的视频在线观看| 在线观看66精品国产| 国产免费现黄频在线看| 国产精品 欧美亚洲| 91字幕亚洲| 人成视频在线观看免费观看| 久久青草综合色| 亚洲国产精品999在线| 久久性视频一级片| 女人爽到高潮嗷嗷叫在线视频| 91麻豆精品激情在线观看国产 | 91成人精品电影| 中文字幕另类日韩欧美亚洲嫩草| 日韩有码中文字幕| 一本大道久久a久久精品| 欧美老熟妇乱子伦牲交| aaaaa片日本免费| 黄色成人免费大全| 琪琪午夜伦伦电影理论片6080| 欧洲精品卡2卡3卡4卡5卡区| 天天躁夜夜躁狠狠躁躁| 亚洲精华国产精华精| 99国产精品99久久久久| 婷婷六月久久综合丁香| 久久久精品国产亚洲av高清涩受| 国产精品国产高清国产av| 国产三级黄色录像| √禁漫天堂资源中文www| 男人舔女人的私密视频| 亚洲五月天丁香| 精品国产国语对白av| 精品国产美女av久久久久小说| 操美女的视频在线观看| 动漫黄色视频在线观看| 亚洲精品中文字幕一二三四区| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| 精品熟女少妇八av免费久了| 亚洲欧美精品综合一区二区三区| videosex国产| 搡老乐熟女国产| 午夜影院日韩av| 国产欧美日韩一区二区精品| 女性被躁到高潮视频| av有码第一页| 黄色丝袜av网址大全| 亚洲伊人色综图| 在线观看免费日韩欧美大片| 欧美人与性动交α欧美软件| 久久久水蜜桃国产精品网| www.精华液| 99久久人妻综合| 久久精品成人免费网站| 成人永久免费在线观看视频| 亚洲久久久国产精品| 婷婷丁香在线五月| 91字幕亚洲| 高清毛片免费观看视频网站 | 啦啦啦 在线观看视频| 色尼玛亚洲综合影院| 91成人精品电影| 99香蕉大伊视频| www日本在线高清视频| 国产免费现黄频在线看| 久久热在线av| 免费少妇av软件| 三级毛片av免费| 天堂俺去俺来也www色官网| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人| 另类亚洲欧美激情| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 亚洲五月婷婷丁香| 色婷婷久久久亚洲欧美| 黄色视频不卡| 91大片在线观看| 精品一区二区三卡| 精品久久久久久电影网| 精品一区二区三区视频在线观看免费 | 午夜老司机福利片| 丰满饥渴人妻一区二区三| 欧美精品亚洲一区二区| 国产精品乱码一区二三区的特点 | 一个人观看的视频www高清免费观看 | 男人的好看免费观看在线视频 | 日日摸夜夜添夜夜添小说| 欧美日韩亚洲高清精品| 一级,二级,三级黄色视频| 亚洲一区二区三区不卡视频| 日韩三级视频一区二区三区| 国产av一区在线观看免费| 国产在线观看jvid| 国产精品 国内视频| 国产伦一二天堂av在线观看| 国产午夜精品久久久久久| 国产人伦9x9x在线观看| 男女下面进入的视频免费午夜 | 叶爱在线成人免费视频播放| 日日夜夜操网爽| 18禁国产床啪视频网站| 成人亚洲精品av一区二区 | 伦理电影免费视频| 美女大奶头视频| 精品久久久久久电影网| 女人被躁到高潮嗷嗷叫费观| 嫁个100分男人电影在线观看| 欧美精品亚洲一区二区| а√天堂www在线а√下载| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 99久久国产精品久久久| 亚洲精品国产区一区二| 高清毛片免费观看视频网站 | 日本免费a在线| 美女福利国产在线| 99在线人妻在线中文字幕| 久久影院123| 黄色女人牲交| 人妻久久中文字幕网| 国产成人系列免费观看| 国产97色在线日韩免费| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩综合在线一区二区| 在线观看免费视频日本深夜| 一个人观看的视频www高清免费观看 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| ponron亚洲| 久久久国产成人免费| 热re99久久精品国产66热6| 午夜免费激情av| 国产精品免费视频内射| 精品福利永久在线观看| 国产精品国产高清国产av| 狠狠狠狠99中文字幕| 成人永久免费在线观看视频| 久久中文字幕一级| 咕卡用的链子| 大陆偷拍与自拍| 久久精品国产综合久久久| 欧美日韩中文字幕国产精品一区二区三区 | 美女扒开内裤让男人捅视频| 久久狼人影院| 亚洲精品国产一区二区精华液| 精品国产一区二区三区四区第35| 成人三级做爰电影| 国产一区在线观看成人免费| 精品日产1卡2卡| 亚洲第一av免费看| 精品午夜福利视频在线观看一区| 丁香欧美五月| 久久草成人影院| 亚洲欧美日韩无卡精品| 亚洲第一欧美日韩一区二区三区| 在线观看一区二区三区激情| 久久草成人影院| 久久国产精品男人的天堂亚洲| 久久精品国产99精品国产亚洲性色 | 动漫黄色视频在线观看| 国产色视频综合| 国产无遮挡羞羞视频在线观看| 香蕉国产在线看| 丝袜在线中文字幕| 亚洲色图 男人天堂 中文字幕| 999久久久精品免费观看国产| 久久精品人人爽人人爽视色| 国产成人av激情在线播放| 午夜福利,免费看| 两个人看的免费小视频| 岛国视频午夜一区免费看| 正在播放国产对白刺激| 巨乳人妻的诱惑在线观看| www国产在线视频色| 最好的美女福利视频网| 欧美黄色淫秽网站| 国产成人系列免费观看| 精品日产1卡2卡| 午夜a级毛片| 国产精品 欧美亚洲| 他把我摸到了高潮在线观看| 国产aⅴ精品一区二区三区波| 亚洲成人国产一区在线观看| 欧美久久黑人一区二区| 深夜精品福利| 一级毛片女人18水好多| 欧美日本中文国产一区发布| 黄色毛片三级朝国网站| 淫秽高清视频在线观看| netflix在线观看网站| 久久久精品欧美日韩精品| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩乱码在线| 亚洲国产精品合色在线| 亚洲国产欧美一区二区综合| 国产精品亚洲一级av第二区| 国产色视频综合| 成人av一区二区三区在线看| 日本一区二区免费在线视频| 妹子高潮喷水视频| 久久精品人人爽人人爽视色| 老司机午夜福利在线观看视频| 18禁黄网站禁片午夜丰满| 国产成人av激情在线播放| 久久久久国产精品人妻aⅴ院| 国产精品国产av在线观看| 午夜精品久久久久久毛片777| 国产精品电影一区二区三区| 18禁观看日本| 亚洲av美国av| 无人区码免费观看不卡| 又黄又爽又免费观看的视频| 午夜视频精品福利| 亚洲精品中文字幕一二三四区| 久久国产精品人妻蜜桃| 国产精品影院久久| 女性生殖器流出的白浆| 亚洲人成伊人成综合网2020| 高清欧美精品videossex| 大型黄色视频在线免费观看| 午夜免费观看网址| 欧美日韩一级在线毛片| www日本在线高清视频| 国产xxxxx性猛交| 亚洲国产看品久久| 久久久久久久久久久久大奶| 亚洲人成网站在线播放欧美日韩| 俄罗斯特黄特色一大片| 成年版毛片免费区| 亚洲黑人精品在线| 亚洲成人免费电影在线观看| 午夜精品久久久久久毛片777| 成人特级黄色片久久久久久久| 中文字幕人妻熟女乱码| 亚洲成人久久性| 日本撒尿小便嘘嘘汇集6| 欧美一区二区精品小视频在线| 精品一区二区三区四区五区乱码| 国产在线观看jvid| 亚洲精品国产色婷婷电影| 一a级毛片在线观看| 国产aⅴ精品一区二区三区波| 啦啦啦 在线观看视频| 最近最新中文字幕大全电影3 | 极品教师在线免费播放| 欧美日韩瑟瑟在线播放| 久久99一区二区三区| 黄色成人免费大全| 夜夜爽天天搞| 日韩精品中文字幕看吧| 啪啪无遮挡十八禁网站| 国产乱人伦免费视频| 中文欧美无线码| 嫩草影视91久久| 老鸭窝网址在线观看| 精品国产超薄肉色丝袜足j| 午夜老司机福利片| 久久国产精品男人的天堂亚洲| 国产区一区二久久| 一区二区日韩欧美中文字幕| 淫妇啪啪啪对白视频| 免费av毛片视频| 看免费av毛片| 精品免费久久久久久久清纯| 久久精品亚洲精品国产色婷小说| 中文字幕人妻丝袜制服| 无限看片的www在线观看| 国产精品爽爽va在线观看网站 | 99在线视频只有这里精品首页| 亚洲精品在线观看二区| 18禁黄网站禁片午夜丰满| 久久国产乱子伦精品免费另类| 法律面前人人平等表现在哪些方面| 国产xxxxx性猛交| 国产在线精品亚洲第一网站| 日韩视频一区二区在线观看| 国产熟女xx| 夜夜爽天天搞| 欧美日韩亚洲高清精品| 99久久国产精品久久久| 一区福利在线观看| 日日夜夜操网爽| 成人免费观看视频高清| 老熟妇乱子伦视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 最近最新中文字幕大全免费视频| 亚洲国产精品合色在线| tocl精华| 丝袜在线中文字幕| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产 | 最新美女视频免费是黄的| 久久精品国产清高在天天线| 夜夜躁狠狠躁天天躁| 亚洲精品美女久久av网站| 成人黄色视频免费在线看| 欧美一区二区精品小视频在线| 中出人妻视频一区二区| 91精品三级在线观看| 91九色精品人成在线观看| 国产亚洲欧美在线一区二区| 国产精品久久久久成人av| 视频在线观看一区二区三区| 亚洲成人免费电影在线观看| 欧美色视频一区免费| 美女国产高潮福利片在线看| 麻豆av在线久日| 中文亚洲av片在线观看爽| 欧美色视频一区免费| 亚洲第一av免费看| 亚洲精品久久午夜乱码| 99精品欧美一区二区三区四区| 看黄色毛片网站| 国产精品久久视频播放| 999久久久精品免费观看国产| 成人国语在线视频| a级片在线免费高清观看视频| 久久久久久久久免费视频了| 精品国产乱码久久久久久男人| 欧美日韩精品网址| 午夜免费激情av| 亚洲美女黄片视频| 男女高潮啪啪啪动态图| 国产免费av片在线观看野外av| 久久 成人 亚洲| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品久久久久久毛片| 亚洲欧美激情综合另类| 老鸭窝网址在线观看| 一进一出抽搐gif免费好疼 | 精品国产一区二区三区四区第35| 亚洲一区高清亚洲精品| 国产精品二区激情视频| www日本在线高清视频| 波多野结衣一区麻豆| 成人国产一区最新在线观看| 色尼玛亚洲综合影院| 一个人免费在线观看的高清视频| 丁香六月欧美| 老鸭窝网址在线观看| 韩国精品一区二区三区| 国产97色在线日韩免费| 国产精品一区二区三区四区久久 | a级毛片在线看网站| 一边摸一边做爽爽视频免费| 这个男人来自地球电影免费观看| 午夜影院日韩av| 黄片播放在线免费| 欧美日韩中文字幕国产精品一区二区三区 | 国产区一区二久久| 国产成年人精品一区二区 | 啪啪无遮挡十八禁网站| 变态另类成人亚洲欧美熟女 | 最近最新中文字幕大全电影3 | 曰老女人黄片| 国内毛片毛片毛片毛片毛片| 老汉色∧v一级毛片| 一级片'在线观看视频| 黄色 视频免费看| av网站在线播放免费| 亚洲久久久国产精品| 欧美激情久久久久久爽电影 | 欧美另类亚洲清纯唯美| 男女做爰动态图高潮gif福利片 | 夜夜躁狠狠躁天天躁| 男女下面插进去视频免费观看| 日韩免费av在线播放| 久久久久精品国产欧美久久久| 成在线人永久免费视频| 亚洲国产精品合色在线| 久久久精品国产亚洲av高清涩受| 黄片小视频在线播放| 最新美女视频免费是黄的| 亚洲中文av在线| 午夜两性在线视频| 久久久久国产一级毛片高清牌| 亚洲美女黄片视频| 一级毛片女人18水好多| 麻豆av在线久日| 亚洲专区中文字幕在线| 国产av一区在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 国产乱人伦免费视频| 真人一进一出gif抽搐免费| 国产高清激情床上av| 成人永久免费在线观看视频| 国产熟女午夜一区二区三区| 国产亚洲精品第一综合不卡| 99精品在免费线老司机午夜| 1024视频免费在线观看| 亚洲成人免费av在线播放| 首页视频小说图片口味搜索| 欧美日本中文国产一区发布| 变态另类成人亚洲欧美熟女 |