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

    北運(yùn)河沉積物中氨氧化微生物的群落特征

    2015-11-17 09:25:58鮑林林陳永娟王曉燕首都師范大學(xué)資源環(huán)境與旅游學(xué)院北京00048首都師范大學(xué)首都圈水環(huán)境研究中心北京00048
    中國環(huán)境科學(xué) 2015年1期
    關(guān)鍵詞:北運(yùn)河古菌拷貝數(shù)

    鮑林林,陳永娟,王曉燕,2*(.首都師范大學(xué)資源環(huán)境與旅游學(xué)院,北京 00048;2.首都師范大學(xué)首都圈水環(huán)境研究中心,北京 00048)

    北運(yùn)河沉積物中氨氧化微生物的群落特征

    鮑林林1,陳永娟1,王曉燕1,2*(1.首都師范大學(xué)資源環(huán)境與旅游學(xué)院,北京 100048;2.首都師范大學(xué)首都圈水環(huán)境研究中心,北京 100048)

    采用T-RFLP、RT-qPCR和克隆測序等分子生物學(xué)技術(shù),以氨單加氧酶基因(amoA)為分子標(biāo)記,研究了北運(yùn)河表層沉積物中氨氧化古菌(AOA)和氨氧化細(xì)菌(AOB)的群落多樣性、豐度、系統(tǒng)發(fā)育及其與環(huán)境因子的響應(yīng)關(guān)系.結(jié)果表明,沉積物中AOB的群落多樣性和豐度均高于AOA,是北運(yùn)河沉積物中氨氧化過程的主要功能微生物.沉積物中氨氧化微生物群落結(jié)構(gòu)沿干流和支流存在明顯的空間分異,而AOA的種類組成空間差異較??;沉積物的氨氮)和硝態(tài)氮)是影響氨氧化微生物群落特征的主要因子,AOB對環(huán)境變化的敏感性更高;AOA和AOB的amoA基因拷貝數(shù)分別為1.32×105~1.91×106copies/g、5.39×105~8.3×106copies/g.閘壩下游沉積物的氨氧化微生物豐度最高.系統(tǒng)發(fā)育分析表明,amoA基因序列多屬于土壤/沉積物分支,較多AOB的克隆序列與土壤亞硝化螺菌屬(Nitrosospira)的類群相似性可達(dá)98%.受污水處理廠退水的影響,部分amoA基因序列與污水處理廠廢水和活性污泥中發(fā)現(xiàn)的類群同源性高.污染物質(zhì)來源、支流匯入和閘壩攔截對河流沉積物氨氧化微生物的群落特征影響顯著.

    沉積物;amoA;群落結(jié)構(gòu);豐度;系統(tǒng)發(fā)育

    每年有超過一半的氮固定由人類活動產(chǎn)生,大量活性氮輸入生物圈,加劇了全球氮循環(huán)的氮負(fù)荷干擾了相應(yīng)微生物的生態(tài)功能[1].氨氧化作用(硝化過程的第一步)作為氮循環(huán)的關(guān)鍵限速步驟,是全球氮循環(huán)的中心環(huán)節(jié),主要由功能微生物氨氧化古菌(AOA)和氨氧化細(xì)菌(AOB)所驅(qū)動[2-3]. AOA和AOB在不同的環(huán)境條件下,兩者相對的群落結(jié)構(gòu)、多樣性、豐度和氨氧化效率優(yōu)勢不盡相同,并且受到pH值、、鹽度和有機(jī)碳等環(huán)境因子的影響[4-6].古菌和細(xì)菌對硝化過程的相對貢獻(xiàn)率,是目前氮循環(huán)微生物生態(tài)學(xué)的研究重點(diǎn),在大多數(shù)生態(tài)環(huán)境如土壤和海洋的氨氧化過程中AOA居主導(dǎo)地位[3,7],但是AOA和AOB在淡水生態(tài)系統(tǒng)尤其是沉積物中的生態(tài)特征、功能差異及多樣化環(huán)境要素的影響機(jī)制仍需要進(jìn)一步的研究[8-9].

    北運(yùn)河上游河段(沙河閘-楊洼閘)主要流經(jīng)北京市北部和東部地區(qū),流域范圍覆蓋了北京市城市化水平最高、人類活動最頻繁的區(qū)域,周圍環(huán)境受人類活動影響而變化劇烈[10].作為北京市最主要的泄洪、排污河道,大量未經(jīng)處理的工農(nóng)業(yè)廢水、城市生活污水以及污水處理廠的退水以點(diǎn)源或非點(diǎn)源的形式排入河道,導(dǎo)致水體氨氮含量嚴(yán)重超標(biāo)[11-12].河道水流緩慢導(dǎo)致沉積物大量淤積,而污染物質(zhì)的沉降、吸附使得沉積物成為內(nèi)源污染源,其中氮循環(huán)功能微生物的活性和氧化還原條件決定了過量氮素的有效轉(zhuǎn)化和輸出[13],對北運(yùn)河的水質(zhì)改善具有重要作用.范改娜等[14]和張建偉等[15]在北運(yùn)河岸邊濕地和溫榆河(北運(yùn)河上游)沉積物中已分別發(fā)現(xiàn),AOA和AOB的普遍存在.本研究利用分子生物學(xué)技術(shù),對比研究了北運(yùn)河沉積物中氨氧化微生物(AOA和AOB)的群落特征和系統(tǒng)發(fā)育,并對其與環(huán)境因子的響應(yīng)關(guān)系進(jìn)行定性和定量的分析,以探討AOA和AOB在城市河流氨氧化過程中的作用和功能差異,為微生物參與氨氧化過程的調(diào)控作用和建立北運(yùn)河生物地球化學(xué)循環(huán)模型提供理論基礎(chǔ),同時也為控制城市河流氮素污染、協(xié)調(diào)水生態(tài)平衡提供科學(xué)依據(jù).

    1 材料與方法

    1.1 樣品采集

    2013年7月,于北運(yùn)河京區(qū)流域范圍主干水系上進(jìn)行實(shí)驗(yàn)樣品的采集.沿途共設(shè)置5個樣點(diǎn),依次編號為1、2、3、4、5(圖1).樣點(diǎn)設(shè)置考慮了水庫、閘壩和支流匯入等影響河流水質(zhì)的主要因素:位于沙河水庫(水系上游區(qū)域)堤壩下游的1號點(diǎn)和楊洼閘(水系出京的最后關(guān)卡)閘門下游的5號點(diǎn),位于水系主要支流(壩河和涼水河)匯入口的3、4號樣點(diǎn),以及位于馬坊橋的2號點(diǎn).每個樣點(diǎn)分別采集上覆水和表層沉積物(3~5cm),沉積物用無菌袋封裝,冷藏保存(4℃),運(yùn)回實(shí)驗(yàn)室后-20℃凍存以備后續(xù)分子實(shí)驗(yàn)使用.

    圖1 采樣點(diǎn)分布示意Fig.1 The distribution of sampling sites 1.沙河水庫壩后;2.馬坊橋;3.壩河口;4.榆林莊橋;5.楊洼閘后

    1.2 樣品理化分析

    上覆水pH值由Mettler Toledo Delta 320pH計(jì)測定,溶解氧濃度(DO)采用取樣現(xiàn)場進(jìn)行硫酸錳和堿性試劑固定的碘量法(GB 7489-87)測定,總有機(jī)碳(TOC)和總氮(TN)使用Multi NC2100(Analytik Jena,德國)測定,氨氮、亞硝氮)和硝態(tài)氮)使用FIAstar流動注射分析儀(FOSS,丹麥)測定.

    1.3 沉積物DNA提取及amoA基因PCR擴(kuò)增

    稱取0.3g沉積物樣品,采用Power Soil DNA試劑盒(MO BIO Laboratories Inc.,美國)進(jìn)行樣品DNA提取,作為PCR擴(kuò)增的DNA模板.AOA和AOB的amoA基因PCR擴(kuò)增反應(yīng)體系為:DNA模版10 μL,前后引物各2μL,2×Taq PCR Master Mix(天根生化科技,北京)25 μL,加雙蒸水(ddH2O)到總體積為50 μL.PCR擴(kuò)增引物和溫度程序設(shè)置如表1所示.擴(kuò)增產(chǎn)物用1%的瓊脂凝膠電泳檢測.用于末端限制性片段長度多態(tài)性技術(shù)(TRFLP)分析的PCR擴(kuò)增產(chǎn)物,對前引物(ArchamoAF, amoA-1F)5’端進(jìn)行FAM熒光標(biāo)記.

    表1 氨氧化微生物amoA基因PCR擴(kuò)增引物和條件Table 1 Primers and procedures for PCR amplification of amoA genes

    1.4 T-RFLP分析和RT-qPCR

    利用T-RFLP進(jìn)行氨氧化微生物群落結(jié)構(gòu)分析[18].按照TIANquick Mini Purification Kit純化試劑盒(天根生化科技,北京)對PCR產(chǎn)物純化回收后,采用限制性內(nèi)切酶HhaⅠ(Takara寶生物工程有限公司,大連)進(jìn)行酶切.反應(yīng)體系(20μL)為:PCR純化產(chǎn)物10μL,HhaⅠ酶1μL,10mol/L Buffer 2μL,ddH2O 7μL.37℃培養(yǎng)箱內(nèi)酶切5h后,70℃水浴滅活15min,酶切產(chǎn)物由上?;倒緶y定.

    采用SYBR Green法,由ABI 7500FAST(Applied Biosystems, 美國)型熒光定量PCR儀進(jìn)行amoA的絕對定量分析(實(shí)時熒光定量PCR,RT-qPCR).樣品amoA基因拷貝數(shù)的熒光定量測試由上海美吉生物醫(yī)藥科技有限公司完成.

    1.5 克隆、測序及系統(tǒng)發(fā)育分析

    用Gel Extraction Kit(OMEGA Bio-tek,美國)對PCR產(chǎn)物進(jìn)行切膠純化,將回收產(chǎn)物與pEASY-T1載體(北京全式金生物技術(shù)有限公司)進(jìn)行連接,轉(zhuǎn)入Trans1-T1Phage Resistant化學(xué)感受態(tài)細(xì)胞,在氨芐青霉素平板上培養(yǎng)后進(jìn)行藍(lán)、白斑篩選,并通過克隆檢測引物進(jìn)行菌液PCR,挑選出陽性克隆子送至北京諾賽基因組研究中心有限公司測序.返回序列經(jīng)NCBI提交獲取序列登陸號(AOA:KF856971-KF857161; AOB: KJ093846-KJ094018)后,利用DOTUR軟件以2%的差異度劃分獨(dú)立操作單元(operational taxonomic unit,OTU)[19],將代表序列在GenBank數(shù)據(jù)庫中進(jìn)行BLAST比對,并選取已發(fā)表的氨氧化微生物amoA序列作為參比序列,使用MEGA 4軟件以鄰位相連法(Neighbor-Joining)構(gòu)建系統(tǒng)發(fā)育樹[20].

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

    利用PRIMER 5.0對整理后的T-RFLP數(shù)據(jù)進(jìn)行微生物群落的香農(nóng)多樣性指數(shù)(Shannon-Weiner, H’)和均勻度指數(shù)(E),以及各樣點(diǎn)的相似性(Similarity)分析.計(jì)算公式如下:

    式中:S為該樣樣點(diǎn)微生物群落T-RF片段總數(shù);Pi= Ni/ N,為第i個T-RF片段的T-RFLP圖譜峰面積(Ni)所占總峰面積(N)的百分比.

    采用SPSS16.0和PRIMER 5.0對數(shù)據(jù)進(jìn)行處理.氨氧化微生物群落與環(huán)境因子之間的關(guān)系采用Canoco for windows 4.5軟件進(jìn)行分析.

    2 結(jié)果

    2.1 沉積物及上覆水的理化指標(biāo)

    北運(yùn)河的河道設(shè)有多道閘壩、橡膠壩等水利措施,水流平緩,沉積物淤積明顯.由表2可見,各樣點(diǎn)上覆水TN濃度均很高,其中以為主,高出Ⅴ類水質(zhì)標(biāo)準(zhǔn)值4~8倍.和的濃度相對較低,閘壩后水體(沙河水庫閘和楊洼閘)的濃度明顯高于其他樣點(diǎn).上覆水pH>7.5偏堿性,TOC濃度從上游到下游各點(diǎn)呈先增加后減少的趨勢.各樣點(diǎn)沉積物的pH值基本上低于上覆水,而沉積物中的與的濃度比明顯高于上覆水中的比值,2號點(diǎn)沉積物的和含量最高.位于庫壩、閘壩下游的1、5號樣點(diǎn),水流較快,沉積物泥沙含量明顯高于其他樣點(diǎn),其TN和TOC含量也明顯低于其他樣點(diǎn).沉積物的TN和TOC存在顯著的正相關(guān)關(guān)系(P<0.01).總體看來,上覆水營養(yǎng)物質(zhì)水平對沉積物具有一定的影響作用,而沉積物中氨氧化微生物對的轉(zhuǎn)化效率可能更高.

    表2 表層沉積物和上覆水的理化性質(zhì)Table 2 Physicochemical properties of sediments and overlying waters

    2.2 沉積物氨氧化微生物的多樣性和豐度

    如表3所示,每個T-RF片段至少代表一種微生物,北運(yùn)河沉積物中共檢出33種不同的AOA和44種不同的AOB,總體來看,北運(yùn)河沉積物AOB的多樣性高于AOA.均勻度指數(shù)和多樣性指數(shù)變化一致,5號點(diǎn)(楊洼閘)的氨氧化微生物種數(shù)最少,且amoA的均勻度和多樣性最低;4號點(diǎn)(榆林莊橋)的AOA多樣性最高,2號點(diǎn)(馬坊橋)的AOB多樣性最高,3、4號點(diǎn)AOA多樣性較高于AOB.

    表3 氨氧化微生物amoA的T-RFs和多樣性Table 3 The terminal fragments of amoA genes and the calculated diversity indices

    圖2 各樣點(diǎn)主要T-RFs的相對含量Fig.2 The relative abundance of the dominant T-RFs in each sediment

    圖2所示為5個樣點(diǎn)沉積物氨氧化微生物主要種類(T-RF百分含量大于6)的組成狀況.AOA的優(yōu)勢菌種比較突出,主要為分布在1、2、5號點(diǎn)的169bp和3、4號點(diǎn)的116bp.AOB優(yōu)勢菌的種類數(shù)較多,包括74,65,474,106bp等,且各樣點(diǎn)之間T-RFs組成的差異性較大,可見,北運(yùn)河沉積物中AOB比AOA具有更高的生物多樣性.利用PRIMER 5.0對各樣點(diǎn)氨氧化微生物的T-RFs組成進(jìn)行Similarity聚類分析,5個沉積物的AOA和AOB群落均被主要聚為兩大類:位于水系干流的1、2、5號樣點(diǎn)和位于支流入?yún)R處的3、4號樣點(diǎn).其中2、5號樣點(diǎn)AOA群落組成的相似性高達(dá)81%,1號與2、5號AOA的相似性約70%,而位于支流的3、4號樣點(diǎn)AOA群落相似性也有57%;而AOB相似性最高的樣點(diǎn)3、4號也僅53%,可見,各樣點(diǎn)AOA的群落組成相似性較高,AOB的空間差異大,群落組成多樣性更高.

    圖3 各樣點(diǎn)氨氧化古菌和細(xì)菌amoA基因拷貝數(shù)Fig.3 Abundances of archaea and bacterial amoA genes in each sediment

    RT-qPCR分析沉積物氨氧化微生物的豐度結(jié)果如圖3所示.AOA和AOB的amoA拷貝數(shù)分別為1.32×105~1.91×106copies/g、5.39×105~8.3×106copies/g,細(xì)菌amoA的拷貝數(shù)是古菌的2~11倍,AOB的豐度高于AOA.1、5號樣點(diǎn)的amoA拷貝數(shù)明顯高于其他樣點(diǎn),且1號點(diǎn)的amoA拷貝數(shù)最高,AOB的豐度也最高,4號點(diǎn)(榆林莊橋)的氨氧化微生物豐度最低,氨氧化微生物的豐度主要受細(xì)菌amoA拷貝數(shù)的影響(R2= 0.99,P < 0.01).細(xì)菌和古菌的amoA占總amoA的比例為極顯著的負(fù)相關(guān)關(guān)系(P < 0.01),AOA、AOB在北運(yùn)河沉積物中的群落生長表現(xiàn)為一定的競爭關(guān)系.

    2.3 環(huán)境因子與微生物群落結(jié)構(gòu)的相關(guān)分析

    利用Canoco選擇合適的分析模型,根據(jù)物種的排序結(jié)果,對AOA優(yōu)勢種與沉積物環(huán)境理化因子的關(guān)系進(jìn)行冗余分析(RDA),AOB優(yōu)勢種與環(huán)境因子的關(guān)系進(jìn)行典范對應(yīng)分析(CCA)分析(圖4).

    圖4 各樣點(diǎn)優(yōu)勢種與沉積物環(huán)境因子的關(guān)系Fig.4 The correlation analysis of the dominant T-RFs with sediment environmental factors

    圖5 古菌amoA基因序列的系統(tǒng)發(fā)育樹Fig.5 Phylogenetic tree of representative archaeal amoA sequences and reference sequences from GenBank

    圖6 細(xì)菌amoA基因序列的系統(tǒng)發(fā)育樹Fig.6 Phylogenetic tree of representative bacterial amoA sequences and reference sequences from GenBank

    如圖4所示,各沉積物樣點(diǎn)受環(huán)境因子影響的AOA和AOB群落結(jié)構(gòu)明顯分為兩類:干流(1、 2、5號樣點(diǎn))和支流(3、4號樣點(diǎn)),與Similarity分析結(jié)果一致.AOA優(yōu)勢種與環(huán)境因子的RDA結(jié)果顯示,排序軸AX1和AX2的解釋度分別為85.5%、8.7%,氨氮占總氮的比例(-0.9927)、氧化態(tài)氮占總氮的比例(-0.8351)和)與AX1的相關(guān)性較高與AX2的相關(guān)性較高.CCA分析沉積物AOB優(yōu)勢種與環(huán)境因子的結(jié)果表明,排序軸AX1、AX2和AX3的解釋度分別為48.3%、20.9%和18.6%,與AX1相關(guān)性高的環(huán)境因子有)、)、)和,與AX2相關(guān)性較高的環(huán)境因子為TOC/TN(0.8198,)與AX3相關(guān)性較高的環(huán)境因子有)、TN(0.6857)和pH(-0.6584).總體而言,沉積物TN中、所占比例和的含量是影響氨氧化微生物群落結(jié)構(gòu)的主要因素,且AOB群落結(jié)構(gòu)受多種環(huán)境因子的影響,其對環(huán)境變化的敏感性比AOA高.

    2.4 沉積物氨氧化微生物的系統(tǒng)發(fā)育分析

    如圖5所示,沉積物中AOA的amoA基因序列在系統(tǒng)發(fā)育樹中可劃分為3個分支,其中52.5%的序列屬于第Ⅰ分支,與海岸帶沉積物及農(nóng)業(yè)土壤AOA的親緣關(guān)系較近;第Ⅱ個分支包括了8.1%的主要來自1號樣點(diǎn)的克隆序列,與崇明島東潮灘沉積物的序列(JQ345811)相似度高達(dá)99%.其余39.5%的克隆序列屬于第Ⅲ個分支,與飲用水和污水處理廠水樣中的古菌amoA序列同源性較高.此外,1號樣點(diǎn)amoA序列在3個分支均有分布,2號樣點(diǎn)主要屬于第Ⅲ個分支,而3號樣點(diǎn)僅分布于第Ⅰ分支.綜上所述,北運(yùn)河沉積物AOA的系統(tǒng)發(fā)育主要屬于土壤/沉積物分支類群,又因其高污染的環(huán)境特點(diǎn),部分序列與污水處理廠類群的親緣性較高,主要包括2號點(diǎn)(馬坊橋)的相關(guān)序列,可能與該處強(qiáng)烈人類活動(集市,多生活固體廢棄物)影響有關(guān).

    北運(yùn)河沉積物中AOB的amoA序列也主要劃分為3個分支(圖6),第Ⅰ個分支包括了45%的序列,主要與土壤amoA的序列親緣關(guān)系較近,該分支主要包括3和5號樣點(diǎn)的克隆序列,其中序列5b30與來自森林土壤[21]的Nitrosospira sp. LT2MFa(AY189145)相似性88%,序列5b10與來自農(nóng)田土壤[22]的Nitrosospira sp. Np39-19(AF042170)相似性高達(dá)98%,該分支與亞硝化螺菌屬(Nitrosospira)類群密切相關(guān);12.4%的序列(來自1號點(diǎn))屬于第Ⅱ分支,與污水處理廠污泥中AOB類群相似性極高;第Ⅲ個分支與河口、湖泊沉積物的AOB類群相似,包括42.6%的測試序列,2號樣點(diǎn)所有代表序列從屬于該分支.

    3 討論

    3.1 群落特征與環(huán)境因子的響應(yīng)關(guān)系

    已有研究表明,在很多環(huán)境中AOA的群落特征和功能作用均優(yōu)于AOB,且AOA的多樣性具有明顯的環(huán)境來源特性,是驅(qū)動氨氧化過程的主要微生物[13,23],如在海洋生態(tài)系統(tǒng)和土壤環(huán)境中古菌amoA基因比細(xì)菌的相對豐度高3000~8000倍[24-25],淡水中AOA比AOB更適于厭氧或低氨氮濃度條件[26-27].但是,在紅樹林沉積物[17]、河流沉積物[28]以及氨氮污染嚴(yán)重的濕地(潮灘)[29]中AOB的群落豐度更具優(yōu)勢,除氮反應(yīng)器活性污泥中的AOB群落多樣性和豐度均高于AOA,且AOB對環(huán)境變化的敏感性比AOA高[30],本文針對amoA基因的T-RFLP和RT-qPCR實(shí)驗(yàn)分析也得到相似的結(jié)果,沉積物中AOB的種類數(shù)多于AOA,其中AOA的優(yōu)勢種類主要為169bp和116bp,AOB的優(yōu)勢種類包括74,56,474,106bp等,總體來看AOB的多樣性高于AOA.沉積物中AOB的群落豐度也高于AOA,細(xì)菌amoA的拷貝數(shù)是古菌的2~11倍,總amoA基因拷貝數(shù)與AOB豐度顯著相關(guān),由此可見,AOB可能是北運(yùn)河沉積物中氨氧化過程的主導(dǎo)微生物.

    北運(yùn)河沉積物中AOB的主要類群(分支Ⅰ,45%的amoA序列)可能屬于亞硝化螺菌屬(Nitrosospira).而在GenBank數(shù)據(jù)庫中沒有發(fā)現(xiàn)純培養(yǎng)的AOA已知菌種與北運(yùn)河沉積物的amoA基因序列有較近發(fā)育關(guān)系.沉積物中大部分的AOA類群與土壤/沉積物來源的菌種相似性較高,AOA的系統(tǒng)發(fā)育可能與奇古菌Group I.1b類群密切相關(guān)[2].此外,部分amoA基因序列與污水處理廠的廢水[39]、活性污泥[40]中的類群有良好的親緣關(guān)系,這與北運(yùn)河接納大量污水處理廠退水有關(guān),通常污水處理廠退水中除了大量的有機(jī)質(zhì)和營養(yǎng)元素外,還帶有很多特有的微生物,對河流的水質(zhì)狀況和微生物系統(tǒng)具有很強(qiáng)的擾動(如微生物物種的殖民)[41-42].

    3.2 氨氧化微生物的空間分異

    本研究涉及的5個沉積物樣點(diǎn),對AOA和AOB群落結(jié)構(gòu)進(jìn)行Similarity聚類分析,結(jié)果可分為1、2、5號和3、4號兩類,其中,3、4號樣點(diǎn)分別位于支流壩河和涼水河的入?yún)R處,污水處理廠退水的影響強(qiáng)烈的主要補(bǔ)給來源為再生水和污水混合類, 1、2、5號樣點(diǎn)位于水系干流上主要來源于再生水水源[10].此外,位于上游河段的1、2號樣點(diǎn)沉積物中的AOB種類數(shù)明顯高于AOA,而來自支流(壩河、涼水河)沉積物中的AOA種類數(shù)更多,到下游干流的5號樣點(diǎn),受支流匯入的影響,沉積物中AOA種類數(shù)明顯多于AOB,可見,水系中不同點(diǎn)位主要污染物質(zhì)的來源不同,氨氧化微生物的群落結(jié)構(gòu)也表現(xiàn)出相應(yīng)的空間分異,且支流匯入對河流沉積物群落組成的影響作用顯著.

    熒光定量分析(圖3)發(fā)現(xiàn)位于河道閘壩后的1號(沙河水庫壩后)和5號(楊洼閘后)樣點(diǎn)的amoA基因拷貝數(shù)明顯高于其他樣點(diǎn),Spearman相關(guān)性分析表明,總的amoA基因豐度與沉積物的TN、TOC含量呈顯著負(fù)相關(guān)(P < 0.05),與含量也存在一定的負(fù)相關(guān)關(guān)系,閘壩后的上覆水和沉積物中TN、TOC含量也明顯低于其他樣點(diǎn),楊洼閘后的濃度也相對較低,這可能是由于閘壩頻繁開啟擾動水體有益沉積物微生物增長,而促進(jìn)了氨氧化微生物的豐度發(fā)育和營養(yǎng)物質(zhì)的轉(zhuǎn)化,同時庫壩內(nèi)相對穩(wěn)定的環(huán)境有利于營養(yǎng)物質(zhì)的沉降、截留[11],Wall等[43]通過對Kaskaskia河上的Shelbyville水庫(美國)內(nèi)營養(yǎng)物質(zhì)轉(zhuǎn)化和微生物作用的研究也表明,河流水庫庫區(qū)對水體營養(yǎng)物質(zhì)的攔截和沉積物中微生物的功能作用,對庫區(qū)下游水生態(tài)系統(tǒng)具有積極作用.北運(yùn)河主河道和支流沿途都設(shè)有多道水閘、橡膠壩等水利調(diào)控設(shè)施,具有一定的攔洪、攔污作用,并明顯地影響著沉積物氨氧化微生物的群落結(jié)構(gòu)和空間分布,以及河流營養(yǎng)物質(zhì)的轉(zhuǎn)化效率.

    4 結(jié)論

    4.1 北運(yùn)河沉積物中AOB的群落多樣性高于AOA.由于支流匯入和不同污染物質(zhì)來源的影響,干流和支流沉積物氨氧化微生物群落結(jié)構(gòu)存在明顯分異,各樣點(diǎn)之間AOA的組成相似性較高,AOB的群落分布空間差異較大.

    4.3 北運(yùn)河沉積物中氨氧化微生物的主要類群屬于土壤/沉積物分支,其中,AOB與土壤亞硝化螺菌屬(Nitrosospira)的相似性較高,受污水處理廠退水影響,部分amoA基因序列與污水處理廠廢水和活性污泥中發(fā)現(xiàn)的類群同源性強(qiáng).可見,污染物質(zhì)來源、支流匯入和閘壩攔截對河流沉積物氨氧化微生物的群落特征具有很大的影響作用.

    [1]Fowler D, Coyle M, Skiba U, et al. The global nitrogen cycle in the twenty-first century [J]. Phil. Trans. R. Soc. B, 2013,368(1621):1-13.

    [2]賀紀(jì)正,張麗梅.氨氧化微生物生態(tài)學(xué)與氮循環(huán)研究進(jìn)展 [J].生態(tài)學(xué)報, 2009,29(1):406-415.

    [3]Francis C A, Beman J M, Kuypers M M M. New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation [J]. The ISME Journal, 2007,1:19-27.

    [4]Liu S, Shen L D, Lou L P, et al. Spatial Distribution and Factors Shaping the Niche Segregation of Ammonia-Oxidizing Microorganisms in the Qiantang River, China [J]. Applied and Environmental Microbiology, 2013,79(13):4065-4071.

    [5]Vissers E W, Anselmetti F S, Bodelier P L E, et al. Temporal and spatial coexistence of archaeal and bacterial amoA genes and gene transcripts in Lake Lucerne [J]. Achaea, 2013,2013(2013):1-11.

    [6]Zheng Y L, Hou L J, Newell S, et al. Community Dynamics and Activity of Ammonia-oxidizing Prokaryotes in Intertidal Sediments of the Yangtze Estuary [J]. Applied Microbiology and Biotechnology, 2014,80(1):408-419.

    [7]郭建華,彭永臻.異養(yǎng)硝化、厭氧氨氧化及氨氧化古菌與新的氮循環(huán) [J]. 環(huán)境科學(xué)學(xué)報, 2008,28(8):1489-1498.

    [8]Fernàndez-Guerra A, Casamayor E O. Habitat-associated phylogenetic community patterns of microbial ammonia oxidizers[J]. PLoS ONE, 2012,7(10):e47330.

    [9]Shen J P, Xu Z H, He J Z. Frontiers in the microbial processes of ammonia oxidation in soils and sediments [J]. Journal of Soils and Sediments, 2014,14(6):1023-1029.

    [10]荊紅衛(wèi),張志剛,郭 婧.北京北運(yùn)河水系水質(zhì)污染特征及污染來源分析 [J]. 中國環(huán)境科學(xué), 2013,33(2):319-327.

    [11]于 洋.北運(yùn)河水體中氨氮的氧化過程及微生物響應(yīng)特征 [D].北京:首都師范大學(xué), 2012:63-70.

    [12]張汪壽,李曉秀,王曉燕,等.北運(yùn)河武清段水污染時空變異特征[J]. 環(huán)境科學(xué)學(xué)報, 2012,32(4):836-846.

    [13]龔 俊,宋延靜,張曉黎.海岸帶沉積物中氮循環(huán)功能微生物多樣性 [J]. 生物多樣性, 2013,21(4):433-444.

    [14]范改娜,祝貴兵,王 雨,等.河流濕地氮循環(huán)修復(fù)過程中的新型功能微生物 [J]. 環(huán)境科學(xué)學(xué)報, 2010,30(8):1558-1563.

    [15]張建偉,孫衛(wèi)玲,邵 軍,等.溫榆河中硝化和反硝化基因的Real-time PCR定量 [J]. 環(huán)境科學(xué)研究, 2013,26(1):64-71.

    [16]鮑士旦.土壤農(nóng)化分析 [M]. 北京:中國農(nóng)業(yè)出版社,2010:42-49.

    [17]Li M, Cao H L, Hong Y G, et al. Spatial distribution and abundances of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in mangrove sediments [J]. Applied Microbiology and Biotechnology, 2011,89(4):1243-1254.

    [18]張 巖,朱 敏,劉煥光,等.利用RFLP分析DO對富積床系統(tǒng)中AOB群落結(jié)構(gòu)的影響 [J]. 中國環(huán)境科學(xué), 2014,34(9):2387-2393.

    [19]Schloss P D, Handelsman J. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness [J]. Applied and Environmental Microbiology,2005,71(3):1501-1506.

    [20]Tamura K, Dudley J, Nei M, et al. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0 [J]. Molecular Biology and Evolution, 2007,24(8):1596-1599.

    [21]Mintie A T, Heichen R S, Cromack K Jr, et al. Ammoniaoxidizing bacteria along meadow-to-forest transects in the Oregon Cascade Mountains [J]. Applied and Environmental Microbiology, 2003,69(6):3129-3136.

    [22]Norton J M, Alzerreca J J, Suwa Y, et al. Diversity of ammonia monooxygenase operon in autotrophic ammonia-oxidizing bacteria [J]. Archives of Microbiology, 2002,177(2):139-149.

    [23]Erguder T H, Boon N, Wittebolle, et al. Environmental factors shaping the ecological niches of ammonia-oxidizing archaea [J]. FEMS Microbiology Reviews, 2009,33(5):855-869.

    [24]賈仲君,翁佳華,林先貴,等.氨氧化古菌的生態(tài)學(xué)研究進(jìn)展 [J].微生物學(xué)報, 2010,50(4):431-437.

    [25]Leininger S, Urich T, Schloter M, et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils [J]. Nature, 2006,442(17):806-809.

    [26]Liu Z H, Huang S B, Sun G P, et al. Diversity and abundance of ammonia-oxidizing archaea in the Dongjiang River [J]. Microbiological Research, 2011,166(5):337-345.

    [27]劉晶靜,吳偉祥,丁 穎,等.氨氧化古菌及其在氮循環(huán)中的重要作用 [J]. 應(yīng)用生態(tài)學(xué)報, 2010,21(8):2154-2160.

    [28]Sun W, Xia C, Xu M, et al. Distribution and abundance of archaeal and bacterial ammonia oxidizers in the sediments of the Dongjiang River, a drinking water supply for Hong Kong [J]. Microbes and Environments, 2013,28(4):457-465.

    [29]Zheng Y L, Hou L J, Liu M, et al. Diversity, abundance, and activity of ammonia-oxidizing bacteria and archaea in Chong ming eastern intertidal sediments [J]. Applied Microbiology and Biotechnology, 2013,97(18):8351-8363.

    [30]Jin T, Zhang T, Yan Q M. Characterization and quantification of ammonia-oxidizing archaea (AOA) and bacteria (AOB) in a nitrogen-removing reactor using T-RFLP and qPCR [J]. Applied Microbiology and Biotechnology, 2010,87(3):1167-1176.

    [31]Magalh?es C M, Joye S B, Moreira R M, et al. Effect of salinity and inorganic nitrogen concentrations on nitrification and denitrification rates in intertidal sediments and rocky biofilms of the Douro River estuary, Portugal [J]. Water Research, 2005,39(9): 1783-1794.

    [32]Verhamme D T, Prosser J I, Nicol G W. Ammonia concentration determines differential growth of ammonia-oxidizing archaea and bacteria in soil microcosms [J]. The ISME Journal, 2011,5:1067-1071.

    [33]Jia Z J, Conrad R. Bacteria rather than Archaea dominate microbial ammonia oxidation in an agricultural soil [J]. Environmental Microbiology, 2009,11(7):1658-1671.

    [34]Di H, Cameron K, Shen J, et al. Nitrification driven by bacteria and not archaea in nitrogen-rich grassland soils [J]. Nature Geoscience, 2009,2(9):621-624.

    [35]Chen X, Zhang L M, Shen J P, et al. Abundance and community structure of ammonia-oxidizing archaea and bacteria in an acid paddy soil. Biology and Fertility Soils [J]. 2011,47(3):323-331.

    [36]Nicol G W, Leininger S, Schleper C, et al. The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria [J]. Environmental Microbiology, 2008,10(11):2966-2978.

    [37]王 超,單保慶.子牙河水系水和沉積物好氧氨氧化微生物分布特征 [J]. 環(huán)境科學(xué)學(xué)報, 2012,32(12):2943-2950.

    [38]Zhao D Y, Luo J, Zeng J, et al. Effects of submerged macrophytes on the abundance and community composition of ammonia-oxidizing prokaryotes in a eutrophic lake [J]. Environmental Science and Pollution Research, 2014,21(1): 389-398.

    [39]Bai Y H, Sun Q H, Wen D H, et al. Abundance of ammoniaoxidizing bacteria and archaea in industrial and domestic wastewater treatment systems [J]. FEMS Microbiol. Ecol., 2012,80(2):323-330.

    [40]Langone M, Yan J, Haaijer S C, et al. Coexistence of nitrifying,anammox and denitrifying bacteria in a sequencing batch reactor[J]. Frontiers in Microbiology, 2014,5(28):1-12.

    [41]Cébron A, Garnier J. Nitrobacter and Nitrospira genera as representatives of nitrite-oxidizing bacteria: detection,quantification and growth along the lower Seine River (France)[J]. Water Research, 2005,39(20):4979-4992.

    [42]Wakelin S A, Colloff M J, Kookana R S. Effect of wastewater treatment plant effluent on microbial function and community structure in the sediment of a freshwater stream with variable seasonal flow [J]. Applied and Environmental Microbiology,2008,74(9):2659-2668.

    [43]Wall L G, Tank J L, Royer T V, et al. Spatial and temporal variability in sediment denitrification within an agriculturally influenced reservoir [J]. Biogeochemistry, 2005,76(1):85-111.

    Diversity and abundance of ammonia-oxidizing prokaryotes in surface sediments in Beiyun River.

    BAO Lin-lin1,CHEN Yong-juan1, WANG Xiao-yan1,2*(1.College of Resources, Environment and Tourism, Capital Normal University,Beijing 100048, China;2.Research Center of Aquatic Environment in the Capital Region, Capital Normal University,Beijing 100048, China). China Environmental Science, 2015,35(1):179~189

    The community diversity, abundance, and phylogeny of ammonia-oxidizing archaea (AOA) and bacteria (AOB)in the surface sediments of Beiyun River were investigated using various molecular techniques targeting the ammonia monooxygenase (amoA) genes. Diversity and abundance of bacteria were higher than those of archaea. AOB were the dominant ammonia-oxidizing microorganisms in Beiyun River. For AOB, sediment samples from main stream and tributary were cluster into two categories. In contrast, AOA communities showed relatively slight site-specific difference. Correlation analysis indicated that concentration of ammonia) and the sum of nitrite and nitrate of sediment were the main factors affecting the community compositions of ammonia-oxidizing prokaryotes, and AOB was much more sensitive to environmental change. The amoA gene copy numbers of AOA and AOB in wet sediment ranged from 1.32×105to 1.91×106copies/g and from 5.39×105to 8.3×106copies/g, respectively. The sediment ammonia-oxidizing prokaryotes in downstream of the dam showed the highest abundance. Phylogenetic analysis revealed that most of the cloned amoA sequences of both AOA and AOB belonged to soil/sediment group. A considerable proportion of AOB sequences were highly affiliated with soil Nitrosospira lineage. Because of the influence from wastewater treatment plant effluent, some amoA sequences were similar to those recovered from activated sludge and wastewater.

    sediment;amoA;community composition;abundance;phylogeny

    X172

    A

    1000-6923(2015)01-0179-11

    鮑林林(1990-),女,四川綿陽人,首都師范大學(xué)資源環(huán)境與旅游學(xué)院碩士研究生,主要從事生態(tài)環(huán)境治理與建設(shè)研究.發(fā)表論文8篇.

    2014-04-28

    國家自然科學(xué)基金項(xiàng)目(41271495);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金聯(lián)合資助項(xiàng)目(20121108110006)

    * 責(zé)任作者, 教授, cnuwxy@sohu.com

    猜你喜歡
    北運(yùn)河古菌拷貝數(shù)
    不同pH和氧氣條件下土壤古菌與海洋古菌的競爭適應(yīng)機(jī)制*
    變油為氣,“榨干”廢棄油田
    大自然探索(2022年5期)2022-07-11 03:10:33
    海洋古菌
    線粒體DNA拷貝數(shù)變異機(jī)制及疾病預(yù)測價值分析
    北運(yùn)河流域水系綜合治理面臨的問題及建議
    胎兒染色體組拷貝數(shù)變異與產(chǎn)前超聲異常的相關(guān)性分析
    北運(yùn)河水環(huán)境保潔設(shè)備及配套設(shè)施研究
    DNA序列拷貝數(shù)變化決定黃瓜性別
    東太平洋海隆深海熱液區(qū)沉積物古菌多樣性分析
    線粒體DNA拷貝數(shù)的研究新進(jìn)展
    好男人在线观看高清免费视频| 九九久久精品国产亚洲av麻豆| av.在线天堂| 黄色丝袜av网址大全| 久久精品91蜜桃| 少妇的逼水好多| 国产又黄又爽又无遮挡在线| 俺也久久电影网| 99久久无色码亚洲精品果冻| 亚洲性夜色夜夜综合| 亚洲精品粉嫩美女一区| 亚洲18禁久久av| av国产免费在线观看| 亚洲自拍偷在线| 深夜a级毛片| 此物有八面人人有两片| 嫩草影视91久久| 国产乱人视频| 亚洲精品亚洲一区二区| 少妇高潮的动态图| 黄色配什么色好看| 两性午夜刺激爽爽歪歪视频在线观看| 黄色欧美视频在线观看| 日本三级黄在线观看| 91精品国产九色| 国内精品宾馆在线| 在现免费观看毛片| 亚洲精华国产精华精| 亚洲精华国产精华精| 欧美zozozo另类| 精品午夜福利在线看| 亚洲精品一卡2卡三卡4卡5卡| 特大巨黑吊av在线直播| 深爱激情五月婷婷| 三级国产精品欧美在线观看| a级一级毛片免费在线观看| 午夜爱爱视频在线播放| 国产精品一及| av福利片在线观看| 校园春色视频在线观看| 免费大片18禁| 国产精品无大码| 国产午夜精品论理片| 男人狂女人下面高潮的视频| 深爱激情五月婷婷| 久久婷婷人人爽人人干人人爱| av天堂中文字幕网| 国产熟女欧美一区二区| 一区二区三区激情视频| 亚洲电影在线观看av| 午夜福利在线观看免费完整高清在 | 一进一出好大好爽视频| 日日撸夜夜添| 国产欧美日韩精品亚洲av| 男女边吃奶边做爰视频| 中文字幕熟女人妻在线| 国产精品自产拍在线观看55亚洲| 69av精品久久久久久| 午夜免费激情av| 99在线视频只有这里精品首页| 禁无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 一级黄片播放器| 如何舔出高潮| 又爽又黄a免费视频| 成人美女网站在线观看视频| 日韩一区二区视频免费看| 日韩欧美在线乱码| 久久久久久久久久成人| 久久精品国产亚洲av涩爱 | 国产高清有码在线观看视频| 岛国在线免费视频观看| 性插视频无遮挡在线免费观看| 亚洲七黄色美女视频| 日韩av在线大香蕉| 一本久久中文字幕| 久久国内精品自在自线图片| av女优亚洲男人天堂| 午夜激情福利司机影院| 两人在一起打扑克的视频| 久久久久性生活片| 人妻少妇偷人精品九色| 熟女人妻精品中文字幕| 夜夜夜夜夜久久久久| 啦啦啦啦在线视频资源| 婷婷丁香在线五月| 内地一区二区视频在线| 亚洲男人的天堂狠狠| 九色国产91popny在线| 久久久久免费精品人妻一区二区| av福利片在线观看| 亚洲成人久久性| 午夜福利高清视频| 自拍偷自拍亚洲精品老妇| 日韩欧美精品v在线| 啪啪无遮挡十八禁网站| videossex国产| 男插女下体视频免费在线播放| 成年人黄色毛片网站| 成人av在线播放网站| 简卡轻食公司| 干丝袜人妻中文字幕| 亚洲图色成人| 国产欧美日韩精品一区二区| 91狼人影院| 一个人看视频在线观看www免费| 伊人久久精品亚洲午夜| 国产精品98久久久久久宅男小说| 亚洲成人免费电影在线观看| 在线观看av片永久免费下载| 国产午夜精品久久久久久一区二区三区 | 中文字幕精品亚洲无线码一区| 亚洲最大成人手机在线| 久久久久九九精品影院| 日韩 亚洲 欧美在线| 九九爱精品视频在线观看| 亚洲在线观看片| 亚洲无线观看免费| 国产熟女欧美一区二区| 久久久久久九九精品二区国产| www日本黄色视频网| 97超视频在线观看视频| 有码 亚洲区| 91在线精品国自产拍蜜月| 俺也久久电影网| 三级毛片av免费| 精品99又大又爽又粗少妇毛片 | 久久精品夜夜夜夜夜久久蜜豆| 国产精品一及| 久久久精品欧美日韩精品| 亚洲av一区综合| 俺也久久电影网| 国产69精品久久久久777片| 午夜久久久久精精品| 亚洲无线观看免费| 亚洲经典国产精华液单| 1000部很黄的大片| 午夜影院日韩av| 美女黄网站色视频| 国产精品女同一区二区软件 | 看十八女毛片水多多多| 亚洲一区二区三区色噜噜| 国产三级在线视频| 日本 av在线| 麻豆国产97在线/欧美| 91在线精品国自产拍蜜月| 精品欧美国产一区二区三| 99久国产av精品| 真人做人爱边吃奶动态| 亚洲,欧美,日韩| 亚洲人成网站在线播放欧美日韩| 国产爱豆传媒在线观看| 美女大奶头视频| 国产男人的电影天堂91| 观看美女的网站| 精品一区二区三区视频在线| 在线看三级毛片| 精品久久久久久久久久免费视频| 久久久久免费精品人妻一区二区| 国内精品美女久久久久久| av国产免费在线观看| 国产亚洲91精品色在线| 国产伦人伦偷精品视频| 久久午夜亚洲精品久久| 国产精品久久久久久久久免| 九色国产91popny在线| 中文资源天堂在线| 少妇高潮的动态图| 久久精品人妻少妇| 国产私拍福利视频在线观看| 国产精品精品国产色婷婷| 国内精品宾馆在线| 91精品国产九色| 亚洲在线自拍视频| 成人国产综合亚洲| 欧美+亚洲+日韩+国产| 草草在线视频免费看| 久久精品影院6| 日本 欧美在线| 欧美丝袜亚洲另类 | 国产av在哪里看| 亚洲一级一片aⅴ在线观看| 国产国拍精品亚洲av在线观看| 69人妻影院| 99riav亚洲国产免费| 久久久久性生活片| 又黄又爽又刺激的免费视频.| 欧美激情国产日韩精品一区| 悠悠久久av| 少妇丰满av| 床上黄色一级片| 中国美女看黄片| 麻豆精品久久久久久蜜桃| 国产大屁股一区二区在线视频| 国产伦在线观看视频一区| 不卡一级毛片| 婷婷六月久久综合丁香| 俄罗斯特黄特色一大片| 最近在线观看免费完整版| 国产v大片淫在线免费观看| 真人做人爱边吃奶动态| 国产精品电影一区二区三区| 网址你懂的国产日韩在线| 日本色播在线视频| 国产91精品成人一区二区三区| 噜噜噜噜噜久久久久久91| 最后的刺客免费高清国语| 又爽又黄无遮挡网站| 看黄色毛片网站| 亚洲一区二区三区色噜噜| 久久99热这里只有精品18| 欧美性猛交黑人性爽| 99久久精品国产国产毛片| 特大巨黑吊av在线直播| 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 成人国产麻豆网| 欧美日韩中文字幕国产精品一区二区三区| 最近中文字幕高清免费大全6 | 国产aⅴ精品一区二区三区波| 国产探花极品一区二区| 亚洲精品亚洲一区二区| 精品一区二区免费观看| 国内毛片毛片毛片毛片毛片| 精品人妻熟女av久视频| 国产在线精品亚洲第一网站| av福利片在线观看| 国产不卡一卡二| 校园人妻丝袜中文字幕| 久久精品国产亚洲网站| 亚洲三级黄色毛片| 国产视频一区二区在线看| 亚洲成人中文字幕在线播放| 亚洲无线在线观看| 99精品久久久久人妻精品| 精品久久久噜噜| 国产高清不卡午夜福利| 97人妻精品一区二区三区麻豆| 九九热线精品视视频播放| 久9热在线精品视频| 欧洲精品卡2卡3卡4卡5卡区| 老司机深夜福利视频在线观看| 日本与韩国留学比较| 国产av不卡久久| 给我免费播放毛片高清在线观看| 91麻豆av在线| 欧美激情国产日韩精品一区| 日韩一本色道免费dvd| 亚洲av美国av| 最后的刺客免费高清国语| 免费黄网站久久成人精品| 日韩国内少妇激情av| 国产精品无大码| 一个人免费在线观看电影| 美女 人体艺术 gogo| 男女那种视频在线观看| 免费观看在线日韩| 日本 av在线| 两个人视频免费观看高清| 婷婷丁香在线五月| 欧美激情国产日韩精品一区| 国产真实伦视频高清在线观看 | 99热精品在线国产| 伊人久久精品亚洲午夜| .国产精品久久| 美女大奶头视频| 亚洲av一区综合| 最后的刺客免费高清国语| 国产成人av教育| 免费黄网站久久成人精品| 国产色婷婷99| 成年女人看的毛片在线观看| 国产91精品成人一区二区三区| 毛片一级片免费看久久久久 | 国产探花极品一区二区| 深夜a级毛片| 日日干狠狠操夜夜爽| 老司机午夜福利在线观看视频| 尾随美女入室| 国产av在哪里看| 超碰av人人做人人爽久久| 久久99热6这里只有精品| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 18禁在线播放成人免费| 99热这里只有是精品在线观看| 午夜免费男女啪啪视频观看 | 三级男女做爰猛烈吃奶摸视频| 色综合亚洲欧美另类图片| 久久亚洲真实| 又爽又黄无遮挡网站| 内射极品少妇av片p| 日韩欧美精品v在线| 最近中文字幕高清免费大全6 | 国产毛片a区久久久久| 成人国产麻豆网| 一区二区三区免费毛片| 99久久成人亚洲精品观看| 国产亚洲精品久久久久久毛片| 精品一区二区免费观看| 亚洲精华国产精华精| 老师上课跳d突然被开到最大视频| 亚洲国产高清在线一区二区三| 亚洲一级一片aⅴ在线观看| 男人的好看免费观看在线视频| 国产精品一区二区免费欧美| 麻豆成人av在线观看| 国产精品久久久久久精品电影| 春色校园在线视频观看| 国产色婷婷99| 一进一出抽搐gif免费好疼| 亚洲av日韩精品久久久久久密| 99热这里只有是精品在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产真实伦视频高清在线观看 | 亚洲图色成人| 欧美日韩综合久久久久久 | 一本久久中文字幕| 最近最新中文字幕大全电影3| 九九在线视频观看精品| 国产男人的电影天堂91| 非洲黑人性xxxx精品又粗又长| 国产精品人妻久久久久久| 中文字幕av在线有码专区| 午夜福利在线观看免费完整高清在 | 99九九线精品视频在线观看视频| 又爽又黄无遮挡网站| 一个人看的www免费观看视频| 亚州av有码| 最近在线观看免费完整版| 波多野结衣巨乳人妻| 国产三级在线视频| 免费搜索国产男女视频| 国产又黄又爽又无遮挡在线| 色哟哟哟哟哟哟| 国产成人福利小说| 99热6这里只有精品| 伊人久久精品亚洲午夜| www日本黄色视频网| 日韩高清综合在线| 国产精品电影一区二区三区| 午夜福利18| 中亚洲国语对白在线视频| 久久久久久久亚洲中文字幕| ponron亚洲| 免费观看人在逋| 中文字幕免费在线视频6| 女同久久另类99精品国产91| 亚洲国产欧美人成| 久久久国产成人精品二区| 国产黄片美女视频| 美女高潮喷水抽搐中文字幕| 最近在线观看免费完整版| 午夜免费男女啪啪视频观看 | 亚洲精品日韩av片在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩东京热| 国产私拍福利视频在线观看| 国产成人aa在线观看| 欧美日韩瑟瑟在线播放| 99热网站在线观看| 在现免费观看毛片| 国产午夜精品论理片| 国产免费av片在线观看野外av| 禁无遮挡网站| 久久久久久久久中文| 一a级毛片在线观看| 国内精品宾馆在线| 在线观看舔阴道视频| av国产免费在线观看| 男人的好看免费观看在线视频| 18+在线观看网站| 午夜福利欧美成人| 麻豆成人午夜福利视频| 国产伦在线观看视频一区| 国产精品久久久久久精品电影| 一进一出抽搐动态| 老司机深夜福利视频在线观看| 日韩欧美三级三区| 色视频www国产| 久久香蕉精品热| 国产精品自产拍在线观看55亚洲| 午夜福利成人在线免费观看| 成年女人看的毛片在线观看| 观看美女的网站| 又黄又爽又刺激的免费视频.| 国产男靠女视频免费网站| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| 久久国产乱子免费精品| a在线观看视频网站| 波多野结衣巨乳人妻| 精品国产三级普通话版| 别揉我奶头~嗯~啊~动态视频| 伦精品一区二区三区| 99热网站在线观看| 少妇人妻精品综合一区二区 | 97超视频在线观看视频| 成人av在线播放网站| 亚洲在线自拍视频| 女同久久另类99精品国产91| 国产精品福利在线免费观看| 精品乱码久久久久久99久播| 亚洲精品影视一区二区三区av| 性欧美人与动物交配| 亚洲国产精品合色在线| 97碰自拍视频| 亚洲国产欧洲综合997久久,| 日日干狠狠操夜夜爽| 三级毛片av免费| 国产综合懂色| 特级一级黄色大片| 日韩精品有码人妻一区| 长腿黑丝高跟| 久久精品国产99精品国产亚洲性色| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 在线播放国产精品三级| 97碰自拍视频| 天美传媒精品一区二区| 99热这里只有是精品50| 亚洲成人精品中文字幕电影| 免费观看的影片在线观看| 欧美+日韩+精品| 88av欧美| 成人亚洲精品av一区二区| 我的老师免费观看完整版| 成人午夜高清在线视频| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| aaaaa片日本免费| 国产视频一区二区在线看| 嫩草影院精品99| 在线播放国产精品三级| 人人妻人人澡欧美一区二区| 日本欧美国产在线视频| 亚洲最大成人中文| 亚洲精品粉嫩美女一区| 精品99又大又爽又粗少妇毛片 | 欧美成人a在线观看| 黄色欧美视频在线观看| 精品久久久久久久久亚洲 | 中文亚洲av片在线观看爽| 日韩欧美免费精品| 欧美一级a爱片免费观看看| 午夜a级毛片| 亚洲欧美日韩高清专用| 免费无遮挡裸体视频| 九九热线精品视视频播放| 又爽又黄a免费视频| 亚洲自拍偷在线| 国产精品一区www在线观看 | 制服丝袜大香蕉在线| 国产在线精品亚洲第一网站| 村上凉子中文字幕在线| 国产精品嫩草影院av在线观看 | 内射极品少妇av片p| avwww免费| 热99re8久久精品国产| 欧美区成人在线视频| 日韩欧美三级三区| 蜜桃亚洲精品一区二区三区| 在线观看av片永久免费下载| 在线观看午夜福利视频| 成年女人毛片免费观看观看9| 精品一区二区三区视频在线| av在线亚洲专区| 国产成人福利小说| 久久精品国产亚洲av天美| 午夜影院日韩av| 国产人妻一区二区三区在| 黄色日韩在线| 亚洲av一区综合| 波多野结衣巨乳人妻| 99热只有精品国产| 亚洲人成网站在线播| 午夜影院日韩av| 国产真实乱freesex| 久久香蕉精品热| 免费人成视频x8x8入口观看| 午夜视频国产福利| 国产精品一区二区三区四区久久| 人妻久久中文字幕网| 波多野结衣高清无吗| 久久久国产成人免费| 免费av不卡在线播放| 国产黄色小视频在线观看| 成人精品一区二区免费| 99久久精品热视频| 男人狂女人下面高潮的视频| 国产极品精品免费视频能看的| 国产一区二区激情短视频| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉精品热| 琪琪午夜伦伦电影理论片6080| 欧美极品一区二区三区四区| 婷婷精品国产亚洲av在线| 1024手机看黄色片| 国产老妇女一区| 性插视频无遮挡在线免费观看| 搞女人的毛片| 日韩强制内射视频| 天堂av国产一区二区熟女人妻| 亚洲国产欧洲综合997久久,| 黄色丝袜av网址大全| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 亚洲真实伦在线观看| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 久久久色成人| 日本黄色视频三级网站网址| 亚洲五月天丁香| 99久久九九国产精品国产免费| 国产成人福利小说| 免费观看精品视频网站| 国产亚洲欧美98| 免费看日本二区| 在线观看午夜福利视频| 精品日产1卡2卡| 黄色一级大片看看| 欧美日韩瑟瑟在线播放| 成人三级黄色视频| 少妇猛男粗大的猛烈进出视频 | 久久精品综合一区二区三区| 美女 人体艺术 gogo| 91在线观看av| 久久精品国产亚洲av香蕉五月| 嫩草影院新地址| 一本精品99久久精品77| 亚洲最大成人中文| 欧美日韩精品成人综合77777| av视频在线观看入口| 国产美女午夜福利| 真实男女啪啪啪动态图| 精品人妻偷拍中文字幕| 欧美性猛交╳xxx乱大交人| 乱码一卡2卡4卡精品| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 亚洲av电影不卡..在线观看| ponron亚洲| 99国产精品一区二区蜜桃av| 中文字幕av在线有码专区| 久久这里只有精品中国| 18禁黄网站禁片午夜丰满| av福利片在线观看| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国内精品美女久久久久久| 18禁在线播放成人免费| 国产综合懂色| 久久人人精品亚洲av| 久久午夜亚洲精品久久| 特级一级黄色大片| 免费高清视频大片| 婷婷丁香在线五月| 国产高清视频在线播放一区| 丝袜美腿在线中文| 少妇的逼水好多| 国产蜜桃级精品一区二区三区| 在线免费观看不下载黄p国产 | 国产男人的电影天堂91| 国产精品久久久久久久久免| 午夜福利高清视频| 少妇猛男粗大的猛烈进出视频 | x7x7x7水蜜桃| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 91麻豆av在线| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 特级一级黄色大片| 男女边吃奶边做爰视频| 亚洲精品456在线播放app | 久久6这里有精品| 少妇的逼水好多| 长腿黑丝高跟| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 国产麻豆成人av免费视频| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| 国产精品一区二区免费欧美| 美女大奶头视频| 狂野欧美白嫩少妇大欣赏| 一区二区三区激情视频| 国产黄色小视频在线观看| 亚洲av中文字字幕乱码综合| 色综合亚洲欧美另类图片| 欧美色欧美亚洲另类二区| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 亚洲成av人片在线播放无| 亚洲av成人av| 亚洲天堂国产精品一区在线| 国产真实伦视频高清在线观看 | 亚洲,欧美,日韩| 久久久久久大精品| 大又大粗又爽又黄少妇毛片口| 国产伦精品一区二区三区视频9| 亚洲欧美激情综合另类| av在线蜜桃| 男人和女人高潮做爰伦理| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 亚洲成av人片在线播放无| 成人特级av手机在线观看| 国产高清三级在线| 久久亚洲精品不卡| 亚洲一区高清亚洲精品| 国产高清视频在线播放一区| 国产亚洲av嫩草精品影院| av天堂中文字幕网| 日本a在线网址|