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

    福州左海湖細(xì)菌群落的季節(jié)性變化特征及典范對(duì)應(yīng)分析

    2016-10-28 01:28:27劉蘭英陳麗華李玥仁
    關(guān)鍵詞:文庫類群克隆

    劉蘭英,呂 新,陳麗華,李玥仁*

    (1.福建省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所,福建 福州 350003;2.福建省精密儀器農(nóng)業(yè)測(cè)試重點(diǎn)實(shí)驗(yàn)室,福建 福州 350003)

    ?

    福州左海湖細(xì)菌群落的季節(jié)性變化特征及典范對(duì)應(yīng)分析

    劉蘭英1,2,呂新1,2,陳麗華1,2,李玥仁1,2*

    (1.福建省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測(cè)技術(shù)研究所,福建福州350003;2.福建省精密儀器農(nóng)業(yè)測(cè)試重點(diǎn)實(shí)驗(yàn)室,福建福州350003)

    為揭示福州左海湖細(xì)菌群落的季節(jié)變化及其與環(huán)境因子的關(guān)系,通過構(gòu)建16S rRNA 基因克隆文庫對(duì)全湖水體細(xì)菌群落進(jìn)行了調(diào)查,應(yīng)用典范對(duì)應(yīng)分析(Canonical correspondence analysis,CCA)探討影響細(xì)菌群落季節(jié)性變化的環(huán)境因素。結(jié)果表明:左海湖細(xì)菌群落組成中β-變形菌門β-Proteobacteria和藍(lán)細(xì)菌門Cyanobacteria為優(yōu)勢(shì)類群,α-變形菌門α-Proteobacteria、擬桿菌門Bacteroidetes和放線菌門Actinobacteria為亞優(yōu)勢(shì)類群。細(xì)菌群落結(jié)構(gòu)具有明顯的季節(jié)性變化特征,春、夏季Cyanobacteria為優(yōu)勢(shì)類群,秋、冬季β-Proteobacteria為優(yōu)勢(shì)類群。春季水體細(xì)菌群落結(jié)構(gòu)主要受pH、葉綠素a(Chla)、透明度(SD)和總氮(TN)的影響,夏季水體受pH、Chla、SD、總磷(TP)影響較大,秋季水體受TP、Chla、TN影響較大,冬季水體則受水溫(WT)、pH、SD影響較大。本研究發(fā)現(xiàn),左海湖細(xì)菌群落結(jié)構(gòu)的季節(jié)動(dòng)態(tài)與水溫及營養(yǎng)狀況密切相關(guān)。

    左海湖;細(xì)菌群落;16S rRNA 基因克隆文庫;季節(jié)變化;典范對(duì)應(yīng)分析

    湖泊生態(tài)系統(tǒng)中細(xì)菌種類繁多且數(shù)量龐大,約105~108個(gè)·mL-1,占湖泊中生物量的絕大部分,是湖泊環(huán)境中食物網(wǎng)和食物鏈的重要組成部分[1- 2]。細(xì)菌對(duì)湖泊生態(tài)系統(tǒng)的穩(wěn)定性具有調(diào)控作用,細(xì)菌群落結(jié)構(gòu)是其整體功能的表現(xiàn),湖泊細(xì)菌群落的組成與特性在一定程度上反映了湖泊生態(tài)系統(tǒng)的環(huán)境狀況[3]。因此研究湖泊細(xì)菌的群落組成、多樣性的季節(jié)差異,對(duì)于了解細(xì)菌在湖泊生態(tài)系統(tǒng)功能演替中的作用是非常必要的。

    對(duì)水體細(xì)菌多樣性的研究已有很多報(bào)道,主要集中于湖泊、水庫、溪流、海洋等環(huán)境[4],尤其是湖泊細(xì)菌多樣性研究越來越受到研究者的關(guān)注。許多研究表明,湖泊營養(yǎng)狀況不同,湖泊細(xì)菌的群落結(jié)構(gòu)組成也各不相同[5-6]。Zhang等[7]借助克隆文庫和基因芯片相結(jié)合的分子生物學(xué)技術(shù),對(duì)青藏高原高山湖泊細(xì)菌多樣性的研究表明,鹽濃度和海拔與該湖泊中的細(xì)菌群落結(jié)構(gòu)呈顯著性相關(guān)。Adams等[8]通過研究北極24個(gè)湖泊和溪流樣品中的細(xì)菌多樣性,發(fā)現(xiàn)水體的溫度在這些生境的細(xì)菌群落組成中起著關(guān)鍵作用。

    左海湖位于福州市西北部,是典型的半封閉式淺水湖泊,是福州市民主要的休閑和娛樂場(chǎng)所之一。自2008年發(fā)現(xiàn)被稱為“左海水怪”的不名果凍狀生物體以來,左海湖的水體質(zhì)量狀況已引起了廣泛的社會(huì)關(guān)注?!八帧钡某霈F(xiàn)可能與左海湖環(huán)境因子發(fā)生了變化有關(guān),而湖泊中微生物群落的動(dòng)態(tài)變化對(duì)湖泊環(huán)境具有較好的指示作用[9]。本研究在左海湖設(shè)置4個(gè)采樣點(diǎn),基于構(gòu)建16S rRNA基因克隆文庫的方法,對(duì)湖區(qū)細(xì)菌群落組成及季節(jié)變化規(guī)律進(jìn)行了研究,并分析細(xì)菌群落組成與水體環(huán)境因子間的關(guān)系,以期為左海湖水環(huán)境質(zhì)量評(píng)價(jià)及污染防治提供理論依據(jù)。

    1 材料與方法

    1.1樣品采集

    分別于2012年2月(冬季)、4月(春季)、8月(夏季)和10月(秋季)在左海湖進(jìn)行水樣采集。設(shè)置4個(gè)采樣點(diǎn),具體采樣位置的經(jīng)緯度信息見表1,每個(gè)取樣站位設(shè)置3個(gè)重復(fù),所有樣品采集完成后立即放入帶有冰塊的保溫箱中運(yùn)回實(shí)驗(yàn)室,用于后續(xù)試驗(yàn)分析。

    1.2左海湖環(huán)境因子的測(cè)定

    采用透明度盤對(duì)水體的透明度(SD)進(jìn)行現(xiàn)場(chǎng)測(cè)定,其他環(huán)境因子均在實(shí)驗(yàn)室內(nèi)完成測(cè)定。水體的pH值通過酸度計(jì)(PE20,梅特勒)測(cè)定;采用90%丙酮萃取熒光分光光度法對(duì)水體的葉綠素a含量進(jìn)行測(cè)定[10];總氮(TN)、總磷(TP)含量和高錳酸鉀指數(shù)(COD)測(cè)定方法均參照王心芳等[11]的報(bào)道。

    1.3水體細(xì)菌的收集與基因組DNA的提取

    采用8 μm孔徑的定性濾紙,將水樣中的大部分顆粒雜質(zhì)及真核生物過濾除去,再通過聚碳酸酯微孔薄膜(Whatman,0.2 μm,47 mm) 進(jìn)行真空抽濾,然后用無菌水洗滌濾膜2~3次,取下濾膜并將其剪碎后裝入1.5 mL離心管中,于-20℃保存?zhèn)溆?。水體細(xì)菌總基因組DNA的提取參考Wu等[12]的報(bào)道。

    1.4細(xì)菌16S rRNA基因克隆文庫構(gòu)建與篩選

    采用細(xì)菌16S rRNA序列通用型引物27f (5′-AGA GTT TGA TCC TGG CTC AG-3′) 和1492r (5′-TAC GGC TAC CTT GTT ACG ACT T-3′)[13]進(jìn)行PCR擴(kuò)增,設(shè)置3組試驗(yàn)重復(fù),PCR擴(kuò)增方法參考Wu等[4]的報(bào)道。反應(yīng)程序:①94℃,3 min;②94℃,30 s;③60℃,1 min;④72℃,60 s。共30個(gè)循環(huán),最后再72℃延伸10 min?;旌? 組平行的PCR擴(kuò)增產(chǎn)物,經(jīng)純化后的PCR產(chǎn)物連接到pMD18-T載體(TaKaRa公司),轉(zhuǎn)化E.coliDH5α感受態(tài)細(xì)胞。通過藍(lán)白斑法篩選陽性克隆,用測(cè)序引物M13f和M13r進(jìn)行菌落PCR鑒定,PCR擴(kuò)增體系和反應(yīng)條件同上。

    1.5酶切分析與測(cè)序

    將驗(yàn)證為陽性重組子的質(zhì)粒DNA進(jìn)行酶切鑒定,所用的限制性內(nèi)切酶為Csp6 I和Hinf I[14],通過2.5%~3.0% 瓊脂糖凝膠電泳分離驗(yàn)證酶切產(chǎn)物。借助Photoshop軟件分析酶切譜型,統(tǒng)計(jì)各泳道的酶切條帶數(shù)和大小,將完全相同的酶切譜型進(jìn)行歸類,并計(jì)算每個(gè)酶切譜型出現(xiàn)的頻率,挑取具有不同譜型的克隆子到生工生物工程(上海)有限公司進(jìn)行測(cè)序。測(cè)序獲得的16S rRNA基因序列通過DNAMAN 8.0軟件進(jìn)行相互比較,將16S rRNA基因序列相似性大于或者等于97% 的劃為一個(gè)操作分類單元 (Operational Taxonomic Units,OTU)。

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

    細(xì)菌多樣性根據(jù)Shannon-Wiener指數(shù)(H′)和Simpson指數(shù)(1/D)[15]進(jìn)行分析,均勻度指數(shù)(E)可以反映群落中各個(gè)物種多度分布的均勻程度,其計(jì)算參考邱小琮等[16]的報(bào)道,采用文庫覆蓋率(Coverage C)評(píng)估所構(gòu)建的文庫對(duì)目標(biāo)環(huán)境微生物多樣性的體現(xiàn)程度[4]。具體計(jì)算公式如下:

    (1)

    (2)

    式中,Pi=ni/N,ni為第i種OTU的克隆數(shù),N為總克隆數(shù)。

    E=H′/Hmax

    (3)

    式中,Hmax為最大Shannon指數(shù),大小為log2S,S為克隆文庫中總OTU種類數(shù)。

    C=(1-n1/N)×100%

    (4)

    式中,n1代表在16SrRNA基因克隆文庫中僅出現(xiàn)一次的OTU數(shù)量,N為總克隆數(shù)。

    通過CanocoforWindows4.5應(yīng)用軟件對(duì)群落組成數(shù)據(jù)和環(huán)境因子數(shù)據(jù)進(jìn)行典范對(duì)應(yīng)分析,探討細(xì)菌群落與環(huán)境因子之間的關(guān)系,分析時(shí)除了pH值,其余數(shù)據(jù)均進(jìn)行l(wèi)og(x+1)轉(zhuǎn)換[17-18],采用Origin8.0軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行分析并作圖。

    2 結(jié)果與分析

    2.1環(huán)境因子季節(jié)動(dòng)態(tài)

    左海湖各環(huán)境因子指標(biāo)呈現(xiàn)季節(jié)性的變化規(guī)律(表1)。湖水各季節(jié)的pH值介于7.1~7.9,均呈弱堿性;水溫(WT)變化幅度為12.8~32.2℃,四季分明;透明度(SD)的變化范圍介于20~48cm,在秋季呈現(xiàn)最低值;化學(xué)需氧量(COD)的變化幅度較小,波動(dòng)范圍介于7.63~10.51mg·L-1;總氮(TN)含量的變化幅度較大,在0.39~1.28mg·L-1之間波動(dòng),其中夏季各采樣點(diǎn)的總氮含量均較高;總磷(TP)含量的波動(dòng)范圍在0.023~0.052mg·L-1,在夏季呈現(xiàn)最高值,冬季呈現(xiàn)最低值。葉綠素a(Chla)含量的變化范圍介于0.012~0.048mg·L-1,在冬季呈現(xiàn)最低值,其余季節(jié)變化規(guī)律不顯著。

    表1 左海湖水體理化指標(biāo)

    注:①S1~S4代表左海湖4個(gè)不同的采樣站位;②表格中的數(shù)值代表平均值±標(biāo)準(zhǔn)差(3組平行試驗(yàn)的標(biāo)準(zhǔn)差)。

    2.2克隆文庫分析

    每個(gè)季節(jié)樣品分別構(gòu)建1個(gè)16S rRNA基因克隆文庫,并在每個(gè)文庫中隨機(jī)挑取120個(gè)克隆子。通過菌落PCR驗(yàn)證后共得到412個(gè)陽性克隆子,其中春季115個(gè),夏季105個(gè),秋季101個(gè),冬季91個(gè)。通過限制性內(nèi)切酶Csp6Ⅰ和HinfⅠ對(duì)4個(gè)季節(jié)文庫中獲得的陽性克隆子進(jìn)行酶切鑒定,共計(jì)得到298個(gè)酶切結(jié)果。

    分別計(jì)算4個(gè)季節(jié)文庫的覆蓋率Coverage(C)、Shannon-Wiener指數(shù)(H′)、Simpson 指數(shù)(1/D)和均勻度指數(shù)(E)(表2),以比較各季節(jié)樣品的細(xì)菌群落多樣性與豐富度。從表2可以看出,各季節(jié)樣品的文庫覆蓋率(C)為69.5% ~79.1%。張勝男等[19]在梁素海浮游細(xì)菌群落結(jié)構(gòu)研究中指出,當(dāng)湖區(qū)克隆文庫的覆蓋率為65.2%~75.4%時(shí),表明文庫基本可代表細(xì)菌群落結(jié)構(gòu)的組成。因此,由表2中的文庫覆蓋率范圍可知,本研究所構(gòu)建的克隆文庫較為理想,可以很好地反映左海湖中細(xì)菌群落的多樣性。

    多樣性分析表明,春、夏、秋、冬各季節(jié)文庫的Shannon-Wiener指數(shù)分別為2.85、3.30、3.29和3.52,其中冬季樣品的HShannon-Wiener指數(shù)和HSimpson指數(shù)均為最高,說明其細(xì)菌群落多樣性最豐富;相比之下,春季樣品的HShannon-Wiener指數(shù)和HSimpson指數(shù)均較低,說明其細(xì)菌群落多樣性較低。從均勻度指數(shù)(E)可以看出,各克隆文庫中細(xì)菌群落的均勻度較為接近,和多樣性指數(shù)變化趨勢(shì)一致。

    表2 左海湖細(xì)菌16S rRNA基因克隆文庫分析結(jié)果

    2.3細(xì)菌群落組成的季節(jié)動(dòng)態(tài)

    將所構(gòu)建克隆文庫的16S rRNA基因片段送上海生工進(jìn)行雙向測(cè)序,測(cè)序結(jié)果中的上、下游序列通過Invitrogen軟件拼接。根據(jù)97%序列相似性原則,將所有季節(jié)文庫中得到的克隆子劃分為169個(gè)OTUs類型。BLAST比對(duì)結(jié)果顯示,大部分OTUs(90%以上)的16S rRNA基因序列與GenBank 中已知序列具有很高的相似度(97%~100%)。本試驗(yàn)中獲得的412個(gè)陽性克隆子代表的169種基因型歸屬為11門一級(jí)的細(xì)菌分類單元,分別來自變形菌門的α-,β-,γ-和δ-亞族、放線菌門Actinobacteria、擬桿菌門Bacteroidetes、藍(lán)細(xì)菌門Cyanobacteria、綠彎菌門Chloroflexi、酸桿菌門Acidobacteria、疣微菌門Verrucomicrobia和厚壁菌門Firmicutes。

    不同細(xì)菌類群組成的季節(jié)性變化較大(圖1),在各季節(jié)中均大量分布的細(xì)菌類群是β-Proteobacteria和Cyanobacteria。其中,β-Proteobacteria在秋、冬季比例較高,在春季比例相對(duì)較低,為秋、冬季文庫的優(yōu)勢(shì)細(xì)菌類群;Cyanobacteria在春、夏季最豐富,尤其在春季,該細(xì)菌類群占了絕對(duì)的優(yōu)勢(shì)地位。α-Proteobacteria在4個(gè)不同季節(jié)中均有發(fā)現(xiàn),春冬季比例較高,夏秋季比例較低。Actinobacteria和Bacteroidetes在不同季節(jié)里都有分布,冬季數(shù)量較多,其他季節(jié)較少。除了上述幾大主要的細(xì)菌類群,一些非優(yōu)勢(shì)類群,如Chloroflexi、Verrucomicrobia、Firmicutes等并未在所有季節(jié)里均有檢測(cè)到,而且所占的比例也較少,這可能與研究中所采用的16S rRNA基因文庫構(gòu)建方法無法很好地獲取環(huán)境中的非優(yōu)勢(shì)細(xì)菌類群有關(guān)。

    2.4細(xì)菌群落與環(huán)境因子的CCA分析

    不同季節(jié)細(xì)菌群落組成與環(huán)境因子之間的關(guān)系如圖2所示。依據(jù)除趨勢(shì)對(duì)應(yīng)分析(DCA),物種分布的最大梯度長度為2.040(<3.0),具有明顯的線性模型特征[20],因此選擇冗余分析(RDA)探討細(xì)菌群落組成與環(huán)境因子之間的相關(guān)性。本研究中環(huán)境因子第1、2主軸間的相關(guān)系數(shù)為0,說明排序結(jié)果是可信的,圖中箭頭所處的象限表示環(huán)境因子與排序軸間的正負(fù)相關(guān)性[21-22]。從RDA分析圖中可以看出,不同季節(jié)細(xì)菌群落受環(huán)境因子影響不同,本研究選取的7 個(gè)環(huán)境因子對(duì)細(xì)菌群落結(jié)構(gòu)都有一定程度的影響。春季水體細(xì)菌的分布主要受Chla、pH、SD和TN的影響,排序軸1與Chla(r= -0.959 1)、pH(r=-0.873 7)、SD(r=-0.814 4)呈負(fù)相關(guān),軸2與TN(r=-0.767 0)呈負(fù)相關(guān);夏季水體細(xì)菌的分布受pH、SD、TP、Chla的影響較大,軸1與pH(r=-0.960 7)、TP(r=-0.769 3)、Chla(r=-0.742 7)呈負(fù)相關(guān),軸2與SD(r=-0.796 1)負(fù)相關(guān);秋季水體細(xì)菌的分布主要受TP、Chla、TN的影響,環(huán)境因子與主軸1的相關(guān)性較大,其中呈正相關(guān)的為TP(r=0.995 3)、TN(r=0.920 8),與Chla(r=-0.941 1)呈負(fù)相關(guān);冬季水體細(xì)菌的分布受WT、pH、SD影響較大,環(huán)境因子與主軸2的相關(guān)性較大,其中呈正相關(guān)的為pH(r=0.985 3),呈負(fù)相關(guān)的為WT(r=-0.983 3)和SD(r=-0.793 4)。

    3 討論與結(jié)論

    與其他已知淡水生境細(xì)菌類群的組成情況[4-5]相似,本研究結(jié)果表明,左海湖水體中常見的細(xì)菌類群為α-變形菌門α-Proteobacteria、β-變形菌門β-Proteobacteria、放線菌門Actinobacteria、擬桿菌門Bacteroidetes和藍(lán)細(xì)菌門Cyanobacteria,其季節(jié)分布特征表現(xiàn)為β-Proteobacteria和Cyanobacteria在所有季節(jié)里所占的比例均最高。春、夏季Cyanobacteria為優(yōu)勢(shì)類群,秋、冬季β-Proteobacteria為優(yōu)勢(shì)類群。β-Proteobacteria和Cyanobacteria通常被認(rèn)為在富營養(yǎng)化水體中占優(yōu)勢(shì)[4, 6],因此在左海湖檢測(cè)到豐富的β-Proteobacteria和Cyanobacteria細(xì)菌類群在某種程度上說明了左海湖水質(zhì)狀況發(fā)生了變化。而優(yōu)勢(shì)類群顯著的季節(jié)性變化,原因可能為春、夏季氣溫高,有利于Cyanobacteria的生長,因此Cyanobacteria數(shù)量通常有所上升,甚至占了絕對(duì)的優(yōu)勢(shì)地位,由于種間競(jìng)爭(zhēng)的加劇,導(dǎo)致β-Proteobacteria數(shù)量有所下降。另外,β-Proteobacteria和Cyanobacteria在季節(jié)性變化過程中存在明顯的此消彼長的規(guī)律,還可能與水體中營養(yǎng)物質(zhì)發(fā)生變化有關(guān)[23],水體中營養(yǎng)物質(zhì)的高低會(huì)導(dǎo)致一些細(xì)菌的生長受到促進(jìn)或抑制。

    本研究對(duì)福州左海湖水體的細(xì)菌群落組成、環(huán)境因子進(jìn)行典范對(duì)應(yīng)分析,結(jié)果表明環(huán)境因子對(duì)不同季節(jié)里細(xì)菌群落組成的影響程度明顯不同。造成左海湖細(xì)菌群落結(jié)構(gòu)變化的環(huán)境因子包括水體的pH、水溫、總氮、總磷和葉綠素a,表明左海湖細(xì)菌群落結(jié)構(gòu)具有明顯的季節(jié)變化,并且與湖泊中的環(huán)境因子有著不可分割的密切關(guān)系,這與Pollet等[24],Du等[25],Wang等[26]的研究結(jié)果類似。氮、磷是水體主要營養(yǎng)污染物[27],通過環(huán)境因子季節(jié)動(dòng)態(tài)分析發(fā)現(xiàn),春、夏、秋季左海湖水的總氮、總磷含量均高于冬季,表明這期間湖水富營養(yǎng)化程度比冬季湖水嚴(yán)重。已有研究表明,富營養(yǎng)化湖泊中碳酸鹽的化學(xué)平衡決定著水體pH的變化,而碳酸鹽體系的化學(xué)平衡又受到藻類種群生長、轉(zhuǎn)化,光合、呼吸作用等過程的影響,藻類初級(jí)生產(chǎn)力與水體的pH相關(guān)性顯著[28-29]。本研究中細(xì)菌群落組成與水體pH的相關(guān)性可能是由不同水環(huán)境條件下藍(lán)細(xì)菌的生長繁殖變化導(dǎo)致的。藍(lán)細(xì)菌是一類光能營養(yǎng)微生物,其數(shù)量在春、夏季最豐富,這可能與春、夏季水體具有更為豐富的有機(jī)營養(yǎng),較高的水體溫度以及較強(qiáng)的光照強(qiáng)度等有關(guān)。另有研究表明葉綠素a含量的高低體現(xiàn)了藻類生物量的多少[30],本研究中我們發(fā)現(xiàn)藍(lán)細(xì)菌在春、夏、秋季都比較豐富,而葉綠素a含量在這3個(gè)季節(jié)里也都比較高。因此,認(rèn)為左海湖細(xì)菌群落與水體葉綠素a含量的相關(guān)性也可能是由不同季節(jié)時(shí)期藍(lán)細(xì)菌的生長繁殖變化導(dǎo)致的。

    本研究通過構(gòu)建16S rRNA基因克隆文庫的方法,探討了福州左海湖中不同季節(jié)細(xì)菌群落結(jié)構(gòu)的組成,該方法對(duì)于細(xì)菌群落多樣性研究有一定幫助,但也存在一定的局限性。該方法主要采用隨機(jī)挑選克隆子的方法,可能會(huì)高估目標(biāo)環(huán)境中優(yōu)勢(shì)細(xì)菌群落的豐富度,不能了解環(huán)境中的非優(yōu)勢(shì)細(xì)菌群落,而那些數(shù)量上不占優(yōu)勢(shì)的細(xì)菌類群可能對(duì)水體環(huán)境也起著重要的作用。近年來流行的高通量測(cè)序技術(shù)可以更為精確的研究環(huán)境中微生物群落的多樣性,已逐漸得到廣泛應(yīng)用,今后研究中可結(jié)合高通量測(cè)序技術(shù)對(duì)左海湖中的細(xì)菌群落組成進(jìn)行更全面認(rèn)識(shí)。

    [1]SUN J G, JIANG R B, REN T Z, et al. Prospect for farm land and water pollution and microorganism repair in China[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2008, 29(1): 41-47.

    [2]LIU H X, SONG X Y, HUANG L M, et al. Progress in the mechanism of modulation on marine bacterial production[J]. Ecological Science, 2008, 27(1): 61-64.

    [3]YANNARELL A C, KENT A D, LAUSTER G H, et al. Temporal patterns in bacterial communities in three temperate lakes of different trophic status[J]. Microb Ecol, 2003, 46: 391-405.

    [4]WU L, GE G, ZHU G F, et al. Diversity and composition of the bacterial community of Poyang Lake (China) as determined by 16S rRNA gene sequence analysis[J]. World J Microbiol Biotechnol, 2012, 28: 233-244.

    [5]LIU Z, HUANG S, SUN G, et al. Phylogenetic diversity, composition and distribution of bacterioplankton community in the Dongjiang River, China[J]. FEMS Microbiol Ecol, 2012, 80(1): 30-44.

    [6]CHEN Z B, ZHOU Z Y, PENG X, et al. Effects of wet and dry seasons on the aquatic bacterial community structure of the Three Gorges Reservoir[J]. World J Microbiol Biotechnol, 2013, 29(5): 841-853.

    [7]ZHANG R, WU Q, PICENO Y M, et al. Diversity of bacterioplankton in contrasting Tibetan lakes revealed by high-density microarray and clone library analysis[J]. FEMS Microbiol Ecol, 2013, 86(2): 277-287.

    [8]ADAMS H E, CRUMP B C, KLING G W. Temperature controls on aquatic bacterial production and community dynamics in arctic lakes and streams[J]. Environmental Microbiology, 2010, 12(5): 1319-1333.

    [9]TAMAKI H, SEKIGUCHI Y, HANADA S, et al. Comparative analysis of bacteria diversity in freshwater sediment of a shallow eutrophic lake by molecular and improved cultivation-based techniques[J]. Applied and Environmental Microbiology, 2005, 71: 2162-2169.

    [10]楊彩根, 宋學(xué)宏, 孫丙耀. 浮游植物葉綠素a含量簡(jiǎn)易測(cè)定方法的比較[J]. 海洋科學(xué), 2007, 31(1): 6-9.

    [11]王心芳, 魏復(fù)盛, 齊文啟. 水和廢水監(jiān)測(cè)分析方法. [M]. 第4版. 北京: 中國環(huán)境科學(xué)出版社, 2002:216-219.

    [12]WU Q L, ZWART G, SCHAUER M, et al. Bacterioplankton 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.

    [13]MWIRICHIA R, COUSIN S, MUIGAI A W, et al. Bacterial Diversity in the Haloalkaline Lake Elmenteita, Kenya[J]. Curr Microbiol, 2011, 62: 209-221.

    [14]白藍(lán), 趙明文, 賈軍偉, 等. 16S rDNA 克隆文庫法探索轉(zhuǎn)基因香石竹對(duì)土壤細(xì)菌群落的影響[J]. 微生物學(xué)通報(bào), 2012, 39(4): 435-447.

    [15]張繼強(qiáng), 陳文業(yè), 康建軍,等. 敦煌西湖鹽化草甸蘆葦群落特征及多樣性沿水分梯度的分布格局[J]. 水土保持通報(bào), 2013, 33(2): 173-176.

    [16]邱小琮, 趙紅雪. 寧夏沙湖浮游植物群落結(jié)構(gòu)及多樣性研究[J]. 水生態(tài)學(xué)雜志, 2011, 32(1): 20-26.

    [17]LIU X, LU X, CHEN Y. The effects of temperature and nutrient ratios on Microcystis blooms in Lake Taihu, China: An 11-year investigation[J]. Harmful Algae, 2011, 10(3): 337-343.

    [18]MARIA J, MORA R, SCRANTON M I, et al. Bacterial community composition in a large marine anoxic basin: a Cariaco Basin time-series survey[J]. FEMS Microbiol Ecol, 2013, 84: 625-639.

    [19]張勝男, 趙吉睿, 張曉軍, 等. 烏梁素海浮游細(xì)菌群落結(jié)構(gòu)及其對(duì)富營養(yǎng)化因子的響應(yīng)[J]. 微生物學(xué)通報(bào),2014, 41(6): 1082-1093.

    [20]鐘晶晶, 劉茂松, 王玉, 等. 太湖流域河流與湖泊間主要水質(zhì)指標(biāo)的空間關(guān)聯(lián)特征[J]. 生態(tài)學(xué)雜志, 2014, 33(8): 2176-2182.

    [21]TER BRAAK C J F. Canonical correspondence analysis:A new eigenvector technique for multivariate direct gradient analysis[J]. Ecology, 1986, 67(5): 1167-1179.

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

    [23]王曉丹, 翟振華, 趙爽, 等. 北京翠湖表流和潛流濕地對(duì)細(xì)菌多樣性的影響[J]. 環(huán)境科學(xué), 2009, 30(1): 280-288.

    [24]POLLET T, TADONLéKé R D, HUMBERT J F. Spatiotemporal changes in the structure and composition of a less-abundant bacterial phylum (Planctomycetes) in two perialpine lakes[J]. Applied and environmental microbiology, 2011, 77(14): 4811-4821.

    [25]DU J, XIAO K, LI L, et al. Temporal and Spatial Diversity of Bacterial Communities in Coastal Waters of the South China Sea[J]. PLoS One, 2013, 8(6): e66968.

    [26]WANG J J, ZHANG Y, LI Z G, et al. Higher seasonal variation of actinobacterial communities than spatial heterogeneity in the surface sediments of Taihu Lake, China[J]. Can J Microbiol, 2013, 59(5): 353-358.

    [27]孫寓姣, 程陳, 丁愛中, 等. 官廳水庫水質(zhì)特征及水體微生物多樣性的響應(yīng)[J]. 中國環(huán)境科學(xué), 2015, 35(5): 1547-1553.

    [28]朱文昌, 陸敏, 石浚哲. 梅梁湖水體浮游植物與環(huán)境因子的關(guān)系[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù), 2010, 22(3): 27-30.

    [29]楊麗標(biāo), 韓小勇, 孫璞, 等. 巢湖藻類組成與環(huán)境因子典范對(duì)應(yīng)分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2011, 30(5): 952-958.

    [30]馮勝, 李定龍, 秦伯強(qiáng). 太湖水華過程中微生物群落的動(dòng)態(tài)變化[J]. 寧波大學(xué)學(xué)報(bào):理工版, 2010, 23(1): 7-12.

    (責(zé)任編輯:林海清)

    Seasonal Changes on and Environmental Factors Affecting Bacterial Community at Zuohai Lake in Fuzhou

    LIU Lan-ying1, 2, Lü Xin1, 2, CHEN Li-hua1, 2, LI Yue-ren1, 2*

    (1.InstituteforAgriculturalStandardsandTestingTechnology,FujianAcademyofAgriculturalScience,Fuzhou,F(xiàn)ujian350003,China; 2.KeyLaboratoryforPrecisionInstrumentTestsinAgriculturalFields,Fuzhou,F(xiàn)ujian350003,China)

    To determine the seasonal changes on and the environmental factors that affected the bacterial community at Zuohai Lake in Fuzhou, a survey was conducted at 4 sampling sites on the lake from February to August in 2012. A genetic library of their 16S rRNA gene clones was established. Canonical correspondence analysis (CCA) was used to examine the correlation between the bacterial diversity and the environmental factors. The results indicated thatβ-ProteobacteriaandCyanobacteriawere the predominate lineages in the water, followed byα-Proteobacteria,Bacteroidetes, andActinobacteria. A significant seasonal variation in composition of the bacterial community was detected. Among various categories of bacteria,Cyanobacteriaappeared most dominant in the spring and summer; whileβ-Proteobacteria, in the autumn and winter. The results of CCA showed that the bacterial community in the water was mainly affected by pH, chlorophyll a (Chla), transparency (SD), and total nitrogen (TN) in the spring; pH, Chla, SD, and total phosphate(TP), in the summer; TP, Chla, and TN in the fall; and, water temperature, pH, and SD, in the winter. Thus, in general, the bacterial growth in the lake closely related to the temperature as well as the nutrients available in the water.

    Zuohai lake; bacterial community; 16S rRNA gene clone library; seasonal changes; canonical correspondence analysis

    2016-03-22初稿;2016-04-30修改稿

    劉蘭英(1987-),女,研究實(shí)習(xí)員,研究方向:環(huán)境微生物多樣性

    李玥仁(1966-),男,研究員,研究方向:微生物多樣性與生物安全(E-mail: yuerenli@yeah.net)

    福建省自然科學(xué)基金項(xiàng)目(2011J01118);福建省科技計(jì)劃項(xiàng)目——省屬公益類科研院所基本科研專項(xiàng)(2015R1025-9)

    Q 938.8; X 172

    A

    1008-0384(2016)07-758-07

    劉蘭英,呂新,陳麗華,等.福州左海湖細(xì)菌群落的季節(jié)性變化特征及典范對(duì)應(yīng)分析[J].福建農(nóng)業(yè)學(xué)報(bào),2016,31(7):758-764.

    LIU L-Y,Lü X,CHEN L-H,et al.Seasonal Changes on and Environmental Factors Affecting Bacterial Community at Zuohai Lake in Fuzhou[J].FujianJournalofAgriculturalSciences,2016,31(7):758-764.

    猜你喜歡
    文庫類群克隆
    克隆狼
    專家文庫
    浙江:誕生首批體細(xì)胞克隆豬
    優(yōu)秀傳統(tǒng)文化啟蒙文庫
    幽默大師(2020年10期)2020-11-10 09:07:22
    關(guān)于推薦《當(dāng)代詩壇百家文庫》入選詩家的啟事
    中華詩詞(2019年1期)2019-11-14 23:33:56
    薏苡種質(zhì)資源ISSR分子標(biāo)記篩選及親緣關(guān)系分析
    專家文庫
    黑石頂自然保護(hù)區(qū)土壤動(dòng)物功能類群對(duì)季節(jié)動(dòng)態(tài)的響應(yīng)
    抗BP5-KLH多克隆抗體的制備及鑒定
    Galectin-7多克隆抗體的制備與鑒定
    只有这里有精品99| 亚洲精品一区蜜桃| 午夜日本视频在线| 精品少妇内射三级| 成人国产av品久久久| 777米奇影视久久| 黑丝袜美女国产一区| 久久精品国产亚洲网站| 成人亚洲精品一区在线观看| 久久久久国产精品人妻一区二区| 国产片特级美女逼逼视频| 成年美女黄网站色视频大全免费 | videos熟女内射| 欧美日韩av久久| 婷婷色综合www| 精品一区在线观看国产| 久久国产乱子免费精品| 国产深夜福利视频在线观看| 在线天堂最新版资源| 欧美精品高潮呻吟av久久| 中文精品一卡2卡3卡4更新| 久久久a久久爽久久v久久| 国产无遮挡羞羞视频在线观看| 大码成人一级视频| 男女免费视频国产| 国产精品人妻久久久久久| 中文字幕久久专区| 91久久精品国产一区二区成人| 日韩人妻高清精品专区| 色5月婷婷丁香| 久久人人爽av亚洲精品天堂| 色婷婷av一区二区三区视频| 如日韩欧美国产精品一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 国产日韩欧美亚洲二区| 多毛熟女@视频| 午夜日本视频在线| 亚洲国产精品999| 黄色视频在线播放观看不卡| 深夜a级毛片| 成人无遮挡网站| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区三区在线 | 国产无遮挡羞羞视频在线观看| 亚洲欧美中文字幕日韩二区| 久久人人爽人人爽人人片va| av在线观看视频网站免费| 高清黄色对白视频在线免费看 | 在线观看一区二区三区激情| 国产精品欧美亚洲77777| 丝袜喷水一区| av有码第一页| av在线app专区| 香蕉精品网在线| 美女视频免费永久观看网站| 久久久久久人妻| 国产精品久久久久成人av| 国产综合精华液| 中国国产av一级| 如日韩欧美国产精品一区二区三区 | videossex国产| 日韩成人av中文字幕在线观看| 欧美精品国产亚洲| 草草在线视频免费看| 精品国产露脸久久av麻豆| 在线观看www视频免费| 国产美女午夜福利| 免费大片18禁| 日本vs欧美在线观看视频 | 久久狼人影院| 久久这里有精品视频免费| 简卡轻食公司| 日韩免费高清中文字幕av| 久久久久久久久大av| 国产午夜精品久久久久久一区二区三区| 久久久久久久久久成人| 亚洲av在线观看美女高潮| 国产中年淑女户外野战色| 精品久久久久久久久av| 久久久亚洲精品成人影院| videossex国产| 天天操日日干夜夜撸| 伊人久久精品亚洲午夜| 观看免费一级毛片| 久久国内精品自在自线图片| 久久午夜综合久久蜜桃| 成人毛片60女人毛片免费| 国产老妇伦熟女老妇高清| 777米奇影视久久| 免费人妻精品一区二区三区视频| 亚洲在久久综合| 在线观看一区二区三区激情| 校园人妻丝袜中文字幕| 亚洲av福利一区| 看非洲黑人一级黄片| 中文字幕人妻熟人妻熟丝袜美| 亚州av有码| 久久久久久人妻| 国产老妇伦熟女老妇高清| 热re99久久国产66热| 黄色毛片三级朝国网站 | 一本久久精品| 亚洲国产精品一区三区| 国产午夜精品久久久久久一区二区三区| 久久毛片免费看一区二区三区| 国产一区二区三区av在线| 精品亚洲乱码少妇综合久久| 国产亚洲5aaaaa淫片| 女性被躁到高潮视频| 大片电影免费在线观看免费| 久久97久久精品| 少妇被粗大猛烈的视频| 只有这里有精品99| 亚洲熟女精品中文字幕| 免费高清在线观看视频在线观看| 色网站视频免费| 如日韩欧美国产精品一区二区三区 | 国产亚洲精品久久久com| 高清不卡的av网站| 日韩亚洲欧美综合| 纯流量卡能插随身wifi吗| 欧美97在线视频| 80岁老熟妇乱子伦牲交| 亚洲三级黄色毛片| 国产午夜精品一二区理论片| 成人综合一区亚洲| 丁香六月天网| 久久精品国产亚洲av涩爱| 99久久精品一区二区三区| 丝袜喷水一区| 女性被躁到高潮视频| 三级国产精品片| a 毛片基地| 91精品伊人久久大香线蕉| 国产精品国产三级国产av玫瑰| av福利片在线| 精品一品国产午夜福利视频| 色视频在线一区二区三区| a级毛片免费高清观看在线播放| 少妇人妻一区二区三区视频| 久久精品久久久久久噜噜老黄| 日本-黄色视频高清免费观看| 色5月婷婷丁香| 欧美xxxx性猛交bbbb| 色婷婷久久久亚洲欧美| 成年人午夜在线观看视频| 成人无遮挡网站| 99久国产av精品国产电影| 青春草国产在线视频| 国产精品久久久久久久久免| 国产视频内射| 一级av片app| 国产精品成人在线| 老女人水多毛片| tube8黄色片| 纯流量卡能插随身wifi吗| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| av国产久精品久网站免费入址| 春色校园在线视频观看| 亚洲国产最新在线播放| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品古装| 国产成人精品一,二区| 美女视频免费永久观看网站| 中文字幕免费在线视频6| 高清在线视频一区二区三区| 亚洲精品一区蜜桃| 国产色婷婷99| 久久人人爽人人片av| 午夜免费鲁丝| 观看美女的网站| 少妇的逼好多水| 亚洲欧美成人综合另类久久久| 老司机影院毛片| 最新的欧美精品一区二区| 高清黄色对白视频在线免费看 | 国产午夜精品一二区理论片| 丰满迷人的少妇在线观看| 一区在线观看完整版| 99热国产这里只有精品6| 一本色道久久久久久精品综合| 国产熟女午夜一区二区三区 | 久久精品夜色国产| 99热6这里只有精品| 亚洲欧美日韩另类电影网站| 亚洲情色 制服丝袜| √禁漫天堂资源中文www| 成人漫画全彩无遮挡| 亚洲精品久久午夜乱码| 性色av一级| 高清视频免费观看一区二区| 又爽又黄a免费视频| 久久国产精品男人的天堂亚洲 | 中国三级夫妇交换| 欧美区成人在线视频| 高清不卡的av网站| 在线观看三级黄色| av女优亚洲男人天堂| 水蜜桃什么品种好| 久久久久久久久久久丰满| 又粗又硬又长又爽又黄的视频| 成人综合一区亚洲| 亚洲色图综合在线观看| 亚洲国产精品一区二区三区在线| 三级国产精品片| 啦啦啦中文免费视频观看日本| 亚洲色图综合在线观看| 青青草视频在线视频观看| h日本视频在线播放| 国产伦精品一区二区三区四那| 高清在线视频一区二区三区| 国产 精品1| 一本—道久久a久久精品蜜桃钙片| 国产精品一区二区三区四区免费观看| 少妇人妻精品综合一区二区| 日韩av免费高清视频| 人人妻人人澡人人爽人人夜夜| 一级黄片播放器| 搡女人真爽免费视频火全软件| 人妻制服诱惑在线中文字幕| 久久精品久久久久久久性| 日韩强制内射视频| 亚洲精品国产成人久久av| 天堂俺去俺来也www色官网| 精品少妇内射三级| 成人18禁高潮啪啪吃奶动态图 | 国产淫片久久久久久久久| 亚洲精品乱久久久久久| 欧美三级亚洲精品| 久久狼人影院| 一级毛片我不卡| 精品久久久久久电影网| 国产黄片美女视频| 午夜精品国产一区二区电影| 国产综合精华液| 少妇人妻一区二区三区视频| 久久久久久久久久成人| 国产无遮挡羞羞视频在线观看| 最近最新中文字幕免费大全7| 日本黄色日本黄色录像| 国产 精品1| 中国国产av一级| 亚洲精品久久久久久婷婷小说| 亚洲三级黄色毛片| 中文乱码字字幕精品一区二区三区| 精品国产一区二区久久| 久久女婷五月综合色啪小说| 蜜桃久久精品国产亚洲av| av国产精品久久久久影院| 午夜福利在线观看免费完整高清在| 国产乱来视频区| 成人黄色视频免费在线看| 国产成人精品福利久久| 狠狠精品人妻久久久久久综合| 国产黄片视频在线免费观看| 少妇人妻一区二区三区视频| 亚洲国产日韩一区二区| 久久精品国产亚洲av天美| 亚洲成人手机| 亚洲欧美成人综合另类久久久| 赤兔流量卡办理| 少妇熟女欧美另类| 久久久国产一区二区| av网站免费在线观看视频| 国产在线男女| 日韩一区二区三区影片| 成人漫画全彩无遮挡| 亚洲欧美一区二区三区黑人 | 欧美高清成人免费视频www| 狂野欧美白嫩少妇大欣赏| 久久ye,这里只有精品| 欧美精品人与动牲交sv欧美| 欧美高清成人免费视频www| 性色av一级| 中国国产av一级| 国产一区二区在线观看av| 你懂的网址亚洲精品在线观看| 国产日韩欧美在线精品| 最近2019中文字幕mv第一页| 91久久精品国产一区二区三区| av网站免费在线观看视频| 男人添女人高潮全过程视频| 亚洲精品第二区| 99九九线精品视频在线观看视频| 最黄视频免费看| 国产深夜福利视频在线观看| 一区二区三区乱码不卡18| 婷婷色综合大香蕉| 久久久久国产网址| 国产在线一区二区三区精| 婷婷色av中文字幕| 91精品伊人久久大香线蕉| 亚洲国产欧美在线一区| 99精国产麻豆久久婷婷| 男女无遮挡免费网站观看| 蜜桃久久精品国产亚洲av| 久久精品夜色国产| 热re99久久国产66热| 国产成人精品婷婷| 日韩av在线免费看完整版不卡| 久久av网站| 涩涩av久久男人的天堂| 99热6这里只有精品| 国产亚洲欧美精品永久| 国产午夜精品久久久久久一区二区三区| 女人久久www免费人成看片| 精品亚洲成a人片在线观看| 丝袜脚勾引网站| 天堂8中文在线网| √禁漫天堂资源中文www| 欧美日韩视频高清一区二区三区二| 男女边吃奶边做爰视频| 视频区图区小说| 国产午夜精品久久久久久一区二区三区| 少妇丰满av| 亚洲av国产av综合av卡| 成人亚洲精品一区在线观看| 美女中出高潮动态图| 亚洲人成网站在线播| 亚洲欧美一区二区三区国产| 日日摸夜夜添夜夜爱| 国产精品欧美亚洲77777| 永久免费av网站大全| 亚洲性久久影院| 精品一区二区三卡| 国产精品蜜桃在线观看| 2021少妇久久久久久久久久久| 如何舔出高潮| 国产在线男女| freevideosex欧美| 深夜a级毛片| 高清av免费在线| 久久久久人妻精品一区果冻| 视频中文字幕在线观看| 久久韩国三级中文字幕| 老司机亚洲免费影院| 18禁在线无遮挡免费观看视频| 亚洲av国产av综合av卡| 免费大片黄手机在线观看| 国产成人精品一,二区| 高清不卡的av网站| 亚洲精品亚洲一区二区| 大香蕉久久网| 国产成人精品一,二区| 三上悠亚av全集在线观看 | 最近手机中文字幕大全| 18禁在线播放成人免费| 日韩中文字幕视频在线看片| 中文天堂在线官网| 老司机影院毛片| 大又大粗又爽又黄少妇毛片口| 人妻系列 视频| 高清欧美精品videossex| 中文字幕av电影在线播放| 精品久久久久久久久av| 在线观看免费日韩欧美大片 | 国产精品国产av在线观看| av在线老鸭窝| 亚洲人成网站在线观看播放| 国产伦理片在线播放av一区| 五月开心婷婷网| 高清午夜精品一区二区三区| 边亲边吃奶的免费视频| 最新的欧美精品一区二区| av福利片在线| 黑人高潮一二区| 欧美少妇被猛烈插入视频| 中国国产av一级| 亚洲成人手机| 曰老女人黄片| 麻豆精品久久久久久蜜桃| 男人和女人高潮做爰伦理| 久久综合国产亚洲精品| 国产精品人妻久久久影院| 99热国产这里只有精品6| 男的添女的下面高潮视频| 精品国产露脸久久av麻豆| 国产在线免费精品| 麻豆成人av视频| 赤兔流量卡办理| 高清在线视频一区二区三区| 特大巨黑吊av在线直播| 最近2019中文字幕mv第一页| 91久久精品国产一区二区成人| 777米奇影视久久| 亚洲一级一片aⅴ在线观看| 国产色爽女视频免费观看| 啦啦啦啦在线视频资源| 女性被躁到高潮视频| 性色avwww在线观看| 久久精品久久精品一区二区三区| 免费大片18禁| 亚洲精品久久午夜乱码| 亚洲精品久久午夜乱码| 久久99热这里只频精品6学生| 国产精品女同一区二区软件| 国产亚洲最大av| 一级a做视频免费观看| 国产亚洲91精品色在线| 国产亚洲欧美精品永久| 2018国产大陆天天弄谢| 一级毛片 在线播放| 少妇人妻久久综合中文| 精品国产乱码久久久久久小说| 最黄视频免费看| 人妻夜夜爽99麻豆av| 制服丝袜香蕉在线| 三级国产精品欧美在线观看| 一级毛片黄色毛片免费观看视频| 免费观看在线日韩| 波野结衣二区三区在线| 成人亚洲欧美一区二区av| 久久久久久久久久人人人人人人| 人体艺术视频欧美日本| 全区人妻精品视频| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人爽人人夜夜| 亚洲成人av在线免费| 国产av码专区亚洲av| 成人黄色视频免费在线看| 少妇裸体淫交视频免费看高清| 久久久久久久国产电影| 高清黄色对白视频在线免费看 | 亚洲精品中文字幕在线视频 | 99久久综合免费| 大片免费播放器 马上看| 亚洲va在线va天堂va国产| 人妻人人澡人人爽人人| 美女cb高潮喷水在线观看| av一本久久久久| 青春草亚洲视频在线观看| 在现免费观看毛片| 精品99又大又爽又粗少妇毛片| 夜夜看夜夜爽夜夜摸| av免费在线看不卡| 成人漫画全彩无遮挡| 国产黄色视频一区二区在线观看| 欧美成人精品欧美一级黄| 国产伦理片在线播放av一区| 高清黄色对白视频在线免费看 | 高清午夜精品一区二区三区| 免费看光身美女| 国产有黄有色有爽视频| 亚洲欧美成人精品一区二区| 夜夜爽夜夜爽视频| 亚洲国产精品999| 久久精品国产亚洲av涩爱| 亚洲在久久综合| 婷婷色综合大香蕉| 国产亚洲最大av| a级毛片免费高清观看在线播放| 我要看日韩黄色一级片| 国产精品一区www在线观看| 丝袜喷水一区| .国产精品久久| 99热网站在线观看| 国产亚洲精品久久久com| 女性生殖器流出的白浆| 日韩人妻高清精品专区| 成人黄色视频免费在线看| 爱豆传媒免费全集在线观看| 91久久精品国产一区二区三区| 国产av国产精品国产| 免费观看的影片在线观看| 热re99久久精品国产66热6| 精品人妻一区二区三区麻豆| av视频免费观看在线观看| h日本视频在线播放| 国产永久视频网站| 久久久久精品久久久久真实原创| 777米奇影视久久| 亚洲成人手机| 欧美3d第一页| 观看免费一级毛片| 国产精品久久久久久久久免| 91精品一卡2卡3卡4卡| 精品人妻偷拍中文字幕| 国产伦精品一区二区三区视频9| 丰满乱子伦码专区| av国产久精品久网站免费入址| 国产日韩欧美视频二区| 一个人免费看片子| 亚洲无线观看免费| 国产成人精品无人区| 亚洲怡红院男人天堂| 少妇猛男粗大的猛烈进出视频| 国产精品蜜桃在线观看| 免费久久久久久久精品成人欧美视频 | 日韩欧美精品免费久久| 三级国产精品片| 人人妻人人爽人人添夜夜欢视频 | 免费不卡的大黄色大毛片视频在线观看| 久久久久网色| 少妇精品久久久久久久| 日韩在线高清观看一区二区三区| 亚洲av日韩在线播放| 亚洲综合色惰| 人妻人人澡人人爽人人| 日日啪夜夜爽| 欧美日韩综合久久久久久| 国产成人91sexporn| 99久久综合免费| 99久久人妻综合| 一级毛片黄色毛片免费观看视频| 哪个播放器可以免费观看大片| 亚洲真实伦在线观看| 大码成人一级视频| 91精品国产九色| 欧美精品一区二区大全| 国产免费视频播放在线视频| 中文字幕免费在线视频6| 国产一区二区三区综合在线观看 | 欧美成人精品欧美一级黄| 夜夜看夜夜爽夜夜摸| 久久97久久精品| 嘟嘟电影网在线观看| 久久毛片免费看一区二区三区| 国产极品粉嫩免费观看在线 | 亚洲伊人久久精品综合| 黄色毛片三级朝国网站 | 制服丝袜香蕉在线| 久久国产亚洲av麻豆专区| 日韩在线高清观看一区二区三区| 91aial.com中文字幕在线观看| 男男h啪啪无遮挡| 秋霞伦理黄片| 国产精品国产三级国产av玫瑰| 国产乱来视频区| av福利片在线观看| 99久久中文字幕三级久久日本| 欧美日韩av久久| 国产欧美亚洲国产| 精品久久久久久久久av| 日本黄色片子视频| 国产一区二区三区综合在线观看 | 免费观看性生交大片5| 黑人巨大精品欧美一区二区蜜桃 | 久久精品国产亚洲av涩爱| 国产高清不卡午夜福利| 视频中文字幕在线观看| 欧美精品高潮呻吟av久久| 日本av免费视频播放| 少妇丰满av| 99热国产这里只有精品6| 国产成人a∨麻豆精品| 人人妻人人看人人澡| 最黄视频免费看| 最近手机中文字幕大全| 久久亚洲国产成人精品v| av免费观看日本| 最近手机中文字幕大全| 在线精品无人区一区二区三| 一区在线观看完整版| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 中文天堂在线官网| 久久久久久久久大av| 26uuu在线亚洲综合色| 美女cb高潮喷水在线观看| 97超碰精品成人国产| 蜜臀久久99精品久久宅男| 亚洲成人手机| 国产色爽女视频免费观看| 亚洲精品久久久久久婷婷小说| 高清不卡的av网站| 亚洲人成网站在线观看播放| 亚洲精品色激情综合| 在线免费观看不下载黄p国产| 亚洲无线观看免费| 亚洲国产av新网站| 欧美日韩av久久| 国产有黄有色有爽视频| 丰满迷人的少妇在线观看| 精品亚洲乱码少妇综合久久| 人妻系列 视频| 99热这里只有是精品在线观看| 九九在线视频观看精品| 狂野欧美激情性xxxx在线观看| 狠狠精品人妻久久久久久综合| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 欧美一级a爱片免费观看看| 婷婷色综合大香蕉| 狂野欧美白嫩少妇大欣赏| 国产白丝娇喘喷水9色精品| 大码成人一级视频| 美女视频免费永久观看网站| 插逼视频在线观看| 亚洲高清免费不卡视频| 午夜福利,免费看| 夫妻午夜视频| 亚洲第一av免费看| 久久国产乱子免费精品| 黄色视频在线播放观看不卡| 国产在线视频一区二区| 人妻制服诱惑在线中文字幕| 一级,二级,三级黄色视频| 美女内射精品一级片tv| 亚洲丝袜综合中文字幕| 欧美日韩在线观看h| 2021少妇久久久久久久久久久| 久久亚洲国产成人精品v| 我的女老师完整版在线观看| 亚洲欧美日韩另类电影网站| 欧美三级亚洲精品| av国产精品久久久久影院| a级一级毛片免费在线观看| 成人特级av手机在线观看| 免费看不卡的av| 久久综合国产亚洲精品| 亚洲久久久国产精品| 免费观看性生交大片5| 99re6热这里在线精品视频| 99九九线精品视频在线观看视频|