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

    龍景湖沉積物的細菌群落垂向分布特征

    2017-06-28 12:36:29牛鳳霞吉芳英沈秋實顏海波重慶大學三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室低碳綠色建筑國際聯(lián)合研究中心重慶400045
    中國環(huán)境科學 2017年6期
    關(guān)鍵詞:中層條帶深層

    牛鳳霞,吉芳英,趙 艮,張 倩,沈秋實,何 強,顏海波 (重慶大學三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室,低碳綠色建筑國際聯(lián)合研究中心,重慶 400045)

    龍景湖沉積物的細菌群落垂向分布特征

    牛鳳霞,吉芳英*,趙 艮,張 倩,沈秋實,何 強,顏海波 (重慶大學三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室,低碳綠色建筑國際聯(lián)合研究中心,重慶 400045)

    為探索不同沉積時間條件下沉積物細菌群落的垂向變化及其在湖泊生態(tài)系統(tǒng)物質(zhì)遷移轉(zhuǎn)化過程中所發(fā)揮作用,對重慶市園博園龍景湖的原有河道(OR)和新形成庫灣(NB)沉積物進行了分層采樣,利用變性梯度凝膠電泳(DGGE)技術(shù)并結(jié)合環(huán)境參數(shù)進行冗余分析(RDA).結(jié)果表明,OR沉積物的細菌多樣性指數(shù)(H)和豐度(S)高于NB,垂向上H、S和均勻度(E)都由表層至深層先減小后增大.系統(tǒng)發(fā)育分析顯示,龍景湖沉積物主要包含了7個門類的細菌,OR表層沉積物以δ-變形菌(Deltaproteobacteria)、擬桿菌(Bacteroidetes)和綠彎菌(Chloroflexi)為主,中層以綠彎菌為主,深層以δ-變形菌和綠彎菌為主;而NB表層和中層的主要菌種與OR深層一致,其他菌種還有Ignavibacteriae、放線菌(Actinobacteria)、厚壁菌 (Firmicutes)和綠菌門(Chlorobi).RDA表明,沉積物中總氮(TN)、總硫(TS)、孔隙水的總有機碳(TOC)濃度以及平均粒徑是影響細菌群落結(jié)構(gòu)和分布的主要因子.通過DGGE技術(shù)所獲得的細菌多與有機質(zhì)的降解有關(guān).OR的TOC、TN含量高于NB,相應(yīng)的沉積物中δ-變形菌、擬桿菌、Ignavibacteriae等與有機質(zhì)降解相關(guān)的細菌種類和數(shù)量較多.

    細菌群落;垂向分布;變性梯度凝膠電泳(DGGE);冗余分析(RDA)

    沉積物中蘊含著豐富的微生物[1-2],它們在生物化學轉(zhuǎn)化過程中發(fā)揮著重要作用,包括有機質(zhì)(OM)的礦化和污染物的降解[4].對微生物至關(guān)重要的底物和電子受體,根據(jù)能量代謝的效率,沿沉積物垂向逐漸被消耗.這些與深度相關(guān)的生物地球化學性質(zhì)梯度能夠為不同的微生物代謝活動提供生態(tài)位[5].研究表明,細菌群落沿水平和垂向有明顯的空間分異[3-16],研究不同區(qū)域沉積物中細菌群落的組成和代謝功能是了解沉積物生物圈演變和發(fā)展的一個關(guān)鍵[17].由于沉積物中僅有 0.25%的微生物能夠進行獨立培養(yǎng),基于 16S rRNA的變性梯度凝膠電泳(DGGE)被廣泛應(yīng)用于沉積物細菌生態(tài)學研究[16,18-19].Koizumi等[5]研究 Biwa湖細菌群落結(jié)構(gòu)時發(fā)現(xiàn),沉積物中Nitrospira-like非常豐富,其相對豐度隨深度遞增;Shao等[18]研究了太湖2個富營養(yǎng)化湖區(qū)的細菌群落組成,結(jié)果表明變形菌是沉積物中的主要細菌門,但不同區(qū)域沉積物中的優(yōu)勢菌群并不相同;史春瀟等[15]分析了甲烷礁鄰區(qū)沉積物中細菌群落結(jié)構(gòu)的垂向特征,結(jié)果顯示細菌群落結(jié)構(gòu)隨深度變化明顯,表層以 δ-變形菌為主,中層以 α-變形菌和δ-變形菌為主,深層以α-變形菌及未培養(yǎng)微生物為主.目前關(guān)于沉積物中細菌群落沿深度變化的研究多集中在海洋[4,7,9]、土壤[20]、濕地[12]和淺水湖泊[18],較少關(guān)注深水湖泊沉積物的細菌群落結(jié)構(gòu)差異及其與生物地球化學性質(zhì)垂向梯度的相關(guān)性[5].深水湖泊通常具有儲水量大、水力停留時間長、水溫分層等特殊性[21].許多環(huán)境變量,如溶解氧(DO)[22]、pH值[23-24]、無機物含量[25]等,都與垂向水溫梯度密切相關(guān)[25].水溫分層使得混合較好的表水層與下層水交換不暢,下水層溶解氧較低甚至出現(xiàn)厭氧狀態(tài)[21].Kristen等[20]的模擬實驗發(fā)現(xiàn),厭氧條件會顯著改變土壤微生物的群落組成.

    重慶市園博園龍景湖是典型的新建河道型深水水庫,2012年由原龍景湖水庫大壩加高蓄水而成[26].龍景湖的水面總面積約 0.67km2,水文水質(zhì)特征表現(xiàn)為水位較深(平均水深11~22m)[27],水流滯緩,水溫分層顯著,底層水體溶解氧含量低等[28].水庫建成初期,龍景湖水質(zhì)總體呈中富營養(yǎng)化狀態(tài)[29].由于龍景湖閘壩攔截蓄水前未進行河道清淤處理和建設(shè)初期施工攔截蓄水形成的新淹沒底部區(qū)域的清庫工作,龍景湖底部的沉積物類型存在差異.整個湖區(qū)沉積物大致可分為 2類:園博園建設(shè)前已存在的水體區(qū)域(原有湖區(qū))和園區(qū)建成后新增加的區(qū)域(新增淹沒區(qū))[29].本研究利用DGGE技術(shù),以龍景湖不同區(qū)域沉積物為研究對象,分析細菌群落結(jié)構(gòu)的垂向分布特征及其與環(huán)境因子的相關(guān)性,以期為湖泊沉積物細菌多樣性及其生態(tài)功能的研究提供基礎(chǔ).

    1 材料與方法

    1.1 采樣點布設(shè)及預(yù)處理

    在龍景湖的原有湖區(qū)的河道和新增淹沒區(qū)的庫灣選取代表性樣點,于2015年10月26日用無擾動柱狀采樣器(Corer 60,Uwitec,Austria)采集沉積物樣品,分別標記為 OR(原有河道, 29°41′00″N,106°33′00″E)和 NB(新淹沒庫灣, 29°40′53″N,106°32′55″E),每個樣點采集3根平行柱樣,樣點位置如圖1所示.現(xiàn)場將沉積物以1cm間隔分割,并測定每層沉積物的氧化還原電位(Eh).測定完成后迅速將泥樣放入滅菌的聚乙烯離心管,用便攜式冰箱低溫運回實驗室.OR共分割10個梯度,NB共分割6個梯度.

    每個樣品分為2份,分別于4℃(理化性質(zhì)測定)和?80℃(DNA提取)保存.采集3根平行柱樣的樣品分別測定理化指標,各項指標取平均值.有機質(zhì)(OM)的測定選取新鮮泥樣,用燒失量(LossOn Ignition, LOI)表示.CT15RT型高速冷凍離心機離心(8000r/min,10min)后獲取的上清液為孔隙水,用于測定 TOC濃度.離心后的泥樣經(jīng)冷凍干燥、研磨過100目篩后用元素分析儀測定沉積物的總硫(TS)和總氮(TN).分別從 3個平均柱樣中提取每層沉積物的細菌 DNA 進行混勻[11,30],用于PCR擴增等后續(xù)步驟.

    1.2 DGGE提取和分析方法

    1.2.1 DNA提取 總 DNA提取采用 Fast DNATMSPIN Kit For Soil試劑盒,步驟依照廠家手冊.以樣品基因組 DNA為模板,采用細菌通用引物341F,926R擴增樣品16S rDNA高變區(qū)序列,引物信息見表1.PCR擴增體系(50μL)為:10×PCR buffer 5μL;dNTP (2.5mmol/L)3.2μL;rTaq(5U/μL) 0.4μL;GC-341F(20μmol/L)1μL;926R(20μmol/L) 1μL;模板DNA 50ng;補雙蒸水至50μL.PCR擴增程序為:94℃預(yù)變性5min;94℃變性40s,55℃復(fù)性40s,72℃延伸 1min,30個循環(huán);最終 72℃延伸10min.PCR產(chǎn)物采用 OMEGA公司 DNA Gel Extraction Kit純化回收.

    表1 PCR引物序列Table 1 Primers of Polymerase Chain Reaction (PCR)

    1.2.2 變性梯度凝膠電泳(DGGE) 取 10 μL PCR的產(chǎn)物進行變性梯度凝膠電泳(DGGE)分析.采用變性梯度為35%~55%、濃度為8%的聚丙烯酰胺凝膠在 1×TAE緩沖液中 80V 60℃下電泳16h.變性梯度凝膠電泳(DGGE)完畢后,采用銀染法染色、步驟包括:固定液固定、Milli-Q純水清洗、銀染液染色、顯色液顯色、終止液終止反應(yīng).

    1.2.3 DGGE圖譜中優(yōu)勢條帶的回收與測序在紫外光照射下,用滅菌的手術(shù)刀切下OR和NB均出現(xiàn)的條帶和部分特征條帶(本研究選取 16條具有代表性的條帶,編號分別為1、5、6、7、8、9、12、13、14、15、17、19、21、22、23、24),并采用 OMEGA 公司的 Poly-Gel DNA Extraction Kit回收目的條帶.以2μL回收產(chǎn)物為模板,341F/926R為引物進行PCR擴增(體系和引物不變).將重新擴增的DNA片段切膠回收、純化后,連接到Pmd18-T載體上,并轉(zhuǎn)化至DH5α感受態(tài)細胞中,篩選陽性克隆,進行序列測定.

    1.2.4 數(shù)據(jù)分析 采用 Quantity one 4軟件對DGGE電泳圖譜中每個樣品的條帶數(shù)目、條帶密度進行數(shù)字化分析,用戴斯系數(shù)表征DGGE圖譜中沉積物不同層的相似程度,一般戴維系數(shù)越大,相似性越高.Shannon多樣性指數(shù)(H)、Pielou均勻度指數(shù)(E)和Patrick豐度指數(shù)(S)的計算公式如下所示:

    式中:pi為樣品中單一條帶的強度在該樣品所有條帶總強度中所占的比率;N為DGGE圖譜單一泳道上所有條帶的豐度;Ni為第 i條帶的豐度;S是某樣品中所有條帶數(shù)目總和.

    將基因測序結(jié)果與GenBank數(shù)據(jù)庫進行比對,匹配同源性最高的序列.所獲得序列與其他序列的相似度≥97%時,被認為隸屬于同一種系.采用DNAstar和Cluster軟件對測序結(jié)果進行分析,下載最相似的菌株序列作為系統(tǒng)發(fā)育樹的參考序列.然后采用MEGA軟件,Neighbor-joining法構(gòu)建系統(tǒng)發(fā)育樹,自展數(shù)(bootstrap)為 1000.使用軟件CANOCO 5.0對沉積物樣品的基本理化指標進行主成分分析(PCA)和冗余分析(RDA).

    2 結(jié)果與分析

    2.1 環(huán)境參數(shù)和沉積物理化性質(zhì)

    龍景湖原有河道(OR)和新形成庫灣(NB)沉積物-水界面上覆水的性質(zhì)見表 2.OR的水深較深,達 16.1m,上覆水溶解氧(DO)濃度較低(1.14mg/L).OR沉積物上覆水的pH值呈弱堿性(7.15),而 NB呈弱酸性(6.60).細菌群落結(jié)構(gòu)的差異與沉積物、土壤中氧化還原電位(Eh)的急劇變化有關(guān)[4,7,31].依據(jù)沉積物中 Eh的垂向變化趨勢,本研究將龍景湖沉積物樣品分為:表層(0~3cm),中層(3~7cm),深層(7~10cm),不同層的理化指標見表3.水平方向上,原有河道(OR)

    表2 龍景湖各采樣點基本水質(zhì)參數(shù)Table 2 Parameters of water quality at each sample site in Longjing Lake

    沉積物的OM含量和孔隙水的TOC濃度均高于新形成庫灣(NB).垂向上,OR沉積物表層至深層的Eh整體呈遞增趨勢,而TN、TS、OM、平均粒徑及孔隙水TOC濃度均逐漸降低;而NB表層沉積物的TN含量略低于中層.NB表層沉積物(0~3cm)的平均粒徑和 TS含量明顯高于 OR,而中層相反,這也說明了近年來兩個區(qū)域的沉積環(huán)境存在較大差異.

    表3 龍景湖沉積物理化性質(zhì)Table 3 Physical and chemical characteristics of the sediments in Longjing Lake

    2.2 不同深度DGGE圖譜分析

    圖2 PCR-DGGE技術(shù)分析原有河道和新形成庫灣沉積物細菌群落垂向分布的指紋圖譜Fig.2 PCR-DGGE profiles of bacterial communities in the sediments of OR and NB Bay

    龍景湖沉積物中不同深度樣品中分離到 10到25條不同的條帶,細菌的DGGE電泳圖譜如圖2所示.原有河道(OR)沉積物較新形成庫灣(NB)具有更高的細菌多樣性,其中條帶數(shù)目最多的層位是OR的表層0~1cm,最少的為NB的中層4~5cm.不同層細菌群落結(jié)構(gòu)差異較大,而同一層內(nèi)不同深度相似性較高[32].

    圖3 龍景湖沉積物細菌多樣性、均勻度和豐度指數(shù)的垂向分布Fig.3 Indexes of H, E and S in the sediments of Longjing Lake along depth

    根據(jù)不同深度沉積物中細菌的豐度(S)、香農(nóng)指數(shù)(H)和均勻度(E)可知(圖 3),OR表層0~1cm的香農(nóng)指數(shù)最高(2.89).OR和NB的香農(nóng)指數(shù)自表層向下均呈現(xiàn)先減小后增大的趨勢. OR中下層(4~8cm)和NB中層(4~5cm)沉積物的均勻度指數(shù)相對較低,反映了中層沉積物中細菌群落結(jié)構(gòu)發(fā)生了變化.

    2.3 細菌群落組成

    本研究所選取的16個條帶大部分與數(shù)據(jù)庫中已鑒定或未培養(yǎng)的細菌序列具有很高的同源性,根據(jù)比對結(jié)果構(gòu)建的系統(tǒng)發(fā)育樹見圖 4.龍景湖沉積物主要涵蓋了 7個類群的細菌,包括變形菌(Proteobacteria,條帶1、12、15、17、22、24)、擬桿菌門(Bacteroidetes,條帶 8、13、19)、綠彎菌門(Chloroflexi,條帶 5、6)、Ignavibacteriae(條帶9、21)、放線菌(Actinobacteria,條帶7)、厚壁菌門(Firmicutes,條帶14)和綠菌門(Chlorobi,條帶23).其中,變形菌主要包含了 β-變形菌(Betaproteobacteria,條 帶 1)和 δ-變 形 菌(Deltaproteobacteria,條帶12、15、17、22、24)2個亞類.

    圖4 龍景湖沉積物中細菌16S rRNA序列系統(tǒng)發(fā)育樹Fig.4 Phylogenetic tree of bacteria 16S rRNA sequence in sediments of Longjing Lake

    龍景湖沉積物中以δ-變形菌和綠彎菌的為主,不同樣點和層位的細菌群落組成差異較大.OR表層沉積物以δ-變形菌、擬桿菌和綠彎菌為主,中層以綠彎菌為主,深層以 δ-變形菌和綠彎菌為主;而NB的表層和中層均以δ-變形菌和綠彎菌為主.擬桿菌門和 Ignavibacteriae門主要出現(xiàn)在OR,但在沉積物中的富集區(qū)域不同.擬桿菌門的條帶8和13主要集中在表層、中層,而條帶19則出現(xiàn)在深層.Ignavibacteriae門的條帶9和21分別集中在沉積物的中層和深層.綠菌門、放線菌、厚壁菌在OR和NB沉積物中都有富集,其中綠菌門(條帶23)富集在OR的表層和 NB的中層,放線菌門(條帶 7)廣泛存在于OR的深層和NB的表層.僅在OR表層0~1cm出現(xiàn)的厚壁菌門(條帶14),在NB沉積物中普遍存在.

    3 討論

    3.1 不同深度細菌群落結(jié)構(gòu)分析

    采用DGGE技術(shù)分析龍景湖原有河道(OR)和新形成庫灣(NB)沉積物的細菌群落結(jié)構(gòu)可知,OR和 NB的主要類群為變形菌,δ-變形菌為主要亞群.太湖[18]、洱海[33]、德國Saidenbach飲用水水庫[1]、南海北部九龍甲烷礁鄰區(qū)[15]沉積物中的主要細菌門也為變形菌,但不同區(qū)域的優(yōu)勢亞群不同.Saidenbach水庫、太湖沉積物的優(yōu)勢亞群分別為β-和δ-變形菌、β-變形菌.洱海表層沉積物以β-變形菌為主,而深層主要為δ-變形菌;甲烷礁鄰區(qū)表層沉積物以δ-變形菌為主,中層以α-、δ-變形菌為主,深層以α-變形菌及未培養(yǎng)微生物為主.δ-變形菌綱是變形菌中唯一主要由厭氧菌組成的分類[34],龍景湖沉積物中的條帶 17和 24分別與 δ-變形菌綱的 Syntrophorhabdus aromaticivorans和Smithella propionica親緣關(guān)系密切(相似性 97%),均為專性厭氧菌.上述兩種菌的代謝類型和在沉積物中的分布不同, Syntrophorhabdus aromaticivorans能夠與產(chǎn)甲烷菌耦合,將苯酚降解為醋酸[35],主要分布在OR的深層,而Smithella propionica是一種互養(yǎng)的、丙酸鹽氧化菌,除了降解丙酸鹽也能夠利用甲酯或丁酸鹽緩慢生長[36],在OR的表層和深層以及NB的中層富集.

    擬桿菌門在高分子有機物的降解過程中發(fā)揮著重要作用[6],DGGE中的條帶8和13分別與擬桿菌門的 Candidatus Aquirestis calciphila和Parasegetibacter luojiensis遺傳距離緊密(相似性>97%).條帶8在OR和NB的表層沉積物中均有出現(xiàn),主要集中在 OR的表層(0~1cm),這是由于Candidatus A. calciphila是淡水浮游生物的重要組成部分,主要出現(xiàn)在靜止內(nèi)陸水體的浮游帶,種群豐度的季節(jié)性變化明顯,在春季出現(xiàn)的大幅度增值[37]可能導(dǎo)致一些浮游生物沉降到了沉積物表層.條帶13主要在沉積物的表層以及OR的中層(5~6cm)富集,這與 Parasegetibacter luojiensis的生理學特征有關(guān).Parasegetibacter luojiensis最早從中國新疆胡楊樹林的土壤中分離出來,為異養(yǎng)型、嚴格好氧菌,擁有完整的滑行-運動蛋白[38],該菌株可能通過滑行遷移[39]到 OR沉積物的中層.條帶 19也屬于擬桿菌門,但與 Prolixibacter bellariivorans的同源性較低(90%),主要存在于OR的深層沉積物中.

    根據(jù)細菌的遺傳學特性,Podosokorskaya等[40]認為,綠菌門和 Ignavibacteriae都起源于擬桿菌門.綠菌門與 Ignavibacteriae的分歧始于綠菌門能夠利用太陽能,代謝類型由異養(yǎng)轉(zhuǎn)變?yōu)樽责B(yǎng),而 Ignavibacteriae依然為化能有機異養(yǎng)細菌,并且演變?yōu)榭衫酶喾N類的含碳底物作為電子受體.龍景湖沉積物所提取的條帶23與綠菌門的Chlorobaculum limnaeum相似性高達99%,普遍分布在 OR的表層沉積物中.Chlorobaculum limnaeum是一種綠硫菌,適宜生活在湖泊,需要利用硫化物和硫或 H2作為光合作用的電子供體[41],最終氧化為硫酸鹽.田濤等[27]在研究龍景湖原水庫沉積物-水界面硫酸鹽的分布時發(fā)現(xiàn),原水庫沉積物中硫酸鹽的濃度隨深度遞增,而表層沉積物中硫酸鹽的濃度較低,這表明 OR表層沉積物中除了含有綠硫菌,可能同時還存在著硫酸鹽還原菌.Sinkko等[6]的研究表明硫酸鹽還原菌(δ-變形菌)在海岸和海洋沉積物(尤其是表層)中很豐富.

    Ignavibacteriae是一個新確立的門[40],與綠菌門的主要差異在于其兼性厭氧、專性有機營養(yǎng)型、無綠色體、不能進行光合作用[42].分別與條帶 9和 21具有較高同源性(相似性<90%)的是Melioribacter roseus和Ignavibacterium album.M. roseus可利用單糖、二糖或多糖進行好氧呼吸、發(fā)酵或還原多種電子供體(如亞硝酸鹽、Fe(III)、As(V)).該菌種的高水解性使其成為富含 OM沉積物中的重要分解者[40],主要分布在龍景湖 OR和NB的中層.而Ignavibacterium album有多種代謝機制,可編碼氧化酶和還原酶,在好氧和缺氧環(huán)境都能進行有機異養(yǎng)代謝[42],可同時分布在龍景湖OR的表層和深層.

    綠彎菌是水合物較少而 OM豐富沉積物的主要類群[43],與條帶 5、6同源性較高的分別為Dehalogenimonas alkenigignens(相似性 88%)和Dehalogenimonas lykanthroporepellens(相似性90%),都屬于 organohalide respiring bacteria (OHRB)[44].OHRB的共性特征為能在厭氧條件下發(fā)生還原性脫氯反應(yīng),利用鹵代化合物作為呼吸作用的電子受體,氫氣作為電子供體[45].鹵代化合物的存在通常與工業(yè)污染相關(guān)[46],其中 1,2-二氯甲烷(1,2-DCP)已經(jīng)是環(huán)境和飲用水質(zhì)的一個威脅,尤其是已受污染地區(qū)[44].迄今為止,僅OHRB可以通過完全脫氯作用將1,2-DCP轉(zhuǎn)化成丙烯或丙烷,減弱污染物的毒性.基于綠彎菌門(條帶5,6)在OR和NB各層沉積物中普遍存在且數(shù)量較多,對龍景湖沉積物中綠彎菌門還原脫鹵反應(yīng)機理的理解有利于環(huán)境中有機污染降解過程的進一步評估[44].

    放線菌門和厚壁菌門等也與 OM降解有密切關(guān)系.放線菌對環(huán)境中有機物和異形生物質(zhì)的降解具有重要作用[47],并且是一類在抗生素和酶制劑產(chǎn)業(yè)中具有重要應(yīng)用價值的微生物[48].條帶7與放線菌門的Nitriliruptor alkaliphilus有一定同源性(83%),Nitriliruptor alkaliphilus能夠分解脂肪腈類[49],在OR和NB沉積物中普遍存在,主要富集在 OR的表層(0~1cm)和深層,張勇等[47]的研究也表明,東太湖表層沉積物中放線菌群落多樣性較深層高.由于大多數(shù)腈類有毒且難降解,關(guān)于龍景湖沉積物中放線菌的研究,對于綠色化工和環(huán)境保護領(lǐng)域具有重要意義.與其他菌種的分布不同,厚壁菌門的條帶14主要富集在NB沉積物中,僅在 OR的表層(0~1cm)出現(xiàn).與條帶 14同 源 性 較 高(相 似 性 96%)的 是Saccharofermentans acetigenes,其生理學特征為專性厭氧,產(chǎn)芽孢,不運動,能夠通過葡萄糖發(fā)酵產(chǎn)生醋酸,但不能降解纖維素,在產(chǎn)甲烷環(huán)境的復(fù)雜OM降解中發(fā)揮重要作用[50].

    通過DGGE技術(shù)在龍景湖沉積物所提取的條帶中,大多數(shù)細菌為嚴格厭氧或兼性厭氧菌,與OM 的降解相關(guān).影響沉積物中細菌群落組成和分布的因素有很多,除了細菌本身的代謝類型,還有水動力作用[9]、有機質(zhì)的來源[33]及組成[24]、沉積物類型[19,51]、pH值、溶解氧(DO)濃度、氧化還原電位(Eh)[7,20]、鹽度[13,23]、營養(yǎng)鹽含量[1,12,23]等.不同沉積環(huán)境下的影響因子及作用程度并不相同.

    3.2 細菌群落與環(huán)境因子的響應(yīng)關(guān)系

    對龍景湖沉積物的細菌群落和理化指標進行了主成分分析(PCA),得到各菌種與環(huán)境變量的相關(guān)關(guān)系(圖5).環(huán)境變量與菌種的夾角<90°時說明二者呈正相關(guān),夾角越小關(guān)系越密切.第一主成分分析表明,細菌群落分布主要與TN相關(guān),其次是Eh.第二方向軸與TS和平均粒徑正相關(guān),而與孔隙水中總TOC濃度負相關(guān).由于理化指標較多,采用預(yù)選方式確定更適合的環(huán)境變量,用來代表細菌群落組成變化的驅(qū)動因子[52].只有環(huán)境變量顯著地提高模型的解釋能力(P<0.05)時,才能在預(yù)選中保留下來.預(yù)選過程中,細菌群落分布的變化主要受TN (P = 0.008)、TS (P = 0.001)、TOC (P = 0.019)和平均粒徑(P = 0.050)影響,分別解釋了細菌DGGE垂向分布中19%、20%、10%和8%的變化.結(jié)果表明細菌群落組成受上述4種環(huán)境因子的影響較大.

    龍景湖(OR和 NB)沉積物中的有機質(zhì)含量較高,尤其是表層,為好氧、厭氧、兼性厭氧細菌提供豐富的營養(yǎng)物質(zhì)和能源[6].本研究中通過DGGE技術(shù)所獲得的龍景湖沉積物中的細菌多與有機質(zhì)的降解有關(guān).Roske等[1]研究Saidenbach飲用水庫中的細菌多樣性和分布時發(fā)現(xiàn),氨濃度、含水率、OM含量是主要的驅(qū)動因子.Edlund等[7]的研究也表明,波羅的海沉積物中的細菌群落主要與有機碳相關(guān),其次受TN和 Eh的影響. Kristen等[20]的研究則表明,厭氧條件會顯著改變土壤微生物的群落組成.龍景湖 OR水深較 NB更深,沉積物中的有機質(zhì)、TN含量更高,δ-變形菌、擬桿菌、Ignavibacteriae等降解有機物的細菌種類和數(shù)量更多,通過有機質(zhì)的降解或形成的缺氧環(huán)境還會直接或間接促進了營養(yǎng)鹽的釋放(如正磷酸鹽和銨等)[6].

    圖5 細菌種類分布與沉積物環(huán)境因子的PCA(a)和RDA(b)分析Fig.5 PCA (a) and RDA (b) analyses between bacterial community and environmental factors of sediment

    對沉積物細菌群落和預(yù)選出的 4種理化指標進行冗余分析(RDA),由圖4可知,厚壁菌(條帶14)主要受沉積物中平均粒徑的影響.NB沉積物的平均粒徑較大,這可能是厚壁菌門普遍分布在NB沉積物中的原因.Colare s等[19]對Rhizophora mangle沉積物中細菌群落組成和環(huán)境變量的RDA分析也表明,沉積物的砂-黏土比例顯著影響沉積物中細菌和古菌的分布.其他可能影響龍景湖沉積物中細菌群落結(jié)構(gòu)和分布的因素為有機質(zhì)的來源[33]及組成[24]、沉積物類型[19,51]、pH值[23-24]等,會在接下來的研究中進一步分析.

    4 結(jié)論

    4.1 龍景湖原有河道(OR)沉積物的細菌群落多樣性較新形成庫灣(NB)的高.OR表層0~1cm沉積物中細菌多樣性最高,香農(nóng)指數(shù)為2.89.表層至深層香農(nóng)指數(shù)呈先下降后上升的趨勢,與豐度和均勻度變化規(guī)律相同.

    4.2 龍景湖沉積物中共發(fā)現(xiàn)7個門類的細菌,細菌群落分布隨深度差異明顯,OR表層沉積物以δ-變形菌、擬桿菌和綠彎菌為主,中層以綠彎菌為主,深層以δ-變形菌和綠彎菌為主;而NB的表層和中層的主要菌種都為δ-變形菌和綠彎菌. 4.3 通過DGGE技術(shù)所獲得的龍景湖沉積物中的細菌多與有機質(zhì)的降解有關(guān).通過 PCA和RDA表明細菌群落結(jié)構(gòu)和分布主要受沉積物中TN、TS含量、孔隙水中TOC濃度和平均粒徑的影響.OR沉積物中的有機質(zhì)、TN含量較NB高,δ-變形菌、擬桿菌、Ignavibacteriae等降解有機物的細菌種類和數(shù)量較多.平均粒徑可能是影響厚壁菌門主要分布在NB沉積物中的原因.

    [1]Roske K, Sachse R, Scheerer C, et al. Microbial diversity and composition of the sediment in the drinking water reservoir Saidenbach (Saxonia, Germany) [J]. Systematic and Applied Microbiology, 2012,35(1):35-44.

    [2]Wang Y, Sheng H F, He Y, et al. Comparison of the Levels of Bacterial Diversity in Freshwater, Intertidal Wetland, and Marine Sediments by Using Millions of Illumina Tags [J]. Applied and Environmental Microbiology, 2012,78(23):8264-8271.

    [3]鮑林林,陳永娟,王曉燕.北運河沉積物中氨氧化微生物的群落特征 [J]. 中國環(huán)境科學, 2015,35(1):179-189.

    [4]Urakawa H, Yoshida T, Nishimura M, et al. Characterization of depthrelated population variation in microbial communities of a coastal marine sediment using 16S rDNA-based approaches and quinone profiling [J]. Environmental Microbiology, 2000,2(5): 542-554.

    [5]Koizumi Y, Kojima H, Fukui M. Characterization of depthrelated microbial community structure in lake sediment by denaturing gradient gel electrophoresis of amplified 16S rDNA and reversely transcribed 16S rRNA fragments [J]. Fems Microbiol Ecol, 2003,46(2):147-157.

    [6]Sinkko H, Lukkari K, Sihvonen L M, et al. Bacteria Contribute to Sediment Nutrient Release and Reflect Progressed Eutrophication-Driven Hypoxia in an Organic-Rich Continental Sea [J]. Plos One, 2013,8(6):e67061.

    [7]Edlund A, Hardeman F, Jansson J K, et al. Active bacterial community structure along vertical redox gradients in Baltic Sea sediment [J]. Environmental Microbiology, 2008,10(8):2051-2063.

    [8]趙興青,楊柳燕,尹大強,等.不同空間位點沉積物理化性質(zhì)與微生物多樣性垂向分布規(guī)律 [J]. 環(huán)境科學, 2008,(12):3537-3545.

    [9]Boer S I, Hedtkamp S I C, van Beusekom J E E, et al. Time- and sediment depth-related variations in bacterial diversity and community structure in subtidal sands [J]. Isme Journal, 2009, 3(7):780-791.

    [10]劉愛菊,王洪海,潘嘉芬,等.孝婦河表層沉積物中重金屬賦存形態(tài)與微生物群落組成 [J]. 中國環(huán)境科學, 2010,30(8):1103-1109.

    [11]杜 萍,劉晶晶,沈李東,等.Biolog和PCR-DGGE技術(shù)解析椒江口沉積物微生物多樣性 [J]. 環(huán)境科學學報, 2012,32(6):1436-1444.

    [12]Liu Y, Zhang J X, Zhao L, et al. Spatial distribution of bacterial communities in high-altitude freshwater wetland sediment [J]. Limnology, 2014,15(3):249-256.

    [13]時 玉,孫懷博,劉勇勤,等.青藏高原淡水湖普莫雍錯和鹽水湖阿翁錯湖底沉積物中細菌群落的垂直分布 [J]. 微生物學通報, 2014,41(11):2379-2387.

    [14]向 興,王紅梅,龔林鋒,等.細菌群落在神農(nóng)架大九湖泥炭蘚與表層沉積物的垂向變化及其生態(tài)意義 [J]. 中國科學:地球科學, 2014,(6):1244-1252.

    [15]史春瀟,雷懷彥,趙 晶,等.南海北部九龍甲烷礁鄰區(qū)沉積物層中垂向細菌群落結(jié)構(gòu)特征研究 [J]. 沉積學報, 2014,32(6):1072-1082.

    [16]沈 爍,楊長明,成水平.PCR-DGGE分析合肥市塘西河表層沉積物細菌群落結(jié)構(gòu)空間分布特征 [J]. 應(yīng)用與環(huán)境生物學報, 2015,21(1):80-87.

    [17]Bukin S V, Pavlova O N, Manakov A Y, et al. The Ability of Microbial Community of Lake Baikal Bottom Sediments Associated with Gas Discharge to Carry Out the Transformation of Organic Matter under Thermobaric Conditions [J]. Frontiers in Microbiology, 2016,7: Article 690.

    [18]Shao K, Gao G, Wang Y, et al. Vertical diversity of sediment bacterial communities in two different trophic states of the eutrophic Lake Taihu, China [J]. Journal of Environmental Sciences-China, 2013,25(6):1186-1194.

    [19]Colares G B, Melo V M M. Relating microbial community structure and environmental variables in mangrove sediments inside Rhizophora mangle L. habitats [J]. Applied Soil Ecology, 2013,64:171-177.

    [20]DeAngelis K M, Silver W L, Thompson A W, et al. Microbial communities acclimate to recurring changes in soil redox potential status [J]. Environmental Microbiology, 2010,12(12):3137-3149.

    [21]Davis M L, Cornwell D A.環(huán)境工程導(dǎo)論 [M]. 清華大學出版社, 2010.

    [22]Koizumi Y, Kojima H, Oguri K, et al. Vertical and temporal shifts in microbial communities in the water column and sediment of saline meromictic Lake Kaiike (Japan), as determined by a 16S rDNA-based analysis, and related to physicochemical gradients [J]. Environmental Microbiology, 2004,6(6):622-637.

    [23]Xiong J B, Liu Y Q, Lin XG, et al. Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau [J]. Environmental Microbiology, 2012,14(9):2457-2466.

    [24]Oni O E, Schmidt F, Miyatake T, et al. Microbial Communities and Organic Matter Composition in Surface and Subsurface Sediments of the Helgoland Mud Area, North Sea [J]. Frontiers in Microbiology, 2015,6:Article 1290.

    [25]De Wever A, Muylaert K, Van der Gucht K, et al. Bacterial community composition in Lake Tanganyika: vertical and horizontal heterogeneity [J]. Applied and environmental microbiology, 2005,71(9):5029-5037.

    [26]吉芳英,顏海波,何 強,等.龍景湖龍景溝匯水區(qū)沉積物-水界面氮形態(tài)空間分布特征 [J]. 中國環(huán)境科學, 2015,35(10):3101-3107.

    [27]田 濤,張代鈞,李玉蓮,等.重慶園博園龍景湖水體硫酸鹽還原及氮化物和TOC的影響[C]//2014中國環(huán)境科學學會學術(shù)年會論文集. 2014:3688-3695.

    [28]張甘林.山地城市景觀深水湖泊湖灣人工強化復(fù)氧技術(shù)示范研究 [D]. 重慶:重慶大學, 2014.

    [29]潘延安,雷 沛,張 洪,等.重慶園博園龍景湖新建初期內(nèi)源氮磷分布特征及擴散通量估算 [J]. 環(huán)境科學, 2014,35(5):1727-1734.

    [30]Hu A Y, Yang X Y, Chen N W, et al. Response of bacterial communities to environmental changes in a mesoscale subtropical watershed, Southeast China [J]. Science of the Total Environment, 2014,472:746-756.

    [31]Pett-Ridge J, Firestone M K. Redox fluctuation structures microbial communities in a wet tropical soil [J]. Applied &Environmental Microbiology, 2005,71(71):6998-7007.

    [32]趙興青,楊柳燕,尹大強,等.太湖沉積物中微生物多樣性垂向分布特征 [J]. 地學前緣, 2008,15(6):177-184.

    [33]Xiong W, Xie P, Wang SR, et al. Sources of organic matter affect depth-related microbial community composition in sediments of Lake Erhai, Southwest China [J]. Journal of Limnology, 2015, 74(2):310-323.

    [34]Thomas S H, Wagner R D, Arakaki A K, et al. The mosaic genome of Anaeromyxobacter dehalogenans strain 2CP-C suggests an aerobic common ancestor to the delta-proteobacteria [J]. PLoS One, 2008,3(5):e2103.

    [35]Qiu Y L, Hanada S, Ohashi A, et al. Syntrophorhabdus aromaticivorans gen. nov., sp. nov., the first cultured anaerobe capable of degrading phenol to acetate in obligate syntrophic associations with a hydrogenotrophic methanogen [J]. Applied and environmental microbiology, 2008,74(7):2051-2058.

    [36]Liu Y, Balkwill D L, Aldrich H C, et al. Characterization of the anaerobic propionate-degrading syntrophs Smithella propionica gen. nov., sp. nov. and Syntrophobacter wolinii [J]. Int J Syst Bacteriol, 1999,49Pt 2:545-556.

    [37]Hahn M W, Schauer M. 'Candidatus Aquirestis calciphila' and 'Candidatus Haliscomenobacter calcifugiens', filamentous, planktonic bacteria inhabiting natural lakes [J]. Int J Syst Evol Microbiol, 2007,57(Pt5):936-940.

    [38]Gomez-Pereira P R, Schuler M, Fuchs B M, et al. Genomic content of uncultured Bacteroidetes from contrasting oceanic provinces in the North Atlantic Ocean [J]. Environmental Microbiology, 2012,14(1):52-66.

    [39]Zhang K, Tang Y, Zhang L, et al. Parasegetibacter luojiensis gen. nov., sp. nov., a member of the phylum Bacteroidetes isolated from a forest soil [J]. Int J Syst Evol Microbiol, 2009,59(Pt 12):3058-3062.

    [40]Podosokorskaya O A, Kadnikov V V, Gavrilov S N, et al. Characterization of Melioribacter roseus gen. nov., sp. nov., a novel facultatively anaerobic thermophilic cellulolytic bacterium from the class Ignavibacteria, and a proposal of a novel bacterial phylum Ignavibacteriae [J]. Environmental Microbiology, 2013, 15(6):1759-1771.

    [41]Imhoff JF. Phylogenetic taxonomy of the family Chlorobiaceae on the basis of 16S rRNA and fmo (Fenna-Matthews-Olson protein) gene sequences [J]. Int J Syst Evol Microbiol, 2003,53(Pt 4):941-951.

    [42]Iino T, Mori K, Uchino Y, et al. Ignavibacterium album gen. nov., sp. nov., a moderately thermophilic anaerobic bacterium isolated from microbial mats at a terrestrial hot spring and proposal of Ignavibacteria classis nov., for a novel lineage at the periphery of green sulfur bacteria [J]. Int J Syst Evol Microbiol, 2010,60(Pt 6):1376-1382.

    [43]Inagaki F, Nunoura T, Nakagawa S, et al. Biogeographical distribution and diversity of microbes in methane hydrate-bearing deep marine sediments, on the Pacific Ocean Margin [J]. Proceedings of the National Academy of Sciences of the United States of America, 2006,103(8):2815-2820.

    [44]Martin-Gonzalez L, Mortan S H, Rosell M, et al. Stable Carbon Isotope Fractionation During 1,2-Dichloropropane-to-Propene Transformation by an Enrichment Culture Containing Dehalogenimonas Strains and a dcpA Gene [J]. Environmental Science & Technology, 2015,49(14):8666-8674.

    [45]Hug L A, Maphosa F, Leys D, et al. Overview of organohaliderespiring bacteria and a proposal for a classification system for reductive dehalogenases [J]. Philos T R Soc B, 2013,368(1616):68-76.

    [46]Kadnikov V V, Mardanov A V, Beletsky A V, et al. Microbial community structure in methane hydrate-bearing sediments of freshwater Lake Baikal [J]. Fems Microbiol Ecol, 2012,79(2):348-358.

    [47]張 勇,沈 吉,王建軍,等.東太湖表層沉積物放線菌群落結(jié)構(gòu)多樣性及其空間異質(zhì)性 [J]. 生態(tài)科學, 2011,30(1):8-13.

    [48]Takahashi Y, Omura S. Isolation of new actinomycete strains for the screening of new bioactive compounds [J]. Journal of General & Applied Microbiology, 2003,49(3):141-154.

    [49]Sorokin D Y, van Pelt S, Tourova TP, et al. Nitriliruptor alkaliphilus gen. nov., sp. nov., a deep-lineage haloalkaliphilic actinobacterium from soda lakes capable of growth on aliphatic nitriles, and proposal of Nitriliruptoraceae fam. nov. and Nitriliruptorales ord. nov [J]. Int J Syst Evol Microbiol, 2009, 59(Pt 2):248-253.

    [50]Chen S, Niu L, Zhang Y. Saccharofermentans acetigenes gen. nov., sp. nov., an anaerobic bacterium isolated from sludge treating brewery wastewater [J]. Int J Syst Evol Microbiol, 2010,60(Pt 12):2735-2738.

    [51]Jackson C R, Weeks A Q. Influence of particle size on bacterial community structure in aquatic sediments as revealed by 16S rRNA gene sequence analysis [J]. Applied and Environmental Microbiology, 2008,74(16):5237-5240.

    [52]Wang X Q, Cui H Y, Shi J H, et al. Relationship between bacterial diversity and environmental parameters during composting of different raw materials [J]. Bioresource Technology, 2015,198:395-402.

    Vertical distribution of bacterial communities in sediments of Longjing Lake


    NIU Feng-xia, JI Fang-ying*, ZHAOGen, ZHANG Qian, SHEN Qiu-shi, HE Qiang, YAN Hai-bo (Key Laboratory of Three Gorges Reservoir Region’s Eco-Environment, Ministry of Education, National Centre for International Research of Low-carbon and Green Buildings, Chongqing University, Chongqing 400045, China). China Environmental Science, 2017,37(6):2322~2331

    This study explored the vertical distribution of bacterial communities in sediments under different deposition time and its role in the migration and transformation of substances in lake ecosystem. Stratified sediment samples were collected respectively from Original River (OR) and New-formed Bay (NB) in Longjing Lake. Redundancy analysis (RDA) were conducted based on Denaturing Gradient Gel Electrophoresis (DGGE) technology and environmental parameters. Results showed that Shannon-Wiener indexes (H) and Richness indexes (S) in OR were higher than those in NB and vertical change trend of H, S and Evenness (E) indexes decreased first and then increased as depth increased. Phylogenetic analysis revealed that Longjing lake mainly included seven phyla of bacteria. Deltaproteobacteria, Bacteroidetes and Chloroflexi dominated in surface layers whereas Chloroflexi dominated in middle layers of OR sediments. Deltaproteobacteria and Chloroflexi were the main phyla in deep layers of both sites. Other phyla included Ignavibacteriae, Actinobacteria, Firmicutes and Chlorobi. RDA showed that the main factors influencing bacterial community composition and distribution were total nitrogen (TN), total sulphur (TS) and the concentrations of total organic carbon (TOC) in pore water, as well as average particle size in sediments. Otherwise, the obtained bacteria were closely related to the degradation of organic matter (OM). That TOC and TN contents in OR were higher than those in NB caused relatively higher bacteria species and abundance of Deltaproteobacteria, Bacteroidetes, and Ignavibacteriae.

    bacterial community;vertical distribution;DGGE;RDA

    X171.5

    A

    1000-6923(2017)06-2322-10

    牛鳳霞(1989-),女,河南開封人,重慶大學博士研究生,主要從事水污染控制研究.發(fā)表論文3篇.

    2016-09-23

    國家水體污染控制與治理科技重大專項(2012ZX07307-001);重慶大學 2015年研究生科研創(chuàng)新項目基金(CYB15041);中央高校基本業(yè)務(wù)費項目(2015CDJXY)

    * 責任作者, 教授, jfy@cqu.edu.cn

    猜你喜歡
    中層條帶深層
    考慮各向異性滲流的重力壩深層抗滑穩(wěn)定分析
    雙麥克斯韋分布下極區(qū)中層塵埃粒子帶電研究
    SAM系統(tǒng)對TDCS數(shù)據(jù)的優(yōu)化處理與深層應(yīng)用
    基于條帶模式GEOSAR-TOPS模式UAVSAR的雙基成像算法
    對“醫(yī)患失去信任”的深層憂慮
    基于 Savitzky-Golay 加權(quán)擬合的紅外圖像非均勻性條帶校正方法
    中國光學(2015年1期)2015-06-06 18:30:20
    企業(yè)中層管理團隊在運營中的困境與作用
    電視節(jié)目低俗化的深層反思
    一種基于MATLAB的聲吶條帶圖像自動拼接算法
    海岸工程(2014年4期)2014-02-27 12:51:28
    從組織支持談中層管理者激勵*
    国产熟女xx| а√天堂www在线а√下载| 亚洲最大成人av| 我要看日韩黄色一级片| 啪啪无遮挡十八禁网站| 亚洲最大成人手机在线| 51国产日韩欧美| 国产精品美女特级片免费视频播放器| 别揉我奶头~嗯~啊~动态视频| 美女 人体艺术 gogo| 嫩草影院入口| 亚洲av电影不卡..在线观看| 18+在线观看网站| 欧美区成人在线视频| 久久欧美精品欧美久久欧美| 久久99热这里只有精品18| 欧美一区二区国产精品久久精品| 免费人成视频x8x8入口观看| 亚洲人成网站在线播| 99精品在免费线老司机午夜| 小说图片视频综合网站| 成人av一区二区三区在线看| 最新中文字幕久久久久| 90打野战视频偷拍视频| 午夜激情欧美在线| 成人一区二区视频在线观看| 黄色一级大片看看| 一进一出抽搐gif免费好疼| 日韩欧美精品v在线| 国产精品久久久久久精品电影| av欧美777| 亚洲性夜色夜夜综合| 97热精品久久久久久| aaaaa片日本免费| 精品人妻一区二区三区麻豆 | 自拍偷自拍亚洲精品老妇| 五月玫瑰六月丁香| 老鸭窝网址在线观看| 国产精品美女特级片免费视频播放器| 精品久久久久久久久久久久久| 看黄色毛片网站| 免费无遮挡裸体视频| 91字幕亚洲| 99久久精品热视频| 欧美性猛交╳xxx乱大交人| 一进一出抽搐动态| 99热精品在线国产| 在线天堂最新版资源| 色在线成人网| 国产黄a三级三级三级人| 久久久久久大精品| 18美女黄网站色大片免费观看| 欧美极品一区二区三区四区| 日韩欧美在线乱码| 波野结衣二区三区在线| 少妇被粗大猛烈的视频| 露出奶头的视频| 亚洲国产欧美人成| 久久午夜福利片| 日本撒尿小便嘘嘘汇集6| 日韩欧美在线二视频| 亚洲色图av天堂| 大型黄色视频在线免费观看| 色吧在线观看| 国产精品久久视频播放| 99久国产av精品| xxxwww97欧美| 亚洲一区二区三区不卡视频| 日韩欧美国产一区二区入口| 亚洲精品粉嫩美女一区| 男插女下体视频免费在线播放| 色视频www国产| 中文字幕人妻熟人妻熟丝袜美| 午夜免费男女啪啪视频观看 | 丰满的人妻完整版| 又黄又爽又免费观看的视频| 丁香六月欧美| 国产成人啪精品午夜网站| 国产伦精品一区二区三区四那| 两性午夜刺激爽爽歪歪视频在线观看| 12—13女人毛片做爰片一| 久久精品影院6| 欧美成人a在线观看| 色精品久久人妻99蜜桃| 久久国产乱子免费精品| 级片在线观看| 99精品久久久久人妻精品| 色播亚洲综合网| 性色av乱码一区二区三区2| 国产亚洲av嫩草精品影院| 亚洲第一欧美日韩一区二区三区| 搡老熟女国产l中国老女人| 91久久精品电影网| 精品一区二区免费观看| 18美女黄网站色大片免费观看| 级片在线观看| 亚洲精品亚洲一区二区| 日本a在线网址| 特大巨黑吊av在线直播| 欧美最新免费一区二区三区 | 国产不卡一卡二| 国产精品一区二区三区四区久久| 亚洲五月婷婷丁香| 别揉我奶头 嗯啊视频| 欧美乱色亚洲激情| 此物有八面人人有两片| 国产精品三级大全| 51国产日韩欧美| 亚洲三级黄色毛片| 在线天堂最新版资源| 男人舔女人下体高潮全视频| 最后的刺客免费高清国语| 99热这里只有精品一区| 国内精品一区二区在线观看| eeuss影院久久| 国产精品自产拍在线观看55亚洲| 成人性生交大片免费视频hd| 亚洲在线观看片| 免费在线观看影片大全网站| 国产一区二区在线av高清观看| 一区二区三区免费毛片| 精品人妻一区二区三区麻豆 | 在线观看美女被高潮喷水网站 | 欧美日韩国产亚洲二区| av中文乱码字幕在线| 国产精品一区二区三区四区久久| 麻豆国产av国片精品| 亚洲人成网站高清观看| 全区人妻精品视频| 精品午夜福利在线看| 性色avwww在线观看| 国产精品久久久久久久久免 | 黄色女人牲交| 男女视频在线观看网站免费| 欧美国产日韩亚洲一区| 黄色配什么色好看| 男女下面进入的视频免费午夜| 性插视频无遮挡在线免费观看| 精品日产1卡2卡| 日韩国内少妇激情av| 免费看a级黄色片| 床上黄色一级片| 毛片女人毛片| 又黄又爽又刺激的免费视频.| 欧美绝顶高潮抽搐喷水| 国产私拍福利视频在线观看| 国产淫片久久久久久久久 | 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久,| 五月玫瑰六月丁香| 亚洲成av人片免费观看| 校园春色视频在线观看| 国产av一区在线观看免费| 国产激情偷乱视频一区二区| 人人妻,人人澡人人爽秒播| 国产黄a三级三级三级人| 国产欧美日韩一区二区精品| 亚洲在线观看片| 欧美日韩综合久久久久久 | 变态另类成人亚洲欧美熟女| 在线观看66精品国产| 欧美在线黄色| 久久午夜福利片| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩东京热| 美女大奶头视频| 一夜夜www| 国产精品伦人一区二区| 激情在线观看视频在线高清| 我的老师免费观看完整版| 变态另类丝袜制服| 别揉我奶头~嗯~啊~动态视频| 免费人成在线观看视频色| 亚洲精品在线美女| 亚洲内射少妇av| 欧美黑人巨大hd| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 亚洲精品色激情综合| 亚洲欧美清纯卡通| 丰满人妻熟妇乱又伦精品不卡| 亚洲三级黄色毛片| 看片在线看免费视频| 99riav亚洲国产免费| 国模一区二区三区四区视频| 免费搜索国产男女视频| 亚洲av免费高清在线观看| 男插女下体视频免费在线播放| 最近在线观看免费完整版| 一边摸一边抽搐一进一小说| 免费搜索国产男女视频| 中国美女看黄片| 欧美成狂野欧美在线观看| 亚洲成a人片在线一区二区| www.熟女人妻精品国产| 亚洲在线观看片| 69av精品久久久久久| 亚洲黑人精品在线| 国模一区二区三区四区视频| 久久6这里有精品| 久久久久九九精品影院| 精品人妻熟女av久视频| 一本精品99久久精品77| av在线蜜桃| 国产欧美日韩精品亚洲av| 国产美女午夜福利| 欧美激情久久久久久爽电影| 国产av在哪里看| 国产精品女同一区二区软件 | 精品一区二区三区视频在线| 成人欧美大片| 18禁黄网站禁片免费观看直播| 亚洲av一区综合| 一级作爱视频免费观看| 在线观看免费视频日本深夜| 国产精品三级大全| 又爽又黄无遮挡网站| 日韩精品中文字幕看吧| 高清日韩中文字幕在线| 色视频www国产| 伦理电影大哥的女人| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 小蜜桃在线观看免费完整版高清| 男人舔女人下体高潮全视频| xxxwww97欧美| 一级黄色大片毛片| 欧美一区二区国产精品久久精品| 日本黄色片子视频| 欧美三级亚洲精品| 国产一区二区亚洲精品在线观看| 亚洲精品456在线播放app | 伊人久久精品亚洲午夜| 在线观看一区二区三区| 欧美激情在线99| 亚洲欧美激情综合另类| 看十八女毛片水多多多| 久久中文看片网| 非洲黑人性xxxx精品又粗又长| 午夜日韩欧美国产| av国产免费在线观看| av女优亚洲男人天堂| 久久久久久久久中文| 少妇人妻一区二区三区视频| 午夜老司机福利剧场| 乱码一卡2卡4卡精品| 国产激情偷乱视频一区二区| 亚洲成人久久性| 欧美日韩中文字幕国产精品一区二区三区| 亚洲天堂国产精品一区在线| 少妇被粗大猛烈的视频| 男人狂女人下面高潮的视频| 神马国产精品三级电影在线观看| 麻豆国产av国片精品| 婷婷精品国产亚洲av| 99国产综合亚洲精品| 高潮久久久久久久久久久不卡| avwww免费| 99热6这里只有精品| 亚洲国产高清在线一区二区三| 搡老妇女老女人老熟妇| 国产亚洲精品av在线| 精品午夜福利在线看| 欧美高清性xxxxhd video| av在线观看视频网站免费| 欧美一区二区精品小视频在线| 99久久无色码亚洲精品果冻| 国产午夜精品论理片| 久久精品影院6| 亚洲人成网站在线播| av中文乱码字幕在线| 1024手机看黄色片| 日本黄色视频三级网站网址| 男女视频在线观看网站免费| 成人国产综合亚洲| 精品免费久久久久久久清纯| 成人永久免费在线观看视频| 90打野战视频偷拍视频| 99热精品在线国产| 欧美精品国产亚洲| 亚洲一区高清亚洲精品| 亚洲经典国产精华液单 | 全区人妻精品视频| 欧美bdsm另类| 成人无遮挡网站| 婷婷丁香在线五月| 成年版毛片免费区| 日韩成人在线观看一区二区三区| 女人十人毛片免费观看3o分钟| 首页视频小说图片口味搜索| 免费一级毛片在线播放高清视频| 亚洲精品在线美女| 亚洲av美国av| 午夜老司机福利剧场| 久久人妻av系列| 精品国产亚洲在线| 国产精品久久久久久精品电影| 亚洲精品久久国产高清桃花| 免费观看精品视频网站| 亚洲乱码一区二区免费版| 日韩av在线大香蕉| 首页视频小说图片口味搜索| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添av毛片 | 国产主播在线观看一区二区| 一二三四社区在线视频社区8| 欧美成人一区二区免费高清观看| 国产美女午夜福利| av在线老鸭窝| 欧美日韩黄片免| 国产精品人妻久久久久久| 身体一侧抽搐| 草草在线视频免费看| 天天躁日日操中文字幕| 亚洲精品456在线播放app | 精品久久久久久久人妻蜜臀av| 亚洲专区中文字幕在线| 久久6这里有精品| 看片在线看免费视频| 婷婷色综合大香蕉| 长腿黑丝高跟| 三级国产精品欧美在线观看| 亚洲,欧美精品.| 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 国产单亲对白刺激| 成人美女网站在线观看视频| 在线天堂最新版资源| 国产精品久久久久久精品电影| 国产人妻一区二区三区在| 99国产精品一区二区蜜桃av| 91麻豆av在线| 亚洲av中文字字幕乱码综合| 国内少妇人妻偷人精品xxx网站| 十八禁网站免费在线| av在线天堂中文字幕| 热99re8久久精品国产| av专区在线播放| 久久国产精品影院| 久久精品久久久久久噜噜老黄 | 久久久久九九精品影院| 免费在线观看成人毛片| 日日夜夜操网爽| av在线观看视频网站免费| 国产精品精品国产色婷婷| 国产v大片淫在线免费观看| 国内精品一区二区在线观看| 国产蜜桃级精品一区二区三区| 97超视频在线观看视频| 麻豆国产av国片精品| 久久精品国产亚洲av涩爱 | 最后的刺客免费高清国语| 看黄色毛片网站| 好男人电影高清在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲aⅴ乱码一区二区在线播放| 如何舔出高潮| 真实男女啪啪啪动态图| 99久久精品一区二区三区| 亚洲精品在线美女| 亚洲av电影在线进入| 很黄的视频免费| 尤物成人国产欧美一区二区三区| 国产一级毛片七仙女欲春2| 老熟妇乱子伦视频在线观看| 久久欧美精品欧美久久欧美| 午夜a级毛片| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点| 亚洲乱码一区二区免费版| 一夜夜www| 直男gayav资源| 我的老师免费观看完整版| 久久精品综合一区二区三区| 精品无人区乱码1区二区| 国产成+人综合+亚洲专区| 婷婷色综合大香蕉| 欧美zozozo另类| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| АⅤ资源中文在线天堂| 欧美极品一区二区三区四区| .国产精品久久| 天堂√8在线中文| 亚洲第一欧美日韩一区二区三区| 亚洲中文日韩欧美视频| 2021天堂中文幕一二区在线观| 91狼人影院| 国产高潮美女av| 嫩草影院精品99| 特大巨黑吊av在线直播| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区在线臀色熟女| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 亚洲自偷自拍三级| 99riav亚洲国产免费| 欧美一级a爱片免费观看看| 成年女人看的毛片在线观看| 美女高潮喷水抽搐中文字幕| 99热6这里只有精品| 国产毛片a区久久久久| 国产国拍精品亚洲av在线观看| 国产精品一区二区免费欧美| 国产不卡一卡二| 国产单亲对白刺激| 此物有八面人人有两片| av女优亚洲男人天堂| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| 日韩国内少妇激情av| 久久久久久久午夜电影| 一本精品99久久精品77| 国产午夜精品久久久久久一区二区三区 | 亚洲熟妇中文字幕五十中出| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 最近在线观看免费完整版| 精品国产亚洲在线| 99riav亚洲国产免费| 亚洲精品一卡2卡三卡4卡5卡| 噜噜噜噜噜久久久久久91| 丰满人妻熟妇乱又伦精品不卡| 国产男靠女视频免费网站| 给我免费播放毛片高清在线观看| 国产成+人综合+亚洲专区| 一a级毛片在线观看| 综合色av麻豆| 99热这里只有是精品50| 每晚都被弄得嗷嗷叫到高潮| 亚洲av电影在线进入| 两人在一起打扑克的视频| 亚洲第一区二区三区不卡| 哪里可以看免费的av片| 国产成人aa在线观看| 少妇的逼水好多| 一区福利在线观看| 99riav亚洲国产免费| 欧美成人免费av一区二区三区| 特级一级黄色大片| 一级黄色大片毛片| 99视频精品全部免费 在线| 在线观看一区二区三区| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 综合色av麻豆| 最新在线观看一区二区三区| 欧美区成人在线视频| 天堂√8在线中文| 国产又黄又爽又无遮挡在线| 国产三级中文精品| 久久中文看片网| 99久久成人亚洲精品观看| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 亚洲人成伊人成综合网2020| 亚洲国产日韩欧美精品在线观看| 久久国产乱子免费精品| or卡值多少钱| 日韩有码中文字幕| av黄色大香蕉| 可以在线观看毛片的网站| 亚洲av成人不卡在线观看播放网| 搡老妇女老女人老熟妇| 高清在线国产一区| 人妻夜夜爽99麻豆av| 国产精品一区二区三区四区久久| 高清毛片免费观看视频网站| www.熟女人妻精品国产| 精品国产亚洲在线| a级毛片a级免费在线| 如何舔出高潮| 三级男女做爰猛烈吃奶摸视频| 美女高潮喷水抽搐中文字幕| 禁无遮挡网站| 可以在线观看毛片的网站| 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成人久久爱视频| 精品国内亚洲2022精品成人| 国产三级中文精品| 国产一区二区三区在线臀色熟女| 99久久精品国产亚洲精品| 在线观看舔阴道视频| 日韩中字成人| 噜噜噜噜噜久久久久久91| 亚洲欧美精品综合久久99| 国产精品野战在线观看| 婷婷精品国产亚洲av在线| 成人午夜高清在线视频| 亚洲av电影不卡..在线观看| 免费av观看视频| 国产主播在线观看一区二区| 简卡轻食公司| 特级一级黄色大片| 色吧在线观看| 国产精华一区二区三区| 国产欧美日韩精品一区二区| 亚洲精品在线观看二区| 欧美成人性av电影在线观看| 亚洲久久久久久中文字幕| 夜夜看夜夜爽夜夜摸| 琪琪午夜伦伦电影理论片6080| 热99re8久久精品国产| 亚洲一区二区三区不卡视频| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 内射极品少妇av片p| 久久久久久九九精品二区国产| 嫩草影院精品99| 99久国产av精品| 99国产精品一区二区三区| 亚洲综合色惰| 欧美色欧美亚洲另类二区| 成人永久免费在线观看视频| 亚洲18禁久久av| 热99re8久久精品国产| 757午夜福利合集在线观看| 黄片小视频在线播放| 两人在一起打扑克的视频| 亚洲第一区二区三区不卡| 性色av乱码一区二区三区2| 国产精品野战在线观看| 大型黄色视频在线免费观看| 久久久精品欧美日韩精品| 51国产日韩欧美| 欧美不卡视频在线免费观看| 在线免费观看的www视频| 在线免费观看不下载黄p国产 | 美女黄网站色视频| 老司机午夜十八禁免费视频| 国产老妇女一区| 观看美女的网站| 脱女人内裤的视频| 久久精品国产清高在天天线| 桃色一区二区三区在线观看| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频| 日日夜夜操网爽| 国产精品av视频在线免费观看| 国产欧美日韩一区二区精品| 变态另类成人亚洲欧美熟女| 中文字幕熟女人妻在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲精品亚洲一区二区| 欧美成人一区二区免费高清观看| 亚洲一区二区三区不卡视频| 日本精品一区二区三区蜜桃| 亚洲成a人片在线一区二区| 国产免费一级a男人的天堂| а√天堂www在线а√下载| 欧美中文日本在线观看视频| 夜夜看夜夜爽夜夜摸| av视频在线观看入口| 一级黄色大片毛片| 欧美成狂野欧美在线观看| 亚洲成人免费电影在线观看| 国产精品亚洲av一区麻豆| 久久精品国产自在天天线| 午夜福利成人在线免费观看| 国产成人福利小说| 日日摸夜夜添夜夜添av毛片 | 国产精品,欧美在线| 日韩欧美在线乱码| 观看美女的网站| 网址你懂的国产日韩在线| 亚洲黑人精品在线| 亚洲中文字幕日韩| 亚洲av免费高清在线观看| 内射极品少妇av片p| 1024手机看黄色片| 此物有八面人人有两片| 亚洲人成伊人成综合网2020| 黄片小视频在线播放| 成人精品一区二区免费| 毛片一级片免费看久久久久 | 亚洲人成伊人成综合网2020| 国产精品一区二区性色av| 亚洲18禁久久av| 看黄色毛片网站| 男女之事视频高清在线观看| 一本综合久久免费| 麻豆成人午夜福利视频| 国产av在哪里看| 久久99热6这里只有精品| av欧美777| 国产aⅴ精品一区二区三区波| 长腿黑丝高跟| 精品国内亚洲2022精品成人| 国产成人福利小说| 亚洲国产日韩欧美精品在线观看| 亚洲欧美清纯卡通| 91字幕亚洲| 亚洲精品456在线播放app | 亚洲内射少妇av| 欧美最黄视频在线播放免费| 国产精品久久视频播放| 人妻夜夜爽99麻豆av| 51国产日韩欧美| 黄片小视频在线播放| 国产精品电影一区二区三区| 欧美性猛交黑人性爽| 国产免费男女视频| 久久天躁狠狠躁夜夜2o2o| 97碰自拍视频| 亚洲av电影不卡..在线观看| 亚洲,欧美精品.| 91在线观看av| 亚洲精华国产精华精| 三级男女做爰猛烈吃奶摸视频| 麻豆国产av国片精品| 少妇人妻一区二区三区视频| 国产aⅴ精品一区二区三区波| 国产精品,欧美在线| 国产精品乱码一区二三区的特点|