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

    夏季達里諾爾湖浮游細菌群落表、底層結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子

    2022-01-13 01:44:38杜蕾李暢游李文寶史小紅楊旭劉晶晶
    生態(tài)科學(xué) 2021年6期
    關(guān)鍵詞:浮游菌門底層

    杜蕾, 李暢游, 李文寶, 史小紅, 楊旭, 劉晶晶

    夏季達里諾爾湖浮游細菌群落表、底層結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子

    杜蕾, 李暢游*, 李文寶, 史小紅, 楊旭, 劉晶晶

    內(nèi)蒙古農(nóng)業(yè)大學(xué)水資源保護與利用自治區(qū)重點實驗室, 呼和浩特 010018

    浮游細菌是湖泊水體物質(zhì)循環(huán)和能量流動過程中的關(guān)鍵要素, 影響著整個生態(tài)系統(tǒng)的遷移轉(zhuǎn)化。為探索中國北方典型內(nèi)陸封閉型湖泊浮游細菌群落結(jié)構(gòu)對環(huán)境因子變動的響應(yīng), 以夏季內(nèi)蒙古達里諾爾湖(以下簡稱“達里湖”)為例, 采集表、底層水共32個樣品, 基于16S rRNA基因測序, 詳細分析了浮游細菌群落結(jié)構(gòu)的表、底層變化特征及其主要影響因素, 結(jié)果顯示: (1)夏季達里湖水體中浮游細菌群落多樣性和均勻度均呈現(xiàn)表層>底層, 豐富度則呈底層>表層。(2)浮游細菌群落以變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍細菌門(Cyanobacteria)和擬桿菌門(Bacteroidetes)四大優(yōu)勢菌門為主要組成部分, 且變形菌門含量在表、底層水體中均占絕對優(yōu)勢: 在底層水中占比46.46 %, 在表層水中占比36.21 %。(3)基于不同環(huán)境因子和浮游細菌優(yōu)勢菌門間的Pearson相關(guān)性分析, 表層水體中總?cè)芙庑怨腆w(TDS)、溶解態(tài)無機磷(DIP)是多數(shù)優(yōu)勢菌門的主要驅(qū)動因子, 而底層水體中的總氮(TN)、葉綠素a(Chl)和水溫(WT)為主要驅(qū)動因子??偟膩碚f, 達里諾爾湖表、底層浮游細菌群落結(jié)構(gòu)與不同環(huán)境因子間存在顯著相關(guān), 而水體環(huán)境因子的表、底層差異是影響多數(shù)優(yōu)勢菌門演化的重要因素。

    浮游細菌; 表、底層變化; 夏季; 群落結(jié)構(gòu); 達里諾爾湖

    0 前言

    湖泊在區(qū)域生態(tài)系統(tǒng)中起著至關(guān)重要的作用, 而湖泊水體中的微生物則是這一生態(tài)系統(tǒng)中的關(guān)鍵組成部分[1-2]。浮游細菌作為主要分解者[3-4], 參與了湖泊水體中的各種物質(zhì)循環(huán)和能量流動過程[5-8]。此外, 水體中的浮游細菌群落結(jié)構(gòu)變化與其水質(zhì)也有著密切的聯(lián)系, 且對水體環(huán)境變化響應(yīng)快速: 不同形態(tài)營養(yǎng)元素的輸入、水體環(huán)境因子的變化都會導(dǎo)致細菌群落結(jié)構(gòu)的改變[9-10]。因此, 分析水體浮游細菌群落結(jié)構(gòu)及優(yōu)勢菌門含量變化可以很好地反映湖泊的水環(huán)境質(zhì)量狀況[11-14]。

    目前, 關(guān)于湖泊水體生態(tài)系統(tǒng)中浮游細菌菌群結(jié)構(gòu)特征及其與環(huán)境因子之間關(guān)系的研究已經(jīng)廣泛展開。如Garcia, De Figueiredo等通過研究季節(jié)變化對湖泊水體環(huán)境的影響, 進而討論了其中細菌群落結(jié)構(gòu)的變化[1,8]; 在國內(nèi), 姜學(xué)霞、呂明姬、沈烽、楊長明等分別對河、湖水體及沉積物中的浮游細菌群落結(jié)構(gòu)特征及影響因素進行了分析[14-17]。結(jié)果顯示不同環(huán)境因子, 如水體溫度(WT)、溶解氧含量(DO)、溶解性總固體含量(TDS)、營養(yǎng)元素(如總磷-TP、總氮-TN等)等, 對浮游細菌的群落結(jié)構(gòu)影響明顯[17,18]。然而, 相對于這些研究區(qū)域, 中國內(nèi)蒙古高原內(nèi)陸封閉型湖泊不僅水文循環(huán)過程影響因素簡單[19], 而且氣候環(huán)境變化存在明顯的冬季冰封期(相對封閉環(huán)境)-夏季非冰封期(相對開放環(huán)境)交替的區(qū)域特點[20, 21]。因此, 在具有重要區(qū)域生態(tài)環(huán)境意義的內(nèi)陸湖泊開展菌群結(jié)構(gòu)變化研究, 既可以為深入探討區(qū)域生態(tài)環(huán)境演變提供理論基礎(chǔ), 又可以為不同氣候條件下的湖泊水體生態(tài)系統(tǒng)演變對比提供基礎(chǔ)材料。

    基于此, 本文將以內(nèi)蒙古赤峰市克什克騰旗境內(nèi)達里諾爾湖(簡稱“達里湖”)為研究區(qū)域, 通過對2018年夏季湖水進行垂向采樣, 分析不同空間部分環(huán)境因子參數(shù)(如WT、DO、TDS、水深、不同形態(tài)營養(yǎng)元素含量等)的變化特征; 同時, 利用16S rRNA基因測序?qū)Σ煌顦悠分械母∮渭毦郝浣Y(jié)構(gòu)進行測試, 分析浮游細菌群落結(jié)構(gòu)組成及垂向變化特征; 進一步結(jié)合Pearson相關(guān)性分析方法研究浮游細菌優(yōu)勢菌門和不同環(huán)境因子間的關(guān)系, 討論影響浮游細菌群落組成及演變的關(guān)鍵因素, 為進一步研究達里湖水環(huán)境變化、水生態(tài)演化過程提供理論基礎(chǔ), 以期為達里諾爾湖水體環(huán)境保護提供一定參考。

    1 研究區(qū)概況

    達里湖(43o12′—43o24′ N, 116o24′—116o56′ E)位于內(nèi)蒙古自治區(qū)赤峰市克什克騰旗境內(nèi), 為封閉式蘇達型半咸水湖, 屬高原內(nèi)陸湖, 湖水無外瀉, 補給源為大氣降水, 河水(亮子河、沙里河、浩來河與貢格爾河)及地下水 (圖1), 其中亮子河與貢格爾河為主要補給河流, 達里湖年入湖水量存在差異, 主要由于年均降水量不同[19]。

    達里湖水面面積平均約為189 km2, 年均面積變化在1—2 km2內(nèi)[19]。湖泊的面積變化造成不同區(qū)域水深出現(xiàn)一定程度波動, 具體體現(xiàn)在從西南到東北逐漸變淺, 湖水最深處約13 m[21]。其中, 西南側(cè)水深變化明顯, 形成了顯著的水下陡坡。而其他區(qū)域湖岸均比較平緩, 湖水深度變化也較平緩[22, 23]。

    注: (a) 取樣點分布示意圖。取樣點依具體位置由北到南依次編號。左上小圖代表內(nèi)蒙古自治區(qū)的年均降水量區(qū)域變化特征; (b)達里湖水深變化概況示意圖, 修改自肖舉樂, 等[22, 23]. 圖中白色數(shù)值代表當(dāng)時水深。由于近年達里湖水位變化, 具體湖水深度存在一定程度的波動。

    Figure 1 Distribution of sampling points and water depth profile of Dali-nor Lake

    2 材料與方法

    2.1 樣品采集

    根據(jù)達里湖的水面面積、水深以及入湖河口的位置, 本文采用GPS全球定位系統(tǒng), 在達里湖內(nèi)共布設(shè)17個取樣點(圖1), 于2018年7月進行了湖水樣品采集工作。樣品采集主要劃分為表層、底層2層, 其中DL-9水深較淺, 僅采集表層水樣, 由于實驗操作失誤丟失DL-5底層水樣, 共采集水樣32個。采集完成的水樣裝入經(jīng)滅菌處理的1 L的聚乙烯取樣瓶中, 并立即密封避光保存在低溫保溫箱中盡快運回實驗室進行過濾處理, 其中500 mL用于環(huán)境因子的測定, 500 mL用于收集細菌, 首先水樣先經(jīng)過5 μm孔徑濾膜過濾, 以去除顆粒雜質(zhì)(濾膜可用以測試附著細菌)[24, 25], 然后在經(jīng)過0.22 μm孔徑濾膜過濾[24, 25], 過濾完成后, 濾膜置于凍融管中保存在-20℃條件下, 用于后續(xù)實驗分析。

    2.2 水體環(huán)境因子測定

    利用YSI Pro多參數(shù)水質(zhì)分析儀現(xiàn)場測定DO、WT、TDS、EC, 塞氏圓盤測定SD、測深鉛錘測定水深。同時, TN、TP、溶解態(tài)無機磷(DIP)、溶解態(tài)總磷(DTP)、葉綠素(Chl)等環(huán)境因子的測定則是在內(nèi)蒙古農(nóng)業(yè)大學(xué)水資源利用與保護自治區(qū)重點實驗室完成, 標(biāo)準(zhǔn)參照《水和廢水監(jiān)測分析方法》(第四版)[26]。

    2.3 DNA提取和高通量測序

    按照FastDNA? Spin Kit for Soil試劑盒的具體實驗操作步驟方法, 提取湖水中微生物基因組DNA。隨后利用微量紫外分光光度計(NanoDrop-2000)測定提取出的DNA濃度和純度后, 選用通用引物515F(5`- GTGCCAGCMGCC-GCGGTAA-3`)和806R (5`-GGACTACVSGGG TATCTAAT-3`)[27]對V4-V5 可變區(qū)進行PCR 擴增。

    擴增程序[28]為95℃預(yù)變性3 min, 27 個循環(huán)(95 ℃變性30 s, 55 ℃退火30 s, 72 ℃延伸30 s), 最后72 ℃延伸10 min (PCR 儀: ABI GeneAmp? 9700 型)。擴增體系為20 uL, 4 uL 5*FastPfu 緩沖液, 2 uL 2.5 mM dNTPs, 0.8 uL 引物(5 uM), 0.4 uL FastPfu 聚合酶; 10 ng DNA 模板。每個樣品3 個重復(fù), 將不同樣品的PCR 擴增產(chǎn)物均一化至10 nmol·L-1后等體積混合, 利用上海美吉生物醫(yī)藥科技有限公司的MiSeq PE300測序儀完成序列測定。

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

    高通量測序所得數(shù)據(jù)保存在NCBI數(shù)據(jù)庫內(nèi), 其生物信息學(xué)分析采用 Qiime軟件進行, 根據(jù)序列的相似度, 將序列歸為多個OTU, OTU 產(chǎn)出后, 統(tǒng)計各個樣品中含有OTU的情況及每個OTU中所含序列的數(shù)目, 得到每個OTU 的分類學(xué)信息[29]。選取相似度在97%條件下的OTU生成預(yù)期的稀釋曲線, 并利用Mothur軟件計算豐富度指數(shù)Chao和Ace, 文庫覆蓋度, 均勻度指數(shù)Heip以及多樣性指數(shù)Simpson 和Shannon指數(shù)進行Alpha 多樣性分析[30]。利用非度量多維尺度分析(NMDS)方法, 對不同深度湖水中浮游細菌群落結(jié)構(gòu)的垂向差異進行分析。利用Pearson相關(guān)分析進行了環(huán)境因子與優(yōu)勢菌群間的相關(guān)性分析。

    3 結(jié)果與分析

    3.1 達里湖水體環(huán)境因子表、底層差異變化

    為探究達里湖水體環(huán)境因子表、底層差異, 對各取樣點的理化指標(biāo)計算均值, 并計算各指標(biāo)表、底均值的變異系數(shù)。Chl的變化最大, 變異系數(shù)為10.889%; pH值的變化最小, 變異系數(shù)為0.007%。除Chl外, 各理化指標(biāo)的變異系數(shù)均較小, 可見達里湖表、底層水體交換過程相對緩慢, 水體理化指標(biāo)變化相對穩(wěn)定[31]。

    表1 表、底層水體環(huán)境因子的平均值及變異系數(shù)*

    注: “*”平均值為每層水體所有取樣點的平均值, 黑體字代表平均值最高。

    在夏季開放條件下的達里湖, 表層水更易受太陽輻射、大氣降水及河流輸入等外界因素的直接影響, 因此導(dǎo)致達里湖表層水的溫度、TDS值較高。夏季適宜的水體溫度導(dǎo)致藻類大量繁殖, 造成表層水體的Chl值高于底層, 而底層水體中TN、DIP、COD環(huán)境因子參數(shù)的平均值高于表層。

    3.2 達里湖浮游細菌多樣性分析

    表、底層水體樣品測序結(jié)果的文庫覆蓋率均高于99.7%(表2), 測序條數(shù)達到10000條以上時稀釋性曲線趨向平穩(wěn)(圖2), 表明測序數(shù)據(jù)量合理, 能夠代表物種的豐富度。采用群落豐富度指數(shù)Chao和Ace、群落多樣性指數(shù)Shannon 和Simpson、均勻度指數(shù)Heip和文庫覆蓋率等對達里湖浮游細菌群落進行評估, 結(jié)果表明樣品具有豐富的群落組成, 其群落多樣性均較高(表2)。多樣性指數(shù)Shannon、Simpson、均勻度指數(shù)Heip均體現(xiàn)表層水體浮游細菌群落多樣性最高, 豐富度指數(shù)Chao 和Ace為底層水體最高, 文庫覆蓋率為表層高于底層。從標(biāo)準(zhǔn)差對比來看, 達里湖表、底層水體的浮游細菌群落豐富度較多樣性差異明顯。

    3.3 浮游細菌群落結(jié)構(gòu)組成的表、底層分布特征

    為具體探究達里湖表、底層水體中浮游細菌群落結(jié)構(gòu)的差異性特征, 采用非度量多維尺度分析(NMDS)方法, 繪制出樣點的二維分布圖(圖3), 其中不同樣點在二維圖中的距離可以反映群落結(jié)構(gòu)的相似程度, 樣點之間的距離越近則表示群落組成結(jié)構(gòu)越相似。

    通過表、底層各個取樣點中細菌群落結(jié)構(gòu)的分析, 可以發(fā)現(xiàn)表層水體各取樣點的NMDS1值主要分布在-0.4—0.2之間, NMDS2值主要分布在-0.2—0.45之間; 底層水各取樣點的NMDS1值主要分布在-0.4—0.3間, NMDS2值主要分布在-0.15—0.35間(圖3)。顯然, 夏季開放條件下達里湖水體浮游細菌菌群結(jié)構(gòu)差異性程度呈現(xiàn)表層>底層的特征。

    基于表、底層水體中浮游細菌群落結(jié)構(gòu)的差異, 對表層和底層水中的浮游細菌在門類水平上的變化進行了統(tǒng)計, 具體分析了達里湖表、底層水體中浮游細菌群類結(jié)構(gòu)特征及差異(圖4, 表3), 變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍細菌門(Cyanobacteria) 3種浮游細菌菌門豐度在表、底層水體中均占優(yōu)勢比例(表3), 其中, 變形菌門的標(biāo)準(zhǔn)差遠高于其他菌門: 變形菌門(Proteobacteria)的豐度最高, 在底層水(47.65%)中高于表層水(36.21%)中; 而放線菌門的豐度在底層水(24.89%)高于表層水(24.66%); 藍細菌門在表層水(11.22%)中高于底層水(8.99%)中。相對地, 擬桿菌門(Bacteroidetes)、軟壁菌門(Tenericutes)、異常球菌-棲熱菌門(Deinococcus-Thermus)、疣微菌門(Verrucomicrobia)、硬壁菌門 (Firmicutes)、浮霉?fàn)罹T(Plancto-mycetes)、酸桿菌門(Acidobacteria)和綠菌門(Chlorobi)等8種門類浮游細菌同樣在表層、底層湖水中的比例存在差異(表3)。

    表2 浮游細菌表、底層分布多樣性分析

    注: 黑體字代表平均值最高。

    圖2 浮游細菌稀釋性曲線

    Figure 2 Rarefaction curve of bacterioplankton community

    圖3 達里湖表、底層水體浮游細菌群落結(jié)構(gòu)的NMDS分析

    Figure 3 NMDS analysis of phytoplankton community structure in surface and bottom waters of Dali-nor Lake

    表3 表、底層湖水門類水平上浮游細菌群落結(jié)構(gòu)比例(%)

    圖4 表、底層水體門類水平浮游細菌群落結(jié)構(gòu)組成

    Figure 4 Composition of phytoplankton community in surface and bottom water bodies on phylum level

    3.4 環(huán)境因子與浮游細菌優(yōu)勢菌門之間的關(guān)系

    浮游細菌作為湖泊水體中主要的微生物種群, 其群落結(jié)構(gòu)變化會對水體環(huán)境因子的垂向變化出現(xiàn)一定程度的響應(yīng)。因此, 利用Pearson相關(guān)分析對達里湖表、底水體中門類水平上的浮游細菌物種豐富度和環(huán)境因子變化間的相互關(guān)系進行了詳細分析, 結(jié)果顯示表、底層水體中優(yōu)勢細菌門類受不同環(huán)境因子的影響程度也不相同(圖5)。

    注: 橫向坐標(biāo)為水體環(huán)境因子, 縱向為物種信息, 中間不同顏色代表Pearson 相關(guān)分析的相關(guān)系數(shù)r值, r值范圍為-1—1. 當(dāng)r值小于零, 則顯示兩者負相關(guān), 當(dāng)r值大于零, 則顯示兩者正相關(guān). "*"代表顯著性檢驗P值<0.05有顯著性差異, "**"代表有極顯著差異。

    Figure 5 Pearson analyzing result of planktonic bacterial flora and environmental factors in Dali-nor Lake

    表層水體中, TDS、DIP作為與浮游細菌優(yōu)勢菌門群落結(jié)構(gòu)相關(guān)性最顯著的環(huán)境因子, 在和不同優(yōu)勢門類浮游細菌間的相關(guān)性結(jié)果上也有體現(xiàn)。異常球菌-棲熱菌門(Deinococcus- Thermus)、硬壁菌門(Firmicutes)與TDS呈極顯著正相關(guān), 而與DIP呈顯著負相關(guān); 變形菌門(Proteobacteria)與這TN、SD、TDS 3個環(huán)境因子呈顯著負相關(guān); pH與擬桿菌門(Bacteroidetes)、軟壁菌門(Tenericutes)呈顯著負相關(guān); 綠菌門(Chlorobi)與DIP呈顯著負相關(guān) (圖5a)。

    底層水體中, 浮游細菌優(yōu)勢菌門群落結(jié)構(gòu)受TN、Chl和WT因子影響最顯著。TN與擬桿菌門(Bacteroidetes)、綠菌門(Chlorobi)、軟壁菌門(Tenericutes)呈顯著負相關(guān); DIP與軟壁菌門(Tenericutes)呈顯著正相關(guān); Chl與放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)、綠菌門(Chlorobi)呈顯著正相關(guān); WT與酸桿菌門(Acidobacteria)、異常球菌-棲熱菌門(Deinococcus- Thermus)和硬壁菌門(Firmicutes)呈顯著正相關(guān); 水深與酸桿菌門(Acidobacteria)、異常球菌-棲熱菌門(Deinococcus-Thermus)呈顯著負相關(guān)(圖5b)。

    4 討論

    本文利用16S rRNA基因測序技術(shù)闡明了達里湖浮游細菌群落的空間變化, 多樣性指數(shù)表現(xiàn)為表層高, 底層低, 而豐富度指數(shù)為表層低, 底層高, 每個樣品在菌群組成上表、底層差異顯著, 這與其中變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、藍細菌門(Cyanobacteria)是達里湖的主要菌門有關(guān), 這一分布特征也和大部分的淡水湖泊浮游細菌的組成相似[32,33]。

    在達里湖夏季水體中, Proteobacteria作為優(yōu)勢菌中數(shù)量最高的菌門, 在達里湖表、底層水體中均占據(jù)明顯優(yōu)勢, 其物種多樣性豐富, 生存方式多樣[34], 這也導(dǎo)致在表、底層湖水中Proteobacteria 豐度變化與湖水理化指標(biāo)參數(shù)的相互關(guān)系復(fù)雜多變; 而相對豐度也相對較高的 Actinobacteria, 其豐度變化往往會受到水體富營養(yǎng)化等因素影響[35], 隨著水深的增加, Actinobacteria與營養(yǎng)元素的相關(guān)性產(chǎn)生差異; 其次是Cyanobacteria, 其生長方式為在表層湖水進行光合作用, 消耗表層水中的營養(yǎng)物質(zhì), 當(dāng)營養(yǎng)物質(zhì)被消耗盡時, 他們便會隨細胞體下沉到底層水體或沉積物中繼續(xù)吸收營養(yǎng)物質(zhì), 進而造成達里湖表層水體到底層水體與氮、磷營養(yǎng)鹽的相關(guān)關(guān)系發(fā)生改變。

    達里湖屬典型內(nèi)陸封閉型湖泊, 常年湖水無外泄, 在夏季湖泊水體處于開放條件下, 表層水體更易受太陽輻射、大氣降水及河流輸入等外界因素的直接影響, 因此導(dǎo)致達里湖表層水的溫度、TDS值較高, 并且由于夏季適宜的水體溫度導(dǎo)致藻類大量繁殖, 造成表層水體的Chl值也較底層偏高, 進而顯示出表層水體的理化性質(zhì)主要受外界因素輸入的影響; 而底層水體中TN、DIP、COD等環(huán)境因子參數(shù)高于表層, 這是由于在湖泊底部, 底層水體與底泥接觸, 而底泥中的氮、磷等污染物釋放造成污染物濃度偏高。隨著環(huán)境因子在表、底層水體的變化, 其浮游細菌群落結(jié)構(gòu)也相應(yīng)發(fā)生改變??偟膩碚f, 湖泊水體中微生物群落的結(jié)構(gòu)特征與其環(huán)境理化性質(zhì)緊密相關(guān)并相互限制影響[36]。本研究發(fā)現(xiàn), 與達里湖浮游細菌呈顯著相關(guān)的環(huán)境因子分別為TN、DIP、Chl、PH、WT、SD、TDS、水深。TN、DIP及水深因子在表、底層均為主要環(huán)境因子。氮、磷元素循環(huán)是各種元素循環(huán)的中心, 微生物參與下的氮、磷循環(huán)是生態(tài)系統(tǒng)中的主要初級生產(chǎn)力, 在水生態(tài)環(huán)境中, 過量的氮、磷會導(dǎo)致水體發(fā)生富營養(yǎng)化。Haukka等[37]在對韋西湖微生物群落結(jié)構(gòu)的分析中發(fā)現(xiàn), 氮濃度與浮游細菌的組成直接相關(guān), 其中Verrucomicrobia能夠?qū)Φ谞I養(yǎng)鹽進行有效的利用, 而Actinobacteria則適宜低濃度氮磷營養(yǎng)環(huán)境。在本研究中達里湖浮游細菌對氮磷營養(yǎng)鹽的響應(yīng)并不明顯, 氮濃度與多數(shù)浮游細菌的關(guān)系為不顯著相關(guān), 且DIP作為對達里湖部分浮游細菌有顯著影響的磷營養(yǎng)鹽, 對Chlorobi與Tenericutes有顯著促進作用, 對Deinococcus-Thermus、Firmicutes呈抑制作用; 李玉華等[38]在對松花湖水體的研究中得出水深為主要環(huán)境因子, 這與本研究中水深變化影響水體環(huán)境因子改變進而影響菌群結(jié)構(gòu)產(chǎn)生差異的結(jié)果相同, 顯示了松花湖與達里湖浮游細菌菌群結(jié)構(gòu)均受到水深變化的影響。夏季開放條件下, 達里湖表層懸浮大量藻類, 使得Chl含量偏高且與底層存在顯著差異, 從Pearson相關(guān)分析中Chl對多數(shù)優(yōu)勢菌門的生長呈顯著促進作用, 除此之外, TDS也對表層的優(yōu)勢菌門具有顯著影響, 這與王歡等對碧流河水庫的研究存在一定差異[39]。顯然, 湖泊水體微生物結(jié)構(gòu)特征除受水體環(huán)境因子的調(diào)控外, 水域所處區(qū)域氣候條件及自身的區(qū)域特征也會對水體浮游細菌的菌群結(jié)構(gòu)及優(yōu)勢菌門的分布存在影響。

    總之, 在開放條件下的達里湖, 表層和底層水體受外界環(huán)境變化影響的程度存在差異, 這進一步導(dǎo)致了水體中環(huán)境因子的含量出現(xiàn)波動, 最終引起浮游細菌群落結(jié)構(gòu)及優(yōu)勢菌門的豐度出現(xiàn)變化。

    5 結(jié)論

    1) 夏季開放條件下, 湖水中變形菌門(Proteo-bacteria)、放線菌門(Actinobacteria)、藍細菌門(Cyanobacteria)和擬桿菌門(Bacteroidetes)為優(yōu)勢菌門, 且Proteobacteria含量在表、底層水體中中均占絕對優(yōu)勢: 底層水占比46.46 %, 表層水中占比36.21 %。表、底層水體中優(yōu)勢菌的占比存在一定差異, 這與水體環(huán)境差異有關(guān)。

    2)表層水體中浮游細菌的群落多樣性、分布均勻程度最高。不過, 底層水體的整體細菌豐富度要高于表層水。

    3) 水深變化引起的不同水體中環(huán)境因子含量的波動成為夏季達里湖浮游細菌菌群結(jié)構(gòu)和優(yōu)勢菌門豐度變化的關(guān)鍵因子, 即不同水深處浮游細菌群落結(jié)構(gòu)及優(yōu)勢菌門與不同水體環(huán)境因子間的相關(guān)性存在差異, 不同環(huán)境因子促進或抑制了相應(yīng)優(yōu)勢菌門的發(fā)育。其中, 表層水體中TDS、DIP對浮游細菌優(yōu)勢菌門影響最顯著, 底層水體中TN、Chl和WT為主要影響因子。整體上, 水體環(huán)境因子的表、底層差異成為影響浮游細菌的優(yōu)勢菌門群落結(jié)構(gòu)的重要因素。

    [1] Garcia S L, Salka I, Grossart H P, et al. Depth-discrete profiles of bacterial communities reveal pronounced spatio-temporal dynamics related to lake stratification[J]. Environmental Microbiology Reports, 2013, 5(4): 549–555.

    [2] Li M, Wei G, Shi W,et al. Distinct distribution patterns of ammonia-oxidizing archaea and bacteria in sediment and water column of the Yellow River estuary[J]. Scientific Reports, 2018, 8(1): 1584.

    [3] 陳禎, 何聃, 任麗娟. 溫度和營養(yǎng)鹽水平對淡水浮游細菌群落結(jié)構(gòu)的潛在影響[J].南京大學(xué)學(xué)報(自然科學(xué)), 2016, 52(4): 590–600.

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

    [5] Hanson C A, Fuhrman J A, Horner-Devine M C, et al. Beyond biogeographic patterns: processes shaping the microbial landscape[J]. Nature Reviews Microbiology, 2012, 10(7): 497–506.

    [6] 閆冰,李俊生,陸晴,等. 北京城市公園綠地土壤微生物群落碳源代謝活性特征[J].環(huán)境科學(xué)研究, 2019, 32(9): 1567–1574.

    [7] Yuan G, He J, Min C, et al. Factors dominating bacterioplankton abundance and production in the Nordic seas and the Chukchi Sea in summer 2012[J]. Acta Oceanologica Sinica, 2017, 36: 153–162.

    [8] De Figueiredo D R, Pereira M J, Correia A. Seasonal modulation of bacterioplankton community at a temperate eutrophic shallow lake[J]. World Journal of Microbiology Biotechnology, 2010, 26(6): 1067–1077.

    [9] Chen K, Allen J, Lu J R. Community Structures of Phytoplankton with Emphasison Toxic Cyanobacteria in an Ohio Inland Lake during Bloom Season[J]. Journal of Water Resource and Protection, 2017, 9(11): 1–29.

    [10] Hahn M W. The microbial diversity of inland waters[J]. Current Opinion in Biotechnology, 2006, 17(3): 256-261.

    [11] Sims A, Zhang Y, Gajaraj S, et al. Toward the development of microbial indicators for wetland assessment[J]. Water Research, 2013, 47(5): 1711–1725.

    [12] Jiang J G, Shen Y F. Development of the microbial communities in Lake Donghu in relation to water quality[J]. Environmental Monitoring and Assessment, 2007, 127(1-3): 227–236.

    [13] 程豹, 望雪, 徐雅倩, 等. 瀾滄江流域浮游細菌群落結(jié)構(gòu)特征及驅(qū)動因子分析[J]. 環(huán)境科學(xué), 2018, 39(8): 3649–3659.

    [14] 姜學(xué)霞, 朱木得, 陳雪曼, 等. 肇慶波海湖細菌群落多樣性研究[J]. 環(huán)境科學(xué)導(dǎo)刊, 2019, 38(5): 20–23.

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

    [16] 沈烽, 趙大勇, 黃睿, 等. 南京玄武湖浮游細菌群落結(jié)構(gòu)的季節(jié)變化及其與環(huán)境因子的關(guān)系[J]. 湖泊科學(xué), 2017, 29(3): 662–669.

    [17] 楊長明, 吳亞瓊, 王育來, 等. 南淝河表層沉積物細菌群落結(jié)構(gòu)特征及驅(qū)動因素[J]. 中國環(huán)境科學(xué), 2018, 38(9): 3552–3561.

    [18] 李文寶, 劉志嬌, 楊旭, 等. 內(nèi)蒙古高原達里諾爾湖夏季水體穩(wěn)定同位素變化特征[J]. 湖泊科學(xué), 2019, 31(2): 539–550.

    [19] 劉曉旭, 李暢游, 李文寶, 等. 冰封期達里諾爾湖同位素與營養(yǎng)鹽分布特征及關(guān)系的定量分析[J]. 湖泊科學(xué), 2015, 27(6): 1159–1167.

    [20] 李衛(wèi)平. 高原典型湖泊營養(yǎng)元素地球化學(xué)循環(huán)與重金屬污染研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2012.

    [21] 王旭陽. 基于3S技術(shù)的達里諾爾湖水深反演研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2017.

    [22] 李文寶, 李暢游, 劉曉旭, 等. 達里諾爾湖水體穩(wěn)定氫、氧同位素組成變化對結(jié)冰過程的響應(yīng)[J]. 地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報), 2015, 40(12): 2081–2090.

    [23] 范佳偉, 肖舉樂, 溫銳林, 等. 內(nèi)蒙古達里湖全新世有機碳氮同位素記錄與環(huán)境演變[J]. 第四紀研究, 2015, 35(4): 856–870.

    [24] Yang C Y, W Q, Simon P N, et al. Distinct Network Interactions in Particle-Associated and Free-Living Bacterial Communities during aBloom in a Plateau Lake[J]. Frontiers in Microbiology, 2017, 8: 1202–1227.

    [25] Angelika R, Herlemann D P R. Particle-Associated Differ from Free-Living Bacteria in Surface Waters of the Baltic Sea[J]. Frontiers in Microbiology, 2015, 6: 1297–1310.

    [26] 國家環(huán)境保護總局. 水和廢水監(jiān)測分析方法(第四版)[M]. 北京: 中國環(huán)境科學(xué)出版社, 2016, 35(1): 139–147.

    [27] Bates S T, Berg-Lyons D, Caporaso J G, et al. Examining the global distribution of dominant archaeal populations in soil[J]. The International Society for Microbial Ecology Journal, 2011, 5(5): 908–917.

    [28] Marti E, Balcázar, José Luis. Use of pyrosequencing to explore the benthic bacterial community structure in a river impacted by wastewater treatment plant discharges[J]. Research in Microbiology, 2014, 165(6): 468–471.

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

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

    [31] 劉志嬌. 達里諾爾湖水動力條件及氫氧穩(wěn)定同位素試驗研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2015.

    [32] 錢瑋, 范竟成, 王金虎, 等. 東太湖地區(qū)濕地水質(zhì)與浮游細菌群落結(jié)構(gòu)關(guān)聯(lián)性研究[J]. 生物技術(shù), 2017, 1: 90–96.

    [33] 劉晨宇, 祝亞楠, 王繼華, 等. 不同季節(jié)尾水排放對河道細菌群落結(jié)構(gòu)的影響[J]. 環(huán)境科學(xué)研究, 2018, 31(5): 893–900.

    [34] 丁軼睿, 李定龍, 張毅敏, 等. 滆湖底泥細菌群落結(jié)構(gòu)及多樣性[J]. 環(huán)境科學(xué)學(xué)報, 2017, 37(5): 1649–1656.

    [35] Allgaier M, Grossart H P. Diversity and Seasonal Dynamics of Actinobacteria Populations in Four Lakes in Northeastern Germany[J]. Applied and Environmental Microbiology, 2006, 72(5): 3489–3497

    [36] 萬甜, 何夢夏, 任杰輝, 等. 渭河流域水體細菌群落的環(huán)境響應(yīng)及生態(tài)功能預(yù)測[J]. 環(huán)境科學(xué), 2019, 40(08): 3588–3595.

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

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

    [39] 王歡, 趙文, 謝在剛, 等. 碧流河水庫細菌群落結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子[J]. 環(huán)境科學(xué), 2018, 39(8): 3660– 3669.

    Surface and bottom characteristics of bacterioplankton community in summer Dali-nor Lake and its key driving factors

    Du Lei, Li Changyou*, Li Wenbao, Shi Xiaohong, Yang Xu, Liu Jingjing

    IMAR Key Laboratory of Water Resources Protecting and Utilization, Inner Mongolia Agricultural University, Hohhot 010018, China

    Bacterioplankton is the key element in the material cycle and energyflow in lake, affecting the migration and transformation of the ecosystem. In this paper, we collect 32 water samples from the surface and bottom water in summer Dali-nor Lake, which locates in Inner Mongolia Plateau and is an inland closed lake. Based on the T-RFLP technology and 16S rRNA sequencing,the community structure of bacterioplankton and its vertical variation characteristics were analyzed in detail. In summer, the diversity and evenness of the bacterioplankton in surface water are higher than bottom water. Meanwhile, bottom water has the highest species richness.On the whole, the Proteobacteria, Actinobacteria, Cyanobacteria, and Bacteroidetes are the dominant phyla in summer Dali-nor Lake. The Proteobacteria is the most dominant phylum in all water layers, with the content of about 46.46% in bottom water, and about 36.21% in surface water, respectively. Furthermore, TDS, DIP in surface water are correlated with most dominant bacteria, TN, Chland WT in the bottom water are the main environmental factors influencing community distribution. In summary, there is a significant correlation between the surface and bottom planktonic bacterial community structure of Dali-norl Lake and different environmental factors, while the water environmental factors in surface and bottom waters have become important factors affecting the evolution of most dominant bacteria.

    bacterioplankton; change of surface and bottom water layers; summer; community structure characteristics; Dali-nor Lake

    10.14108/j.cnki.1008-8873.2021.06.002

    X172

    A

    1008-8873(2021)06-013-08

    杜蕾, 李暢游, 李文寶, 等. 夏季達里諾爾湖浮游細菌群落表、底層結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子[J]. 生態(tài)科學(xué), 2021, 40(6): 13–20.

    Du Lei, Li Changyou, Li Wenbao, et al. Surface and bottom characteristics of bacterioplankton community in summer Dali-nor Lake and its key driving factors[J]. Ecological Science, 2021, 40(6): 13–20.

    2020-05-11;

    2020-06-20

    國家自然科學(xué)基金地區(qū)項目(52160021, 51869020); 國家重點研發(fā)計劃—國際合作專項(2017YFE0114800); 內(nèi)蒙古自治區(qū)自然科學(xué)基金項目(2019MS05032); 內(nèi)蒙古自治區(qū)科技攻關(guān)項目(2020GG0009); 內(nèi)蒙古自治區(qū)高等學(xué)校“青年科技英才支持計劃”項目(NJYT-20-A14)

    杜蕾(1996—), 女, 內(nèi)蒙古赤峰市人, 碩士研究生, 主要研究方向為湖泊水環(huán)境演變及修復(fù), E-mail: 18847163393@163.com

    通信作者:李暢游, E-mail: nndlichangyou@163.com

    猜你喜歡
    浮游菌門底層
    航天企業(yè)提升采購能力的底層邏輯
    浞河浮游藻類的調(diào)查研究與水質(zhì)評價
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    《浮游》
    流行色(2019年10期)2019-12-06 08:13:26
    饑餓與重攝食對河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進展
    妊娠期糖尿病腸道菌群變化及臨床價值分析
    浮游藻類在水體PAHs富集與降解中的研究
    圍墾對椒江口夏季浮游植物群落結(jié)構(gòu)和多樣性的影響
    回到現(xiàn)實底層與悲憫情懷
    小說林(2014年5期)2014-02-28 19:51:47
    国产精品久久久av美女十八| 国产精品一国产av| av卡一久久| 日本av手机在线免费观看| 国产精品免费大片| 黄色怎么调成土黄色| 视频在线观看一区二区三区| av网站免费在线观看视频| 日韩不卡一区二区三区视频在线| 欧美精品人与动牲交sv欧美| 亚洲欧美成人综合另类久久久| 欧美另类一区| av电影中文网址| 熟女人妻精品中文字幕| 亚洲伊人色综图| 日韩在线高清观看一区二区三区| 免费观看av网站的网址| 国产国语露脸激情在线看| 精品午夜福利在线看| 久久99热这里只频精品6学生| 91精品伊人久久大香线蕉| 国产片内射在线| av黄色大香蕉| 高清av免费在线| 最新中文字幕久久久久| 午夜免费观看性视频| 丝瓜视频免费看黄片| 成人无遮挡网站| 校园人妻丝袜中文字幕| 久久人妻熟女aⅴ| 午夜福利视频精品| 黑人欧美特级aaaaaa片| 一级毛片 在线播放| 视频区图区小说| 日韩av不卡免费在线播放| 日韩人妻精品一区2区三区| 亚洲av欧美aⅴ国产| 哪个播放器可以免费观看大片| 在线观看三级黄色| 热re99久久国产66热| 国产精品熟女久久久久浪| www日本在线高清视频| 日日摸夜夜添夜夜爱| 国产成人av激情在线播放| 成人黄色视频免费在线看| 国产黄色免费在线视频| 免费女性裸体啪啪无遮挡网站| 插逼视频在线观看| 你懂的网址亚洲精品在线观看| 老司机影院毛片| 人体艺术视频欧美日本| 免费看不卡的av| 人妻一区二区av| 国产有黄有色有爽视频| 日韩一本色道免费dvd| 五月伊人婷婷丁香| 中国三级夫妇交换| 免费av中文字幕在线| 国产免费视频播放在线视频| 青春草国产在线视频| 18+在线观看网站| a级毛片在线看网站| 黑人猛操日本美女一级片| 丝袜脚勾引网站| 美女中出高潮动态图| 成人毛片60女人毛片免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产av成人精品| 老女人水多毛片| 插逼视频在线观看| √禁漫天堂资源中文www| 国产精品久久久久久精品电影小说| 女人被躁到高潮嗷嗷叫费观| 黄色一级大片看看| 人妻 亚洲 视频| 色哟哟·www| 少妇人妻久久综合中文| 久久99蜜桃精品久久| 日本色播在线视频| 成人毛片60女人毛片免费| 视频在线观看一区二区三区| 91aial.com中文字幕在线观看| 插逼视频在线观看| 日产精品乱码卡一卡2卡三| 搡老乐熟女国产| 建设人人有责人人尽责人人享有的| 日韩视频在线欧美| 亚洲精品aⅴ在线观看| 亚洲av.av天堂| 综合色丁香网| 日韩av在线免费看完整版不卡| 高清黄色对白视频在线免费看| 亚洲欧美成人综合另类久久久| 亚洲,欧美精品.| 久久久久久久亚洲中文字幕| 午夜影院在线不卡| av黄色大香蕉| 黄网站色视频无遮挡免费观看| 久久久亚洲精品成人影院| 国产日韩欧美在线精品| 极品少妇高潮喷水抽搐| 亚洲伊人色综图| 日韩视频在线欧美| 18禁观看日本| videos熟女内射| 免费人成在线观看视频色| 丝袜脚勾引网站| 欧美精品av麻豆av| 精品一区二区免费观看| 男女免费视频国产| 成人毛片a级毛片在线播放| 国产精品 国内视频| 寂寞人妻少妇视频99o| 最新的欧美精品一区二区| 黄片无遮挡物在线观看| 天天操日日干夜夜撸| 精品国产一区二区三区久久久樱花| 天堂中文最新版在线下载| 美女大奶头黄色视频| 在线观看免费日韩欧美大片| 日本vs欧美在线观看视频| 久久人人爽人人爽人人片va| 亚洲人与动物交配视频| 国产毛片在线视频| 制服诱惑二区| av在线播放精品| 国产亚洲精品久久久com| 日日摸夜夜添夜夜爱| 一本久久精品| 人妻人人澡人人爽人人| 色婷婷久久久亚洲欧美| 麻豆精品久久久久久蜜桃| 麻豆精品久久久久久蜜桃| 国产在线一区二区三区精| 欧美亚洲日本最大视频资源| 丝瓜视频免费看黄片| 国产色爽女视频免费观看| 国产成人精品一,二区| 另类亚洲欧美激情| 欧美激情极品国产一区二区三区 | 欧美日韩综合久久久久久| 亚洲精品久久久久久婷婷小说| 看非洲黑人一级黄片| 搡女人真爽免费视频火全软件| 亚洲天堂av无毛| 我的女老师完整版在线观看| 国产精品人妻久久久影院| 美女xxoo啪啪120秒动态图| 黄色毛片三级朝国网站| 国产片特级美女逼逼视频| 99久久人妻综合| 伦理电影大哥的女人| 国产极品粉嫩免费观看在线| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 亚洲人成网站在线观看播放| 国产伦理片在线播放av一区| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 国产精品久久久久成人av| 亚洲av在线观看美女高潮| 亚洲av电影在线进入| 国产极品粉嫩免费观看在线| 午夜激情av网站| 一区二区日韩欧美中文字幕 | 亚洲精品久久成人aⅴ小说| 在线精品无人区一区二区三| 少妇被粗大猛烈的视频| 免费观看无遮挡的男女| 亚洲国产成人一精品久久久| 国产熟女欧美一区二区| 国产麻豆69| 嫩草影院入口| a级毛色黄片| 一个人免费看片子| 日本黄大片高清| 欧美精品国产亚洲| 精品一区二区三区视频在线| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 欧美日韩视频精品一区| 视频区图区小说| 日韩大片免费观看网站| 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品电影小说| 亚洲国产成人一精品久久久| 国产免费现黄频在线看| 99久久中文字幕三级久久日本| 久久影院123| 视频中文字幕在线观看| 女的被弄到高潮叫床怎么办| 捣出白浆h1v1| 国产一区二区在线观看日韩| 曰老女人黄片| 亚洲少妇的诱惑av| 大片免费播放器 马上看| 一级a做视频免费观看| 国产精品久久久久成人av| 欧美精品高潮呻吟av久久| 大陆偷拍与自拍| 亚洲欧美一区二区三区黑人 | 纵有疾风起免费观看全集完整版| av片东京热男人的天堂| 最后的刺客免费高清国语| 国产深夜福利视频在线观看| 黄色配什么色好看| 美女主播在线视频| 色视频在线一区二区三区| 秋霞在线观看毛片| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人 | 色哟哟·www| 又黄又粗又硬又大视频| 久久久久久久国产电影| 高清视频免费观看一区二区| 久久久久久久久久久久大奶| 2018国产大陆天天弄谢| 亚洲国产最新在线播放| 欧美成人午夜免费资源| 欧美最新免费一区二区三区| 日韩制服骚丝袜av| 久久人人爽人人爽人人片va| 9热在线视频观看99| 精品一品国产午夜福利视频| 满18在线观看网站| 免费久久久久久久精品成人欧美视频 | 波多野结衣一区麻豆| 黄色怎么调成土黄色| 亚洲精品中文字幕在线视频| 岛国毛片在线播放| 欧美日韩精品成人综合77777| 美女xxoo啪啪120秒动态图| 一边摸一边做爽爽视频免费| 夜夜爽夜夜爽视频| 中文天堂在线官网| 精品久久国产蜜桃| 亚洲av免费高清在线观看| 亚洲第一av免费看| 1024视频免费在线观看| 一级毛片电影观看| 国产黄色视频一区二区在线观看| 精品人妻在线不人妻| 91午夜精品亚洲一区二区三区| 成人国语在线视频| 精品卡一卡二卡四卡免费| 91成人精品电影| 90打野战视频偷拍视频| 久久久久精品久久久久真实原创| 毛片一级片免费看久久久久| 国产乱人偷精品视频| 日本黄大片高清| 欧美国产精品一级二级三级| 国产精品99久久99久久久不卡 | 夜夜骑夜夜射夜夜干| 高清视频免费观看一区二区| 激情五月婷婷亚洲| 男人添女人高潮全过程视频| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品古装| 午夜福利网站1000一区二区三区| 一本久久精品| 亚洲av电影在线进入| √禁漫天堂资源中文www| 超碰97精品在线观看| 成人亚洲欧美一区二区av| 久久人人爽人人爽人人片va| 亚洲欧美日韩另类电影网站| 精品少妇黑人巨大在线播放| 久久精品国产自在天天线| 成人无遮挡网站| 成年女人在线观看亚洲视频| 婷婷色麻豆天堂久久| 欧美另类一区| 国产日韩欧美视频二区| 赤兔流量卡办理| 日本与韩国留学比较| 日日摸夜夜添夜夜爱| 国产日韩欧美亚洲二区| 男女免费视频国产| 两个人看的免费小视频| 欧美 日韩 精品 国产| 夫妻性生交免费视频一级片| 久久久国产一区二区| 国产 精品1| 国产毛片在线视频| 欧美bdsm另类| 国语对白做爰xxxⅹ性视频网站| 国产免费又黄又爽又色| 男女边摸边吃奶| 久久久久久久精品精品| videos熟女内射| 亚洲精品视频女| 午夜福利影视在线免费观看| 最近中文字幕2019免费版| av在线播放精品| 宅男免费午夜| 男女国产视频网站| 欧美变态另类bdsm刘玥| 欧美另类一区| 另类亚洲欧美激情| 国产片内射在线| 青春草视频在线免费观看| 我的女老师完整版在线观看| av黄色大香蕉| 亚洲av.av天堂| 成年av动漫网址| 亚洲美女黄色视频免费看| 水蜜桃什么品种好| 满18在线观看网站| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 色5月婷婷丁香| 免费av中文字幕在线| 亚洲第一av免费看| 亚洲综合色网址| 中文字幕亚洲精品专区| 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| 亚洲成人一二三区av| 亚洲av中文av极速乱| 欧美日韩av久久| 亚洲av日韩在线播放| 美国免费a级毛片| 久久久久久久久久久免费av| 丝袜美足系列| 色94色欧美一区二区| 亚洲欧美日韩卡通动漫| 伦理电影大哥的女人| 91在线精品国自产拍蜜月| 亚洲国产精品999| 久久久精品94久久精品| 免费少妇av软件| 色婷婷av一区二区三区视频| 亚洲av在线观看美女高潮| 26uuu在线亚洲综合色| 国产无遮挡羞羞视频在线观看| 2022亚洲国产成人精品| 麻豆精品久久久久久蜜桃| 亚洲精品自拍成人| 最近的中文字幕免费完整| 欧美精品一区二区大全| 亚洲精品国产av成人精品| 夜夜爽夜夜爽视频| 99热国产这里只有精品6| 免费观看a级毛片全部| xxxhd国产人妻xxx| 丁香六月天网| 亚洲av国产av综合av卡| 亚洲高清免费不卡视频| 久热这里只有精品99| 日产精品乱码卡一卡2卡三| 高清不卡的av网站| 高清在线视频一区二区三区| a级毛片在线看网站| 91精品伊人久久大香线蕉| 久久这里有精品视频免费| 精品国产一区二区久久| 国产精品一区www在线观看| 午夜福利,免费看| 日韩av免费高清视频| 精品一区二区三卡| 少妇的逼好多水| 国产xxxxx性猛交| 亚洲国产精品一区二区三区在线| 成年人午夜在线观看视频| 高清黄色对白视频在线免费看| 日本欧美国产在线视频| 久久久久久久久久久久大奶| 99热全是精品| 午夜av观看不卡| 亚洲国产欧美日韩在线播放| 波多野结衣一区麻豆| 丝袜在线中文字幕| 亚洲av成人精品一二三区| 男人舔女人的私密视频| 亚洲国产av影院在线观看| 国产乱人偷精品视频| 99香蕉大伊视频| 热re99久久精品国产66热6| av播播在线观看一区| 欧美3d第一页| 国产极品粉嫩免费观看在线| 97在线人人人人妻| 久久久久视频综合| 午夜免费观看性视频| 欧美另类一区| 欧美+日韩+精品| 综合色丁香网| 色94色欧美一区二区| 午夜91福利影院| 国产高清不卡午夜福利| 国产精品久久久久成人av| 男女免费视频国产| 侵犯人妻中文字幕一二三四区| 国产极品天堂在线| 如何舔出高潮| 亚洲综合色惰| 亚洲三级黄色毛片| 欧美另类一区| 亚洲在久久综合| 日韩成人av中文字幕在线观看| 欧美日本中文国产一区发布| 一二三四中文在线观看免费高清| 黄色一级大片看看| 综合色丁香网| 午夜免费观看性视频| 一级毛片电影观看| 久久久久精品人妻al黑| 大片免费播放器 马上看| 香蕉丝袜av| 美女国产高潮福利片在线看| 亚洲一区二区三区欧美精品| 亚洲国产欧美日韩在线播放| 精品一区二区免费观看| 亚洲国产看品久久| 成人18禁高潮啪啪吃奶动态图| 亚洲成av片中文字幕在线观看 | 在线观看国产h片| 久久av网站| 国产av国产精品国产| 免费av不卡在线播放| 国产在视频线精品| 国产永久视频网站| 亚洲精品456在线播放app| 精品国产乱码久久久久久小说| 国产乱人偷精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 午夜视频国产福利| 亚洲精品久久午夜乱码| 波野结衣二区三区在线| 男人爽女人下面视频在线观看| 亚洲激情五月婷婷啪啪| 国产一级毛片在线| 桃花免费在线播放| 久久婷婷青草| 亚洲欧洲精品一区二区精品久久久 | 一级片'在线观看视频| 在线观看www视频免费| 制服丝袜香蕉在线| 好男人视频免费观看在线| 捣出白浆h1v1| 欧美性感艳星| 久热这里只有精品99| 自拍欧美九色日韩亚洲蝌蚪91| av黄色大香蕉| 亚洲国产av影院在线观看| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 亚洲伊人色综图| 啦啦啦在线观看免费高清www| tube8黄色片| 我的女老师完整版在线观看| 高清在线视频一区二区三区| 26uuu在线亚洲综合色| 9热在线视频观看99| 插逼视频在线观看| 精品卡一卡二卡四卡免费| 免费大片黄手机在线观看| 观看美女的网站| kizo精华| 久久ye,这里只有精品| 久热这里只有精品99| 国产国拍精品亚洲av在线观看| 婷婷成人精品国产| 欧美精品一区二区免费开放| 免费黄色在线免费观看| 欧美精品人与动牲交sv欧美| 亚洲欧美中文字幕日韩二区| 日本-黄色视频高清免费观看| 一本久久精品| 亚洲一区二区三区欧美精品| 男人舔女人的私密视频| 午夜激情久久久久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 天天躁夜夜躁狠狠久久av| 精品久久蜜臀av无| 高清在线视频一区二区三区| 最新中文字幕久久久久| 飞空精品影院首页| 一区在线观看完整版| 性色av一级| 日本色播在线视频| 中文字幕亚洲精品专区| 最新中文字幕久久久久| 久久久久久人妻| 日本午夜av视频| 国产一级毛片在线| 日本午夜av视频| 热re99久久精品国产66热6| 亚洲av国产av综合av卡| 99热6这里只有精品| 亚洲av中文av极速乱| 最近的中文字幕免费完整| 亚洲精品乱码久久久久久按摩| www日本在线高清视频| 最近最新中文字幕大全免费视频 | 中文字幕免费在线视频6| 80岁老熟妇乱子伦牲交| 亚洲av福利一区| 亚洲丝袜综合中文字幕| 久久久久精品久久久久真实原创| 亚洲国产精品国产精品| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| 国产av精品麻豆| 老司机亚洲免费影院| 伊人亚洲综合成人网| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av涩爱| 毛片一级片免费看久久久久| 插逼视频在线观看| 亚洲国产精品一区二区三区在线| 久久精品夜色国产| 91aial.com中文字幕在线观看| 在线亚洲精品国产二区图片欧美| www.色视频.com| 精品一区二区三卡| av天堂久久9| 伦精品一区二区三区| 亚洲情色 制服丝袜| 十八禁网站网址无遮挡| 午夜福利视频精品| 成年美女黄网站色视频大全免费| 国产精品一区www在线观看| 91国产中文字幕| 亚洲精品国产av蜜桃| 草草在线视频免费看| 一级爰片在线观看| av一本久久久久| 免费看光身美女| 国语对白做爰xxxⅹ性视频网站| 中文天堂在线官网| 超碰97精品在线观看| 国产精品不卡视频一区二区| 午夜激情久久久久久久| 18+在线观看网站| 午夜免费鲁丝| 欧美日韩综合久久久久久| 成年人午夜在线观看视频| 一区二区三区精品91| 久久久久久久国产电影| 日韩大片免费观看网站| 精品国产一区二区三区四区第35| 欧美少妇被猛烈插入视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女脱内裤让男人舔精品视频| av女优亚洲男人天堂| 国产精品偷伦视频观看了| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 一级爰片在线观看| 国产片内射在线| www.av在线官网国产| 午夜91福利影院| 22中文网久久字幕| 韩国高清视频一区二区三区| 欧美精品一区二区免费开放| 免费黄网站久久成人精品| 爱豆传媒免费全集在线观看| 蜜桃在线观看..| 欧美日韩视频高清一区二区三区二| 国产精品.久久久| 亚洲精品av麻豆狂野| 美女脱内裤让男人舔精品视频| 亚洲人与动物交配视频| 日本色播在线视频| 在线观看人妻少妇| 一级a做视频免费观看| 欧美人与善性xxx| 日韩欧美精品免费久久| 永久网站在线| 人妻人人澡人人爽人人| 满18在线观看网站| 麻豆精品久久久久久蜜桃| 纵有疾风起免费观看全集完整版| 亚洲人与动物交配视频| 亚洲熟女精品中文字幕| 丁香六月天网| 黑人猛操日本美女一级片| 亚洲成av片中文字幕在线观看 | 亚洲精品国产色婷婷电影| 欧美+日韩+精品| 日韩免费高清中文字幕av| 视频在线观看一区二区三区| 亚洲精品国产av蜜桃| 精品午夜福利在线看| 国产又色又爽无遮挡免| av国产精品久久久久影院| 在线观看人妻少妇| √禁漫天堂资源中文www| 国产探花极品一区二区| av国产久精品久网站免费入址| 纵有疾风起免费观看全集完整版| 黄色视频在线播放观看不卡| 91国产中文字幕| 久热久热在线精品观看| 久久精品国产自在天天线| 国产日韩欧美亚洲二区| 日本wwww免费看| 天天影视国产精品| 欧美性感艳星| 久久久国产精品麻豆| 婷婷色综合www| 这个男人来自地球电影免费观看 | 精品福利永久在线观看| 婷婷成人精品国产| 亚洲一级一片aⅴ在线观看| 考比视频在线观看| √禁漫天堂资源中文www| 精品人妻熟女毛片av久久网站| 成人无遮挡网站| 视频在线观看一区二区三区|