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

    長(zhǎng)江口及鄰近海域表層海水細(xì)菌多樣性及群落結(jié)構(gòu)*

    2015-04-10 05:46:40劉曉輝王健鑫王帥兵樊英萍俞凱成劉明華
    海洋與湖沼 2015年6期
    關(guān)鍵詞:長(zhǎng)江口文庫(kù)類(lèi)群

    劉曉輝 王健鑫 王帥兵 樊英萍 俞凱成 蔣 然 劉明華

    (浙江海洋學(xué)院海洋微生物分子生態(tài)與應(yīng)用實(shí)驗(yàn)室 舟山 316022)

    海洋約占地球表面積的 71%,浩瀚的海洋中有著豐富的生物資源,而海洋微生物以其數(shù)量大、種類(lèi)多、分布廣及適應(yīng)能力強(qiáng)等特點(diǎn)在其中占據(jù)舉足輕重的地位。特定環(huán)境下的微生物群落結(jié)構(gòu)不僅可以反映微生物所處的生境特點(diǎn),同時(shí)也能說(shuō)明微生物在生態(tài)系統(tǒng)中的功能及其在生物地球化學(xué)循環(huán)中所起的作用(Lu et al,2008),因此通過(guò)分析自然環(huán)境中的細(xì)菌群落結(jié)構(gòu)來(lái)研究生態(tài)學(xué)、環(huán)境學(xué)及全球變化,有著重要的現(xiàn)實(shí)意義。

    長(zhǎng)江是我國(guó)第一大河,平均年入海徑流量約9.28×1011m3,占輸入東??倧搅髁康?4.4%,在陸源物質(zhì)向東海的輸送中起著重要作用(張東聲,2011)。大量營(yíng)養(yǎng)鹽、污染物等輸入長(zhǎng)江口及東海,對(duì)河口和沿岸海域的水文、沉積過(guò)程、地貌變化、生物生態(tài)等產(chǎn)生重要影響,另外長(zhǎng)江口及鄰近海域還受到臺(tái)灣暖流、浙江近岸流和黃海沿岸流等水系的共同作用(Ning et al,2011),強(qiáng)烈的季風(fēng)氣候通過(guò)海氣相互作用也影響東海的環(huán)流和溫鹽結(jié)構(gòu),這種多水系的匯合、復(fù)雜的水文狀況和日益頻繁的人類(lèi)活動(dòng)干擾,使該海域不僅是我國(guó)海岸帶陸-海相互作用研究的關(guān)鍵水域(寧修仁等,2004),也是典型近岸生態(tài)系統(tǒng)研究的重要區(qū)域。

    近年來(lái),關(guān)于長(zhǎng)江口及鄰近海域微生物多樣性的研究日益增多,Jiao等(2007)應(yīng)用16S rRNA基因文庫(kù)的方法研究了秋季長(zhǎng)江口和東海陸架區(qū)的微生物多樣性,劉敏等(2008)采用PCR-DGGE技術(shù)對(duì)長(zhǎng)江口外低氧區(qū)的細(xì)菌群落組成進(jìn)行了分析,曾永輝(2008)和潘洛安(2005)也分別研究了長(zhǎng)江口及鄰近海域浮游細(xì)菌多樣性和分布特征,但總體由于采樣站點(diǎn)、季節(jié)、方法的差異,對(duì)該海域細(xì)菌組成及群落結(jié)構(gòu)仍然缺乏比較全面的認(rèn)識(shí)。

    本文采用DAPI染色計(jì)數(shù)法對(duì)長(zhǎng)江口及鄰近海域10個(gè)站點(diǎn)表層海水中的細(xì)菌生物量進(jìn)行統(tǒng)計(jì),并結(jié)合16S rRNA基因克隆文庫(kù)法,構(gòu)建系統(tǒng)發(fā)育樹(shù),研究了長(zhǎng)江口及鄰近海域細(xì)菌群落結(jié)構(gòu)及其多樣性,為了解海洋微生物多樣性提供了分子生物學(xué)方面的基礎(chǔ)數(shù)據(jù),也為研究海洋細(xì)菌與生態(tài)系統(tǒng)之間的關(guān)系提供一定的參考。

    1 材料與方法

    1.1 樣品采集

    2014年7月,搭乘“浙???號(hào)”科學(xué)考察船,利用 SBE 32型采水器(海鳥(niǎo)公司),在長(zhǎng)江口及鄰近海域共采集 10個(gè)站點(diǎn)表層海水樣品,其采樣站點(diǎn)分布如圖1所示。

    圖1 長(zhǎng)江口及鄰近海域采樣站點(diǎn)分布圖Fig.1 The sampling sites of the Changjiang River estuary and adjacent areas

    1.1.1 DAPI樣品采集 采集表層(水面下2m)海水30mL,經(jīng)100μm的篩絹過(guò)濾后,分裝至10mL的無(wú)菌凍存管中,立即加濃度為 1%的戊二醛固定液進(jìn)行固定,黑暗固定 15分鐘后放入低溫冰箱的冷凍格(–20°C)中保存。

    1.1.2 環(huán)境總DNA提取樣品采集 采集表層海水1000mL,裝于滅菌的聚乙烯塑料瓶中。樣品先用100μm的篩絹過(guò)濾,再經(jīng)孔徑為3μm和0.22μm的混合纖維素濾膜(Millipore公司)過(guò)濾,將過(guò)濾后的0.22μm濾膜放入滅菌的 5mL凍存管中,置?20°C冰箱中保存,航次結(jié)束后放置于實(shí)驗(yàn)室?80°C超低溫冰箱中長(zhǎng)期保存。

    1.2 DAPI染色計(jì)數(shù)

    取1—2mL待測(cè)水樣,加入終濃度為10μg/mL的DAPI染色劑,室溫下暗房染色30 min,加入含0.2μm黑色聚碳酸酯膜的過(guò)濾器,濾干后取出濾膜,放在干凈的載玻片上展平,在濾膜中間滴一滴香柏油,置于熒光顯微鏡(100倍油鏡)下觀察并隨機(jī)選擇10個(gè)視野進(jìn)行計(jì)數(shù),只計(jì)算清晰可見(jiàn)且為藍(lán)光的菌體。根據(jù)公式 BN=A×S1/(S2×V)計(jì)算每毫升細(xì)菌數(shù),其中 BN為每毫升細(xì)菌的數(shù)量,A為視野中細(xì)菌的平均數(shù),S1為濾膜的面積,S2為顯微鏡的視野面積,V為過(guò)濾水樣的體積(趙海萍等,2007)。

    1.3 環(huán)境總DNA提取和16S rDNA擴(kuò)增

    利用 FastPrep快速核酸提取儀(MP Biomedicals公司)和FastDNA spin kit for soil試劑盒進(jìn)行海水樣品總DNA提取,核酸蛋白檢測(cè)儀(Bio-Rad公司)測(cè)定DNA濃度和純度,純化后用于PCR擴(kuò)增。擴(kuò)增采用引物 341f (5'-CCTACGGGAGGCAGCAG-3')和 907r(5'-CCGTCAATTCCTTTGAGTTT-3')。50 μL 反應(yīng)體系為: 10×buffer 5 μL,dNTPs 0.5 μL,Primer 341f 1μL,Primer 907r 1 μL,H2O 41 μL,Taq 酶(TaKaRa)0.5 μL,模板DNA 1 μL。PCR擴(kuò)增條件為: 95°C預(yù)變性5 min,94°C 變性 1 min,56°C 退火 30 s,72°C 延伸 2 min,循環(huán)35次,最后72°C延伸10 min。

    1.4 細(xì)菌克隆文庫(kù)構(gòu)建及序列分析

    PCR回收產(chǎn)物與pMD-18T vector(TaKaRa公司)在16°C連接過(guò)夜,將連接產(chǎn)物轉(zhuǎn)化到大腸桿菌DH5α感受態(tài)細(xì)胞(TaKaRa公司),經(jīng)藍(lán)白斑篩選,挑選陽(yáng)性克隆,重新擴(kuò)增插入片斷,將含有合適大小插入片段的克隆送蘇州金唯智生物科技有限公司進(jìn)行測(cè)序。獲得序列后,首先利用Usearch軟件(http: //drive5.com/usearch/)去除嵌合體,再應(yīng)用 BLASTn 程序搜索相似性序列,進(jìn)行系統(tǒng)發(fā)育分析。采用ClustalX (Version 1.8)對(duì)序列進(jìn)行比對(duì)分析,通過(guò)MEGA 5軟件構(gòu)建系統(tǒng)發(fā)育樹(shù)(Kumar et al,2004),采用Neighbor-Joining建樹(shù)方法,建樹(shù)結(jié)果進(jìn)行 Bootstrap1000次系統(tǒng)檢驗(yàn),利用PHYLIP軟件包中DNASIS程序計(jì)算距離矩陣,利用 DOTUR軟件將相似性>97%的序列歸為一個(gè)OTU,計(jì)算多樣性指數(shù)并繪制稀釋曲線(Schloss et al,2005)。

    2 結(jié)果與分析

    2.1 DAPI染色計(jì)數(shù)

    DAPI染色計(jì)數(shù)結(jié)果如表1所示,其中B03站點(diǎn)總菌數(shù)(1.48×106cells/mL)最高,其次為 A05站點(diǎn)(1.23×106cells/mL),而 A04 站點(diǎn)的總菌數(shù)(1.16×105cells/mL)最低,10個(gè)站點(diǎn)的平均值為 4.92×105cells/mL,細(xì)菌豐度總體較高,但各站點(diǎn)之間細(xì)菌豐度也存在明顯差異,細(xì)菌豐度的高值區(qū)主要集中于舟山外海一帶,而低值區(qū)主要集中于長(zhǎng)江口區(qū)域,整體變化趨勢(shì)為從長(zhǎng)江口至外海區(qū)域細(xì)菌豐度逐漸增加。

    2.2 細(xì)菌多樣性分析

    2.2.1 細(xì)菌 16S rDNA 文庫(kù)分析 每個(gè)站點(diǎn)選取110個(gè)左右陽(yáng)性克隆進(jìn)行測(cè)序,片段長(zhǎng)度為 550—700bp;利用 DOTUR軟件,將相似性>97%的克隆子歸為一個(gè)OTU,各站點(diǎn)的克隆文庫(kù)基本情況及多樣性指數(shù)如表2所示,細(xì)菌群落的稀釋曲線如圖2所示。

    表2 長(zhǎng)江口及鄰近海域表層海水細(xì)菌多樣性指數(shù)一覽表Tab.2 Bacterial diversity indices of the surface seawater in Changjiang River estuary and adjacent areas

    Chao1和ACE指數(shù)用于計(jì)算群落的物種豐富度,在不同算法的計(jì)算下它們有著相同的變化趨勢(shì)。由表2可知,豐富度指數(shù)Chao 1與ACE指數(shù)同時(shí)顯示A02站點(diǎn)擁有最高豐富度,C03站點(diǎn)次之,豐富度最低的為C04站點(diǎn)。Shannon和Simposon指數(shù)用于計(jì)算群落的多樣性,由于計(jì)算方法的差異它們?cè)跀?shù)值上的變化趨勢(shì)是相反的,但兩者均顯示 A02站點(diǎn)細(xì)菌多樣性最高,C03站點(diǎn)次之,C04站點(diǎn)最低。結(jié)合各站點(diǎn)的群落豐富度和多樣性指數(shù),可以看出 A02站點(diǎn)的細(xì)菌群落多樣性最高。由稀釋性曲線可知,A02和C03站點(diǎn)的稀釋曲線仍然有上升的趨勢(shì),表明該站點(diǎn)中細(xì)菌群落多樣性很高,克隆子數(shù)量還未完全覆蓋其多樣性,這與其豐富度和多樣性指數(shù)結(jié)果是一致的,其余站點(diǎn)的稀釋性曲線均趨于平穩(wěn),基本上能反映真實(shí)環(huán)境中細(xì)菌的群落結(jié)構(gòu)情況。

    2.2.2 細(xì)菌克隆文庫(kù)物種類(lèi)群分析 通過(guò)BLASTn搜索相似序列,可將10個(gè)站點(diǎn)共1041個(gè)克隆子分為 11個(gè)細(xì)菌門(mén)類(lèi),分別為變形菌門(mén)(Proteobacteria)、擬桿菌門(mén)(Bacteroidetes)、放線菌門(mén)(Actinobacteria)、藍(lán)細(xì)菌門(mén)(Cyanobacteria)、厚壁菌門(mén)(Firmicutes)、綠彎菌門(mén)(Chloroflexi)、酸桿菌門(mén)(Acidobacteria)、浮霉菌門(mén)(Planctomycetes)、纖維桿菌門(mén)(Fibrobacteres)、疣微菌門(mén)(Verrucomicrobia)和柔膜菌門(mén)(Tenericutes),另有 2.02%的細(xì)菌分類(lèi)地位不確定。變形菌門(mén)又分為 α-變形菌綱(Alphaproteobacteria)、β-變形菌綱(Betaproteobacteria)、γ-變形菌綱(Gammaproteobacteria)、δ-變形菌綱(Deltaproteobacteria)和ε-變形菌綱(Epsilonproteobacteria)五個(gè)綱。各細(xì)菌類(lèi)群在每個(gè)站點(diǎn)克隆文庫(kù)中的比例如圖3所示。

    圖2 十個(gè)文庫(kù)的稀釋性曲線Fig.2 The rarefaction curve of the 10 clone libraries

    圖3 細(xì)菌類(lèi)群在不同采樣點(diǎn)克隆文庫(kù)中的比例Fig.3 The ratios of bacterial groups in clone library of different sites

    綜合 10個(gè)站點(diǎn)文庫(kù)的結(jié)果可以發(fā)現(xiàn),變形菌門(mén)為研究區(qū)域中最大的優(yōu)勢(shì)菌群(占 60.61%),其中 α-變形菌占 50.62%,β-變形菌占 4.32%,γ-變形菌占4.13%,δ-變形菌占 1.44%,ε-變形菌占 0.10%;僅次于變形菌門(mén)的是擬桿菌門(mén)(占15.18%);第三大優(yōu)勢(shì)菌群為放線菌門(mén)(占14.79%)。

    對(duì)各站點(diǎn)文庫(kù)進(jìn)行比較分析,可見(jiàn)除 B03站點(diǎn)外,變形菌門(mén)是其余 9個(gè)站點(diǎn)文庫(kù)的主要優(yōu)勢(shì)類(lèi)群,其中C04站點(diǎn)文庫(kù)比例最高(占90.66%),變形菌門(mén)中又以 α-變形菌為主要類(lèi)群,其余類(lèi)群(特別是 δ-變形菌和ε-變形菌)所占比例相對(duì)較少,ε-變形菌僅在A02站點(diǎn)有發(fā)現(xiàn)且含量很低(僅占 0.93%)。作為第二大優(yōu)勢(shì)類(lèi)群的擬桿菌,僅在B03站點(diǎn)克隆文庫(kù)中的含量超過(guò)變形菌門(mén)(占43.16%),成為該站點(diǎn)文庫(kù)的主要類(lèi)群,其它站點(diǎn)均低于變形菌門(mén)。放線菌門(mén)在10個(gè)站點(diǎn)克隆文庫(kù)中均有分布,其中 A02站點(diǎn)放線菌含量最高(占35.51%),與該站點(diǎn)的變形菌門(mén)含量相當(dāng)。

    除了A01站點(diǎn)未見(jiàn)相應(yīng)克隆,藍(lán)細(xì)菌門(mén)在各文庫(kù)的分布比較均勻,但含量不高,最高值出現(xiàn)在 C05站點(diǎn)(占12.26%)。浮霉菌門(mén)只在A02和C03站點(diǎn)文庫(kù)中有分布,并且含量極低(分別為0.93%和0.92%)??拷L(zhǎng)江口的A01和A02站點(diǎn)文庫(kù)含有其它站點(diǎn)沒(méi)有的厚壁菌門(mén),其含量分別占各自文庫(kù)的8.16%和3.74%。綠彎菌門(mén)主要分布于A斷面(A01、A02、A03和A04站點(diǎn)),其中 A02文庫(kù)中綠彎菌門(mén)含量最高(占 3.74%)。柔膜菌門(mén)僅出現(xiàn)在 C05站點(diǎn)文庫(kù)中,且僅占 0.94%。A04站點(diǎn)文庫(kù)含有其它站點(diǎn)文庫(kù)沒(méi)有的纖維桿菌門(mén)(占1.98%)。酸桿菌門(mén)和疣微菌門(mén)只在 A02站點(diǎn)的克隆文庫(kù)中出現(xiàn),含量分別為2.80%和1.87%。同時(shí)分類(lèi)地位不確定的類(lèi)群在各站點(diǎn)文庫(kù)中都有極少量的分布。

    2.3 細(xì)菌系統(tǒng)發(fā)育分析

    2.3.1 變形菌門(mén)系統(tǒng)發(fā)育分析 使用MEGA5.0軟件,對(duì)各站點(diǎn)典型變形菌序列與BLASTn比對(duì)相似序列構(gòu)建系統(tǒng)發(fā)育樹(shù)(圖4),克隆子相似菌株分離源的分析結(jié)果如表3所示。

    由系統(tǒng)發(fā)育樹(shù)和相似序列生境分析結(jié)果可知,α-變形菌綱克隆子的相似序列主要來(lái)自表層海水、富營(yíng)養(yǎng)化淡水水域,個(gè)別源自深海熱液口;β-變形菌綱克隆子的相似序列主要來(lái)自淡水湖泊和海水環(huán)境;γ-變形菌綱克隆子的相似序列主要來(lái)自表層水體、海洋沉積物和土壤環(huán)境;δ-變形菌綱克隆子相似序列的分離源主要是墨西哥灣石油泄漏海區(qū)和西班牙海域碳?xì)浠衔镂廴镜某练e物,推測(cè)該菌可能與石油烴污染環(huán)境相關(guān);ε-變形菌則在水體環(huán)境中分布較少,唯一的基因型A02B21與湖泊沉積物中的ε-變形菌同源性達(dá)99%。SAR類(lèi)群是海洋中比較典型的細(xì)菌類(lèi)群,也是影響海洋碳循環(huán)的重要因素,由圖表可知,SAR11的相似菌株來(lái)自美國(guó)普吉特海灣表層海水和墨西哥灣石油污染海域;SAR406的分離源則分別來(lái)自沿海表層和深海的海水環(huán)境,說(shuō)明該類(lèi)群細(xì)菌對(duì)于壓力、光照等環(huán)境具有很強(qiáng)的適應(yīng)性。

    圖4 根據(jù)細(xì)菌16S rDNA基因序列構(gòu)建的變形菌系統(tǒng)發(fā)育樹(shù)Fig.4 The phylogenetic tree of Proteobacteria according to 16S rDNA gene sequences

    2.3.2 其它細(xì)菌類(lèi)群系統(tǒng)發(fā)育分析 除變形菌門(mén)外,各站點(diǎn)其它典型細(xì)菌序列與BLASTn比對(duì)相似序列所構(gòu)建的系統(tǒng)發(fā)育樹(shù)如圖5所示,克隆子相似菌株分離源的分析結(jié)果見(jiàn)表4。

    由系統(tǒng)發(fā)育樹(shù)的結(jié)果可知,變形菌門(mén)以外的其它代表性細(xì)菌克隆,其相似序列大多來(lái)自海水環(huán)境,少數(shù)來(lái)源于沉積物和淡水環(huán)境。其中擬桿菌門(mén)(黃桿菌綱)克隆子的相似序列主要來(lái)自于河口和海灣等生境;放線菌門(mén)、綠彎菌門(mén)、厚壁菌門(mén)等細(xì)菌克隆子的相似菌株則大多來(lái)源于受芳香烴、石油、多環(huán)芳烴等污染的海域、土壤或沉積物;酸桿菌門(mén)克隆子的相似序列來(lái)自太湖東北部爆發(fā)微囊藻水華的水樣;疣微菌門(mén)的相似序列則來(lái)自夏威夷群島沿岸海水;藍(lán)細(xì)菌門(mén)克隆子相似序列的分布比較廣泛,主要來(lái)自各海域的表層海水;浮霉菌門(mén)克隆子相似序列的生境來(lái)自西南太平洋勞盆地深海熱液區(qū)沉積物;另外分布較少的纖維桿菌門(mén)和柔膜菌門(mén),其克隆子相似序列分別來(lái)自瑞典爆發(fā)水華的淡水湖泊和美國(guó)佐治亞州Sapelo島沿岸表層海水。在分類(lèi)地位不明確的細(xì)菌中,基因型B02B83(1個(gè)克隆子)、C03B1(4個(gè)克隆子)、和C04B90(1個(gè)克隆子)在發(fā)育樹(shù)中處于同一個(gè)單獨(dú)的分支,其相似序列主要來(lái)自我國(guó)南海近岸表層海水、阿根廷Espejo湖水以及廣東湛江灣海水。

    表3 變形菌克隆子序列比對(duì)及相似菌株分離源分析Tab.3 The taxonomic affiliation of Proteobacteria and the closest species origins

    3 討論

    3.1 長(zhǎng)江口及鄰近海域表層海水細(xì)菌豐度

    從長(zhǎng)江口到舟山外海的區(qū)域是受人類(lèi)活動(dòng)影響最大的區(qū)域,每年有著大量的陸源有機(jī)物質(zhì)輸入,而該區(qū)域又是臺(tái)灣暖流和長(zhǎng)江沖淡水交匯的區(qū)域(唐曉暉等,2004),復(fù)雜的水文地理?xiàng)l件使得該區(qū)域的細(xì)菌豐度同時(shí)受眾多環(huán)境因素的影響。趙三軍等(2003)研究發(fā)現(xiàn)秋季長(zhǎng)江口外海域表層海水的細(xì)菌豐度為3.97×105—1.162×106cells/mL 之間,平均值為 7.19×105cells/mL;劉晶晶等(2011)調(diào)查了夏季長(zhǎng)江口及鄰近海域浮游細(xì)菌的數(shù)量,結(jié)果為 6.92×105—5.54×106cells/mL之間;王新等(2010)對(duì)春季東海赤潮高發(fā)區(qū)的浮游細(xì)菌豐度調(diào)查結(jié)果為 5.85×104—9.26×105cells/mL之間;上述研究結(jié)果與本次調(diào)查結(jié)果(細(xì)菌豐度在 1.16×105—1.48×106cells/mL 之間,平均值為 4.92×105cells/mL)都有一定差異,原因主要可能是采樣季節(jié)不同以及環(huán)境特征的改變。

    圖5 根據(jù)細(xì)菌16S rDNA基因序列構(gòu)建的其它分支細(xì)菌系統(tǒng)發(fā)育樹(shù)Fig.5 The phylogenetic tree of other bacteria according to 16S rDNA gene sequences

    就細(xì)菌豐度的空間變化而言,本研究發(fā)現(xiàn)細(xì)菌豐度分別在A05、B03、C05等靠近外海的站點(diǎn)形成高值區(qū),在相對(duì)離岸較近的 A01—A04和 C03—C04站點(diǎn)海域形成低值區(qū),細(xì)菌豐度的總體空間變化趨勢(shì)是外海較近岸高,其主要原因可能是外海區(qū)域的幾個(gè)站點(diǎn)位于營(yíng)養(yǎng)鹽豐富、光照充足的高生產(chǎn)力海域,浮游植物的繁盛促進(jìn)了細(xì)菌的大量生長(zhǎng);相反,在長(zhǎng)江口近岸海域,由于水團(tuán)之間混合所產(chǎn)生的高濁度和低透光率影響了浮游植物生長(zhǎng),從而形成海洋細(xì)菌低密度區(qū)。季倩(2008)研究認(rèn)為長(zhǎng)江口靠近最大混濁帶異養(yǎng)細(xì)菌豐度較低,由長(zhǎng)江口向外海其豐度顯著增加;寧修仁等(2004)也認(rèn)為在長(zhǎng)江口及鄰近區(qū)域的沖淡水區(qū)中部以及沖淡水和海水交匯處是低細(xì)菌豐度和活性區(qū)域,上述報(bào)道與本次的研究結(jié)果基本一致。本次研究中B02站點(diǎn)雖然處在A03和C03站點(diǎn)相近經(jīng)度區(qū)域,但細(xì)菌生物量處于較高水平,原因可能是 B02站點(diǎn)位于擁有大量天然海藻場(chǎng)和人工貽貝養(yǎng)殖場(chǎng)的枸杞島海域,養(yǎng)殖區(qū)域的存在對(duì)細(xì)菌豐度的影響較大,有研究表明浮游植物和海藻光合作用產(chǎn)物以及部分海藻降解產(chǎn)生的腐殖質(zhì)促進(jìn)了浮游細(xì)菌的新陳代謝(V?r?s et al,2003),同時(shí)貝類(lèi)養(yǎng)殖的殘餌、糞便所分解的有機(jī)質(zhì)含量較高,在該海域獨(dú)特的水文地理?xiàng)l件下可能導(dǎo)致水體中細(xì)菌豐度偏高。

    表4 其它細(xì)菌克隆子序列比對(duì)及相似菌株分離源分析Tab.4 The taxonomic affiliation of other bacteria and the closest species origins

    3.2 長(zhǎng)江口及鄰近海域細(xì)菌多樣性及群落結(jié)構(gòu)

    近年來(lái)國(guó)內(nèi)外有許多針對(duì)長(zhǎng)江口海域表層海水細(xì)菌多樣性及群落結(jié)構(gòu)的相關(guān)研究,如李和陽(yáng)等(2012)研究了長(zhǎng)江口典型站位浮游細(xì)菌類(lèi)群的組成特征,結(jié)果表明變形菌門(mén)(α-,β-和 γ-變形菌綱)為該海域優(yōu)勢(shì)菌群;Feng等(2009)對(duì)長(zhǎng)江口及鄰近海域海水和沉積物樣品進(jìn)行較全面的研究,結(jié)果顯示變形菌門(mén)是該海洋環(huán)境中含量最豐富的種類(lèi)(占克隆數(shù)的72.9%),其中α-變形菌綱和γ-變形菌綱為該海域海水中的最大優(yōu)勢(shì)菌群;Sekiguchi等(2002)分析了長(zhǎng)江口特定站位細(xì)菌的遺傳多樣性,發(fā)現(xiàn)優(yōu)勢(shì)菌為 α-變形菌和γ-變形菌。從以上研究結(jié)果可以發(fā)現(xiàn)變形菌門(mén)是長(zhǎng)江口及鄰近海域表層海水中的優(yōu)勢(shì)菌群,作者在本實(shí)驗(yàn)中檢測(cè)到變形菌門(mén)含量達(dá)到總克隆數(shù)的60.61%,是本次研究海域表層海水中的最大優(yōu)勢(shì)類(lèi)群,這與以往研究結(jié)果基本一致。變形菌門(mén)作為本實(shí)驗(yàn)的優(yōu)勢(shì)類(lèi)群,在長(zhǎng)江口及鄰近海域生態(tài)系統(tǒng)中發(fā)揮著重要作用,尤其是含量豐富、分布廣泛的α-變形菌綱(占總克隆數(shù)的50.62%)。目前對(duì)α-變形菌的研究發(fā)現(xiàn),其在低價(jià)硫化物的轉(zhuǎn)化過(guò)程中起間接作用,常以H2S作為代謝中的電子供體(李濤等,2008);多數(shù)α-變形菌還具有廣泛的耐鹽性,并且利用 DOC等碳源的能力很強(qiáng)(Sj?stedt et al,2012)。本實(shí)驗(yàn)中有5個(gè)站點(diǎn)的 α-變形菌含量高達(dá) 50%以上,說(shuō)明相應(yīng)環(huán)境特征非常適合該細(xì)菌類(lèi)群的生長(zhǎng),同時(shí)也驗(yàn)證了從長(zhǎng)江口到外海區(qū)域鹽度和營(yíng)養(yǎng)鹽急劇變化條件下 α-變形菌仍為本次研究海域的優(yōu)勢(shì)菌群。β-變形菌是世界范圍內(nèi)淡水湖泊和河口區(qū)域的優(yōu)勢(shì)菌群(Mueller-Spitz et al,2009;Cottrell et al,2005),但本次研究中β-變形菌含量?jī)H占總文庫(kù)的4.32%,這可能與長(zhǎng)江口沖淡水的影響有關(guān)。γ-變形菌是海水和沉積物環(huán)境中的常見(jiàn)細(xì)菌類(lèi)群,與光能營(yíng)養(yǎng)、趨堿性、氨氧化和硫降解等過(guò)程相關(guān)(Bowman et al,2005;Ravenschlag et al,1999),目前較多報(bào)道顯示γ-變形菌在沉積物環(huán)境中屬于優(yōu)勢(shì)類(lèi)群(Bowman et al,2003;Inagaki et al,2003),但在海水環(huán)境中不是優(yōu)勢(shì)菌群,我們的研究結(jié)果也驗(yàn)證了這一觀點(diǎn)。δ-變形菌常在硫循環(huán)以及鐵等金屬的代謝中扮演重要角色(Xu et al,2008),本次研究發(fā)現(xiàn) δ-變形菌克隆與原油泄漏海域的脫硫桿菌相似度高,脫硫桿菌是一種能降解有機(jī)物(尤其是石油烴類(lèi))的硫酸鹽還原細(xì)菌(Musat et al,2009),這預(yù)示著 δ-變形菌含量豐富的站點(diǎn)可能受到石油烴的污染。ε-變形菌雖然能廣泛利用環(huán)境中的氫、硫化物、硝酸鹽、氧氣等多種物質(zhì)作為電子受體或供體(Takai et al,2003;Nakagawa et al,2005),但在海水中的分布極少,本次調(diào)查也僅在A02站點(diǎn)有發(fā)現(xiàn)。

    擬桿菌門(mén)、放線菌門(mén)和藍(lán)細(xì)菌門(mén)是本次研究分布比較廣并且比例也較高的類(lèi)群。本次實(shí)驗(yàn)發(fā)現(xiàn)的擬桿菌主要包括黃桿菌綱和鞘脂桿菌綱兩大類(lèi),其中黃桿菌綱與藻類(lèi)的赤潮和降解過(guò)程有關(guān)(Abell et al,2005),而鞘脂桿菌綱重要類(lèi)群為噬胞菌屬(Cytophaga)可以降解纖維素(Li et al,2011),上述擬桿菌門(mén)含量較高與該海域初級(jí)生產(chǎn)力水平有關(guān)。放線菌是一類(lèi)在海洋中廣泛存在且豐度較高的類(lèi)群,有研究表明放線菌在有機(jī)污染物的降解過(guò)程當(dāng)中起著非常重要的作用(李和陽(yáng)等,2012),同時(shí)被普遍認(rèn)為是生物活性物質(zhì)和海洋藥物開(kāi)發(fā)的重要微生物來(lái)源(Bull et al,2005),這表明長(zhǎng)江口及鄰近海域可能蘊(yùn)藏著豐富的放線菌及相關(guān)生物活性產(chǎn)物資源。藍(lán)細(xì)菌是廣泛分布于海洋中能進(jìn)行放氧光合作用的超微型浮游植物(孫晟等,2003),是海洋初級(jí)生產(chǎn)力和微食物網(wǎng)中重要的組成成分。而在本次調(diào)查中,多個(gè)站點(diǎn)均發(fā)現(xiàn)藍(lán)細(xì)菌且其數(shù)量不少,可以看出藍(lán)細(xì)菌對(duì)該海域生態(tài)系統(tǒng)的穩(wěn)定性與多樣性有重要意義。

    作者在實(shí)驗(yàn)中有發(fā)現(xiàn)一些分布范圍和比例都較小的類(lèi)群,包括厚壁菌門(mén)、綠彎菌門(mén)、酸桿菌門(mén)、浮霉菌門(mén)、疣微菌門(mén)和纖維桿菌門(mén)。厚壁菌門(mén)中的芽孢桿菌主要存在于A01和A02站點(diǎn),Vieira等(2008)指出厚壁菌豐富的地方可能表明附近海區(qū)的污染,這可能與A01和A02兩個(gè)站點(diǎn)處在近河口位置有關(guān)。綠彎菌是除變形菌和硝化螺旋菌外,第三類(lèi)能通過(guò)硝酸鹽的氧化而獲得能量的一類(lèi)細(xì)菌(Sorokin et al,2012),作者發(fā)現(xiàn)綠彎菌僅存在于近海區(qū)域受陸源輸入影響而營(yíng)養(yǎng)鹽豐富的幾個(gè)站點(diǎn)當(dāng)中,因此推測(cè)綠彎菌傾向于生活在營(yíng)養(yǎng)鹽豐富的海域中。酸桿菌一般由外界環(huán)境所引入,A02站點(diǎn)位于長(zhǎng)江口受地表徑流影響較大,檢測(cè)到的酸桿菌很可能為陸源菌,另外Barns等(2007)曾報(bào)道過(guò)酸桿菌門(mén)類(lèi)群多存在于重金屬污染等酸性較強(qiáng)的環(huán)境中,這說(shuō)明長(zhǎng)江口區(qū)域可能存在一定程度酸化的傾向。浮霉菌門(mén)細(xì)菌屬于兼性厭氧細(xì)菌,并且這類(lèi)細(xì)菌可能與深海有機(jī)質(zhì)的礦化有關(guān)(S?rensen et al,2005),本研究結(jié)果顯示浮霉菌的相似序列來(lái)源于深海熱液區(qū)。

    3.3 長(zhǎng)江口及鄰近海域特殊細(xì)菌類(lèi)群及分布

    本次實(shí)驗(yàn)還發(fā)現(xiàn)了以往同海域研究中鮮有報(bào)道的柔膜菌門(mén)和SAR類(lèi)群。目前關(guān)于柔膜菌門(mén)的報(bào)道并不多,僅侯梅鋒等(2011)在青海湖底泥中有發(fā)現(xiàn),主要是一類(lèi)耐鹽菌,另有Beazley等(2012)在墨西哥灣原油泄漏的海岸鹽沼中發(fā)現(xiàn)大量柔膜菌的存在。本次研究在 C05站點(diǎn)發(fā)現(xiàn)少量的柔膜菌門(mén),具體原因不詳。SAR11是最早發(fā)現(xiàn)的一類(lèi)廣泛分布于亞熱帶海域的浮游細(xì)菌類(lèi)群,研究表明SAR11細(xì)菌能夠調(diào)控海洋中相當(dāng)數(shù)量DOM的流通,在C、N等元素的循環(huán)中起著重要的作用(Rappé et al,2002),另外海洋細(xì)菌生理代謝機(jī)制研究方面的重要發(fā)現(xiàn)之一是 SAR11菌群需要外源還原態(tài)硫來(lái)維持生長(zhǎng)(Tripp et al,2008)。對(duì)SAR406類(lèi)群的研究顯示,它是海洋最低含氧帶(Oxygen Minimum Zone,OMZ)中的優(yōu)勢(shì)類(lèi)群(Wright et al,2012),可能與缺氧環(huán)境有關(guān);Crump等(2007)在切薩皮克灣缺氧底層水中發(fā)現(xiàn)高豐度的 SAR406類(lèi)群,類(lèi)似于 OMZ的細(xì)菌群落構(gòu)成。本次研究發(fā)現(xiàn) A02站點(diǎn)同時(shí)有SAR11和SAR406類(lèi)群的分布,可能與該站點(diǎn)恰好位于長(zhǎng)江口外缺氧區(qū),有機(jī)質(zhì)和營(yíng)養(yǎng)鹽含量豐富有關(guān)。

    4 結(jié)論

    本實(shí)驗(yàn)通過(guò)DAPI染色計(jì)數(shù)和構(gòu)建16S rRNA基因文庫(kù)的方法對(duì)長(zhǎng)江口及鄰近海域10個(gè)站點(diǎn)表層海水中細(xì)菌的生物量及多樣性進(jìn)行研究。結(jié)果發(fā)現(xiàn)長(zhǎng)江口及鄰近海域表層海水中細(xì)菌的豐度總體較高,但各站點(diǎn)間存在一定差異。海水中細(xì)菌類(lèi)群主要包括變形菌門(mén)、擬桿菌門(mén)、放線菌門(mén)、藍(lán)細(xì)菌門(mén)、厚壁菌門(mén)、綠彎菌門(mén)、酸桿菌門(mén)、浮霉菌門(mén)、纖維桿菌門(mén)、疣微菌門(mén)和柔膜菌門(mén)等 11個(gè)類(lèi)群,其中 α-變形菌綱所占比例最高,其次為擬桿菌門(mén)和放線菌門(mén),其它類(lèi)群的分布比例和范圍都較小,同時(shí)還發(fā)現(xiàn)了廣受關(guān)注的SAR類(lèi)群和海洋環(huán)境中少見(jiàn)的柔膜菌門(mén)。研究表明某些細(xì)菌類(lèi)群可能成為該海域生態(tài)系統(tǒng)受污染的指示生物,此外未確定分類(lèi)地位的細(xì)菌可能是一些未被發(fā)現(xiàn)的細(xì)菌種類(lèi),有待進(jìn)一步調(diào)查研究。以上結(jié)果充分說(shuō)明了長(zhǎng)江口及鄰近海域表層海水細(xì)菌群落結(jié)構(gòu)具有豐富的多樣性和復(fù)雜性。

    王 新,李志江,鄭天凌,2010.海洋浮游細(xì)菌在東海赤潮高發(fā)區(qū)的分布與活性.環(huán)境科學(xué),31(2): 287—295

    寧修仁,史君賢,蔡昱明等,2004.長(zhǎng)江口和杭州灣海域生物生產(chǎn)力鋒面及其生態(tài)學(xué)效應(yīng).海洋學(xué)報(bào),26(6): 96—106

    劉 敏,王子峰,朱開(kāi)玲等,2008.應(yīng)用PCR-DGGE技術(shù)分析長(zhǎng)江口低氧區(qū)的細(xì)菌群落組成.高技術(shù)通訊,18(6):650—656

    劉晶晶,杜 萍,曾江寧等,2011.夏季長(zhǎng)江口浮游細(xì)菌和浮游病毒的分布特征及環(huán)境制約因素.海洋學(xué)研究,29(3):118—124

    孫 晟,肖 天,岳海東,2003.秋季與春季東、黃海藍(lán)細(xì)菌(Synechococcus spp.)生態(tài)分布特點(diǎn).海洋與湖沼,34(2):161—168

    李 濤,王 鵬,汪品先,2008.南海南部陸坡表層沉積物細(xì)菌和古菌多樣性.微生物學(xué)報(bào),48(3): 323—329

    李和陽(yáng),陳明霞,焦念志,2012.長(zhǎng)江口典型站位浮游細(xì)菌類(lèi)群的組成特征研究.海洋學(xué)報(bào),34(1): 183—188

    張東聲,2011.長(zhǎng)江口及其鄰近海域微生物的多樣性和生態(tài)分布特征研究.杭州: 浙江大學(xué)博士學(xué)位論文,13—15

    季 倩,2008.長(zhǎng)江口鄰近海域水體及底棲微微型浮游生物研究.上海: 華東師范大學(xué)碩士學(xué)位論文,93—95

    趙三軍,肖 天,岳海東,2003.秋季東、黃海異養(yǎng)細(xì)菌(Heterotrophic Bacteria)的分布特點(diǎn).海洋與湖沼,34(3):295—305

    趙海萍,李清雪,陶建華,2007.海洋細(xì)菌熒光顯微計(jì)數(shù)法及其應(yīng)用.河北工程大學(xué)學(xué)報(bào)(自然科學(xué)版),24(1): 57—60

    侯梅鋒,何士龍,李 棟等,2011.連云港海底底泥及青海湖底泥細(xì)菌多樣性研究.環(huán)境科學(xué),32(9): 2681—2688

    唐曉暉,王 凡,2004.長(zhǎng)江口鄰近海域夏,冬季水文特征分析.海洋科學(xué)集刊,46: 42—66

    曾永輝,2008.典型海洋環(huán)境中浮游細(xì)菌多樣性及環(huán)境適應(yīng)機(jī)制的研究.廈門(mén): 廈門(mén)大學(xué)博士學(xué)位論文,34—41

    潘洛安,2005.東海陸架及南海北部鋒區(qū)微生物群落結(jié)構(gòu)初探.上海: 華東師范大學(xué)碩士學(xué)位論文,109—113

    Abell G C J,Bowman J P,2005.Colonization and community dynamics of class Flavobacteria on diatom detritus in experimental mesocosms based on Southern Ocean seawater.FEMS Microbiology Ecology,53(3): 379—391

    Barns S M,Cain E C,Sommerville L et al,2007.Acidobacteria phylum sequences in uranium-contaminated subsurface sediments greatly expand the known diversity within the phylum.Applied and Environmental Microbiology,73(9): 3113—3116

    Beazley M J,Martinez R J,Rajan S et al,2012.Microbial community analysis of a coastal salt marsh affected by the Deepwater Horizon oil spill.PLoS One,7(7): e41305

    Bowman J P,McCammon S A,Dann A L,2005.Biogeographic and quantitative analyses of abundant uncultivated γ-proteobacterial clades from marine sediment.Microbial Ecology,49(3): 451—460

    Bowman J P,McCuaig R D,2003.Biodiversity,community structural shifts,and biogeography of prokaryotes within Antarctic continental shelf sediment.Applied and Environmental Microbiology,69(5): 2463—2483

    Bull A T,Stach J E M,Ward A C et al,2005.Marine actinobacteria: perspectives,challenges,future directions.Antonie van Leeuwenhoek,87(1): 65—79

    Cottrell M T,Waidner L A,Yu L et al,2005.Bacterial diversity of metagenomic and PCR libraries from the Delaware River.Environmental Microbiology,7(12): 1883—1895

    Crump B C,Peranteau C,Beckingham B et al,2007.Respiratory succession and community succession of bacterioplankton in seasonally anoxic estuarine waters.Applied and Environmental Microbiology,73(21): 6802—6810

    Feng B W,Li X R,Wang J H et al,2009.Bacterial diversity of water and sediment in the Changjiang estuary and coastal area of the East China Sea.FEMS Microbiology Ecology,70(2): 236—248

    Inagaki F,Suzuki M,Takai K et al,2003.Microbial communities associated with geological horizons in coastal subseafloor sediments from the Sea of Okhotsk.Applied and Environmental Microbiology,69(12): 7224—7235

    Jiao N Z,Zhang Y,Zeng Y H et al,2007.Ecological anomalies in the East China Sea: Impacts of the three gorges dam?.Water Research,41(6): 1287—1293

    Kumar S,Tamura K,Nei M,2004.MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment.Briefings in Bioinformatics,5(2): 150—163

    Li Y,Lee C G,Watanabe T et al,2011.Identification of microbial communities that assimilate substrate from root cap cells in an aerobic soil using a DNA-SIP approach.Soil Biology and Biochemistry,43(9): 1928—1935

    Lu J R,Santo Domingo J W,Lamendella R et al,2008.Phylogenetic diversity and molecular detection of bacteria in gull feces.Applied and Environmental Microbiology,74(13):3969—3976

    Mueller-Spitz S R,Goetz G W,McLellan S L,2009.Temporal and spatial variability in nearshore bacterioplankton communities of Lake Michigan.FEMS Microbiology Ecology,67(3): 511—522

    Musat F,Galushko A,Jacob J et al,2009.Anaerobic degradation of naphthalene and 2-methylnaphthalene by strains of marine sulfate-reducing bacteria.Environmental Microbiology,11(1):209—219

    Nakagawa S,Takai K,Inagaki F et al,2005.Distribution,phylogenetic diversity and physiological characteristics of epsilon-Proteobacteria in a deep-sea hydrothermal field.Environmental Microbiology,7(10): 1619—1632

    Ning X,Lin C,Su J et al,2011.Long-term changes of dissolved oxygen,hypoxia,and the responses of the ecosystems in the East China Sea from 1975 to 1995.Journal of Oceanography,67(1): 59—75

    Rappé M S,Connon S A,Vergin K L et al,2002.Cultivation of the ubiquitous SAR11 marine bacterioplankton clade.Nature,418(6898): 630—633

    Ravenschlag K,Sahm K,Pernthaler J et al,1999.High bacterial diversity in permanently cold marine sediments.Applied and Environmental Microbiology,65(9): 3982—3989

    Schloss P D,Handelsman J,2005.Introducing DOTUR,a computer program for defining operational taxonomic units and estimating species richness.Applied and Environmental Microbiology,71(3): 1501—1506

    Sekiguchi H,Koshikawa H,Hiroki M et al,2002.Bacterial distribution and phylogenetic diversity in the Changjiang estuary before the construction of the Three Gorges Dam.Microbial Ecology,43(1): 82—91

    Sj?stedt J,Koch-Schmidt P,Pontarp M et al,2012.Recruitment of members from the rare biosphere of marine bacterioplankton communities after an environmental disturbance.Applied and Environmental Microbiology,78(5): 1361—1369

    S?rensen K B,Canfield D E,Teske A P et al,2005.Community composition of a hypersaline endoevaporitic microbial mat.Applied and Environmental Microbiology,71(11): 7352—7365

    Sorokin D Y,Lücker S,Vejmelkova D et al,2012.Nitrification expanded: discovery,physiology and genomics of a nitrite-oxidizing bacterium from the phylum Chloroflexi.The ISME journal,6(12): 2245—2256

    Takai K,Inagaki F,Nakagawa S et al,2003.Isolation and phylogenetic diversity of members of previously uncultivated ε-Proteobacteria in deep-sea hydrothermal fields.FEMS Microbiology Letters,218(1): 167—174

    Tripp H J,Kitner J B,Schwalbach M S et al,2008.SAR11 marine bacteria require exogenous reduced sulphur for growth.Nature,452(7188): 741—744

    Vieira R P,Gonzalez A M,Cardoso A M et al,2008.Relationships between bacterial diversity and environmental variables in a tropical marine environment,Rio de Janeiro.Environmental Microbiology,10(1): 189—199

    V?r?s L,Balogh K V,Koncz E et al,2003.Phytoplankton and bacterioplankton production in a reed-covered water body.Aquatic Botany,77(2): 99—110

    Wright J J,Konwar K M,Hallam S J,2012.Microbial ecology of expanding oxygen minimum zones.Nature Reviews Microbiology,10(6): 381—394

    Xu H X,Wu M,Wang X G et al,2008.Bacterial diversity in deep-sea sediment from northeastern Pacific Ocean.Acta Ecologica Sinica,28(2): 479—485

    猜你喜歡
    長(zhǎng)江口文庫(kù)類(lèi)群
    專(zhuān)家文庫(kù)
    優(yōu)秀傳統(tǒng)文化啟蒙文庫(kù)
    幽默大師(2020年10期)2020-11-10 09:07:22
    關(guān)于推薦《當(dāng)代詩(shī)壇百家文庫(kù)》入選詩(shī)家的啟事
    薏苡種質(zhì)資源ISSR分子標(biāo)記篩選及親緣關(guān)系分析
    專(zhuān)家文庫(kù)
    黑石頂自然保護(hù)區(qū)土壤動(dòng)物功能類(lèi)群對(duì)季節(jié)動(dòng)態(tài)的響應(yīng)
    長(zhǎng)江口鳳鱭繁殖群體的年齡結(jié)構(gòu)和生長(zhǎng)特性
    長(zhǎng)江口及其鄰近水域仔稚魚(yú)種類(lèi)組成及分布特征
    2011年春夏季黃海和東海微型浮游動(dòng)物類(lèi)群組成及其攝食的研究
    風(fēng)場(chǎng)對(duì)夏季長(zhǎng)江口外東北海域低鹽水團(tuán)的作用
    特级一级黄色大片| 欧美区成人在线视频| av女优亚洲男人天堂| 七月丁香在线播放| 91av网一区二区| 黄片无遮挡物在线观看| 精品酒店卫生间| 亚洲成人av在线免费| 大话2 男鬼变身卡| 午夜日本视频在线| 久久久午夜欧美精品| 免费观看在线日韩| 亚洲人成网站在线观看播放| 久久精品国产鲁丝片午夜精品| 免费在线观看成人毛片| 日日干狠狠操夜夜爽| 舔av片在线| 欧美性猛交黑人性爽| 国产亚洲91精品色在线| 我要看日韩黄色一级片| 美女内射精品一级片tv| 亚洲精品影视一区二区三区av| 天美传媒精品一区二区| 国产精品爽爽va在线观看网站| 日韩精品青青久久久久久| 在线观看av片永久免费下载| 国产精品一区二区三区四区免费观看| 日韩,欧美,国产一区二区三区 | 人体艺术视频欧美日本| 国产白丝娇喘喷水9色精品| 中文字幕av在线有码专区| 啦啦啦观看免费观看视频高清| 午夜精品国产一区二区电影 | 联通29元200g的流量卡| 夜夜爽夜夜爽视频| 天天一区二区日本电影三级| 看黄色毛片网站| 免费人成在线观看视频色| 国产一级毛片在线| 干丝袜人妻中文字幕| 99热全是精品| 免费黄色在线免费观看| 一级毛片我不卡| 毛片女人毛片| 观看美女的网站| 99久久中文字幕三级久久日本| 精品久久久久久电影网 | 国产精品美女特级片免费视频播放器| 少妇裸体淫交视频免费看高清| 岛国在线免费视频观看| 亚洲欧美日韩卡通动漫| 亚洲av电影在线观看一区二区三区 | 精品一区二区三区人妻视频| 亚洲,欧美,日韩| 久久草成人影院| 国产老妇伦熟女老妇高清| 亚洲av免费高清在线观看| 精品一区二区三区人妻视频| 亚洲欧美精品综合久久99| 欧美成人一区二区免费高清观看| 欧美一区二区国产精品久久精品| 波多野结衣高清无吗| 欧美日本视频| 成年免费大片在线观看| 一夜夜www| 亚洲aⅴ乱码一区二区在线播放| 少妇熟女欧美另类| 欧美另类亚洲清纯唯美| 久久久精品大字幕| 国产午夜精品论理片| 人体艺术视频欧美日本| 一区二区三区高清视频在线| 深夜a级毛片| 国产成人91sexporn| 亚洲精品久久久久久婷婷小说 | 麻豆国产97在线/欧美| 青青草视频在线视频观看| 欧美性感艳星| 18禁裸乳无遮挡免费网站照片| 精品人妻熟女av久视频| 晚上一个人看的免费电影| 天堂网av新在线| 亚洲在久久综合| 成人高潮视频无遮挡免费网站| 在线观看66精品国产| 天天躁夜夜躁狠狠久久av| 深夜a级毛片| 亚洲欧洲日产国产| 男人和女人高潮做爰伦理| 级片在线观看| 夫妻性生交免费视频一级片| 久久热精品热| 黄片wwwwww| 国内揄拍国产精品人妻在线| 最新中文字幕久久久久| 亚洲综合色惰| av在线播放精品| 精品久久久久久久末码| 99热6这里只有精品| 亚洲精品乱码久久久v下载方式| 国产又黄又爽又无遮挡在线| 国产精品人妻久久久久久| 国产精品国产高清国产av| 欧美三级亚洲精品| 自拍偷自拍亚洲精品老妇| 麻豆国产97在线/欧美| 亚洲av电影在线观看一区二区三区 | 精品人妻熟女av久视频| 久久精品国产亚洲av天美| 国产午夜精品论理片| 国产一级毛片七仙女欲春2| 国产真实伦视频高清在线观看| 精品国产三级普通话版| 国产一区亚洲一区在线观看| 噜噜噜噜噜久久久久久91| 青青草视频在线视频观看| 一个人观看的视频www高清免费观看| 成人无遮挡网站| 插阴视频在线观看视频| 国产精华一区二区三区| 国产精品一区www在线观看| 插阴视频在线观看视频| 欧美色视频一区免费| 一本一本综合久久| 99久国产av精品国产电影| 99在线视频只有这里精品首页| 18禁裸乳无遮挡免费网站照片| 男插女下体视频免费在线播放| 六月丁香七月| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 99久久精品国产国产毛片| 久久久久网色| 欧美xxxx性猛交bbbb| 亚洲国产成人一精品久久久| 午夜精品一区二区三区免费看| 成人午夜精彩视频在线观看| 亚洲最大成人中文| 老司机影院毛片| 永久免费av网站大全| 日本与韩国留学比较| 五月伊人婷婷丁香| АⅤ资源中文在线天堂| 天美传媒精品一区二区| 一级毛片我不卡| av在线蜜桃| 国产午夜福利久久久久久| 亚洲精品亚洲一区二区| 日本与韩国留学比较| 国产精品一区二区三区四区久久| 九色成人免费人妻av| 听说在线观看完整版免费高清| 青春草亚洲视频在线观看| 岛国在线免费视频观看| 亚洲欧美清纯卡通| 国产麻豆成人av免费视频| 久久久a久久爽久久v久久| 欧美潮喷喷水| 免费一级毛片在线播放高清视频| 国产成人一区二区在线| 夜夜看夜夜爽夜夜摸| 欧美又色又爽又黄视频| 91aial.com中文字幕在线观看| 国产伦精品一区二区三区视频9| 免费不卡的大黄色大毛片视频在线观看 | 少妇猛男粗大的猛烈进出视频 | 亚洲人成网站在线观看播放| 欧美人与善性xxx| 欧美区成人在线视频| 中文天堂在线官网| 色噜噜av男人的天堂激情| 久久久久久久久大av| 成人特级av手机在线观看| 成人午夜高清在线视频| 99在线视频只有这里精品首页| 国产69精品久久久久777片| 亚洲国产精品国产精品| 亚洲人成网站高清观看| 成人午夜精彩视频在线观看| .国产精品久久| 久久久亚洲精品成人影院| 亚洲av.av天堂| 欧美zozozo另类| 亚洲精品国产av成人精品| 国产男人的电影天堂91| 99久久无色码亚洲精品果冻| 欧美极品一区二区三区四区| 18+在线观看网站| 久久精品人妻少妇| 色综合亚洲欧美另类图片| 久久久久九九精品影院| 成人美女网站在线观看视频| 久久久亚洲精品成人影院| 国产精品女同一区二区软件| 女的被弄到高潮叫床怎么办| 国产精品国产高清国产av| 久久久久久伊人网av| 国语对白做爰xxxⅹ性视频网站| 别揉我奶头 嗯啊视频| 一级黄色大片毛片| 黑人高潮一二区| 国产人妻一区二区三区在| 久久这里有精品视频免费| 亚洲精品亚洲一区二区| 亚洲精品国产成人久久av| 99热6这里只有精品| 不卡视频在线观看欧美| 免费av毛片视频| 亚洲婷婷狠狠爱综合网| av在线老鸭窝| 国产一级毛片七仙女欲春2| 国产一区二区在线av高清观看| 亚洲一区高清亚洲精品| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看 | 亚洲真实伦在线观看| 日韩成人av中文字幕在线观看| 91精品国产九色| 亚洲成人av在线免费| 色哟哟·www| 午夜激情福利司机影院| 国产伦理片在线播放av一区| 精华霜和精华液先用哪个| 亚洲国产精品成人综合色| 直男gayav资源| 看十八女毛片水多多多| 免费av毛片视频| 欧美色视频一区免费| 看片在线看免费视频| 亚洲欧美日韩无卡精品| 色播亚洲综合网| 精品人妻一区二区三区麻豆| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 女的被弄到高潮叫床怎么办| 熟女电影av网| 国产精品国产三级国产av玫瑰| av国产久精品久网站免费入址| 真实男女啪啪啪动态图| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 日韩强制内射视频| 少妇熟女欧美另类| 黄色一级大片看看| 婷婷色综合大香蕉| 久久久久九九精品影院| 秋霞伦理黄片| 亚洲精品自拍成人| 2021天堂中文幕一二区在线观| 听说在线观看完整版免费高清| 亚洲一区高清亚洲精品| h日本视频在线播放| 最近最新中文字幕免费大全7| 日韩成人av中文字幕在线观看| 直男gayav资源| 精品国产露脸久久av麻豆 | 午夜爱爱视频在线播放| 久久这里有精品视频免费| 日韩人妻高清精品专区| 精品免费久久久久久久清纯| av线在线观看网站| 亚洲av日韩在线播放| 久久久久久伊人网av| 午夜爱爱视频在线播放| 91狼人影院| 欧美bdsm另类| 免费观看性生交大片5| 韩国高清视频一区二区三区| 成人一区二区视频在线观看| 我的老师免费观看完整版| 国产精品一区二区三区四区久久| 亚洲欧美清纯卡通| 真实男女啪啪啪动态图| 人体艺术视频欧美日本| 美女国产视频在线观看| 能在线免费看毛片的网站| 精品人妻一区二区三区麻豆| 中文天堂在线官网| 国产精品精品国产色婷婷| 蜜桃亚洲精品一区二区三区| 国产高清不卡午夜福利| 日韩中字成人| 日本一二三区视频观看| 免费一级毛片在线播放高清视频| 又粗又爽又猛毛片免费看| av在线天堂中文字幕| 欧美极品一区二区三区四区| av天堂中文字幕网| 国产男人的电影天堂91| 亚洲欧美一区二区三区国产| 欧美日韩国产亚洲二区| eeuss影院久久| 亚洲欧美日韩卡通动漫| 高清在线视频一区二区三区 | 亚洲av男天堂| 日韩亚洲欧美综合| 国产伦精品一区二区三区四那| 蜜臀久久99精品久久宅男| 欧美激情在线99| 亚洲在线自拍视频| 亚洲精品久久久久久婷婷小说 | 中文在线观看免费www的网站| 五月伊人婷婷丁香| 91精品国产九色| 午夜福利在线在线| 日韩 亚洲 欧美在线| 小说图片视频综合网站| 亚洲国产精品成人综合色| 伦精品一区二区三区| 综合色av麻豆| 国产午夜精品一二区理论片| 国产一区二区亚洲精品在线观看| 亚洲丝袜综合中文字幕| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 国产单亲对白刺激| 日韩欧美精品v在线| 久久人人爽人人片av| 一区二区三区高清视频在线| 成人av在线播放网站| 国产午夜福利久久久久久| 一区二区三区免费毛片| 国语自产精品视频在线第100页| 免费观看人在逋| 亚洲精华国产精华液的使用体验| 一级毛片电影观看 | 精品一区二区免费观看| 亚洲图色成人| 草草在线视频免费看| 如何舔出高潮| 在线观看一区二区三区| 午夜福利成人在线免费观看| 中文欧美无线码| 国产麻豆成人av免费视频| 国产午夜福利久久久久久| 最后的刺客免费高清国语| 黄色配什么色好看| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 我要看日韩黄色一级片| 中国美白少妇内射xxxbb| 亚洲国产欧美人成| 人人妻人人看人人澡| 97在线视频观看| 成人欧美大片| 亚洲精品国产成人久久av| 免费一级毛片在线播放高清视频| a级毛片免费高清观看在线播放| 国产精品电影一区二区三区| 色哟哟·www| 欧美97在线视频| 亚洲国产欧美在线一区| 久久精品熟女亚洲av麻豆精品 | 中国国产av一级| 联通29元200g的流量卡| 国国产精品蜜臀av免费| 午夜视频国产福利| 成人漫画全彩无遮挡| 91久久精品电影网| 最近的中文字幕免费完整| 亚洲欧美精品自产自拍| 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 欧美色视频一区免费| 三级国产精品欧美在线观看| 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 男人和女人高潮做爰伦理| 国产一级毛片七仙女欲春2| 久久国内精品自在自线图片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产日韩欧美精品在线观看| 久久综合国产亚洲精品| 久久精品国产亚洲av涩爱| 99九九线精品视频在线观看视频| 亚洲性久久影院| av在线天堂中文字幕| av在线亚洲专区| 免费观看人在逋| 搞女人的毛片| 69av精品久久久久久| 国产白丝娇喘喷水9色精品| 女人十人毛片免费观看3o分钟| 天天一区二区日本电影三级| 毛片一级片免费看久久久久| 波野结衣二区三区在线| 国产精华一区二区三区| 午夜福利网站1000一区二区三区| 亚洲欧美日韩高清专用| 日韩一本色道免费dvd| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 99热6这里只有精品| 国产av在哪里看| 美女高潮的动态| 美女大奶头视频| 国产精品蜜桃在线观看| 搡女人真爽免费视频火全软件| av线在线观看网站| 亚洲精品色激情综合| 午夜精品一区二区三区免费看| 国产av不卡久久| 一级毛片aaaaaa免费看小| 色哟哟·www| 亚洲精华国产精华液的使用体验| 午夜亚洲福利在线播放| 白带黄色成豆腐渣| 美女cb高潮喷水在线观看| 久久久色成人| 我的老师免费观看完整版| 国产精品一区www在线观看| 国内少妇人妻偷人精品xxx网站| 少妇人妻精品综合一区二区| 国产精品福利在线免费观看| 国产91av在线免费观看| 99热网站在线观看| 久久久精品大字幕| 成年女人永久免费观看视频| 国产伦在线观看视频一区| 免费观看性生交大片5| 国产白丝娇喘喷水9色精品| 日韩欧美精品免费久久| 亚洲精品aⅴ在线观看| 中文字幕精品亚洲无线码一区| 日日啪夜夜撸| 久久久久性生活片| 免费观看a级毛片全部| 亚洲不卡免费看| 亚洲国产精品专区欧美| 大又大粗又爽又黄少妇毛片口| АⅤ资源中文在线天堂| 国产精品人妻久久久久久| 日韩成人av中文字幕在线观看| 伊人久久精品亚洲午夜| 亚洲五月天丁香| av福利片在线观看| 麻豆成人午夜福利视频| 只有这里有精品99| 免费黄网站久久成人精品| 最新中文字幕久久久久| 久久久久性生活片| 亚洲国产高清在线一区二区三| av视频在线观看入口| av专区在线播放| 国产成人免费观看mmmm| 国产亚洲一区二区精品| 国产精品一区二区性色av| 国模一区二区三区四区视频| 成人高潮视频无遮挡免费网站| 国产高清国产精品国产三级 | 免费观看精品视频网站| 日本午夜av视频| 亚洲成人久久爱视频| 亚洲av二区三区四区| 中文字幕熟女人妻在线| 亚洲在线观看片| 国产精品久久久久久精品电影| 日本黄色片子视频| 身体一侧抽搐| 一夜夜www| 日韩视频在线欧美| 亚洲美女视频黄频| 女人被狂操c到高潮| 国产女主播在线喷水免费视频网站 | 欧美又色又爽又黄视频| 国内精品宾馆在线| 在线观看66精品国产| 嫩草影院新地址| 欧美高清成人免费视频www| 国产成年人精品一区二区| 免费看av在线观看网站| 久久热精品热| 一区二区三区乱码不卡18| 国产伦一二天堂av在线观看| 国产高清视频在线观看网站| 精品99又大又爽又粗少妇毛片| 欧美又色又爽又黄视频| 1024手机看黄色片| 色网站视频免费| 99国产精品一区二区蜜桃av| 日韩欧美在线乱码| 国产精品三级大全| 在线播放国产精品三级| 91午夜精品亚洲一区二区三区| 能在线免费看毛片的网站| 欧美人与善性xxx| 91在线精品国自产拍蜜月| 一级黄色大片毛片| 最近的中文字幕免费完整| 亚洲成人中文字幕在线播放| 99热这里只有是精品50| 欧美性感艳星| 国产av码专区亚洲av| 人妻少妇偷人精品九色| 久久国产乱子免费精品| 日本一本二区三区精品| 嫩草影院精品99| 国产精品一区二区在线观看99 | 精品一区二区三区视频在线| 成人三级黄色视频| 丝袜喷水一区| 欧美性猛交╳xxx乱大交人| 国语自产精品视频在线第100页| 麻豆成人午夜福利视频| 婷婷色综合大香蕉| 欧美97在线视频| 观看美女的网站| 色网站视频免费| 三级毛片av免费| 少妇被粗大猛烈的视频| 成人美女网站在线观看视频| 日本免费a在线| 一级黄色大片毛片| 少妇熟女欧美另类| 久久久久久久久大av| 国产又色又爽无遮挡免| 好男人视频免费观看在线| 天天躁夜夜躁狠狠久久av| 亚洲av电影不卡..在线观看| 一边亲一边摸免费视频| 少妇熟女aⅴ在线视频| 久久久色成人| www日本黄色视频网| 国产成人精品一,二区| 在线观看美女被高潮喷水网站| 中文亚洲av片在线观看爽| 99热全是精品| 高清日韩中文字幕在线| 欧美日本视频| 最近手机中文字幕大全| 中文欧美无线码| 真实男女啪啪啪动态图| 免费大片18禁| av.在线天堂| 久久久国产成人精品二区| 国产私拍福利视频在线观看| 中文资源天堂在线| 日韩高清综合在线| 亚洲精品久久久久久婷婷小说 | 婷婷色综合大香蕉| 啦啦啦韩国在线观看视频| 亚洲av熟女| 国产精品综合久久久久久久免费| 午夜精品国产一区二区电影 | 午夜福利高清视频| 久久久久久久国产电影| 亚洲国产欧美人成| 老司机影院毛片| 人妻系列 视频| 免费电影在线观看免费观看| 高清视频免费观看一区二区 | 免费观看精品视频网站| 一二三四中文在线观看免费高清| 中文字幕av成人在线电影| 色尼玛亚洲综合影院| 中文字幕熟女人妻在线| 亚洲国产高清在线一区二区三| 久久久国产成人精品二区| 亚洲国产欧洲综合997久久,| 国产精品国产高清国产av| 亚洲成人久久爱视频| 日本wwww免费看| 在线播放无遮挡| 美女cb高潮喷水在线观看| 亚洲精品日韩av片在线观看| 久久精品熟女亚洲av麻豆精品 | 国产人妻一区二区三区在| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| 性色avwww在线观看| 只有这里有精品99| 亚洲在线观看片| 丝袜美腿在线中文| 久久欧美精品欧美久久欧美| 国产 一区 欧美 日韩| 搡老妇女老女人老熟妇| 三级国产精品欧美在线观看| 亚洲人成网站在线播| 丰满乱子伦码专区| 在线天堂最新版资源| 嘟嘟电影网在线观看| 精品久久国产蜜桃| 欧美日韩综合久久久久久| 又爽又黄无遮挡网站| 欧美精品国产亚洲| h日本视频在线播放| 免费黄色在线免费观看| 日韩欧美三级三区| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 在线免费观看的www视频| 亚洲在久久综合| 国产精品美女特级片免费视频播放器| 99热网站在线观看| 不卡视频在线观看欧美| 日韩 亚洲 欧美在线| 亚洲av.av天堂| 日韩国内少妇激情av| 亚洲人与动物交配视频| 国产精品,欧美在线| 日日摸夜夜添夜夜爱| 男女啪啪激烈高潮av片| 观看免费一级毛片| 两个人视频免费观看高清| 午夜日本视频在线| 国产成人a区在线观看| 久久韩国三级中文字幕| 欧美日本亚洲视频在线播放| 精品国内亚洲2022精品成人| 欧美成人a在线观看| 中文欧美无线码| 国产伦精品一区二区三区四那| 午夜免费激情av| 国产精品无大码|