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

    三份南海島礁珊瑚砂樣品中可培養(yǎng)細(xì)菌多樣性

    2022-03-31 01:02:50李存崔林青楊紅強(qiáng)龍麗娟田新朋
    熱帶海洋學(xué)報(bào) 2022年2期
    關(guān)鍵詞:放線菌島礁類(lèi)群

    李存, 崔林青, 楊紅強(qiáng), 龍麗娟, 田新朋

    1. 中國(guó)科學(xué)院熱帶海洋生物資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室(中國(guó)科學(xué)院南海海洋研究所), 廣東 廣州 510301;

    2. 中國(guó)科學(xué)院邊緣海與大洋地質(zhì)重點(diǎn)實(shí)驗(yàn)室(中國(guó)科學(xué)院南海海洋研究所), 廣東 廣州 510301;

    3. 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(廣州), 廣東 廣州 511458;

    4. 中國(guó)科學(xué)院大學(xué), 北京 100049

    珊瑚礁是典型的生物礁, 珊瑚礁生態(tài)系統(tǒng)是全球物種多樣性最高、資源最豐富的生態(tài)系統(tǒng)之一, 被譽(yù)為“海洋中的熱帶雨林”(龍麗娟 等,2019)。島礁生態(tài)系統(tǒng)中的微生物至關(guān)重要, 它不僅參與珊瑚礁生態(tài)系統(tǒng)的化學(xué)循環(huán)和物質(zhì)轉(zhuǎn)化,還與珊瑚礁生態(tài)系統(tǒng)中各種生物的健康息息相關(guān)(周進(jìn) 等, 2014), 但目前對(duì)島礁環(huán)境微生物的研究相對(duì)較少, 因此島礁環(huán)境中微生物資源有待深入的研究。

    微生物在地球能量流動(dòng)和物質(zhì)循環(huán)中扮演了重要角色, 也是珊瑚礁生態(tài)系統(tǒng)中重要成員, 開(kāi)展微生物多樣性研究將有助于了解微生物在珊瑚礁生態(tài)系統(tǒng)中的功能, 對(duì)珊瑚礁生態(tài)系統(tǒng)微生物資源的開(kāi)發(fā)和利用有促進(jìn)作用。近年來(lái)為了對(duì)純培養(yǎng)微生物多樣性進(jìn)行研究, 許多學(xué)者通過(guò)改進(jìn)培養(yǎng)基配方以及優(yōu)化培養(yǎng)條件, 顯著提高了分離得到的可培養(yǎng)微生物的多樣性, 并獲得了更多新的微生物類(lèi)群。孫創(chuàng)等(2021)采用改良的 2216E固體培養(yǎng)基、R2A 固體培養(yǎng)基、MBM 固體培養(yǎng)基、TCBS 固體培養(yǎng)基和改良的2216E 液體培養(yǎng)基對(duì)西太平洋海水中的微生物進(jìn)行分離純化, 而且獲得了多株新分類(lèi)單元。Xian 等(2020)通過(guò)微生物網(wǎng)絡(luò)分析, 發(fā)現(xiàn)溫單胞菌在群落中可以將復(fù)雜底物轉(zhuǎn)化為綠彎菌等難培養(yǎng)微生物生長(zhǎng)必需的小分子底物, 于是利用菌株TepidimonasSYSU G00190W 的上清液改良培養(yǎng)基, 實(shí)現(xiàn)了熱泉生境中未培養(yǎng)綠彎菌的定向分離培養(yǎng), 進(jìn)一步代謝組學(xué)研究表明, 培養(yǎng)基上清液中含有很多小分子有機(jī)基質(zhì), 可作為未培養(yǎng)細(xì)菌生長(zhǎng)的潛在營(yíng)養(yǎng)物質(zhì)。何媛秋等(2020)通過(guò)采用不同培養(yǎng)溫度、鹽度、pH、樣品稀釋倍數(shù)和營(yíng)養(yǎng)濃度條件對(duì)南海沉積物樣品進(jìn)行可培養(yǎng)細(xì)菌的多樣性研究, 發(fā)現(xiàn)寡營(yíng)養(yǎng)、20 倍的樣品稀釋倍數(shù)、pH 6 和0.05%以下的有機(jī)營(yíng)養(yǎng)物濃度更有利于稀有細(xì)菌類(lèi)群的分離培養(yǎng)。熊盈盈等(2021)總結(jié)前人對(duì)微生物分離培養(yǎng)的經(jīng)驗(yàn), 提出針對(duì)未培養(yǎng)的環(huán)境微生物培養(yǎng)基和培養(yǎng)條件改良方法, 比如添加抑制劑、樣品抽提液或菌液、信號(hào)分子、電子供受體以及維生素、生物酶或蛋白因子等物質(zhì), 或用其他凝固劑代替瓊脂,延長(zhǎng)培養(yǎng)時(shí)間以及極限稀釋培養(yǎng)法等。

    珊瑚礁樣品采集難度較大, 非常珍貴, 為獲得更多種類(lèi)的可培養(yǎng)微生物, 本研究采用5 種不同營(yíng)養(yǎng)成分的培養(yǎng)基, 并添加不同的微量元素, 對(duì)3 個(gè)島礁陸域珊瑚砂樣品進(jìn)行可培養(yǎng)微生物的分離; 同時(shí)分析了不同培養(yǎng)基條件下島礁可培養(yǎng)細(xì)菌的多樣性, 島礁可培養(yǎng)細(xì)菌的群落組成, 以及使用PICRUSt2 軟件對(duì)不同島礁微生物群落功能進(jìn)行了合理預(yù)測(cè), 以展示其生態(tài)功能狀況。

    1 材料與方法

    1.1 材料

    1.1.1 珊瑚砂樣品

    本次試驗(yàn)的南海島礁珊瑚砂樣品從有人島礁永暑礁(Y11-S)以及無(wú)人島礁扁參礁(BS2-S)和貝殼礁(BK-S)采集獲得, 使用無(wú)菌采樣勺取珊瑚砂表層0~3cm 的新鮮沉積砂樣, 裝入無(wú)菌塑料保藏袋, 放于-20℃保藏備用。

    1.1.2 分離培養(yǎng)基

    試驗(yàn)采用5 種寡營(yíng)養(yǎng)培養(yǎng)基, 即5% 2216E固體培養(yǎng)基、 5% AIA (DifocMactinomycete isolation agar)固體培養(yǎng)基、5% R2A 固體培養(yǎng)基、SN 固體培養(yǎng)基和10% 菌液瓊脂培養(yǎng)基(表1)。同時(shí)單獨(dú)添加不同的微量元素(終濃度為0.05%),分別為醋酸鈉、二甲基巰基丙酸(DMSP)、半胱氨酸、甲硫氨酸、?;撬帷C糠N培養(yǎng)基中瓊脂終濃度為15. 0g·L–1。

    表1 不同培養(yǎng)基的類(lèi)型及其成分Tab 1 Components of selective isolation media used in this study

    1.1.3 試驗(yàn)試劑與儀器

    細(xì)菌 16S rRNA 基因擴(kuò)增的通用引物 27F(5′-AGAGTTTGATCCTGGCTCAG-3′)和1492R (5′-TACGGCTACCTTGTTACGACTT-3′)由生工生物工程(上海)股份有限公司合成。5% Chelex 100 試劑(5g Chelex 100 Resin (Bio-rad, USA)溶解于100mL 滅菌超純水); Taq DNA 聚合酶(北京全式金); PCR 儀(Eppendorf, German); 凝膠成像系統(tǒng)(Bio-rad, USA);小型臺(tái)式高速離心機(jī)和恒溫箱(ESCO, Singapore)。

    1.2 菌株的分離、純化與保藏

    稱取0.5g 珊瑚砂樣品, 使用滅菌的無(wú)菌海水稀釋20 倍, 手動(dòng)用力震蕩搖動(dòng)3min 后取150μL 稀釋后的樣品懸濁液涂布于固體培養(yǎng)基表面, 每份樣品每種培養(yǎng)基涂布兩個(gè)平板。將涂布好的分離培養(yǎng)基倒置培養(yǎng)于28℃恒溫培養(yǎng)箱中, 連續(xù)培養(yǎng)一個(gè)月后挑取培養(yǎng)基中所有細(xì)菌形態(tài)的單菌落, 利用平板劃線法在2216E 固體培養(yǎng)基(MA)上進(jìn)行分離純化, 將純化后的菌種轉(zhuǎn)入體積質(zhì)量30%甘油管中, 放置于-80℃冰箱中保藏。

    1.3 16S rRNA 基因擴(kuò)增測(cè)序分析

    細(xì)菌總 DNA 通過(guò) 5% Chelex-100 試劑提取(Walsh et al, 1991), 采用細(xì)菌通用引物27F 和1492R進(jìn)行16S rRNA 基因PCR 擴(kuò)增。獲得的PCR 產(chǎn)物經(jīng)瓊脂糖凝膠電泳檢測(cè)后, 送廣州天一輝遠(yuǎn)測(cè)序公司進(jìn)行基因測(cè)序。測(cè)序結(jié)果采用 SeqMan 5.0(Swindell et al, 1997)軟件進(jìn)行拼接, 去除質(zhì)量不佳的測(cè)序堿基后, 將獲得的有效序列(近1500bp)提交至EzBioCloud (Yoon et al, 2017)和NCBI (https://www.ncbi.nlm.nih.gov/)進(jìn)行在線BLAST 比對(duì), 比對(duì)結(jié)果按照 98%作為區(qū)分物種的標(biāo)準(zhǔn)(Yarza et al,2014), 若與已知菌種相似度小于等于98%, 則認(rèn)為該菌株為潛在新種。

    1.4 系統(tǒng)進(jìn)化分析

    根據(jù)序列比對(duì)結(jié)果, 選取同源性較高的模式菌株的 16S rRNA 基因序列作為參比序列, 使用CLUSTAL_W 軟件進(jìn)行多序列比對(duì), 使用MEGA-X(Kumar et al, 2018)軟件, 結(jié)合Kimura2-parameter 模型估算系統(tǒng)進(jìn)化矩陣, 采用鄰近法(neighbor-joining,N-J)(Saitou et al, 1987)進(jìn)行聚類(lèi)分析, 設(shè)置1000 次重采樣構(gòu)建系統(tǒng)發(fā)育樹(shù), 使用iTOL (Letunic et al,2019)在線網(wǎng)站對(duì)系統(tǒng)發(fā)育樹(shù)進(jìn)行可視化。

    1.5 多樣性分析與群落功能預(yù)測(cè)

    統(tǒng)計(jì)3 份樣品中獲得的細(xì)菌種類(lèi)和豐度, 制作樣品-細(xì)菌物種豐度矩陣表, 使用R 語(yǔ)言vegan 包計(jì)算樣品的α多樣性, 包括Shannon 多樣性指數(shù)、Simpson 指數(shù)、invsimpson 多樣性指數(shù)和Chao1 指數(shù)。使用PICRUSt2 軟件(Douglas et al, 2020)基于16S rRNA 基因數(shù)據(jù)對(duì)群落功能進(jìn)行預(yù)測(cè), 參考京都基因與基因組百科全書(shū)(kyoto encyclopedia of genes and genomes, KEGG)數(shù)據(jù)庫(kù), 得到直系同源(KEGG orthology)功能的豐度預(yù)測(cè)表及KEGG 代謝途徑豐度表。

    2 結(jié)果與分析

    2.1 細(xì)菌多樣性

    通過(guò)菌株分離純化, 本研究共從126 個(gè)未污染的分離平板獲得349 株細(xì)菌, 經(jīng)16S rRNA 基因序列相似性比對(duì)分析, 顯示這些新分離的菌株分別與目前134 個(gè)已知種親緣關(guān)系相近, 這些菌株分別屬于4 門(mén)、6 綱、26 目、43 科、73 屬。4 個(gè)分離到的細(xì)菌門(mén)中放線菌門(mén)(Actinobacteria)的菌株占所有分離菌株的60%, 為優(yōu)勢(shì)菌門(mén), 其他依次為變形菌門(mén)(Proteobacteria, 28%)、厚壁菌門(mén)(Firmicutes, 11%)和擬桿菌門(mén)(Bacteroidetes, 2%)。分離到較多的屬級(jí)類(lèi)群為鏈霉菌屬(Streptomyces, 46 株)、糖多孢菌屬(Saccharopolyspora, 40 株) 和產(chǎn)微球莖菌屬(Microbulbifer, 37 株)。圖1 展示了純培養(yǎng)菌株在不同分類(lèi)等級(jí)中細(xì)菌類(lèi)群所占的比例。圖2 展示了134 個(gè)物種(選取同種菌株中相似度最高的作為代表菌株)與其最相似已知物種菌株的16S rRNA 基因組列構(gòu)建的系統(tǒng)發(fā)育關(guān)系。對(duì)獲得的4 個(gè)門(mén)級(jí)類(lèi)群純培養(yǎng)菌株的多樣性分別進(jìn)行描述。

    圖1 分離獲得的純培養(yǎng)細(xì)菌在各分類(lèi)等級(jí)中的分布情況a. 在不同門(mén)級(jí)、綱級(jí)、目級(jí)水平下細(xì)菌類(lèi)群比例; b. 在不同屬水平下細(xì)菌類(lèi)群比例(標(biāo)注了菌株中豐度前25 的屬)Fig. 1 Distribution of isolated pure culture strains in each classification grade. (a) Proportion of pure culture strains at different phylum, class and order levels, and (b) proportion of pure culture strains in the top 25 genera

    圖2 采用鄰接法基于16S rRNA 基因序列構(gòu)建的134 個(gè)代表菌株系統(tǒng)發(fā)育樹(shù)Fig. 2 Neighbor-joining phylogenetic tree of 134 representative strains reconstructed based on the 16S rRNA gene sequences

    2.1.1 放線菌門(mén)

    分離得到的放線菌門(mén)細(xì)菌共計(jì)210 株, 屬于放線菌綱(Actinomycetia)的15 目22 科37 屬74 個(gè)菌種, 分離菌株數(shù)量最多的為鏈霉菌和糖多孢菌屬, 其中包含鏈霉菌屬下的11 個(gè)菌種、糖多孢菌屬下的3 個(gè)菌種, 其中Saccharopolyspora spongiae(27 株)、Streptomycesdiacarni(20 株)、Saccharopolyspora tripterygii(12 株)、Streptomyces albiaxialis(11 株) 、Nocardioides iriomotensis(11 株)為分離到的主要優(yōu)勢(shì)菌株, 其他優(yōu)勢(shì)屬還包括微桿菌屬(Microbacterium)、類(lèi)諾卡氏屬(Nocardioides) 、Barrientosiimonas、 布勞氏菌屬(Prauserella)、短桿菌屬(Brevibacterium)、迪茨氏菌屬(Dietzia)、兩面神菌屬(Janibacter)、羅斯氏菌屬(Rothia)、紅球菌屬(Rhodococcus)、戈登氏菌屬(Gordonia)等。

    2.1.2 變形菌門(mén)

    本研究分離得到的變形菌門(mén)菌株共計(jì) 96 株,涵蓋了3 綱(α-變形菌綱、β-變形菌綱、γ-變形菌綱)、8 目(Burkholderiales、Caulobacterales、Cellvibrionales、Lysobacterales 、 Rhizobiales 、 Rhodobacterales 、Rhodospirillales、Sphingomonadales)、15 科、23 屬、33 種。在綱水平上優(yōu)勢(shì)類(lèi)群為α-變形菌綱(53 株)和γ- 變形菌綱(42 株) 。 前者中以紅細(xì)菌目(Caulobacterales)、根瘤菌目(Rhizobiales)和柄桿菌目(Caulobacterales)為主, 后者則以Cellvibrionales 為主。獲得菌株中最多的物種為Microbulbifer variabilis(34 株)、Ruegeria arenilitoris(11 株)、Pararhizobium haloflavum(9 株)。

    2.1.3 厚壁菌門(mén)和擬桿菌門(mén)

    本研究中獲取的37 株厚壁菌門(mén)菌株分別屬于芽孢桿菌目的4 科11 屬, 以芽孢桿菌屬(Bacillus)、葡萄球菌屬(Staphylococcus) 、 短芽孢桿菌屬(Brevibacillus)為主, 優(yōu)勢(shì)菌株為Bacillus altitudinis(5 株)。另外, 擬桿菌門(mén)僅分離得到6 株, 分別屬于土壤中華微桿菌(Sinomicrobium soli, 3 株)和黃色簡(jiǎn)浩然菌(Haoranjiania flava, 3 株)。

    2.2 3 份珊瑚砂樣品中可培養(yǎng)細(xì)菌多樣性的比較

    本研究所采集的3 份樣品中可培養(yǎng)細(xì)菌的主要優(yōu)勢(shì)類(lèi)群不同, 選擇Shannon、Simpson、Invsimpson、Chao1 和ACE 指數(shù)對(duì)3 份樣品的細(xì)菌群落多樣性進(jìn)行評(píng)估。由表2 和圖3 可知, 樣品BK-S 和BS2-S的多樣性較高, 而Y11-S 樣品的多樣性較為單一。BK-S 樣品中獲得的80 個(gè)菌種分別屬于4 門(mén)、5 綱、21 目、32 科、48 屬, 優(yōu)勢(shì)類(lèi)群為放線菌綱、α-變形菌綱和γ-變形菌綱。BS2-S 樣品中獲得的64 個(gè)物種分別屬于4 門(mén)、7 綱、20 目、32 科、43 屬, 優(yōu)勢(shì)類(lèi)群為放線菌綱、α-變形菌綱和芽孢桿菌綱。Y11-S樣品中獲得的20 個(gè)物種分別屬于2 門(mén)、2 綱、6 目、7 科、10 屬, 其優(yōu)勢(shì)類(lèi)群為放線菌綱的鏈霉菌屬。雖然采集3 個(gè)樣品的水深均為10m 左右, 采集樣品的水溫、透光性等環(huán)境參數(shù)也相似, 但樣品可培養(yǎng)細(xì)菌多樣性存在明顯差異, 樣品間細(xì)菌多樣性差異是否與人類(lèi)活動(dòng)有關(guān), 仍需進(jìn)一步驗(yàn)證。

    表2 3 份珊瑚砂樣品中細(xì)菌群落多樣性指數(shù)Tab. 2 Diversity index of pure cultured strains isolated from three coral reef samples

    圖3 3 份樣品(Y11、BK-S、BS2-S)中共有菌種數(shù)量Fig. 3 Numbers of same species discovered in three samples of Y11-S, BK-S and BS2-S

    2.3 不同培養(yǎng)基分離到的可培養(yǎng)細(xì)菌多樣性比較

    本研究使用5 種培養(yǎng)基獲得的最優(yōu)勢(shì)門(mén)級(jí)類(lèi)群都為放線菌門(mén), 最優(yōu)勢(shì)綱為放線菌綱。在5 種培養(yǎng)基中, 5% 2216、SN、10%菌液瓊脂培養(yǎng)基分離得到的細(xì)菌多樣性較高, 除放線菌綱外, 3 種培養(yǎng)基中都分離到了芽孢桿菌綱、α-變形菌綱和γ-變形菌綱。在SN 培養(yǎng)基中還獲得了3 株黃桿菌綱菌種, 在10%菌液瓊脂培養(yǎng)基中還獲得了3 株Chitinophagia 綱和1 株β-變形菌綱菌種。在這種水平上, 5 種培養(yǎng)基分離得到細(xì)菌的多樣性由高到低依次為5% 2216 培養(yǎng)基(61 種)、10%菌液瓊脂培養(yǎng)基(59 種)、SN 培養(yǎng)基(53 種)、5% AIA 培養(yǎng)基(9 種)、5% R2A 培養(yǎng)基(8種)(圖4)。5% 2216 培養(yǎng)基和10%菌液瓊脂培養(yǎng)基中含有豐富多樣的營(yíng)養(yǎng)物質(zhì), 其營(yíng)養(yǎng)成分配比更接近于海洋環(huán)境, 如碳氮硫濃度、鹽度, 可以為多種類(lèi)型的菌株提供生長(zhǎng)所需的營(yíng)養(yǎng)環(huán)境。分析添加的微量元素的分離效果, 含有醋酸鈉的培養(yǎng)基獲得了更多的菌株數(shù)量; 同時(shí)發(fā)現(xiàn)與已知物種Pseudactinotalea terrae、Actinoplanes sediminis、Tessaracoccus flavescens、Achromobacter piechaudii等相近的潛在新分類(lèi)單元僅在含有醋酸鈉的分離培養(yǎng)基中獲得, 與Micromonospora aurantiaca相近的潛在新分類(lèi)單元僅在含有甲硫氨酸的分離培養(yǎng)基中獲得, 與Actinomadura meyerae相近的潛在新分類(lèi)單元僅在含有DMSP 的培養(yǎng)基中獲得。

    圖4 不同培養(yǎng)基獲得的可培養(yǎng)細(xì)菌在不同分類(lèi)水平上的多樣性Fig. 4 Bacterial diversity in different classification grades isolated from different selective isolation media

    2.4 潛在新物種信息

    以16S rRNA 基因相似度大于98% 時(shí)作為同一物種, 本試驗(yàn)中共分離得到56 株潛在新種, 與目前有效發(fā)表的最相近菌種的16S rRNA 基因相似度范圍在89.9%~98.0%, 其中89.9%~95.0% 的有7 株,95.0%~97.0%的有21 株, 97.0%~98.0%的有28 株;去重復(fù)后為18 個(gè)潛在新種, 見(jiàn)表3。圖5 為56 株潛在新種與已知物種的系統(tǒng)進(jìn)化發(fā)育樹(shù), 其中31 株潛在新分類(lèi)單元屬于放線菌綱, 19 株屬于α-變形菌綱, 3 株屬于芽孢桿菌綱, 3 株屬于Chitinophagia 綱。獲得的31 株潛在放線菌綱新種分別屬于短桿菌屬(Brevibacterium)、短狀桿菌屬(Brachybacterium)、Barrientosiimonas屬 、 賴 氨 酸 單 胞 菌 屬(Lysinimonas)、微桿菌屬(Microbacterium)、花金龜菌屬(Protaetiibacter)、迪茨氏菌屬(Dietzia)、戈登氏菌屬(Gordonia)、類(lèi)諾卡氏菌屬(Nocardioides)、普勞氏菌屬(Prauserella)、阮氏菌屬(Ruania)和鏈霉菌屬(Streptomyces)。對(duì)潛在新種的來(lái)源培養(yǎng)基進(jìn)行分析,發(fā)現(xiàn)在5% 2216 培養(yǎng)基(19 株)和SN 培養(yǎng)基(19 株)中獲得了最多的潛在新種。圖6 展示了部分潛在新物種的平板菌落形態(tài)情況。

    圖5 不同培養(yǎng)基分離得到的潛在新種的系統(tǒng)發(fā)育樹(shù)Fig. 5 Phylogenetic tree of potential new taxa obtained from different media

    圖6 不同潛在新種在2216E 培養(yǎng)基上培養(yǎng)14 天后的菌落平板a. SCSIO 63409 (Inquilinus limosus); b. SCSIO 63263 (Brevibacterium casei); c. SCSIO 63671 (Mesorhizobium oceanicum); d. SCSIO 63659 (Nocardioides iriomotensis); e. SCSIO 63385 (Gordonia terrae); f. SCSIO 63505 (Cohnella cellulosilytica); g. SCSIO 63658(Prauserella rugose); h. SCSIO 68054 (Oceanobacillus iheyensis)Fig. 6 Pure culture morphology of potential new taxa on 2216E medium after incubation for 14 days. (a) SCSIO 63409(Inquilinus limosus); (b) SCSIO 63263 (Brevibacterium casei); (c) SCSIO 63671 (Mesorhizobium oceanicum); (d) SCSIO 63659 (Nocardioides iriomotensis); (e) SCSIO 63385 (Gordonia terrae); (f) SCSIO 63505 (Cohnella cellulosilytica); (g)SCSIO 63658 (Prauserella rugose); (h) SCSIO 68054 (Oceanobacillus iheyensis)

    表3 潛在新菌種去重復(fù)信息表Tab. 3 Potential new taxa information after removing the duplication

    2.5 PICRUSt2 功能預(yù)測(cè)分析

    基于PICRUSt2 軟件, 進(jìn)行菌群功能分析, 得到不同樣品細(xì)菌的功能預(yù)測(cè)信息。利用KEGG 數(shù)據(jù)庫(kù)對(duì)測(cè)序數(shù)據(jù)進(jìn)行比對(duì), 所有樣品分析的功能涉及8 類(lèi)生物代謝通路: 代謝(metabolism)、功能蛋白(protein families)、遺傳信息處理(genetic information processing)、環(huán)境信息處理(environmental information processing)、細(xì)胞轉(zhuǎn)化(cellular processes)、生物體系統(tǒng)(organismal systems)、人類(lèi)疾病(human diseases)和未分類(lèi)(Unclassified)(圖7), 這8 個(gè)代謝通路主要體現(xiàn)了細(xì)菌菌群的生長(zhǎng)繁殖、物質(zhì)代謝與轉(zhuǎn)運(yùn)、能量代謝與傳遞、環(huán)境互作等方面的功能信息。分析結(jié)果顯示, 代謝和編碼功能蛋白是主要通路, 豐度分別為40.64%~43.0%和23.18%~23.83%。進(jìn)一步對(duì)預(yù)測(cè)基因二級(jí)功能層級(jí)進(jìn)行分析, 共發(fā)現(xiàn)碳水化合物代謝(carbohydrate metabolism)、氨基酸代謝(amino acid metabolism)、轉(zhuǎn)錄(transcription)、膜轉(zhuǎn)運(yùn) (membrane transport) 、 信 號(hào) 傳 導(dǎo) (signal transduction)、細(xì)胞生長(zhǎng)和死亡(cell growth and death)和能量代謝(energy metabolism)等38 個(gè)子功能。分析發(fā)現(xiàn)3 個(gè)樣品的群落功能組成大體一致, 在碳水化合物代謝、氨基酸代謝、編碼遺傳信息處理蛋白、編碼信號(hào)轉(zhuǎn)導(dǎo)和細(xì)胞過(guò)程蛋白等通路具有更高的豐度。PICRUSt2 功能預(yù)測(cè)分析顯示出有人島礁和無(wú)人島礁細(xì)菌群體在其生態(tài)功能方面并無(wú)差異, 這說(shuō)明島礁周邊水下環(huán)境細(xì)菌生態(tài)功能并未因?yàn)閸u礁陸域上人為活動(dòng)而受到影響。

    圖7 不同來(lái)源樣品的PICRUSt2 功能預(yù)測(cè)熱圖Fig. 7 Function prediction heatmap of PICRUSt2 for three different samples

    3 討論

    珊瑚砂是以珊瑚碎屑為主并有石灰藻、有孔蟲(chóng)、棘皮動(dòng)物碎片、貝殼碎片等組成的砂質(zhì)沉積物, 顆粒大小不等, 規(guī)格多種多樣, 其最大的特點(diǎn)是表面微孔豐富, 適宜細(xì)菌生存。本研究采用純培養(yǎng)方式對(duì)南海島礁珊瑚砂環(huán)境可培養(yǎng)細(xì)菌多樣性進(jìn)行研究,共獲得349 株純培養(yǎng)菌株。鑒定結(jié)果顯示這些菌株分別屬于4 門(mén)、6 綱、26 目、43 科、73 屬、134 種,包括18 個(gè)潛在新種, 其中放線菌門(mén)是最優(yōu)勢(shì)的細(xì)菌類(lèi)群, 占所有分離菌株的60%, 原因可能與3 個(gè)樣品營(yíng)養(yǎng)貧瘠有關(guān); 其他優(yōu)勢(shì)類(lèi)群依次為變形菌門(mén)(28%)、厚壁菌門(mén)(11%)和擬桿菌門(mén)(2%)。分離得到較多的屬級(jí)類(lèi)群為鏈霉菌屬(46 株)、糖多孢菌屬(40株)和產(chǎn)微球莖菌屬(37 株)。分析發(fā)現(xiàn)樣品BK-S 和BS2-S 中的細(xì)菌多樣性較高, 而Y11 樣品的細(xì)菌多樣性較為單一。樣品BK-S 和BS2-S 來(lái)源于無(wú)人島礁, Y11 樣品來(lái)源的島礁上有人類(lèi)活動(dòng), 對(duì)細(xì)菌群落功能的預(yù)測(cè)分析結(jié)果顯示有人島礁和無(wú)人島礁細(xì)菌群體在其生態(tài)功能方面并無(wú)差異, 這說(shuō)明島礁周邊水下環(huán)境細(xì)菌生態(tài)功能并未因?yàn)閸u礁陸域上人為活動(dòng)而受到影響, 因此Y11 樣品來(lái)源島礁周?chē)h(huán)境中的細(xì)菌多樣性低的原因仍然未知。Sch?ttner 等(2011)對(duì)紅海亞喀巴灣東北部淺岸礁的碳酸鹽砂和硅酸鹽砂樣品進(jìn)行免培養(yǎng)細(xì)菌多樣性的研究, 發(fā)現(xiàn)珊瑚砂樣品中細(xì)菌的群落組成主要為γ-變形菌綱、放線菌綱、α-變形菌綱、δ-變形菌綱、厚壁菌門(mén)、擬桿菌門(mén)等。韓敏敏等(2020)使用6 種普適性培養(yǎng)基對(duì)Khai島和Pathiu 島采集的2 份珊瑚礁沉積物樣品進(jìn)行分離培養(yǎng), 獲得的可培養(yǎng)細(xì)菌類(lèi)群為變形菌門(mén)(54.55%)、厚壁菌門(mén)(31.82%)、放線菌門(mén)(9.09%)和擬桿菌門(mén)(4.54%)。對(duì)比已有的對(duì)珊瑚砂細(xì)菌多樣性的研究, 本研究獲得的可培養(yǎng)類(lèi)群與珊瑚砂中主要細(xì)菌類(lèi)群一致, 并且通過(guò)優(yōu)化改良培養(yǎng)基, 從樣品中獲取了更多種類(lèi)的菌株, 有助于更好的反映珊瑚砂樣品中可培養(yǎng)細(xì)菌的多樣性。

    本研究將常用的2216E、R2A、AIA 等基礎(chǔ)培養(yǎng)基稀釋到5%, 以及使用SN、10%菌液瓊脂等寡營(yíng)養(yǎng)培養(yǎng)基對(duì)樣品進(jìn)行分離, 從得到的純培養(yǎng)菌株多樣性來(lái)看, 寡營(yíng)養(yǎng)培養(yǎng)基珊瑚砂樣品具有更好的分離效果。根據(jù)文獻(xiàn)報(bào)道, 在營(yíng)養(yǎng)豐富的培養(yǎng)基中獲得的微生物的數(shù)量并非最多, 相反培養(yǎng)基中高濃度的營(yíng)養(yǎng)物會(huì)影響寡營(yíng)養(yǎng)微生物細(xì)胞生長(zhǎng), 而用于培養(yǎng)寡營(yíng)養(yǎng)環(huán)境中微生物的培養(yǎng)基中有機(jī)碳水平不應(yīng)高于50mg·L-1(Jannasch et al, 1959; Kadota, 1981;Kjelleberg et al, 1985)。Cho 等(2004)在分析太平洋及遠(yuǎn)洋樣品時(shí), 發(fā)現(xiàn)有44 株γ-變形菌無(wú)法在經(jīng)典培養(yǎng)基上形成菌落, 而其中7 株菌卻可以在10% R2A 瓊脂培養(yǎng)基上形成菌落。K o c h(2 0 0 1)認(rèn)為當(dāng)寡營(yíng)養(yǎng)微生物從營(yíng)養(yǎng)貧乏的生態(tài)環(huán)境中突然轉(zhuǎn)移到富營(yíng)養(yǎng)環(huán)境中時(shí), 微生物無(wú)法快速適應(yīng)環(huán)境變化, 可能會(huì)在生長(zhǎng)過(guò)程中產(chǎn)生有毒物質(zhì), 例如有氧代謝形成的有毒產(chǎn)物累積或細(xì)胞內(nèi)自由基過(guò)快積累, 導(dǎo)致細(xì)胞損傷, 影響生長(zhǎng)。同時(shí), 大量非代謝物質(zhì)的轉(zhuǎn)運(yùn)攝取會(huì)導(dǎo)致細(xì)胞壁、細(xì)胞質(zhì)生長(zhǎng)不均衡, 膜蛋白比例過(guò)高阻礙細(xì)胞質(zhì)膜生長(zhǎng), 細(xì)胞內(nèi)代謝失調(diào), 影響細(xì)胞增殖機(jī)制。

    環(huán)境中有些微生物具有乙酰輔酶 A 連接酶(Acetate-CoA ligase, EC 6.2.1.1), 能利用乙酸鹽作為唯一碳源, 使乙酸鹽生成乙酰輔酶A, 而乙酰輔酶A 是細(xì)胞生化代謝過(guò)程中重要的中間產(chǎn)物, 可以轉(zhuǎn)化成多種細(xì)胞所需要的營(yíng)養(yǎng)物質(zhì)。硫元素是生物體必需的營(yíng)養(yǎng)元素, 在海洋生態(tài)系統(tǒng)中硫酸鹽、硫化物儲(chǔ)存量巨大, 不同價(jià)態(tài)硫化合物之間的轉(zhuǎn)化主要由代謝功能多樣的微生物等來(lái)完成, 環(huán)境中的硫元素對(duì)某些微生物的生長(zhǎng)至關(guān)重要。因此, 本研究在寡營(yíng)養(yǎng)培養(yǎng)基的基礎(chǔ)上, 添加乙酸鈉、含硫氨基酸、DMSP 等小分子物質(zhì)來(lái)增加從環(huán)境中獲取微生物的幾率, 從結(jié)果來(lái)看, 與以往使用經(jīng)典培養(yǎng)基進(jìn)行環(huán)境樣品分離相比, 本研究獲得的菌株在多樣性和潛在新種數(shù)量等方面都有顯著的提高。

    猜你喜歡
    放線菌島礁類(lèi)群
    肺曲霉菌合并肺放線菌感染一例
    薏苡種質(zhì)資源ISSR分子標(biāo)記篩選及親緣關(guān)系分析
    體系作戰(zhàn)條件下島礁作戰(zhàn)中輔助決策問(wèn)題研究
    黑石頂自然保護(hù)區(qū)土壤動(dòng)物功能類(lèi)群對(duì)季節(jié)動(dòng)態(tài)的響應(yīng)
    基于OODA過(guò)程的島礁防空CGF模型
    南大西洋深海沉積物中可培養(yǎng)放線菌的多樣性
    2011年春夏季黃海和東海微型浮游動(dòng)物類(lèi)群組成及其攝食的研究
    黃花蒿內(nèi)生放線菌A5次生代謝產(chǎn)物分離鑒定
    近35年來(lái)熱帶風(fēng)暴對(duì)我國(guó)南海島礁的影響分析
    龍灣自然保護(hù)區(qū)森林土壤動(dòng)物群落多樣性及功能類(lèi)群
    成人国产综合亚洲| 不卡av一区二区三区| 久久久国产欧美日韩av| 午夜免费鲁丝| 亚洲国产中文字幕在线视频| 老司机午夜福利在线观看视频| 欧美色视频一区免费| 国产一区二区激情短视频| 此物有八面人人有两片| 美女国产高潮福利片在线看| 欧美大码av| 两人在一起打扑克的视频| 国产精品久久久久久人妻精品电影| 国产精品二区激情视频| 黄色女人牲交| 亚洲真实伦在线观看| 日韩成人在线观看一区二区三区| 亚洲国产欧美日韩在线播放| 久久精品亚洲精品国产色婷小说| av电影中文网址| 在线观看免费日韩欧美大片| 露出奶头的视频| 色综合亚洲欧美另类图片| 色哟哟哟哟哟哟| 欧美乱码精品一区二区三区| 国产成人精品久久二区二区免费| 欧美黄色淫秽网站| a级毛片a级免费在线| 久久国产乱子伦精品免费另类| 亚洲五月色婷婷综合| 国产伦在线观看视频一区| 50天的宝宝边吃奶边哭怎么回事| 久久久精品国产亚洲av高清涩受| 国产高清有码在线观看视频 | 成人亚洲精品一区在线观看| 国产野战对白在线观看| 亚洲av五月六月丁香网| 欧美性长视频在线观看| 色综合欧美亚洲国产小说| 精品无人区乱码1区二区| 国产精品久久视频播放| 黄频高清免费视频| 90打野战视频偷拍视频| 一二三四社区在线视频社区8| 日韩精品免费视频一区二区三区| 亚洲专区国产一区二区| 国产伦一二天堂av在线观看| 中文字幕av电影在线播放| 亚洲在线自拍视频| 亚洲 欧美一区二区三区| 亚洲第一电影网av| 国产成人av教育| 亚洲国产欧美一区二区综合| 国产一区二区激情短视频| 亚洲精品美女久久av网站| 色综合亚洲欧美另类图片| 色在线成人网| 久久久久久人人人人人| 91老司机精品| 99久久国产精品久久久| cao死你这个sao货| 啦啦啦 在线观看视频| 欧美在线黄色| 日韩免费av在线播放| 一边摸一边做爽爽视频免费| 亚洲av中文字字幕乱码综合 | 欧美日韩瑟瑟在线播放| 高潮久久久久久久久久久不卡| 婷婷丁香在线五月| 欧美日韩福利视频一区二区| 午夜福利在线观看吧| 国产精品爽爽va在线观看网站 | 宅男免费午夜| 亚洲国产高清在线一区二区三 | 中文亚洲av片在线观看爽| 国产精品乱码一区二三区的特点| 久久久国产成人精品二区| 日本 欧美在线| 嫩草影院精品99| av有码第一页| 久久精品国产清高在天天线| 自线自在国产av| 欧美亚洲日本最大视频资源| 老汉色∧v一级毛片| www日本在线高清视频| netflix在线观看网站| 十分钟在线观看高清视频www| 欧美乱妇无乱码| 一级片免费观看大全| 性色av乱码一区二区三区2| 国产欧美日韩精品亚洲av| 一a级毛片在线观看| 麻豆久久精品国产亚洲av| 国产精品1区2区在线观看.| 超碰成人久久| 亚洲国产欧美网| 特大巨黑吊av在线直播 | 搞女人的毛片| 久久精品国产综合久久久| 国产成人av激情在线播放| 日韩大尺度精品在线看网址| 黄色丝袜av网址大全| 一卡2卡三卡四卡精品乱码亚洲| 两性夫妻黄色片| 在线观看免费日韩欧美大片| 国产av一区二区精品久久| 久久人人精品亚洲av| 国产精品影院久久| 国产一级毛片七仙女欲春2 | 亚洲精品国产精品久久久不卡| 中文字幕人成人乱码亚洲影| 人人妻,人人澡人人爽秒播| 999精品在线视频| 亚洲欧美精品综合一区二区三区| 亚洲午夜理论影院| 亚洲中文日韩欧美视频| 精品国内亚洲2022精品成人| 国产久久久一区二区三区| 日本在线视频免费播放| 亚洲欧美精品综合久久99| 美女扒开内裤让男人捅视频| 欧美日韩福利视频一区二区| 久久九九热精品免费| 亚洲五月婷婷丁香| 午夜福利免费观看在线| 哪里可以看免费的av片| 成在线人永久免费视频| 色综合站精品国产| 久久久精品欧美日韩精品| 免费女性裸体啪啪无遮挡网站| 亚洲黑人精品在线| 国产成人影院久久av| 在线观看免费午夜福利视频| 黄色片一级片一级黄色片| 国产日本99.免费观看| 国产极品粉嫩免费观看在线| 欧美人与性动交α欧美精品济南到| 国产麻豆成人av免费视频| 国产成人av教育| 欧美人与性动交α欧美精品济南到| 好男人在线观看高清免费视频 | 中文字幕精品免费在线观看视频| 亚洲最大成人中文| 国产精品一区二区精品视频观看| 亚洲国产高清在线一区二区三 | 精华霜和精华液先用哪个| 一区福利在线观看| 亚洲激情在线av| 亚洲国产欧洲综合997久久, | 亚洲av电影不卡..在线观看| 色综合欧美亚洲国产小说| 亚洲欧美精品综合久久99| 婷婷丁香在线五月| 美女扒开内裤让男人捅视频| 日本三级黄在线观看| 精品一区二区三区av网在线观看| 真人做人爱边吃奶动态| 国产亚洲精品综合一区在线观看 | svipshipincom国产片| 国产aⅴ精品一区二区三区波| 色尼玛亚洲综合影院| 久久久久精品国产欧美久久久| xxx96com| 亚洲欧美日韩高清在线视频| 午夜免费鲁丝| 精品久久久久久久人妻蜜臀av| 亚洲国产欧美网| 久久久久免费精品人妻一区二区 | 在线观看www视频免费| 免费在线观看影片大全网站| 亚洲成人久久性| 午夜福利在线在线| 亚洲自拍偷在线| 久久99热这里只有精品18| 国产一区二区激情短视频| 国产伦人伦偷精品视频| 视频在线观看一区二区三区| 黄片大片在线免费观看| 好男人在线观看高清免费视频 | 精品福利观看| 国产爱豆传媒在线观看 | av中文乱码字幕在线| 日韩欧美一区视频在线观看| 亚洲,欧美精品.| 国产精品免费一区二区三区在线| 女生性感内裤真人,穿戴方法视频| 俄罗斯特黄特色一大片| 男女之事视频高清在线观看| 一边摸一边抽搐一进一小说| 女警被强在线播放| 可以在线观看的亚洲视频| 欧美乱妇无乱码| 久久午夜综合久久蜜桃| 黄色a级毛片大全视频| 欧美 亚洲 国产 日韩一| 午夜日韩欧美国产| 好看av亚洲va欧美ⅴa在| 动漫黄色视频在线观看| 久久久久亚洲av毛片大全| 亚洲一区二区三区不卡视频| 国产精品一区二区免费欧美| 高潮久久久久久久久久久不卡| 美国免费a级毛片| 在线观看66精品国产| 久久久久久九九精品二区国产 | 亚洲精品中文字幕在线视频| 国产一区二区三区视频了| 天天躁夜夜躁狠狠躁躁| www日本黄色视频网| 精品乱码久久久久久99久播| 久久久水蜜桃国产精品网| 听说在线观看完整版免费高清| 亚洲七黄色美女视频| 男人的好看免费观看在线视频 | 亚洲五月天丁香| 桃红色精品国产亚洲av| 99精品久久久久人妻精品| 日本一区二区免费在线视频| 特大巨黑吊av在线直播 | 国产又爽黄色视频| 欧美午夜高清在线| 天堂动漫精品| 欧美日韩中文字幕国产精品一区二区三区| 人人妻人人澡欧美一区二区| 国产亚洲精品第一综合不卡| 女同久久另类99精品国产91| 午夜福利高清视频| 日韩精品免费视频一区二区三区| 成人亚洲精品一区在线观看| 人成视频在线观看免费观看| 午夜精品在线福利| 久9热在线精品视频| 99精品久久久久人妻精品| 成熟少妇高潮喷水视频| 最近最新中文字幕大全电影3 | 国产精品野战在线观看| 久久精品亚洲精品国产色婷小说| 久久婷婷成人综合色麻豆| 亚洲精品在线观看二区| 亚洲成人国产一区在线观看| 国产精品精品国产色婷婷| 香蕉国产在线看| 久久精品国产99精品国产亚洲性色| 中文在线观看免费www的网站 | 日韩欧美国产在线观看| 欧美乱码精品一区二区三区| 国产亚洲av高清不卡| 此物有八面人人有两片| 美女高潮喷水抽搐中文字幕| 亚洲一码二码三码区别大吗| 最近最新免费中文字幕在线| 亚洲人成电影免费在线| 中文字幕精品亚洲无线码一区 | 欧美一级a爱片免费观看看 | 真人做人爱边吃奶动态| 岛国在线观看网站| 老汉色∧v一级毛片| 一级a爱片免费观看的视频| 久久精品国产亚洲av香蕉五月| 嫩草影院精品99| 最近在线观看免费完整版| 久久久久九九精品影院| 国产单亲对白刺激| 黄片播放在线免费| 欧美乱码精品一区二区三区| 动漫黄色视频在线观看| 久久国产精品男人的天堂亚洲| 丰满的人妻完整版| www日本黄色视频网| 亚洲 国产 在线| 亚洲成人久久性| 国产视频一区二区在线看| 一二三四在线观看免费中文在| svipshipincom国产片| 日本在线视频免费播放| 久久伊人香网站| 在线观看日韩欧美| 免费看美女性在线毛片视频| or卡值多少钱| 在线永久观看黄色视频| 老汉色∧v一级毛片| 精品国产超薄肉色丝袜足j| 色哟哟哟哟哟哟| 男人舔女人的私密视频| 一进一出抽搐gif免费好疼| 18禁观看日本| 国产99白浆流出| 伊人久久大香线蕉亚洲五| 色哟哟哟哟哟哟| 国产亚洲av高清不卡| 国产爱豆传媒在线观看 | 真人做人爱边吃奶动态| 国产亚洲av高清不卡| 亚洲片人在线观看| 久久久精品欧美日韩精品| 青草久久国产| 精品久久久久久成人av| 国产麻豆成人av免费视频| 欧美中文日本在线观看视频| 精品久久久久久久人妻蜜臀av| 动漫黄色视频在线观看| 丰满的人妻完整版| 一个人观看的视频www高清免费观看 | 亚洲中文字幕日韩| 精品国产一区二区三区四区第35| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| 午夜久久久久精精品| 黄色成人免费大全| 亚洲第一av免费看| 天天添夜夜摸| 女警被强在线播放| 欧美最黄视频在线播放免费| 国产99白浆流出| 久久天堂一区二区三区四区| 欧美在线一区亚洲| 成年女人毛片免费观看观看9| 黑人操中国人逼视频| 天天躁夜夜躁狠狠躁躁| 久久久水蜜桃国产精品网| 香蕉国产在线看| 久久久国产成人精品二区| 人人妻人人看人人澡| 自线自在国产av| 日韩三级视频一区二区三区| 操出白浆在线播放| 黑人欧美特级aaaaaa片| 在线观看免费视频日本深夜| 国产精品香港三级国产av潘金莲| 国产成人精品无人区| 精品久久蜜臀av无| 白带黄色成豆腐渣| 一区福利在线观看| 最近最新免费中文字幕在线| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 成人av一区二区三区在线看| 国产午夜精品久久久久久| 视频在线观看一区二区三区| 午夜成年电影在线免费观看| 777久久人妻少妇嫩草av网站| 午夜a级毛片| 亚洲国产欧美一区二区综合| 国产激情久久老熟女| 女同久久另类99精品国产91| 1024香蕉在线观看| а√天堂www在线а√下载| 99精品久久久久人妻精品| 久久青草综合色| www.自偷自拍.com| 亚洲九九香蕉| 亚洲男人的天堂狠狠| 欧美日韩精品网址| 精品乱码久久久久久99久播| 久久久久九九精品影院| 亚洲在线自拍视频| 精品电影一区二区在线| 国产三级黄色录像| 国产精品美女特级片免费视频播放器 | 亚洲一区二区三区不卡视频| 女人高潮潮喷娇喘18禁视频| av在线天堂中文字幕| 少妇的丰满在线观看| 欧美绝顶高潮抽搐喷水| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 午夜成年电影在线免费观看| 老司机在亚洲福利影院| 大型黄色视频在线免费观看| 国产在线观看jvid| 免费av毛片视频| 国产精品免费一区二区三区在线| 国产欧美日韩精品亚洲av| 日韩欧美免费精品| www.www免费av| 欧美色欧美亚洲另类二区| 熟女少妇亚洲综合色aaa.| 中文字幕最新亚洲高清| 国产伦一二天堂av在线观看| 精品国产亚洲在线| 国产又爽黄色视频| 日韩三级视频一区二区三区| 亚洲中文av在线| 中国美女看黄片| 久久这里只有精品19| 精品久久久久久久久久久久久 | 亚洲成人久久性| 亚洲精品中文字幕一二三四区| 欧美黄色淫秽网站| 精品第一国产精品| 人人妻人人澡欧美一区二区| 狠狠狠狠99中文字幕| 中文字幕久久专区| 国产精品香港三级国产av潘金莲| 日韩有码中文字幕| 亚洲五月色婷婷综合| 国产亚洲精品久久久久5区| xxx96com| 久久国产精品男人的天堂亚洲| 别揉我奶头~嗯~啊~动态视频| 欧美日韩精品网址| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 亚洲国产精品合色在线| 88av欧美| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 亚洲第一av免费看| 高清毛片免费观看视频网站| 成人一区二区视频在线观看| av中文乱码字幕在线| 中亚洲国语对白在线视频| 高清在线国产一区| 丝袜在线中文字幕| 日韩欧美免费精品| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线| 啪啪无遮挡十八禁网站| 亚洲欧美日韩高清在线视频| 国产男靠女视频免费网站| 欧美激情高清一区二区三区| 国产激情久久老熟女| 老汉色∧v一级毛片| 久久精品国产综合久久久| 午夜免费成人在线视频| 午夜亚洲福利在线播放| 欧美国产日韩亚洲一区| 欧美亚洲日本最大视频资源| 久久久久久久精品吃奶| 美国免费a级毛片| 婷婷六月久久综合丁香| 欧美日韩一级在线毛片| 色婷婷久久久亚洲欧美| 国产蜜桃级精品一区二区三区| svipshipincom国产片| 热re99久久国产66热| 熟妇人妻久久中文字幕3abv| 两性夫妻黄色片| 男女那种视频在线观看| 精品人妻1区二区| 国产成人影院久久av| 熟妇人妻久久中文字幕3abv| 国产91精品成人一区二区三区| 中文字幕精品亚洲无线码一区 | 每晚都被弄得嗷嗷叫到高潮| 免费看十八禁软件| 女性被躁到高潮视频| 欧美激情极品国产一区二区三区| 午夜久久久久精精品| 日韩成人在线观看一区二区三区| 麻豆av在线久日| 一本大道久久a久久精品| 国产一区二区在线av高清观看| 伦理电影免费视频| 亚洲电影在线观看av| 欧美成人性av电影在线观看| 可以在线观看毛片的网站| 国产又爽黄色视频| 久久香蕉激情| 日韩精品中文字幕看吧| 琪琪午夜伦伦电影理论片6080| 国产精品 欧美亚洲| 男人舔奶头视频| 欧美成狂野欧美在线观看| 国产激情偷乱视频一区二区| 国产精品久久久av美女十八| 国产精华一区二区三区| 无限看片的www在线观看| 国产91精品成人一区二区三区| e午夜精品久久久久久久| 看黄色毛片网站| 国产爱豆传媒在线观看 | 最好的美女福利视频网| 久久精品91无色码中文字幕| 精品高清国产在线一区| 国产av一区二区精品久久| 日韩一卡2卡3卡4卡2021年| 欧美不卡视频在线免费观看 | 久久精品aⅴ一区二区三区四区| 18禁黄网站禁片午夜丰满| 久久国产精品影院| 日韩欧美免费精品| 成人一区二区视频在线观看| 亚洲第一av免费看| 欧美色视频一区免费| 超碰成人久久| av有码第一页| xxx96com| 久久亚洲精品不卡| 国产麻豆成人av免费视频| 色老头精品视频在线观看| 国产精品精品国产色婷婷| 国产精品永久免费网站| 老司机靠b影院| 亚洲aⅴ乱码一区二区在线播放 | 级片在线观看| 国产在线观看jvid| 最近在线观看免费完整版| 国产精品一区二区精品视频观看| 嫩草影院精品99| 欧美一级a爱片免费观看看 | 搡老岳熟女国产| bbb黄色大片| 婷婷亚洲欧美| 听说在线观看完整版免费高清| 亚洲国产精品久久男人天堂| 此物有八面人人有两片| 香蕉国产在线看| 曰老女人黄片| av在线天堂中文字幕| 可以在线观看毛片的网站| 国产三级在线视频| 99精品久久久久人妻精品| 精品少妇一区二区三区视频日本电影| 亚洲九九香蕉| 中文字幕av电影在线播放| 亚洲成人国产一区在线观看| 久久久国产成人精品二区| 好男人电影高清在线观看| 麻豆av在线久日| 在线十欧美十亚洲十日本专区| 视频区欧美日本亚洲| 亚洲人成网站高清观看| 一进一出好大好爽视频| 亚洲人成77777在线视频| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av| 亚洲成人国产一区在线观看| 欧美成人免费av一区二区三区| 日本五十路高清| 国产精华一区二区三区| 白带黄色成豆腐渣| 久久中文字幕一级| 高潮久久久久久久久久久不卡| 国产高清激情床上av| av片东京热男人的天堂| 好男人电影高清在线观看| 国产色视频综合| 午夜两性在线视频| 日本 欧美在线| 国产精品国产高清国产av| 51午夜福利影视在线观看| 色播在线永久视频| 午夜免费观看网址| 人妻久久中文字幕网| 色在线成人网| aaaaa片日本免费| 久久精品人妻少妇| 久久香蕉激情| 女同久久另类99精品国产91| 国产精品久久久久久人妻精品电影| 一本大道久久a久久精品| 一个人免费在线观看的高清视频| 亚洲 欧美一区二区三区| 亚洲,欧美精品.| 高清在线国产一区| 久热这里只有精品99| 亚洲av电影不卡..在线观看| 成人三级黄色视频| 国产视频内射| 12—13女人毛片做爰片一| 欧美成狂野欧美在线观看| 一边摸一边做爽爽视频免费| 男女视频在线观看网站免费 | 久久伊人香网站| 精品国产一区二区三区四区第35| 十分钟在线观看高清视频www| 免费看a级黄色片| 欧美色欧美亚洲另类二区| 此物有八面人人有两片| 在线观看免费午夜福利视频| 亚洲精品国产区一区二| 男人的好看免费观看在线视频 | 日本熟妇午夜| 51午夜福利影视在线观看| 不卡av一区二区三区| 免费观看人在逋| 国产成+人综合+亚洲专区| 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 国产亚洲欧美在线一区二区| 欧美最黄视频在线播放免费| 波多野结衣高清作品| 午夜免费鲁丝| 十八禁网站免费在线| 亚洲 欧美一区二区三区| 啦啦啦 在线观看视频| 中文字幕久久专区| 国产欧美日韩精品亚洲av| 人妻丰满熟妇av一区二区三区| 亚洲av成人av| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 欧美性猛交黑人性爽| 欧美日本视频| 一个人免费在线观看的高清视频| 色播在线永久视频| 天天添夜夜摸| 好看av亚洲va欧美ⅴa在| av在线天堂中文字幕| 久久精品亚洲精品国产色婷小说| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 国产精品av久久久久免费| 欧美 亚洲 国产 日韩一| 两个人视频免费观看高清| 欧美不卡视频在线免费观看 | 性色av乱码一区二区三区2| 好看av亚洲va欧美ⅴa在| 精品免费久久久久久久清纯| 超碰成人久久| 日日干狠狠操夜夜爽| 日韩精品中文字幕看吧| 久久精品91无色码中文字幕|