• <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ì)
    少妇猛男粗大的猛烈进出视频 | 黄色一级大片看看| 日韩精品青青久久久久久| 国产蜜桃级精品一区二区三区| 免费大片18禁| 在线看三级毛片| 国产又黄又爽又无遮挡在线| 亚洲色图av天堂| av天堂在线播放| 黑人高潮一二区| 国模一区二区三区四区视频| 亚洲国产精品久久男人天堂| 日本成人三级电影网站| 久久精品国产清高在天天线| 99riav亚洲国产免费| 国产精品久久久久久久久免| 内射极品少妇av片p| 日本免费一区二区三区高清不卡| 国产成年人精品一区二区| 一区二区三区免费毛片| 天堂影院成人在线观看| 日本黄色视频三级网站网址| 国产一区亚洲一区在线观看| 亚洲欧美日韩高清专用| 国产精品电影一区二区三区| 欧美+亚洲+日韩+国产| 久久精品影院6| 日本欧美国产在线视频| 午夜精品一区二区三区免费看| 熟女电影av网| 免费高清视频大片| 久99久视频精品免费| 日本一本二区三区精品| 国产免费一级a男人的天堂| 我的女老师完整版在线观看| 久久午夜亚洲精品久久| 国产精品野战在线观看| 日日啪夜夜撸| 午夜亚洲福利在线播放| 人妻久久中文字幕网| 特级一级黄色大片| 国产一区二区在线观看日韩| 久久天躁狠狠躁夜夜2o2o| 男女啪啪激烈高潮av片| 男女边吃奶边做爰视频| 亚洲成人久久爱视频| 久久精品国产亚洲网站| 国产精品久久视频播放| 别揉我奶头~嗯~啊~动态视频| a级毛片免费高清观看在线播放| 美女xxoo啪啪120秒动态图| 一区二区三区四区激情视频 | 成人av在线播放网站| 麻豆精品久久久久久蜜桃| 欧美另类亚洲清纯唯美| 国产在线精品亚洲第一网站| 亚洲丝袜综合中文字幕| 熟妇人妻久久中文字幕3abv| 国产白丝娇喘喷水9色精品| 亚洲欧美精品自产自拍| 午夜影院日韩av| 高清日韩中文字幕在线| 少妇猛男粗大的猛烈进出视频 | 成人综合一区亚洲| 日韩,欧美,国产一区二区三区 | 欧美色欧美亚洲另类二区| 人妻久久中文字幕网| 久久久久国产精品人妻aⅴ院| 91av网一区二区| 一夜夜www| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| av卡一久久| 国产久久久一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 国产伦精品一区二区三区视频9| 一个人观看的视频www高清免费观看| 中文字幕免费在线视频6| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看 | 晚上一个人看的免费电影| 国产伦一二天堂av在线观看| 亚洲自拍偷在线| 在线观看美女被高潮喷水网站| 亚洲成a人片在线一区二区| 日韩成人伦理影院| 看非洲黑人一级黄片| 丝袜美腿在线中文| av在线老鸭窝| 亚洲成人久久性| 国产黄a三级三级三级人| 麻豆成人午夜福利视频| 国产精品亚洲美女久久久| 成人三级黄色视频| 在线观看午夜福利视频| 中文字幕久久专区| 色在线成人网| 超碰av人人做人人爽久久| 99久国产av精品| 日韩精品中文字幕看吧| 亚州av有码| 欧美色视频一区免费| a级一级毛片免费在线观看| 欧美+亚洲+日韩+国产| 成人漫画全彩无遮挡| 成人鲁丝片一二三区免费| 欧美激情国产日韩精品一区| 日韩制服骚丝袜av| 国产淫片久久久久久久久| 美女被艹到高潮喷水动态| 91午夜精品亚洲一区二区三区| 亚洲性夜色夜夜综合| 成人鲁丝片一二三区免费| 日韩国内少妇激情av| 国内揄拍国产精品人妻在线| 日本欧美国产在线视频| 熟妇人妻久久中文字幕3abv| 亚洲综合色惰| 欧美色欧美亚洲另类二区| 成人高潮视频无遮挡免费网站| 欧美bdsm另类| 亚洲欧美日韩高清专用| 国产精品乱码一区二三区的特点| 国产精品一区二区三区四区久久| 夜夜夜夜夜久久久久| 国产成人freesex在线 | 赤兔流量卡办理| 亚洲自拍偷在线| 一个人观看的视频www高清免费观看| 一本精品99久久精品77| 久久欧美精品欧美久久欧美| 成人午夜高清在线视频| 日韩大尺度精品在线看网址| av在线观看视频网站免费| 校园人妻丝袜中文字幕| 久久久久精品国产欧美久久久| 成年女人看的毛片在线观看| 免费大片18禁| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 一级毛片我不卡| 中文字幕久久专区| 色噜噜av男人的天堂激情| 春色校园在线视频观看| 色5月婷婷丁香| 久久精品夜夜夜夜夜久久蜜豆| 久久天躁狠狠躁夜夜2o2o| 联通29元200g的流量卡| 长腿黑丝高跟| 99久久无色码亚洲精品果冻| 亚洲av成人精品一区久久| 日韩一本色道免费dvd| 日本免费a在线| 日本在线视频免费播放| 日产精品乱码卡一卡2卡三| 岛国在线免费视频观看| 成人欧美大片| 97人妻精品一区二区三区麻豆| 亚洲五月天丁香| 成人综合一区亚洲| 最新中文字幕久久久久| 最近的中文字幕免费完整| 内射极品少妇av片p| 免费av毛片视频| 国产精品美女特级片免费视频播放器| 一个人看的www免费观看视频| 日本精品一区二区三区蜜桃| 晚上一个人看的免费电影| 亚洲电影在线观看av| 久久久精品欧美日韩精品| 我要搜黄色片| 国产精品一区二区免费欧美| 亚洲国产精品sss在线观看| 成人无遮挡网站| 国产色婷婷99| 深夜a级毛片| 乱人视频在线观看| 久久精品夜色国产| 国产成人freesex在线 | 久久久国产成人精品二区| 一个人看的www免费观看视频| 少妇的逼好多水| 性欧美人与动物交配| 三级男女做爰猛烈吃奶摸视频| 99在线人妻在线中文字幕| 一级av片app| 在线a可以看的网站| 国产一区二区在线av高清观看| 男人狂女人下面高潮的视频| 午夜福利在线观看免费完整高清在 | 超碰av人人做人人爽久久| 一本精品99久久精品77| 精品久久久久久久久亚洲| 久久国内精品自在自线图片| 久久亚洲国产成人精品v| 国产精品久久久久久精品电影| 精品人妻熟女av久视频| 久久鲁丝午夜福利片| 日本精品一区二区三区蜜桃| 亚洲不卡免费看| 欧美一区二区精品小视频在线| 亚洲av第一区精品v没综合| 亚洲在线自拍视频| 国产精品av视频在线免费观看| 亚洲国产色片| 国产毛片a区久久久久| 亚洲国产欧洲综合997久久,| 欧美另类亚洲清纯唯美| 久久久久久九九精品二区国产| 一区二区三区四区激情视频 | 日本免费一区二区三区高清不卡| 18禁裸乳无遮挡免费网站照片| 在线播放无遮挡| 91麻豆精品激情在线观看国产| a级毛色黄片| 高清午夜精品一区二区三区 | 干丝袜人妻中文字幕| 国产视频一区二区在线看| 国产精品爽爽va在线观看网站| 国产av一区在线观看免费| 蜜臀久久99精品久久宅男| 欧美xxxx黑人xx丫x性爽| 久久久久久久久中文| 免费av毛片视频| 国产精品一二三区在线看| 精品不卡国产一区二区三区| av在线播放精品| 免费人成在线观看视频色| 欧美日韩精品成人综合77777| av女优亚洲男人天堂| 成年版毛片免费区| 中国国产av一级| 亚洲天堂国产精品一区在线| 国产午夜精品久久久久久一区二区三区 | 久久精品国产亚洲av天美| 亚洲激情五月婷婷啪啪| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| 亚洲最大成人手机在线| 在线a可以看的网站| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 搡女人真爽免费视频火全软件 | 免费观看在线日韩| 亚洲av第一区精品v没综合| 中文字幕免费在线视频6| 在线观看免费视频日本深夜| 一a级毛片在线观看| 亚洲欧美中文字幕日韩二区| 老熟妇乱子伦视频在线观看| 中国国产av一级| 天堂网av新在线| 黄片wwwwww| 国产成人aa在线观看| 亚洲专区国产一区二区| а√天堂www在线а√下载| 免费一级毛片在线播放高清视频| 国产精品国产三级国产av玫瑰| 中出人妻视频一区二区| ponron亚洲| 一级毛片我不卡| 日本熟妇午夜| 91麻豆精品激情在线观看国产| 好男人在线观看高清免费视频| 国产探花在线观看一区二区| 日韩欧美精品v在线| 亚洲最大成人手机在线| 亚洲人成网站在线播放欧美日韩| 91在线精品国自产拍蜜月| 联通29元200g的流量卡| 亚洲自偷自拍三级| 精品人妻一区二区三区麻豆 | 国产真实乱freesex| 在线a可以看的网站| 搡老妇女老女人老熟妇| 国产精品不卡视频一区二区| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕| 日本 av在线| 国产日本99.免费观看| 日本-黄色视频高清免费观看| 岛国在线免费视频观看| 不卡视频在线观看欧美| 最近视频中文字幕2019在线8| 亚洲人成网站高清观看| 精品乱码久久久久久99久播| 国内久久婷婷六月综合欲色啪| 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站在线播放欧美日韩| videossex国产| 日日摸夜夜添夜夜添av毛片| 久久久精品94久久精品| 真人做人爱边吃奶动态| 国产男人的电影天堂91| 观看免费一级毛片| 丰满人妻一区二区三区视频av| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 狂野欧美白嫩少妇大欣赏| 精品一区二区三区视频在线观看免费| 在线观看午夜福利视频| 精品不卡国产一区二区三区| 97热精品久久久久久| 久久久久国产网址| 亚洲欧美日韩卡通动漫| 国产毛片a区久久久久| 在线观看美女被高潮喷水网站| 麻豆一二三区av精品| 日本精品一区二区三区蜜桃| 精品久久久噜噜| 亚洲av.av天堂| 精品日产1卡2卡| 国产男靠女视频免费网站| 亚洲欧美中文字幕日韩二区| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| h日本视频在线播放| 一个人免费在线观看电影| 91在线观看av| 我要看日韩黄色一级片| 欧美高清成人免费视频www| 免费高清视频大片| 联通29元200g的流量卡| 精品久久国产蜜桃| 亚洲av中文字字幕乱码综合| 久久久久久大精品| 久久久欧美国产精品| 女生性感内裤真人,穿戴方法视频| 久久久精品大字幕| 伊人久久精品亚洲午夜| 少妇裸体淫交视频免费看高清| 赤兔流量卡办理| 国产av不卡久久| 真实男女啪啪啪动态图| 国产精品免费一区二区三区在线| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 蜜桃亚洲精品一区二区三区| av中文乱码字幕在线| 狂野欧美白嫩少妇大欣赏| av.在线天堂| 国产精品女同一区二区软件| 午夜久久久久精精品| 国产高清视频在线播放一区| 国产精品久久视频播放| 国产精品亚洲一级av第二区| 高清毛片免费观看视频网站| 91狼人影院| 99热只有精品国产| 少妇熟女欧美另类| 精品人妻偷拍中文字幕| 可以在线观看的亚洲视频| 国产精品一区二区性色av| 亚洲四区av| 成年女人看的毛片在线观看| 亚洲国产精品成人久久小说 | 久久精品影院6| 久久久精品大字幕| 精品久久久久久久久av| 午夜视频国产福利| 激情 狠狠 欧美| 国产精品野战在线观看| 一进一出好大好爽视频| 国产视频内射| 日本黄色视频三级网站网址| 嫩草影院入口| 国产伦在线观看视频一区| av免费在线看不卡| 亚洲人成网站在线播| 黄片wwwwww| 午夜精品国产一区二区电影 | 久久中文看片网| 久久精品国产鲁丝片午夜精品| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲av嫩草精品影院| 九九在线视频观看精品| 天美传媒精品一区二区| 日韩大尺度精品在线看网址| 3wmmmm亚洲av在线观看| 中文字幕免费在线视频6| 午夜福利18| 亚洲人与动物交配视频| 亚洲专区国产一区二区| 一夜夜www| 久久鲁丝午夜福利片| 日韩av在线大香蕉| 少妇高潮的动态图| 99riav亚洲国产免费| av.在线天堂| 久久久色成人| 久久欧美精品欧美久久欧美| 久久草成人影院| 男人舔女人下体高潮全视频| 级片在线观看| av免费在线看不卡| 啦啦啦韩国在线观看视频| 亚洲国产欧美人成| 91狼人影院| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 男女视频在线观看网站免费| 99视频精品全部免费 在线| 99riav亚洲国产免费| 日本-黄色视频高清免费观看| 1024手机看黄色片| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 日韩中字成人| 尤物成人国产欧美一区二区三区| 国产精品人妻久久久影院| 国内精品久久久久精免费| 少妇的逼好多水| 看黄色毛片网站| 免费大片18禁| 人人妻人人澡人人爽人人夜夜 | 大香蕉久久网| 午夜免费激情av| 成年女人看的毛片在线观看| 最近2019中文字幕mv第一页| 亚洲国产精品sss在线观看| 色av中文字幕| 国产私拍福利视频在线观看| 人妻丰满熟妇av一区二区三区| 久久久久久国产a免费观看| 午夜日韩欧美国产| 亚洲av第一区精品v没综合| 婷婷精品国产亚洲av| 春色校园在线视频观看| 男女啪啪激烈高潮av片| 老司机午夜福利在线观看视频| 精品久久久久久久久亚洲| 久久精品国产鲁丝片午夜精品| 黑人高潮一二区| 欧美潮喷喷水| 欧美性感艳星| 人人妻人人澡人人爽人人夜夜 | 少妇高潮的动态图| 午夜福利视频1000在线观看| 日韩欧美 国产精品| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 蜜桃亚洲精品一区二区三区| 亚洲国产欧洲综合997久久,| 99久久精品一区二区三区| 久久久久久久午夜电影| 欧美人与善性xxx| 欧美激情在线99| 成人亚洲欧美一区二区av| 精品欧美国产一区二区三| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 97碰自拍视频| 午夜精品国产一区二区电影 | 一进一出好大好爽视频| 国产探花在线观看一区二区| 乱人视频在线观看| 亚洲av二区三区四区| 国产亚洲精品久久久com| 一进一出抽搐gif免费好疼| 免费在线观看成人毛片| 一a级毛片在线观看| 丝袜喷水一区| 淫秽高清视频在线观看| 99久久九九国产精品国产免费| 午夜免费激情av| 中文字幕av成人在线电影| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 淫妇啪啪啪对白视频| 97在线视频观看| 亚洲最大成人中文| 久久久久九九精品影院| 久久久午夜欧美精品| 极品教师在线视频| 亚洲精品色激情综合| 亚洲成av人片在线播放无| 日韩av在线大香蕉| 麻豆久久精品国产亚洲av| 桃色一区二区三区在线观看| 黄片wwwwww| 精品一区二区三区av网在线观看| 欧美一区二区亚洲| 九九在线视频观看精品| 秋霞在线观看毛片| 一区福利在线观看| 国产精品电影一区二区三区| 国产高清有码在线观看视频| 欧美三级亚洲精品| 麻豆av噜噜一区二区三区| 日韩欧美一区二区三区在线观看| 一级毛片aaaaaa免费看小| 在线观看美女被高潮喷水网站| 免费av不卡在线播放| 国产欧美日韩一区二区精品| 午夜影院日韩av| 精品久久久久久久久久久久久| 97在线视频观看| 中文字幕av在线有码专区| av在线天堂中文字幕| 亚洲专区国产一区二区| 午夜激情福利司机影院| 99久久精品热视频| 小蜜桃在线观看免费完整版高清| 99久久无色码亚洲精品果冻| 日本免费一区二区三区高清不卡| 午夜精品国产一区二区电影 | 国产精品爽爽va在线观看网站| 亚洲av五月六月丁香网| 午夜爱爱视频在线播放| 在线观看一区二区三区| 美女被艹到高潮喷水动态| 69人妻影院| 久久6这里有精品| 免费观看在线日韩| 性欧美人与动物交配| 精品乱码久久久久久99久播| 国产伦在线观看视频一区| 亚洲va在线va天堂va国产| 亚洲精品在线观看二区| 直男gayav资源| 久久精品综合一区二区三区| 亚洲色图av天堂| 熟女人妻精品中文字幕| 国产不卡一卡二| 欧美潮喷喷水| 日本免费a在线| 成年女人看的毛片在线观看| 免费看av在线观看网站| 五月玫瑰六月丁香| 国产精品乱码一区二三区的特点| 国产片特级美女逼逼视频| 日韩欧美精品v在线| 99热这里只有是精品50| 我要搜黄色片| 色哟哟哟哟哟哟| 精品日产1卡2卡| 九色成人免费人妻av| 大香蕉久久网| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 久久久久久久久中文| 97热精品久久久久久| 蜜臀久久99精品久久宅男| 国产aⅴ精品一区二区三区波| 国产精品爽爽va在线观看网站| 成人亚洲欧美一区二区av| 在线看三级毛片| 美女黄网站色视频| 亚洲国产日韩欧美精品在线观看| 日韩欧美精品免费久久| 日本五十路高清| 十八禁国产超污无遮挡网站| 一级毛片久久久久久久久女| 色哟哟哟哟哟哟| 精品日产1卡2卡| 男女边吃奶边做爰视频| 亚洲精品一卡2卡三卡4卡5卡| 一级黄片播放器| 我的女老师完整版在线观看| 国产精品99久久久久久久久| 国内精品宾馆在线| www日本黄色视频网| 国产伦精品一区二区三区四那| 欧美在线一区亚洲| 男人舔奶头视频| 国产大屁股一区二区在线视频| 乱人视频在线观看| av福利片在线观看| 麻豆乱淫一区二区| 深夜精品福利| 婷婷亚洲欧美| 韩国av在线不卡| 美女黄网站色视频| 看免费成人av毛片| 熟女人妻精品中文字幕| 丰满的人妻完整版| 欧美人与善性xxx| 国产真实乱freesex| 日韩制服骚丝袜av| 国产精品亚洲美女久久久| 亚洲欧美日韩高清专用| 搡女人真爽免费视频火全软件 | 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点| 可以在线观看的亚洲视频| 国产精品久久久久久精品电影| 欧美人与善性xxx| 成人av在线播放网站| 丰满的人妻完整版| 特级一级黄色大片| 国产精品美女特级片免费视频播放器| 免费av毛片视频| 一个人看视频在线观看www免费| 精品欧美国产一区二区三| 在线观看一区二区三区| 美女被艹到高潮喷水动态| 午夜老司机福利剧场| 少妇人妻一区二区三区视频| 99热精品在线国产| 久久久久国内视频| 欧美日韩在线观看h| 99久久精品一区二区三区| 俺也久久电影网| 国产一区二区在线观看日韩| 在线国产一区二区在线| 亚洲av免费高清在线观看| 久久亚洲精品不卡| 亚洲成av人片在线播放无| 精品欧美国产一区二区三| 久久99热这里只有精品18| 一个人看视频在线观看www免费|