• <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ì)鈷的干擾
    中文字幕人成人乱码亚洲影| 亚洲最大成人中文| 国产伦人伦偷精品视频| 久久久国产精品麻豆| 国产精品一区二区三区四区久久| 国产欧美日韩精品亚洲av| 亚洲激情在线av| 国产精品久久久久久精品电影| 欧美黑人精品巨大| 岛国视频午夜一区免费看| 亚洲国产欧洲综合997久久,| 国产精品,欧美在线| 变态另类丝袜制服| 精品熟女少妇八av免费久了| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 久久99热这里只有精品18| av视频在线观看入口| 最新在线观看一区二区三区| 一本大道久久a久久精品| 男人舔奶头视频| 久久精品影院6| 国产亚洲av嫩草精品影院| svipshipincom国产片| 欧美激情久久久久久爽电影| 国产不卡一卡二| 免费观看人在逋| 美女 人体艺术 gogo| av视频在线观看入口| 成人av在线播放网站| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看影片大全网站| 久久久久久久午夜电影| 99久久精品国产亚洲精品| 国产激情欧美一区二区| 97碰自拍视频| 久久久久久大精品| 亚洲免费av在线视频| 欧美高清成人免费视频www| 男人舔女人的私密视频| 亚洲中文av在线| 在线观看日韩欧美| 欧美+亚洲+日韩+国产| 亚洲熟妇中文字幕五十中出| 国产欧美日韩一区二区精品| 波多野结衣巨乳人妻| 亚洲专区中文字幕在线| 极品教师在线免费播放| 女人被狂操c到高潮| 成人永久免费在线观看视频| 国产精品久久久av美女十八| 亚洲黑人精品在线| 欧美大码av| 成人特级黄色片久久久久久久| 在线永久观看黄色视频| 久久午夜亚洲精品久久| 波多野结衣巨乳人妻| 日韩精品青青久久久久久| 久久久久久国产a免费观看| 午夜精品久久久久久毛片777| 99热这里只有是精品50| 日本黄大片高清| 男女床上黄色一级片免费看| 99久久国产精品久久久| 熟妇人妻久久中文字幕3abv| 精华霜和精华液先用哪个| 黄频高清免费视频| 亚洲国产精品久久男人天堂| 欧美三级亚洲精品| 国产精品九九99| 久99久视频精品免费| 可以在线观看毛片的网站| 一区二区三区高清视频在线| 国产一区二区在线av高清观看| 久久久国产成人免费| 国内精品久久久久精免费| 亚洲中文字幕一区二区三区有码在线看 | 不卡av一区二区三区| 无遮挡黄片免费观看| 国产成年人精品一区二区| 亚洲熟女毛片儿| 亚洲精品中文字幕一二三四区| 淫秽高清视频在线观看| 99热只有精品国产| 天堂动漫精品| 久久久久久大精品| 男女之事视频高清在线观看| 欧美黑人巨大hd| 国产熟女午夜一区二区三区| 欧美乱妇无乱码| 国产69精品久久久久777片 | av片东京热男人的天堂| 日本三级黄在线观看| 91在线观看av| 国内精品一区二区在线观看| 又粗又爽又猛毛片免费看| 最近视频中文字幕2019在线8| 国产真实乱freesex| 国产1区2区3区精品| 国产精品久久久久久亚洲av鲁大| 亚洲av熟女| 午夜成年电影在线免费观看| 欧美日韩黄片免| 国产成年人精品一区二区| 一级作爱视频免费观看| netflix在线观看网站| 国产三级在线视频| 丝袜美腿诱惑在线| bbb黄色大片| 欧美一级毛片孕妇| 日本黄大片高清| 亚洲国产精品999在线| 午夜视频精品福利| 在线视频色国产色| 免费电影在线观看免费观看| 亚洲av熟女| 中文字幕精品亚洲无线码一区| 两人在一起打扑克的视频| 好男人在线观看高清免费视频| 欧美精品啪啪一区二区三区| 男女午夜视频在线观看| 视频区欧美日本亚洲| 九九热线精品视视频播放| 久久中文字幕一级| 麻豆av在线久日| 午夜老司机福利片| 午夜福利成人在线免费观看| 国产激情欧美一区二区| 老熟妇乱子伦视频在线观看| 精品国产超薄肉色丝袜足j| 一区福利在线观看| 久久久精品欧美日韩精品| 一本一本综合久久| 精品久久蜜臀av无| 日本一本二区三区精品| 日韩国内少妇激情av| ponron亚洲| 十八禁网站免费在线| cao死你这个sao货| 色av中文字幕| 99久久久亚洲精品蜜臀av| a级毛片a级免费在线| 老司机午夜福利在线观看视频| 国产精品98久久久久久宅男小说| 久久99热这里只有精品18| 精品乱码久久久久久99久播| 制服人妻中文乱码| 亚洲av成人一区二区三| 成人av在线播放网站| 国产av麻豆久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 国产成人系列免费观看| 国产精品久久久人人做人人爽| 久久久久久国产a免费观看| 人成视频在线观看免费观看| 精品无人区乱码1区二区| 一级毛片女人18水好多| 白带黄色成豆腐渣| 黄色女人牲交| 亚洲人成77777在线视频| tocl精华| 一本精品99久久精品77| 成人一区二区视频在线观看| 国产99白浆流出| 久久精品人妻少妇| 日韩欧美三级三区| 听说在线观看完整版免费高清| 午夜精品在线福利| 欧美国产日韩亚洲一区| 婷婷精品国产亚洲av| av国产免费在线观看| 国产精品九九99| 琪琪午夜伦伦电影理论片6080| 中文字幕精品亚洲无线码一区| 色av中文字幕| 亚洲av片天天在线观看| 久久久久久九九精品二区国产 | 亚洲一码二码三码区别大吗| 亚洲欧美激情综合另类| 校园春色视频在线观看| 国产精品98久久久久久宅男小说| 亚洲国产精品sss在线观看| 欧美黄色片欧美黄色片| 欧美日韩亚洲综合一区二区三区_| 久久久精品国产亚洲av高清涩受| 国产1区2区3区精品| 女同久久另类99精品国产91| 91字幕亚洲| 国产亚洲精品久久久久5区| 制服诱惑二区| 亚洲第一欧美日韩一区二区三区| 叶爱在线成人免费视频播放| 亚洲专区字幕在线| 色哟哟哟哟哟哟| 日韩欧美一区二区三区在线观看| 极品教师在线免费播放| 中文字幕精品亚洲无线码一区| 91在线观看av| 不卡一级毛片| 精品国产美女av久久久久小说| 男女之事视频高清在线观看| 可以在线观看的亚洲视频| 91九色精品人成在线观看| 美女大奶头视频| 一本久久中文字幕| 日日摸夜夜添夜夜添小说| 亚洲精品在线观看二区| 十八禁网站免费在线| 精品欧美国产一区二区三| 精品熟女少妇八av免费久了| 日本在线视频免费播放| 中文字幕久久专区| 又黄又爽又免费观看的视频| 三级国产精品欧美在线观看 | 啦啦啦免费观看视频1| 久久精品aⅴ一区二区三区四区| 亚洲精品久久国产高清桃花| 老司机福利观看| 欧美精品啪啪一区二区三区| 黄色 视频免费看| 亚洲 欧美 日韩 在线 免费| 精品午夜福利视频在线观看一区| 亚洲男人天堂网一区| 国产黄a三级三级三级人| 久久精品91无色码中文字幕| 99精品久久久久人妻精品| 亚洲中文日韩欧美视频| 国产高清视频在线播放一区| 午夜免费观看网址| 中国美女看黄片| 中出人妻视频一区二区| 免费在线观看视频国产中文字幕亚洲| 国产精品影院久久| av有码第一页| 1024手机看黄色片| 亚洲欧美日韩高清在线视频| 精品不卡国产一区二区三区| 国产精品久久久久久精品电影| 亚洲精华国产精华精| 精品久久久久久成人av| 精品福利观看| 免费av毛片视频| 给我免费播放毛片高清在线观看| 国产免费男女视频| 午夜福利免费观看在线| 啦啦啦观看免费观看视频高清| 男女下面进入的视频免费午夜| 欧美成人一区二区免费高清观看 | 变态另类成人亚洲欧美熟女| 国产精品免费视频内射| 精华霜和精华液先用哪个| 久久中文字幕人妻熟女| 久久久久性生活片| 国产熟女xx| 亚洲专区中文字幕在线| 天天添夜夜摸| 精品一区二区三区四区五区乱码| 国产激情偷乱视频一区二区| 12—13女人毛片做爰片一| 好看av亚洲va欧美ⅴa在| 国产三级在线视频| 两个人视频免费观看高清| 国产精品久久电影中文字幕| 免费在线观看视频国产中文字幕亚洲| 国产91精品成人一区二区三区| 狂野欧美激情性xxxx| 黄色 视频免费看| 最近最新免费中文字幕在线| 岛国在线观看网站| a级毛片a级免费在线| 亚洲第一欧美日韩一区二区三区| 很黄的视频免费| 十八禁人妻一区二区| 久久香蕉激情| 可以免费在线观看a视频的电影网站| 国产成人精品久久二区二区免费| 日本黄色视频三级网站网址| 国内揄拍国产精品人妻在线| 久久亚洲真实| 国产片内射在线| 亚洲精品在线美女| 亚洲国产欧洲综合997久久,| 久99久视频精品免费| netflix在线观看网站| 97碰自拍视频| 日韩欧美免费精品| 久久草成人影院| 又紧又爽又黄一区二区| 久久99热这里只有精品18| 18美女黄网站色大片免费观看| 国产又黄又爽又无遮挡在线| 午夜视频精品福利| 免费观看人在逋| 男女那种视频在线观看| 99久久久亚洲精品蜜臀av| 国模一区二区三区四区视频 | 999精品在线视频| 1024视频免费在线观看| 亚洲av五月六月丁香网| 久久久久久久精品吃奶| 日韩精品中文字幕看吧| 男女之事视频高清在线观看| 搞女人的毛片| 日本免费a在线| 看片在线看免费视频| 日韩欧美三级三区| 两个人看的免费小视频| 国产精品,欧美在线| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| www国产在线视频色| 亚洲精品在线观看二区| 在线看三级毛片| 老熟妇仑乱视频hdxx| 在线观看免费日韩欧美大片| 全区人妻精品视频| 日本一二三区视频观看| 精品国产亚洲在线| 成年人黄色毛片网站| 在线观看免费视频日本深夜| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| 听说在线观看完整版免费高清| 国产1区2区3区精品| 长腿黑丝高跟| 可以在线观看毛片的网站| 99精品欧美一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| www日本黄色视频网| 首页视频小说图片口味搜索| 女人高潮潮喷娇喘18禁视频| 国产97色在线日韩免费| 天天躁狠狠躁夜夜躁狠狠躁| 可以免费在线观看a视频的电影网站| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美98| 99国产极品粉嫩在线观看| 麻豆成人午夜福利视频| 成人av一区二区三区在线看| 麻豆成人午夜福利视频| 18禁观看日本| 一进一出抽搐动态| 国产三级中文精品| 非洲黑人性xxxx精品又粗又长| 日韩欧美 国产精品| 亚洲欧美日韩高清专用| ponron亚洲| 国产真实乱freesex| 黄色 视频免费看| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 欧美色视频一区免费| 亚洲一区二区三区色噜噜| 成人手机av| 精品国产亚洲在线| 欧美人与性动交α欧美精品济南到| 欧美大码av| 黄色女人牲交| 欧美久久黑人一区二区| 老汉色av国产亚洲站长工具| 天堂av国产一区二区熟女人妻 | 国产精品 欧美亚洲| 99re在线观看精品视频| 国产精品 欧美亚洲| 99久久精品热视频| 中文在线观看免费www的网站 | 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 国产激情欧美一区二区| 亚洲av电影在线进入| 久久久久久久久中文| 淫秽高清视频在线观看| 成人永久免费在线观看视频| 国产高清有码在线观看视频 | 亚洲狠狠婷婷综合久久图片| 特级一级黄色大片| 亚洲黑人精品在线| 欧美成人午夜精品| 香蕉国产在线看| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频| 两个人的视频大全免费| 国产成人av教育| 日韩精品青青久久久久久| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 老司机午夜福利在线观看视频| www日本黄色视频网| 亚洲av电影不卡..在线观看| 国产免费av片在线观看野外av| 高清在线国产一区| 亚洲18禁久久av| 啦啦啦韩国在线观看视频| 久久香蕉国产精品| 搡老熟女国产l中国老女人| 成在线人永久免费视频| 天天一区二区日本电影三级| 久久天躁狠狠躁夜夜2o2o| av福利片在线观看| 久久午夜亚洲精品久久| 97碰自拍视频| 精品久久蜜臀av无| 国产精品久久久人人做人人爽| 欧美大码av| 亚洲国产欧洲综合997久久,| 精品高清国产在线一区| 国产亚洲精品久久久久5区| 757午夜福利合集在线观看| 一个人观看的视频www高清免费观看 | 亚洲av成人一区二区三| 51午夜福利影视在线观看| 中文字幕高清在线视频| 黄片小视频在线播放| 亚洲成人久久爱视频| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 亚洲精品久久国产高清桃花| 久久国产精品人妻蜜桃| 中文字幕人成人乱码亚洲影| 美女 人体艺术 gogo| 亚洲中文字幕日韩| or卡值多少钱| 亚洲人成网站高清观看| 18禁黄网站禁片免费观看直播| 欧美乱妇无乱码| 国产精品一区二区三区四区免费观看 | 欧美日韩福利视频一区二区| 欧美乱色亚洲激情| 免费看美女性在线毛片视频| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 欧美一级毛片孕妇| 国产欧美日韩一区二区三| 久久久精品大字幕| 国产精品自产拍在线观看55亚洲| 欧美激情久久久久久爽电影| 黄片大片在线免费观看| 午夜福利高清视频| 久久久国产成人精品二区| 久久精品夜夜夜夜夜久久蜜豆 | 国产在线观看jvid| 日韩欧美国产在线观看| av有码第一页| 日本熟妇午夜| 男人的好看免费观看在线视频 | 日本 欧美在线| 欧美日韩亚洲综合一区二区三区_| 亚洲人成网站高清观看| 97人妻精品一区二区三区麻豆| 成人av在线播放网站| 欧美精品亚洲一区二区| 日本免费a在线| 中文字幕av在线有码专区| 亚洲激情在线av| 色综合欧美亚洲国产小说| 91字幕亚洲| 久久性视频一级片| 国产精品国产高清国产av| 黄色视频,在线免费观看| 在线播放国产精品三级| 色在线成人网| 两性夫妻黄色片| 国产精品98久久久久久宅男小说| 国产精品自产拍在线观看55亚洲| 精品国产乱子伦一区二区三区| 国产不卡一卡二| 成人国语在线视频| 国内少妇人妻偷人精品xxx网站 | 女警被强在线播放| 久久人人精品亚洲av| 亚洲avbb在线观看| 一二三四社区在线视频社区8| 真人做人爱边吃奶动态| 久久久久久久精品吃奶| 老汉色∧v一级毛片| 亚洲欧美激情综合另类| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清在线视频| 日韩欧美免费精品| 熟女电影av网| 在线看三级毛片| 欧美黑人精品巨大| 日日爽夜夜爽网站| 日日干狠狠操夜夜爽| 一区二区三区激情视频| 久久久精品国产亚洲av高清涩受| 国内揄拍国产精品人妻在线| 脱女人内裤的视频| 在线观看美女被高潮喷水网站 | 精品免费久久久久久久清纯| 日韩中文字幕欧美一区二区| 精品一区二区三区四区五区乱码| 久久精品aⅴ一区二区三区四区| 脱女人内裤的视频| 狠狠狠狠99中文字幕| 在线观看www视频免费| 久久久久性生活片| 他把我摸到了高潮在线观看| 亚洲中文日韩欧美视频| 国产精品 国内视频| 午夜视频精品福利| 午夜福利成人在线免费观看| 亚洲精品美女久久av网站| 成人永久免费在线观看视频| 男插女下体视频免费在线播放| 人人妻,人人澡人人爽秒播| 成熟少妇高潮喷水视频| 变态另类丝袜制服| 亚洲午夜精品一区,二区,三区| tocl精华| 巨乳人妻的诱惑在线观看| 淫妇啪啪啪对白视频| 极品教师在线免费播放| 中文亚洲av片在线观看爽| 老司机午夜十八禁免费视频| 国产精品久久久久久精品电影| 一区二区三区国产精品乱码| 日本 欧美在线| 男人舔女人的私密视频| 午夜成年电影在线免费观看| 国产黄a三级三级三级人| 国产免费男女视频| 亚洲精品久久国产高清桃花| 午夜福利高清视频| 狂野欧美激情性xxxx| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 中文字幕熟女人妻在线| 亚洲午夜精品一区,二区,三区| 午夜福利视频1000在线观看| 亚洲真实伦在线观看| 成人亚洲精品av一区二区| 亚洲成人国产一区在线观看| 成人精品一区二区免费| а√天堂www在线а√下载| 成人特级黄色片久久久久久久| bbb黄色大片| 亚洲成人中文字幕在线播放| 久久国产精品人妻蜜桃| 亚洲av电影在线进入| 国产精品一区二区三区四区免费观看 | 国产黄a三级三级三级人| 91老司机精品| 嫩草影院精品99| 亚洲欧美精品综合一区二区三区| 99精品久久久久人妻精品| 天天躁狠狠躁夜夜躁狠狠躁| 制服人妻中文乱码| a级毛片a级免费在线| 国产精品亚洲美女久久久| 中亚洲国语对白在线视频| 成人av一区二区三区在线看| 成年女人毛片免费观看观看9| 他把我摸到了高潮在线观看| 99国产精品99久久久久| 亚洲aⅴ乱码一区二区在线播放 | 18禁裸乳无遮挡免费网站照片| 99久久久亚洲精品蜜臀av| 又爽又黄无遮挡网站| 亚洲成人国产一区在线观看| 好看av亚洲va欧美ⅴa在| 18美女黄网站色大片免费观看| 午夜福利高清视频| 精品日产1卡2卡| 中文字幕人妻丝袜一区二区| 在线观看舔阴道视频| 夜夜爽天天搞| 国产爱豆传媒在线观看 | 欧美黑人精品巨大| 黑人欧美特级aaaaaa片| 欧美性猛交黑人性爽| 又黄又爽又免费观看的视频| 1024手机看黄色片| 丁香欧美五月| 黄色a级毛片大全视频| 日韩欧美免费精品| 一区二区三区激情视频| 亚洲成a人片在线一区二区| 国产精品久久久久久精品电影| 亚洲av成人一区二区三| 国产精品久久视频播放| 欧美又色又爽又黄视频| 国产精品久久电影中文字幕| 国产一区二区在线av高清观看| 又黄又爽又免费观看的视频| 99在线视频只有这里精品首页| 国产在线观看jvid| 午夜精品一区二区三区免费看| 99国产综合亚洲精品| 白带黄色成豆腐渣| 日本一区二区免费在线视频| 亚洲精品久久成人aⅴ小说| 久久久久久久久中文| 母亲3免费完整高清在线观看| 国产99白浆流出| 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 亚洲中文字幕日韩| 亚洲av熟女| 婷婷亚洲欧美| 欧美黑人巨大hd| 日本成人三级电影网站| 欧美成人一区二区免费高清观看 | 看免费av毛片| 亚洲国产欧美人成| 亚洲午夜理论影院| 日本免费a在线| 成人国产综合亚洲| 色精品久久人妻99蜜桃| 国产亚洲精品久久久久久毛片| 精品欧美国产一区二区三| 两个人的视频大全免费| 中文字幕精品亚洲无线码一区| 久久国产乱子伦精品免费另类|