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

    南水北調(diào)中線干渠(河南段)浮游細(xì)菌群落組成及影響因素

    2017-06-05 08:42:09陳兆進(jìn)陳海燕李玉英魯開杰趙海軍朱靜亞胡蘭群
    中國環(huán)境科學(xué) 2017年4期
    關(guān)鍵詞:方城渠首浮游

    陳兆進(jìn),陳海燕,李玉英*,黃 進(jìn),魯開杰,趙海軍,李 冰,朱靜亞,胡蘭群

    ?

    南水北調(diào)中線干渠(河南段)浮游細(xì)菌群落組成及影響因素

    陳兆進(jìn)1,2,陳海燕2,3,李玉英1,2*,黃 進(jìn)2,3,魯開杰4,趙海軍2,3,李 冰2,3,朱靜亞1,2,胡蘭群2,3

    (1.南陽師范學(xué)院農(nóng)業(yè)工程學(xué)院,河南省南水北調(diào)中線水源區(qū)生態(tài)安全重點(diǎn)實(shí)驗(yàn)室,河南南陽473061;2.南水北調(diào)中線水源區(qū)水安全河南省協(xié)同創(chuàng)新中心,河南南陽473061;3.南水北調(diào)中線渠首環(huán)境監(jiān)測(cè)應(yīng)急中心,河南淅川474475;4.方城縣南水北調(diào)辦公室,河南方城 473200)

    為了解南水北調(diào)中線工程通水后干渠河南段浮游細(xì)菌群落組成及影響因素,于2016年7月對(duì)中線干渠河南段5個(gè)典型生態(tài)點(diǎn)位表層水樣進(jìn)行采集.水質(zhì)監(jiān)測(cè)結(jié)果表明干渠水質(zhì)總體較好,除總氮(TN)和高錳酸鹽指數(shù)(CODMn)外,其他各項(xiàng)指標(biāo)符合I類水標(biāo)準(zhǔn)要求.渠首、方城和新鄭3點(diǎn)位水體處于中營養(yǎng)狀態(tài),溫縣、安陽2點(diǎn)位水體處于輕度富營養(yǎng)狀態(tài).采用高通量測(cè)序技術(shù)對(duì)浮游細(xì)菌群落組成進(jìn)行研究,結(jié)果表明主要由放線菌門 (Actinobacteria)、變形菌門 (Proteobacteria)、浮霉菌門(Planctomycetes) 、擬桿菌門 (Bacteroidetes)、藍(lán)藻門 (Cyanobacteria)、疣微菌門 (Verrucomicrobia)、裝甲菌門(Armatimonadetes)等25個(gè)門和407個(gè)屬組成,5點(diǎn)位樣品具有豐富的群落組成,浮游細(xì)菌群落多樣性排序?yàn)?溫縣>新鄭>方城>渠首>安陽.對(duì)浮游細(xì)菌群落組成與環(huán)境因子的關(guān)系進(jìn)行典范對(duì)應(yīng)分析(CCA),結(jié)果表明: TN、Chl a、NH4+-N、COD、CODMn能顯著影響浮游細(xì)菌群落,pH值和TP對(duì)其也有一定程度影響.

    干渠;浮游細(xì)菌;群落結(jié)構(gòu);高通量測(cè)序;環(huán)境因子

    南水北調(diào)中線工程(以下簡稱“中線工程”)是一項(xiàng)跨流域、跨多省市、長距離的特大型調(diào)水工程,擔(dān)負(fù)著北京、天津、石家莊、鄭州等數(shù)十座城市保障供水的重大任務(wù).中線工程總干渠輸水形式以明渠為主,局部布置管涵[1].開放式渠道沒法完全避免外界污染物的侵入,如污染物由空氣進(jìn)入水體以及部分渠段地下水內(nèi)排與地下水污染物的滲漏等[2-3].因此,開展干渠河南段水質(zhì)監(jiān)測(cè)對(duì)確保南水北調(diào)中線供水水質(zhì)安全具有重要意義.

    浮游細(xì)菌是水生態(tài)系統(tǒng)的重要組成部分,在有機(jī)物降解及營養(yǎng)物循環(huán)過程中具有重要作用,而其組成變化對(duì)水環(huán)境也起到明顯的指示作用[4-5].因此,全面了解水生態(tài)系統(tǒng)中浮游細(xì)菌群落多樣性、分布特征及其在生態(tài)系統(tǒng)中的功能和作用,對(duì)于管理和維護(hù)水生態(tài)環(huán)境具有深遠(yuǎn)的意義[6].由于只有不到1%的微生物可以純培養(yǎng),限制了自然界中絕大部分微生物的認(rèn)識(shí).近年來,高通量測(cè)序技術(shù)(High throughput sequencing)的出現(xiàn),大大加速了宏基因組學(xué)的發(fā)展[7].該技術(shù)操作簡單、成本低廉,能快速高通量得到特定的 DNA 片段,較之前的分子生物學(xué)方法能更全面地展示生物群落的組成結(jié)構(gòu),極大地推動(dòng)了生物多樣性研究[8-13].

    中線工程干渠水質(zhì)狀況直接關(guān)系到受水區(qū)居民的飲水安全.本課題組已對(duì)中線工程調(diào)水前后不同時(shí)期庫體和干渠南陽段水體浮游藻類、浮游動(dòng)物進(jìn)行了監(jiān)測(cè)[14],但有關(guān)中線工程干渠浮游細(xì)菌群落研究尚未報(bào)道.本研究選取干渠河南段5個(gè)典型站位,對(duì)水體中總氮(TN)、總磷(TP)、化學(xué)需氧量(COD)、葉綠素a含量(Chl a)等水質(zhì)指標(biāo)進(jìn)行監(jiān)測(cè),同時(shí)采用高通量測(cè)序技術(shù)對(duì)浮游細(xì)菌群落結(jié)構(gòu)和多樣性進(jìn)行分析,研究其群落組成和分布特征.結(jié)合典范對(duì)應(yīng)分析(CCA)進(jìn)一步探討細(xì)菌豐度、群落結(jié)構(gòu)與環(huán)境因子間的關(guān)系,以期為中線工程干渠水環(huán)境保護(hù)提供參考依據(jù).

    1 材料與方法

    1.1 樣品采集與水質(zhì)理化指標(biāo)測(cè)定

    根據(jù)中線工程河南段干渠的地理位置特征,選取從南到北沿程的屬于陶岔渠首(樣品編號(hào)Q)、南陽市的方城獨(dú)樹(F)、鄭州市的新鄭十里鋪(X)、焦作市溫縣陳家溝(W)和安陽市安陽縣洪河屯(A)5個(gè)點(diǎn)位為采樣點(diǎn),具體信息見表1.參照《水和廢水監(jiān)測(cè)分析方法》(第四版)[15]于2016年7月用柱狀采水器采集表層(深度50cm)水樣2L,所有點(diǎn)位均取平行樣3組. 600mL水樣用于提取細(xì)菌總 DNA,置于預(yù)先滅菌容器中存于冰盒,運(yùn)回實(shí)驗(yàn)室后進(jìn)行檢測(cè)分析,剩余水樣用于水質(zhì)指標(biāo)測(cè)定.參照《水和廢水監(jiān)測(cè)分析方法》(第四版)[15]測(cè)定水樣pH值、總氮(TN)、氨氮(NH4+-N)、總磷(TP)、化學(xué)需氧量(COD)、高錳酸鹽指數(shù)(CODMn)以及葉綠素a含量(Chl a)7項(xiàng)理化指標(biāo).

    表1 采樣點(diǎn)分布

    1.2 樣品總 DNA提取

    600mL表層新鮮水樣經(jīng)20μm篩絹預(yù)過濾至滅菌燒杯中,以除去大型浮游植物和浮游動(dòng)物.預(yù)過濾后的水樣經(jīng)0.22μm無菌微孔濾膜過濾,用以收集浮游細(xì)菌,將濾膜剪碎置于50mL無菌離心管中.按照 Omega Water DNA Kit (Omega, USA)的說明,提取水體中的DNA.將提取得到的細(xì)菌總DNA通過微量紫外分光光度計(jì)(NanoDrop?ND-1000, Wilmington, DE, USA)測(cè)定DNA濃度和純度.

    1.3 高通量測(cè)序

    采用通用引物338F(5′-ACTCCTACGGGA- GGCAGCA-3′)和806R(5′-GGACTACHVGG- GTWTCTAAT-3′)對(duì)浮游細(xì)菌16S rRNA基因的V3~V4區(qū)擴(kuò)增,修飾后的通用引物含有不同的Tag標(biāo)簽用以區(qū)分不同樣品.PCR擴(kuò)增體系為20μL,其中含5×FastPfu Buffer 4μL、2.5mmol/L dNTPs 2μL、Forward Primer(5μmol/L) 0.8μL、Reverse Primer(5μmol/L) 0.8μL、FastPfu Polymerase 0.4μL、DNA模板10ng.補(bǔ)ddH2O至20μL.PCR擴(kuò)增的反應(yīng)條件為: 94℃,5min;30×(94℃,30s;54℃,30s;72℃,45s);72℃,10min.每個(gè)樣品3個(gè)重復(fù),利用上海美吉生物醫(yī)藥科技有限公司的MiSeq PE300測(cè)序儀(Illumina Inc..San Diego. CA.USA)完成序列測(cè)定.

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

    1.4.1 高通量數(shù)據(jù)分析 高通量數(shù)據(jù)的生物信息學(xué)分析采用 Qiime進(jìn)行,根據(jù)序列的相似度,將序列歸為多個(gè)OTU,OTU產(chǎn)出后,統(tǒng)計(jì)各個(gè)樣品含有OTU情況及每個(gè)OTU中含有序列的數(shù)目,得到每個(gè)OTU的分類學(xué)信息[16].選取相似度在97%條件下的OTU生成預(yù)期的稀釋曲線,并利用軟件mothur計(jì)算豐富度指數(shù)Chao1和ACE,覆蓋度指數(shù)以及多樣性指數(shù)Simpson和Shannon指數(shù)進(jìn)行Alpha多樣性分析[17].利用PCoA、Correlation聚類法分別進(jìn)行數(shù)據(jù)處理、細(xì)菌群落分布、主成分分析和聚類分析.采用Linear Discriminant Analysis(LDA) coupled with effect size measurements(LEfSe)在線工具尋找組與組之間有統(tǒng)計(jì)學(xué)差異的生物標(biāo)志物[18].

    1.4.2 典范對(duì)應(yīng)分析 典范對(duì)應(yīng)分析(CCA)能將可能相關(guān)的多個(gè)環(huán)境因子一起進(jìn)行分析,可以很好地反映種群或群落與環(huán)境的關(guān)系.本文采用CANOCO軟件(版本4.5)進(jìn)行浮游細(xì)菌和環(huán)境因子的CCA分析,排序結(jié)果用物種-環(huán)境因子關(guān)系的雙序圖表示[19].在主軸1和主軸2構(gòu)成的平面中,箭頭表示環(huán)境因子,向量長短代表著相應(yīng)環(huán)境因子在主軸中的作用,箭頭所處象限表示環(huán)境因子與排序軸間相關(guān)性的正負(fù).

    1.4.3 方差分析 不同采樣點(diǎn)位間的環(huán)境因子等數(shù)據(jù)采用 SPSS19.0軟件中的單因素方差分析,顯著性水平設(shè)定為<0.05,相關(guān)圖表制作在 Excel中完成.

    2 結(jié)果與分析

    2.1 水質(zhì)狀態(tài)評(píng)價(jià)

    表2 各采樣點(diǎn)水體主要理化參數(shù)

    注 :同一列不同小寫字母,表示處理在<0.05水平達(dá)到顯著.

    2016年7月干渠河南段渠首、方城等5點(diǎn)位水質(zhì)的理化監(jiān)測(cè)結(jié)果表明各水體的水質(zhì)總體較好,除TN (除渠首外,其他4點(diǎn)位超過1.00mg/L,為III類地表水標(biāo)準(zhǔn))和CODMn外(渠首樣品超過3.00mg/L,為II類地表水標(biāo)準(zhǔn)),其他各項(xiàng)指標(biāo)符合《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB38382- 2002)[20]I類水標(biāo)準(zhǔn)要求(表2).按照《地表水資源質(zhì)量評(píng)價(jià)技術(shù)規(guī)程》(SL395-2007)[21]進(jìn)行水庫營養(yǎng)狀態(tài)的評(píng)價(jià),以TN、TP、CODMn和Chl a 4個(gè)參數(shù)作為水質(zhì)營養(yǎng)狀態(tài)參數(shù),渠首、方城等5點(diǎn)位TN賦分值均介于60~70, TP、CODMn和Chl a的賦分值介于40~50之間.渠首、方城和新鄭3點(diǎn)位營養(yǎng)狀態(tài)指數(shù)EI分別為47.5、50.0和50.0,為中營養(yǎng)狀態(tài).溫縣、安陽樣品EI指數(shù)均為52.5,為輕度富營養(yǎng)狀態(tài)(表2).

    2.2 浮游細(xì)菌稀釋性曲線

    稀釋性曲線可以反映樣品的取樣深度,被用來評(píng)價(jià)測(cè)序量是否足以覆蓋所有類群.渠首、方城等5點(diǎn)位浮游細(xì)菌稀釋性曲線隨測(cè)序條數(shù)的增加物種豐富度呈現(xiàn)前期增加后期趨向平緩.5點(diǎn)位每個(gè)樣本的測(cè)序曲線在測(cè)序條數(shù)達(dá)到15000條以上時(shí)物種數(shù)都基本趨向平穩(wěn),說明此時(shí)測(cè)序數(shù)據(jù)量合理,真實(shí)環(huán)境中細(xì)菌群落結(jié)構(gòu)的置信度較高,能夠代表物種的豐富度.

    2.3 浮游細(xì)菌群落多樣性

    高通量測(cè)序表明渠首、方城等5點(diǎn)位樣品具有豐富的群落組成,其群落多樣性均較高(表3).5組樣品文庫覆蓋率均在98.90%以上,能夠反映該區(qū)域浮游細(xì)菌群落的種類和結(jié)構(gòu).溫縣和新鄭樣品OTU數(shù)、豐富度指數(shù)Chao1、Shannon指數(shù)高于其他3組樣品,同時(shí)Simpson指數(shù)和文庫覆蓋率低于其他組樣品,其中OTU數(shù)、文庫覆蓋率、Shannon指數(shù)和Simpson指數(shù)差異達(dá)到顯著水平(0.05),表明溫縣和新鄭樣品浮游細(xì)菌群落多樣性高于其他樣品(表3).渠首、方城和安陽3組樣品浮游細(xì)菌群落多樣性指數(shù)有所差異,但豐富度指數(shù)、多樣性指數(shù)等評(píng)估指標(biāo)未達(dá)顯著水平.綜合OTU數(shù)、豐富度指數(shù)Chao1和ACE以及多樣性指數(shù)Simpson和Shannon指數(shù),5組樣品浮游細(xì)菌群落多樣性排序?yàn)?溫縣>新鄭>方城>渠首>安陽.

    表3 浮游細(xì)菌群落多樣性評(píng)估表

    注 :同一列不同小寫字母,表示處理在0.05水平達(dá)到顯著.

    2.4 浮游細(xì)菌多樣性的PCoA分析(主坐標(biāo)分析)

    利用Qiime 軟件對(duì)渠首、方城等5點(diǎn)位樣品相關(guān)群落進(jìn)行了PCoA分析.如圖1所示,橫坐標(biāo) PC1貢獻(xiàn)度為50.67%,縱坐標(biāo) PC2貢獻(xiàn)度為14.27%,5點(diǎn)位中渠首、方城樣品聚集在PCoA分析圖右上側(cè),新鄭、溫縣樣品聚集在PCoA分析圖下側(cè),安陽樣品分布于PCoA分析圖左上側(cè). PCoA分析結(jié)果表明不同樣品來源浮游細(xì)菌群落結(jié)構(gòu)有所差異,其中渠首、方城浮游細(xì)菌群落組成較為相似,地理距離相近的新鄭、溫縣樣品浮游細(xì)菌群落組成也較為相似,安陽樣品與其他組樣品群落組成有所差異.后續(xù)使用基于unweighted unifrac的UPGMA 方法對(duì)渠首等5點(diǎn)位樣品浮游細(xì)菌群落構(gòu)成的相似性進(jìn)行聚類分析,同樣可以有效的將不同地理分布樣品實(shí)現(xiàn)組間聚類(圖2).

    2.5 浮游細(xì)菌群落組成

    高通量測(cè)序結(jié)果表明,15組實(shí)驗(yàn)樣品OTU平均數(shù)525(表3),經(jīng)過分析主要為細(xì)菌的 25個(gè)門,包括放線菌門(Actinobacteria)、變形菌門(Proteobacteria)、浮霉菌門(Planctomycetes)、擬桿菌門 (Bacteroidetes)、藍(lán)藻門(Cyanobacteria)、疣微菌門(Verrucomicrobia)、裝甲菌門(Armatimonadetes)、芽單胞菌門(Gemmatim- onadetes)等(圖3).

    其中渠首樣品包含14個(gè)門,方城、新鄭樣品包含20個(gè)門,溫縣樣品包含21個(gè)門,安陽樣品包含17個(gè)門.浮游細(xì)菌主要分布于407個(gè)屬,其中渠首樣品包括296個(gè)屬,方城樣品包括303個(gè)屬,新鄭樣品包含339個(gè)屬,溫縣樣品包括342個(gè)屬,安陽樣品包括293個(gè)屬.屬于放線菌門、變形菌門、浮霉菌門、擬桿菌門、藍(lán)藻門、擬桿菌門和疣微菌門的序列總和占全部序列的97.28%~98.34%,這些細(xì)菌為優(yōu)勢(shì)種群(圖3).其中,放線菌門為渠首樣品最優(yōu)勢(shì)種群,占總?cè)郝涞?0.64%,同時(shí)為方城、新鄭、溫縣和安陽樣品第二優(yōu)勢(shì)菌群,占比分別為20.31%、30.80%、20.57%和18.86%.變形菌門為方城、新鄭、溫縣和安陽樣品最優(yōu)勢(shì)種群和渠首樣品第二優(yōu)勢(shì)菌群,分別占總?cè)郝涞?1.70%、45.55%、50.28%、58.40%和34.03%.

    2.6 不同樣品差異細(xì)菌分析

    為了進(jìn)一步確定各組細(xì)菌種屬的豐富度差異,采用在線統(tǒng)計(jì)工具LEfSe來尋找宏基因組生物標(biāo)志物.該方法是建立在統(tǒng)計(jì)學(xué)差異、生物學(xué)一致性和相同檢驗(yàn)效能估計(jì)基礎(chǔ)之上的,通過分析高通量測(cè)序數(shù)據(jù)在每個(gè)樣品中的分類信息,在生物學(xué)上提供不同級(jí)別的可能的生物標(biāo)志物[18]. LEfSe法計(jì)算出不同樣品中的差異菌屬,所有菌在門、綱、目、科、屬水平的差異信息用餅形圖表示(圖4).在門的水平上,渠首樣品中3個(gè)門(放線菌門、芽單胞菌門和WCHB1-60)、方城樣品中2個(gè)門(變形菌門和裝甲菌門)、新鄭樣品中2個(gè)門(TM6)、溫縣樣品中1個(gè)門(酸桿菌門)和安陽樣品中1個(gè)門(擬桿菌門)的細(xì)菌存在顯著差異.在屬的水平上,渠首樣品中放線菌門的hgcI_clade等,變形菌門的、、LD28等,藍(lán)藻門的等9個(gè)門共37個(gè)屬的細(xì)菌存在顯著差異;方城樣品中變形菌門的、、、等,裝甲菌門的,放線菌門的,芽單胞菌門的等6個(gè)門共25個(gè)屬的細(xì)菌存在顯著差異;新鄭樣品中變形菌門的、、、等,放線菌門的、等,厚壁菌門的,擬桿菌門的、、等6個(gè)門共27個(gè)屬的細(xì)菌存在顯著差異;溫縣樣品中酸桿菌門的,擬桿菌門的、、等,變形菌門的、、、、等7個(gè)門共68個(gè)屬的細(xì)菌存在顯著差異;安陽樣品中放線菌門的,變形菌門的、、、等,擬桿菌門的等4個(gè)門共27個(gè)屬的細(xì)菌存在顯著差異.

    2.7 浮游細(xì)菌群落組成與環(huán)境因子的相關(guān)分析

    采用CCA分析水體環(huán)境因子與LEfSe分析篩選出的主要細(xì)菌之間的關(guān)系.由圖5可知,所分析的環(huán)境因子中Chl a、NH4+-N、TN、pH和TP與第1排序軸(AX1)正相關(guān),其中Chl a、NH4+-N、TN相關(guān)系數(shù)較大,呈顯著正相關(guān)(<0.05). COD、CODMn與AX1負(fù)相關(guān),其中CODMn達(dá)顯著水平. pH、COD、CODMn、NH4+-N和TP與第2排序軸 (AX2)正相關(guān),其中COD、CODMn相關(guān)系數(shù)都較大,達(dá)顯著水平(<0.05).TN和Chl a與第2排序軸 (AX2)負(fù)相關(guān), TN相關(guān)性系數(shù)為-0.78,呈顯著性負(fù)相關(guān).從第1、第2排序軸的相關(guān)性分析可以看出,TN、Chl a、NH4+-N、COD、CODMn能顯著影響浮游細(xì)菌群落,pH和TP對(duì)其也有一定影響.

    CODMn、COD在CCA分析圖中箭頭方向和長度較為相似,與TN箭頭相反,所以CODMn、COD與細(xì)菌群落相關(guān)性分析結(jié)果相似,與TN相反.相關(guān)性分析表明,10個(gè)門46個(gè)屬的細(xì)菌與CODMn、COD和TN呈顯著性相關(guān),其中變形菌門占23個(gè)屬,為最大比例組成.變形菌門中Betaproteobacteria和Deltaproteobacteria的細(xì)菌與CODMn、COD呈顯著正相關(guān),與TN顯著負(fù)相關(guān),包括、、、、OM27clade等10個(gè)屬細(xì)菌.于此類似,疣微菌門,擬桿菌門的、、藍(lán)藻門的、放線菌門的、、hgcI clade等均與CODMn、COD呈顯著正相關(guān), 與TN顯著負(fù)相關(guān).放線菌門的、與CODMn、COD呈顯著負(fù)相關(guān),與TN顯著正相關(guān).Chl a與NH4+-N在CCA分析圖中箭頭方向和長度較為相似,所以與細(xì)菌群落相關(guān)性分析也相似.相關(guān)性分析表明,11個(gè)門112個(gè)屬的細(xì)菌與Chl a和NH4+-N呈顯著性相關(guān),其中變形菌門和擬桿菌門為主要組成,分別占72和15屬.變形菌門中Alphaproteobacteria和Gammaproteobacteria的絕大多數(shù)細(xì)菌與Chl a和NH4+-N呈顯著正相關(guān),包括等34個(gè)屬細(xì)菌.浮霉菌門的SM1A02、的細(xì)菌與Chl a和NH4+-N呈顯著正相關(guān).藍(lán)藻門的芽單胞菌門的裝甲菌門的以及放線菌門的hgcI clade、、等與Chl a和NH4+-N呈顯著負(fù)相 關(guān).

    3 討論

    3.1 中線干渠河南段水質(zhì)測(cè)定

    丹江口水庫和干渠水質(zhì)狀況直接關(guān)系到受水區(qū)居民的飲水安全.長期監(jiān)測(cè)表明,丹江流域存在農(nóng)業(yè)面源污染,其中TN明顯超標(biāo),造成丹江口庫區(qū)水體的TN超標(biāo)現(xiàn)象嚴(yán)重[22-24].為了解通水干渠河南段水質(zhì)狀況及影響水質(zhì)的主要污染因子,本研究測(cè)定了通水后第二年豐水期(2016年7月)干渠河南段5個(gè)典型站位水質(zhì).結(jié)果表明TN (除渠首外,干渠4點(diǎn)位超過1.00mg/L,為III類地表水標(biāo)準(zhǔn))和CODMn(渠首樣品超過3.00mg/L,為II類地表水標(biāo)準(zhǔn))為庫區(qū)主要污染因子,該結(jié)果與朱媛媛等[23]關(guān)于丹江口庫區(qū)水質(zhì)主要受TN及CODMn等指標(biāo)的影響,且TN為丹江口庫區(qū)水質(zhì)主要限制因子的研究結(jié)果基本一致.

    以TN、TP、CODMn和Chl a 4個(gè)參數(shù)對(duì)水庫營養(yǎng)狀態(tài)進(jìn)行評(píng)價(jià),渠首和方城點(diǎn)位營養(yǎng)狀態(tài)指數(shù)EI為47.5和50.0,為中營養(yǎng)狀態(tài),與王晨溪等[14]關(guān)于中線工程調(diào)水前后不同時(shí)期丹江口庫區(qū)和干渠南陽段水質(zhì)監(jiān)測(cè)結(jié)果基本一致.干渠其他3點(diǎn)位營養(yǎng)狀態(tài)指數(shù)EI分別為50.0(新鄭,中營養(yǎng)狀態(tài),)和52.5(溫縣、安陽, 輕度富營養(yǎng)狀態(tài)),該結(jié)果表明,隨著輸水過程水體富營養(yǎng)化程度逐漸增強(qiáng).中線工程總干渠輸水形式以明渠為主,水體富營養(yǎng)化是明渠污染的主要形式.研究表明水體能依靠生物降解、重力沉淀作用凈化水質(zhì),郭慶園[25]對(duì)南水北調(diào)中線工程京石段水質(zhì)監(jiān)測(cè)表明,水體中NH4+-N經(jīng)過水體自凈和沿程降解作用(硝化作用),含量降低趨勢(shì)明顯.本實(shí)驗(yàn)結(jié)果表明干渠河南段隨著輸水過程水體富營養(yǎng)化程度逐漸增強(qiáng),分析可能有以下原因:(1)良好的水生生態(tài)系統(tǒng),對(duì)水質(zhì)凈化具有重要作用.比如,泥沙在天然水體的凈化過程中起著非常重要的作用,張繼育等[26]采集南水北調(diào)京石段干渠底泥模擬對(duì)水體鉻的去除.中線干渠河南段通水時(shí)間較短,可能尚未形成良好的水生生態(tài)系統(tǒng)(水生植物、水生動(dòng)物、底泥等),對(duì)水質(zhì)凈化能力較弱.(2)開放式渠道沒法完全避免外界污染物的侵入,如污染物由空氣進(jìn)入水體以及部分渠段地下水內(nèi)排與地下水污染物的滲漏等,這些均可能會(huì)引起水質(zhì)下降[2-3].本實(shí)驗(yàn)只是初步對(duì)干渠某一時(shí)間點(diǎn)進(jìn)行了監(jiān)測(cè),針對(duì)干渠水質(zhì)變化的確切原因和時(shí)空規(guī)律,需要對(duì)干渠水體及其水生生態(tài)系統(tǒng)等進(jìn)行長期監(jiān)測(cè)和系統(tǒng)研究.

    3.2 中線干渠河南段浮游細(xì)菌群落組成

    浮游細(xì)菌是水生態(tài)系統(tǒng)的重要組成部分,本課題組曾于2016年枯水期(1月)對(duì)丹江口庫區(qū)浮游細(xì)菌群落進(jìn)行測(cè)定,并分析了庫區(qū)浮游細(xì)菌群落與水溫、Chla、COD、BOD、DO等諸多生態(tài)因子的關(guān)系[27],但關(guān)于中線干渠浮游細(xì)菌群落組成及其影響因素研究鮮見報(bào)道,相關(guān)工作亟待開展.為了研究通水后中線干渠浮游細(xì)菌群落組成,本實(shí)驗(yàn)采用Illumina公司MiSeq測(cè)序平臺(tái)對(duì)渠首、方城等干渠河南段5點(diǎn)位浮游細(xì)菌群落組成進(jìn)行了研究,發(fā)現(xiàn)其主要由放線菌門、變形菌門、浮霉菌門、擬桿菌門、藍(lán)藻門、疣微菌門、裝甲菌門、芽單胞菌門等25個(gè)門組成,hgcI clade、CL500-29marine group、MWH-UniP1aquatic group、unclassified、等407個(gè)屬組成,表現(xiàn)出群落組成的豐富性,該結(jié)果全面地展示了干渠河南段的浮游細(xì)菌組成結(jié)構(gòu).同時(shí)比較發(fā)現(xiàn),干渠浮游細(xì)菌群落組成與之前研究的鄱陽湖、太湖、滇池等淡水湖泊[10,28-30]以及長江、南明河等河流[13,31-32]及其河口[33]等水體典型的細(xì)菌組成類群相似,這些典型的淡水細(xì)菌大都隸屬于變形菌門、藍(lán)細(xì)菌門、擬桿菌門、放線菌門和疣微菌門等[6,29],這些細(xì)菌在湖泊、河流的元素循環(huán)和能量流動(dòng)中起到重要作用.本次實(shí)驗(yàn)結(jié)果揭示了中線干渠河南段的浮游細(xì)菌種群類別,為中線工程水質(zhì)保護(hù)提供重要的依據(jù).

    3.3 浮游細(xì)菌群落與環(huán)境因子的關(guān)系

    研究表明,水體中的環(huán)境因子,包括形態(tài)特征(例如:水體大小和深度)、物理化學(xué)特征(例如:溫度、pH、鹽度、無機(jī)營養(yǎng)鹽等)、有機(jī)質(zhì)的濃度和類型以及食物網(wǎng)的組成和物種間的相互作用等均能影響浮游細(xì)菌的群落組成[6].本實(shí)驗(yàn)采用在線統(tǒng)計(jì)工具LEfSe分析各組細(xì)菌種屬的豐富度差異,同時(shí)對(duì)差異浮游細(xì)菌群落與理化指標(biāo)進(jìn)行CCA分析.結(jié)果表明,TN、NH4+-N、Chl a、COD、CODMn能顯著影響浮游細(xì)菌群落.丹江流域存在農(nóng)業(yè)面源污染,其中總氮明顯超標(biāo)[22,34],渠首TN超過0.5mg/L,為II類地表水標(biāo)準(zhǔn),干渠點(diǎn)位樣品TN顯著高于渠首樣品(<0.05),其中溫縣、安陽樣品NH4+-N差異達(dá)極顯著水平(< 0.01),表明隨著輸水過程水體TN濃度逐漸增強(qiáng)(表2).在各種水生態(tài)環(huán)境中,溶解無機(jī)氮是異養(yǎng)細(xì)菌重要的氮源,同時(shí)過量的氮可能導(dǎo)致水體富營養(yǎng)化,造成水質(zhì)惡化.Haukka等[35]報(bào)道在淡水生態(tài)系統(tǒng)中添加總氮和磷酸鹽后,浮游細(xì)菌的群落結(jié)構(gòu)發(fā)生變化,這種變化可能與生態(tài)系統(tǒng)中添加營養(yǎng)物質(zhì)后引起的水體富營養(yǎng)化有關(guān).干渠河南段水質(zhì)監(jiān)測(cè)結(jié)果表明渠首、方城和新鄭點(diǎn)位為中營養(yǎng)狀態(tài),溫縣和安陽為輕度富營養(yǎng)狀態(tài),隨著輸水過程水體富營養(yǎng)化程度逐漸增強(qiáng),從而影響干渠浮游細(xì)菌群落組成.

    4 結(jié)論

    4.1 2016年7月中線干渠河南段渠首、方城等5點(diǎn)位監(jiān)測(cè)結(jié)果表明水體的水質(zhì)總體較好,除TN和CODMn外,其他各項(xiàng)指標(biāo)符合I類水標(biāo)準(zhǔn)要求.以TN、TP、CODMn和Chl a 4個(gè)參數(shù)作為水質(zhì)營養(yǎng)狀態(tài)參數(shù),渠首、方城和新鄭3點(diǎn)位水體處于中營養(yǎng)狀態(tài),溫縣、安陽2點(diǎn)位水體處于輕度富營養(yǎng)狀態(tài).

    4.2 高通量測(cè)序表明丹江口水庫具有豐富的浮游細(xì)菌群落組成,主要由放線菌門、變形菌門、浮霉菌門、擬桿菌門、藍(lán)藻門、裝甲菌門和疣微菌門等25個(gè)門和407個(gè)屬的細(xì)菌組成.浮游細(xì)菌群落多樣性排序?yàn)?溫縣>新鄭>方城>渠首>安陽.

    4.3 浮游細(xì)菌與環(huán)境因子CCA分析表明,TN、NH4+-N、Chl a、COD、CODMn能顯著影響浮游細(xì)菌群落,pH和TP對(duì)其也有一定程度影響.

    [1] 國家計(jì)委,水利部.南水北調(diào)工程總體規(guī)劃 [R]. 北京:國家發(fā)展計(jì)劃委員會(huì), 2002.

    [2] 王興偉,陳家軍.地下水內(nèi)排對(duì)南水北調(diào)中線干渠水質(zhì)的影響 [J]. 南水北調(diào)與水利科技, 2015,13(5):858-861.

    [3] 高國軍.南水北調(diào)中線工程北京段水質(zhì)分析及其預(yù)測(cè)研究 [D]. 北京:北京林業(yè)大學(xué), 2016.

    [4] 黃 藝,舒中亞.基于浮游細(xì)菌生物完整性指數(shù)的河流生態(tài)系統(tǒng)健康評(píng)價(jià)-以滇池流域?yàn)槔?[J]. 環(huán)境科學(xué), 2013,34(8):3010-3018.

    [5] Newton R J, Jones S E, Eiler A, et al. A guide to the natural history of freshwater lake bacteria [J]. Microbiology and Molecular Biology Reviews, 2011,75(1):14-49.

    [6] 任麗娟,何 聃,邢 鵬,等.湖泊水體細(xì)菌多樣性及其生態(tài)功能研究進(jìn)展 [J]. 生物多樣性, 2013,21(4):421-432.

    [7] van Dijk E L, Auger H, Jaszczyszyn Y, et al. Ten years of next-generation sequencing technology [J]. Trends in Genetics, 2014,30(9):418-426.

    [8] Fan L M, Barry K, Hu G D, et al. Bacterioplankton community analysis in tilapia ponds by Illumina high-throughput sequencing [J]. World Journal of Microbiology and Biotechnology, 2015,32(1):1-11.

    [9] 白 潔,劉小沙,侯 瑞,等.南海南部海域浮游細(xì)菌群落特征及影響因素研究 [J]. 中國環(huán)境科學(xué), 2014,34(11):2950-2957.

    [10] 寇文伯,黃正云,張 杰,等.鄱陽湖湖泊細(xì)菌群落組成及結(jié)構(gòu)-以松門山為例 [J]. 生態(tài)學(xué)報(bào), 2015,35(23):7608-7614.

    [11] Yu L Y, Zhang W J, Liu L M, et al. Determining microeukaryotic plankton community around xiamen island, southeast china, using Illumina MiSeq and PCR-DGGE techniques [J]. PLoS ONE, 2015,10(5):e0127721.

    [12] Liu L M, Yang J, Yu Z, et al. The biogeography of abundant and rare bacterioplankton in the lakes and reservoirs of China [J]. ISME Journal, 2015:2068–2077.

    [13] 毛 莉,張 明,白 芳,等.三峽庫區(qū)支流藍(lán)藻水華對(duì)浮游細(xì)菌群落組成的影響 [J]. 水生生物學(xué)報(bào), 2016,40(3):609-614.

    [14] 王晨溪,朱靜亞,牛其愷,等.丹江口水庫和干渠南陽段微型生物群落的周期性變化 [J]. 安徽師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016,39(2):150-156.

    [15] 國家環(huán)境保護(hù)總局.水和廢水監(jiān)測(cè)分析方法(第四版) [M]. 北京:中國環(huán)境科學(xué)出版社, 2002:239-284.

    [16] Caporaso J G, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data [J]. Nature Methods, 2010,7(5):335-336.

    [17] Schloss P D, Westcott S L, Ryabin T, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities [J]. Applied and Environmental Microbiology, 2009,75(23):7537-7541.

    [18] Segata N, Izard J, Waldron L, et al. Metagenomic biomarker discovery and explanation [J]. Genome Biology, 2011,12(6):1-18.

    [19] Ter Braak C J, Smilauer P. CANOCO reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5) [M]. Microcomputer Power, Ithaca, NY, 2002.

    [20] GB/3838-2002 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)[S].

    [21] SL395-2007 地表水資源質(zhì)量評(píng)價(jià)技術(shù)規(guī)程[S].

    [22] 徐國策,李占斌,李 鵬,等.丹江中游典型小流域土壤總氮的空間分布 [J]. 地理學(xué)報(bào), 2012,67(11):1547-1555.

    [23] 朱媛媛,田進(jìn)軍,李紅亮,等.丹江口水庫水質(zhì)評(píng)價(jià)及水污染特征 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2016,35(1):139-147.

    [24] 陳 靜,丁衛(wèi)東,焦 飛,等.丹江口水庫總氮含量較高的調(diào)查分析 [J]. 中國環(huán)境監(jiān)測(cè), 2005,21(3):54-57.

    [25] 郭慶園.南水北調(diào)京石段水質(zhì)遷移轉(zhuǎn)化規(guī)律研究 [D]. 青島理工大學(xué), 2011.

    [26] 張繼育,王劭然,張大磊.南水北調(diào)中線工程京石段干渠底泥吸附鉻研究 [J]. 青島理工大學(xué)學(xué)報(bào), 2013,34(6):60-65.

    [27] 陳兆進(jìn),丁傳雨,朱靜亞,等.丹江口水庫枯水期浮游細(xì)菌群落組成及影響因素研究 [J]. 中國環(huán)境科學(xué), 2017,37(1):336-344.

    [28] Niu Y, Shen H, Chen J, et al. Phytoplankton community succession shaping bacterioplankton community composition in Lake Taihu, China [J]. Water Research, 2011,45(14):4169-4182.

    [29] Dai Y, Yang Y Y, Wu Z, et al. Spatiotemporal variation of planktonic and sediment bacterial assemblages in two plateau freshwater lakes at different trophic status [J]. Applied Microbiology and Biotechnology, 2016,100(9):4161-4175.

    [30] 呂明姬,汪杰,范 錚,等.滇池浮游細(xì)菌群落組成的空間分布特征及其與環(huán)境因子的關(guān)系 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2011,31(2):299-306.

    [31] 唐 婧,徐小蓉,商傳禹,等.南明河城區(qū)河段細(xì)菌多樣性與環(huán)境因子的關(guān)系 [J]. 微生物學(xué)報(bào), 2015,55(8):1050-1059.

    [32] 王 鵬,陳 波,李傳瓊,等.贛江南昌段豐水期細(xì)菌群落特征 [J]. 中國環(huán)境科學(xué), 2016,36(8):2453-2462.

    [33] 闞金軍,孫 軍.河口細(xì)菌群落多樣性及其控制因素:以切薩皮克灣為例 [J]. 生物多樣性, 2011,19(6):770-778.

    [34] 肖春艷,俐 武,趙同謙,等.南水北調(diào)中線源頭區(qū)蓄水前土壤氮磷分布特征 [J]. 中國環(huán)境科學(xué), 2013,33(10):1814-1820.

    [35] Haukka K, Kolmonen E, Hyder R, et al. Effect of nutrient loading on bacterioplankton community composition in lake Mesocosms [J]. Microbial Ecology, 2006,51(2):137-146.

    Community structure and influencing factors of bacterioplankton in the Main Cancel of the Mid-line Project of South-to-North Water Division in sections of Henan Province.

    CHEN Zhao-jin1,2, CHEN Hai-yan2,3, LI Yu-ying1,2*, HUANG Jin2,3, LU Kai-jie4, ZHAO Hai-jun2,3,LI Bing2,3, ZHU Jing-ya1,2,HU Lan-qun2,3

    (1.Key Laboratory of Ecological Security for Water Source Region of Mid-line Project of South-to-North Diversion Project of Henan Province, College of Agricultural Engineering, Nanyang Normal University, Nanyang 473061, China;2.Henan Collaborative Innovation Center of Water Security for Water Source Region of Mid-line Project of South-to-North Diversion Project, Nanyang 473061, China;3.Emergency Centre for Environmental Monitoring of Mid-line Project of South to North Water Division, Xichuan 474475, China;4.Office of South to North Water Diversion of Fangcheng County, Fangcheng 473200, China).

    To determine which factors may have an impact on the community composition of bacterioplankton, 15surface water samples were collected from the main cancel of the Mid-line Project of South-to-North Water Division in 5sections of Henan province, in July 2016, and 7essential environmental factors were evaluated. The results showed that water quality in the main cancel was strongly influenced by total nitrogen (TN) and the permanganate index (CODMn). Water system in Qushou, Fangcheng, and Xinzheng sections was under the mesotrophic state, and in Wenxian and Anyang sections, under the light eutrophic state. High-throughput sequencing was used to analyze distribution characteristics of the community structure, diversity of bacterioplankton. Phylogenetic analysis based on 16S ribosomal DNA sequences revealed that bacterioplankton could be subdivided into 25major phylogenetic and 407genus groups. The dominant phylogenetic groups included Actinobacteria, Proteobacteria, Planctomycetes, Bacteroidetes, Cyanobacteria, Verrucomicrobia, and Armatimonadetes. Diversity of the bacterioplankton community in the main cancel of the Mid-line Project of South-to-North Water Division could be ranked in the following order: Wenxian >Xinzheng > Fangcheng > Qushou > Anyang. The results of canonical correspondence analysis (CCA) of the relation between the bacterioplankton and 7 environmental factors showed that TN, ammonia nitrogen (NH4+-N), chlorophyll a (Chl a), chemical oxygen demand (COD), and the CODMnwere the main environmental factors affecting the distribution characteristics of the bacterioplankton community.

    main cancel;bacterioplankton;community structure;high-throughput sequencing;environmental factor

    172

    A

    1000-6923(2017)04-1505-09

    2016-09-18

    國家自然科學(xué)基金項(xiàng)目(41601332);河南省教育廳高等學(xué)校重點(diǎn)科研項(xiàng)目(16A210012);南陽師范學(xué)院高層次人才科研啟動(dòng)費(fèi)資助項(xiàng)目(ZX2014069)

    陳兆進(jìn)(1985-),男,安徽合肥人,講師,博士,主要從事土水污染控制、修復(fù)與資源化利用研究.發(fā)表論文10余篇.

    * 責(zé)任作者, 教授, lyying200508@163.com

    , 2017,37(4):1505~1513

    猜你喜歡
    方城渠首浮游
    昆侖渠首改建水力模型試驗(yàn)方案研究
    浞河浮游藻類的調(diào)查研究與水質(zhì)評(píng)價(jià)
    從楚方城到萬里長城
    聚力碧水攻堅(jiān) 護(hù)航美麗方城
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    ROARINGWILD壹方城店
    人大專題詢問之“方城樣本”
    應(yīng)用三維動(dòng)力有限元法對(duì)新疆葉爾羌河中游渠首工程閘室抗震穩(wěn)定性分析
    浮游藻類在水體PAHs富集與降解中的研究
    黑龍江省灌區(qū)渠首現(xiàn)狀及發(fā)展趨勢(shì)
    亚洲久久久国产精品| 另类精品久久| 蜜桃国产av成人99| 亚洲人成网站在线观看播放| 成人免费观看视频高清| 这个男人来自地球电影免费观看 | 免费观看性生交大片5| 日本欧美国产在线视频| 精品亚洲乱码少妇综合久久| 色94色欧美一区二区| 亚洲精品久久成人aⅴ小说 | 国产一区有黄有色的免费视频| a级毛片在线看网站| 伦理电影免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费现黄频在线看| 精品亚洲乱码少妇综合久久| a级毛色黄片| 国产在视频线精品| av天堂久久9| 日本wwww免费看| 全区人妻精品视频| 日韩精品免费视频一区二区三区 | 91午夜精品亚洲一区二区三区| 蜜桃在线观看..| 国产精品久久久久久久久免| 色视频在线一区二区三区| 国产一区二区在线观看日韩| 五月开心婷婷网| 精品亚洲乱码少妇综合久久| 99久久精品国产国产毛片| 欧美精品国产亚洲| 高清黄色对白视频在线免费看| 日本91视频免费播放| 王馨瑶露胸无遮挡在线观看| 美女国产视频在线观看| 日韩av不卡免费在线播放| 亚洲国产欧美日韩在线播放| 久久女婷五月综合色啪小说| 蜜桃国产av成人99| 免费大片黄手机在线观看| 韩国av在线不卡| 最近中文字幕高清免费大全6| 亚洲色图 男人天堂 中文字幕 | 精品久久久久久久久av| 国产亚洲欧美精品永久| 日本爱情动作片www.在线观看| 99久久综合免费| 老司机影院成人| 韩国av在线不卡| 久久久久久久大尺度免费视频| 欧美 日韩 精品 国产| 精品亚洲乱码少妇综合久久| 午夜影院在线不卡| 国产成人午夜福利电影在线观看| 久久久久久久久久成人| 精品酒店卫生间| 人体艺术视频欧美日本| 日本免费在线观看一区| 制服丝袜香蕉在线| 久久免费观看电影| 蜜臀久久99精品久久宅男| 国产黄频视频在线观看| 少妇丰满av| 亚洲精品色激情综合| 国产精品人妻久久久影院| 2018国产大陆天天弄谢| 色哟哟·www| 欧美性感艳星| 啦啦啦视频在线资源免费观看| 国语对白做爰xxxⅹ性视频网站| 女性被躁到高潮视频| 一边摸一边做爽爽视频免费| 99热这里只有是精品在线观看| 久久久久久久大尺度免费视频| 国产高清三级在线| 日本欧美视频一区| h视频一区二区三区| 日本黄色日本黄色录像| 国产一区二区三区综合在线观看 | 亚洲av在线观看美女高潮| 久久久a久久爽久久v久久| 午夜免费观看性视频| 欧美+日韩+精品| 日韩熟女老妇一区二区性免费视频| 成年美女黄网站色视频大全免费 | av在线app专区| 少妇人妻精品综合一区二区| 免费黄网站久久成人精品| 久久久久久久亚洲中文字幕| 两个人免费观看高清视频| 人妻制服诱惑在线中文字幕| 最近中文字幕2019免费版| 亚洲欧美精品自产自拍| 久久免费观看电影| 成人黄色视频免费在线看| 一本—道久久a久久精品蜜桃钙片| 免费av中文字幕在线| 国产高清有码在线观看视频| 大片免费播放器 马上看| videossex国产| 亚洲精品日本国产第一区| 国产精品99久久99久久久不卡 | 最近2019中文字幕mv第一页| 一级黄片播放器| 免费黄网站久久成人精品| 美女中出高潮动态图| 18+在线观看网站| 日日摸夜夜添夜夜添av毛片| 五月开心婷婷网| 寂寞人妻少妇视频99o| av网站免费在线观看视频| 免费高清在线观看视频在线观看| 久久99蜜桃精品久久| 九九久久精品国产亚洲av麻豆| 啦啦啦在线观看免费高清www| 国产日韩一区二区三区精品不卡 | 99热全是精品| 亚洲美女搞黄在线观看| 日韩视频在线欧美| 国产成人一区二区在线| 欧美激情 高清一区二区三区| 免费av不卡在线播放| 成人毛片a级毛片在线播放| 久久久午夜欧美精品| 国产午夜精品久久久久久一区二区三区| 国产精品一国产av| 国产伦精品一区二区三区视频9| 爱豆传媒免费全集在线观看| 亚洲欧洲精品一区二区精品久久久 | 九九久久精品国产亚洲av麻豆| 日本wwww免费看| 一区二区三区精品91| 极品人妻少妇av视频| 亚洲av电影在线观看一区二区三区| 性色avwww在线观看| 麻豆乱淫一区二区| 九九爱精品视频在线观看| 男人添女人高潮全过程视频| 久久久久久久久久成人| 久久99一区二区三区| 26uuu在线亚洲综合色| 自拍欧美九色日韩亚洲蝌蚪91| 视频中文字幕在线观看| 久久国产亚洲av麻豆专区| 国产精品久久久久成人av| 99热这里只有是精品在线观看| 亚洲美女搞黄在线观看| 99久久精品一区二区三区| 欧美激情国产日韩精品一区| 九九在线视频观看精品| 亚洲一区二区三区欧美精品| 国产精品秋霞免费鲁丝片| 午夜激情久久久久久久| 精品少妇久久久久久888优播| 亚洲av免费高清在线观看| 91精品国产九色| 久热这里只有精品99| 成人亚洲欧美一区二区av| 久久久久久久久久久久大奶| 国产一区二区在线观看日韩| 黑人猛操日本美女一级片| 日韩欧美精品免费久久| 丰满饥渴人妻一区二区三| 99久久人妻综合| 在线精品无人区一区二区三| 母亲3免费完整高清在线观看 | 9色porny在线观看| 一本大道久久a久久精品| 纯流量卡能插随身wifi吗| 美女中出高潮动态图| 色94色欧美一区二区| 久久99精品国语久久久| 亚洲欧洲国产日韩| 久久久久人妻精品一区果冻| 亚洲高清免费不卡视频| 51国产日韩欧美| 久久久a久久爽久久v久久| 少妇高潮的动态图| 亚洲成色77777| 久久免费观看电影| 欧美bdsm另类| 精品久久国产蜜桃| 18禁在线无遮挡免费观看视频| 午夜激情福利司机影院| 久久精品国产鲁丝片午夜精品| 精品99又大又爽又粗少妇毛片| 国产淫语在线视频| 成年人免费黄色播放视频| 九九在线视频观看精品| 26uuu在线亚洲综合色| 国产精品欧美亚洲77777| 久久午夜福利片| 如何舔出高潮| 亚洲精品一二三| 午夜福利在线观看免费完整高清在| 少妇人妻 视频| 最近中文字幕高清免费大全6| 最近的中文字幕免费完整| 少妇精品久久久久久久| 国产在线视频一区二区| 国产极品天堂在线| 久久人人爽人人爽人人片va| 18禁在线播放成人免费| 人妻少妇偷人精品九色| 伊人久久国产一区二区| 国产精品一区二区三区四区免费观看| 大香蕉久久成人网| 男女国产视频网站| 丝袜脚勾引网站| 欧美一级a爱片免费观看看| 免费av不卡在线播放| 国产 一区精品| 久久ye,这里只有精品| 午夜激情av网站| 青春草视频在线免费观看| 成人亚洲欧美一区二区av| 日韩中字成人| 久久综合国产亚洲精品| 男人爽女人下面视频在线观看| 香蕉精品网在线| 亚洲精华国产精华液的使用体验| 中文字幕人妻丝袜制服| 国产极品粉嫩免费观看在线 | xxxhd国产人妻xxx| 熟女人妻精品中文字幕| av在线app专区| 精品一区二区三卡| 久久国产精品大桥未久av| 在线观看免费日韩欧美大片 | 中文字幕制服av| 少妇被粗大的猛进出69影院 | 亚洲欧美一区二区三区国产| 精品99又大又爽又粗少妇毛片| 亚洲,欧美,日韩| 国产精品 国内视频| 少妇人妻久久综合中文| a级毛片免费高清观看在线播放| 国产亚洲精品久久久com| av黄色大香蕉| 人妻制服诱惑在线中文字幕| 亚洲av国产av综合av卡| 国产男女内射视频| 曰老女人黄片| 狠狠精品人妻久久久久久综合| 欧美激情 高清一区二区三区| 亚洲精品视频女| 日韩中文字幕视频在线看片| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 欧美精品高潮呻吟av久久| 久久久久国产网址| 亚洲熟女精品中文字幕| 日韩成人av中文字幕在线观看| 欧美精品亚洲一区二区| 夜夜爽夜夜爽视频| 国产精品一区二区三区四区免费观看| 国产黄频视频在线观看| 国产片内射在线| 国产精品99久久99久久久不卡 | 国产亚洲av片在线观看秒播厂| av有码第一页| 国产精品偷伦视频观看了| 夫妻午夜视频| 简卡轻食公司| 日韩视频在线欧美| 国产免费福利视频在线观看| 国产成人91sexporn| 在线观看国产h片| 午夜激情av网站| 久久久a久久爽久久v久久| 国产毛片在线视频| 久久热精品热| 亚洲国产欧美在线一区| 久久 成人 亚洲| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 亚洲av中文av极速乱| av一本久久久久| 夜夜骑夜夜射夜夜干| 国产乱人偷精品视频| 国产在线免费精品| 3wmmmm亚洲av在线观看| xxxhd国产人妻xxx| 如何舔出高潮| 亚洲国产日韩一区二区| 欧美精品一区二区免费开放| 在线观看国产h片| 51国产日韩欧美| 考比视频在线观看| 母亲3免费完整高清在线观看 | 91久久精品国产一区二区三区| 美女内射精品一级片tv| a级毛色黄片| 成人综合一区亚洲| 天堂中文最新版在线下载| 自拍欧美九色日韩亚洲蝌蚪91| 老司机亚洲免费影院| 亚洲欧美色中文字幕在线| 亚洲欧美精品自产自拍| 欧美精品国产亚洲| 亚洲怡红院男人天堂| 亚洲精品av麻豆狂野| 亚洲国产成人一精品久久久| 伊人久久精品亚洲午夜| 亚洲av国产av综合av卡| 黑人猛操日本美女一级片| 水蜜桃什么品种好| 成人午夜精彩视频在线观看| 日韩av不卡免费在线播放| 丰满乱子伦码专区| 亚洲国产色片| 午夜视频国产福利| 日韩av在线免费看完整版不卡| 亚洲欧洲精品一区二区精品久久久 | 午夜福利视频在线观看免费| 最近中文字幕2019免费版| 99热网站在线观看| 美女cb高潮喷水在线观看| 午夜福利视频在线观看免费| 你懂的网址亚洲精品在线观看| 十分钟在线观看高清视频www| 亚洲综合精品二区| 国产精品麻豆人妻色哟哟久久| 亚洲经典国产精华液单| 99热网站在线观看| 欧美日韩国产mv在线观看视频| videos熟女内射| 日本黄大片高清| 国产亚洲欧美精品永久| 色婷婷久久久亚洲欧美| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 精品久久国产蜜桃| 韩国av在线不卡| 国产毛片在线视频| 国产极品粉嫩免费观看在线 | 欧美日韩亚洲高清精品| 亚洲av二区三区四区| 日韩三级伦理在线观看| 熟女人妻精品中文字幕| 久久久久久久久久成人| 日本av免费视频播放| 亚洲色图综合在线观看| 精品人妻一区二区三区麻豆| 国产熟女欧美一区二区| 国产伦精品一区二区三区视频9| 国内精品宾馆在线| 国产成人午夜福利电影在线观看| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 999精品在线视频| 夫妻性生交免费视频一级片| xxx大片免费视频| 少妇 在线观看| 久久精品国产亚洲av涩爱| 亚洲欧美成人综合另类久久久| 免费观看的影片在线观看| 赤兔流量卡办理| 如日韩欧美国产精品一区二区三区 | 街头女战士在线观看网站| 丰满乱子伦码专区| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 国产亚洲午夜精品一区二区久久| 五月天丁香电影| 亚洲无线观看免费| 亚洲美女视频黄频| 国内精品宾馆在线| 80岁老熟妇乱子伦牲交| 免费观看a级毛片全部| 久久久久久久精品精品| 精品人妻一区二区三区麻豆| 久久99一区二区三区| 男女边吃奶边做爰视频| 国产精品国产三级国产专区5o| 久久鲁丝午夜福利片| 草草在线视频免费看| 最后的刺客免费高清国语| 国产色爽女视频免费观看| 美女脱内裤让男人舔精品视频| 777米奇影视久久| 国产精品熟女久久久久浪| 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 好男人视频免费观看在线| 亚洲熟女精品中文字幕| 黄片播放在线免费| 欧美成人午夜免费资源| 国产免费一级a男人的天堂| 一边亲一边摸免费视频| 成人二区视频| 老司机影院毛片| 国产精品三级大全| 欧美亚洲日本最大视频资源| 欧美最新免费一区二区三区| xxx大片免费视频| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| xxx大片免费视频| 少妇熟女欧美另类| 欧美日本中文国产一区发布| 亚州av有码| 天堂俺去俺来也www色官网| 老司机影院毛片| 大码成人一级视频| 少妇熟女欧美另类| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩精品成人综合77777| 久久韩国三级中文字幕| 久久狼人影院| 国产精品成人在线| 亚洲精品av麻豆狂野| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲av福利一区| 妹子高潮喷水视频| 男女免费视频国产| 91在线精品国自产拍蜜月| 黑人巨大精品欧美一区二区蜜桃 | 波野结衣二区三区在线| 国产高清国产精品国产三级| 九九爱精品视频在线观看| 一级,二级,三级黄色视频| 亚洲,一卡二卡三卡| 久久热精品热| 啦啦啦中文免费视频观看日本| 青春草视频在线免费观看| 26uuu在线亚洲综合色| 少妇熟女欧美另类| 欧美三级亚洲精品| 国产女主播在线喷水免费视频网站| 国产永久视频网站| 久久韩国三级中文字幕| 考比视频在线观看| 久久久久视频综合| 亚洲精品国产av蜜桃| 亚洲精品久久久久久婷婷小说| 我要看黄色一级片免费的| 国产在线免费精品| 亚洲av电影在线观看一区二区三区| 大又大粗又爽又黄少妇毛片口| 日韩精品有码人妻一区| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 美女内射精品一级片tv| 97超碰精品成人国产| 国产成人免费观看mmmm| 如日韩欧美国产精品一区二区三区 | 中文字幕最新亚洲高清| 国产精品久久久久久久电影| 午夜福利在线观看免费完整高清在| 狠狠婷婷综合久久久久久88av| 欧美亚洲 丝袜 人妻 在线| 精品国产乱码久久久久久小说| 日日撸夜夜添| 欧美性感艳星| 在线观看美女被高潮喷水网站| 观看美女的网站| 自拍欧美九色日韩亚洲蝌蚪91| 又大又黄又爽视频免费| 久久午夜综合久久蜜桃| 国产精品熟女久久久久浪| 99热6这里只有精品| 999精品在线视频| 极品少妇高潮喷水抽搐| 高清毛片免费看| 亚洲精品成人av观看孕妇| 能在线免费看毛片的网站| 天堂俺去俺来也www色官网| 秋霞伦理黄片| 在线看a的网站| 九九久久精品国产亚洲av麻豆| 日韩强制内射视频| 日韩免费高清中文字幕av| 久久精品国产亚洲网站| 国产亚洲一区二区精品| 国产免费又黄又爽又色| 在线精品无人区一区二区三| 中文字幕人妻丝袜制服| 美女主播在线视频| av一本久久久久| 下体分泌物呈黄色| 久久久a久久爽久久v久久| 人妻一区二区av| 观看av在线不卡| 高清av免费在线| 能在线免费看毛片的网站| 欧美丝袜亚洲另类| 中文字幕久久专区| 午夜福利,免费看| 亚洲无线观看免费| av网站免费在线观看视频| 男女边吃奶边做爰视频| 人妻人人澡人人爽人人| 老女人水多毛片| 久久久久久久久久人人人人人人| 男人操女人黄网站| 久久综合国产亚洲精品| 另类精品久久| 80岁老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 精品久久久久久久久亚洲| 免费观看无遮挡的男女| 中文字幕精品免费在线观看视频 | 熟妇人妻不卡中文字幕| 麻豆乱淫一区二区| 超碰97精品在线观看| 男女啪啪激烈高潮av片| 成年美女黄网站色视频大全免费 | av一本久久久久| 黑人猛操日本美女一级片| 国产日韩欧美亚洲二区| 国产精品99久久久久久久久| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区国产| 大又大粗又爽又黄少妇毛片口| 寂寞人妻少妇视频99o| 亚洲精品国产色婷婷电影| 精品久久久噜噜| 中文欧美无线码| 我要看黄色一级片免费的| 欧美xxⅹ黑人| 69精品国产乱码久久久| 日本av手机在线免费观看| 99精国产麻豆久久婷婷| 日本av免费视频播放| 夫妻午夜视频| 国产在视频线精品| 欧美日韩亚洲高清精品| 日韩一区二区视频免费看| 欧美最新免费一区二区三区| 蜜臀久久99精品久久宅男| 亚洲精品乱久久久久久| 亚洲色图综合在线观看| 极品少妇高潮喷水抽搐| 全区人妻精品视频| 狂野欧美激情性xxxx在线观看| 欧美丝袜亚洲另类| 久久鲁丝午夜福利片| av在线观看视频网站免费| 久久综合国产亚洲精品| 久久久精品区二区三区| 少妇精品久久久久久久| 亚洲欧美成人精品一区二区| 欧美精品一区二区大全| 亚洲av福利一区| videossex国产| a 毛片基地| 美女中出高潮动态图| 久久久久久久久久久免费av| 精品国产一区二区久久| 五月天丁香电影| 亚洲人成网站在线观看播放| 寂寞人妻少妇视频99o| 制服诱惑二区| 亚洲欧洲精品一区二区精品久久久 | 三级国产精品片| 美女视频免费永久观看网站| 欧美变态另类bdsm刘玥| 久久99蜜桃精品久久| 18禁在线播放成人免费| 黄色毛片三级朝国网站| 亚洲情色 制服丝袜| 韩国高清视频一区二区三区| 欧美人与善性xxx| 99国产精品免费福利视频| 国产在线免费精品| 亚洲国产精品专区欧美| 春色校园在线视频观看| 国产精品国产av在线观看| 两个人的视频大全免费| 人人澡人人妻人| 性色av一级| 伊人久久精品亚洲午夜| 久久久久人妻精品一区果冻| 插阴视频在线观看视频| 午夜福利视频在线观看免费| 中文天堂在线官网| 大码成人一级视频| 亚洲美女搞黄在线观看| 成人亚洲精品一区在线观看| 亚洲精品色激情综合| a级毛片黄视频| 新久久久久国产一级毛片| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| 婷婷色综合www| 欧美激情 高清一区二区三区| 日本午夜av视频| 熟女人妻精品中文字幕| 精品人妻熟女毛片av久久网站| 国产一区亚洲一区在线观看| 人妻系列 视频| 亚洲av二区三区四区| 久久久国产一区二区| 国产一级毛片在线| 成年美女黄网站色视频大全免费 | 午夜av观看不卡| 91精品伊人久久大香线蕉| 夜夜爽夜夜爽视频| 街头女战士在线观看网站| 免费观看的影片在线观看| 日韩精品免费视频一区二区三区 | 大话2 男鬼变身卡| 在线看a的网站| 狂野欧美激情性xxxx在线观看| 新久久久久国产一级毛片| 欧美bdsm另类| 51国产日韩欧美| 老司机影院毛片| 国产国拍精品亚洲av在线观看| 晚上一个人看的免费电影| 人成视频在线观看免费观看| 国产高清国产精品国产三级|