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

    鄱陽湖不同區(qū)域沉積物細菌群落結(jié)構(gòu)、功能變化及其與環(huán)境因子的關(guān)系

    2017-06-05 14:21:21寇文伯于昊天劉亞軍馬燕天
    環(huán)境科學研究 2017年4期
    關(guān)鍵詞:結(jié)構(gòu)功能

    金 笑, 寇文伯, 于昊天, 劉亞軍, 馬燕天, 吳 蘭

    南昌大學生命科學學院, 鄱陽湖環(huán)境與資源利用教育部重點實驗室, 江西 南昌 330031

    鄱陽湖不同區(qū)域沉積物細菌群落結(jié)構(gòu)、功能變化及其與環(huán)境因子的關(guān)系

    金 笑, 寇文伯, 于昊天, 劉亞軍, 馬燕天, 吳 蘭*

    南昌大學生命科學學院, 鄱陽湖環(huán)境與資源利用教育部重點實驗室, 江西 南昌 330031

    為了解鄱陽湖湖區(qū)內(nèi)具有不同形態(tài)、功能的沉積物微生物群落的空間差異,采用T- RFLP(末端限制性片段長度多態(tài)性)和Biolog- Eco微平板法,對松門山(主湖區(qū))、饒河(入湖口)、南磯山(碟形湖)和白沙洲(湖汊)樣地表層沉積物進行細菌遺傳多樣性和微生物代謝功能多樣性進行分析.結(jié)果表明:鄱陽湖具有不同形態(tài)、功能的樣地間細菌群落結(jié)構(gòu)及微生物功能多樣性存在協(xié)同的空間異質(zhì)性,顯示饒河樣地與其他樣地之間差異顯著(P<0.05).基于T- RFLP結(jié)果,饒河樣地與其他樣地細菌群落結(jié)構(gòu)的差異主要來源于中等相對豐度(0.20%~2.00%)的TRFs片段.此外,饒河樣地微生物代謝活性顯著低于其他樣地(P<0.05),主要是因為饒河樣地微生物對多聚物類和碳水化合物類碳源的利用能力顯著降低.基于微生物群落與環(huán)境因子的冗余分析(RDA)結(jié)果顯示,營養(yǎng)物質(zhì)(有機質(zhì)、TOC、TN)和重金屬(Zn)是影響細菌群落結(jié)構(gòu)的主要因子,而水深、pH、速效氮〔w(NH4+-N)、w(NO3--N)〕、重金屬〔w(Zn)、w(Cu)、w(Cd)〕是影響微生物代謝多樣性的顯著因子.研究顯示,重金屬(Zn)對鄱陽湖沉積物細菌群落結(jié)構(gòu)及微生物功能多樣性具有顯著影響.研究結(jié)果可為探索大型淡水湖泊底棲微生物群落結(jié)構(gòu)和功能的空間異質(zhì)性以及探討微生物群落結(jié)構(gòu)與功能對環(huán)境因素的響應(yīng)提供理論依據(jù).

    鄱陽湖; Biolog-Eco微平板法; T-RFLP分析; 微生物代謝功能; 細菌群落結(jié)構(gòu)

    近年來,隨著對淡水資源需求量的增加以及對水資源的不合理利用,淡水湖泊的水環(huán)境安全日益成為國內(nèi)外學者十分關(guān)注的環(huán)境問題[1- 2].在淡水湖泊生態(tài)系統(tǒng)中,表層沉積物是微生物聚居和生源物質(zhì)賦存的主要介質(zhì),而微生物可以通過同化、異化等代謝過程來影響沉積物中營養(yǎng)鹽及重金屬的分布、轉(zhuǎn)化、利用等[3- 5].因此,湖泊沉積物微生物群落可作為水環(huán)境變化和演替的重要標志.研究淡水湖泊沉積物中的微生物群落對環(huán)境因子的響應(yīng)規(guī)律,并揭示其結(jié)構(gòu)與功能的聯(lián)系,對了解淡水湖泊水環(huán)境中生源要素的循環(huán)具有重要作用.

    目前,有諸多學者進行了依賴于地理距離的沉積物細菌群落差異比較,這種差異主要受到了氣候、地質(zhì)背景、水文特征以及來自于湖體及集水區(qū)人為因素的影響[6- 8],水體環(huán)境的變化和外來物質(zhì)的輸入都會影響沉積物中的細菌群落結(jié)構(gòu)[9].如Macalady等[10]對美國加利福尼亞州受到無機汞污染的克里爾萊克湖中3個不同區(qū)域進行研究,發(fā)現(xiàn)其沉積物細菌群落結(jié)構(gòu)與采樣區(qū)域密切相關(guān),在高汞區(qū)域演替出的類脫硫桿菌微生物在沉積物汞甲基化過程中發(fā)揮了重要作用.此外,研究表明,營養(yǎng)物質(zhì)會對微生物代謝功能活性起到促進作用[11],而在重金屬污染區(qū)域進行的研究中,重金屬濃度升高對微生物代謝功能活性的作用呈現(xiàn)先促進后抑制的規(guī)律[12].然而,以往對細菌群落的研究大多是建立在單一類型的環(huán)境因素對微生物群落的影響方面,對于具有復(fù)雜性生境的大型淡水湖泊來說,單一類型環(huán)境因素的改變并不足以反映微生物對環(huán)境的響應(yīng)特征,因而其研究還遠遠不足.除此之外,依據(jù)杜萍等的研究[13- 15],在環(huán)境改變的條件下,微生物群落結(jié)構(gòu)與功能的關(guān)系分別呈現(xiàn)協(xié)同、功能優(yōu)先于結(jié)構(gòu)以及結(jié)構(gòu)優(yōu)先于功能變化的特征.由此可見,湖泊底泥微生物群落結(jié)構(gòu)和功能如何響應(yīng)環(huán)境變化并沒有得到一致性的答案.

    鄱陽湖地處江西省北部、長江中下游南岸,是我國第一大淡水湖.鄱陽湖湖河形態(tài)多樣,除了主湖區(qū)以及出入湖河流外,還有許多人控或天然湖汊和碟形湖等.作為長江洪水的天然調(diào)蓄場所,鄱陽湖同樣也是生物多樣性的重要區(qū)域,為人類提供了豐富的淡水以及魚類資源,是廣大湖區(qū)人民賴以生存發(fā)展的重要基礎(chǔ)[16].該研究選取鄱陽湖湖區(qū)范圍內(nèi)4個形態(tài)及功能差異明顯的樣地——松門山(主湖區(qū))、南磯山(碟形湖)、饒河(入湖口)、白沙洲(湖汊),利用T-RFLP(末端限制性片段長度多態(tài)性)和Biolog-Eco微平板技術(shù),從群落結(jié)構(gòu)和代謝活性兩方面進行鄱陽湖沉積物微生物群落結(jié)構(gòu)和功能的空間變化研究,探索湖泊底棲微生物與環(huán)境因子的響應(yīng)關(guān)系,以期為進一步探討生態(tài)環(huán)境的地域變化提供理論依據(jù).

    1 材料與方法

    1.1 樣品采集與處理

    于2014年3月在鄱陽湖4個樣地(見圖1)進行沉積物微生物群落結(jié)構(gòu)及功能的調(diào)查研究,4個樣地分別為松門山(SMS,29°12′N、116°11′E)、南磯山(NJS,28°55′N、116°16′E)、饒河(RH,29°01′N、116°29′E)和白沙洲(BSZ,29°10′N、116°37′E).松門山樣地地處鄱陽湖南主湖區(qū),湖面開闊,水流平穩(wěn);南磯山樣地為碟形湖,環(huán)境相對封閉,枯水期與主湖區(qū)分離,有大量的漁業(yè)養(yǎng)殖活動;饒河樣地為入湖口,湖水流速快,并且受到上游來水影響顯著;白沙洲樣地為湖汊,長期受到旅游業(yè)等人為影響.每個樣地設(shè)置3個平行采樣點,分別分布在湖(河)中央以及向兩側(cè)延伸200~500 m位置.沉積物樣品的采集采用抓斗式采泥器,同時進行水深等原位參數(shù)的測定.樣品采集后,立即放入無菌的聚乙烯自封袋中,于冰上放置并迅速運回實驗室.將每份樣品一分為二:一份于4 ℃下保存,用于理化指標以及微生物群落功能多樣性的測定;另一份裝入無菌離心管中保存于-80 ℃冰箱,用于細菌群落結(jié)構(gòu)分析.樣品采集及處理過程盡可能在無菌條件下進行.

    樣地:1—松門山(主湖區(qū));2—南磯山(碟形湖);3—饒河(入湖口);4—白沙洲(湖汊).圖1 鄱陽湖樣地分布Fig.1 Sampling sites in the Poyang Lake

    1.2 理化參數(shù)測定

    樣品pH以及AFDM(去灰分干質(zhì)量,用以表征有機質(zhì)含量)、w(NO3--N)、w(NH4+-N)、w(TP)、w(TN)、w(TOC)測定方法見文獻[17];w(Cu)、w(Zn)、w(Pb)、w(Cd)采用原子吸收分光光度法[18](AA800型圓原子吸收分光光譜儀,美國PE公司)進行測定.

    1.3 沉積物微生物總DNA的提取及T- RFLP分析

    稱取-80 ℃下保存的沉積物樣品0.3 g,采用LU等[19]所運用的方法進行沉積物總DNA的提取.采用通用引物338F(5′-ACTCCTACGGGAGGCWGC- 3′)和1392R(5′-ACGGGCGGTGTGTACA- 3′)[20]進行16S rDNA PCR擴增,上游引物5′端進行6-FAM熒光標記(生工生物工程(上海)股份有限公司).

    將PCR擴增產(chǎn)物用限制性內(nèi)切酶HaeⅢ進行酶切.酶切產(chǎn)物經(jīng)純化(Cycle-Pure Kit,OMEGA)后,進行STR微衛(wèi)星檢測(生工生物工程(上海)股份有限公司).

    T-RFLP分析過程[20]:經(jīng)STR微衛(wèi)星檢測得到TRFs圖譜,舍去低于50 bp和高于600 bp長度的片段.計算每一種片段的相對豐度,并舍去各樣品中相對豐度小于0.5%的片段.

    1.4 Biolog-Eco微平板技術(shù)

    對沉積物細菌功能多樣性的分析方法見文獻[21],用AWCD(average well color development,平均顏色變化率)來評判微生物群落的碳源利用能力,指示微生物群落的代謝活性,其計算公式:

    式中:Ri為第i個非對照孔的吸光值;R0為對照孔的吸光值;31表示有31種碳源.

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

    根據(jù)TRFs在不同樣品圖譜中的分布及豐度,采用PC- ORD 5.0軟件進行Heatmap分析,并利用Venny 2.1.0(http:bioinfogp.cnb.csic.estoolsvennyindex.html)進行韋恩圖的繪制.采用CANOCO 5.0對T- RFLP分析結(jié)果及培養(yǎng)72 h的碳源利用結(jié)果結(jié)合環(huán)境因子進行冗余分析(RDA).采用PRIMER 5.0進行分類組間細菌群落結(jié)構(gòu)及功能差異的ANOSIM分析,并對引起差異的主要類群及碳源進行Simper分析.采用SPSS 20.0軟件進行環(huán)境參數(shù)的差異性檢驗及碳源與環(huán)境因子的相關(guān)性分析,顯著性水平設(shè)置為0.05.

    2 結(jié)果與分析

    2.1 鄱陽湖沉積物理化性質(zhì)及重金屬含量

    鄱陽湖不同樣地沉積物理化因子及重金屬含量見表1.結(jié)果顯示,除w(TOC)和w(Pb)外,沉積物其他理化因子及重金屬質(zhì)量分數(shù)在不同樣地間均差異顯著(P<0.05).鄱陽湖沉積物pH為5.67~6.63,屬弱酸性.各營養(yǎng)參數(shù)均在松門山樣地最低,與碳、氮相關(guān)的參數(shù)均在南磯山樣地最高,w(TP)在饒河樣地最高.4種重金屬含量均在饒河樣地最高,松門山和南磯山樣地最低.利用地積累指數(shù)法并參照鄱陽湖土壤重金屬背景值[22]對重金屬含量進行生態(tài)污染評價,結(jié)果顯示,Pb在4個樣地中均無污染,其他3種重金屬的污染程度表現(xiàn)為Cu>Zn>Cd.

    2.2 細菌群落多樣性及其與環(huán)境因子的相關(guān)性

    對鄱陽湖4個樣地沉積物樣品進行微衛(wèi)星檢測,共獲得的TRFs數(shù)目為52~67,表明鄱陽湖沉積物細菌群落結(jié)構(gòu)均具有較高的多樣性.

    為了分析鄱陽湖沉積物細菌群落結(jié)構(gòu)的空間變化,進行了樣地之間的Heatmap聚類分析(見圖2)及ANOSIM分析.結(jié)果表明,沉積物細菌群落在4個樣地間差異顯著(P<0.05).所有樣品聚為兩支,饒河樣地與其他樣地細菌群落結(jié)構(gòu)相似性最小,單獨為一支;其他樣地聚為一支,其中南磯山樣地與白沙洲樣地細菌群落結(jié)構(gòu)相似性最高,約為37%.由此可見,鄱陽湖沉積物中饒河入湖口樣地細菌群落結(jié)構(gòu)顯著不同于其他樣地,而南磯山碟形湖和白沙洲湖汊樣地中的細菌群落結(jié)構(gòu)更為相似.

    通過對鄱陽湖沉積物中各種TRFs片段的相對豐度進行計算,得到細菌群落組成的韋恩圖(見圖3).由圖3可知,共有24種TRFs廣泛分布于鄱陽湖沉積物樣品中,占所有TRFs片段的87%,這24種TRFs中的8種(TRF77、TRF600、TRF356、TRF572、TRF270、TRF370、TRF537和TRF288)在4個樣地中均處于優(yōu)勢地位(相對豐度均大于2%,下同).此外,TRF251(1.20%)、TRF394(0.70%)、TRF555(0.80%)、TRF589(1.40%)均分布于南磯山、白沙洲和松門山樣地卻不存在于饒河樣地中,同時TRF108(0.07%)、TRF127(0.04%)、TRF444(0.04%)、TRF497(0.04%)4種TRFs片段只存在饒河樣地而不存在其他樣地中.進而進行Simper分析以了解造成饒河樣地與其他樣地細菌群落差異的主要物種組成,結(jié)果顯示,TRF573(1.6%)、TRF588(0.96%)、TRF576(1.20%)、TRF589(1.47%)、TRF571(1.92%)起到了關(guān)鍵作用,關(guān)鍵作用,并且在饒河樣地的相對豐度顯著高于其他樣地(P<0.05).綜上可知,饒河樣地與其他樣地細菌群落結(jié)構(gòu)的差異主要來源于某些處于中等相對豐度(0.20%~2.00%)片段在樣地之間的豐度差異,而不是片段的獨有或缺失.

    注: 不同的數(shù)字代表TRFs片段個數(shù).圖3 基于T- RFLP結(jié)果的韋恩圖Fig.3 Venn Diagram based on the T-RFLP analysis

    為了解鄱陽湖沉積物環(huán)境因子與細菌群落的關(guān)系,進行了冗余分析,結(jié)果見圖4.由Monte Carlo檢驗可知,w(TN)(F=3.2,P=0.002)、AFDM(F=2.1,P=0.006)、w(TOC)(F=16.8,P=0.022)和w(Zn)(F=1.9,P=0.042)是對細菌群落具有顯著影響的環(huán)境因子,并且冗余分析排序圖中前兩個排序軸共解釋了34.38%的變異,表明除了已知環(huán)境因子外,可能還有其他因素對鄱陽湖沉積物細菌群落的分布具有較大影響.在以上四種參數(shù)所確定的二維空間內(nèi),各樣地區(qū)域相互獨立.值得注意的是,饒河樣地細菌群落結(jié)構(gòu)主要受w(Zn)的影響;與此相比,南磯山樣地細菌群落結(jié)構(gòu)主要受營養(yǎng)參數(shù)〔AFDM、w(TOC)和w(TN)〕的影響.

    注: RH-A、 RH-B、 RH-C分別代表饒河樣地的3個平行采樣點, 其他依次類推.圖4 基于T- RFLP的細菌群落與環(huán)境因子的冗余分析Fig.4 RDA analysis of bacterial communities and physiochemical parameters based on T-RFLP

    2.3 細菌功能多樣性及其與環(huán)境因素的相關(guān)性

    為了解沉積物微生物群落功能多樣性,利用Biolog-Eco微平板對微生物利用碳源的能力進行分析,每隔12 h測定其AWCD值,繪制其隨時間變化的動態(tài)曲線,結(jié)果如圖5所示.由圖5可見,鄱陽湖沉積物微生物群落利用碳源的AWCD值隨著培養(yǎng)時間的增加而逐漸升高,0~36 h各樣地的AWCD值較低,表明處于代謝延遲期,>36~96 h為對數(shù)生長期,>96 h后為代謝穩(wěn)定期.結(jié)果顯示,饒河樣地微生物碳源利用能力顯著低于其他樣地(P<0.05),并且南磯山、松門山和白沙洲沉積物微生物群落的代謝能力之間無顯著差異(P>0.05).

    圖5 不同樣地微生物群落的AWCDFig.5 AWCD of the sediment microbial community in different sites

    注: 不同小寫字母代表不同樣地間差異達到P<0.05顯著水平.圖6 培養(yǎng)72 h時不同樣地土壤微生物群落對6種不同類型碳源的利用能力Fig.6 The ability of the sediment microbial community to utilize six types of carbon source after 72 hours′ cultivation

    Biolog-Eco微平板中含有六大類碳源,共計31種.通過分析沉積物微生物對6種不同碳源的利用能力,能夠反映微生物的功能代謝類群,培養(yǎng)72 h時的測定結(jié)果如圖6所示.由圖6可見,總體上鄱陽湖沉積物細菌群落對六大類碳源的利用能力表現(xiàn)為胺類>氨基酸類>酚酸類>羧酸類>多聚物類>碳水化合物類.不同樣地間微生物對六大類碳源的利用能力差異顯著(P<0.05),并表現(xiàn)出不同的變化規(guī)律.微生物群落利用碳水化合物類、羧酸類、氨基酸類、胺類碳源的AWCD值均在松門山樣地最高,而其利用多聚物類和酚酸類碳源的AWCD值在白沙洲樣地最高,說明松門山樣地和白沙洲樣地微生物群落具有最強的碳源利用水平.微生物群落利用多聚物類和碳水化合物類碳源的AWCD值在饒河樣地最低,而其利用羧酸類、酚酸類、氨基酸類和胺類碳源的AWCD值在南磯山樣地最低.整體上,鄱陽湖沉積物微生物群落對胺類、氨基酸類、酚酸類、羧酸類碳源的利用高于多聚物類和碳水化合物類碳源,但是饒河樣地微生物群落對多聚物類和碳水化合物類的利用相比其他樣地顯著降低.

    為了進一步了解鄱陽湖沉積物環(huán)境因子與微生物群落功能多樣性的關(guān)系,利用培養(yǎng)72 h后微生物對各種碳源利用的AWCD值進行冗余分析,結(jié)果如圖7所示.結(jié)果表明,w(NH4+-N)(F=3.5,P=0.016)、pH(F=3.4,P=0.014)、w(NO3--N)(F=3.4,P=0.01)、水深(F=3.2,P=0.018)、w(Zn)(F=2.7,P=0.032)、w(Cu)(F=2.3,P=0.028)、w(Cd)(F=2.3,P=0.038)是影響細菌群落功能多樣性的主要因子,前兩軸共解釋了60.61%的變異.

    × 饒河樣地 ● 南磯山樣地 ▲ 松門山樣地 ■ 白沙洲樣地圖7 基于Biolog- Eco微平板技術(shù)的細菌群落功能與環(huán)境因子的冗余分析Fig.7 Redundancy analysis of bacterial communities function and physiochemical parameters based on Biolog-Eco plate

    饒河樣地細菌群落功能主要受到w(Cu)、w(Zn)、w(Cd)的影響,而南磯山樣地主要受到w(NO3--N)和w(NH4+-N)的影響,松門山樣地主要受到水深和pH的影響,而白沙洲樣地受各種參數(shù)影響均較?。鶕?jù)Pearson相關(guān)性分析(數(shù)據(jù)未列出)可知,微生物對多聚物類和碳水化合物類碳源的利用主要與重金屬含量呈顯著負相關(guān)(P<0.05),故多聚物類和碳水化合物類可能是引起饒河樣地與其他樣地微生物代謝功能差異的主要代謝底物.

    3 討論

    3.1 鄱陽湖沉積物細菌群落結(jié)構(gòu)及微生物功能多樣性的空間變化

    該研究結(jié)果表明,鄱陽湖沉積物細菌群落結(jié)構(gòu)及功能多樣性在各樣地間具有較高的空間異質(zhì)性,并且群落結(jié)構(gòu)和功能呈現(xiàn)協(xié)同變化規(guī)律,均表現(xiàn)為饒河樣地與其他樣地之間差異顯著.Franklin等[23]在2.5~11 m的較小尺度下得到微生物物種相似度隨距離增加而顯著降低的空間格局;而王佳等[24]在20~200 km的相對較大尺度下,對渾河沉積物進行研究時也發(fā)現(xiàn),從上游到下游沉積物細菌群落結(jié)構(gòu)具有明顯的區(qū)域分布特征,呈現(xiàn)由單一向復(fù)雜演替的趨勢;此外,在其他一些研究中,微生物的某些物種或為了適應(yīng)環(huán)境的變化而表現(xiàn)出不同的功能模式,使得結(jié)構(gòu)的變化延遲于功能的變化[14],或因環(huán)境改變而演替出其他具有相同功能的微生物類群,使得微生物功能變化延遲于結(jié)構(gòu)變化[15].因此,淡水湖泊表層沉積物微生物群落結(jié)構(gòu)和功能變化是否協(xié)同發(fā)生、如何協(xié)同發(fā)生及發(fā)生機制還需要更多的研究證明.

    針對基于T-RFLP得到的細菌群落結(jié)構(gòu)結(jié)果進行Simper分析可知,在對饒河樣地與其他樣地的差異中起主導(dǎo)作用的是中等相對豐度(0.20%~2.00%)的TRFs片段,而沉積物中的優(yōu)勢和稀有TRFs片段對細菌群落結(jié)構(gòu)變化的影響較小[9].結(jié)果顯示,鄱陽湖沉積物中等相對豐度的微生物類群可能能夠執(zhí)行一些更加專有的功能,而優(yōu)勢類群可能因含有更多的功能冗余而無法在樣品差異中得以體現(xiàn)[25].秦紅靈等[26]運用T-RFLP技術(shù)對不同土地利用方式下土壤細菌群落結(jié)構(gòu)進行研究時也發(fā)現(xiàn),除了優(yōu)勢片段外,某些相對數(shù)量較少的TRFs片段受土地利用方式的影響也較大.表明在對細菌群落結(jié)構(gòu)進行空間異質(zhì)性的研究過程中,對各種豐度類群的變化關(guān)注不能忽視.

    Biolog-Eco微平板技術(shù)可以有效地評價不同空間水平上同一環(huán)境區(qū)系中菌群功能的變化特征[27].饒河地區(qū)微生物對多聚物類及碳水化合物類碳源的利用能力顯著低于其他樣地,表明這兩類碳源可能是導(dǎo)致鄱陽湖沉積物空間代謝差異的主要碳源.有研究表明[28],在具有不同重金屬污染程度的礦區(qū)土壤中,其微生物群落的碳源代謝差異是由碳水化合物類碳源導(dǎo)致的.也有研究[29]表明,低劑量Cu添加能提高微生物群落對多聚物類碳源的轉(zhuǎn)化與利用的能力.由此可見,微生物對多聚物和碳水化合物類碳源利用的變化可以作為沉積物重金屬含量變化的靈敏、有效的生物學指標.

    3.2 環(huán)境因子對湖泊沉積物細菌群落結(jié)構(gòu)與微生物代謝功能的影響

    營養(yǎng)物質(zhì)〔AFDM、w(TN)、w(TOC)〕及重金屬w(Zn)是細菌群落結(jié)構(gòu)的顯著影響因子,而水深、pH、速效氮〔w(NH4+-N)、w(NO3--N)〕、重金屬〔w(Zn)、w(Cu)、w(Cd)〕是微生物代謝功能多樣性的顯著影響因子.并且w(Zn)是造成饒河樣地與其他樣地微生物群落結(jié)構(gòu)(解釋量16.1%)和功能(解釋量21.4%)差異的共同環(huán)境因子.已有研究[30- 31]證實,較低的w(Zn)能夠促進微生物生長,而較高的w(Zn)則會對微生物生長產(chǎn)生毒害作用.此外,該研究表明,重金屬w(Cu)、w(Cd)對細菌群落結(jié)構(gòu)的影響不顯著,但卻顯著影響微生物代謝功能,這可能是微生物對重金屬長期適應(yīng)的結(jié)果[32].同時,饒河樣地的復(fù)合高重金屬含量對微生物代謝有一定程度的抑制作用,與TENG等[33]對Ag尾礦區(qū)重金屬復(fù)合污染土壤微生物群落研究結(jié)果相符.該研究中,w(TP)也呈現(xiàn)出饒河樣地顯著高于其他樣地的特征,但通過Monte Carlo檢驗,w(TP)并不是影響其細菌群落的顯著影響因子,由此可見,在重金屬和磷同時污染時,淡水湖泊沉積物微生物的群落結(jié)構(gòu)和功能對重金屬變化更為敏感[34].

    4 結(jié)論

    a) 鄱陽湖沉積物中,除w(TOC)和w(Pb)外,所測定的其他理化參數(shù)及重金屬含量在形態(tài)及功能不同的樣地間均具有顯著的空間差異.

    b) 鄱陽湖沉積物細菌群落結(jié)構(gòu)及微生物的碳源利用能力在各樣地間呈現(xiàn)協(xié)同變化特征.其中饒河樣地的細菌群落與其他樣地差異最大,在細菌群落結(jié)構(gòu)上主要體現(xiàn)在中等相對豐度TRFs片段在不同樣地間的豐度差異;而饒河樣地微生物的碳源代謝能力顯著低于其他樣地,主要是由于微生物對多聚物類和碳水化合物類碳源的利用能力顯著降低所致.

    c) 鄱陽湖沉積物細菌群落結(jié)構(gòu)及功能多樣性受環(huán)境因子的影響較大,其中細菌群落結(jié)構(gòu)多樣性受到營養(yǎng)物質(zhì)〔AFDM、w(TN)、w(TOC)〕及重金屬w(Zn)的顯著影響,而微生物的碳源利用功能多樣性受到水深、pH、速效氮〔w(NH4+-N)、w(NO3--N)〕、重金屬〔w(Zn)、w(Cu)、w(Cd)〕的顯著影響.

    [1] 寧淼,葉文虎.我國淡水湖泊的水環(huán)境安全及其保障對策研究[J].北京大學學報:自然科學版,2009(5):848- 854. NING Miao,YE Huwen.Analysis and guaranteeing countermeasures for water environment security of Chinese freshwater lakes[J].Acta Scientiarum Naturalium Universitatis Pekinensis,2009(5):848- 854.

    [2] AZIZULLAH A,KHATTAK M N K,RICHTER P,etal.Water pollution in Pakistan and its impact on public health:a review[J].Environment International,2011,37(2):479- 497.

    [3] SPRING S,SCHULZE R,OVERMANN J O,etal.Identification and characterization of ecologically significant prokaryotes in the sediment of fresh water lakes:molecular and cultivation studies[J].FEMS Microbiology,2000,24:573- 590.

    [4] WANG Guanghua,LIU Junjie,QI Xiaoning,etal.Effects of fertilization on bacterial community structure and function in a black soil of Dehui region estimated by Biolog and PCR-DGGE methods[J].Acta Ecologica Sinica,2008,28(1):220- 226.

    [5] GAD H,WACHENDORF C,JOERGENSEN R G.Response of maize and soil microorganisms to decomposing poplar root residues after shallow or homogenous mixing into soil[J].Journal of Plant Nutrition and Soil Science,2015,178(3):507- 514.

    [6] ZWIRGLMAIER K,KEIZ K,ENGEL M,etal.Seasonal and spatial patterns of microbial diversity along a trophic gradient in the interconnected lakes of the Osterseen Lake District,Bavaria[J].Frontiers in Microbiology,2015,6:1168.

    [7] YANNARELL A C,KENT A D,etal.Temporal patterns in bacterial communities in three temperate lakes of different trophic status[J].Microbial Ecology,2003,46:391- 405.

    [8] SONG H,LI Z,DU B,etal.Bacterial communities in sediments of the shallow Lake Dongping in China[J].Journal of Applied Microbiology,2012,112(1):79- 89.

    [9] 宋洪寧,杜秉海,張明巖,等.環(huán)境因素對東平湖沉積物細菌群落結(jié)構(gòu)的影響[J].微生物學報,2010,50(8):1065- 1071. SONG Hongning,DU Binghai,ZHANG Mingyan,etal.Effect of environmental factors on bacterial community in Lake Dongping sediment[J].Acta Microbiologica Sinica,2012,50(8):1065- 1071.

    [10] MACALADY J L,MACK E E,NELSON D C,etal.Sediment microbial community structure and mercury methylation in mercury-polluted Clear Lake,California[J].Applied and Environmental Microbiology,2000,66(4):1479- 1488.

    [11] 秦旭,高文勝,周文杰,等.控釋肥對蘋果園土壤微生物碳源利用能力的影響[J].果樹學報,2014(5):809- 814. QIN Xu,GAO Wensheng,ZHOU Wenjie,etal.Effects of controlled release fertilizer on microbial community function in apple orchard soil[J].Journal of Fruit Science,2014(5):809- 814.

    [12] 郭星亮,谷潔,陳智學,等.銅川煤礦區(qū)重金屬污染對土壤微生物群落代謝和酶活性的影響[J].應(yīng)用生態(tài)學報,2012,23(3):798- 806. GUO Xingliang,GU Jie,CHEN Zhixue,etal.Effects of heavy metals pollution on soil microbial communities metabolism and soil enzyme activities in coal mining area of Tongchuan,Shanxi Province of northwest China[J].Chinese Journal of Applied Ecology,2012,23(3):798- 806.

    [13] 杜萍,劉晶晶,沈李東,等.Biolog和 PCR-DGGE 技術(shù)解析椒江口沉積物微生物多樣性[J].環(huán)境科學學報,2012,32(6):1436- 1444. DU Ping,LIU Jingjing,SHEN Lidong,etal.Diversity of microorganisms in sediments of the Jiaojiang Estuary as eatimated by Biolog and PCR-DGGE[J].Acta Scientiae Circumstantiae,2012,32(6):1436- 1444.

    [14] YAN Qingyun,BI Yonghong,DENG Ye,etal.Impacts of the Three Gorges Dam on microbial structure and potential function[J].Scientific Reports,2015,5:1- 9.

    [15] LEAR G,BELLAMY J,CASE B S,etal.Fine-scale spatial patterns in bacterial community composition and function within freshwater ponds[J].The ISME Journal,2014,8(8):1715- 1726.

    [16] 蔡其華.健康長江與生態(tài)鄱陽湖:在長江流域湖泊的保護與管理研討會上的主題報告[J].人民長江,2009,40(21):1- 4.

    [17] 張杰,余潮,王自海,等.不同植被群落表層土壤中細菌群落多樣性[J].環(huán)境科學研究,2013,26(8):866- 872. ZHANG Jie,YU Chao,WANG Zihai,etal.Diversity of bacterial communities in surface soils under different vegetation communities[J].Research of Environmental Sciences,2013,26(8):866- 872.

    [18] 田林鋒,胡繼偉,羅桂林,等.貴州百花湖沉積物重金屬穩(wěn)定性及潛在生態(tài)風險性研究[J].環(huán)境科學學報,2012,32(4):885- 894. TIAN Linfeng,HU Jiwei,LUO Guilin,etal.Ecological risk and stability of heavy metals in sediments from Lake Baihua in Guizhou Province[J].Acta Scientiae Circumstantiae,2012,32(4):885- 894.

    [19] LU S M,LIAO M J,ZHANG M,etal.A rapid DNA extraction method for quantitative real-time PCR amplification from fresh water sediment[J].The Journal of Food,Agriculture and Environment,2012,10(34):1252- 1255.

    [20] 劉芳鵬,許光勤,劉倩純,等.南磯山濕地典型植物群落土壤細菌群落結(jié)構(gòu)變化[J].濕地科學,2015,2015(4):444- 450. LIU Fangpeng,XU Guangqin,LIU Qianchun,etal.Variation of structure of bacterial communities in soils with typical plant communities in Nanjishan Wetlands[J].Wetland Science,2015,2015(4):444- 450.

    [21] ZHANG Haihan,HUANG Tinglin,LIU Tingting.Sediment enzyme activities and microbial community diversity in an oligotrophic drinking water reservoir,eastern China[J].Plos One,2013,8(10):e78571.

    [22] 胡春華,李鳴,夏穎.鄱陽湖表層沉積物重金屬污染特征及潛在生態(tài)風險評價[J].江西師范大學學報:自然科學版,2011,35(4):427- 430. HU Chunhua,LI Ming,XIA Ying.The pollution characteristics and potential ecological risk assessment of heavy metals in the surface sediments of Poyang Lake[J].Journal of Jiangxi Normal University(Natural Science),2011,35(4):427- 430.

    [23] FRANKLIN R B,MILLS A L.Multi- scale variation in spatial heterogeneity for microbial community structure in an eastern Virginia agricultural field[J].FEMS Microbiology Ecology,2003,44(3):335- 346.

    [24] 王佳,彭劍峰,宋永會,等.渾河底泥細菌與古菌群落結(jié)構(gòu)空間變化特征研究[J].環(huán)境科學學報,2016,36(1):92- 99. WANG Jia,PENG Jianfeng,SONG Yonghui,etal.Study on the spatial distributions of bacterial and archaeal communities in the sediments of Hun River[J].Acta Scientiae Circumstantiae,2016,36(1):92- 99.

    [25] CHAMBERS L G,GUEVARA R,BOYER J N,etal.Effects of salinity and inundation on microbial community structure and function in a mangrove peat soil[J].Wetlands,2016,36(2):361- 371.

    [26] 秦紅靈,袁紅朝,張慧,等.紅壤坡地利用方式對土壤細菌群落結(jié)構(gòu)的影響[J].土壤學報,2011,48(3):594- 602. QIN Hongling,YUAN Hongchao,ZHANG Hui,etal.Soil bacteria community structure in upland red soil in relation to land use pattern[J].Acta Pedologica Sinica,2011,48(3):594- 602.

    [27] 王強,戴九蘭,吳大千,等.微生物生態(tài)研究中基于BIOLOG方法的數(shù)據(jù)分析[J].生態(tài)學報,2010,30(3):817- 823. WANG Qiang,DAI Jiulan,WU Daqian,etal.Statistical analysis of data from Biolog method in the study of microbial ecology[J].Acta Ecologica Sinica,2010,30(3):817- 823.

    [28] 孫玉青,張莘,吳照祥,等.雄黃礦區(qū)不同砷污染土壤中微生物群落結(jié)構(gòu)及碳源利用特征[J].環(huán)境科學學報,2015,35(11):3669- 3678. SUN Yuqing,ZHANG Zi,WU Zhaoxiang,etal.Soil microbial community structure and carbon source metabolic diversity in the Realgar mining area[J],Acta Scientiae Circumstantiae,2015,35(11):3669- 3678.

    [29] GUO Xingliang,GU Jie,GAO Hua,etal.Effects of Cu on metabolisms and enzyme activities of microbial communities in the process of composting[J].Bioresource Technology,2012,108:140- 148.

    [30] 林海,崔軒,董穎博,等.銅尾礦庫重金屬Cu,Zn對細菌群落結(jié)構(gòu)的影響[J].中國環(huán)境科學,2014(12):3182- 3188. LIN Hai,CUI Xuan,DONG Yingbo,etal.Impact on bacterial community structure of heavy mental of Cu and Zn in copper mine tailings[J].China Environmental Science,2014(12):3182- 3188.

    [31] KELLY J J,TATE R L.Effects of heavy metal contamination and remediation on soil microbial communities in the vicinity of a zinc smelter[J].Journal of Environmental Quality,1998,27(3):609- 617.

    [32] AVIDANO L,GAMALERO E,COSSA G P,etal.Characterization of soil health in an Italian polluted site by using microorganisms as bioindicators[J].Applied Soil Ecology,2005,30(1):21- 33.

    [33] TENG Ying,HUANG Changyong,LUO Yongming,etal.Microbial activities and functional diversity of community in soils polluted with PB- ZN- AG mine tailings[J].Acta Pedologica Sinica,2004,41(1):113- 119

    [34] KELLY J J,TATE R L.Effects of heavy metal contamination and remediation on soil microbial communities in the vicinity of a zinc smelter[J].Journal of Environmental Quality,1998,27(3):609- 617.

    Environmental Factors Influencing the Spatial Distribution of Sediment Bacterial Community Structure and Function in Poyang Lake

    JIN Xiao, KOU Wenbo, YU Haotian, LIU Yajun, MA Yantian, WU Lan*

    School of Life Sciences, Key Laboratory of Poyang Lake Environment and Resource, Ministry of Education,Nanchang University, Nanchang 330031, China

    This study aimed to explore the spatial distribution patterns of bacterial community structure and metabolic function in the surface sediments of Poyang Lake using the T-RFLP and Biolog-Ecoplate.Sediment samples were selected from four sites:Songmenshan(SMS),Nanjishan(NJS),Baishazhou(BSZ)and Raohe(RH).The results demonstrated that both the bacterial community structure and function from RH were significantly different from the other sites(P<0.05),suggesting a similar spatial variation in both bacterial community and microbial metabolic functions.The T-RFLP fragments with a relative abundance between 0.20%-2.00% mainly contributed to the differences of bacterial community composition among the four sites.Meanwhile,the average well color development(AWCD)used to characterize the microbial carbon metabolic activity of the Raohe samples(AWCD=0.65)was remarkably lower than the other sites(P<0.05);further analysis attributed this decrease to the reduced utilization of polymers and carbohydrates.Redundancy analysis(RDA)results showed that nutrient parameters(i.e.,organic matter,TN and TOC)and heavy metal Zn obviously affected the bacterial community structure,while water depth,pH,available nitrogen(i.e.,NH4+-N,NO3--N)and heavy metals(i.e.,Zn,Cu and Cd)had significant effects on the microbial carbon metabolic function.Overall,heavy metal Zn was the most important factor that shaped the bacterial community structure and microbial carbon metabolic function of lake sediments.This study may give insight into the spatial heterogeneity and environmental effects on sediment microbes in large freshwater lakes.

    Poyang Lake; Biolog- Ecoplate; T- RFLP; microbial metabolic function; bacterial community structure

    2016- 10- 08

    2016- 11- 29

    國家自然科學基金項目(31360127,31060082,31660027)

    金笑(1990-),女,山東煙臺人,ncusk724@hotmail.com.

    *責任作者,吳蘭(1969-),女,江西萍鄉(xiāng)人,教授,博士,博導(dǎo),主要從事環(huán)境微生物學研究,wl690902@hotmail.com

    X172

    1001- 6929(2017)04- 0529- 08

    A

    10.13198j.issn.1001- 6929.2017.01.51

    金笑,寇文伯,于昊天,等.鄱陽湖不同區(qū)域沉積物細菌群落結(jié)構(gòu)、功能變化及其與環(huán)境因子的關(guān)系[J].環(huán)境科學研究,2017,30(4):529-536.

    JIN Xiao,KOU Wenbo,YU Haotian,etal.Environmental factors influencing the spatial distribution of sediment bacterial community structure and function in Poyang Lake[J].Research of Environmental Sciences,2017,30(4):529-536.

    猜你喜歡
    結(jié)構(gòu)功能
    也談詩的“功能”
    中華詩詞(2022年6期)2022-12-31 06:41:24
    《形而上學》△卷的結(jié)構(gòu)和位置
    哲學評論(2021年2期)2021-08-22 01:53:34
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    關(guān)于非首都功能疏解的幾點思考
    懷孕了,凝血功能怎么變?
    媽媽寶寶(2017年2期)2017-02-21 01:21:24
    “簡直”和“幾乎”的表達功能
    論《日出》的結(jié)構(gòu)
    創(chuàng)新治理結(jié)構(gòu)促進中小企業(yè)持續(xù)成長
    中西醫(yī)結(jié)合治療甲狀腺功能亢進癥31例
    国产午夜福利久久久久久| av视频在线观看入口| 国产精品一区www在线观看| 亚洲av电影不卡..在线观看| 中国美白少妇内射xxxbb| 1000部很黄的大片| 亚洲18禁久久av| 97热精品久久久久久| 免费看日本二区| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 在线免费十八禁| 国产在线精品亚洲第一网站| 成年女人看的毛片在线观看| 天堂动漫精品| 国产精品嫩草影院av在线观看| 99久久成人亚洲精品观看| 国内精品一区二区在线观看| 村上凉子中文字幕在线| a级毛片免费高清观看在线播放| 日韩高清综合在线| 久久久久久久亚洲中文字幕| 国产熟女欧美一区二区| 国产高清视频在线观看网站| 成人性生交大片免费视频hd| 亚洲中文日韩欧美视频| 老熟妇仑乱视频hdxx| 国产男人的电影天堂91| 性色avwww在线观看| 亚洲美女搞黄在线观看 | 少妇猛男粗大的猛烈进出视频 | 搡老熟女国产l中国老女人| 22中文网久久字幕| 国产精品av视频在线免费观看| 97超视频在线观看视频| 婷婷六月久久综合丁香| av国产免费在线观看| 99热网站在线观看| av免费在线看不卡| 国产av不卡久久| 欧美性猛交╳xxx乱大交人| 亚洲七黄色美女视频| 午夜福利在线观看吧| 国产麻豆成人av免费视频| 国产av麻豆久久久久久久| 免费在线观看影片大全网站| 噜噜噜噜噜久久久久久91| 最近的中文字幕免费完整| 国产av麻豆久久久久久久| ponron亚洲| 国产黄a三级三级三级人| 别揉我奶头~嗯~啊~动态视频| 国产色婷婷99| 欧美又色又爽又黄视频| 一夜夜www| 欧美成人精品欧美一级黄| 日韩国内少妇激情av| 美女内射精品一级片tv| 国产 一区 欧美 日韩| 成人亚洲欧美一区二区av| 九九久久精品国产亚洲av麻豆| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 欧美区成人在线视频| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 国产中年淑女户外野战色| 国产激情偷乱视频一区二区| 午夜福利在线在线| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 国内精品一区二区在线观看| 国产精品国产三级国产av玫瑰| 卡戴珊不雅视频在线播放| 亚洲图色成人| 插逼视频在线观看| 亚洲国产欧美人成| 少妇熟女欧美另类| 国产精品亚洲美女久久久| 蜜臀久久99精品久久宅男| 久久综合国产亚洲精品| 久久精品国产鲁丝片午夜精品| 国产视频一区二区在线看| 美女黄网站色视频| 精品一区二区三区视频在线| 日日撸夜夜添| 嫩草影院入口| 亚洲人成网站高清观看| 不卡一级毛片| 九九热线精品视视频播放| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 精品国产三级普通话版| 春色校园在线视频观看| 99热精品在线国产| 国产午夜精品久久久久久一区二区三区 | 久久久久久大精品| 国产v大片淫在线免费观看| 大香蕉久久网| 亚州av有码| 亚洲欧美日韩高清专用| 熟妇人妻久久中文字幕3abv| 国产精品一及| 国产精品美女特级片免费视频播放器| 中国美白少妇内射xxxbb| 亚洲成人精品中文字幕电影| 欧美成人一区二区免费高清观看| 啦啦啦观看免费观看视频高清| 大型黄色视频在线免费观看| 免费大片18禁| av黄色大香蕉| 不卡一级毛片| 免费观看人在逋| 午夜精品国产一区二区电影 | 日韩人妻高清精品专区| а√天堂www在线а√下载| 午夜精品国产一区二区电影 | 人人妻人人澡欧美一区二区| 中国美女看黄片| 国产成人a区在线观看| 91久久精品电影网| 最近中文字幕高清免费大全6| 中文字幕免费在线视频6| 插逼视频在线观看| 乱人视频在线观看| 亚洲第一电影网av| 亚洲成人久久性| a级毛片a级免费在线| 97超级碰碰碰精品色视频在线观看| 久久鲁丝午夜福利片| 久久久久久久久大av| 午夜精品在线福利| 内射极品少妇av片p| 日本 av在线| 午夜福利视频1000在线观看| 看黄色毛片网站| 深夜a级毛片| 男女视频在线观看网站免费| av专区在线播放| 欧美日韩在线观看h| 卡戴珊不雅视频在线播放| 又黄又爽又刺激的免费视频.| 成人毛片a级毛片在线播放| 中文字幕熟女人妻在线| 中文字幕人妻熟人妻熟丝袜美| 成人av一区二区三区在线看| 18禁在线播放成人免费| 欧美高清性xxxxhd video| av.在线天堂| 成年女人永久免费观看视频| 欧美高清性xxxxhd video| 国产极品精品免费视频能看的| 天堂动漫精品| 欧美另类亚洲清纯唯美| 国产爱豆传媒在线观看| 一个人观看的视频www高清免费观看| 精品久久久久久久久久免费视频| 色在线成人网| 久久久色成人| 成人高潮视频无遮挡免费网站| 一级a爱片免费观看的视频| 国产成人福利小说| 在线观看美女被高潮喷水网站| 国模一区二区三区四区视频| 欧美色欧美亚洲另类二区| 天堂影院成人在线观看| 欧美人与善性xxx| 中文字幕av成人在线电影| 搡女人真爽免费视频火全软件 | 欧美成人精品欧美一级黄| 久久久久国产网址| 男女做爰动态图高潮gif福利片| 丝袜美腿在线中文| 夜夜夜夜夜久久久久| 观看免费一级毛片| 久久久久国产精品人妻aⅴ院| 国产三级中文精品| 免费电影在线观看免费观看| 久久6这里有精品| 在线免费观看不下载黄p国产| 亚洲第一电影网av| 欧美+日韩+精品| 国产三级中文精品| 国产国拍精品亚洲av在线观看| 国内精品宾馆在线| 色5月婷婷丁香| 男女做爰动态图高潮gif福利片| 成人av在线播放网站| 精品少妇黑人巨大在线播放 | 2021天堂中文幕一二区在线观| 欧美zozozo另类| 淫妇啪啪啪对白视频| 国产午夜精品久久久久久一区二区三区 | 久久久久性生活片| 麻豆久久精品国产亚洲av| 免费观看人在逋| 丝袜喷水一区| avwww免费| 麻豆国产av国片精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品一及| 国产精品伦人一区二区| 国产精品综合久久久久久久免费| 成人无遮挡网站| 啦啦啦韩国在线观看视频| 少妇被粗大猛烈的视频| 国产中年淑女户外野战色| 亚洲经典国产精华液单| 一个人免费在线观看电影| 精品久久久久久久久久免费视频| 国产高清激情床上av| 可以在线观看的亚洲视频| 成人高潮视频无遮挡免费网站| 天堂√8在线中文| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 日韩高清综合在线| 黄色配什么色好看| 久久久精品欧美日韩精品| 尤物成人国产欧美一区二区三区| 深夜a级毛片| 欧美色视频一区免费| 看免费成人av毛片| 国产一区二区三区在线臀色熟女| 久久久精品大字幕| 久久精品国产自在天天线| 精品人妻熟女av久视频| 插逼视频在线观看| 夜夜爽天天搞| 精品少妇黑人巨大在线播放 | 香蕉av资源在线| 桃色一区二区三区在线观看| 国产精品久久久久久久电影| 日韩av在线大香蕉| 内射极品少妇av片p| 长腿黑丝高跟| 真实男女啪啪啪动态图| av国产免费在线观看| 人妻少妇偷人精品九色| 在线天堂最新版资源| 日韩强制内射视频| 欧美一级a爱片免费观看看| 日韩中字成人| 国产精品一区www在线观看| av天堂中文字幕网| 69人妻影院| 黄色欧美视频在线观看| 香蕉av资源在线| 色视频www国产| 精品一区二区三区人妻视频| 男人的好看免费观看在线视频| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 成人鲁丝片一二三区免费| 12—13女人毛片做爰片一| 嫩草影院入口| 婷婷亚洲欧美| 一夜夜www| 啦啦啦韩国在线观看视频| 婷婷六月久久综合丁香| 校园春色视频在线观看| 国产一区二区三区av在线 | 男女那种视频在线观看| 亚洲国产精品sss在线观看| 最近在线观看免费完整版| 国语自产精品视频在线第100页| 日韩强制内射视频| 麻豆一二三区av精品| 亚洲综合色惰| 内地一区二区视频在线| 91av网一区二区| 在线免费观看不下载黄p国产| 国产精品久久久久久久电影| 亚洲精品影视一区二区三区av| 日韩在线高清观看一区二区三区| 床上黄色一级片| 亚洲精品乱码久久久v下载方式| 国产三级在线视频| 少妇人妻精品综合一区二区 | 波多野结衣高清无吗| 午夜老司机福利剧场| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 国产aⅴ精品一区二区三区波| 人人妻人人澡欧美一区二区| 一级黄片播放器| 亚洲一区二区三区色噜噜| 毛片一级片免费看久久久久| 久久精品国产亚洲av涩爱 | 国产麻豆成人av免费视频| 国产熟女欧美一区二区| 国产精品久久视频播放| 国语自产精品视频在线第100页| 99九九线精品视频在线观看视频| 在线国产一区二区在线| 亚洲熟妇熟女久久| 日本一本二区三区精品| 国产私拍福利视频在线观看| 直男gayav资源| av国产免费在线观看| 热99re8久久精品国产| 永久网站在线| 国产精华一区二区三区| 亚洲成人av在线免费| 欧美激情在线99| 国产精品不卡视频一区二区| 久久久久久久午夜电影| 波多野结衣巨乳人妻| 日本精品一区二区三区蜜桃| 精品一区二区三区av网在线观看| 国产精品免费一区二区三区在线| 尤物成人国产欧美一区二区三区| 人人妻人人澡欧美一区二区| 综合色丁香网| 久久久久性生活片| 99热精品在线国产| 深夜a级毛片| 国产中年淑女户外野战色| 不卡一级毛片| 国产成人影院久久av| 国产毛片a区久久久久| 嫩草影院精品99| 亚洲丝袜综合中文字幕| 色哟哟·www| 99热这里只有是精品50| 在线免费观看的www视频| 99热这里只有精品一区| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 91午夜精品亚洲一区二区三区| 在线免费十八禁| 在线天堂最新版资源| 国产成人a∨麻豆精品| 久久久久国产网址| 岛国在线免费视频观看| 人妻制服诱惑在线中文字幕| 国产高清有码在线观看视频| 国产精品99久久久久久久久| 欧美日韩乱码在线| 国产探花极品一区二区| 天美传媒精品一区二区| 日韩国内少妇激情av| 超碰av人人做人人爽久久| 看片在线看免费视频| 亚洲国产欧洲综合997久久,| 成人一区二区视频在线观看| 精品久久久久久久久av| 亚洲人成网站在线观看播放| 91久久精品国产一区二区三区| 亚洲在线自拍视频| 国产男人的电影天堂91| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看| 国产色婷婷99| 男女边吃奶边做爰视频| 国产成人a∨麻豆精品| 大香蕉久久网| 久久天躁狠狠躁夜夜2o2o| 国产精品国产高清国产av| 两个人视频免费观看高清| 成人永久免费在线观看视频| 精品免费久久久久久久清纯| av在线亚洲专区| 噜噜噜噜噜久久久久久91| 真人做人爱边吃奶动态| 欧美性猛交╳xxx乱大交人| 午夜日韩欧美国产| 成人漫画全彩无遮挡| 91在线观看av| 99九九线精品视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 91在线精品国自产拍蜜月| 伊人久久精品亚洲午夜| 久久6这里有精品| 婷婷亚洲欧美| 美女内射精品一级片tv| 免费观看在线日韩| 在现免费观看毛片| 国产精品综合久久久久久久免费| 最近2019中文字幕mv第一页| 亚洲精品亚洲一区二区| 中文字幕精品亚洲无线码一区| 欧美+亚洲+日韩+国产| 中文字幕精品亚洲无线码一区| 在线观看av片永久免费下载| 国产日本99.免费观看| av天堂在线播放| 三级毛片av免费| 久久欧美精品欧美久久欧美| 在线看三级毛片| ponron亚洲| 国产真实乱freesex| 在线观看av片永久免费下载| 韩国av在线不卡| 国产精品福利在线免费观看| 极品教师在线视频| 热99在线观看视频| 最近最新中文字幕大全电影3| 99久久中文字幕三级久久日本| 亚洲av成人精品一区久久| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 亚洲国产欧美人成| 国产伦在线观看视频一区| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲av涩爱 | 如何舔出高潮| 草草在线视频免费看| 国产乱人偷精品视频| 国产精品国产三级国产av玫瑰| av在线老鸭窝| 精品久久久久久久久av| 欧美日本视频| 日日干狠狠操夜夜爽| 人人妻人人澡人人爽人人夜夜 | 一个人看视频在线观看www免费| 男女之事视频高清在线观看| 精品福利观看| 美女cb高潮喷水在线观看| 国产极品精品免费视频能看的| 别揉我奶头~嗯~啊~动态视频| 日本免费一区二区三区高清不卡| 美女高潮的动态| 精华霜和精华液先用哪个| 国产久久久一区二区三区| 香蕉av资源在线| 国内精品久久久久精免费| АⅤ资源中文在线天堂| 亚洲无线观看免费| 亚洲国产精品成人综合色| 在线观看66精品国产| 亚洲欧美清纯卡通| 12—13女人毛片做爰片一| 久久人人爽人人片av| 午夜福利在线在线| 精品99又大又爽又粗少妇毛片| 99久久精品热视频| 欧美国产日韩亚洲一区| 国产精品一区二区三区四区久久| 欧美最黄视频在线播放免费| 乱人视频在线观看| 网址你懂的国产日韩在线| 久久韩国三级中文字幕| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 国产一区二区三区av在线 | 国产精品不卡视频一区二区| 国产高清有码在线观看视频| 欧美日韩综合久久久久久| 久久久色成人| 午夜福利在线在线| 九九爱精品视频在线观看| 国产私拍福利视频在线观看| 赤兔流量卡办理| 深爱激情五月婷婷| 成人亚洲欧美一区二区av| 国产亚洲精品久久久久久毛片| 97热精品久久久久久| 亚洲成人久久性| 少妇熟女欧美另类| 成人国产麻豆网| 中文在线观看免费www的网站| АⅤ资源中文在线天堂| av在线亚洲专区| 毛片一级片免费看久久久久| 亚洲欧美精品自产自拍| 免费在线观看成人毛片| av国产免费在线观看| 中文字幕人妻熟人妻熟丝袜美| 免费高清视频大片| 特级一级黄色大片| 久久久久国内视频| 天天躁夜夜躁狠狠久久av| 亚洲一区高清亚洲精品| 日韩欧美一区二区三区在线观看| 99热精品在线国产| 我要搜黄色片| 国产熟女欧美一区二区| 国产麻豆成人av免费视频| 日本色播在线视频| 国产又黄又爽又无遮挡在线| or卡值多少钱| 成年版毛片免费区| 成人av在线播放网站| 亚洲国产精品sss在线观看| 色尼玛亚洲综合影院| 欧美精品国产亚洲| 亚洲性久久影院| 欧美zozozo另类| 亚洲国产高清在线一区二区三| 美女 人体艺术 gogo| 深夜a级毛片| 乱码一卡2卡4卡精品| 男女边吃奶边做爰视频| 亚洲最大成人中文| 国产亚洲精品久久久com| 亚洲真实伦在线观看| 中文字幕免费在线视频6| 国产成年人精品一区二区| 亚洲人成网站高清观看| 久久久久久久久久黄片| 级片在线观看| 久久久成人免费电影| 直男gayav资源| 国内揄拍国产精品人妻在线| 性欧美人与动物交配| 国产男靠女视频免费网站| 亚洲成人久久爱视频| 国产高潮美女av| 成人性生交大片免费视频hd| 亚洲欧美精品自产自拍| 最近手机中文字幕大全| 色哟哟哟哟哟哟| 午夜免费激情av| 精品一区二区三区视频在线观看免费| 18禁在线无遮挡免费观看视频 | av在线播放精品| 看免费成人av毛片| 最近中文字幕高清免费大全6| 人妻制服诱惑在线中文字幕| 校园人妻丝袜中文字幕| av在线亚洲专区| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 在线观看66精品国产| 精品一区二区三区视频在线| 日本色播在线视频| 精品久久久久久久久久免费视频| 欧美一区二区亚洲| 免费观看在线日韩| 在线看三级毛片| av免费在线看不卡| 人人妻人人看人人澡| 卡戴珊不雅视频在线播放| 午夜福利高清视频| 如何舔出高潮| 国产精品三级大全| 看免费成人av毛片| 熟女人妻精品中文字幕| 婷婷色综合大香蕉| 欧美一级a爱片免费观看看| 日本a在线网址| 在线观看免费视频日本深夜| 97超级碰碰碰精品色视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产成人福利小说| 国产探花极品一区二区| 99在线视频只有这里精品首页| 国产精品福利在线免费观看| 亚洲av成人精品一区久久| 我要搜黄色片| 永久网站在线| 精品日产1卡2卡| 精品午夜福利在线看| 成人高潮视频无遮挡免费网站| 毛片女人毛片| 搡老岳熟女国产| 国产色婷婷99| 欧美激情在线99| 欧美+亚洲+日韩+国产| 99热这里只有精品一区| 91久久精品国产一区二区三区| 欧美不卡视频在线免费观看| 热99re8久久精品国产| 男女做爰动态图高潮gif福利片| 亚洲av中文字字幕乱码综合| 一边摸一边抽搐一进一小说| 精品人妻偷拍中文字幕| 尤物成人国产欧美一区二区三区| 舔av片在线| 亚洲第一区二区三区不卡| 亚洲精品456在线播放app| 亚洲内射少妇av| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx性猛交bbbb| 日日撸夜夜添| 亚洲,欧美,日韩| 亚洲欧美日韩高清专用| 久久精品国产自在天天线| 看黄色毛片网站| 美女cb高潮喷水在线观看| 少妇人妻一区二区三区视频| 欧美潮喷喷水| 亚洲婷婷狠狠爱综合网| 亚洲真实伦在线观看| 又黄又爽又刺激的免费视频.| av在线播放精品| 日韩人妻高清精品专区| 乱码一卡2卡4卡精品| 综合色丁香网| 插逼视频在线观看| 欧美成人免费av一区二区三区| 丰满乱子伦码专区| 99视频精品全部免费 在线| 99久国产av精品| 国产一区二区亚洲精品在线观看| 国产精品亚洲一级av第二区| 美女高潮的动态| 色哟哟·www| 亚洲成人av在线免费| 国产精品久久久久久久久免| 丰满乱子伦码专区| 国产亚洲精品久久久久久毛片| 中出人妻视频一区二区| 亚洲国产欧美人成| 国产精品久久久久久久电影| 99久国产av精品| or卡值多少钱| 成年女人看的毛片在线观看| 亚洲美女搞黄在线观看 | 熟女电影av网| 国产精品久久久久久久久免|