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

    太湖不同湖區(qū)冬季沉積物細(xì)菌群落多樣性

    2018-03-02 05:24:18薛銀剛江曉棟耿金菊滕加泉謝文理陳心一常州大學(xué)環(huán)境與安全工程學(xué)院江蘇常州64江蘇省環(huán)境保護(hù)水環(huán)境生物監(jiān)測重點實驗室常州市環(huán)境監(jiān)測中心江蘇常州00污染控制與資源化研究國家重點實驗室南京大學(xué)環(huán)境學(xué)院江蘇南京00
    中國環(huán)境科學(xué) 2018年2期
    關(guān)鍵詞:河口區(qū)湖心湖區(qū)

    薛銀剛,劉 菲,江曉棟,耿金菊,滕加泉,謝文理,張 皓,陳心一 (.常州大學(xué)環(huán)境與安全工程學(xué)院,江蘇 常州 64;.江蘇省環(huán)境保護(hù)水環(huán)境生物監(jiān)測重點實驗室,常州市環(huán)境監(jiān)測中心,江蘇 常州 00;.污染控制與資源化研究國家重點實驗室,南京大學(xué)環(huán)境學(xué)院,江蘇 南京 00)

    太湖是我國第三大淺層淡水湖,富營養(yǎng)化是其目前存在的主要環(huán)境問題[1].湖泊沉積層對富營養(yǎng)化非常敏感[2],沉積物細(xì)菌群落在氮、磷等元素的礦化和循環(huán)中發(fā)揮重要作用,影響著淡水湖泊營養(yǎng)物質(zhì)的生物化學(xué)循環(huán)[3-4].營養(yǎng)物質(zhì)的可利用性使得沉積物細(xì)菌群落在不同營養(yǎng)水平湖區(qū)之間存在差異.

    以往的研究利用T-RFLP[5-6]和PCRDGGE[7]等方法對太湖浮游細(xì)菌群落特征已經(jīng)展開了大量的工作.隨著分子生物學(xué)測序的發(fā)展,高通量測序技術(shù)被證實是識別微生物群落整體概況的高效工具[8],并被廣泛應(yīng)用于土壤、生物膜和沉積物等復(fù)雜環(huán)境的微生物群落結(jié)構(gòu)研究[9-12].目前,關(guān)于太湖微生物群落的研究大多集中于北部藻型湖區(qū)[13-14]和浮游細(xì)菌多樣性[5-7],本文基于高通量測序?qū)Σ煌瑺I養(yǎng)水平湖區(qū)冬季沉積物細(xì)菌群落多樣性展開研究.

    太湖自2007年藍(lán)藻大面積暴發(fā)后,相關(guān)部門采取了一系列治理措施,如控源截污、生態(tài)修復(fù)、引水稀釋和藍(lán)藻打撈等.對于藍(lán)藻水華的防控主要集中在春夏兩季對藍(lán)藻的大規(guī)模打撈,對于早期的藍(lán)藻防控工作相對缺乏.本文以太湖不同營養(yǎng)水平湖區(qū)表層沉積物為研究對象,利用高通量測序技術(shù)分析各個點位與不同湖區(qū)的沉積物細(xì)菌多樣性,解析不同湖區(qū)沉積物細(xì)菌群落的空間差異;通過環(huán)境因子相關(guān)性分析,研究水環(huán)境因子對沉積物細(xì)菌群落的影響.通過對冬季沉積物細(xì)菌群落結(jié)構(gòu)進(jìn)行分析,以期為太湖沉積物細(xì)菌生物化學(xué)循環(huán)和實施早期藍(lán)藻防控工作提供基礎(chǔ)數(shù)據(jù).

    1 材料與方法

    1.1 樣品采集與理化分析

    采樣時間為2015年11月,采集的樣品為表層沉積物,采樣點為太湖不同營養(yǎng)水平湖區(qū)(圖1).采樣點依次是沙塘港(STG)、閭江口(LJK)、小灣里(XWL)、椒山(JS)、平臺山(PTS)、西山西(XSX)、漁洋山(YYS)、胥湖心(XHX)、七都(QD).其中STG、LJK和XWL位于藻類主導(dǎo)的河口區(qū)(HK);JS、PTS和XSX位于受外部因素影響最小的湖心區(qū)(HX);YYS、XHX和QD位于水生植物主導(dǎo)的水草區(qū)(SC).使用便攜式柱狀采泥器(HYDRO-BIOS, Germany)進(jìn)行沉積物樣品的采集.在運(yùn)輸過程中,所有的樣品均放置在冷藏箱中保存,并在送到實驗室后立即進(jìn)行預(yù)處理.將采集的沉積物樣品離心后得到的固態(tài)污泥保存于-20℃條件下,用于后續(xù)分析測試.

    利用YSI水質(zhì)分析儀(美國,YSI 6600)現(xiàn)場測定上覆水的水溫和pH值.每個采樣點采集3份水樣帶回實驗室測定TN、NH3-N、、、TP和水質(zhì)指標(biāo),測定方法參照《水和廢水監(jiān)測分析方法》(第四版)[15].

    圖1 太湖采樣點位置Fig.1 Location of sampling sites in Lake Taihu

    1.2 DNA提取和PCR擴(kuò)增

    使用FastDNA Spin Kit for Soil試劑盒(MP bio, USA)按照試劑盒說明步驟提取基因組DNA,使用NanoDrop2000超微量蛋白質(zhì)核酸分析儀(Thermo Fisher Scientific, USA)測定提取出的DNA濃度和純度后,將DNA樣品于-20℃中保存,用于后續(xù)的PCR擴(kuò)增等分析.

    將提取好的DNA產(chǎn)物利用16S rDNA V1~V2區(qū)擴(kuò)增引物進(jìn)行PCR擴(kuò)增,引物序列(5'-3'):正向引物AGAGTTTGATYMTGGCTCAG,反向引物TGCTGCCTCCCGTAGGAGT.PCR產(chǎn)物以等密度比混合后使用Qiagen Gel Extraction Kit(Qiagen, Germany)純化.

    1.3 高通量測序

    純化后的產(chǎn)物送至江蘇中宜金大分析檢測有限公司,使用Illumina(USA)的Hiseq平臺進(jìn)行高通量測序和分析.經(jīng)TruSeq?DNA PCR-Free樣品制備試劑盒(Illumina, USA)生成測序文庫.使用Qubit @ 2.0熒光計(Thermo Fisher Scientific,USA)和Agilent Bioanalyzer 2100系統(tǒng)(Agilent,USA)評估測序文庫質(zhì)量.最后,在Illumina HiSeq 2500平臺(Illumina, USA)上進(jìn)行測序,產(chǎn)生了250bp的配對末端讀數(shù).

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

    測序得到的原始數(shù)據(jù),存在一定比例的干擾數(shù)據(jù),為了使信息分析的結(jié)果更加準(zhǔn)確、可靠,首先對原始數(shù)據(jù)進(jìn)行拼接、過濾,得到有效數(shù)據(jù).然后基于有效數(shù)據(jù)進(jìn)行OTUs聚類和物種分類分析.根據(jù)OTUs聚類結(jié)果,對每個OTU的代表序列做物種注釋,得到對應(yīng)的物種信息和基于物種的豐度分布情況.對于每個代表性的序列,使用基于RDP分類器算法的GreenGene數(shù)據(jù)庫來注釋分類信息[16].

    為了研究不同OTU的系統(tǒng)發(fā)育關(guān)系,以及不同樣本中優(yōu)勢種的差異,使用MUSCLE軟件進(jìn)行多序列比對(3.8.31版本);基于QIIME計算樣品香農(nóng)指數(shù)(Shannon)、辛普森指數(shù)(Simpson)、文庫覆蓋率(Coverage)和系統(tǒng)發(fā)育距離,使用平均聯(lián)動的非加權(quán)組平均法的聚類分析方法(UPGMA)來解釋距離矩陣構(gòu)建系統(tǒng)發(fā)育樹[17];R語言程序(2.15.3版)繪制韋恩圖;Origin(8.5版本)繪制屬水平優(yōu)勢物種豐度柱狀圖;優(yōu)勢物種差異顯著性利用SPSS(17.0版本)軟件進(jìn)行分析,顯著性水平設(shè)定為P<0.05;結(jié)合Pearson相關(guān)性分析和CANOCO軟件(4.5版本)來進(jìn)行優(yōu)勢細(xì)菌與水環(huán)境因子的冗余分析(RDA)[18].

    2 結(jié)果與分析

    2.1 水樣理化指標(biāo)

    太湖冬季不同湖區(qū)上覆水的理化指標(biāo)如表1,可以看出各個湖區(qū)的營養(yǎng)水平為河口區(qū)>湖心區(qū)>水草區(qū).太湖水體整體呈弱堿性,水溫在15℃左右.TN含量在不同湖區(qū)之間差異明顯,TP含量在不同湖區(qū)之間差異性較弱.水草區(qū)的NH3-N含量較高;河口區(qū)和湖心區(qū)的TN、TP、、和含量較高.

    表1 太湖不同湖區(qū)冬季水質(zhì)理化指標(biāo)Table 1 Physio-chemical indices of water in different zones of Lake Taihu in winter

    2.2 OTU聚類與多樣性分析

    各采樣點的OTU聚類結(jié)果與多樣性指數(shù)見表2,測序文庫的覆蓋率均在99%以上,表明各個采樣點的測序結(jié)果完全可以代表樣本的真實情況.每個采樣點的有效序列數(shù)(用于OTU聚類等后續(xù)分析的有效數(shù)據(jù))均在20000以上,本次測序產(chǎn)生的OTUs是將所有樣品的有效序列進(jìn)行統(tǒng)計,以97%的一致性將序列聚類而成[19],然后對OTUs的代表序列進(jìn)行物種注釋.稀釋曲線(圖2)表明STG、LJK和XWL在有效測序數(shù)達(dá)到4500時,曲線已保持平緩狀態(tài);其余采樣點當(dāng)測序有效序列達(dá)到8000時,曲線漸進(jìn)平緩;表明用于OTU聚類的測序數(shù)滿足要求[20],且聚類產(chǎn)生的OTUs可以用于后續(xù)物種注釋.稀釋曲線反映出各個采樣點聚類產(chǎn)生的OTUs為YYS>XHX>QD>JS>XSX>PTS>XWL>STG>LJK,這與表2所述基本一致.

    基于多序列比對計算得出各采樣點的Shannon指數(shù)和Simpson指數(shù),結(jié)果表明,各采樣點的Shannon指數(shù)為YYS>QD>XHX>JS>XWL>XSX>PTS>STG>LJK;Simpson指數(shù)為YYS>QD>XWL>XHX>XSX>JS>PTS>STG>LJK.總體來看,Shannon指數(shù)和Simpson指數(shù)為水草區(qū)>湖心區(qū)>河口區(qū).

    表2 各點位OTU聚類與多樣性分析Table 2 Clustering and diversity analysis of each sampling site

    圖2 稀釋曲線Fig.2 Rarefaction curve

    Venn圖用于展示可視化數(shù)據(jù)集之間的重疊和差異,圖3反映了河口區(qū)、湖心區(qū)和水草區(qū)樣品的OTUs組成相似性以及數(shù)據(jù)重疊情況.不同營養(yǎng)水平湖區(qū)均具有其特有的細(xì)菌種類,3個湖區(qū)OTUs共計1006,共有的OTUs為78.湖心區(qū)和水草區(qū)共有的OTU數(shù)最多,其次為河口區(qū)和湖心區(qū)共有的OTU數(shù),河口區(qū)與水草區(qū)共有的OTU數(shù)最少.

    圖3 OTU分布韋恩圖Fig.3 Venn diagram denoting the distribution of OTUs

    描述湖區(qū)到湖區(qū)物種組成的β多樣性對于了解各種空間尺度上物種多樣性的模式至關(guān)重要[21],可以考察不同湖區(qū)的生境多樣性,有利于了解整個太湖區(qū)域的生物生產(chǎn)力差異.選用距離指標(biāo)來衡量不同湖區(qū)樣品間的相異系數(shù)(圖4).水草區(qū)與湖心區(qū)相異系數(shù)較小(0.214),水草區(qū)與河口區(qū)以及河口區(qū)與湖心區(qū)相異系數(shù)較大,分別為0.650和0.613.

    圖4 Beta多樣性指數(shù)Fig.4 Beta diversity index

    2.3 沉積物細(xì)菌群落結(jié)構(gòu)

    高通量測序結(jié)果表明,河口區(qū)、湖心區(qū)和水草區(qū)冬季沉積物樣品聚類產(chǎn)生的OTUs分別為378、669和641,經(jīng)分析主要隸屬于25個門,51個綱,96個屬.在門水平上,優(yōu)勢類群主要包括綠彎菌門(Chloroflexi)、變形菌門(Proteobacteria)、酸桿菌門(Acidobacteria)、藍(lán)藻門(Cyanobacteria)和擬桿菌門(Bacteroidetes)等(圖5).系統(tǒng)發(fā)育樹表明位于同一營養(yǎng)水平湖區(qū)沉積物的細(xì)菌組成較為相似.就河口區(qū)而言,豐度高于10%的優(yōu)勢類群為Proteobacteria(29.53%)、硝化螺旋菌門(Nitrospirae,16.47%)和Chloroflexi(13.80%).除Chloroflexi(34.78%)和Proteobacteria(16.16%)外,Cyanobacteria在湖心區(qū)豐度高達(dá)19.81%.Chloroflexi同時也是水草區(qū)豐度最高的優(yōu)勢類群,占比約為總細(xì)菌含量的1/3,其次為Proteobacteria(17.64%)、Acidobacteria(12.62%)和Cyanobacteria(10.32%).另外,Bacteroidetes、綠菌門(Chlorobi)、放線菌門(Actinobacteria)、芽單胞菌門(Gemmatimonadetes)和浮霉菌門(Planctomycetes)等優(yōu)勢類群也均有檢出,豐度均高于1%.

    圖5可以看出各個點位優(yōu)勢類群的豐度存在明顯差別,Proteobacteria和Chloroflexi是太湖各個點位沉積物中相對豐度最高的兩個類群.9個采樣點占比第一的沉積物細(xì)菌分別為STG(Proteobacteria,37.48%)、LJK(Proteobacteria,32.99%)、XWL(Proteobacteria,18.12%)、JS(Proteobacteria,34.44%)、PTS(Chloroflexi,36.67%)、XSX(Chloroflexi,48.15%)、XHX(Chloroflexi,30.58%)、QD(Chloroflexi,37.80%)、YYS(Proteobacteria,28.99%).其中,Chloroflexi中豐度較高的是厭氧繩菌綱(Anaerolineae,11.29%)和KD4-96(7.48%);Proteobacteria中豐度較高的是β-變形菌綱(Betaproteobacteria,3.03%)、γ-變形菌綱(Gammaproteobacteria,9.48%)和δ-變形菌綱(Deltaproteobacteria,3.73%).

    圖5 UPGMA聚類樹與門水平優(yōu)勢沉積物細(xì)菌Fig.5 UPGMA cluster tree of dominant bacteria in winter sediment of Lake Taihu at the phyla level

    選取屬水平排名前10的優(yōu)勢細(xì)菌重新進(jìn)行相對豐度的計算,以更顯著地表明不同湖區(qū)優(yōu)勢沉積物細(xì)菌的豐度差異性(圖6).其中歸屬于Others的細(xì)菌主要有蛭弧菌屬(Bdellovibrio)、厭氧粘細(xì)菌屬(Anaeromyxobacter)、Blvii28_wastewater-sludge_group以及分枝桿菌屬(Mycobacterium)等(圖中未展示). unidentified_Chloroplast、聚球藻屬(Synechococcus)、微囊藻屬(Microcystis)和unidentified_ Nitrosomonadaceae是太湖冬季沉積物中優(yōu)勢類群.其中,河口區(qū)的主要優(yōu)勢類群為Microcystis、unidentified_Nitrosomonadaceae和unidentified_Chloroplast,平均豐度分別達(dá)到26.24%、24.14%和18.87%;這3類細(xì)菌也是湖心區(qū)豐度排名前3的優(yōu)勢類群;而水草區(qū)最具優(yōu)勢的類群為Synechococcus,平均豐度達(dá)到31.65%.值得注意的是,Microcystis在每個采樣點均有檢出,其中JS的Microcystis含量最高,僅次于JS的是LJK.unidentified_Chloroplast、Synechococcus、unidentified_Nitrospiraceae、Blastocatella和Sva0081_sediment_groups分別在XSX、QD、XWL、PTS、STG豐度最高.

    選取5類門水平優(yōu)勢細(xì)菌、3類綱水平優(yōu)勢細(xì)菌與4類屬水平優(yōu)勢細(xì)菌與不同湖區(qū)進(jìn)行差異顯著性分析(表3).單因素方差分析表明,不同優(yōu)勢細(xì)菌在湖區(qū)之間的差異顯著性不同.2類門水平優(yōu)勢細(xì)菌Nitrospirae(P=0.033)和Acidobacteria (P=0.042)、2類綱水平優(yōu)勢細(xì)菌Anaerolineae(P=0.043)和Nitrospira(P=0.033)在河口區(qū)、湖心區(qū)和水草區(qū)之間差異顯著.3類門水平優(yōu)勢細(xì)菌(Chloroflexi、Proteobacteria和Cyanobacteria)、1類綱水平優(yōu)勢細(xì)菌(Gammaproteobacteria)和4類屬水平優(yōu)勢細(xì)菌(unidentified_Chloroplast、Synechococcus、Microcystis和unidentified_Nitrospiraceae)在不同湖區(qū)之間差異不顯著(P>0.05).

    圖6 屬水平優(yōu)勢沉積物細(xì)菌Fig.6 Dominant bacteria in sediment at the genus level

    表3 優(yōu)勢細(xì)菌與不同湖區(qū)差異顯著性Table 3 Significance of difference of dominant bacteria among different lake zones

    2.4 沉積物細(xì)菌與水環(huán)境因子的相關(guān)性

    圖7展示了水環(huán)境營養(yǎng)鹽含量、水溫以及pH值與屬水平細(xì)菌的冗余分析結(jié)果,第一主軸對屬水平細(xì)菌群落方差變化的解釋量為62.40%,第二主軸的解釋量為18.80%.其中pH值與第一主軸呈顯著正相關(guān),與第二主軸呈顯著負(fù)相關(guān);、水溫、TN、TP、、和NH3-N與第一主軸呈顯著負(fù)相關(guān),與第二主軸呈顯著正相關(guān).可見,上層水環(huán)境營養(yǎng)鹽水平、水溫以及pH值均對太湖冬季表層沉積物細(xì)菌群落有一定影響.

    unidentified_Xanthomonadales、Ignavibacterium和Desulfatiglans與TN、TP、、和NH3-N具有較強(qiáng)正相關(guān)性,相關(guān)系數(shù)均大于0.9;與pH值呈負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)絕對值大于0.99.Synechococcus和Blastocatella與pH值呈正相關(guān)性關(guān)系,與TN、TP、、和NH3-N呈負(fù)相關(guān)關(guān)系. Sva0081_sediment_group、CL500-3和uidentified_Nitrospiraceae受環(huán)境因子影響較小,NH3-N是其主要影響因子.Microcystis顯示出與水溫和較強(qiáng)的正相關(guān)關(guān)系,相關(guān)系數(shù)均大于0.99.

    圖7 優(yōu)勢細(xì)菌與水環(huán)境因子的冗余分析Fig.7 Redundancy analysis of dominant bacteria and water environmental factors

    3 討論

    3.1 太湖不同湖區(qū)沉積物細(xì)菌群落特征

    太湖水域較廣闊,由于其水深相對較淺,外加風(fēng)浪的作用促使沉積層物質(zhì)與水體物質(zhì)交換頻繁[22],引起沉積層理化性質(zhì)與營養(yǎng)程度的變化,進(jìn)而影響微生物的多樣性.為了研究不同湖區(qū)沉積物細(xì)菌群落差異性和生境多樣性,本文選用α多樣性和β多樣性對太湖不同營養(yǎng)水平湖區(qū)進(jìn)行分析.結(jié)果表明,湖心區(qū)和水草區(qū)的細(xì)菌群落組成較為相似;通過Shannon指數(shù)和Simpson指數(shù)對每個樣品中物種組成的豐富度和均勻度兩方面進(jìn)行考察[23],結(jié)果表明,水草區(qū)物種組成的均勻度和豐富度均最高,其次為湖心區(qū),河口區(qū)最低;之前已有研究報道表明營養(yǎng)水平的增加會導(dǎo)致細(xì)菌群落多樣性顯著下降[24].造成太湖3個湖區(qū)群落多樣性差異可能是由于水草區(qū)主要由水生植物主導(dǎo),富營養(yǎng)化程度較低,水質(zhì)較好;湖心區(qū)地理位置上遠(yuǎn)離河口區(qū),受外界影響小,富營養(yǎng)化程度較河口區(qū)低;由于周邊農(nóng)業(yè)、水產(chǎn)養(yǎng)殖業(yè)、旅游業(yè)和城市化的發(fā)展,河口區(qū)長期接納農(nóng)業(yè)化學(xué)品、土壤顆粒、工業(yè)廢水和生活污水的侵蝕,營養(yǎng)輸入過多,水質(zhì)較差.營養(yǎng)通量的改變可能導(dǎo)致生態(tài)系統(tǒng)發(fā)生巨大變化,使得不同營養(yǎng)區(qū)的細(xì)菌群落組成各異.Bacteroidetes、Chlorobi、unidentified_Chloroplast和unidentified_Nitrospiraceae在富營養(yǎng)區(qū)的豐度更高,而Actinobacteria在貧營養(yǎng)區(qū)的檢出率更高,劉芳華等[25]也發(fā)現(xiàn)Actinobacteria在太湖貧營養(yǎng)區(qū)的含量相對較高,營養(yǎng)可用性被認(rèn)為是改變湖泊沉積物中微生物群落的主要因素[26].營養(yǎng)水平的不同使Nitrospirae、Acidobacteria、Anaerolineae和Nitrospira在3個湖區(qū)之間具有顯著的差異.同時,不同營養(yǎng)區(qū)還存在特有的細(xì)菌種類,如Sva0081_ sediment_group僅在水草區(qū)的3個采樣點有檢出.

    沉積物細(xì)菌具有種類豐富以及代謝活性高的特點,是湖泊中完成轉(zhuǎn)化營養(yǎng)物質(zhì)工作的承擔(dān)者[27].太湖冬季沉積物中主要的優(yōu)勢類群為Chloroflexi、Proteobacteria、Acidobacteria和Cyanobacteria等.高光等[28]通過對太湖水生植物主導(dǎo)和藻類主導(dǎo)湖區(qū)冬季和夏季的沉積物細(xì)菌群落進(jìn)行比較,發(fā)現(xiàn)Proteobacteria在兩類湖區(qū)豐度均最高,除此之外,Cyanobacteria僅在藻類主要湖區(qū)的樣品中有檢出,水生植物主導(dǎo)湖區(qū)的優(yōu)勢類群為Acidobacteria、Bacteroidetes.該研究團(tuán)隊對太湖不同湖區(qū)春季沉積物細(xì)菌群落的研究結(jié)果表明,梅梁灣和水草區(qū)的主要優(yōu)勢細(xì)菌為Chloroflexi、Proteobacteria;Proteobacteria和Bacteroidetes在湖心區(qū)豐度較高[25].可見,Proteobacteria和Chloroflexi是太湖沉積物中長期存在的優(yōu)勢類群.大多數(shù)的Proteobacteria對于可降解的有機(jī)基質(zhì)具有很強(qiáng)的代謝能力,在富營養(yǎng)化湖泊占據(jù)主導(dǎo)地位并不意外.Chloroflexi是涉及有機(jī)污染物生物降解的細(xì)菌,太湖的富營養(yǎng)狀態(tài)使這類細(xì)菌具有明顯的生存優(yōu)勢.趙大勇[29]等利用T-RFLP方法對太湖3個不同營養(yǎng)水平湖區(qū)(梅梁灣、胥口灣、東太湖)夏季表層沉積物的細(xì)菌群落開展研究,發(fā)現(xiàn)Gammaproteobacteria廣泛存在于3個湖區(qū)中,Deltaproteobacteria在梅梁灣豐度較高,Betaproteobacteria為東太湖含量較多,這與本次研究結(jié)果相類似;但其僅在梅梁灣區(qū)域檢出Chloroflexi,而本次研究表明Chloroflexi分布于太湖不同湖區(qū)所有采樣點,可能是由于高通量測序的檢測范圍寬、定量準(zhǔn)確以及測序深度高等優(yōu)點.

    近幾年來,盡管太湖的營養(yǎng)污染得到一定的削減,但藍(lán)藻的繁殖頻率和持續(xù)時間都有所增加[30].早前的報道表明,全球氣候的變化導(dǎo)致地表水溫升高,可能會在藍(lán)藻繁殖中起重要作用[31].針對富營養(yǎng)化和氣候變化促進(jìn)藍(lán)藻有害增殖和擴(kuò)張的兩個過程,O'Neil等[32]研究了淡水、河口和海洋生態(tài)系統(tǒng)的代表藍(lán)藻屬與營養(yǎng)和氣候變化之間的關(guān)系,結(jié)果表明藍(lán)藻屬對氮磷等營養(yǎng)元素具有極高的競爭力;并且在更高溫度下,其生長速度更快.藍(lán)藻水華暴發(fā)往往經(jīng)歷4個階段:湖底休眠、水體復(fù)蘇、上浮積聚和死亡分解.很多藍(lán)藻細(xì)胞在冬季會沉到底泥里,而到春季會復(fù)蘇到水體,且藍(lán)藻大規(guī)模積聚的時間提前,相關(guān)部門在春季便開始對太湖藍(lán)藻實施大規(guī)模打撈,但藍(lán)藻并未減少,且2017年藍(lán)藻打撈量是2016年同期的5.5倍[33].本次通過高通量測序技術(shù)在太湖9個采樣點均檢出藍(lán)藻,在3個湖區(qū)占比均較高,以湖心區(qū)含量最高,其次為水草區(qū),河口區(qū)豐度最低.研究表明,沉積物中的營養(yǎng)物質(zhì)會為藍(lán)藻復(fù)蘇及生長繁殖提供營養(yǎng)[34],致使春季太湖水體中藍(lán)藻大規(guī)模富集.可見,下一年春季太湖水體中藍(lán)藻會大規(guī)模積聚,可能以湖心區(qū)藍(lán)藻覆蓋面積更大.設(shè)想如若將冬季沉積物中的藻種消除,便可能克服后續(xù)大規(guī)模的藍(lán)藻打撈工作,從而實現(xiàn)早期藍(lán)藻水華的防控.

    3.2 沉積物細(xì)菌群落與環(huán)境因子的關(guān)系

    營養(yǎng)鹽的輸入可能對接收水質(zhì)產(chǎn)生深遠(yuǎn)的影響,同時上覆水的理化性質(zhì)也可對沉積物細(xì)菌群落結(jié)構(gòu)有一定的影響.營養(yǎng)狀況可能對淡水湖泊沉積細(xì)菌群落的形成起重要作用[35].參考《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB38382-2002)[36]Ⅲ類水標(biāo)準(zhǔn)要求,河口區(qū)和湖心區(qū)的TN和TP含量都已超標(biāo),由氮磷含量反應(yīng)出太湖此時已呈現(xiàn)營養(yǎng)狀態(tài).冗余分析結(jié)果表明水環(huán)境因子顯著地影響著沉積物細(xì)菌群落.Winters等[37]發(fā)現(xiàn)硝酸鹽、銨鹽和磷酸鹽等環(huán)境因子均可影響勞倫森大湖沉積物細(xì)菌群落組成.通過考察不同污染源對太湖春夏秋冬四季沉積物細(xì)菌群落的影響,Jing等[38]發(fā)現(xiàn)銨鹽和有機(jī)物顯著影響細(xì)菌群落組成,4種不同的污染源(污水處理廠,城市,郊區(qū),農(nóng)業(yè))對太湖沉積物細(xì)菌群落影響不同,其中污水處理廠最可能導(dǎo)致沉積物生物同質(zhì)化.

    Microcystis在太湖9個采樣點均有檢出,上覆水的水溫和NO3--N是冬季沉積物中Microcystis的主要影響因子.Liu等[39]對太湖進(jìn)行了為期11年的調(diào)查工作,發(fā)現(xiàn)Microcystis與NOx-N和水溫呈正相關(guān),在含氮化合物濃度下降時,水溫升高是微囊藻水華暴發(fā)的主要推動力.李朋富教授團(tuán)隊對太湖梅梁灣、貢湖灣和湖心區(qū)進(jìn)行為期6個月的研究工作,發(fā)現(xiàn)從春季到秋季Microcystis含量持續(xù)升高[40].本次Microcystis在太湖不同湖區(qū)均有檢出,當(dāng)水溫持續(xù)上升后,Microcystis可能會迅速增長以致大規(guī)模積聚.秦伯強(qiáng)等[41]發(fā)現(xiàn)產(chǎn)毒的Microcystis主要存在于太湖營養(yǎng)充足的湖區(qū),而無毒的Microcystis則是營養(yǎng)受限湖區(qū)的優(yōu)勢細(xì)菌.本次研究中,Microcystis在河口區(qū)和湖心區(qū)豐度很高,這與河口區(qū)和湖心區(qū)的營養(yǎng)鹽含量密不可分;同時河口區(qū)與湖心區(qū)的水溫也高于水草區(qū).對于存在于河口區(qū)和湖心區(qū)等營養(yǎng)水平較高湖區(qū)的Microcystis,應(yīng)給予極大關(guān)注.通過找出Microcystis等有害生物的制約因子,以控制環(huán)境制約因子達(dá)到對有害微囊藻水華的防控.

    4 結(jié)論

    4.1 營養(yǎng)程度不同,湖區(qū)物種豐富度和均勻度也不同,3個湖區(qū)的物種豐富度和均勻度為水草區(qū)>湖心區(qū)>河口區(qū).

    4.2 不同湖區(qū)優(yōu)勢細(xì)菌群落組成存在差異,Nitrospirae、Acidobacteria、Anaerolineae和Nitrospira占比在河口區(qū)、湖心區(qū)和水草區(qū)之間差異顯著(P<0.05).冬季太湖沉積物細(xì)菌主要類群為Chloroflexi(22.18%)、Proteobacteria(21.11%)、Cyanobacteria(13.37%)和Nitrospirae(10.17%);在屬分類水平上,主要的優(yōu)勢類群為unidentified_Chloroplast(4.78%)和Microcystis(3.81%).

    4.3 冬季太湖沉積物中全面檢出Cyanobacteria,且含量豐富.隸屬于Cyanobacteria的Microcystis在太湖不同湖區(qū)9個采樣點也均有檢出,NO3--N和水溫是Microcystis的主要影響因子,可為早期藍(lán)藻和有害微囊藻水華防控工作提供基礎(chǔ)數(shù)據(jù).

    [1] Qin B Q, Zhu G W, Gao G A, et al. A drinking water crisis in Lake Taihu, China: linkage to climatic variability and lake management [J]. Environmental Management, 2010,45(1):105-112.

    [2] Ye W J, Liu X L, Lin S Q, et al. The vertical distribution of bacterial and archaeal communities in the water and sediment of Lake Taihu [J]. Fems Microbiology Ecology, 2009,70(2):263-276.

    [3] Yang Z, Zhang M, Shi X L, et al. Nutrient reduction magnifies the impact of extreme weather on cyanobacterial bloom formation in large shallow Lake Taihu (China) [J]. Water Research, 2016,103:302-310.

    [4] 鄭艷玲,侯立軍,陸 敏,等.崇明東灘夏冬季表層沉積物細(xì)菌多樣性研究 [J]. 中國環(huán)境科學(xué), 2012,32(2):300-310.

    [5] 顧婷婷,孔繁翔,譚 嘯,等.越冬和復(fù)蘇時期太湖水體藍(lán)藻群落結(jié)構(gòu)的時空變化 [J]. 生態(tài)學(xué)報, 2011,31(1):21-30.

    [6] Tang X M, Gao G A, Chao J Y, et al. Dynamics of organicaggregate-associated bacterial communities and related environmental factors in Lake Taihu, a large eutrophic shallow lake in China [J]. Limnology and Oceanography, 2010,55(2):469-480.

    [7] Shi L M, Cai Y F, Kong F X, et al. Specific association between bacteria and buoyant Microcystis, colonies compared with other bulk bacterial communities in the eutrophic Lake Taihu, China [J].Environmental Microbiology Reports, 2012,4(6):669-678.

    [8] 陳兆進(jìn),陳海燕,李玉英,等.南水北調(diào)中線干渠(河南段)浮游細(xì)菌群落組成及影響因素 [J]. 中國環(huán)境科學(xué), 2017,37(4):1505-1513.

    [9] Chao Y Q, Mao Y P, Wang Z P, et al. Diversity and functions of bacterial community in drinking water biofilms revealed by high-throughput sequencing [J]. Scientific Reports, 2015,5:10044.

    [10] Sun Y J, Wang T Y, Peng X W, et al. Bacterial community compositions in sediment polluted by perfluoroalkyl acids(PFAAs) using Illumina high-throughput sequencing [J].Environmental Science and Pollution Research, 2016,23(11):10556-10565.

    [11] 劉 洋,曾全超,黃懿梅.基于454高通量測序的黃土高原不同喬木林土壤細(xì)菌群落特征 [J]. 中國環(huán)境科學(xué), 2016,36(11):3487-3494.

    [12] 白 潔,李海艷,張 健,等.黃海西北部沉積物中細(xì)菌群落16S rDNA多樣性解析 [J]. 中國環(huán)境科學(xué), 2009,29(12):1277-1284.

    [13] Shao K Q, Gao G, Qin B Q, et al. Structural changes of the sediment bacterial community in Meiliang Bay of Lake Taihu,China, during typhoon Morakot, 2009 [J]. Journal of Basic Microbiology, 2012,52(1):99-103.

    [14] Tang X M, Li L L, Shao K Q, et al. Pyrosequencing analysis of free-living and attached bacterial communities in Meiliang Bay,Lake Taihu, a large eutrophic shallow lake in China [J]. Canadian Journal of Microbiology, 2014,61(1):22-31.

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

    [16] Desantis T Z, Hugenholtz P, Larsen N, et al. Greengenes: a chimera-checked 16S rRNA gene database and workbench compatible in ARB [J]. Applied & Environmental Microbiology,2006,72(7):5069-5072.

    [17] 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.

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

    [19] Huse S H, Welch D M, Morrison H G, et al. Ironing out the wrinkles in the rare biosphere through improved OTU clustering[J]. Environmental Microbiology, 2010,12(7):1889-1898.

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

    [21] Zhou J, Jiang Y H, Deng Y, et al. Random sampling process leads to overestimation of β-diversity of microbial communities [J].Mbio, 2013,4(3):e00324-13.

    [22] 逄 勇,顏潤潤,余鐘波,等.風(fēng)浪作用下的底泥懸浮沉降及內(nèi)源釋放量研究 [J]. 環(huán)境科學(xué), 2008,29(9):2456-2464.

    [23] Strecker A L, Arnott S E, Yan N D, et al. Variation in the response of crustacean zooplankton species richness and composition to the invasive predator Bythotrephes longimanus [J]. Canadian Journal of Fisheries and Aquatic Sciences, 2006,63(9):2126-2136.

    [24] Sand-Jensen K, Borum J. Interactions among phytoplankton,periphyton, and macrophytes in temperate freshwaters and estuaries [J]. Aquatic Botany, 1991,41(1-3):137-175.

    [25] Liu F H, Lin G H, Gao G, et al. Bacterial and archaeal assemblages in sediments of a large shallow freshwater lake, Lake Taihu, as revealed by denaturing gradient gel electrophoresis [J].Journal of Applied Microbiology, 2009,106(3):1022-1032.

    [26] Nam Y D, Sung Y, Chang H W, et al. Characterization of the depth-related changes in the microbial communities in Lake Hovsgol sediment by 16S rRNA gene-based approaches [J].Journal of Microbiology, 2008,46(2):125-136.

    [27] Madsen E L. Microorganisms and their roles in fundamental biogeochemical cycles [J]. Current Opinion in Biotechnology,2011,22(3):456-464.

    [28] Shao K Q, Gao G, Qin B Q, et al. Comparing sediment bacterial communities in the macrophyte-dominated and algae-dominated areas of eutrophic Lake Taihu, China [J]. Canadian Journal of Microbiology, 2011,57(4):263-272.

    [29] Huang R, Zhao D Y, Jiang C L, et al. Heterogeneity of bacterial community compositions in surface sediments of three lake zones in Lake Taihu [J]. Proceedings of the National Academy of Sciences, India Section B: Biological Sciences, 2015,85(2):1-10.[30] Yang Z, Zhang M, Shi X L, et al. Nutrient reduction magnifies the impact of extreme weather on cyanobacterial bloom formation in large shallow Lake Taihu (China) [J]. Water Research, 2016,103:302-310.

    [31] Paerl H W, Huisman J. Blooms like it hot [J]. Science,2008,320(5872):57-58.

    [32] O’Neil J M, Davis T W, Burford M A, et al. The rise of harmful cyanobacteria blooms: The potential roles of eutrophication and climate change [J]. Harmful Algae, 2012,14:313-334.

    [33] 宜興日報.太湖藍(lán)藻打撈量是去年同期的5.5倍 [EB/OL].http://www.yixing.gov.cn/zxzx/show-13593314.heml2017-05-16/2017-06-26.

    [34] 陳 楠.太湖沉積物微生物群落組成與物質(zhì)循環(huán)及藍(lán)藻爆發(fā)的相關(guān)性 [D]. 北京:中國農(nóng)業(yè)大學(xué), 2015.

    [35] Yu D, Yang Y, Zhen W, 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.

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

    [37] Winters A D, Marsh T L, Brenden T O, et al. Molecular characterization of bacterial communities associated with sediments in the Laurentian Great Lakes [J]. Journal of Great Lakes Research, 2014,40(3):640-645.

    [38] Jing W, Yi L, Wang P F, et al. Response of bacterial community compositions to different sources of pollutants in sediments of a tributary of Taihu Lake, China [J]. Environmental Science and Pollution Research International, 2016,23(14):13886-13894.

    [39] Liu X, Lu X H, Chen Y W. The effects of temperature and nutrient ratios on Microcystis blooms in Lake Taihu, China: an 11-year investigation [J]. Harmful Algae, 2011,10(3):337-343.

    [40] Wang X Y, Sun M J, Wang J M, et al. Microcystis genotype succession and related environmental factors in Lake Taihu during cyanobacterial blooms [J]. Microbial Ecology, 2012,64(4):986-999.

    [41] Otten T G, Xu H, Qin B Q, et al. Spatiotemporal patterns and ecophysiology of toxigenic microcystis blooms in Lake Taihu,China: implications for water quality management [J].Environmental Science and Technology, 2016,46(6):3480-3488.

    猜你喜歡
    河口區(qū)湖心湖區(qū)
    白天鵝
    漢江河口區(qū)四明灘段四明上護(hù)岸坍岸修復(fù)技術(shù)
    大通湖區(qū)河蟹產(chǎn)業(yè)發(fā)展綜述
    白天鵝
    幼兒園(2018年15期)2018-10-15 19:40:36
    生活在湖區(qū)
    海峽旅游(2018年4期)2018-06-01 11:20:00
    湖區(qū)航道風(fēng)速預(yù)警監(jiān)測點布設(shè)研究
    江西建材(2018年4期)2018-04-10 12:37:24
    “傳承—信仰—未來”教育主題的設(shè)計者
    ——記山東省東營市河口區(qū)河安小學(xué)校長劉波
    教書育人(2017年23期)2017-09-26 06:47:50
    尋找你的伙伴
    湖心哈哈鏡
    閩江下游及河口區(qū)表層沉積物多種微量元素特征
    亚洲国产欧美一区二区综合| 91麻豆av在线| xxxhd国产人妻xxx| 一区二区三区乱码不卡18| 欧美精品亚洲一区二区| 动漫黄色视频在线观看| 国产野战对白在线观看| 国产三级黄色录像| 高清毛片免费观看视频网站 | 丰满饥渴人妻一区二区三| 国产高清视频在线播放一区| 国产成人欧美| 精品人妻在线不人妻| 91成年电影在线观看| 80岁老熟妇乱子伦牲交| av不卡在线播放| 日韩欧美一区视频在线观看| 男人舔女人的私密视频| 天天影视国产精品| www.自偷自拍.com| 久久精品亚洲熟妇少妇任你| 黑人欧美特级aaaaaa片| 一级片'在线观看视频| 香蕉国产在线看| 大香蕉久久成人网| 亚洲av国产av综合av卡| 天天操日日干夜夜撸| www.999成人在线观看| 一级毛片精品| 两性午夜刺激爽爽歪歪视频在线观看 | 成年女人毛片免费观看观看9 | 国产av又大| 婷婷丁香在线五月| 久久天堂一区二区三区四区| 亚洲视频免费观看视频| 久久国产精品大桥未久av| 久久中文字幕人妻熟女| 成年人黄色毛片网站| 自拍欧美九色日韩亚洲蝌蚪91| 欧美一级毛片孕妇| 国产男女超爽视频在线观看| 最近最新免费中文字幕在线| 热99re8久久精品国产| 精品一区二区三区四区五区乱码| 一区二区av电影网| 国产亚洲午夜精品一区二区久久| 多毛熟女@视频| 久久精品亚洲精品国产色婷小说| 999久久久国产精品视频| 精品少妇黑人巨大在线播放| 母亲3免费完整高清在线观看| 香蕉国产在线看| 亚洲免费av在线视频| 青草久久国产| 亚洲九九香蕉| 视频在线观看一区二区三区| 精品久久久久久电影网| 国产亚洲欧美精品永久| 咕卡用的链子| 国产在视频线精品| av线在线观看网站| 美女高潮到喷水免费观看| 99国产综合亚洲精品| 欧美精品一区二区免费开放| 国产男女内射视频| 少妇猛男粗大的猛烈进出视频| 国产精品久久久久久精品电影小说| 纯流量卡能插随身wifi吗| 日韩欧美一区视频在线观看| 亚洲欧美一区二区三区黑人| 老司机影院毛片| 18在线观看网站| 亚洲精品乱久久久久久| 亚洲国产欧美网| 中文亚洲av片在线观看爽 | 欧美一级毛片孕妇| 午夜激情久久久久久久| 搡老熟女国产l中国老女人| 亚洲三区欧美一区| 国产欧美日韩一区二区三| 日韩三级视频一区二区三区| 男女边摸边吃奶| 电影成人av| 国产单亲对白刺激| 国产成人免费无遮挡视频| 一本久久精品| 免费在线观看日本一区| 黄色怎么调成土黄色| 啦啦啦中文免费视频观看日本| 亚洲色图av天堂| 97人妻天天添夜夜摸| 久久精品国产a三级三级三级| 蜜桃在线观看..| 国产成人精品在线电影| 国产黄色免费在线视频| 免费久久久久久久精品成人欧美视频| 亚洲国产成人一精品久久久| 考比视频在线观看| 精品国产乱码久久久久久小说| 欧美国产精品一级二级三级| 中文字幕av电影在线播放| 国产精品久久久久久精品古装| 国产精品国产av在线观看| 最近最新中文字幕大全免费视频| 国产成人免费观看mmmm| 精品福利观看| 99国产精品免费福利视频| 国产欧美日韩一区二区精品| 国产精品电影一区二区三区 | 18禁国产床啪视频网站| 国产国语露脸激情在线看| 丝袜美足系列| 老熟妇仑乱视频hdxx| www.999成人在线观看| 中文欧美无线码| 国产高清国产精品国产三级| 又紧又爽又黄一区二区| 宅男免费午夜| 中文欧美无线码| 精品第一国产精品| 黄色 视频免费看| 久久天堂一区二区三区四区| 久久九九热精品免费| 9色porny在线观看| 亚洲熟女毛片儿| 一进一出好大好爽视频| a在线观看视频网站| 成人18禁在线播放| 亚洲三区欧美一区| 考比视频在线观看| 黄片大片在线免费观看| 精品国产乱子伦一区二区三区| 亚洲精品成人av观看孕妇| 人人妻,人人澡人人爽秒播| 久久精品人人爽人人爽视色| 两人在一起打扑克的视频| 热re99久久国产66热| 午夜精品国产一区二区电影| 久久久国产一区二区| 九色亚洲精品在线播放| 欧美变态另类bdsm刘玥| 久久香蕉激情| 精品卡一卡二卡四卡免费| 久久人人97超碰香蕉20202| 亚洲精品粉嫩美女一区| 99热国产这里只有精品6| 午夜成年电影在线免费观看| 婷婷成人精品国产| 91成年电影在线观看| 国产一区二区在线观看av| 亚洲成人手机| 无遮挡黄片免费观看| 免费女性裸体啪啪无遮挡网站| 日本vs欧美在线观看视频| 亚洲av片天天在线观看| 黄色 视频免费看| 在线天堂中文资源库| 亚洲色图 男人天堂 中文字幕| 久久午夜综合久久蜜桃| 国产无遮挡羞羞视频在线观看| 18禁国产床啪视频网站| 欧美日韩视频精品一区| 一级毛片女人18水好多| 乱人伦中国视频| 午夜福利在线免费观看网站| 国产精品秋霞免费鲁丝片| 亚洲成人手机| 午夜福利视频在线观看免费| 天堂中文最新版在线下载| 极品人妻少妇av视频| av国产精品久久久久影院| 老熟妇仑乱视频hdxx| av在线播放免费不卡| 国产精品 国内视频| 国产高清国产精品国产三级| 12—13女人毛片做爰片一| 精品少妇黑人巨大在线播放| 久久久久久亚洲精品国产蜜桃av| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 又大又爽又粗| 最新美女视频免费是黄的| 精品国产亚洲在线| 黄色视频不卡| 18禁裸乳无遮挡动漫免费视频| 国产老妇伦熟女老妇高清| 日韩精品免费视频一区二区三区| 日韩人妻精品一区2区三区| 成人国语在线视频| 成年人午夜在线观看视频| 巨乳人妻的诱惑在线观看| 夜夜爽天天搞| 两个人看的免费小视频| 亚洲欧美激情在线| 国产伦理片在线播放av一区| 亚洲天堂av无毛| 真人做人爱边吃奶动态| 亚洲第一av免费看| 91国产中文字幕| 国产成人免费无遮挡视频| 不卡一级毛片| 视频区欧美日本亚洲| 免费不卡黄色视频| www日本在线高清视频| 美女午夜性视频免费| 一级片免费观看大全| 搡老乐熟女国产| 一区二区三区激情视频| 久久亚洲精品不卡| 成年人黄色毛片网站| 国产精品电影一区二区三区 | 欧美成人午夜精品| 国产三级黄色录像| 中文字幕色久视频| 午夜两性在线视频| 成年人免费黄色播放视频| 国产成人av教育| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 免费观看av网站的网址| 欧美久久黑人一区二区| videos熟女内射| 日韩大码丰满熟妇| 国产在视频线精品| 视频区欧美日本亚洲| 国产高清videossex| 老熟女久久久| 这个男人来自地球电影免费观看| 在线十欧美十亚洲十日本专区| 精品人妻熟女毛片av久久网站| 狂野欧美激情性xxxx| 蜜桃国产av成人99| 在线看a的网站| av国产精品久久久久影院| 久久精品熟女亚洲av麻豆精品| 最黄视频免费看| 曰老女人黄片| 黄色视频,在线免费观看| 国产高清国产精品国产三级| 精品视频人人做人人爽| 亚洲成国产人片在线观看| 国产一区二区三区视频了| 国产成人精品无人区| 国产99久久九九免费精品| 国产不卡一卡二| 操出白浆在线播放| 精品国产超薄肉色丝袜足j| 亚洲欧美一区二区三区黑人| 欧美日韩国产mv在线观看视频| 国产一卡二卡三卡精品| 国产激情久久老熟女| 狠狠精品人妻久久久久久综合| 久久99一区二区三区| 黄色a级毛片大全视频| 精品国产一区二区三区四区第35| 亚洲国产av影院在线观看| 亚洲熟女毛片儿| 午夜福利乱码中文字幕| 少妇裸体淫交视频免费看高清 | 亚洲少妇的诱惑av| 中文字幕人妻熟女乱码| 国产日韩欧美视频二区| 久久天躁狠狠躁夜夜2o2o| 欧美日韩国产mv在线观看视频| 深夜精品福利| av网站在线播放免费| 黄色 视频免费看| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| 久久久久久久国产电影| 波多野结衣一区麻豆| 亚洲色图综合在线观看| 国产真人三级小视频在线观看| 国产黄频视频在线观看| 一二三四在线观看免费中文在| 亚洲性夜色夜夜综合| 亚洲精品久久午夜乱码| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 少妇猛男粗大的猛烈进出视频| 1024香蕉在线观看| 午夜成年电影在线免费观看| 国产不卡一卡二| 色婷婷av一区二区三区视频| 亚洲av成人不卡在线观看播放网| 大香蕉久久成人网| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 一级a爱视频在线免费观看| 久久中文字幕人妻熟女| 在线观看舔阴道视频| 精品人妻熟女毛片av久久网站| 黄色丝袜av网址大全| 亚洲国产毛片av蜜桃av| 别揉我奶头~嗯~啊~动态视频| 少妇 在线观看| 免费观看av网站的网址| 一本大道久久a久久精品| 成人免费观看视频高清| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品古装| 乱人伦中国视频| 悠悠久久av| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 免费少妇av软件| 久久久久国内视频| 日日爽夜夜爽网站| 国产日韩欧美在线精品| 日本av免费视频播放| 久久精品亚洲熟妇少妇任你| 国产一区二区三区综合在线观看| 亚洲精品国产一区二区精华液| 纯流量卡能插随身wifi吗| 最近最新中文字幕大全免费视频| xxxhd国产人妻xxx| 国产日韩欧美在线精品| 久久九九热精品免费| 精品卡一卡二卡四卡免费| 丁香欧美五月| 这个男人来自地球电影免费观看| 国产精品久久久久久精品电影小说| 国产亚洲精品久久久久5区| 99国产精品99久久久久| 色94色欧美一区二区| 99精品在免费线老司机午夜| 99精品欧美一区二区三区四区| 色尼玛亚洲综合影院| 精品国产一区二区三区四区第35| 大片电影免费在线观看免费| 国产不卡一卡二| 国产精品亚洲av一区麻豆| 日韩制服丝袜自拍偷拍| 国产成人av教育| 老熟妇仑乱视频hdxx| 免费看十八禁软件| 亚洲中文字幕日韩| 黄片大片在线免费观看| 亚洲精品国产区一区二| 成人黄色视频免费在线看| 国产亚洲一区二区精品| 两性夫妻黄色片| 黄网站色视频无遮挡免费观看| 电影成人av| 成年版毛片免费区| av超薄肉色丝袜交足视频| 国产精品熟女久久久久浪| 肉色欧美久久久久久久蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品一区二区三区| 亚洲九九香蕉| 欧美 亚洲 国产 日韩一| 国产精品自产拍在线观看55亚洲 | 韩国精品一区二区三区| aaaaa片日本免费| 亚洲专区中文字幕在线| 国产一区二区在线观看av| 日韩欧美国产一区二区入口| 操美女的视频在线观看| videos熟女内射| 99久久精品国产亚洲精品| 99国产极品粉嫩在线观看| 精品国产一区二区久久| 人人妻,人人澡人人爽秒播| 国产一区二区在线观看av| 乱人伦中国视频| 极品人妻少妇av视频| 亚洲精品在线观看二区| 久久久久视频综合| 国产精品香港三级国产av潘金莲| 国产男女内射视频| 蜜桃在线观看..| 制服人妻中文乱码| 一本一本久久a久久精品综合妖精| 制服人妻中文乱码| 老汉色av国产亚洲站长工具| 中文字幕高清在线视频| 国产精品自产拍在线观看55亚洲 | 免费黄频网站在线观看国产| 777米奇影视久久| 国产一区有黄有色的免费视频| 亚洲精品乱久久久久久| 亚洲国产毛片av蜜桃av| 亚洲国产中文字幕在线视频| 国产欧美日韩一区二区三| 制服人妻中文乱码| 亚洲熟妇熟女久久| 精品国产国语对白av| 天天添夜夜摸| www日本在线高清视频| 69av精品久久久久久 | 久久午夜亚洲精品久久| 狠狠精品人妻久久久久久综合| 黄色片一级片一级黄色片| 男人舔女人的私密视频| 成人国产av品久久久| 无人区码免费观看不卡 | av网站在线播放免费| av在线播放免费不卡| 国产精品欧美亚洲77777| 欧美变态另类bdsm刘玥| 午夜免费鲁丝| 人人妻人人澡人人爽人人夜夜| 国产成人精品久久二区二区免费| 国产伦人伦偷精品视频| 亚洲欧美精品综合一区二区三区| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区视频在线观看免费 | 不卡一级毛片| 久久这里只有精品19| av天堂在线播放| 日本wwww免费看| 精品人妻1区二区| av片东京热男人的天堂| 在线观看一区二区三区激情| 国产精品免费大片| 日韩欧美一区二区三区在线观看 | 久久国产亚洲av麻豆专区| 国产亚洲欧美在线一区二区| 亚洲国产欧美网| 黄频高清免费视频| 超碰成人久久| 久久久久国内视频| 午夜视频精品福利| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产99精品国产亚洲性色 | 国产精品欧美亚洲77777| 一个人免费看片子| 亚洲人成电影免费在线| 曰老女人黄片| 交换朋友夫妻互换小说| 亚洲专区字幕在线| 91麻豆精品激情在线观看国产 | 日韩 欧美 亚洲 中文字幕| 国产欧美日韩综合在线一区二区| 色在线成人网| 色尼玛亚洲综合影院| 亚洲精品中文字幕在线视频| 99精品久久久久人妻精品| 亚洲 国产 在线| 老司机午夜十八禁免费视频| 久久精品熟女亚洲av麻豆精品| av天堂久久9| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 欧美在线黄色| 麻豆国产av国片精品| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻丝袜一区二区| 国产成人av激情在线播放| a级片在线免费高清观看视频| 国产亚洲午夜精品一区二区久久| 国产在线一区二区三区精| 亚洲专区中文字幕在线| 一区二区av电影网| av在线播放免费不卡| 久久人妻av系列| 美女高潮喷水抽搐中文字幕| 99久久99久久久精品蜜桃| 亚洲国产精品一区二区三区在线| 亚洲三区欧美一区| 亚洲人成伊人成综合网2020| 成年动漫av网址| 黄片播放在线免费| 黄色视频,在线免费观看| 欧美激情久久久久久爽电影 | 成人av一区二区三区在线看| 国产97色在线日韩免费| 日韩欧美一区二区三区在线观看 | 夜夜骑夜夜射夜夜干| 纵有疾风起免费观看全集完整版| 日韩中文字幕欧美一区二区| 国产成人一区二区三区免费视频网站| 亚洲美女黄片视频| svipshipincom国产片| 一本一本久久a久久精品综合妖精| 亚洲成人手机| 久久99一区二区三区| 天堂俺去俺来也www色官网| 国产成人av教育| 美女午夜性视频免费| 国产一区二区 视频在线| 色94色欧美一区二区| 色老头精品视频在线观看| 大香蕉久久成人网| 女人高潮潮喷娇喘18禁视频| 亚洲欧洲日产国产| 在线天堂中文资源库| 超碰成人久久| 国产色视频综合| 十八禁高潮呻吟视频| 夫妻午夜视频| 欧美日本中文国产一区发布| 99国产精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 嫩草影视91久久| 亚洲综合色网址| 欧美激情 高清一区二区三区| 午夜日韩欧美国产| 老鸭窝网址在线观看| 婷婷成人精品国产| 精品久久蜜臀av无| 亚洲,欧美精品.| 99九九在线精品视频| 最近最新免费中文字幕在线| 999久久久国产精品视频| 老司机福利观看| 人妻一区二区av| 亚洲国产欧美网| 中亚洲国语对白在线视频| 久久久久久久国产电影| 国产激情久久老熟女| 在线天堂中文资源库| 老司机亚洲免费影院| tocl精华| 一个人免费看片子| 91老司机精品| 国产有黄有色有爽视频| 国产亚洲一区二区精品| 亚洲精品乱久久久久久| 日本精品一区二区三区蜜桃| 欧美乱码精品一区二区三区| 十八禁人妻一区二区| 嫩草影视91久久| 免费av中文字幕在线| 国产欧美亚洲国产| 手机成人av网站| 欧美人与性动交α欧美精品济南到| 亚洲欧美色中文字幕在线| 天堂8中文在线网| 免费女性裸体啪啪无遮挡网站| 精品国产亚洲在线| 久久久久视频综合| 国产欧美日韩精品亚洲av| 一本大道久久a久久精品| 国产色视频综合| 欧美精品高潮呻吟av久久| 人人澡人人妻人| 国内毛片毛片毛片毛片毛片| 妹子高潮喷水视频| 日韩中文字幕视频在线看片| 深夜精品福利| 在线十欧美十亚洲十日本专区| 久9热在线精品视频| 精品国产一区二区三区四区第35| 窝窝影院91人妻| 满18在线观看网站| 日韩免费av在线播放| 多毛熟女@视频| 一个人免费在线观看的高清视频| 亚洲精品成人av观看孕妇| 亚洲成a人片在线一区二区| 亚洲,欧美精品.| 欧美人与性动交α欧美软件| 波多野结衣一区麻豆| 午夜福利欧美成人| 亚洲色图综合在线观看| 夜夜爽天天搞| 日日夜夜操网爽| 亚洲av欧美aⅴ国产| 精品一区二区三区视频在线观看免费 | 欧美亚洲日本最大视频资源| 最黄视频免费看| 国产野战对白在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线观看jvid| av网站在线播放免费| 国产一卡二卡三卡精品| 亚洲精品久久午夜乱码| 91成年电影在线观看| 91九色精品人成在线观看| 免费女性裸体啪啪无遮挡网站| 下体分泌物呈黄色| 日本一区二区免费在线视频| 国产不卡av网站在线观看| 夜夜骑夜夜射夜夜干| 如日韩欧美国产精品一区二区三区| 午夜福利视频在线观看免费| 久久青草综合色| a级毛片在线看网站| 精品亚洲成a人片在线观看| 青青草视频在线视频观看| 亚洲一码二码三码区别大吗| 亚洲国产成人一精品久久久| av一本久久久久| 99国产综合亚洲精品| 欧美亚洲日本最大视频资源| 法律面前人人平等表现在哪些方面| 午夜福利免费观看在线| 国产亚洲av高清不卡| 精品久久蜜臀av无| 女性生殖器流出的白浆| 亚洲中文字幕日韩| 人妻久久中文字幕网| 精品久久久久久久毛片微露脸| 久久久精品区二区三区| 高清av免费在线| 天天添夜夜摸| 久久婷婷成人综合色麻豆| 欧美国产精品一级二级三级| 叶爱在线成人免费视频播放| 亚洲国产欧美日韩在线播放| 丰满人妻熟妇乱又伦精品不卡| 免费观看人在逋| 亚洲全国av大片| 两个人看的免费小视频| 纵有疾风起免费观看全集完整版| 人人妻,人人澡人人爽秒播| 高潮久久久久久久久久久不卡| cao死你这个sao货| av线在线观看网站| 亚洲欧美日韩高清在线视频 | 老汉色∧v一级毛片|