• <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ì)鈷的干擾
    中文欧美无线码| 9热在线视频观看99| 国产亚洲一区二区精品| 少妇被粗大的猛进出69影院| 电影成人av| 久久草成人影院| 久久青草综合色| 在线观看免费视频日本深夜| 精品国产美女av久久久久小说| 性色av乱码一区二区三区2| 如日韩欧美国产精品一区二区三区| 国内久久婷婷六月综合欲色啪| 欧美性长视频在线观看| 国产不卡一卡二| 国内毛片毛片毛片毛片毛片| 两人在一起打扑克的视频| 免费一级毛片在线播放高清视频 | 亚洲精品国产一区二区精华液| 黑人操中国人逼视频| 国产精品综合久久久久久久免费 | 国产无遮挡羞羞视频在线观看| 国产蜜桃级精品一区二区三区 | 精品国产超薄肉色丝袜足j| 80岁老熟妇乱子伦牲交| 日韩人妻精品一区2区三区| 成人三级做爰电影| 岛国在线观看网站| 在线观看www视频免费| 亚洲精品一卡2卡三卡4卡5卡| 黄频高清免费视频| 久久人妻福利社区极品人妻图片| 最近最新免费中文字幕在线| 高清av免费在线| 韩国av一区二区三区四区| 日韩免费av在线播放| 亚洲熟妇中文字幕五十中出 | 久久久精品国产亚洲av高清涩受| 9191精品国产免费久久| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 两性夫妻黄色片| 黑人猛操日本美女一级片| 国产人伦9x9x在线观看| 亚洲熟妇熟女久久| 亚洲九九香蕉| 午夜福利免费观看在线| 五月开心婷婷网| 亚洲av日韩在线播放| 大型av网站在线播放| 久久精品人人爽人人爽视色| 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻丝袜一区二区| 99精品欧美一区二区三区四区| 亚洲视频免费观看视频| 国产在线一区二区三区精| 两个人免费观看高清视频| 欧美日本中文国产一区发布| 亚洲av日韩在线播放| 欧美色视频一区免费| 亚洲第一欧美日韩一区二区三区| 欧美大码av| netflix在线观看网站| 一二三四在线观看免费中文在| 老司机靠b影院| videos熟女内射| 亚洲性夜色夜夜综合| 久久这里只有精品19| 亚洲国产欧美网| 十八禁高潮呻吟视频| 国产精品一区二区在线观看99| 久久久精品国产亚洲av高清涩受| 一本综合久久免费| 一级a爱片免费观看的视频| 欧美激情高清一区二区三区| 精品午夜福利视频在线观看一区| 麻豆成人av在线观看| 丰满的人妻完整版| 黄片大片在线免费观看| 免费一级毛片在线播放高清视频 | 亚洲国产欧美网| 女人高潮潮喷娇喘18禁视频| 久久中文字幕人妻熟女| 高清av免费在线| 国产精品99久久99久久久不卡| 他把我摸到了高潮在线观看| 久久九九热精品免费| 黄色丝袜av网址大全| 麻豆成人av在线观看| 免费在线观看完整版高清| 欧美日韩成人在线一区二区| 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 人人妻人人澡人人爽人人夜夜| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 欧美+亚洲+日韩+国产| 国产男女超爽视频在线观看| av网站在线播放免费| 久热爱精品视频在线9| 午夜免费观看网址| 一进一出抽搐动态| 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| cao死你这个sao货| 午夜免费鲁丝| 日本五十路高清| 香蕉丝袜av| 一级片'在线观看视频| 手机成人av网站| 国产av一区二区精品久久| 中亚洲国语对白在线视频| 国产片内射在线| 天堂俺去俺来也www色官网| 99热只有精品国产| 欧美日韩视频精品一区| 在线观看免费高清a一片| av天堂在线播放| 黄频高清免费视频| 国产精品永久免费网站| 国产无遮挡羞羞视频在线观看| 女人久久www免费人成看片| 日韩一卡2卡3卡4卡2021年| 高清欧美精品videossex| 精品一区二区三卡| 在线观看免费日韩欧美大片| 色综合婷婷激情| 亚洲国产精品sss在线观看 | 美女 人体艺术 gogo| 精品久久久精品久久久| 久久天堂一区二区三区四区| 日本vs欧美在线观看视频| 无限看片的www在线观看| 国产高清激情床上av| 欧美亚洲 丝袜 人妻 在线| 日日夜夜操网爽| 黑人巨大精品欧美一区二区蜜桃| 啪啪无遮挡十八禁网站| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美大码av| 三上悠亚av全集在线观看| 国产精品自产拍在线观看55亚洲 | 高清视频免费观看一区二区| 国产精品亚洲av一区麻豆| 老熟妇乱子伦视频在线观看| 亚洲,欧美精品.| 99热国产这里只有精品6| netflix在线观看网站| 国产不卡一卡二| 久久精品亚洲av国产电影网| 看片在线看免费视频| 黄网站色视频无遮挡免费观看| 亚洲国产毛片av蜜桃av| 亚洲欧美激情在线| 久久久久精品国产欧美久久久| 脱女人内裤的视频| 老熟女久久久| 免费看十八禁软件| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看| 首页视频小说图片口味搜索| 建设人人有责人人尽责人人享有的| 久久精品人人爽人人爽视色| 国产高清国产精品国产三级| 久久国产乱子伦精品免费另类| 夜夜爽天天搞| 少妇猛男粗大的猛烈进出视频| 国产精品久久久av美女十八| 男人操女人黄网站| 欧美日韩av久久| 91精品三级在线观看| 欧美丝袜亚洲另类 | 制服人妻中文乱码| 成人18禁高潮啪啪吃奶动态图| 中亚洲国语对白在线视频| 法律面前人人平等表现在哪些方面| 精品国产国语对白av| 热re99久久精品国产66热6| av在线播放免费不卡| 在线天堂中文资源库| 叶爱在线成人免费视频播放| 久久久精品国产亚洲av高清涩受| 黄片播放在线免费| 欧美不卡视频在线免费观看 | 久久精品国产综合久久久| ponron亚洲| 欧美日本中文国产一区发布| 十八禁人妻一区二区| 两人在一起打扑克的视频| 色婷婷久久久亚洲欧美| 19禁男女啪啪无遮挡网站| 色综合欧美亚洲国产小说| 不卡一级毛片| 18禁美女被吸乳视频| 久久人妻av系列| av线在线观看网站| 青草久久国产| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜一区二区| 欧美成人免费av一区二区三区 | 久久久久久久久久久久大奶| 精品国产国语对白av| 精品国产美女av久久久久小说| 精品无人区乱码1区二区| 激情在线观看视频在线高清 | 欧美国产精品va在线观看不卡| 美女高潮喷水抽搐中文字幕| 久久精品国产99精品国产亚洲性色 | 狠狠婷婷综合久久久久久88av| 色尼玛亚洲综合影院| 亚洲专区中文字幕在线| 一区二区三区精品91| 日韩免费高清中文字幕av| 深夜精品福利| 免费久久久久久久精品成人欧美视频| 97人妻天天添夜夜摸| 久久久久国产精品人妻aⅴ院 | 亚洲精品一卡2卡三卡4卡5卡| 午夜福利欧美成人| 亚洲av熟女| 高清av免费在线| 9热在线视频观看99| 国产免费av片在线观看野外av| 亚洲五月天丁香| 中文字幕最新亚洲高清| 亚洲熟女毛片儿| 99精品欧美一区二区三区四区| tocl精华| 国产精品电影一区二区三区 | 很黄的视频免费| 91在线观看av| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av | 另类亚洲欧美激情| 丝袜人妻中文字幕| 久久国产精品人妻蜜桃| 精品一区二区三区视频在线观看免费 | 欧美乱码精品一区二区三区| 九色亚洲精品在线播放| 免费黄频网站在线观看国产| 丝袜美腿诱惑在线| 亚洲免费av在线视频| 亚洲一区二区三区不卡视频| 欧美在线黄色| 国产成人影院久久av| 久久人人爽av亚洲精品天堂| 女人被躁到高潮嗷嗷叫费观| 国产欧美亚洲国产| 国产国语露脸激情在线看| netflix在线观看网站| 一级a爱片免费观看的视频| 麻豆av在线久日| 国产成人av激情在线播放| 亚洲国产欧美日韩在线播放| 国产亚洲欧美在线一区二区| 黄色丝袜av网址大全| 免费黄频网站在线观看国产| 免费在线观看亚洲国产| 少妇粗大呻吟视频| 久久精品成人免费网站| 欧美激情久久久久久爽电影 | 午夜福利,免费看| 少妇的丰满在线观看| 婷婷成人精品国产| 国产精品成人在线| 侵犯人妻中文字幕一二三四区| 精品国产一区二区三区久久久樱花| 人妻丰满熟妇av一区二区三区 | 成人国语在线视频| 成人特级黄色片久久久久久久| 女性生殖器流出的白浆| 麻豆国产av国片精品| 精品国产美女av久久久久小说| 夜夜爽天天搞| 亚洲avbb在线观看| 不卡一级毛片| 国产有黄有色有爽视频| 久久久久视频综合| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 别揉我奶头~嗯~啊~动态视频| 在线十欧美十亚洲十日本专区| 久久天堂一区二区三区四区| 欧美日韩一级在线毛片| 在线国产一区二区在线| 国产成人系列免费观看| 老鸭窝网址在线观看| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲人成77777在线视频| 波多野结衣一区麻豆| 午夜免费鲁丝| 日本五十路高清| 亚洲av欧美aⅴ国产| 精品高清国产在线一区| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 黄色毛片三级朝国网站| 久久久国产欧美日韩av| 人成视频在线观看免费观看| 我的亚洲天堂| 亚洲av成人不卡在线观看播放网| 亚洲七黄色美女视频| 国产精品欧美亚洲77777| 超碰97精品在线观看| 老司机影院毛片| 在线观看午夜福利视频| 国产精品久久久久成人av| av有码第一页| 99国产精品99久久久久| 久久久国产精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 精品人妻在线不人妻| 狂野欧美激情性xxxx| 丝瓜视频免费看黄片| 99香蕉大伊视频| 999久久久精品免费观看国产| a级毛片在线看网站| 欧美精品啪啪一区二区三区| 精品国产乱码久久久久久男人| 色综合婷婷激情| 大香蕉久久网| 国精品久久久久久国模美| 国产精品二区激情视频| 免费观看精品视频网站| 日韩三级视频一区二区三区| 日日摸夜夜添夜夜添小说| 深夜精品福利| 精品国内亚洲2022精品成人 | 国产精品欧美亚洲77777| 午夜福利欧美成人| 夫妻午夜视频| 国产欧美日韩一区二区三| 美女扒开内裤让男人捅视频| 久热这里只有精品99| 欧美成人免费av一区二区三区 | 精品无人区乱码1区二区| 三上悠亚av全集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜视频精品福利| 国产成人av教育| 久久久久久久国产电影| x7x7x7水蜜桃| av线在线观看网站| 久久久久久亚洲精品国产蜜桃av| 欧美精品高潮呻吟av久久| 国产精品亚洲av一区麻豆| videos熟女内射| 夜夜夜夜夜久久久久| 日韩欧美国产一区二区入口| aaaaa片日本免费| 自线自在国产av| 99久久国产精品久久久| 国产深夜福利视频在线观看| 国产区一区二久久| 午夜福利视频在线观看免费| 亚洲一区高清亚洲精品| 国产xxxxx性猛交| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 亚洲综合色网址| 精品国产一区二区三区四区第35| 夫妻午夜视频| 欧美激情 高清一区二区三区| 国产日韩一区二区三区精品不卡| 精品人妻熟女毛片av久久网站| 久久热在线av| 99re6热这里在线精品视频| 国产精品一区二区免费欧美| 日韩三级视频一区二区三区| 丝袜美足系列| 久久人妻熟女aⅴ| 久久精品亚洲熟妇少妇任你| 国产精品 国内视频| 9热在线视频观看99| 国产精品 国内视频| 变态另类成人亚洲欧美熟女 | 国产成人精品无人区| 色老头精品视频在线观看| 亚洲黑人精品在线| 一级黄色大片毛片| 如日韩欧美国产精品一区二区三区| 婷婷成人精品国产| 国产成人一区二区三区免费视频网站| av线在线观看网站| 可以免费在线观看a视频的电影网站| 在线观看一区二区三区激情| 日韩欧美在线二视频 | 日日爽夜夜爽网站| 国产男女超爽视频在线观看| 欧美国产精品va在线观看不卡| 大型av网站在线播放| 97人妻天天添夜夜摸| 这个男人来自地球电影免费观看| 一区二区三区激情视频| 99热国产这里只有精品6| 日韩制服丝袜自拍偷拍| 久久精品国产99精品国产亚洲性色 | 中文字幕最新亚洲高清| 十八禁网站免费在线| 久久 成人 亚洲| 久久精品人人爽人人爽视色| 国产免费男女视频| 中亚洲国语对白在线视频| 国产亚洲精品久久久久久毛片 | 极品少妇高潮喷水抽搐| 久久久久国产一级毛片高清牌| avwww免费| 久久ye,这里只有精品| 高清视频免费观看一区二区| 一区二区三区国产精品乱码| 男女之事视频高清在线观看| 九色亚洲精品在线播放| 国产av又大| 国产精品综合久久久久久久免费 | 国产有黄有色有爽视频| 成年版毛片免费区| 欧美精品一区二区免费开放| 国产深夜福利视频在线观看| 精品午夜福利视频在线观看一区| 午夜免费观看网址| 免费一级毛片在线播放高清视频 | 岛国毛片在线播放| 无限看片的www在线观看| cao死你这个sao货| 精品国产超薄肉色丝袜足j| 免费观看人在逋| 久久精品亚洲精品国产色婷小说| 亚洲av欧美aⅴ国产| 国产亚洲精品久久久久久毛片 | 国产在线观看jvid| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 亚洲 欧美一区二区三区| 成人免费观看视频高清| 热99国产精品久久久久久7| 免费人成视频x8x8入口观看| 精品国产超薄肉色丝袜足j| 国产精品美女特级片免费视频播放器 | 咕卡用的链子| 亚洲黑人精品在线| 婷婷成人精品国产| 女人精品久久久久毛片| 一区福利在线观看| xxxhd国产人妻xxx| 少妇裸体淫交视频免费看高清 | 男人舔女人的私密视频| 精品久久蜜臀av无| 久久天堂一区二区三区四区| 国产1区2区3区精品| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院 | 国产在线一区二区三区精| 欧美一级毛片孕妇| 免费在线观看完整版高清| 美女福利国产在线| 黄频高清免费视频| 天堂俺去俺来也www色官网| 女警被强在线播放| 午夜福利欧美成人| 50天的宝宝边吃奶边哭怎么回事| 黄色视频,在线免费观看| 在线十欧美十亚洲十日本专区| 在线观看午夜福利视频| 天天躁夜夜躁狠狠躁躁| 大片电影免费在线观看免费| 午夜福利,免费看| 中文字幕高清在线视频| 免费人成视频x8x8入口观看| 亚洲国产欧美日韩在线播放| 国产单亲对白刺激| 久久中文看片网| 热re99久久国产66热| 高清黄色对白视频在线免费看| 成人特级黄色片久久久久久久| 国产国语露脸激情在线看| 久久国产精品男人的天堂亚洲| 制服人妻中文乱码| 国产一区二区三区在线臀色熟女 | 日韩欧美一区二区三区在线观看 | 亚洲,欧美精品.| 男男h啪啪无遮挡| 国产成人精品无人区| 亚洲精品在线美女| 欧美黄色淫秽网站| 日韩人妻精品一区2区三区| 欧美亚洲日本最大视频资源| 久久久国产成人精品二区 | 日韩大码丰满熟妇| 乱人伦中国视频| 欧美中文综合在线视频| 欧美激情极品国产一区二区三区| 757午夜福利合集在线观看| 久久性视频一级片| 亚洲少妇的诱惑av| 久久ye,这里只有精品| 国产精品.久久久| 久久精品亚洲熟妇少妇任你| 欧美乱妇无乱码| av天堂在线播放| 在线播放国产精品三级| av有码第一页| 97人妻天天添夜夜摸| av福利片在线| 久久精品熟女亚洲av麻豆精品| 精品少妇久久久久久888优播| 黄色视频,在线免费观看| 一级a爱视频在线免费观看| 国产欧美日韩综合在线一区二区| 一级a爱片免费观看的视频| 黄片播放在线免费| 男女之事视频高清在线观看| 精品一区二区三区视频在线观看免费 | 国产色视频综合| 精品亚洲成a人片在线观看| 国产人伦9x9x在线观看| 黑丝袜美女国产一区| 国产精品1区2区在线观看. | 热re99久久精品国产66热6| 亚洲免费av在线视频| 欧美老熟妇乱子伦牲交| 搡老熟女国产l中国老女人| 久久精品91无色码中文字幕| 好男人电影高清在线观看| 亚洲avbb在线观看| 国产精品.久久久| 777久久人妻少妇嫩草av网站| tocl精华| 精品国产一区二区三区四区第35| 窝窝影院91人妻| 人成视频在线观看免费观看| 天天躁夜夜躁狠狠躁躁| 色尼玛亚洲综合影院| 午夜福利免费观看在线| 女同久久另类99精品国产91| 夜夜夜夜夜久久久久| 在线免费观看的www视频| 咕卡用的链子| 久久九九热精品免费| 国产精品久久久av美女十八| 脱女人内裤的视频| 欧美在线一区亚洲| 欧美+亚洲+日韩+国产| 高清在线国产一区| 久久中文字幕人妻熟女| 妹子高潮喷水视频| 老司机福利观看| 老汉色av国产亚洲站长工具| 一边摸一边抽搐一进一出视频| 国产极品粉嫩免费观看在线| 在线永久观看黄色视频| 香蕉久久夜色| 亚洲色图 男人天堂 中文字幕| 久久国产精品男人的天堂亚洲| 日日夜夜操网爽| av超薄肉色丝袜交足视频| 乱人伦中国视频| 国产日韩一区二区三区精品不卡| 狠狠婷婷综合久久久久久88av| 丝袜在线中文字幕| 法律面前人人平等表现在哪些方面| 久久久国产精品麻豆| 亚洲精品成人av观看孕妇| 性少妇av在线| 亚洲成av片中文字幕在线观看| 女性被躁到高潮视频| 女警被强在线播放| 1024香蕉在线观看| 国产在视频线精品| 欧美人与性动交α欧美精品济南到| 午夜福利欧美成人| 国产成人精品久久二区二区免费| 欧美日韩成人在线一区二区| 天天影视国产精品| 少妇猛男粗大的猛烈进出视频| 丝袜人妻中文字幕| 欧美精品高潮呻吟av久久| 亚洲精品av麻豆狂野| 精品福利永久在线观看| 午夜免费成人在线视频| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品一区二区三区在线| 伊人久久大香线蕉亚洲五| 最近最新中文字幕大全电影3 | 精品久久久精品久久久| 亚洲色图 男人天堂 中文字幕| 在线观看免费视频网站a站| 亚洲成人手机| a在线观看视频网站| 亚洲第一av免费看| 国产99久久九九免费精品| 亚洲精华国产精华精| 国产日韩欧美亚洲二区| 精品国内亚洲2022精品成人 | 又黄又爽又免费观看的视频| 三级毛片av免费| 人妻一区二区av| 9热在线视频观看99| 亚洲成a人片在线一区二区| 捣出白浆h1v1| 国产精品一区二区精品视频观看| 夜夜夜夜夜久久久久| 在线天堂中文资源库| 国产亚洲av高清不卡| 精品午夜福利视频在线观看一区| 色94色欧美一区二区| 日韩中文字幕欧美一区二区| 大片电影免费在线观看免费| 日韩视频一区二区在线观看| 精品一区二区三区视频在线观看免费 | 精品无人区乱码1区二区|