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

    東湖表層沉積物中氨氧化古菌和氨氧化細(xì)菌豐度及多樣性研究

    2018-03-30 03:33:09潘彥羽代嫣然王飛華
    水生生物學(xué)報(bào) 2018年2期
    關(guān)鍵詞:東湖沉積物氨氮

    潘彥羽 代嫣然 王飛華 梁 威

    (1. 中國(guó)科學(xué)院水生生物研究所淡水生態(tài)與生物技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 武漢 430072; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049)

    隨著人口的增長(zhǎng)以及經(jīng)濟(jì)的快速發(fā)展, 大量氮磷等營(yíng)養(yǎng)元素及有機(jī)污染物未經(jīng)處理即匯入淡水生態(tài)系統(tǒng)[1], 致使其生態(tài)系統(tǒng)的生態(tài)服務(wù)功能日漸退化甚至消失。2015年《中國(guó)環(huán)境統(tǒng)計(jì)年報(bào)》顯示全國(guó)廢水中氨氮排放量2.3×109kg, 其中城鎮(zhèn)生活污水中氨氮排放量占氨氮排放總量的58.3%, 農(nóng)業(yè)源氨氮排放量占排放總量的31.6%, 工業(yè)廢水氨氮排放量占氨氮排放總量的9.4%[2]。水體中氮素含量過高不僅會(huì)極大促進(jìn)藻類的生長(zhǎng)繁殖, 導(dǎo)致水體富營(yíng)養(yǎng)化, 還會(huì)毒害其中的水生動(dòng)物[3,4]。

    氮循環(huán)具有復(fù)雜的變化過程, 并且與微生物多樣性和群落結(jié)構(gòu)密切相關(guān)。這個(gè)過程包括: 固氮作用、氨化作用、硝化作用和反硝化作用[5]。其中氨氧化是硝化作用的第一步, 也是限速步驟[6]。長(zhǎng)期以來, 氨氧化細(xì)菌(Ammonia-oxidizing bacterial,AOB) 被認(rèn)為是氨氧化作用的主要承擔(dān)者[7], 然而,近年來隨著分子生物學(xué)技術(shù)的發(fā)展和研究的深入,新一類具有氨氧化能力的微生物——氨氧化古菌(Ammonia-oxidizing archaea, AOA) 從環(huán)境中被分離和鑒定出來。2005年Konneke等[8]從海洋沉積物中分離了一株氨氧化古菌。2006年Leininger等[9]在歐洲中性土壤中發(fā)現(xiàn)了AOA的廣泛存在, 并且amoA基因豐度大于AOB。氨氧化過程中兩個(gè)關(guān)鍵酶分別是氨單加氧酶和羥胺氧化還原酶。氨單加氧酶是由amoA、amoB、amoC三個(gè)基因編碼的α、β、γ三個(gè)亞基組成。目前, 編碼α亞基的amoA基因被廣泛用于生態(tài)系統(tǒng)中AOA和AOB的研究[10]。大量研究表明, AOA和AOB存在于很多生態(tài)系統(tǒng)中,如土壤[11]、海洋沉積物[12]、湖泊沉積物[13]和活性污泥[14]等, 但是AOA和AOB在淡水生態(tài)系統(tǒng)尤其是沉積物中的生態(tài)特征、功能差異及多樣化環(huán)境要素的影響機(jī)制仍需要進(jìn)一步研究[15]。

    本研究以典型城市湖泊——東湖為例, 通過監(jiān)測(cè)東湖表層沉積物的理化參數(shù), 采用熒光定量PCR和高通量測(cè)序等分子生物學(xué)技術(shù), 以功能基因amoA為分子標(biāo)記, 對(duì)AOA和AOB進(jìn)行定量和定性分析。探討在不同環(huán)境條件下表層沉積物中AOA和AOB的amoA基因的數(shù)量及多樣性, 分析底質(zhì)和水質(zhì)中各種理化因子對(duì)它的影響, 力求為湖泊氮循環(huán)研究提供理論基礎(chǔ)。

    1 材料與方法

    1.1 采樣地點(diǎn)

    本研究選取東湖子湖——郭鄭湖作為研究對(duì)象, 設(shè)置3個(gè)采樣點(diǎn)(圖 1), 1號(hào)點(diǎn)(30°33′4.71″N 114°21′36.68″E); 2號(hào)點(diǎn)(30°33′49.02″N 114°23′3.52″E);3號(hào)點(diǎn)(30°33′6.20″N 114°23′31.73″E)。

    1.2 樣品采集

    分別于2016年1月、4月、7月和10月采集東湖表層沉積物樣品(12個(gè)樣品, 0—20 cm)。將沉積物樣品充分混勻后分兩份置于無菌塑料袋中密封,4℃保存運(yùn)回實(shí)驗(yàn)室, 一份自然風(fēng)干過200目篩用于理化性質(zhì)的測(cè)定, 另一份放于-20℃冰凍保存用于后續(xù)分子實(shí)驗(yàn)。

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

    采用標(biāo)準(zhǔn)方法對(duì)樣品進(jìn)行分析, 其中pH采用pH計(jì)(Mettler Toledo, Switzerland)測(cè)定; 氨氮(-N)、亞硝酸氮(-N)、硝酸鹽氮(-N)采用氯化鉀溶液提取-分光光度法測(cè)定; 底泥總氮(TN)、總磷(TP)測(cè)定按照《土壤農(nóng)化分析》(第三版)所述方法進(jìn)行[16]; 總有機(jī)碳(TOC)采用總有機(jī)碳分析儀(Elementar, Germany)測(cè)定。每個(gè)樣品設(shè)置3個(gè)平行。

    1.4 沉積物中AOA和AOB豐度的測(cè)定

    沉積物樣品DNA提取稱取0.3 g沉積物樣品, 采用PowerSoil DNA試劑盒(MO BIO Laboratories Inc., USA)進(jìn)行樣品DNA提取, 提取產(chǎn)物用1%的瓊脂凝膠電泳檢測(cè)。

    圖 1 東湖采樣點(diǎn)Fig. 1 Sampling sites in Donghu Lake

    絕對(duì)定量標(biāo)準(zhǔn)曲線制作參照He等[17]報(bào)道的方法, 將提取的土壤樣品的微生物全基因組DNA混合并作為PCR擴(kuò)增的模板, 擴(kuò)增引物和程序如表 1所示。擴(kuò)增體系為: 10×PCR buffer (Mg2+)5 μL, dNTP (2.5 mmol/L) 4 μL, 正反向引物(10 mmol/L)各1 μL, Taq酶(2.5U)0.25 μL, DNA模板2 μL, 加ddH2O補(bǔ)足至50 μL。使用Easy Pure QuickGel Extraction Kit (Trans)純化試劑盒對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行切膠純化, T載體與PCR回收產(chǎn)物連接后轉(zhuǎn)入JM109感受態(tài)細(xì)胞, 進(jìn)行藍(lán)白斑篩選。選取白色克隆, 采用菌落PCR鑒定陽性克隆, 所用引物為載體通用引物M13F和M13R。將挑取出的重組質(zhì)粒進(jìn)行測(cè)序。測(cè)序結(jié)果經(jīng)NCBI比對(duì), AOA與泉古菌(JQ345886.1) amoA基因的同源性高達(dá)99%, AOB與亞硝化單胞菌(KC769053.1) amoA基因的同源性高達(dá)99%。因此, 可作為絕對(duì)熒光定量PCR分析的標(biāo)準(zhǔn)DNA。使用EasyPure Plasmid Miniprep Kit(Trans)提取重組質(zhì)粒DNA, 經(jīng)核酸定量?jī)x(Nanodrop 2000, USA) 測(cè)定濃度, AOA的amoA基因重組質(zhì)粒濃度為44.6 ng/μL, AOB的amoA基因重組質(zhì)粒濃度為37.5 ng/μL。根據(jù)阿伏伽德羅常數(shù)計(jì)算出amoA基因的拷貝數(shù), AOA為2.11×1010copies/μL,AOB為1.77×1010copies/μL。將AOA重組質(zhì)粒按10倍梯度稀釋, 用于標(biāo)準(zhǔn)曲線繪制。

    熒光定量PCR采用SYBR Green法進(jìn)行熒光定量PCR, 所用引物參見表 1, 擴(kuò)增體系為20 μL,其中, 2×Trans Start Tip Green qPCRSuperMix 10 μL,濃度為10 μmol/L的AOA、AOB正反向引物各0.4 μL,DNA模板1 μL, 用ddH2O補(bǔ)足至20 μL。采用Roche LightCycler? 480擴(kuò)增儀進(jìn)行定量, 其反應(yīng)程序如表 1所示。實(shí)驗(yàn)設(shè)置陰性對(duì)照, 并與稀釋好的6個(gè)標(biāo)準(zhǔn)樣品進(jìn)行定量擴(kuò)增, 得到標(biāo)準(zhǔn)曲線, 每個(gè)樣品設(shè)置3次重復(fù)。以基線(背景)熒光信號(hào)標(biāo)準(zhǔn)差的10倍作為閾值, 擴(kuò)增效率>80%, 溶解曲線為單一峰。

    1.5 沉積物中AOA和AOB群落結(jié)構(gòu)分析

    HiSeq平臺(tái)高通量測(cè)序針對(duì)氨氧化微生物的功能基因amoA, 采用表 1中的引物進(jìn)行擴(kuò)增。PCR產(chǎn)物用2%瓊脂糖凝膠電泳檢測(cè), 引物擴(kuò)增出來的AOA amoA基因片段長(zhǎng)度在600—700 bp, AOB amoA基因片段長(zhǎng)度在400—500 bp, 使用AxyPrep-DNA凝膠回收試劑盒切膠回收PCR產(chǎn)物, TrisHCl洗脫; 2%瓊脂糖電泳檢測(cè)。參照電泳初步定量結(jié)果, 將PCR產(chǎn)物用QuantiFluoruantTM藍(lán)色熒光定量系統(tǒng)進(jìn)行檢測(cè)定量, 之后按照每個(gè)樣本的測(cè)序量要求, 進(jìn)行相應(yīng)比例的混合。通過Hiseq 2500平臺(tái)(IlluminaSanDiego, USA)進(jìn)行測(cè)序, 測(cè)序長(zhǎng)度250 bp(凌恩生物, 上海)。

    表 1 氨氧化微生物amoA基因PCR擴(kuò)增引物及反應(yīng)條件Tab. 1 Primers and PCR reaction conditions for amoA gene

    AOA和AOB amoA基因多樣性分析利用Cluster及Mothur軟件對(duì)序列進(jìn)行降噪, 并進(jìn)行聚類,獲得不同序列之間的距離矩陣, 轉(zhuǎn)化為相似度, 將相似度超過80%的序列歸為一個(gè)操作分類單元(Operational taxonomic unit, OTU)。通過稀釋曲線比較不同樣品細(xì)菌群落多樣性差異。用Mothur計(jì)算微生物群落的多樣性和豐度指數(shù), 包括Shannon指數(shù)、Simpson指數(shù)、ACE指數(shù)、Chao指數(shù)和文庫覆蓋率C[18]。

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

    挑選豐度最高的前40個(gè)OTUs的代表核酸序列與數(shù)據(jù)庫中已知序列進(jìn)行比對(duì)。采用Neighbor-Joining建樹方法利用MEGA6.0軟件構(gòu)建AOA和AOB的amoA基因系統(tǒng)發(fā)育樹, 建樹結(jié)果進(jìn)行1000次Bootstrap檢驗(yàn)。

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

    用SPSS 16.0統(tǒng)計(jì)分析軟件進(jìn)行方差分析及多因素相關(guān)分析。amoA基因豐度進(jìn)行了以10為底的對(duì)數(shù)轉(zhuǎn)換。氨氧化微生物群落結(jié)構(gòu)與環(huán)境因子之間的關(guān)系采用Canoco軟件進(jìn)行典范對(duì)應(yīng)分析, 確定顯著影響AOA或AOB群落結(jié)構(gòu)的環(huán)境因子, 顯著水平設(shè)置為P<0.05和極其顯著P<0.01。

    2 結(jié)果

    2.1 沉積物理化指標(biāo)

    東湖沉積物基本理化特征如表 2所示。溫度隨季節(jié)變化明顯, 其中夏季最高(28.6—29.3℃), 冬季溫度最低(6.91—7.14℃)。沉積物呈弱堿性(7.0—7.5), 1號(hào)位至3號(hào)位點(diǎn), 沉積物pH呈現(xiàn)逐漸降低的趨勢(shì)。沉積物TN含量介于1.18—2.31 g/kg (干重), 其中冬季、夏季總氮含量較低, 春季、秋季TN含量較高。沉積物總磷含量介于0.58—1.15 g/kg(干重); 沉積物-N濃度介于26.8—60.9 mg/kg(干重);-N (-N+-N)濃度介于11.1—50.4 mg/kg (干重)。TOC濃度介于27.1—49.2 g/kg(干重), 其中2號(hào)點(diǎn)TP、-N及TOC含量均較高,表明2號(hào)位點(diǎn)污染最為嚴(yán)重, 可能是由于2013年對(duì)1號(hào)點(diǎn)區(qū)域進(jìn)行了清淤使得1號(hào)點(diǎn)受污染程度有所好轉(zhuǎn)。

    2.2 沉積物氨氧化微生物的豐度

    AOA和AOB的定量分析結(jié)果顯示(圖 2), 所有樣品中均檢測(cè)到AOA和AOB的存在, 且AOA的豐度高于AOB, AOA與AOB的比值介于1.10—-73.6。AOA的豐度介于1.50×108—8.48×109拷貝/g沉積物(干重), 其中冬季AOA較多, 介于2.83×109—8.48×109拷貝/g沉積物(干重); 夏季最低為1.74×108—1.95×109拷貝/g沉積物(干重)。AOB的豐度介于4.87×107—2.40×109拷貝/g沉積物(干重),4個(gè)季節(jié)AOB豐度差異不顯著。

    表 2 沉積物理化性質(zhì)Tab. 2 Physical and chemical properties in sediments

    將東湖AOA和AOB的amoA基因豐度與沉積物的理化因子做Pearson相關(guān)性分析(表 3), 結(jié)果表明AOA的豐度與溫度呈顯著負(fù)相關(guān)(P<0.05), 沉積物溫度越高, AOA的豐度越少; AOB豐度與理化因子相關(guān)性不顯著(P>0.05)。

    2.3 沉積物氨氧化微生物多樣性

    對(duì)Shannon指數(shù)、Simpson指數(shù)、OTU數(shù)目的Chao指數(shù)和ACE指數(shù)進(jìn)行了分析, 具體結(jié)果見表4—6。利用Mothur進(jìn)行數(shù)據(jù)的過濾分析, 12個(gè)沉積物樣品中共得到高質(zhì)量的AOA amoA基因序列906010條, 平均每個(gè)樣品的序列條數(shù)為37750條。在所有的12個(gè)沉積物樣品中, 按照80%的序列相似度進(jìn)行OTU劃分。東湖沉積物樣品的OTU數(shù)量介于5—16, Shannon指數(shù)介于0.03—0.05。比較東湖不同季節(jié)不同采樣的多樣性指數(shù)發(fā)現(xiàn), 冬季樣品多樣性高于其他季節(jié), 2號(hào)位點(diǎn)多樣性最低。

    圖 2 沉積物中AOA和AOB的豐度Fig. 2 The abundance of AOA and AOB in sediments

    表 3 AOA和AOB豐度與沉積物理化指標(biāo)相關(guān)性分析Tab. 3 Correlation analysis between the abundance of AOA and AOB and physical and chemical indicators in sediments

    同時(shí), 我們對(duì)沉積物AOB amoA基因進(jìn)行測(cè)序(表 5), AOB的文庫覆蓋率均達(dá)到了較高的水平(>99.9%), 這表明了所建立的amoA基因的文庫涵蓋了沉積物中幾乎所有的amoA基因類型。不同采樣點(diǎn)的樣品OTU數(shù)量差異顯著, 沉積物樣品AOB的OTU數(shù)量介于42—102。Shannon指數(shù)最高的是DSP1(1.99), 最低的是DAU2(1.34)。利用不同的指數(shù)來估算樣本中微生物多樣性時(shí), 各指數(shù)反應(yīng)的結(jié)果略有偏差, 本研究中DSP3 Chao指數(shù)最高, 但其Shannon指數(shù)為1.87, 最大為1.99??赡苁歉髦笖?shù)的計(jì)算方法不同, 并且群落多樣性受群落組成和群落分布情況的共同影響。

    對(duì)氨氧化微生物的多樣性指標(biāo)與沉積物的理化性質(zhì)進(jìn)行相關(guān)性分析(表 6), 結(jié)果發(fā)現(xiàn)東湖沉積物中AOA OTU數(shù)量與TOC呈極顯著負(fù)相關(guān)(P<0.01), 而AOB OTU數(shù)量與pH呈顯著負(fù)相關(guān)(P<0.05)。此外, AOA Shannon指數(shù)和AOB Shannon指數(shù)也與TOC呈極顯著負(fù)相關(guān)(P<0.01), 說明沉積物TOC越高, 群落多樣性越低。2號(hào)位點(diǎn)相對(duì)于其他位點(diǎn)TOC含量最高, 所以AOA和AOB的多樣性均較低。

    2.4 沉積物氨氧化微生物的群落結(jié)構(gòu)分析

    高通量測(cè)序解析氨氧化微生物群落的組成和多樣性結(jié)果如圖 3—5所示??傮w來看東湖沉積物中AOA親緣關(guān)系較近的既有來自水體的, 也有來自土壤和沉積物的, AOA序列可以分為Nitrosopumilus、Nitrososphaera以及其他。各個(gè)季節(jié)各個(gè)采樣點(diǎn)之間AOA群落組成有一定的差異, 但是Nitrosopumilus均占絕對(duì)優(yōu)勢(shì)(99.8%), 僅在1號(hào)位點(diǎn)的冬季、春季及夏季檢測(cè)到了Nitrosospharea存在, 含量很低。AOB amoA基因序列分為來自土壤和淡水環(huán)境的Nitrosomonas, Nitrosospira以及相當(dāng)一部分處于發(fā)育樹外枝的序列。在東湖不同季節(jié)不同采樣點(diǎn)采集的12個(gè)樣品中, Nitrosomonas占55.27%, 其他占44.70%, Nitrosospira含量?jī)H占0.03%。

    2.5 環(huán)境因子對(duì)AOA和AOB群落結(jié)構(gòu)的典范對(duì)應(yīng)分析

    利用Canoco將AOA和AOB的優(yōu)勢(shì)種與環(huán)境因子進(jìn)行典范對(duì)應(yīng)(CCA)分析(圖 6)。東湖AOA群落結(jié)構(gòu)與環(huán)境因子排序軸AX1和AX2的解釋度分別為80.6%和9.14%。TOC (-0.9999)、TP (-0.9872)與AX1的相關(guān)性較高;-N (0.6545)與AX2的相關(guān)性較高。AOB優(yōu)勢(shì)種與環(huán)境因子的結(jié)果表明, 排序軸AX1的解釋度為88.5%, 排序軸AX2的解釋度為2.70%。與AX1相關(guān)性高的環(huán)境因子有TP(-0.9984)、TN (-0.9798)、TOC (-0.9933), 與AX2軸相關(guān)性較高的環(huán)境因子為溫度(-0.9631)和pH (0.7244), 總體而言, TN是影響東湖沉積物AOA群落結(jié)構(gòu)的主要因素, TOC和TP是影響AOB群落結(jié)構(gòu)的主要因素。

    表 4 沉積物中AOA的多樣性指數(shù)Tab. 4 Diversity index of AOA in sediments

    表 5 沉積物中AOB的多樣性指數(shù)Tab. 5 Diversity index of AOB in Sediments

    表 6 AOA和AOB的多樣性指數(shù)與沉積物理化指標(biāo)相關(guān)性分析Tab. 6 Correlation analysis between the diversity index of AOA and AOB and physical and chemical indicators in sediments

    3 討論

    已有研究表明, AOA比AOB更適合厭氧或低氨氮濃度條件[19], 在海洋生態(tài)系統(tǒng)[20]和土壤環(huán)境中[21]AOA amoA基因比AOB的豐度高3個(gè)數(shù)量級(jí)。在淡水生態(tài)系統(tǒng)中, Hou等[22]分析了太湖和巢湖不同富營(yíng)養(yǎng)程度點(diǎn)位中沉積物AOA與AOB豐度, 發(fā)現(xiàn)在中營(yíng)養(yǎng)采樣點(diǎn)中AOA多于AOB; Yang等[23]發(fā)現(xiàn)在洱海的沉積物中, AOA數(shù)量多于AOB。這與本研究在東湖沉積物中AOA和AOB豐度的研究結(jié)果類似, 攜帶amoA基因的微生物并不足以說明其參與氨氧化過程, 因此無法僅僅以基因豐度含量高低確定AOA和AOB在氮循環(huán)中的具體貢獻(xiàn), 但目前還沒有發(fā)現(xiàn)不能進(jìn)行氨氧化過程的此類微生物[24],AOA可能在東湖生態(tài)系統(tǒng)氨氧化中占主要作用。

    湖泊沉積物中的AOA和AOB受多種環(huán)境因子的影響, 現(xiàn)有研究發(fā)現(xiàn)AOA豐度受pH[25]、-N[26]、-N[27]、TN[26]、TP[27]及TOC[28]的影響,AOB的豐度受-N[29]、-N[26]、pH[26]的影響。在本研究中, 東湖沉積物AOA豐度在不同季節(jié)存在顯著差異, 其數(shù)量與溫度呈顯著負(fù)相關(guān)。Lu等[30]研究農(nóng)業(yè)池塘沉積物AOA群落分布時(shí), 也得到相似的結(jié)論。Hou等[31]研究了東湖不同富營(yíng)養(yǎng)化程度子湖郭鄭湖和團(tuán)湖AOA群落分布, 發(fā)現(xiàn)AOA豐度與沉積物中有機(jī)物的含量呈負(fù)相關(guān), AOA豐度對(duì)環(huán)境因子變化敏感。本研究中考察了季節(jié)尺度, 方差分析表明不同季節(jié)之間TN、TP、-N均有顯著差異(P<0.05), 使得溫度對(duì)AOA數(shù)量的影響更為明顯。pH[32]被認(rèn)為是影響AOA和AOB amoA基因數(shù)量的一個(gè)重要環(huán)境因子, 但是在本研究中, pH對(duì)沉積物中AOA和AOB的amoA基因數(shù)量無顯著影響。這些不一致的研究結(jié)果表明AOA和AOB對(duì)于氨氧化過程中的貢獻(xiàn)與環(huán)境條件密切相關(guān)。

    利用高通量測(cè)序分析了東湖沉積物AOA和AOB的群落結(jié)構(gòu), 在AOA和AOB中均檢測(cè)到部分序列(其他)與目前GenBank中已經(jīng)探明的amoA序列相似性較低, 可能是未分離或培養(yǎng)的區(qū)別于已知菌種的新型氨氧化微生物。東湖沉積物AOA Shannon指數(shù)介于0.03—0.05, 東湖AOB Shannon指數(shù)介于1.34—1.99, 高于AOA, 在淡水濕地[33]、金山湖[27]及太湖[26]沉積物中AOB多樣性均較AOA高,AOA多樣性較AOB低, 使得優(yōu)勢(shì)屬的相對(duì)豐度和優(yōu)勢(shì)度更加明顯。李虎等[34]研究了浙江甌江沉積物中AOA及AOB組成, 發(fā)現(xiàn)AOA的shannon指數(shù)在2.64—2.89, AOB的shannon指數(shù)在2.83—3.56。多樣性顯著高于東湖, 推測(cè)是由于東湖是湖泊生態(tài)系統(tǒng), 環(huán)境條件相對(duì)封閉單一, 因此生物多樣性也較低。河流生態(tài)系統(tǒng)具有流動(dòng)性, 易受周邊環(huán)境影響,造成了較高的AOA、AOB多樣性。Bollmann等[28]研究了不同營(yíng)養(yǎng)程度湖泊沉積物的AOA組成, 結(jié)果表明湖泊的營(yíng)養(yǎng)程度是影響氨氧化微生物群落組成的重要因素, 富營(yíng)養(yǎng)湖泊Erie沉積物AOA以Nitrososphaera為主, 寡營(yíng)養(yǎng)湖泊Superior以Nitrosopumilus為主。東湖沉積物AOA主要為Nitrosopumilus (99.8%), 說明東湖郭鄭湖的污染程度相對(duì)較輕, 沉積物TN是影響東湖沉積物AOA群落結(jié)構(gòu)的主要因素。對(duì)于AOB, Nitrosomonas青睞于低氨氮環(huán)境而Nitrosospira更適應(yīng)高氨氮的污染環(huán)境[35], 東湖沉積物AOB主要為Nitrosomonas(55.3%),AOB群落結(jié)構(gòu)主要受TOC和TP的影響, 與Yang等[23]在滇池的研究結(jié)論類似。有研究表明AOA和AOB的群落結(jié)構(gòu)多樣性受總磷的影響[36], 但Liu等[33]在淡水濕地的研究發(fā)現(xiàn)AOA和AOB的群落結(jié)構(gòu)多樣性和總磷沒有顯著相關(guān)性, AOA的多樣性和有效磷含量顯著正相關(guān)。到目前為止, AOA及AOB在含磷酸鹽生態(tài)位中的優(yōu)勢(shì)尚未確定。應(yīng)進(jìn)一步研究磷酸鹽水平與AOA和AOB的存在與活性之間的關(guān)系[37]。

    圖 6 環(huán)境因子對(duì)AOA和AOB群落結(jié)構(gòu)影響CCA分析Fig. 6 CCA analysis of the effects of environmental factors on AOA and AOB community structure

    4 結(jié)論

    在東湖沉積物中, AOA豐度比AOB高, 沉積物中氨氧化反應(yīng)可能由AOA驅(qū)動(dòng)。東湖沉積物樣品中AOA豐度與T呈顯著負(fù)相關(guān), 但AOB豐度變化不明顯。沉積物中Nitrosopumilus為優(yōu)勢(shì)AOA菌, 沉積物TN是影響東湖沉積物AOA群落結(jié)構(gòu)的主要因素; Nitrosomonas為優(yōu)勢(shì)AOB菌, 東湖沉積物AOB群落結(jié)構(gòu)主要受TOC和TP的影響。

    參考文獻(xiàn):

    [1]Fowler D, Coyle M, Skiba U, et al. The global nitrogen cycle in the twenty-first century [J]. Philosophical Transactions of the Royal Society B, 2013, 368(1621): 1—13

    [2]Ministry of Environmental Protection. Annual Report on Environmental Statistics. 2015 [環(huán)保部. 2015年環(huán)境統(tǒng)計(jì)年報(bào)]

    [3]Arauzo M. Harmful effects of un-ionised ammonia on the zooplankton community in a deep waste treatment pond[J]. Water Research, 2003, 37(5): 1048—1054

    [4]Pauer J J, Auer M T. Nitrification in the water column and sediment of a hypereutrophic lake and adjoining river system [J]. Water Research, 2000, 34(4): 1247—1254

    [5]Liu G H, Wu H B. Diversity of ammonia-oxidizing archaea in the surface sediments of the northern continental slope of the South China Sea [J]. Tropical Geography,2016, 36(1): 108—113 [劉國(guó)輝, 吳后波. 南海北部陸坡表層沉積物氨氧化古菌多樣性初探. 熱帶地理, 2016,36(1): 108—113]

    [6]Ke X B, Angel R, Lu Y H, et al. Niche differentiation of ammonia oxidizers and nitrite oxidizers in rice paddy soil[J]. Environmental Microbiology, 2013, 15(8): 2275—2292

    [7]Chen Z B, Xia T Y, Wang D K, et al. Spatial and temporal differences of anammox bacterial distribution in Dianchi Lake sediments [J]. Acta Hydrobiologica Sinica,2016, 40(2): 425—430 [陳澤斌, 夏體淵, 王定康, 等. 滇池沉積物中厭氧氨氧化細(xì)菌分布的時(shí)空差異. 水生生物學(xué)報(bào), 2016, 40(2): 425—430]

    [8]Konneke M, Bernhard A, de la Torre J R, et al. Isolation of an autotrophic ammonia-oxidizing marine archaeon [J].Nature, 2005, 437(22): 543—546

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

    [10]Petersen D G, Blazewicz S J, Firestone M, et al. Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska [J]. Environmental Microbiology, 2012, 14(4): 993—1008

    [11]Han P, Li M, Gu J D. Biases in community structures of ammonia/ammonium-oxidizing microorganisms caused by insufficient DNA extractions from Baijiang soil revealed by comparative analysis of coastal wetland sediment and rice paddy soil [J]. Applied Microbiology and Biotechnology, 2013, 97(19): 8741—8756

    [12]Park S J, Ghai R, Martain-Cuadrado A B, et al. Genomes of two new ammonia-oxidizing archaea enriched from deep marine sediments [J]. PLoS One, 2014, 9(5): e96449

    [13]Zhang J Z, Yang Y Y, et al. Distribution of sediment bacterial and archaealcommunities in plateau freshwater [J].Environmental Biotechnology, 2015, 99(8): 3291—3302

    [14]Limpiyakorn T, Furhacker M, Haberl R, et al. amoA-encoding archaea in wastewater treatment plants: a review[J]. Applied Microbiology and Biotechnology, 2013,97(4): 1425—1439

    [15]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

    [16]Bao S D. Analysis of Soil Aggregation [M]. Beijing:China Agricultural Press. 2010, 123—142 [鮑士旦. 土壤農(nóng)化分析. 北京: 中國(guó)農(nóng)業(yè)出版社. 2010, 123—142]

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

    [18]Zhuang L J, Xia C, Tian Q, et al. Community structure of anaerobic ammonia oxidation bacteria in subsurface soil from riparian zone to land in Baiyangdian Lake using high-throughput sequencing technology [J]. Acta Scientiae Circumstantiae, 2017, 37(1): 261—271 [莊林杰, 夏超,田晴, 等. 高通量測(cè)序技術(shù)研究典型湖泊岸邊陸向深層土壤中厭氧氨氧化細(xì)菌的群落結(jié)構(gòu). 環(huán)境科學(xué)學(xué)報(bào),2017, 37(1): 261—271]

    [19]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

    [20]Zhong J J, Weng J H, Lin X G, et al. Microbial ecology of archaeal ammonia oxidation-a review [J]. Acta Microbiologica Sinca, 2010, 50(4): 431—437 [賈仲君, 翁佳華,林先貴, 等. 氨氧化古菌的生態(tài)學(xué)研究進(jìn)展. 微生物學(xué)報(bào), 2010, 50(4): 431—437]

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

    [22]Hou J, Song C L, Cao X Y, et al. Shifts between ammonia oxidizing bacteria and archaea in relation to nitrification potential across trophic gradients in two large Chinese lakes (Lake Taihuand Lake Chaohu) [J]. Water Research, 2013, 47(12): 2285—2296

    [23]Yang Y Y, Zhang J X, Zhao Q, et al. Sediment Ammonia-Oxidizing Microorganisms in two plateau freshwater lakes at different trophic states [J]. Microbiology of Aquatic Systems, 2016, 71(5): 257—265

    [24]Zhou L L, Zhu G B, Wang S Y, et al. Abundance, biodiversity and contribution to ammonia -oxidizing archaea in littoral sediments of Dongting Lake [J]. Acta Microbiologica Sinca, 2013, 33(6): 1741—1747 [周磊榴, 祝貴兵,王衫允, 等. 洞庭湖岸邊帶沉積物氨氧化古菌的豐度、多樣性及對(duì)氨氧化的貢獻(xiàn). 環(huán)境科學(xué)學(xué)報(bào), 2013, 33(6):1741—1747]

    [25]Wang X Y, Wang C, Bao L L, et al. Abundance and community structure of ammonia-oxidizing microorganisms in reservoir sediment and adjacent soils [J]. Applied Microbiology and Biotechnology, 2014, 98(3): 1883—1892

    [26]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(8): 389—398

    [27]Liu B, Li Y M, Zhang J P, et al. Abundance and diversity of ammonia-oxidizing microorganisms in the sediments of Jinshan Lake [J]. Current Microbiology, 2014, 69(7):751—757

    [28]Wu Y C, Xiang Y, Wang J J, et al. Heterogeneity of archaeal and bacterial ammonia oxidizing communities in Lake Taihu, China [J]. Environmental Microbiology Reports, 2010, 2(33): 569—576

    [29]Bollmann A, Bullerjahn G S, Mckay R M. Abundance and diversity of ammonia-oxidizing archaea and bacteria in sediments of trophic end members of the Laurentian Great Lakes, Erie and Superior [J]. PLoS One, 2014, 9(5):e97068

    [30]Lu S M, Liao M J, Xie C X, et al. Seasonal dynamic of ammonia-oxidizing microorganisms in freshwater aquaculture ponds [J]. Annals of Microbiology, 2015,65(5): 651—657

    [31]Hou J, Cao X Y, Song C L, et al. Predominance of ammonia-oxidizing archaea and nirk-gene-bearing denitrifiers among ammonia-oxidizing and denitrifying populations in sediments of a large urban eutrophic lake(Lake Donghu [J]. Canadian Journal of Microbiology, 2013,59(6): 456—464

    [32]Ying J Y, Li X X, Wang N, et al. Contrasting effects of nitrogen forms and soil pH on ammonia oxidizing microorganisms and their responses to long-term nitrogen fertilization in a typical steppe ecosystem [J]. Soil Biology & Biochemistry, 2017, 107(12): 10—18

    [33]Liu Y, Zhang J X, Zhang X L, et al. Depth-related changes of sediment ammonia-oxidizing microorganisms in a high-altitude freshwater wetland [J]. Environmental Biotechnology, 2014, 98(3): 5697—5707

    [34]Li H, Huang F Y, Su J Q, et al. Distribution and diversity of ammonium-oxidizing archaea and ammonium-oxidizing bacteria in surface sediments of Oujiang River [J].Environmental Science, 2015, 36(12): 4459—4666 [李虎,黃福義, 蘇建強(qiáng), 等. 浙江省甌江氨氧化古菌和氨氧化細(xì)菌分布及多樣性特征. 環(huán)境科學(xué), 2015, 36(12): 4459—4666]

    [35]Wang Y F, Feng Y Y, Gu J D, et al. Seasonal dynamics of ammonia/ammonium-oxidizing prokaryotes in oxic and anoxic wetland sediments of subtropical coastal mangrove [J]. Applied Microbiology and Biotechnology,2013, 97(3): 7919—7934

    [36]Wu Y C, Xiang Y, He J Z, et al. Heterogeneity of archaeal and bacterial ammonia oxidizing communities in Lake Taihu, China [J]. Environmental Microbiology Reports,2010, 2(4): 569—576

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

    猜你喜歡
    東湖沉積物氨氮
    懸浮物對(duì)水質(zhì)氨氮測(cè)定的影響
    化工管理(2022年14期)2022-12-02 11:43:52
    晚更新世以來南黃海陸架沉積物源分析
    改進(jìn)型T-S模糊神經(jīng)網(wǎng)絡(luò)的出水氨氮預(yù)測(cè)
    云南化工(2021年8期)2021-12-21 06:37:36
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    東湖一角
    氧化絮凝技術(shù)處理高鹽高氨氮廢水的實(shí)驗(yàn)研究
    游東湖
    間位芳綸生產(chǎn)廢水氨氮的強(qiáng)化處理及工程實(shí)踐
    討論用ICP-AES測(cè)定土壤和沉積物時(shí)鈦對(duì)鈷的干擾
    精品不卡国产一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 中文精品一卡2卡3卡4更新| 久久99蜜桃精品久久| 国产成人freesex在线| 亚洲精品456在线播放app| 日韩成人av中文字幕在线观看| 免费电影在线观看免费观看| 色哟哟哟哟哟哟| 欧美色欧美亚洲另类二区| 欧美性感艳星| 黑人高潮一二区| www.色视频.com| 国产午夜精品久久久久久一区二区三区| 悠悠久久av| 久久精品91蜜桃| 久久亚洲国产成人精品v| 久久久久久久亚洲中文字幕| 婷婷色av中文字幕| 久99久视频精品免费| 午夜福利在线观看吧| 超碰av人人做人人爽久久| 成人漫画全彩无遮挡| 日韩三级伦理在线观看| 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av天美| 韩国av在线不卡| 99久久无色码亚洲精品果冻| 爱豆传媒免费全集在线观看| 国产三级中文精品| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区三区四区免费观看| 精品国内亚洲2022精品成人| 蜜桃久久精品国产亚洲av| 深夜精品福利| 日本熟妇午夜| 色尼玛亚洲综合影院| 久99久视频精品免费| 国产成人a∨麻豆精品| 亚洲欧洲日产国产| 少妇丰满av| 午夜福利高清视频| 白带黄色成豆腐渣| 国产欧美日韩精品一区二区| 欧美成人免费av一区二区三区| 一本一本综合久久| 日日干狠狠操夜夜爽| 亚洲熟妇中文字幕五十中出| 91av网一区二区| 麻豆久久精品国产亚洲av| 国产成人freesex在线| 99精品在免费线老司机午夜| 老司机福利观看| 1024手机看黄色片| 69人妻影院| 91在线精品国自产拍蜜月| 一个人免费在线观看电影| 精品午夜福利在线看| 日本色播在线视频| 亚洲精品乱码久久久久久按摩| 亚洲一区二区三区色噜噜| 十八禁国产超污无遮挡网站| 免费看av在线观看网站| 亚洲精品影视一区二区三区av| 亚洲七黄色美女视频| 日日干狠狠操夜夜爽| 男女视频在线观看网站免费| 欧美一区二区亚洲| a级毛片免费高清观看在线播放| 中文字幕熟女人妻在线| 国产黄色小视频在线观看| 听说在线观看完整版免费高清| 26uuu在线亚洲综合色| 日本成人三级电影网站| 中国美白少妇内射xxxbb| kizo精华| 精品少妇黑人巨大在线播放 | 国产成人福利小说| 久久精品国产自在天天线| 在线观看一区二区三区| 国产高清激情床上av| 在线观看一区二区三区| 男女下面进入的视频免费午夜| 一区福利在线观看| 成人毛片a级毛片在线播放| 三级男女做爰猛烈吃奶摸视频| 中文资源天堂在线| 国产伦理片在线播放av一区 | 国产白丝娇喘喷水9色精品| 婷婷亚洲欧美| 成人午夜精彩视频在线观看| 亚洲国产色片| 在线免费观看不下载黄p国产| 国产精品一区www在线观看| 国产麻豆成人av免费视频| 久久久久网色| 青青草视频在线视频观看| 床上黄色一级片| av又黄又爽大尺度在线免费看 | 国产探花极品一区二区| 尾随美女入室| 免费av观看视频| 国产 一区 欧美 日韩| 老司机福利观看| 久久精品国产鲁丝片午夜精品| 国产一级毛片在线| 国产片特级美女逼逼视频| 亚洲婷婷狠狠爱综合网| 我要搜黄色片| 九九热线精品视视频播放| 国产精品久久久久久精品电影| 最新中文字幕久久久久| 国内揄拍国产精品人妻在线| 国产亚洲精品久久久久久毛片| av在线亚洲专区| 亚洲在线自拍视频| av专区在线播放| 久久亚洲国产成人精品v| 欧美最黄视频在线播放免费| 观看免费一级毛片| 丝袜喷水一区| 小蜜桃在线观看免费完整版高清| 内射极品少妇av片p| 99热这里只有是精品50| 久久精品夜色国产| 婷婷色综合大香蕉| 村上凉子中文字幕在线| 99久久精品热视频| av又黄又爽大尺度在线免费看 | 国产久久久一区二区三区| 晚上一个人看的免费电影| 又爽又黄a免费视频| 天堂影院成人在线观看| 国产又黄又爽又无遮挡在线| 赤兔流量卡办理| 国产精品国产高清国产av| 国产精品一区二区在线观看99 | 边亲边吃奶的免费视频| 精品不卡国产一区二区三区| 欧美区成人在线视频| 亚洲精品乱码久久久久久按摩| 精华霜和精华液先用哪个| 国产精品久久久久久精品电影| 亚洲精品亚洲一区二区| 欧美日韩综合久久久久久| 九草在线视频观看| 精品久久久久久久人妻蜜臀av| 国内精品久久久久精免费| 赤兔流量卡办理| 又爽又黄a免费视频| 国产69精品久久久久777片| 免费看av在线观看网站| 国产一级毛片在线| 亚洲美女搞黄在线观看| 99久久成人亚洲精品观看| 精品久久久久久久久久免费视频| 亚洲在线自拍视频| 啦啦啦啦在线视频资源| av在线播放精品| 久久精品久久久久久噜噜老黄 | 天天一区二区日本电影三级| 国产伦一二天堂av在线观看| 免费av观看视频| 免费av观看视频| 亚洲国产日韩欧美精品在线观看| 少妇的逼水好多| a级一级毛片免费在线观看| 成人欧美大片| 国内精品久久久久精免费| 青春草视频在线免费观看| 联通29元200g的流量卡| 午夜精品一区二区三区免费看| 99热这里只有是精品50| 日韩视频在线欧美| 可以在线观看毛片的网站| 在线国产一区二区在线| 国产亚洲精品av在线| 三级经典国产精品| 欧美精品一区二区大全| 国产精品乱码一区二三区的特点| 熟妇人妻久久中文字幕3abv| 搞女人的毛片| 国产白丝娇喘喷水9色精品| 亚洲精品亚洲一区二区| 国内揄拍国产精品人妻在线| 丰满人妻一区二区三区视频av| 国国产精品蜜臀av免费| 我要搜黄色片| 国产亚洲精品久久久com| 热99re8久久精品国产| 午夜爱爱视频在线播放| 亚洲人与动物交配视频| a级毛色黄片| 99热这里只有精品一区| 国内精品久久久久精免费| 亚洲国产精品合色在线| 国产精品一区二区性色av| 麻豆久久精品国产亚洲av| 午夜福利高清视频| 18禁黄网站禁片免费观看直播| 国产私拍福利视频在线观看| 久久人人爽人人爽人人片va| 精品免费久久久久久久清纯| 成人性生交大片免费视频hd| 特大巨黑吊av在线直播| 久久精品国产鲁丝片午夜精品| 看黄色毛片网站| 美女脱内裤让男人舔精品视频 | 久久久久久九九精品二区国产| 日韩精品青青久久久久久| 日韩欧美一区二区三区在线观看| 嫩草影院精品99| 卡戴珊不雅视频在线播放| 国产淫片久久久久久久久| 毛片女人毛片| 久久精品夜色国产| 老女人水多毛片| 国产伦一二天堂av在线观看| 97超碰精品成人国产| 日韩一区二区视频免费看| 国产探花极品一区二区| 日本成人三级电影网站| 三级经典国产精品| 欧美日韩综合久久久久久| 久久欧美精品欧美久久欧美| 精品国产三级普通话版| 变态另类丝袜制服| 欧美zozozo另类| 国产精华一区二区三区| 午夜精品国产一区二区电影 | 精品久久久久久久久av| 天堂√8在线中文| 深爱激情五月婷婷| 男插女下体视频免费在线播放| 久久精品国产亚洲网站| 精品99又大又爽又粗少妇毛片| 免费观看精品视频网站| 大香蕉久久网| 美女脱内裤让男人舔精品视频 | 国产午夜福利久久久久久| 亚洲成人久久性| 精华霜和精华液先用哪个| 亚洲精品国产成人久久av| 亚洲av第一区精品v没综合| 欧美xxxx黑人xx丫x性爽| 18禁裸乳无遮挡免费网站照片| 韩国av在线不卡| 亚洲精品乱码久久久v下载方式| 国产激情偷乱视频一区二区| 欧美最新免费一区二区三区| 国产女主播在线喷水免费视频网站 | 最近最新中文字幕大全电影3| 成人永久免费在线观看视频| 两个人的视频大全免费| 日日啪夜夜撸| 此物有八面人人有两片| 日韩制服骚丝袜av| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 嘟嘟电影网在线观看| 亚洲av中文字字幕乱码综合| 欧美三级亚洲精品| 国产高清激情床上av| 深爱激情五月婷婷| 久久午夜亚洲精品久久| 中出人妻视频一区二区| 欧美高清成人免费视频www| 国产精品一区二区在线观看99 | 久久这里只有精品中国| av在线天堂中文字幕| 99在线视频只有这里精品首页| 夜夜夜夜夜久久久久| 午夜福利成人在线免费观看| 亚洲在线自拍视频| 婷婷精品国产亚洲av| 男女啪啪激烈高潮av片| 91av网一区二区| 91麻豆精品激情在线观看国产| 中文字幕熟女人妻在线| 狠狠狠狠99中文字幕| 少妇人妻精品综合一区二区 | 午夜免费男女啪啪视频观看| 国产亚洲精品久久久久久毛片| 欧美日韩精品成人综合77777| 色综合站精品国产| 久久久久九九精品影院| 亚洲一区高清亚洲精品| 赤兔流量卡办理| av在线蜜桃| 国产 一区精品| 在线免费观看的www视频| 国产精品爽爽va在线观看网站| 18禁裸乳无遮挡免费网站照片| 亚洲一区高清亚洲精品| 精品国内亚洲2022精品成人| 人人妻人人澡人人爽人人夜夜 | 国产成人精品婷婷| 免费人成视频x8x8入口观看| 在线免费观看不下载黄p国产| 国产又黄又爽又无遮挡在线| 免费看光身美女| 欧美+亚洲+日韩+国产| 久久久久久久久久黄片| 69av精品久久久久久| 亚洲一区高清亚洲精品| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 亚洲激情五月婷婷啪啪| 在线免费观看的www视频| 久久久久久大精品| 高清毛片免费观看视频网站| 婷婷色av中文字幕| 国产伦精品一区二区三区四那| 午夜福利高清视频| 淫秽高清视频在线观看| 麻豆av噜噜一区二区三区| 亚洲精华国产精华液的使用体验 | 久久久久久大精品| 一本一本综合久久| 国产伦精品一区二区三区四那| 久久精品人妻少妇| 亚洲无线在线观看| 国产精品一区www在线观看| 亚洲人成网站在线播| 国产伦精品一区二区三区四那| 欧美激情久久久久久爽电影| 在现免费观看毛片| 免费看日本二区| 搞女人的毛片| 黄片wwwwww| 色播亚洲综合网| 色视频www国产| 精品国内亚洲2022精品成人| 插逼视频在线观看| 久久这里只有精品中国| 亚洲精品国产av成人精品| 亚洲天堂国产精品一区在线| 黄色配什么色好看| 99热网站在线观看| 久久久色成人| 欧美一级a爱片免费观看看| 国产精品一区二区在线观看99 | 婷婷六月久久综合丁香| 一边亲一边摸免费视频| 久久午夜福利片| 天堂影院成人在线观看| 久久精品91蜜桃| 97在线视频观看| 99热6这里只有精品| 亚洲三级黄色毛片| 日本黄色视频三级网站网址| 91aial.com中文字幕在线观看| 狂野欧美激情性xxxx在线观看| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 边亲边吃奶的免费视频| 亚洲激情五月婷婷啪啪| 国产不卡一卡二| or卡值多少钱| 美女 人体艺术 gogo| 在线播放无遮挡| 国产一区二区三区av在线 | 国产亚洲av片在线观看秒播厂 | 国产精品久久久久久精品电影| 亚洲精品国产成人久久av| 最后的刺客免费高清国语| 丝袜喷水一区| 久久热精品热| 亚洲精品粉嫩美女一区| 日韩av在线大香蕉| 欧美极品一区二区三区四区| www.av在线官网国产| 91av网一区二区| 狠狠狠狠99中文字幕| 特大巨黑吊av在线直播| 亚洲欧美日韩东京热| www日本黄色视频网| 国产黄片视频在线免费观看| 国产v大片淫在线免费观看| 看免费成人av毛片| 久久精品久久久久久噜噜老黄 | 亚洲精品粉嫩美女一区| 国产美女午夜福利| 久久这里有精品视频免费| 午夜福利高清视频| 亚洲精品国产成人久久av| 久久人人爽人人片av| 亚洲无线观看免费| 亚洲成人精品中文字幕电影| 超碰av人人做人人爽久久| 如何舔出高潮| 久久久久国产网址| 国产真实伦视频高清在线观看| 最近视频中文字幕2019在线8| 成年免费大片在线观看| 一级毛片aaaaaa免费看小| 神马国产精品三级电影在线观看| 亚洲欧洲国产日韩| 中文字幕熟女人妻在线| 插阴视频在线观看视频| 国产在线精品亚洲第一网站| 亚洲精品自拍成人| 一区二区三区四区激情视频 | av天堂在线播放| 国产精品日韩av在线免费观看| 国产日本99.免费观看| 成人高潮视频无遮挡免费网站| 国产精品国产三级国产av玫瑰| 熟女人妻精品中文字幕| 国产精品乱码一区二三区的特点| 少妇丰满av| 久久久精品94久久精品| 中文字幕精品亚洲无线码一区| 日韩欧美在线乱码| 22中文网久久字幕| 久久韩国三级中文字幕| 久久这里有精品视频免费| 3wmmmm亚洲av在线观看| 一本一本综合久久| 久久久精品大字幕| 我的老师免费观看完整版| 99九九线精品视频在线观看视频| 蜜桃久久精品国产亚洲av| 成年女人看的毛片在线观看| 中国国产av一级| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 黑人高潮一二区| 美女高潮的动态| 亚洲人成网站在线播放欧美日韩| 亚洲内射少妇av| 久久久久网色| 久久亚洲国产成人精品v| 久久久成人免费电影| 午夜激情福利司机影院| 少妇熟女aⅴ在线视频| 日韩欧美三级三区| 色视频www国产| 午夜福利在线观看免费完整高清在 | 国产一区二区在线av高清观看| 99热这里只有是精品在线观看| 亚洲在久久综合| 精品一区二区免费观看| 久久九九热精品免费| 国产伦一二天堂av在线观看| 97超视频在线观看视频| 国产精品一二三区在线看| 91狼人影院| 亚洲精品日韩av片在线观看| 成人午夜精彩视频在线观看| 最近最新中文字幕大全电影3| 国产女主播在线喷水免费视频网站 | 夜夜爽天天搞| 午夜免费男女啪啪视频观看| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 日本黄大片高清| 91久久精品国产一区二区成人| 尤物成人国产欧美一区二区三区| 亚洲第一区二区三区不卡| 激情 狠狠 欧美| 亚洲第一电影网av| 人人妻人人澡欧美一区二区| 国产精华一区二区三区| 亚洲最大成人中文| 久久久精品94久久精品| av在线老鸭窝| 日韩中字成人| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 亚洲七黄色美女视频| 免费av观看视频| 变态另类丝袜制服| 人妻久久中文字幕网| 亚洲成人久久性| 成人无遮挡网站| 成人永久免费在线观看视频| 一边亲一边摸免费视频| 尤物成人国产欧美一区二区三区| 欧美日本亚洲视频在线播放| 春色校园在线视频观看| 搞女人的毛片| 国产精品一区二区三区四区久久| 能在线免费观看的黄片| 91精品一卡2卡3卡4卡| 五月伊人婷婷丁香| 日本欧美国产在线视频| 不卡视频在线观看欧美| 美女黄网站色视频| 美女高潮的动态| av天堂中文字幕网| 日韩制服骚丝袜av| 免费观看人在逋| 日韩欧美国产在线观看| 色综合色国产| 久久久久久伊人网av| 夜夜夜夜夜久久久久| av天堂中文字幕网| 又粗又硬又长又爽又黄的视频 | 精品一区二区三区视频在线| 麻豆成人av视频| 日韩三级伦理在线观看| 夜夜夜夜夜久久久久| 99视频精品全部免费 在线| 成人特级av手机在线观看| 免费黄网站久久成人精品| 欧美精品一区二区大全| 久久韩国三级中文字幕| 亚洲国产日韩欧美精品在线观看| 国产精品免费一区二区三区在线| 性欧美人与动物交配| 女同久久另类99精品国产91| 我要看日韩黄色一级片| 级片在线观看| 只有这里有精品99| 精品一区二区三区视频在线| 欧美最新免费一区二区三区| 搡老妇女老女人老熟妇| 看片在线看免费视频| 又粗又硬又长又爽又黄的视频 | 国产成人午夜福利电影在线观看| 国产精品一及| 亚洲久久久久久中文字幕| 最好的美女福利视频网| 亚洲国产精品成人久久小说 | 色播亚洲综合网| 精品少妇黑人巨大在线播放 | 欧美性猛交黑人性爽| 亚洲在线观看片| av免费在线看不卡| 校园春色视频在线观看| 国产一区二区三区av在线 | 91午夜精品亚洲一区二区三区| 午夜激情欧美在线| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 高清午夜精品一区二区三区 | 色5月婷婷丁香| 国产av不卡久久| 国产 一区 欧美 日韩| 最新中文字幕久久久久| 美女高潮的动态| 日韩av不卡免费在线播放| 男人和女人高潮做爰伦理| 国产高清不卡午夜福利| 99国产精品一区二区蜜桃av| 国产午夜精品论理片| 欧美性猛交黑人性爽| 免费观看精品视频网站| 我的女老师完整版在线观看| 国产91av在线免费观看| 成人三级黄色视频| a级毛色黄片| 亚洲精品亚洲一区二区| 精品一区二区三区视频在线| 全区人妻精品视频| 中国美女看黄片| 老熟妇乱子伦视频在线观看| 综合色丁香网| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 日韩大尺度精品在线看网址| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 国产不卡一卡二| 99久国产av精品| 91aial.com中文字幕在线观看| 深夜a级毛片| 秋霞在线观看毛片| 欧美+亚洲+日韩+国产| 国产精品一区www在线观看| 成人毛片60女人毛片免费| 日本欧美国产在线视频| 日本成人三级电影网站| 欧美成人a在线观看| 婷婷精品国产亚洲av| 欧美xxxx黑人xx丫x性爽| 在线观看66精品国产| 美女内射精品一级片tv| 在线观看66精品国产| 欧美性猛交黑人性爽| 成人午夜精彩视频在线观看| 麻豆国产97在线/欧美| av专区在线播放| 麻豆国产97在线/欧美| 免费一级毛片在线播放高清视频| 97在线视频观看| 亚洲国产日韩欧美精品在线观看| 青青草视频在线视频观看| 成人高潮视频无遮挡免费网站| 少妇猛男粗大的猛烈进出视频 | 欧美+亚洲+日韩+国产| 欧美色欧美亚洲另类二区| 中文欧美无线码| 女的被弄到高潮叫床怎么办| 少妇丰满av| 国产乱人偷精品视频| 97超视频在线观看视频| 精品久久久久久久久久免费视频| 国产亚洲91精品色在线| 亚洲三级黄色毛片| 亚洲精品粉嫩美女一区| 国产精品野战在线观看| 欧美一区二区国产精品久久精品| a级一级毛片免费在线观看| 日韩强制内射视频| 欧美+日韩+精品| 女人十人毛片免费观看3o分钟| 亚洲色图av天堂| 最后的刺客免费高清国语| 久久精品91蜜桃| 亚洲欧美中文字幕日韩二区| 在线观看av片永久免费下载| 国产日本99.免费观看| 久久国产乱子免费精品|