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

    深海熱液噴口金屬硫化物中甲烷菌的多樣性研究

    2012-10-23 03:01:32王淑芳錢媛媛
    海洋科學(xué) 2012年6期
    關(guān)鍵詞:古菌胡安文庫

    王淑芳, 陳 月 錢媛媛

    (1. 淮海工學(xué)院 海洋學(xué)院, 江蘇 連云港 222005; 2. 淮海工學(xué)院 江蘇省海洋生物技術(shù)重點建設(shè)實驗室, 江蘇 連云港 222005)

    深海熱液噴口金屬硫化物中甲烷菌的多樣性研究

    王淑芳1,2, 陳 月1, 錢媛媛1

    (1. 淮海工學(xué)院 海洋學(xué)院, 江蘇 連云港 222005; 2. 淮海工學(xué)院 江蘇省海洋生物技術(shù)重點建設(shè)實驗室, 江蘇 連云港 222005)

    對兩個來自太平洋洋脊胡安·德富卡的兩個熱液煙囪的金屬硫化物4136-2和4148-B1進行甲烷菌甲基輔酶M還原酶的編碼基因mcrA的序列擴增, 構(gòu)建克隆文庫并進行分子進化分析。結(jié)果表明在這個甲烷富集的熱液噴口周圍含有豐富的甲烷產(chǎn)生菌, 沒有任何甲烷氧化菌存在。兩個硫化物樣品的甲烷產(chǎn)生菌種類完全不同。在4136-2硫化物中的甲烷產(chǎn)生菌都與熱液口的高溫環(huán)境有關(guān)系, 主要屬于甲烷球菌目的甲烷暖球菌(Methanocaldococcus), 少部分屬于甲烷火菌目甲烷嗜高熱菌屬的坎氏甲烷嗜高熱菌(Methanopyrus kandleri)。與這兩個屬的可分離菌株的氨基酸同源性為89%~97%, 核苷酸同源性高達92%~100%。4148-B1硫化物中發(fā)現(xiàn)的一類疑似甲基輔酶M還原酶的編碼序列, 它們與已知的甲烷菌mcrA序列核苷酸同源性為69%~72%, 氨基酸同源性僅為43%~47%。這可能是由于4148-B1來自于正在噴發(fā)的超高溫?zé)嵋簢娍谙嚓P(guān)。由于與已發(fā)表的甲烷菌克隆子或菌株同源性較低, 有可能是熱液口特有的以前未發(fā)現(xiàn)的甲烷菌。

    胡安·德富卡; 金屬硫化物; 甲烷產(chǎn)生菌; mcrA

    20世紀(jì)70年代, 科學(xué)家在東太平洋深海區(qū)發(fā)現(xiàn)的海底熱液口及其周圍生機盎然的生物群。這是生命科學(xué)史上一次具有深遠意義的重大發(fā)現(xiàn), 也是人類認識自然的一大進步。研究熱液口環(huán)境下的微生物群落結(jié)構(gòu)對探索生命起源具有重要的基礎(chǔ)研究意義。

    在太平洋洋脊胡安·德富卡(Juan de Fuca Ridge,JDF)熱液場 Endeavour部分, 甲烷濃度高達 1.8~3.4 mmol/L, 二氧化碳濃度高達11.6~18.2 mmol/L, 氫氣濃度高達 50 mmol/L[1]。遠遠高于 EPR、SJdF和Guanymas的熱液口。在東太平洋樣脊一次爆發(fā)過程中,“A”煙囪口中的氫氣從 25 mmol/kg升高到 45 mmol/kg, 甲烷的濃度從 0.13 mmol/kg升高到 0.19 mmol/kg[2]。Gorda熱液噴發(fā)后的甲烷濃度達到背景海水的16倍, 濃度從2 nmol/L到100 nmol/L, 與其他熱液場穩(wěn)定的熱液口環(huán)境中的濃度相似[3]。熱液流體一般含有2~20 mmol/kg的二氧化碳。另外在非常熱的流體中還檢測到約為 300 mmol/kg的二氧化碳[4,5]。豐富的甲烷含量為甲烷氧化菌的生命活動提供必要的物質(zhì)能量, 豐富的氫和二氧化碳為甲烷產(chǎn)生菌的生命活動提供必要的物質(zhì)和能量。生物可以利用這些氣體作為能源維持生長和代謝活動[6]。

    Endeavour 熱液系統(tǒng)最重要的特點是熱液流體具有高濃度的甲烷和氨根離子[7]。熱液流體中甲烷的碳同位素組成最高而且出現(xiàn)在海底系統(tǒng)中沒有沉積物覆蓋的位置, 同時這些同位素數(shù)據(jù)與高達背景值900倍的甲烷含量相呼應(yīng), 說明Endeavour熱液場海底有豐富的甲烷產(chǎn)生菌[8]。能夠利用二氧化碳產(chǎn)生甲烷的超嗜熱的甲烷菌已經(jīng)從東太平洋脊(EPR)13°N,21°N[9-11]和胡安·德富卡熱液場的 Main Endeavour附近的熱液煙囪中分離得到[12]。本文采用甲烷菌的甲基輔酶M還原酶α亞單位編碼基因(mcrA)作為分子標(biāo)記檢測甲烷產(chǎn)生古菌和甲烷氧化古菌在胡安·德富卡熱液煙囪的金屬硫化物4148-B1及4136-2中甲烷循環(huán)中發(fā)揮的作用。

    1 材料與方法

    1.1 菌株、培養(yǎng)基與試劑

    大腸桿菌 (Escherichia coli) DH5α 為本實驗室保存; pMD18-T Vector購自寶生物工程(大連)有限公司(TaKaRa Biotechnology (Dalian) Co., Ltd.); 氨芐青霉素 (Amp)購自南京創(chuàng)瑞公司; PCR引物由上海英俊生物技術(shù)有限公司合成; DNA回收試劑盒購自上海生工公司的; 培養(yǎng)基配制參照 Maniatis等[13]的分子克隆方法。

    1.2 采樣點特征描述及金屬硫化物巖石樣品的采集

    本試驗所用樣品是2005年秋季, 由中國和美國兩國科學(xué)家在胡安·德富卡洋脊(Juan de Fuca Ridge或JDF)熱液區(qū)由Atlantis/Alvin科考船采集的熱液煙囪樣品。采樣噴口位點如圖1所示。樣品4148是在Main Endeavour熱液場(MEF)采集的一個正在噴發(fā)的煙囪體, 可以清晰地觀察到熱液通道的內(nèi)中外各層, 以及各層礦物不同的顏色, 中心溫度估計超過350℃。煙囪結(jié)構(gòu)在超凈工作臺把煙囪體從煙囪通道到外壁分成三層, 4148-B1是煙囪的最外層, 估計溫度在 100℃。硫化物樣品 4136-2采自 MEF北面的Clam Bed熱液場。Clam Bed熱液場是一個大的熱液彌散流活動場。這里觀察到有微生物席、螃蟹和管狀蠕蟲的存在, 但沒有發(fā)現(xiàn)以前這里的標(biāo)志性生物貽貝(Clam), 探測器顯示此處最高溫度是29.2℃。采集的硫化物樣品迅速冷卻, 放入-20℃冰箱保存直到進行DNA提取。在保存過程中沒有發(fā)現(xiàn)樣品有外觀的變化。

    1.3 金屬硫化物樣品中總DNA的提取

    硫化物巖石樣品用已滅菌的鑿子和小刀從大的樣品上分離下來, 用研缽研磨。DNA的提取方法參照zhou等[14]的方法。30 mL的離心管加入5g樣品和13.5 mL DNA提取液(100 mmol/L Tris-HCl [pH 8.0], 100 mmol/L EDTA鈉鹽[pH8.0], 100 mmol/L磷酸鈉 [pH 8.0], 1.5 mol/L NaCl, 1% CTAB)和 100μL蛋白酶K(10 g/L)在37℃溫浴30 min。再加入1.5 mL 20% SDS, 在65℃溫浴2 h, 溫浴過程中每隔15~20 min混勻一次。溫浴結(jié)束離心并收集上清, 加入氯仿/異戊醇(24:1)進行抽提兩次。收集上清, 加入 2/3體積的異丙醇進行 DNA沉淀, 離心得到的沉淀就是DNA的粗提物。用1%的瓊脂糖檢測提取的總DNA的大小。粗提物用上海生工的UNIQ-10柱式PCR產(chǎn)物回收試劑盒純化。純化后的DNA放入-80℃冰箱保存。

    1.4 mcrA基因PCR擴增

    mcrA擴增參照前人的方法[15]用通用引物 ME1(5′-GCM ATG CAR ATH GGW ATG TC-3′)和 ME2(5′-TCA TKG CRT AGT TDG GRT AGT-3′)擴增環(huán)境基因組的mcr基因。反映程序如下: 94℃變性 3 min;再接94℃40 s, 49℃1 min, 72℃1 min, 35個循環(huán); 最后72℃延伸10 min。擴增反應(yīng)體系:DNA (約20~50 ng/μL) 1 μL, Taq 聚合酶 Buffer(10×)5μL,dNTP(20 mmol/L) 4 μL,引物(10 pmol/uL)各 1 μL, Mg2+(25 mmol/L)3~5 μL, Taq 酶(1U/μL) 3μL, 補水至 50 μL。擴增2~3管重復(fù), 產(chǎn)物用1%的瓊脂糖凝膠電泳檢測并回收。

    1.5 mcrA基因的文庫構(gòu)建及陽性菌篩選

    對PCR產(chǎn)物進行膠回收以純化mcrA基因片段,然后用pMD18-T克隆試劑盒(TaKaRa)按照操作說明連接到pMD18-T載體并轉(zhuǎn)化到DH5α大腸桿菌感受態(tài)細胞,從而構(gòu)建 mcrA基因的克隆文庫。37℃培養(yǎng)12~16 h后, 用 pMD18-T載體的通用引物 UPI(5′-GGAAACAGCTATGACCATG-3′)和 UPII(5′-TTGGGTAACGCCAGGGT-3′)進行菌落 PCR, 檢測插入片段大小。檢測到插入片段在960 bp左右的克隆子為陽性。隨機挑選不同數(shù)目的克隆子進行測序。測序數(shù)目根據(jù)多樣性指數(shù)的計算結(jié)果進行, 以文庫的庫容(Coverage)[16]達到85%以上為標(biāo)準(zhǔn)。

    1.6 mcrA陽性克隆子測序及序列分析

    引物M13+/-的PCR產(chǎn)物用Applied Biosystems 3730xl 測序儀進行序列核苷酸測序。核苷酸序列用DNAMAN 6.0蛋白翻譯功能翻譯成氨基酸序列。氨基酸序列和核苷酸序列在 NCBI上的網(wǎng)上數(shù)據(jù)庫進行同源序列比對, 查找同源性最高的克隆子和菌株并對所有氨基酸序列用ClustalX 1.83進行比對。最終截取78個序列的173個氨基酸用Mega 3.1(http://www.megasoftware.net/index.html)進行分子進化樹的構(gòu)建, 并以1000次的步長和臨位法進行計算。Dotur[17](http://www.plantpath.wisc.edu/fac/joh/dotur.html)用來做稀釋曲線和Chao1非參數(shù)豐度度評估及西普森指數(shù)(Simpson)和香農(nóng)指數(shù)(Shanoon index)。庫容(Coverage)用Good[16]的計算公式進行手工計算。最小分類操作單元(OTUs)用來定義不同的克隆子群,本文以古菌mcrA氨基酸大于等于 98%作為分類最小分類操作單元的標(biāo)準(zhǔn)。

    圖1 Endeavour段Main Endeavour熱液場采樣噴口和樣品圖Fig. 1 Hydrothermal field vent and sample of the endeavour segment main endeavour

    2 結(jié)果

    口的兩個金屬硫化物巖石中提取的環(huán)境總 DNA進行PCR擴增, 得到約760 bp的片段。克隆到DH5α菌株后, 4136-2隨機挑選34個陽性克隆子, 4148-B1隨機挑選27個陽性克隆子, 將獲得的61個克隆子進行測序。

    2.1 甲烷菌甲基輔酶M編碼基因mcrA的文庫構(gòu)建

    采用mcrA的通用引物ME1和ME2對來源于熱液

    2.2 甲烷菌甲基輔酶 M 編碼基因mcrA序列分析及多樣性評估

    獲得的 61個氨基酸序列在 NCBI(http://www.ncbi.nlm.nih.gov/)中進行比對, 搜索最同源的克隆子和可培養(yǎng)菌株, 結(jié)果如表1。61個序列中的43個氨基酸序列與甲烷古菌甲烷輔酶 M的α亞基同源, 同源性在 43%~97%之間。同時這 43個序列的核苷酸序列與已發(fā)表的甲烷甲基輔酶M的序列的同源性達到69%~100%。另外18個序列最同序列不是甲烷古菌的mcrA蛋白, 說明在擴增過程中出現(xiàn)了比較嚴重的非特異性擴增。在胡安·德富卡熱液口的金屬硫化物中甲烷菌都屬于廣古菌中的甲烷產(chǎn)生菌。4136-2文庫中 97.06%的甲烷菌來自在甲烷球菌(Methanococcales), 2.94%甲烷菌來自甲烷嗜高熱菌(Methanopyrales)。4148-B1與已知的克隆子及可培養(yǎng)菌株氨基酸同源性只有 43%~47%, 核苷酸同源性為69%~72%, 屬于未知的甲烷古菌。

    表1 胡安·德富卡熱液口甲烷菌mcrA序列NCBI數(shù)據(jù)庫比對分析Tab. 1 Analysis of mcrA sequences from JDF vent

    本試驗采用自動分析軟件 Dotur對兩個硫化物巖石中的的甲烷古菌進行多樣性評估, 結(jié)果如表2。Shannon-Wiener指數(shù), Simpson's指數(shù)以及Chao1[18-19]無偏估計, 文庫覆蓋率等幾個指數(shù)用來評估構(gòu)建的文庫的多樣性, 以評估文庫對取樣位點甲烷古菌多樣性的反應(yīng)程度。43個mcrA序列以≥98%的氨基酸同源性分為7個OUTs, 兩個OTUs來自4148-B1, 5個OTUs來自4136-2。4136-2和4148-B1mcrA基因構(gòu)建的克隆文庫覆蓋率分別為 94.1%和 88.9%,Shannon-Wiener指數(shù)分別為0.68和0.78, Chao1無偏估計值為5.33和2.00。這表明在胡安·德富卡熱液口甲烷菌的多樣性并不豐富, 試驗中構(gòu)建的兩個文庫得到的 7個 OTUs完全可以代表這兩個位點的甲烷古菌的多樣性分布。比較而言, 4148-B1文庫的H’及1/D較4136-2大, 表明4148-B1文庫的多樣性較4136-2豐富。稀釋曲線(Rarefaction curves)分析(圖2)支持以上結(jié)果。

    圖 2 熱液口金屬硫化物巖石中甲烷古菌 mcrA基因多樣性的稀釋曲線Fig. 2 Rarefaction curves of mcrA gene from vent sulfide rock

    表2 胡安·德富卡熱液口甲烷菌mcrA基因克隆文庫的多樣性評估Tab. 2 Measurement of mcrA clone library from JDF vent

    2.3 甲烷菌甲基輔酶 M 編碼基因mcrA系統(tǒng)進化分析

    圖3 胡安·德富卡熱液口甲烷菌mcrA基因編碼的氨基酸的系統(tǒng)發(fā)育樹Fig. 3 Phylogentic tree of mcrA gene from JDF vent系統(tǒng)發(fā)育樹采用鄰接法構(gòu)建; 粗體為本研究中獲得的甲烷菌mcrA基因的克隆子, 命名為4136-2和4148-B1分別表示克隆子來自于硫化物樣品4136-2和4148-B1; 分支上的數(shù)值代表經(jīng)1000次計算后的置信度值; 比例標(biāo)尺代表氨基酸5%的置換概率The tree was constructed by Neighbor-Joining Method; the bold ones were got in this study, names 4136-2 and 4148-B1 expressing clones from sulfide rock 4136-2 and 4148-B1, respectively; the values on the branches expressing 1000 bootstrap and the scale expressing 5% exchange rate

    NCBI比對結(jié)果為甲烷古菌mcrA的43個序列,用MAGE 3.1構(gòu)建分子進化樹, 如圖3所示。4136-2的mcrA序列與已分離的嗜熱菌有最大同源關(guān)系, 在89%~97%。與克隆子的同源性較菌株低, 在70%~72%。4136-2硫化物中的甲烷菌群落97.06%來自在甲烷球菌(Methanococcales), 2.94%來自甲烷嗜高熱菌(Methanopyrales)。4148-B1序列沒有找到氨基酸同源性超過47%的mcrA序列。核苷酸比對可以找到同源性在 69%~72%的同源序列, 說明這些克隆子序列屬于甲烷菌的甲基輔酶M的編碼基因。多樣性類型4136-2-8與Methanopyrus kandleri[20](GenBank序列號: AE010358)有最大同源關(guān)系, 同源性為95%。多樣性類型 4136-2-14與Methanococcus voltae[21](GenBank序列號: X07793)有最近同源關(guān)系,同源性為 89%。與這兩個序列同源關(guān)系最近的克隆子是來源于酸性泥炭的克隆子 FinME_09[22](GenBank序列號: DQ099689), 同源性分別為71%和70%。4136-2-1與 Methanocaldococcus jannaschii[23](GenBank序列號: F64405)的同源性為 97%。4136-2-20和4136-2-24與Methanocaldococcus jannaschii DSM 2661[23](GenBank 序列號: L77117)有最大同源關(guān)系, 同源性分別為93%和97%。與這三個序列最同源的克隆子 novmcr4[24](GenBank序列號:AF525516)來自富營養(yǎng)湖, 與三個序列的同源性在71%~72%間。樣品4148-B1的mcrA序列4148-B1-3與來源于富營養(yǎng)湖的克隆子novmcr2[24](GenBank序列號: AF525515)同源關(guān)系最近, 為 47%。與之同源關(guān)系最近的菌株是嗜熱菌 Methanopyrus kandleri AV19[25](GenBank序列號: AM114193)氨基酸同源性為 43%, 核苷酸同源性為 72%。另一類型克隆子4148-B1-59與來源于富含甲烷的水稻土壤的克隆子RC-I[25](GenBank 序 列 號 : NC_009464)及Methanosphaera stadtmanae DSM 3091[26](GenBank序列號: NC_007681)的氨基酸同源性最高, 都是43%。核苷酸比對時, 與RC-I的同源性為69%。

    3 分析與討論

    3.1 熱液口金屬硫化物甲烷菌類型分析

    甲烷菌甲基輔酶M編碼基因mcrA序列的系統(tǒng)發(fā)育樹分析已經(jīng)表明在胡安·德富卡熱液場熱液口硫化物巖石中甲烷古菌的類型比較豐富, 存在不同類型的甲烷產(chǎn)生菌及一類未知的甲烷菌。溫度較低的硫化物4136-2中的甲烷產(chǎn)生菌與已分離培養(yǎng)的甲烷產(chǎn)生菌有最近的同源關(guān)系。說明胡安·德富卡熱液口與其他熱液口一樣, 其附近環(huán)境中甲烷菌是熱液口附近比較常見的菌種。但熱液通道外壁的硫化物 4148-B1距離噴口較近, 其中的甲烷菌與已知的菌完全不同。

    甲烷產(chǎn)生菌是一類絕對厭氧菌, 是產(chǎn)生甲烷的生物來源[27]。甲烷球菌(Methanococcales)生活在中溫、嗜熱或極端嗜熱的環(huán)境中, 通過還原CO2和水生成甲烷產(chǎn)生能量供其活動[28]。胡安·德富卡熱液場周圍環(huán)境含有豐富的二氧化碳和氫氣, 為甲烷產(chǎn)生菌的代謝和活動提供物質(zhì)和能量[6,29-30]。本試驗中克隆子數(shù)最大的一個多樣性類型 4136-2-1(占總數(shù)的82.4%)的 97%的同源菌株 Methanocaldococcus jannaschii(GenBank序列號: F64405)最初分離自太平洋EPR熱液場的白煙囪附近的沉積物中, 發(fā)表在《Science》雜志上。另外一類多樣性類型屬于甲烷火菌(Methanopyrus)。與本試驗中得到的序列最近親緣關(guān)系的甲烷火菌是 Methanopyrus kandleri AV19(GenBank序列號: AE010358)。它是一株超嗜熱菌, 分離自加利福尼亞海灣二氧化碳和氫氣豐富的黑煙囪中, 與其他產(chǎn)甲烷生物一樣還原二氧化碳和水生成甲烷。

    3.2 金屬硫化物巖石中甲烷菌類型與環(huán)境的關(guān)系

    4136-2和 4148-B1兩個金屬硫化物樣品的甲烷菌甲基輔酶M編碼基因mcrA反映的甲烷菌類型反映了胡安·德富卡熱液口的甲烷菌的類型和變化。4136-2樣品中得到的甲烷菌全部是甲烷產(chǎn)生古菌,且與高溫來源的甲烷球菌及甲烷火菌有最近的同源關(guān)系。這反映了4136-2樣品采集自熱液噴口附近, 并且噴口附近含有豐富的二氧化碳和氫氣, 這與胡安·德富卡的文獻資料報道的相同。4148-B1樣品的克隆子雖然與已知的克隆子序列同源性不高, 但也與熱液口或水稻田等甲烷豐富的環(huán)境相關(guān)。其中一類與嗜熱菌的核苷酸同源性達到 72%, 氨基酸同源性為47%。由于序列很新, 更多的研究需要繼續(xù)。

    致謝: 本試驗是在國家海洋局第三海洋研究所海洋生物遺傳資源重點實驗室完成。實驗期間受到實驗室各位老師的指導(dǎo), 在此表示感謝。

    [1] Childress C, Fisher R. The biology of hydrothermal vent animals: physiology, biochemistry, and autotrophic symbioses, Oceanogr[J]. Mar Biol Annu Rev, 1992,30:337-441.

    [2] Lilley M D, Olson E J. Methane and hydrogen in active hydrothermal systems[J]. Annu Gold Conf, 2001, 11:3682.

    [3] Lilley M D. Real-time mapping of hydrothermal plumes on the endeavour segment of the Juan de Fuca[J]. Eos Trans Am Geophys Union, 1995,76: 420.

    [4] Lilley M D, Lupton J E, Damm K L V. Volatiles in the 9±N hydrothermal system: a comparison of 1991 and 1992 data[J]. Eos Trans Am Geophys Union, Fall Meet Suppl,1992. 73(43): 524.

    [5] Butterfield D A, Massoth G J. Geochemistry of the North Cleft Segment vent fluids-temporal changes in chlorinity and their possible relation to recent volcanism[J]. J Geophys Res, 1994, 99: 4951-4968.

    [6] McCollom T M. Geochemical constraints on primary productivity in submarine hydrothermal vent plumes[J].Deep-Sea Res I, 2000, 47: 85-101.

    [7] Lilley M D. Anomalous CH4and NH4+concentrations at an unsedimented mid-ocean-ridge hydrothermal system[J].Nature, 1993, 364: 45-47.

    [8] Lilley M D. Anomalous CH4and NH4+concentrations at an unsedimented mid-ocean-ridge hydrothermal system[J].Nature, 1993, 346: 45-47.

    [9] Jones W J, Paynter M J B, Gupta R. Characterization of Methanococcus maripaludis sp. nov., a new methanogen isolated from salt marsh sediment[J]. Arch Microbiol,1983,153: 91-97.

    [10] Jeanthon C. Methanococcus vulcanius sp. nov., a novel hyperthermophilic methanogen isolated from East Pacific Rise, and identification of Methanococcus sp. DSM 4213T as Methanococcus fervens sp. nov[J]. Int J Syst Bacteriol,1999, 49: 583-589.

    [11] Jeanthon C. Methanococcus infernus sp. nov., a novel hyperthermophilic lithotrophic methanogen isolated from a deep-sea hydrothermal vent[J]. Int J Syst Bacteriol, 1998,48: 913-919.

    [12] Summit M. Ecology, physiology, and phylogeny of subseafloor thermophiles from mid-ocean ridge environments[M]. Seattle :University of Washington, 2000.

    [13] Maniatis T, Fritsch E F, Sambrook J. Molecular Cloning:A Laboratory Manual[M]. New York: Cold Spring Harbour Laboratory Press, 1982: 440-468.

    [14] Zhou J, Bruns M A, Tiedje J M. DNA recovery from soils of diverse composition[J]. Appl Environ Microbiol, 1996,62(2): 316-322.

    [15] Hales B A. Isolation and identification of methanogen-specific DNA from blanket bog peat by PCR amplification and sequence analysis[J]. Appl Environ Microbiol,1996,62(2): 668-675.

    [16] Good I L. The population frequencies of species and the estimation of population parameters[J]. Biometrika, 1953,40: 237-264.

    [17] Schloss P D, Handelsman J. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness[J]. Appl Environ Microbiol, 2005, 71(3): 1501-1506.

    [18] Chao A. Non-parametric estimation of the number of classes in a population[J]. Scand J Stat, 1984, 11:265-270.

    [19] Chao A, Lee S M. Estimationg the number of classes via sample coverage[J]. J Am Stat Assoc, 1992, 87: 210-217.

    [20] Slesarev A I. The complete genome of hyperthermophile Methanopyrus kandleri AV19 and monophyly of archaeal methanogens[J]. Proc Natl Acad Sci U S A, 2002, 99(7):4644-4649.

    [21] Klein A. Comparative analysis of genes encoding methyl coenzyme M reductase in methanogenic bacteria[J]. Mol Gen Genet, 1988, 213(2-3): 409-420.

    [22] Metje M, Frenzel P. Effect of temperature on anaerobic ethanol oxidation and methanogenesis in acidic peat from a northern wetland[J]. Appl Environ Microbiol, 2005,71(12): p. 8191-8200.

    [23] Bult C J. Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii[J]. Science, 1996,273(5278): 1058-1073.

    [24] Earl J. Analysis of methanogen diversity in a hypereutrophic lake using PCR-RFLP analysis of mcr sequences[J].Microb Ecol, 2003, 46(2): 270-278.

    [25] Erkel C. Genome of rice cluster I archaea——the key methane producers in the rice rhizosphere[J]. Science,2006, 313(5785): 370-372.

    [26] Fricke W F. The Genome sequence of methanosphaera stadtmanae reveals why this human intestinal archaeon is restricted to methanol and H2for methane formation and ATP synthesis[J]. J Bacteriol, 2006, 188(2): 642-658.

    [27] Singh N. Isolation and characterization of methylotrophic methanogens from anoxic marine sediments in Skan Bay,Alaska: description of Methanococcoides alaskense sp.nov., and emended description of Methanosarcina baltica[J]. Int J Syst Evol Microbiol, 2005,55: 2531-2538.

    [28] Euzeby J P, Tindall B J. Nomenclatural type of orders:corrections necessary according to rules 15 and 21a of the bacteriological code (1990 Revision), and designation of appropriate nomenclatural types of classes and subclasses[J]. Int J Syst Evol Microbiol, 2001, 51: 725-727.

    [29] McCollom T M, Shock E L. Geochemical constraints on chemolithoautotrophic metabolism by microorganisms in seafloor hydrothermal systems[J]. Geochim Cosmochim Acta, 1997,61(20): 4375-4391.

    [30] Delaney J R. The quantum event of oceanic crustal accretion: Impacts of diking at mid-ocean ridges[J]. Science,1998, 281(5374): 222-230.

    Studies on the biodiversity of methanogens from hydrothermal vent chimney

    WANG Shu-fang1,2, CHEN Yue1, QIAN Yuan-yuan1
    (1. Huaihai Institute Technology, Oceanography Institute, Lianyungang 222005, China; 2. Key Laboratory of Marine Biotechnology in Jiangsu Province, Huaihai Institute of Technology, Lianyungang 222005, China)

    Jun.,25,2011

    juan de fuca (JDF); sulfide; methogen; mcrA

    The chimney sulfides (named 4136-2 and 4148-B1) obtained from Juan de Fuca Ridge were used as samples, aiming at mcrA gene amplification, clone library construction and phylogenetic tree analysis. The results showed that there were only CH4production microbes but no CH4oxidation microbes in the samples and the species in the two sulfides were different. The methanogen from 4136-2 sulfide rock were linked to a high temperature environment. These methanogen belong to Methanocaldococcus and Methanopyrus kandleri. The similarities of protein sequences between 4136-2 clones and isolated Methanocaldococcus and Methanopyrus kandleri strains were from 89% to 97%, and The similarities of mucleotide sequences between 4136-2 clones and isolated Methanocaldococcus and Methanopyrus kandleri strains were from 92% to 100%. The mcrA sequences from 4148-B1 sulfide rock had low identities with sequences deposited in GenBank. The similarities of mucleotide sequences identities of 4148-B1 clones and published mcrA sequences were form 69% to 72%, and the similarities of protein sequences were form 43% to 37%. This also linked to its super high temperature environment and may suggest there were un-discovered methanogen living in the vent.

    Q93.3

    A

    1000-3096(2012)06-0031-08

    2011-06-25;

    2011-12-11

    淮海工學(xué)院引進人才基金項目(KQ09021)

    王淑芳(1976-), 女, 河北安國人, 講師, 博士。從事微生物資源研究。電話: 0518-85895427, E-mail: wangshufang2001@hotmail.com

    (本文編輯:梁德海)

    猜你喜歡
    古菌胡安文庫
    不同pH和氧氣條件下土壤古菌與海洋古菌的競爭適應(yīng)機制*
    變油為氣,“榨干”廢棄油田
    大自然探索(2022年5期)2022-07-11 03:10:33
    海洋古菌
    專家文庫
    優(yōu)秀傳統(tǒng)文化啟蒙文庫
    幽默大師(2020年10期)2020-11-10 09:07:22
    關(guān)于推薦《當(dāng)代詩壇百家文庫》入選詩家的啟事
    中華詩詞(2019年1期)2019-11-14 23:33:56
    專家文庫
    變成什么好
    東太平洋海隆深海熱液區(qū)沉積物古菌多樣性分析
    飛天超人——胡安!
    久久久a久久爽久久v久久| 亚洲欧美成人综合另类久久久 | 青春草视频在线免费观看| 欧美成人一区二区免费高清观看| 国产精品爽爽va在线观看网站| 国产成人精品久久久久久| 日韩,欧美,国产一区二区三区 | 26uuu在线亚洲综合色| 3wmmmm亚洲av在线观看| 国产精品不卡视频一区二区| 晚上一个人看的免费电影| 午夜免费激情av| 变态另类丝袜制服| 亚洲欧美日韩无卡精品| 免费av不卡在线播放| 国产色婷婷99| 建设人人有责人人尽责人人享有的 | 午夜精品在线福利| av在线天堂中文字幕| av免费在线看不卡| 久久99精品国语久久久| 亚洲最大成人中文| 欧美潮喷喷水| 国产精品久久电影中文字幕| 国产片特级美女逼逼视频| 免费搜索国产男女视频| 精品国产三级普通话版| 99久久精品一区二区三区| 免费搜索国产男女视频| 国产黄a三级三级三级人| 日韩精品青青久久久久久| 九九在线视频观看精品| 又爽又黄a免费视频| 成人av在线播放网站| 一个人观看的视频www高清免费观看| 成人漫画全彩无遮挡| 黑人高潮一二区| 少妇人妻一区二区三区视频| 人妻夜夜爽99麻豆av| 国产三级在线视频| 久久久欧美国产精品| 别揉我奶头 嗯啊视频| 成年免费大片在线观看| 成人无遮挡网站| 欧美性感艳星| 久久久久久久午夜电影| 久久这里只有精品中国| 亚洲第一区二区三区不卡| 国产精品1区2区在线观看.| 内射极品少妇av片p| 国产在线男女| 观看免费一级毛片| 国产老妇伦熟女老妇高清| 嫩草影院新地址| 日本五十路高清| 色综合色国产| a级一级毛片免费在线观看| 高清毛片免费看| 国产精品伦人一区二区| 亚洲欧美一区二区三区国产| 大香蕉97超碰在线| av国产久精品久网站免费入址| 日本爱情动作片www.在线观看| 18禁在线播放成人免费| av免费在线看不卡| eeuss影院久久| 日日啪夜夜撸| 精品熟女少妇av免费看| 伦精品一区二区三区| 一级爰片在线观看| 99久国产av精品| 水蜜桃什么品种好| 亚洲欧美中文字幕日韩二区| 精品一区二区三区视频在线| 国产精品久久久久久av不卡| 国产精品精品国产色婷婷| 国产精品一区二区三区四区久久| 午夜a级毛片| 成人三级黄色视频| 亚洲精品久久久久久婷婷小说 | 亚洲无线观看免费| 亚洲在线观看片| 在线免费观看不下载黄p国产| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久久av| 神马国产精品三级电影在线观看| 国产免费又黄又爽又色| 亚洲欧美一区二区三区国产| 国产人妻一区二区三区在| 三级毛片av免费| 能在线免费观看的黄片| 欧美色视频一区免费| 大香蕉97超碰在线| 午夜激情欧美在线| 免费观看的影片在线观看| 69人妻影院| 99九九线精品视频在线观看视频| 国产国拍精品亚洲av在线观看| 九九久久精品国产亚洲av麻豆| 免费看a级黄色片| 久久久久久九九精品二区国产| 亚洲第一区二区三区不卡| 国内精品美女久久久久久| 黄色日韩在线| 黄色配什么色好看| 看黄色毛片网站| 亚洲无线观看免费| 搞女人的毛片| 不卡视频在线观看欧美| 国产精品乱码一区二三区的特点| 中文字幕人妻熟人妻熟丝袜美| 老司机影院成人| 国产精品久久久久久av不卡| 超碰av人人做人人爽久久| 亚洲av成人av| 色综合亚洲欧美另类图片| 国产精品一及| 亚洲国产欧洲综合997久久,| 爱豆传媒免费全集在线观看| 久久人人爽人人爽人人片va| 日本猛色少妇xxxxx猛交久久| 综合色av麻豆| 国产成人a∨麻豆精品| 免费观看性生交大片5| 丝袜美腿在线中文| 亚洲aⅴ乱码一区二区在线播放| 午夜亚洲福利在线播放| 久久久精品欧美日韩精品| 久久久久久久午夜电影| 高清av免费在线| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩卡通动漫| 久久久精品94久久精品| 男人的好看免费观看在线视频| 乱系列少妇在线播放| 熟女人妻精品中文字幕| 久久久色成人| 久久久久网色| 国产美女午夜福利| 久久精品国产亚洲网站| 免费看av在线观看网站| 欧美高清性xxxxhd video| 级片在线观看| 欧美一级a爱片免费观看看| 亚洲三级黄色毛片| 国产精品电影一区二区三区| 国产激情偷乱视频一区二区| 国产精品人妻久久久久久| 尾随美女入室| 久久久久精品人妻al黑| 久久久久国产网址| 国产欧美亚洲国产| 亚洲成色77777| 日本免费在线观看一区| 日本黄大片高清| 日韩av在线免费看完整版不卡| 汤姆久久久久久久影院中文字幕| 国产 一区精品| 国产一区二区在线观看日韩| 精品久久久久久电影网| av片东京热男人的天堂| 欧美精品国产亚洲| 亚洲精品国产色婷婷电影| 久久久久人妻精品一区果冻| 精品少妇久久久久久888优播| 久久精品久久精品一区二区三区| 国精品久久久久久国模美| 久久97久久精品| 久久青草综合色| 久久午夜福利片| 美女国产高潮福利片在线看| 欧美3d第一页| 精品一品国产午夜福利视频| 国产一区二区激情短视频 | 搡女人真爽免费视频火全软件| 18禁裸乳无遮挡动漫免费视频| 色吧在线观看| 一边亲一边摸免费视频| 国产视频首页在线观看| 亚洲 欧美一区二区三区| 男男h啪啪无遮挡| 老司机亚洲免费影院| 成人二区视频| 国产亚洲精品第一综合不卡 | 国产极品粉嫩免费观看在线| 各种免费的搞黄视频| 美女xxoo啪啪120秒动态图| 久久青草综合色| 高清欧美精品videossex| 国产精品国产三级国产专区5o| 国产无遮挡羞羞视频在线观看| 国产深夜福利视频在线观看| 在线观看人妻少妇| 久久狼人影院| www.av在线官网国产| 一级毛片黄色毛片免费观看视频| 亚洲精品视频女| 日韩精品有码人妻一区| 色吧在线观看| 国产高清国产精品国产三级| 国精品久久久久久国模美| 婷婷色综合www| 如何舔出高潮| 精品国产一区二区三区四区第35| 亚洲人成网站在线观看播放| 久久女婷五月综合色啪小说| 久久久久久久久久久久大奶| 亚洲国产av新网站| 狂野欧美激情性xxxx在线观看| 美女福利国产在线| 热99国产精品久久久久久7| 欧美丝袜亚洲另类| 高清毛片免费看| 日韩在线高清观看一区二区三区| 成人毛片60女人毛片免费| 免费高清在线观看视频在线观看| 日韩大片免费观看网站| 久久ye,这里只有精品| av国产精品久久久久影院| 国产精品久久久久成人av| 五月天丁香电影| 蜜桃国产av成人99| 国产在线一区二区三区精| 亚洲欧美成人精品一区二区| 高清黄色对白视频在线免费看| 欧美亚洲 丝袜 人妻 在线| 22中文网久久字幕| 日本与韩国留学比较| 国产一区二区在线观看日韩| 一本—道久久a久久精品蜜桃钙片| 18禁裸乳无遮挡动漫免费视频| 中文精品一卡2卡3卡4更新| 亚洲欧洲精品一区二区精品久久久 | 日日啪夜夜爽| 天天影视国产精品| 国产成人精品在线电影| 制服诱惑二区| 有码 亚洲区| 国产免费一区二区三区四区乱码| 国产日韩欧美视频二区| 国产又爽黄色视频| 色网站视频免费| 日韩人妻精品一区2区三区| 香蕉国产在线看| 在线观看三级黄色| 最近中文字幕高清免费大全6| 欧美 亚洲 国产 日韩一| 一级黄片播放器| 少妇精品久久久久久久| 亚洲,欧美,日韩| 九九在线视频观看精品| 久久99精品国语久久久| 亚洲av男天堂| 丰满迷人的少妇在线观看| 性色avwww在线观看| 国产深夜福利视频在线观看| 亚洲精品美女久久av网站| 亚洲经典国产精华液单| 中文欧美无线码| 国产视频首页在线观看| 在线 av 中文字幕| 巨乳人妻的诱惑在线观看| 成年人免费黄色播放视频| 两性夫妻黄色片 | 国产精品女同一区二区软件| 欧美 亚洲 国产 日韩一| 亚洲国产欧美日韩在线播放| 午夜视频国产福利| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 国产成人a∨麻豆精品| 少妇熟女欧美另类| h视频一区二区三区| 国产在线视频一区二区| 最近中文字幕2019免费版| 精品人妻偷拍中文字幕| 精品一区二区三区四区五区乱码 | 亚洲av国产av综合av卡| 丝袜人妻中文字幕| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦中文免费视频观看日本| 韩国av在线不卡| 久久狼人影院| 纵有疾风起免费观看全集完整版| 日韩制服骚丝袜av| 黄色怎么调成土黄色| 97人妻天天添夜夜摸| 久久久久久久亚洲中文字幕| 人体艺术视频欧美日本| 精品久久蜜臀av无| 丝袜美足系列| 国产精品一二三区在线看| 精品国产一区二区三区四区第35| 老司机影院毛片| 国产精品偷伦视频观看了| 亚洲内射少妇av| av国产久精品久网站免费入址| 国产无遮挡羞羞视频在线观看| 日韩精品免费视频一区二区三区 | 蜜臀久久99精品久久宅男| 十分钟在线观看高清视频www| 久久99热这里只频精品6学生| 熟女人妻精品中文字幕| 国产亚洲av片在线观看秒播厂| 在线天堂中文资源库| 在线观看免费日韩欧美大片| 久久人人爽人人爽人人片va| 最后的刺客免费高清国语| 亚洲国产欧美日韩在线播放| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 日韩制服骚丝袜av| 亚洲精品久久午夜乱码| 少妇熟女欧美另类| 国产xxxxx性猛交| 亚洲成色77777| 亚洲熟女精品中文字幕| 国产视频首页在线观看| 精品一区二区三区视频在线| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 亚洲色图 男人天堂 中文字幕 | 亚洲国产色片| 日韩欧美一区视频在线观看| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 三级国产精品片| 一区二区日韩欧美中文字幕 | av免费在线看不卡| 欧美日韩国产mv在线观看视频| 在线亚洲精品国产二区图片欧美| 捣出白浆h1v1| 久久精品久久久久久噜噜老黄| 毛片一级片免费看久久久久| 飞空精品影院首页| 不卡视频在线观看欧美| 日本av免费视频播放| 久久热在线av| 亚洲成人一二三区av| 五月开心婷婷网| 国产熟女午夜一区二区三区| 最近中文字幕高清免费大全6| 超碰97精品在线观看| 伦精品一区二区三区| 男人操女人黄网站| 精品卡一卡二卡四卡免费| 精品人妻一区二区三区麻豆| 中文字幕人妻丝袜制服| 日韩三级伦理在线观看| 国产亚洲精品久久久com| 亚洲国产欧美日韩在线播放| 免费人成在线观看视频色| 日韩av不卡免费在线播放| 亚洲国产精品一区三区| 久久久久久久国产电影| 夫妻性生交免费视频一级片| 又大又黄又爽视频免费| 嫩草影院入口| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 精品久久国产蜜桃| 亚洲五月色婷婷综合| 尾随美女入室| 久久精品国产自在天天线| 99久久综合免费| 天天躁夜夜躁狠狠躁躁| 精品酒店卫生间| 人人妻人人澡人人看| 在线 av 中文字幕| 亚洲中文av在线| 成人18禁高潮啪啪吃奶动态图| 日本91视频免费播放| 日本黄色日本黄色录像| 免费播放大片免费观看视频在线观看| 日本欧美国产在线视频| 国产伦理片在线播放av一区| 国产麻豆69| 丰满少妇做爰视频| 国产一区有黄有色的免费视频| 成年美女黄网站色视频大全免费| 一级片'在线观看视频| 久久久久国产精品人妻一区二区| 国产成人一区二区在线| 日韩精品有码人妻一区| 亚洲婷婷狠狠爱综合网| 日本欧美视频一区| 91午夜精品亚洲一区二区三区| 男女高潮啪啪啪动态图| 国产精品久久久久久精品古装| 亚洲国产日韩一区二区| 亚洲综合色惰| 亚洲欧美中文字幕日韩二区| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 最近的中文字幕免费完整| 国产伦理片在线播放av一区| 日韩精品有码人妻一区| a 毛片基地| videosex国产| 国产亚洲最大av| 又黄又爽又刺激的免费视频.| 久久99精品国语久久久| 91在线精品国自产拍蜜月| 国产色爽女视频免费观看| 久久免费观看电影| 成人影院久久| 丝袜脚勾引网站| 久久久欧美国产精品| 亚洲精品久久久久久婷婷小说| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 中文字幕人妻丝袜制服| 国国产精品蜜臀av免费| 成人二区视频| 肉色欧美久久久久久久蜜桃| √禁漫天堂资源中文www| 另类精品久久| 国产男女超爽视频在线观看| 国产精品免费大片| 久久久久久久国产电影| 亚洲国产日韩一区二区| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 最近中文字幕2019免费版| 亚洲欧洲日产国产| 国产成人欧美| 亚洲中文av在线| 草草在线视频免费看| 免费观看无遮挡的男女| videos熟女内射| 中文天堂在线官网| 欧美性感艳星| 国产亚洲一区二区精品| 免费观看av网站的网址| 国产亚洲最大av| 男人舔女人的私密视频| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区三区| 丝袜在线中文字幕| 男女免费视频国产| 国产成人av激情在线播放| 乱码一卡2卡4卡精品| 亚洲精品视频女| 午夜福利影视在线免费观看| 欧美人与善性xxx| 各种免费的搞黄视频| 超碰97精品在线观看| 少妇高潮的动态图| 精品人妻偷拍中文字幕| 色哟哟·www| 18禁观看日本| 在线看a的网站| 日本爱情动作片www.在线观看| 国产又爽黄色视频| freevideosex欧美| 国产1区2区3区精品| 国产精品麻豆人妻色哟哟久久| 亚洲一区二区三区欧美精品| 国产1区2区3区精品| 在线观看一区二区三区激情| 看非洲黑人一级黄片| 欧美 亚洲 国产 日韩一| 99re6热这里在线精品视频| 一本色道久久久久久精品综合| 久热久热在线精品观看| 制服丝袜香蕉在线| 日本爱情动作片www.在线观看| 国产精品久久久av美女十八| 91aial.com中文字幕在线观看| 欧美 亚洲 国产 日韩一| 国产一区二区激情短视频 | 久久人人爽人人片av| 人人妻人人爽人人添夜夜欢视频| 色婷婷av一区二区三区视频| 大香蕉97超碰在线| 80岁老熟妇乱子伦牲交| 最黄视频免费看| 国产免费一区二区三区四区乱码| 日韩大片免费观看网站| 九色亚洲精品在线播放| 亚洲人成77777在线视频| 亚洲欧美日韩卡通动漫| 国产又色又爽无遮挡免| 五月玫瑰六月丁香| 亚洲精品中文字幕在线视频| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| 欧美变态另类bdsm刘玥| 久久久久久伊人网av| 亚洲精品成人av观看孕妇| 在线天堂最新版资源| 亚洲色图综合在线观看| 另类亚洲欧美激情| 日韩三级伦理在线观看| 自线自在国产av| 在线观看免费日韩欧美大片| 欧美日韩av久久| 色网站视频免费| 两个人看的免费小视频| 国产免费又黄又爽又色| 99香蕉大伊视频| 最近2019中文字幕mv第一页| 国产免费现黄频在线看| 免费av中文字幕在线| av线在线观看网站| 母亲3免费完整高清在线观看 | 国产一区二区激情短视频 | 美女中出高潮动态图| 亚洲欧美成人综合另类久久久| 麻豆精品久久久久久蜜桃| 亚洲av欧美aⅴ国产| 久久精品国产综合久久久 | 丰满乱子伦码专区| 国产色婷婷99| 亚洲av男天堂| 中文字幕免费在线视频6| freevideosex欧美| 亚洲av免费高清在线观看| 精品酒店卫生间| 免费av中文字幕在线| 黑人巨大精品欧美一区二区蜜桃 | 国产成人a∨麻豆精品| 国产伦理片在线播放av一区| 欧美精品国产亚洲| 新久久久久国产一级毛片| 五月伊人婷婷丁香| 麻豆乱淫一区二区| 男的添女的下面高潮视频| 丰满少妇做爰视频| 在线天堂中文资源库| 精品第一国产精品| 人妻 亚洲 视频| 最近最新中文字幕免费大全7| 男人舔女人的私密视频| 五月伊人婷婷丁香| 国产精品一二三区在线看| 美女大奶头黄色视频| 热99国产精品久久久久久7| 中文字幕另类日韩欧美亚洲嫩草| 卡戴珊不雅视频在线播放| 精品少妇黑人巨大在线播放| 一级毛片电影观看| 欧美97在线视频| 国产精品国产三级国产av玫瑰| 久久久久久久国产电影| 亚洲精品自拍成人| 欧美日韩av久久| 在线天堂最新版资源| 中国国产av一级| 国产又色又爽无遮挡免| av又黄又爽大尺度在线免费看| 国产日韩欧美亚洲二区| 国产国语露脸激情在线看| 激情视频va一区二区三区| 亚洲精品久久久久久婷婷小说| a级毛片黄视频| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 深夜精品福利| 久久久久人妻精品一区果冻| 国产午夜精品一二区理论片| 国产精品免费大片| 丝袜喷水一区| 99热网站在线观看| av黄色大香蕉| 七月丁香在线播放| 午夜福利视频在线观看免费| 爱豆传媒免费全集在线观看| 欧美精品高潮呻吟av久久| 亚洲精品色激情综合| 青春草视频在线免费观看| 日韩av免费高清视频| 秋霞伦理黄片| 国产亚洲午夜精品一区二区久久| 亚洲国产毛片av蜜桃av| 亚洲一级一片aⅴ在线观看| 99久久人妻综合| 如日韩欧美国产精品一区二区三区| 久久久久网色| 成年美女黄网站色视频大全免费| 最后的刺客免费高清国语| videosex国产| 亚洲国产最新在线播放| 亚洲人成77777在线视频| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 伦精品一区二区三区| 天堂8中文在线网| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 亚洲,欧美精品.| 亚洲av电影在线观看一区二区三区| 亚洲精华国产精华液的使用体验| 中国国产av一级| a级片在线免费高清观看视频| 国国产精品蜜臀av免费| 欧美变态另类bdsm刘玥| 成人影院久久| 亚洲在久久综合| 久久热在线av| 久久久久精品人妻al黑| 日本av手机在线免费观看| 久久久久久久久久成人| 亚洲精品久久午夜乱码| 久久久久精品人妻al黑| 国产成人av激情在线播放| 日韩熟女老妇一区二区性免费视频| 欧美国产精品va在线观看不卡| 欧美精品一区二区大全| 九色亚洲精品在线播放| 亚洲精品自拍成人| 国产日韩欧美视频二区| 大香蕉久久成人网| 亚洲精品国产av蜜桃|