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

    湖泛發(fā)生過程中水體理化性質(zhì)及細(xì)菌群落的變化

    2017-12-29 01:49:58馮勝袁斌王博雯
    生態(tài)科學(xué) 2017年6期
    關(guān)鍵詞:放線菌菌門太湖

    馮勝, 袁斌, 王博雯

    ?

    湖泛發(fā)生過程中水體理化性質(zhì)及細(xì)菌群落的變化

    馮勝1,*, 袁斌1, 王博雯2

    1. 常州大學(xué)環(huán)境與安全工程學(xué)院, 江蘇 常州 213164 2. 中國科學(xué)院南京地理與湖泊研究所, 湖泊與環(huán)境國家重點實驗室, 江蘇 南京 210008

    運用分子生物學(xué)技術(shù)及多元統(tǒng)計分析的方法, 跟蹤湖泛爆發(fā)過程中細(xì)菌群落的動態(tài)及環(huán)境因子變化。結(jié)果表明: 湖泛爆發(fā)過程中水質(zhì)明顯惡化, DO(溶解氧)、SD(透明度)、pH(酸堿度)顯著降低, 其中DO含量均低于0.51 mg·L-1, 水體中TN(總氮)、TP(總磷)等營養(yǎng)鹽的含量遠(yuǎn)遠(yuǎn)高于正常水體中的含量。細(xì)菌豐度和藻類數(shù)量明顯升高, 其中, 湖泛中心區(qū)的細(xì)菌豐度上升到107cells·L-1, 相比未發(fā)生湖泛的區(qū)域提高了一個數(shù)量級。運用高通量測序技術(shù), 分析比較非湖泛區(qū)及不同程度湖泛區(qū)的細(xì)菌群落結(jié)構(gòu)的差異。結(jié)果表明: 湖泛最嚴(yán)重的區(qū)域, α-、γ-變形菌綱(alpha-, gamma-Proteobacteria)分別由未發(fā)生湖泛時的4.10%和13.16%降至0.59%和1.91%; 放線菌門(Actinobacteria)所占比例由未發(fā)生湖泛區(qū)域的19.04%降至0.92%; 而擬桿菌門(Bacteriodetes)所占的比例則由1.12%增加至36.32%。此外, 還檢測到一些特殊生境的細(xì)菌類群, 如。通過對湖泛事件中細(xì)菌群落與關(guān)鍵環(huán)境因子的變化規(guī)律研究, 為湖泛問題的深入研究提供一定科學(xué)依據(jù)。

    水環(huán)境因子; 細(xì)菌豐度; 細(xì)菌群落結(jié)構(gòu); 湖泛

    1 前言

    太湖是我國五大淡水湖之一, 自20世紀(jì)90年代起, 太湖湖西水域、梅梁湖、貢湖曾多次出現(xiàn)湖泛現(xiàn)象[1-2]。湖泛是藻類大量暴發(fā)、積聚和死亡后, 在適宜的氣象、水文條件下與底泥中的有機物發(fā)生生物、化學(xué)反應(yīng), 在缺氧或厭氧條件下釋放硫化物、甲烷和二甲基三硫等硫醚類物質(zhì), 從而在富營養(yǎng)化湖泊中形成伴有惡臭的褐黑色水團(tuán)的現(xiàn)象[3]。繼大規(guī)模重度藍(lán)藻水華之后, “湖泛”問題已成為我國淡水富營養(yǎng)化淺水湖泊面臨的又一嚴(yán)峻生態(tài)災(zāi)害[4]。目前針對湖泛的研究涉及湖泊水文[5]、氣象[6]、物理[7]、化學(xué)[8-9]、生物[10]等多個領(lǐng)域, 大量交叉研究的結(jié)果顯示: 藍(lán)藻和某些快速繁殖水生植物的堆積與污染有機質(zhì)的輸入是湖泛發(fā)生的物質(zhì)基礎(chǔ), 適宜的氣象水文條件(強光照、高溫、低風(fēng)速、水體垂直熱分層等)是湖泛的觸發(fā)因素, 而以厭氧條件下快速分解有機質(zhì)的梭菌、硫酸鹽還原菌、產(chǎn)甲烷菌(古菌)以及甲烷氧化菌為標(biāo)志的微生物是湖泛形成的生物機制。溶解氧作為聯(lián)系湖泊水文、理化和生物過程的關(guān)鍵因子, 成為誘導(dǎo)湖泛形成和消散的核心因素[11-12]。由此可見: 與湖泛發(fā)生時的環(huán)境特征(如低溶解氧, 低pH, 高有機質(zhì), 高總磷、總氮)相對應(yīng)的是其簡化的食物網(wǎng)結(jié)構(gòu)和特殊的微生物類群, 細(xì)菌作為湖泊生態(tài)系統(tǒng)中食物鏈和食物網(wǎng)的重要組成部分, 驅(qū)動湖泊生態(tài)系統(tǒng)中絕大多數(shù)生物活性元素的形態(tài)轉(zhuǎn)化和地球化學(xué)循環(huán)[13], 在湖泛發(fā)生期間及后期的生態(tài)系統(tǒng)恢復(fù)等方面都具有關(guān)鍵作用。細(xì)菌是湖泊生態(tài)系統(tǒng)中的主要分解者, 在物質(zhì)循環(huán)和能量流動中起著重要作用。其數(shù)量和群落結(jié)構(gòu)的改變都會顯著影響整個微食物網(wǎng)的結(jié)構(gòu)及系統(tǒng)中物質(zhì)循環(huán)和能量流動, 對于細(xì)菌群落結(jié)構(gòu)特征的研究受到了廣泛關(guān)注[14-15]。李玉華等[16]采用聚合酶鏈?zhǔn)椒磻?yīng)-變性凝膠梯度電泳(PCR-DGGE)對松花湖水體中微生物群落結(jié)構(gòu)的多樣性進(jìn)行了研究; 沈燕等[17]選用針對α、β、γ變形菌門和擬桿菌門的幾類探針, 采用原位雜交的方法(FISH)研究了太湖中主要變形菌門細(xì)菌的優(yōu)勢類群、空間分布差異及其與營養(yǎng)鹽、Chl-a濃度等環(huán)境因素的相互影響。吳慶龍等[18]運用原位熒光雜交、變形梯度凝膠電泳(DGGE)和反向線點雜交(RLB)等多種分子生物學(xué)方法對青藏高原16個不同鹽度湖區(qū)進(jìn)行了鹽度和緯度對細(xì)菌群落結(jié)構(gòu)影響的研究。目前使用最多的分子生物學(xué)方法主要為DGGE和T-RFLP, 有些研究也會再結(jié)合Sanger測序得到一些特別的序列圖譜的堿基序列, 建立克隆文庫與數(shù)據(jù)庫進(jìn)行對比[19-22]; 而高通量測序技術(shù), 無需構(gòu)建文庫, 避免了克隆的誤差, 能在數(shù)百萬個點上同時閱讀測序, 擁有較好的定量功能和性價比[23-24]。本研究應(yīng)用高通量測序技術(shù)研究湖泛發(fā)生過程中細(xì)菌群落組成的動態(tài)變化, 以探尋湖泛事件過程中細(xì)菌群落的變化規(guī)律, 為湖泛問題的研究提供一定科學(xué)依據(jù)。

    2 材料與方法

    2.1 樣點布設(shè)

    湖泛是一種因藻類和水草等大量聚積、死亡所引發(fā)的湖泊水體突發(fā)性黑臭現(xiàn)象。典型湖泛區(qū)會出現(xiàn)硫化氫、甲烷等氣體逸出水面, 大面積黑水、底泥上泛, 出現(xiàn)如藻華(algal bloom) 相似開放(bloom)的狀態(tài)。湖泛的生消過程大致可分為4個階段, 即(藻類或水草)生物聚積階段、(水體)耗氧缺氧階段、(湖泛)爆發(fā)成災(zāi)階段, 以及(水體)復(fù)氧消退階段。根據(jù)太湖站2000年至2015年近16年的年報報道[25-28], 太湖主要發(fā)生湖泛的區(qū)域在: 貢湖灣北岸沿岸帶。根據(jù)常規(guī)監(jiān)測和現(xiàn)場情況, 在此區(qū)域設(shè)立5個采樣點(1#-5#), 見圖1。其中1#(N:120.286025; E:31.419991) 為未發(fā)生湖泛區(qū)域, 設(shè)為對照。其余點位坐標(biāo)分別為: 2#(N:120.285939; E:31.423140)、3#(N:120.246714; E:31.412227)、4#(N:120.282677; E:31.420650)、5#(N: 120.278729; E:31.418672)。

    2.2 樣品采集

    用采水器采集每個點位表層水樣(約湖面以下50 cm), 置于預(yù)先用酸浸泡過的塑料桶中, 運回實驗室, 在12 h內(nèi)完成理化參數(shù)測定。另取15 mL水樣裝入預(yù)先滅過菌的離心管中, 加入經(jīng)濾膜(0.22 μm)過濾后的甲醛溶液, 使甲醛終濃度為2%(V/V)。隨后將樣品置于冰盒中, 帶回實驗室用于測定細(xì)菌的豐度群落結(jié)構(gòu)。

    圖1 采樣點位圖

    2.3 水體理化參數(shù)的測定

    水溫(T)、透明度(SD)、電導(dǎo)率(EC)、DO和pH值在現(xiàn)場用YSI6600多參數(shù)水質(zhì)分析儀測定??偟?TN)、溶解性總氮(DTN)、總磷(TP)、溶解性總磷(DTP)、Chl-a等參數(shù)的測定參考湖泊富營養(yǎng)化調(diào)查規(guī)范[29]。

    2.4 細(xì)菌豐度、藻類生物量的測定

    由于天然水環(huán)境中只有0.001%—3%的微生物可以被分離培養(yǎng)[30], 本研究利用熒光染料DAPI (4¢,6-diamidino-2-phenylindole)與表面熒光顯微鏡結(jié)合, 直接計數(shù)環(huán)境中細(xì)菌的方法(epifluores-cence direct counting method, FDC)[31-32]。浮游植物個體體積小, 直接稱重較困難, 且其細(xì)胞比重多接近于1。藻類生物量一般按體積來換算, 可用形態(tài)相近似的幾何體積公式計算細(xì)胞體積, 細(xì)胞體積的毫升數(shù)相當(dāng)于細(xì)胞重量的克數(shù), 這樣體積值(μm3)可直接換算為重量值(109μm3≈1 mg鮮藻重)。

    2.5 DNA的提取、高通量測序

    取200 mL原水經(jīng)0.22 μm孔徑的聚碳酸酯膜(Millipore, Cork, Ireland)過濾后收集濾膜, 用事先滅過菌的剪刀剪碎, 放入滅菌后的2 mL離心管中, 置于-80 ℃的超低溫冰箱中保存至DNA提取。微生物樣品DNA采用試劑盒FastDNA SPIN Kit for Soil(MP Biomedicals)進(jìn)行提取并于-20 ℃下保存, 高通量測序?qū)嶒炍猩虾E缮Z生物科技有限公司完成。

    3 結(jié)果與討論

    3.1 理化參數(shù)與細(xì)菌豐度、藻類數(shù)量的對比

    湖泛突發(fā)區(qū)域, 溶解氧(DO)、透明度(SD)、pH值與湖泛未發(fā)生區(qū)域1#點相比顯著降低。其中2#-5#區(qū)域DO濃度均低于0.51 mg·L-1, SD在湖泛最嚴(yán)重點位處僅為0.13 m。隨著采樣點位向湖泛中心靠近, Chl-a、電導(dǎo)率也都在不斷增加, 水溫也略微升高。其中, 湖泛最嚴(yán)重點位處(5#), Chl-a含量是未發(fā)生湖泛區(qū)域(1#)的近5倍(圖2)。這是由于湖泛形成過程中, 藻類大量腐敗, 水體透明度急劇下降, Chl-a升高[33]。由于水面以下缺乏光照, 水面下大量藻體細(xì)胞的呼吸作用、新陳代謝以及有機質(zhì)的協(xié)同作用消耗了大量的氧氣, 短時間內(nèi)溶解氧就降至接近于零, 藻類的大量分解使得湖泛區(qū)水溫升高、電導(dǎo)率增加并降低水體pH值[34-36]。

    圖3中可以明顯的看到隨著湖泛的不斷加劇, 由湖泛外圍向中心, N、P等營養(yǎng)鹽的含量不斷增加。在湖泛最為嚴(yán)重的5#點TN含量達(dá)到14.09 mg·L-1, 是太湖正常水體TN含量的7倍左右; TP含量達(dá)到0.51 mg·L-1, 比對照點的TP含量高出0.40 mg·L-1左右, 且NH4+、溶解性總氮(DTN)、PO43-和DTP均高于正常值。這是由于水華藍(lán)藻在快速生長繁殖過程中, 群體內(nèi)積聚了大量的營養(yǎng)鹽, 當(dāng)“湖泛”發(fā)生時, 一方面, 由于大量藍(lán)藻的腐敗使其體內(nèi)大量聚集的 N、P營養(yǎng)鹽被迅速釋放到水體中; 而另一方面, “湖泛”中的厭氧環(huán)境會促進(jìn)NH4+-N、PO43?-P從沉積物表層向上覆水?dāng)U散[37]。因此, 這些氮、磷營養(yǎng)鹽的變化主要是由藻類劇烈分解導(dǎo)致形成水體厭氧環(huán)境及大量釋放有機物造成的[38-39]。

    圖2 野外湖泛現(xiàn)場數(shù)據(jù)

    圖 3 野外湖泛的理化參數(shù)

    由圖4可以看出: 細(xì)菌豐度隨著湖泛程度的加劇也由106cells·L-1上升到107cells·L-1, 增加了一個數(shù)量級。藻類數(shù)量在湖泛區(qū)均明顯高于對照點, 但湖泛區(qū)的藻類數(shù)量并沒有隨著湖泛的加劇呈現(xiàn)不斷上升趨勢, 而是大致相等。

    3.2 Miseq宏基因組測序結(jié)果及分析

    5個樣品共得到40萬個序列, 經(jīng)優(yōu)化分析后獲得了10127個operational taxonomic unit (OTUs), 樣品的覆蓋度達(dá)到了95.7%—98.45%, 根據(jù)數(shù)據(jù)繪制了用來表示樣本豐富度與均勻度的幾種多樣性指數(shù)圖(圖5)。

    群落生態(tài)學(xué)中研究微生物多樣性, 通過單樣品的多樣性分析(Alpha多樣性)可以反映微生物群落的豐度和多樣性, 包括通過一系列統(tǒng)計學(xué)分析指數(shù)估計環(huán)境群落的物種豐度和多樣性[40]。Chao指數(shù)和Ace指數(shù)是計算菌群豐度的指數(shù), Chao是用chao1算法估計樣品中所含OTU數(shù)目的指數(shù), 在生態(tài)學(xué)中常用來估計物種總數(shù), 由Chao[41]最早提出的, Ace也是由Chao提出的生態(tài)學(xué)中估計物種總數(shù)的常用指數(shù)之一[42]。Simpson由Edward Hugh Simpson[43]提出, 在生態(tài)學(xué)中常用來定量描述一個區(qū)域的生物多樣性。Simpson指數(shù)值越大, 說明群落多樣性越低。Shannon與Simpson多樣性指數(shù)都常用于反映alpha多樣性指數(shù), 而Shannon值越大, 則說明群落多樣性越高[44-45]。在圖5的多樣性指數(shù)中可以看出, 同樣是湖泛區(qū)的3#物種豐度低于其他點位, 多樣性也低于其它點位。同樣是湖泛區(qū), 不同區(qū)域之間的細(xì)菌群落結(jié)構(gòu)也存在著差異。

    圖4 野外湖泛各個點細(xì)菌豐度及藻類生物量對比圖

    圖5 各個樣本的多樣性指數(shù)

    高通量數(shù)據(jù)的分析結(jié)果顯示: 5個樣本中的序列主要為變形菌門(3367個OTUs)的α-、β-和γ-亞綱(α-,β-,and γ-subdivision of Proteobacteria), 擬桿菌門(Bacteriodetes)(996個OUTs), 放線菌門(Actinoba-cteria)(434個OTUs), 厚壁菌門(Firmicutes)(1379個OTUs)、OD1和一些未分類的門。除細(xì)菌外, 還有少量的古菌分別屬于泉古菌門、廣古菌門和一些未分類的門。進(jìn)一步分析顯示, 屬于α-proteobacteria變形菌門的主要有短波單胞菌屬()和鞘氨醇單胞菌()等。屬于β-Proteo-bacteria的主要有食酸菌屬()、噬氫菌屬()等一些未分類的屬。屬于γ-Proteo-bacteria變形菌的主要有不動桿菌屬()、假單胞菌屬()等等。屬于放線菌門(Actinobacteria)下的優(yōu)勢種為微桿菌屬()、分支桿菌屬()、類諾卡氏菌屬()、和一些未分類的菌屬。屬于擬桿菌門(Bacteriodetes)的屬主要有擬桿菌()、理研科的未分類的一些屬和黃桿菌科下的一些屬。屬于厚壁菌門(Firmicutes)的主要有庫特氏菌屬()、球菌科的一些屬、微桿菌屬()、動性微菌屬()、乳桿菌屬()、梭菌屬()和醋桿菌屬()。

    3.3 湖泛不同區(qū)域的細(xì)菌優(yōu)勢種的差異分析

    由圖6可見: 1#點位處葉綠素最低, 相反湖泛最嚴(yán)重的5#葉綠素最高, 各個餅圖表示測序結(jié)果在門水平上的百分比。目前, 在湖泊水體中共發(fā)現(xiàn)21個典型的淡水細(xì)菌門類, 其中變形菌門(Proteoba-cteria)、藍(lán)細(xì)菌門(Cyanobacteria)、擬桿菌門(Bacteroi-detes)、放線菌門(Actinobacteria)和疣微菌門(Verru-comicrobia)是最主要的5個門類[46]。在這次的高通量檢測結(jié)果主要有變形菌門(Proteobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)、擬桿菌門(Bacteroidetes)、放線菌門(Actinobacteria)和OD1。各個樣品的群落結(jié)果呈現(xiàn)顯著差異(表1)。

    從表1中可以看出, 變形菌門(Proteobacteria)里的α-變形菌綱(α-proteobacteria), γ-變形菌綱(γ-Pro-teo-bacteria)、放線菌門(Actinobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)在未發(fā)生湖泛的區(qū)域與發(fā)生湖泛的區(qū)域相比OTU數(shù)量在不斷減少, 而β-proteo-bacteria、擬桿菌門(Bacteroidetes)的OTU數(shù)量則在發(fā)生湖泛的區(qū)域呈現(xiàn)增加的趨勢。α-變形菌綱(α-Proteobacteria)所占比例由未發(fā)生湖泛時的4.10%降至湖泛最嚴(yán)重區(qū)域的0.59%, 進(jìn)一步分析數(shù)據(jù)發(fā)現(xiàn)α-變形菌綱(α-Proteobacteria)下的短波單胞菌屬()和鞘氨醇單胞菌屬()明顯減少。其中, 鞘氨醇單胞菌()可以降解芳香化合物亦可將戊糖、己糖及二糖轉(zhuǎn)變成酸。但它屬于專性需氧型微生物, 所以可能是由于湖泛發(fā)生區(qū)域的厭氧環(huán)境特征使得這種菌種數(shù)量減少。

    另外一個不同區(qū)域變化最為顯著的是放線菌門(Actinobacteria), 放線菌可以分解許多有機物, 包括芳香化合物、石蠟、橡膠、纖維素、木質(zhì)等復(fù)雜化合物和一些氰等毒性強的化合物。在自然界物質(zhì)循環(huán)中, 更在污水及有機固體廢物的生物處理中有積極的作用, 還能促使土壤形成團(tuán)粒結(jié)構(gòu)而改善土壤。在湖泛發(fā)生的區(qū)域放線菌的數(shù)量并沒有增加而是顯著的減少了, 未發(fā)生湖泛的對照點1#中放線菌的OTU比例為19.04%, 而在湖泛最為嚴(yán)重的5#點中僅為0.92%, 降低了近20倍, 造成原因可能是大多數(shù)放線菌是好氧的且放線菌最適生長pH一般為微堿性。而隨著湖泛的加重, 水體中溶氧極低且pH值也在減小。例如, 在這次發(fā)現(xiàn)的放線菌門(Actinobacteria)下的微桿菌屬()、桿菌屬()、節(jié)桿菌屬 ()和類諾卡氏菌()基本上都是好氧菌且不抗酸, 都隨著湖泛的發(fā)生加劇呈現(xiàn)減少的變化趨勢。

    圖6 細(xì)菌群落組成與葉綠素的關(guān)系

    表1 各個樣品中主要幾種淡水菌種的所占比例(%)

    圖7、8中也列出了一些在湖泛發(fā)生的不同區(qū)域中變化較為顯著的菌種, 在湖泛發(fā)生的水體中增加的菌群大都是厭氧菌, 例如能夠分解多糖的擬桿菌門(Bacteroidetes)是大多是厭氧菌, 與未發(fā)生湖泛的水體中(1.12%)相比, 湖泛區(qū)水體中擬桿菌門(Bacte-roidetes)的OTU比例都不同程度地大幅增加, 最高的達(dá)到36.32%。厚壁菌門(Firmicutes)整體OTU的數(shù)量在湖泛發(fā)生區(qū)域水體中是下降的, 而屬于它的厭氧菌種韋榮氏球菌()和球菌消化鏈球菌屬()的OTU數(shù)量卻是增加的。其中球菌消化鏈球菌屬()在厭氧條件下能夠產(chǎn)酸, 但較適合偏堿性環(huán)境生長。這也就能解釋了該菌種在湖泛發(fā)生區(qū)域的最外圍湖泛不是特別嚴(yán)重的2#點中OTU數(shù)量最高, 但在湖泛最為嚴(yán)重的中心5#點OTU數(shù)量含量最少。因為隨著湖泛越嚴(yán)重其pH值越低, 相比之下, 還是同樣厭氧環(huán)境但偏堿性的2#點適合其生長。但隨著球菌消化鏈球菌屬()不斷的繁殖并分解產(chǎn)酸, 也直接改變了其所在的環(huán)境。從圖2對湖泛監(jiān)測的現(xiàn)場數(shù)據(jù)的分析中可以看出, 隨著湖泛的加劇, 水體中的電導(dǎo)率等指標(biāo)都是增大的, 在高通量的數(shù)據(jù)分析中也發(fā)現(xiàn), 在湖泛較為嚴(yán)重的5#點出現(xiàn)了一些在高鹽地帶出現(xiàn)的菌種, 例如屬于擬桿菌門(Ba-cteroidetes)的。

    在對各個樣品中主要細(xì)菌種類的分析顯示, 5個樣品共發(fā)現(xiàn)放線菌門(Actinobacteria)的有434種菌屬, 在未發(fā)生湖泛的1#點占了395種, 而湖泛最為嚴(yán)重區(qū)域比例升高的菌種的5#僅發(fā)現(xiàn)72個屬, 同樣OTU數(shù)量在減少的α-變形菌綱(α-Proteoba-cteria), γ-變形菌綱(γ-Proteobacteria)在湖泛最嚴(yán)重的5#種類的多樣性也減少, 而擬桿菌門(Bacteroi-detes)和β-變形菌綱(β-Proteobacteria)在發(fā)生湖泛的區(qū)域水體中的種類的多樣性升高。但是對于主要的細(xì)菌種類的分析, 各個樣品即湖泛發(fā)生時不同區(qū)域都沒有達(dá)到在5個樣品所發(fā)現(xiàn)的主要幾種細(xì)菌種類的總數(shù), 這說明湖泛的發(fā)生不僅僅是減少了一些細(xì)菌種類, 也出現(xiàn)了一些未發(fā)生湖泛區(qū)域水體中在不存在的細(xì)菌種類。

    圖 7 湖泛區(qū)域比例降低的菌種

    圖 8 湖泛區(qū)域比例升高的菌種

    4 小結(jié)

    本研究的主要目的是通過比較湖泛發(fā)生過程不同區(qū)域之間及其與未發(fā)生湖泛區(qū)域間水體中理化指標(biāo)的改變、細(xì)菌群落結(jié)構(gòu)的差異, 來了解湖泛對水體中原核微生物豐度及群落結(jié)構(gòu)的影響。不同區(qū)域可以看作是湖泛發(fā)生的不同嚴(yán)重程度, 湖泛中心即是湖泛最為嚴(yán)重的區(qū)域。各個不同區(qū)域水體中細(xì)菌群落的組成可以看做是湖泛各個階段的群落演替過程。與未發(fā)生湖泛區(qū)域水體相比, 湖泛區(qū)主要有以下特點:

    1) 水質(zhì)明顯惡化, DO、透明度(SD)、pH值顯著降低, 其中DO含量均低于0.51 mg·L-1, 屬于厭氧環(huán)境。水體中TN、TP等營養(yǎng)鹽的含量分別是正常水體含量的7倍和5倍左右。

    2) 細(xì)菌豐度和藻類數(shù)量都明顯升高, 其中, 湖泛中心的細(xì)菌豐度上升到107cells·L-1高出對照水體中細(xì)菌豐度一個數(shù)量級。

    3) 變形菌門的α-、β-和γ-綱(α-,β-, and γ- Proteobacteria)、擬桿菌門(Bacteriodetes), 放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)等幾種太湖主要的淡水菌種的比例都發(fā)生了顯著變化。α-、γ-變形菌綱分別由未發(fā)生湖泛區(qū)的4.10%和13.16%降至0.59%和1.91%; 放線菌門(Actinobacteria)所占比例由未發(fā)生湖泛區(qū)域的19.04%降至0.92%; 而擬桿菌門(Bacteriodetes)所占的比例則由1.12%增加至36.32%。

    4) 在湖泛區(qū)域水體中的厭氧環(huán)境下, 不再適合一些專性好氧菌, 而是出現(xiàn)了一些新的產(chǎn)生酸性代謝產(chǎn)物的厭氧菌,,在湖泛區(qū)的群落組成中發(fā)現(xiàn)。

    [1] 陸桂華, 張建華, 馬倩. 太湖湖泛成因及防控關(guān)鍵技術(shù)[C]//水利部水文局、國際水文計劃(IHP)中國國家委員會、國際水文科學(xué)協(xié)會(IAHS)中國國家委員會、中國水利學(xué)會水文專業(yè)委員會. 中國水文學(xué)術(shù)討論會論文集. 北京: 科學(xué)出版社, 2012: 867–874.

    [2] 張軍毅, 朱冰川. 水華、藍(lán)藻水華、湖泛及其與富營養(yǎng)化關(guān)系的探討[C]//中國環(huán)境科學(xué)學(xué)會. 中國環(huán)境科學(xué)學(xué)會學(xué)術(shù)年會論文集. 北京: 北京航空航天大學(xué)出版社, 2013: 3993–3995.

    [3] 陸桂華, 馬倩. 太湖水域“湖泛”及其成因研究[J]. 水科學(xué)進(jìn)展, 2009, 20(3): 438–442.

    [4] 申秋實, 周麒麟, 邵世光, 等. 太湖草源性“湖泛”水域沉積物營養(yǎng)鹽釋放估算[J]. 湖泊科學(xué), 2014, 26(2): 177–184.

    [5] 申秋實, 邵世光, 王兆德, 等. 風(fēng)浪條件下太湖藻源性湖泛的消退及其水體恢復(fù)進(jìn)程[J]. 科學(xué)通報, 2012, 57(12): 1060–1066.

    [6] 王成林, 張詠, 張寧紅, 等. 太湖藻源性“湖泛”形成機制的氣象因素分析[J]. 環(huán)境科學(xué), 2011, 32(2): 401–408.

    [7] 梅森, 段洪濤, 張玉超, 等. 浮游植物降解過程中的水體光學(xué)吸收特性變化研究[J]. 中國環(huán)境科學(xué), 2014, 34(5): 1258–1267.

    [8] SHEN Qiushi, ZHOU Qilin, SHANG Jingge, et al. Beyond hypoxia: occurrence and characteristics of black blooms due to the decomposition of the submerged plant Pota-mogeton crispus in a shallow lake[J]. Journal of Environmental Sciences, 2014, 26(2): 281–288.

    [9] 劉國鋒, 申秋實, 張雷, 等. 藻源性黑水團(tuán)環(huán)境效應(yīng): 對水沉積物界面氮磷變化的驅(qū)動作用[J]. 中國環(huán)境科學(xué), 2010, 31(12): 2917–2924.

    [10] LI Huabing, XING Peng, WU Qinglong. Characterization of the bacterial community composition in a hypoxic zone induced byblooms in Lake Taihu, China[J]. FEMS Microbiology Ecology, 2012, 79(3): 773–784.

    [11] BODAMER B L, BRIDGEMAN T B. Experimental dead zones: two designs for creating oxygen gradients in aquatic ecological studies[J]. Limnology and Oceanography- Methods, 2014, 12(12): 441–454.

    [12] 邢鵬, 胡萬婷, 吳瑜凡, 等. 淺水湖泊湖泛(黑水團(tuán))中的微生物生態(tài)學(xué)研究進(jìn)展[J]. 湖泊科學(xué), 2015, 27(4): 567– 574.

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

    [14] 柴麗紅, 崔曉龍, 彭謙, 等. 青海兩鹽湖細(xì)菌多樣性研究[J]. 微生物學(xué)報, 2004, 44 (3): 271-275.

    [15] VREDE K. Nutrient and temperature limitation of bacterioplankton growth in temperate lakes[J]. Microbial ecology, 2005, 49 (2): 245–256.

    [16] 李玉華, 許其功, 趙越. 松花湖水體中不同空間分布的細(xì)菌群落結(jié)構(gòu)分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2013, 32(4): 764–770.

    [17] 沈燕. 基于FISH方法檢測太湖水體中細(xì)菌優(yōu)勢類群的研究[D]. 南京: 南京農(nóng)業(yè)大學(xué), 2009: 16– 21.

    [18] WU Qiuliang, ZWART G, SCHAUER M, et al. Bacterio-plankton community composition along a salinity gradient of sixteen high-mountain lakes located on the Tibetan Plateau, China[J]. Applied and Environmental Microbiology, 2006, 72 (8): 5478–5485.

    [19] 王興春, 楊致榮, 王敏, 等. 高通量測序技術(shù)及其應(yīng)用[J]. 中國生物工程雜志, 2012, 32 (1): 109-114.

    [20] 夏圍圍, 賈仲君. 高通量測序和DGGE分析土壤微生物群落的技術(shù)評價[J]. 微生物學(xué)報, 2014, 54(12): 1489– 1499.

    [21] 尹承苗, 王功帥, 李園園, 等. 連作蘋果園土壤真菌的T-RFLP分析[J]. 生態(tài)學(xué)報, 2014, 34(4): 837–846.

    [22] 孫宇晶, 劉暢, 龍彥, 等. ADx-ARMS法與PCR-Sanger測序法檢測非小細(xì)胞肺癌微小標(biāo)本EGFR基因突變的比較[J]. 臨床檢驗雜志, 2014, 1(7): 485–489.

    [23] TANG Xiangming, XIE Guijuan, SHAO Keqiang, et al. Influence of salinity on the bacterial community compo-sition in Lake Bosten, a large oligosaline lake in arid northwestern China[J]. Applied and Environmental Microbi-ology, 2012, 78 (13): 4748–4751.

    [24] 李靖宇, 張琇. 騰格里沙漠不同生物土壤結(jié)皮微生物多樣性分析[J]. 生態(tài)科學(xué), 2017, 36(3): 36–42.

    [25] 陸桂華, 馬倩. 2009年太湖水域“湖泛”監(jiān)測與分析[J]. 湖泊科學(xué), 2010, 22(4): 481–487.

    [26] 王成林, 張寧紅, 張詠, 等. 基于氣象條件的太湖湖泛預(yù)警研究[J]. 環(huán)境監(jiān)控與預(yù)警, 2010, 2(5): 1–4.

    [27] 曾誠. 太湖藻源性“湖泛”發(fā)生的過程以及環(huán)境影響研究[D]. 武漢: 華中農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 2015: 23–26.

    [28] 李旭文, 牛志春, 姜晟, 等. 環(huán)境衛(wèi)星CCD影像在太湖湖泛暗色水團(tuán)監(jiān)測中的應(yīng)用[J]. 環(huán)境監(jiān)控與預(yù)警, 2012, 4(3): 1–9.

    [29] 金相燦, 劉鴻亮, 屠清瑛, 等. 中國湖泊富營養(yǎng)化[M]. 北京: 中國環(huán)境科學(xué)出版社, 1990: 64.

    [30] AMANN R I. LUDWIG W, SCLEIFER K H. Phylogenetic identification and in situ detection of individual microbial cells without cultivation[J]. Microbiological Reviews, 1995, 59 (1): 143–169.

    [31] PORTER KG, FEIG Y S. The use of DAPI for identifying and counting aquatic microflora[J]. Limnology and Oceanography, 1980, 25 (5): 943–948.

    [32] 馮勝, 高光, 朱廣偉, 等. 基于16S rDNA-DGGE和FDC技術(shù)對富營養(yǎng)化湖泊不同生態(tài)修復(fù)工程區(qū)細(xì)菌群落結(jié)構(gòu)的研究[J]. 應(yīng)用與環(huán)境生物學(xué)報, 2007, 13 (4): 35–540.

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

    [34] 陳荷生. 太湖宜興近岸水域“湖泛”現(xiàn)象初析[J]. 水利水電科技進(jìn)展, 2011, 31(4): 33–37

    [35] 齊敏, 孫小雪, 鄧緒偉, 等. 太湖不同形態(tài)異味物質(zhì)含量, 相互關(guān)系及其與環(huán)境因子關(guān)系的探討[J]. 湖泊科學(xué), 2013, 25(1): 31–38.

    [36] 尚麗霞, 柯凡, 李文朝, 等. 高密度藍(lán)藻厭氧分解過程與污染物釋放實驗研究[J]. 湖泊科學(xué), 2013, 25(1): 47–54.

    [37] XIE Liqiang, XIE Ping, LI Sixin, et al. The low TN: TP ratio, a cause or a result ofblooms[J]. Water Research, 2003, 37(9): 2073–2080.

    [38] 阮曉紅, 張瑛, 黃林楠, 等. 微生物在濕地氮循環(huán)系統(tǒng)的效應(yīng)分析[J]. 水資源保護(hù), 2005, 20(6): 1-7.

    [39] 郭建華, 彭永臻. 異養(yǎng)硝化, 厭氧氨氧化及古菌氨氧化與新的氮循環(huán)[J]. 環(huán)境科學(xué)學(xué)報, 2008, 28(8): 1489– 1498.

    [40] 李祎, 鄭偉, 鄭天凌. 海洋微生物多樣性及其分子生態(tài)學(xué)研究進(jìn)展[J]. 微生物學(xué)通報, 2013, 40(4): 655– 668.

    [41] CHAO A. Nonparametric estimation of the number of classes in a population[J]. Scandinavian Journal of statistics, 1984, 11(4): 265–270.

    [42] ASMA Z, SYLVIE C, LAURENT C, et al. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil suppleme-ntation of diets[J]. FEMS Microbiology Ecology, 2013, 83(2): 504–514.

    [43] 趙三軍, 肖天, 李洪波, 等. 膠州灣異養(yǎng)細(xì)菌及大腸菌群的分布及對陸源污染的指示[J]. 海洋與湖沼, 2005, 36 (6): 541–547.

    [44] GUO Wei, LI Ying, WANG Lizhi, et al. Evaluation of composition and individual variability of rumen microbiota in yaks by 16S rRNA high-throughput sequencing techno-logy[J]. Anaerobe, 2015, 34: 74–79.

    [45] PITTA D W, PINCHAK W E, DOWD S E, et al. Rumen bacterial diversity dynamics associated w ith changing from Bermudagrass hay to grazed winter wheat diets[J]. Microbial Ecology, 2010, 59(3): 511–522.

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

    馮勝, 袁斌, 王博雯. 湖泛發(fā)生過程中水體理化性質(zhì)及細(xì)菌群落的變化[J]. 生態(tài)科學(xué), 2017, 36(6): 25-34.

    FENGsheng, YUANbin, WANG bowen. Changes in physical and chemical properties of water and bacterial community when black bloom occurs[J]. Ecological Science, 2017, 36(6): 24-34.

    Changes in physical and chemical properties of water and bacterial community when black bloom occurs

    FENGsheng1,*, YUANbin1, WANG bowen2

    1. School of Environmental & safety Engineering, Changzhou University, Changzhou, Jiangsu 213164, China 2. State Key Laboratory of Lake Science and Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing, JiangSu 210008, China

    Molecular biology techniques and multivariate statistical analysis method were used to evaluate relationships between bacterial abundance and environmental variables in the areas when black bloom happened regularly. We also investigated the changes of spatial bacterial community during black bloom happened. Results showed that: DO(dissolved oxygen), SD(secchi disk) and pH(potential of hydrogen) decreased significantly in the black bloom field, particularly, DO decreased as low as 0.51 mg·L-1, which indicated an anaerobic environment in the water column. In addition, black bloom increased nutrient concentrations (total nitrogen, total phosphorus), phytoplankton biomass and bacterial abundance. The bacterial abundance in the central of black booms area exceeded 107cells·L-1, which improved an order of magnitude higher than other areas without black bloom; bacterial community structures were significantly different among the areas: alpha-, beta-, and gamma-Proteobacteria, Bacteriodetes, Actinobacteria. Namely, along black bloom, α-, γ-proteobacteriadecreased from 4.10% and 13.16% to 0.59% and 1.91% respectively; anddeceased from 19.04% to 0.92%. In contrast, Bacteriodetes increased from 1.12% to 36.32%. In addition, some special habitats bacteria were also detected in our samples likeand. Through studying the change rules of bacterial community and key environmental factors when black bloom occurs, we can provide some scientific instruction for further research of black bloom.

    aquatic environmental factor; bacterial abundance; bacterial community structure; black bloom

    10.14108/j.cnki.1008-8873.2017.06.004

    X524

    A

    1008-8873(2017)06-025-10

    2016-07-19;

    2016-09-15

    國家自然科學(xué)基金(41371446); 國家自然科學(xué)基金(41271498)

    馮勝(1974—), 男, 江蘇常州人, 博士, 副教授, 主要從事湖泊治理研究, E-mail: shfeng@cczu.edu.cn

    E-mail: shfeng@cczu.edu.cn.

    猜你喜歡
    放線菌菌門太湖
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    肺曲霉菌合并肺放線菌感染一例
    饑餓與重攝食對河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價值分析
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    太湖一角
    南大西洋深海沉積物中可培養(yǎng)放線菌的多樣性
    91久久精品国产一区二区成人| 26uuu在线亚洲综合色| 国产精品成人在线| 国产国语露脸激情在线看| 日韩精品有码人妻一区| 成人国语在线视频| 亚洲精品456在线播放app| 国产精品无大码| 亚洲在久久综合| 免费黄频网站在线观看国产| 精品国产一区二区三区久久久樱花| 亚洲第一av免费看| 国产成人a∨麻豆精品| 各种免费的搞黄视频| 高清视频免费观看一区二区| 大片免费播放器 马上看| av播播在线观看一区| 在线观看免费日韩欧美大片 | 夜夜爽夜夜爽视频| 大香蕉97超碰在线| av天堂久久9| 校园人妻丝袜中文字幕| 成人漫画全彩无遮挡| 午夜免费鲁丝| 蜜臀久久99精品久久宅男| 天天躁夜夜躁狠狠久久av| av播播在线观看一区| 中文字幕人妻丝袜制服| 久久99一区二区三区| 边亲边吃奶的免费视频| 街头女战士在线观看网站| 97超碰精品成人国产| 国产男女超爽视频在线观看| 女性生殖器流出的白浆| 免费观看av网站的网址| 秋霞伦理黄片| 久久久久久久亚洲中文字幕| 在线观看免费视频网站a站| 成人影院久久| 国产成人av激情在线播放 | 国产国语露脸激情在线看| 波野结衣二区三区在线| 视频在线观看一区二区三区| 亚洲国产精品999| 久久人妻熟女aⅴ| 亚洲精品美女久久av网站| 国产乱来视频区| 久热久热在线精品观看| 在线 av 中文字幕| 男人操女人黄网站| 日韩精品有码人妻一区| 乱码一卡2卡4卡精品| 一区二区日韩欧美中文字幕 | av.在线天堂| 高清黄色对白视频在线免费看| 一区二区三区精品91| 国产毛片在线视频| 久久久久国产网址| 国产精品女同一区二区软件| 婷婷色麻豆天堂久久| 黄片无遮挡物在线观看| 18+在线观看网站| 国产成人精品福利久久| 亚洲中文av在线| 久久久午夜欧美精品| 国产黄频视频在线观看| 哪个播放器可以免费观看大片| 国产色爽女视频免费观看| 亚洲无线观看免费| 午夜影院在线不卡| 不卡视频在线观看欧美| 国产淫语在线视频| 亚洲精品456在线播放app| 久久久精品94久久精品| 国产av国产精品国产| 亚洲精品国产av成人精品| 婷婷色综合www| 水蜜桃什么品种好| 黄色视频在线播放观看不卡| 狠狠精品人妻久久久久久综合| 国产爽快片一区二区三区| 国产日韩欧美亚洲二区| 夜夜爽夜夜爽视频| 久久韩国三级中文字幕| 国产 精品1| 成人二区视频| 日韩亚洲欧美综合| 国产老妇伦熟女老妇高清| 毛片一级片免费看久久久久| 九九久久精品国产亚洲av麻豆| 免费看光身美女| 国产av精品麻豆| .国产精品久久| 三级国产精品片| 飞空精品影院首页| 热re99久久国产66热| 亚洲精品国产色婷婷电影| 日产精品乱码卡一卡2卡三| 国产精品成人在线| 亚洲国产av影院在线观看| 免费观看在线日韩| 亚洲国产日韩一区二区| 精品人妻一区二区三区麻豆| 日韩成人伦理影院| 久久精品国产亚洲av涩爱| 精品人妻熟女毛片av久久网站| 晚上一个人看的免费电影| 日本欧美视频一区| 国产精品一区二区在线观看99| 国产黄频视频在线观看| 色吧在线观看| av国产精品久久久久影院| 最近2019中文字幕mv第一页| 韩国av在线不卡| 午夜激情av网站| 亚洲欧美日韩另类电影网站| 亚洲av成人精品一二三区| 久久久精品94久久精品| 国产男人的电影天堂91| 制服丝袜香蕉在线| 国产精品成人在线| 日韩一本色道免费dvd| 久久亚洲国产成人精品v| 久久女婷五月综合色啪小说| 乱人伦中国视频| 精品卡一卡二卡四卡免费| 亚洲精品日本国产第一区| 丰满少妇做爰视频| 一本色道久久久久久精品综合| 精品一区二区免费观看| 一级a做视频免费观看| 国产又色又爽无遮挡免| 日本黄色日本黄色录像| 大陆偷拍与自拍| 特大巨黑吊av在线直播| av专区在线播放| 亚洲精华国产精华液的使用体验| 久久久久久久久久久久大奶| 中文天堂在线官网| 免费av不卡在线播放| 欧美精品一区二区免费开放| 又大又黄又爽视频免费| av又黄又爽大尺度在线免费看| 黑人欧美特级aaaaaa片| 2021少妇久久久久久久久久久| 特大巨黑吊av在线直播| 欧美精品高潮呻吟av久久| 亚洲成色77777| 我的女老师完整版在线观看| 午夜福利影视在线免费观看| 一个人免费看片子| 考比视频在线观看| 成年美女黄网站色视频大全免费 | 久久久久久久亚洲中文字幕| 男女边摸边吃奶| 国产黄片视频在线免费观看| 久久精品国产鲁丝片午夜精品| 免费观看av网站的网址| 97在线视频观看| 久久99热6这里只有精品| 亚洲成色77777| av又黄又爽大尺度在线免费看| 亚洲人与动物交配视频| 国产探花极品一区二区| 九九久久精品国产亚洲av麻豆| 国内精品宾馆在线| 久久久欧美国产精品| 成人手机av| 亚洲精品av麻豆狂野| 免费黄频网站在线观看国产| 男女啪啪激烈高潮av片| 久久久久久久精品精品| 新久久久久国产一级毛片| 国产成人精品无人区| 久久影院123| 国产在线一区二区三区精| 日本黄大片高清| 精品一区在线观看国产| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 亚洲精品成人av观看孕妇| 国产探花极品一区二区| 女人精品久久久久毛片| 水蜜桃什么品种好| 国产成人freesex在线| 九九久久精品国产亚洲av麻豆| 中文乱码字字幕精品一区二区三区| 黄色欧美视频在线观看| 大香蕉久久成人网| 亚洲婷婷狠狠爱综合网| 亚洲丝袜综合中文字幕| 欧美精品一区二区免费开放| 成人无遮挡网站| 午夜免费观看性视频| 国产黄频视频在线观看| 精品卡一卡二卡四卡免费| 男女无遮挡免费网站观看| 免费人妻精品一区二区三区视频| 91精品伊人久久大香线蕉| 国产色爽女视频免费观看| 亚洲一区二区三区欧美精品| 高清不卡的av网站| 国产成人精品一,二区| 最新的欧美精品一区二区| 亚洲欧美日韩卡通动漫| 亚洲国产毛片av蜜桃av| 久久久久久久久久久久大奶| 久久青草综合色| 精品亚洲成国产av| 这个男人来自地球电影免费观看 | 全区人妻精品视频| 男女边吃奶边做爰视频| 国产淫语在线视频| 日产精品乱码卡一卡2卡三| 男男h啪啪无遮挡| 午夜影院在线不卡| 一区二区三区精品91| 好男人视频免费观看在线| av专区在线播放| 婷婷色综合大香蕉| 婷婷色综合www| 熟女av电影| 亚洲欧美成人综合另类久久久| 一区二区三区乱码不卡18| 自线自在国产av| 国产在线免费精品| 久久人人爽人人爽人人片va| 啦啦啦啦在线视频资源| 人人妻人人爽人人添夜夜欢视频| 男的添女的下面高潮视频| 中文字幕最新亚洲高清| 亚洲精品第二区| 在线观看一区二区三区激情| 男女啪啪激烈高潮av片| 美女福利国产在线| 男女无遮挡免费网站观看| 91aial.com中文字幕在线观看| 久久午夜福利片| 成人午夜精彩视频在线观看| 最黄视频免费看| 日韩亚洲欧美综合| 国产一区二区三区综合在线观看 | av不卡在线播放| 蜜桃在线观看..| 久久久久久久久大av| 韩国av在线不卡| 亚洲成人av在线免费| 中国国产av一级| 老女人水多毛片| 国产黄频视频在线观看| 老司机亚洲免费影院| 日韩,欧美,国产一区二区三区| 亚洲综合色网址| 亚洲国产成人一精品久久久| 天天躁夜夜躁狠狠久久av| 午夜老司机福利剧场| 国产毛片在线视频| 国产男女内射视频| 国产在线一区二区三区精| 蜜臀久久99精品久久宅男| 国产免费一区二区三区四区乱码| 成人午夜精彩视频在线观看| 国产一区二区在线观看av| 在线免费观看不下载黄p国产| 久久久久久久久久人人人人人人| 制服人妻中文乱码| 中文字幕久久专区| 国产 一区精品| 国产精品久久久久久精品电影小说| 国产精品一区二区三区四区免费观看| 一级a做视频免费观看| 蜜桃久久精品国产亚洲av| 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 赤兔流量卡办理| 十分钟在线观看高清视频www| 国产精品一二三区在线看| 免费观看在线日韩| 亚洲欧美日韩卡通动漫| 两个人免费观看高清视频| 日韩精品有码人妻一区| 亚洲五月色婷婷综合| 免费观看av网站的网址| 亚洲精品av麻豆狂野| 国产精品国产av在线观看| 成人国语在线视频| 国产精品久久久久久久久免| 天天躁夜夜躁狠狠久久av| av国产久精品久网站免费入址| 亚洲精华国产精华液的使用体验| 成人国语在线视频| 午夜精品国产一区二区电影| 亚洲av国产av综合av卡| 国产精品久久久久成人av| 久久久久国产网址| 男女高潮啪啪啪动态图| 精品亚洲成a人片在线观看| 中文欧美无线码| 美女国产高潮福利片在线看| 亚洲人成77777在线视频| 久久午夜综合久久蜜桃| 免费少妇av软件| 哪个播放器可以免费观看大片| a 毛片基地| 波野结衣二区三区在线| 国产综合精华液| 欧美97在线视频| 日本wwww免费看| 丝袜脚勾引网站| 99精国产麻豆久久婷婷| 啦啦啦中文免费视频观看日本| 99久久精品一区二区三区| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 日韩人妻高清精品专区| 一本大道久久a久久精品| 国产在线免费精品| 成人漫画全彩无遮挡| 中文字幕久久专区| 一级毛片 在线播放| 久热久热在线精品观看| 黄片无遮挡物在线观看| 97精品久久久久久久久久精品| 国产成人精品久久久久久| 最后的刺客免费高清国语| 欧美日本中文国产一区发布| 亚洲欧洲精品一区二区精品久久久 | 少妇的逼水好多| 亚洲av福利一区| 极品少妇高潮喷水抽搐| 26uuu在线亚洲综合色| 欧美亚洲日本最大视频资源| 亚洲激情五月婷婷啪啪| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 91精品伊人久久大香线蕉| 亚洲,欧美,日韩| 大片电影免费在线观看免费| 国产日韩欧美亚洲二区| 亚洲国产精品专区欧美| 熟妇人妻不卡中文字幕| 五月伊人婷婷丁香| 欧美丝袜亚洲另类| 日本与韩国留学比较| 日韩 亚洲 欧美在线| 寂寞人妻少妇视频99o| 久久午夜福利片| 精品酒店卫生间| 久久97久久精品| 国产精品99久久久久久久久| 中文字幕最新亚洲高清| 人妻 亚洲 视频| 一级片'在线观看视频| 一个人免费看片子| 国产成人91sexporn| 99热这里只有精品一区| 黄片无遮挡物在线观看| 一级黄片播放器| 国产精品 国内视频| 精品一区二区三区视频在线| 纵有疾风起免费观看全集完整版| 久久国内精品自在自线图片| 高清午夜精品一区二区三区| 99久久综合免费| www.av在线官网国产| 黑人欧美特级aaaaaa片| 亚洲精华国产精华液的使用体验| 多毛熟女@视频| 中文字幕久久专区| av又黄又爽大尺度在线免费看| av黄色大香蕉| 最新中文字幕久久久久| 久久精品久久久久久久性| 国产黄色免费在线视频| 欧美xxxx性猛交bbbb| 边亲边吃奶的免费视频| 美女国产高潮福利片在线看| 伊人亚洲综合成人网| 色吧在线观看| 久久精品国产鲁丝片午夜精品| 精品久久久久久电影网| 精品亚洲成a人片在线观看| 亚洲av.av天堂| 久久久久久久久久久丰满| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 国产精品成人在线| 欧美丝袜亚洲另类| 91aial.com中文字幕在线观看| 国产精品久久久久久精品古装| 亚洲精品,欧美精品| 99热6这里只有精品| 伊人亚洲综合成人网| 简卡轻食公司| 亚洲av欧美aⅴ国产| 日本欧美国产在线视频| 亚洲欧美一区二区三区黑人 | 免费观看性生交大片5| 国产成人一区二区在线| 国产高清国产精品国产三级| 国产又色又爽无遮挡免| 制服人妻中文乱码| 欧美激情 高清一区二区三区| 国产亚洲最大av| 亚洲人成网站在线观看播放| 黑人高潮一二区| 国产日韩欧美在线精品| 色哟哟·www| 一级毛片aaaaaa免费看小| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 国产男人的电影天堂91| 青春草视频在线免费观看| 欧美精品国产亚洲| 日本欧美视频一区| 中文精品一卡2卡3卡4更新| 亚洲欧美色中文字幕在线| 亚洲精品乱久久久久久| 天天影视国产精品| 国产国拍精品亚洲av在线观看| 午夜视频国产福利| 久久久精品94久久精品| 熟女av电影| 美女xxoo啪啪120秒动态图| 亚洲欧美成人精品一区二区| 18禁在线播放成人免费| 一个人看视频在线观看www免费| 黄片播放在线免费| 成年人免费黄色播放视频| 久久久久久久精品精品| 欧美激情极品国产一区二区三区 | 国产深夜福利视频在线观看| 国产一区二区在线观看av| 极品少妇高潮喷水抽搐| 日韩三级伦理在线观看| 男人添女人高潮全过程视频| 有码 亚洲区| 国产精品人妻久久久久久| 天美传媒精品一区二区| 精品酒店卫生间| 一区二区三区免费毛片| 国产av码专区亚洲av| 日韩免费高清中文字幕av| 欧美三级亚洲精品| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 建设人人有责人人尽责人人享有的| 街头女战士在线观看网站| 丝袜在线中文字幕| 又大又黄又爽视频免费| 在线观看人妻少妇| 一级毛片电影观看| 高清视频免费观看一区二区| 国产免费视频播放在线视频| 夫妻性生交免费视频一级片| 国产在视频线精品| 永久免费av网站大全| 精品人妻熟女av久视频| 人妻系列 视频| 国产无遮挡羞羞视频在线观看| 亚洲三级黄色毛片| 成人国产麻豆网| 国产精品蜜桃在线观看| 国产极品天堂在线| 国产精品一区www在线观看| 永久网站在线| 人妻人人澡人人爽人人| 久久精品人人爽人人爽视色| 国产av国产精品国产| 麻豆成人av视频| 欧美日韩视频高清一区二区三区二| 97超碰精品成人国产| 久久久精品94久久精品| 日韩一区二区视频免费看| 亚洲一级一片aⅴ在线观看| 飞空精品影院首页| 久久99精品国语久久久| 在线免费观看不下载黄p国产| 这个男人来自地球电影免费观看 | 香蕉精品网在线| 日本vs欧美在线观看视频| 亚洲精品一区蜜桃| 国产伦精品一区二区三区视频9| 亚洲av不卡在线观看| 日韩制服骚丝袜av| 日本91视频免费播放| 九九久久精品国产亚洲av麻豆| 在线亚洲精品国产二区图片欧美 | 久久亚洲国产成人精品v| 亚洲欧美日韩卡通动漫| 国产熟女欧美一区二区| 人妻系列 视频| av福利片在线| 高清在线视频一区二区三区| 夫妻午夜视频| 国产成人精品福利久久| 热99国产精品久久久久久7| 只有这里有精品99| 国产在线一区二区三区精| 自拍欧美九色日韩亚洲蝌蚪91| 人妻一区二区av| 精品国产国语对白av| 午夜福利网站1000一区二区三区| 另类精品久久| 亚州av有码| 久久久欧美国产精品| 国产av国产精品国产| 中文字幕人妻丝袜制服| 久久99热6这里只有精品| 亚洲欧美成人综合另类久久久| 大又大粗又爽又黄少妇毛片口| 黑人欧美特级aaaaaa片| 日本免费在线观看一区| 亚洲精品成人av观看孕妇| 最新中文字幕久久久久| av专区在线播放| xxxhd国产人妻xxx| 久久99一区二区三区| 妹子高潮喷水视频| 久久久久视频综合| 久久亚洲国产成人精品v| 丰满少妇做爰视频| 国产黄色免费在线视频| 97在线人人人人妻| av天堂久久9| 国产高清不卡午夜福利| 久久精品夜色国产| 草草在线视频免费看| 在线观看免费视频网站a站| 国精品久久久久久国模美| 欧美xxⅹ黑人| 91精品国产九色| 久久久久视频综合| 伊人久久精品亚洲午夜| 免费日韩欧美在线观看| 成人黄色视频免费在线看| 久久精品国产a三级三级三级| 男女国产视频网站| 色94色欧美一区二区| 少妇 在线观看| 看非洲黑人一级黄片| 免费人成在线观看视频色| 综合色丁香网| 高清黄色对白视频在线免费看| 国产精品久久久久久久久免| 黑人高潮一二区| 久久精品国产亚洲av涩爱| 99久久精品国产国产毛片| 街头女战士在线观看网站| 日本91视频免费播放| 久久综合国产亚洲精品| 中文字幕人妻丝袜制服| 国产一级毛片在线| 日韩av免费高清视频| 久久久久国产网址| 午夜激情av网站| 久久精品国产自在天天线| 日韩一区二区视频免费看| 丝袜美足系列| 天天躁夜夜躁狠狠久久av| 亚洲精品乱码久久久v下载方式| 国产亚洲午夜精品一区二区久久| .国产精品久久| 亚洲精品第二区| 女人精品久久久久毛片| 99国产精品免费福利视频| 99久久精品国产国产毛片| av免费观看日本| 男人添女人高潮全过程视频| 在线观看免费高清a一片| 国产精品 国内视频| 亚洲国产精品国产精品| 国产亚洲最大av| 国产精品久久久久久久电影| 热re99久久精品国产66热6| 自拍欧美九色日韩亚洲蝌蚪91| 日韩,欧美,国产一区二区三区| 午夜激情av网站| 一区二区av电影网| 免费看不卡的av| videossex国产| 国产午夜精品久久久久久一区二区三区| 99热这里只有精品一区| 男女高潮啪啪啪动态图| 高清午夜精品一区二区三区| 热re99久久精品国产66热6| 久久人人爽人人爽人人片va| 亚洲成色77777| 国产日韩欧美视频二区| av又黄又爽大尺度在线免费看| 欧美日本中文国产一区发布| 简卡轻食公司| 日日撸夜夜添| 99国产精品免费福利视频| 亚洲精品日本国产第一区| 69精品国产乱码久久久| 久热这里只有精品99| 美女福利国产在线| 男女边摸边吃奶| 丁香六月天网| av在线播放精品| 日韩一本色道免费dvd| 国产日韩欧美亚洲二区| 在线观看美女被高潮喷水网站| 欧美精品一区二区大全| 美女视频免费永久观看网站| 免费久久久久久久精品成人欧美视频 | av在线app专区| 天堂8中文在线网| 天天操日日干夜夜撸| 亚洲精品一区蜜桃| 黄色怎么调成土黄色| 母亲3免费完整高清在线观看 | av电影中文网址|