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

    西南印度洋真光層海水中固氮細(xì)菌多樣性

    2015-01-05 01:15:26謝尚微楊俊毅張東聲王春生
    海洋學(xué)研究 2015年3期
    關(guān)鍵詞:藍(lán)細(xì)菌固氮親緣

    謝尚微,楊俊毅*,張東聲,王春生

    (1. 國家海洋局 第二海洋研究所, 浙江 杭州 310012; 2. 國家海洋局 海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點實驗室,浙江 杭州 310012)

    西南印度洋真光層海水中固氮細(xì)菌多樣性

    謝尚微1,2,楊俊毅*1,2,張東聲1,2,王春生1,2

    (1. 國家海洋局 第二海洋研究所, 浙江 杭州 310012; 2. 國家海洋局 海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點實驗室,浙江 杭州 310012)

    本文通過構(gòu)建克隆文庫和基因測序的方法研究了西南印度洋真光層海水中固氮細(xì)菌nifH基因的多樣性。從構(gòu)建的2個nifH基因克隆文庫中共獲得76條有效序列,其中46條來自CTD13-30 m文庫,分屬10個OTUs;30條來自CTD13-125 m文庫,分屬8個OTUs。系統(tǒng)發(fā)育分析結(jié)果顯示,研究站位nifH基因序列主要分布于Cluster I和Cluster III兩個分支,其中Cluster I中包含藍(lán)細(xì)菌和變形菌兩個分支,藍(lán)細(xì)菌以Group B為優(yōu)勢類群,并未獲得束毛藻和Group A的nifH基因序列。此外還有少數(shù)nifH基因序列分布于Cluster II??傮w來看,西南印度洋固氮生物基因與大西洋的親緣關(guān)系更近;固氮生物的多樣性較為豐富,受環(huán)境條件的影響,群落結(jié)構(gòu)與其它熱帶、亞熱帶寡營養(yǎng)海域具有明顯的不同。

    固氮生物;nifH基因;多樣性;西南印度洋

    0 引言

    海洋生物固氮是海洋新氮的重要來源,據(jù)估計全球海洋生物固氮量可達(dá)200 Tg N·a-1,這幾乎與陸源生物固氮的總和持平[1]。在開放性大洋環(huán)境中,生物固氮對真光層新氮輸入的貢獻(xiàn)可達(dá)50%左右[2-3],是海洋氮循環(huán)中必不可少的組成部分[4]。

    長期以來關(guān)于海洋固氮生物的研究主要集中在束毛藻Trichodesmium這一絲狀、無異型細(xì)胞的藍(lán)藻上,它們是最早發(fā)現(xiàn)的海洋固氮生物,廣泛分布在全球熱帶和亞熱帶海洋,被認(rèn)為是海洋生物固氮的主要貢獻(xiàn)者[5]。然而,基于束毛藻的生物固氮量并不能滿足海洋氮預(yù)算的平衡[6],本世紀(jì)以來,在開放性寡營養(yǎng)大洋環(huán)境中發(fā)現(xiàn)了大量的單細(xì)胞固氮細(xì)菌,研究證明這些單細(xì)胞細(xì)菌在海洋氮循環(huán)中發(fā)揮著重要作用,其對生物固氮的貢獻(xiàn)甚至可能超過束毛藻[7-8]。它們主要包括單細(xì)胞藍(lán)細(xì)菌和變形菌[9-10]。其中單細(xì)胞固氮藍(lán)細(xì)菌可分為Group A和Group B兩個分支[7-8,11-13]:Group A細(xì)胞內(nèi)缺乏產(chǎn)氧的光合系統(tǒng)II和具有固碳功能的加氧酶,有利于其在白天進(jìn)行固氮作用,在北大西洋和北太平洋都是重要的固氮生物;Group B以Crocosphaerawatsonii(WH8501)為代表種,具有藍(lán)細(xì)菌典型的光合作用機(jī)制,但其種群內(nèi)基因組的多樣性極低[14]。與藍(lán)細(xì)菌相比,變形菌也是一類重要的固氮生物,且分布更加廣泛,具有比藍(lán)細(xì)菌更高的多樣性[9]。除此之外,海洋環(huán)境中的固氮生物還包括厭氧菌(如梭狀芽孢桿菌ClostridiumPrazmowski和硫酸鹽還原菌Sulfate-Reducing Bacteria等)和古菌等,它們的多樣性及其對海洋生物固氮的貢獻(xiàn)也不容忽視[9,15-16]。

    由于多數(shù)固氮生物還不能通過分離培養(yǎng)獲得純種菌株,目前更多的是利用分子生物學(xué)技術(shù),以nifH基因為標(biāo)記物,從分子水平研究固氮生物的多樣性。CHIEN和ZINDER[17]在固氮酶基因系統(tǒng)發(fā)育分析的基礎(chǔ)上將固氮生物分成4個大的分支:Cluster I、II、III、IV。海洋環(huán)境中的固氮生物多樣性豐富,來自所有4個分支的固氮生物在海洋中均有發(fā)現(xiàn),其中Cluster I是最主要的海洋固氮生物分支,包含了藍(lán)細(xì)菌、α-變形菌、β-變形菌和γ-變形菌,在數(shù)量和多樣性上占有絕對優(yōu)勢。Cluster II包含少數(shù)變形菌和古菌,Cluster III中的固氮生物以嚴(yán)格的厭氧微生物為主,Cluster IV以產(chǎn)甲烷古菌為主[9]。

    隨著分子生物學(xué)技術(shù)的進(jìn)步,我們對海洋固氮生物的多樣性、系統(tǒng)發(fā)育和分布范圍有了全新的認(rèn)識。例如,高通量測序的結(jié)果表明,大洋表層海水固氮生物的優(yōu)勢類群可能是變形菌而不是藍(lán)細(xì)菌[18];在傳統(tǒng)觀點中不適合固氮生物生存的極地環(huán)境中也發(fā)現(xiàn)了nifH基因[19]。這些研究成果進(jìn)一步證明了固氮生物在海洋氮循環(huán)中的重要地位。

    盡管已經(jīng)得到各國學(xué)者的重視,但受采樣和研究方法等條件限制,針對固氮細(xì)菌多樣性的研究有限,主要集中在大西洋和太平洋[10,14,20],印度洋的相關(guān)報道則相當(dāng)匱乏,僅有個別報道[21-22]。西南印度洋的營養(yǎng)鹽水平較低、水文環(huán)境條件復(fù)雜[23-24],受西風(fēng)帶影響,海水溫度較低、風(fēng)力較大,并不是傳統(tǒng)觀點上海洋固氮生物的理想生境[5]。然而,據(jù)估計印度洋生物固氮可達(dá)25 Tg/a[2],其單位面積的固氮速率甚至高于大西洋[25],關(guān)于固氮生物的研究與其在該海域的重要性有明顯差距。本文對西南印度洋真光層海水固氮生物進(jìn)行研究,揭示了該海域固氮生物的多樣性和分布特征,為進(jìn)一步評估海洋固氮生物的分布及其對新氮的貢獻(xiàn)提供科學(xué)資料。

    1 材料與方法

    1.1 研究海區(qū)和采樣分析方法

    樣品由“大洋一號”科學(xué)考察船2010年1月采集自西南印度洋(21V-CTD13站位:37.789 4°S、49.636 7°E,圖1)。使用Seabird 911Plus采水器分別采集30 m層和125 m層海水水樣。每層取4 L水樣,用200 μm篩絹預(yù)過濾后,過濾于孔徑為0.2 μm的核孔濾膜上,濾膜保存于-20 ℃,帶回實驗室分析。

    圖1 西南印度洋固氮生物采樣站位Fig.1 Sampling station in the southwest Indian Ocean

    1.2 DNA提取

    參照LANGLOIS et al[12]的方法,使用QIAGEN DNA提取試劑盒(DNeasy Plant mini kit)提取濾膜上的環(huán)境總DNA。

    1.3 nifH基因擴(kuò)增與克隆

    以從濾膜上提取的環(huán)境總DNA為模板,對nifH基因進(jìn)行巢式PCR,兩輪引物分別為:第一輪引物[26]:nifH1: 5’-TTYTAYGGNAARGGNGG-3’,nifH2: 5’-ATRTTRTTNGCNGCRTA-3’;第二輪引物[27]:nifH3: 5’-TGYGAYCCNAARGCNG A-3’,nifH4: 5’-ADNGCCATCATYTCNCC-3’。50 μL PCR反應(yīng)體系:Mg2+(80 nmol/L)、dNTPs(400 μmol/L)、上下游引物(1 μmol/L)、Taq酶(2.5 U)、10×PCR緩沖液 5 μL、DNA模板2 μL,無菌水補(bǔ)充至50 μL。PCR反應(yīng)條件如下:94 ℃變性4 min;30個循環(huán):94 ℃變性60 s,57 ℃退火30 s,72 ℃延伸40 s;72 ℃延伸10 min。第一輪PCR產(chǎn)物作為DNA模板進(jìn)行第二輪PCR反應(yīng),兩輪PCR反應(yīng)條件相同。PCR產(chǎn)物使用1%瓊脂糖凝膠電泳檢驗。

    nifH基因的PCR產(chǎn)物經(jīng)凝膠回收試劑盒純化后,連接到pMD20-T Vector(TaKaRa)上,轉(zhuǎn)入大腸桿菌感受態(tài)細(xì)胞(E.coli JM109 Competent cells, TaKaRa),涂布到含有氨芐青霉素(Ampicillin)/IPTG/X-Gal的LB固體培養(yǎng)基上,37 ℃培養(yǎng)過夜,挑取白色克隆子,4 ℃保存。利用T7和SP6進(jìn)行菌落PCR,篩選插入目的片段的陽性克隆子并直接利用菌液進(jìn)行DNA序列測定,構(gòu)建克隆文庫。測序反應(yīng)由上海生工生物工程公司完成。

    1.4 序列分析

    將測序所得序列在NCBI中進(jìn)行同源性比對,挑選與序列親緣關(guān)系最相似序列,利用CHECK-CHIMERA軟件分析檢查去除嵌合體[28]。使用Clustal X 1.83和Mega 4.0分析多重對比結(jié)果,構(gòu)建系統(tǒng)發(fā)育樹。將nifH基因序列提交GenBank數(shù)據(jù)庫,登記號為KF740730—KF740805。

    1.5 環(huán)境參數(shù)分析

    溫度和鹽度數(shù)據(jù)由Seabird 911Plus CTD直接測得。用于葉綠素a分析的樣品采集自0、10、20、30、50、75、100、125、150、175和200 m層,每層水樣分別取500 mL,依次用孔徑為20、2和0.2 μm的濾膜(Millipore NY2004700、TTTP04700、GTTP04700)分級過濾,分析測定細(xì)胞粒徑大于20、2~20和0.2~2 μm三個粒級葉綠素a質(zhì)量濃度,各粒級葉綠素a質(zhì)量濃度之和為總?cè)~綠素a質(zhì)量濃度。葉綠素a測定采用熒光法[29]。

    2 結(jié)果

    2.1 溫度、鹽度及葉綠素a質(zhì)量濃度的垂直分布

    調(diào)查站位200 m以淺水層溫、鹽垂直分布較均勻,水溫在16~19 ℃之間,隨水深增大而逐漸降低;鹽度在35.5~35.6之間,垂直分布呈單峰形態(tài),高值在水深40 m左右(圖2)。

    圖2 CTD13站位溫、鹽垂直分布Fig.2 Vertical distribution of temperature and salinity on CTD13 station

    葉綠素a質(zhì)量濃度的垂直分布呈典型的單峰形態(tài)(圖3),總?cè)~綠素a質(zhì)量濃度最大值(0.336 μg/L)位于100 m層,30 m層和125 m層的總?cè)~綠素a質(zhì)量濃度相近,分別為0.106 μg/L和0.134 μg/L。微微型浮游植物(Pico,0.2~2 μm)對總?cè)~綠素a質(zhì)量濃度的貢獻(xiàn)占絕對優(yōu)勢,平均達(dá)到70.2%,垂直分布特征也與總?cè)~綠素a質(zhì)量濃度一致;其次是微型浮游植物(Nano,2~20 μm),占總?cè)~綠素a質(zhì)量濃度的27.4%;小型浮游植物(Micro,>20 μm)占總?cè)~綠素a質(zhì)量濃度比例最低,僅2.4%。

    圖3 CTD13站位各粒級葉綠素a質(zhì)量濃度垂直分布Fig.3 Vertical distribution of Chlorophyll a concentration of different size on CTD13 station

    2.2 nifH基因克隆文庫構(gòu)建和多樣性分析

    構(gòu)建CTD13-30 m和CTD13-125 m兩個nifH基因文庫,對獲得的陽性克隆進(jìn)行測序,共獲得76條有效的nifH基因,其中CTD13-30 m文庫46條,CTD13-125 m文庫30條。以基因相似性≥97%為一個OTU統(tǒng)計,CTD13-30 m文庫獲得10個OTUs,CTD13-125 m文庫獲得8個OTUs,稀釋度曲線(圖4)顯示這2個文庫均已趨于飽和。

    系統(tǒng)發(fā)育分析結(jié)果顯示,CTD13-30 m和CTD13-125 m兩個文庫中的nifH基因分屬于Cluster I、Cluster II和Cluster III這3個已知的分支。CTD13-30 m文庫中的nifH基因在3個分支中均有分布,CTD13-125 m文庫中的nifH基因僅分布于Cluster I和Cluster III兩個分支(圖5和圖6)。

    屬于Cluster I的nifH基因大多數(shù)與固氮藍(lán)細(xì)菌親緣關(guān)系較近,其中25條nifH基因?qū)儆贕roup B分支,該分支以Crocosphaerawatsonii為代表,是寡營養(yǎng)海域單細(xì)胞固氮藍(lán)細(xì)菌的一個重要分支[7]。另有3條nifH基因不屬于藍(lán)細(xì)菌nifH基因已知的分支,CTD13-30 m-56與色球藻屬Chroococcus和從現(xiàn)代海洋疊層石中獲取的未培養(yǎng)nifH基因親緣關(guān)系較近,CTD13-125 m-21與來自中國南海海水中的未培養(yǎng)nifH基因親緣關(guān)系較近,CTD13-125 m-13則與來自土壤和活性污泥等環(huán)境中的未培養(yǎng)nifH基因親緣關(guān)系較近(圖5)。4條nifH基因與γ變形菌的親緣關(guān)系較近,來自CTD13-125 m文庫的1條nifH基因與Vibriodiazotrophicus親緣關(guān)系較近,來自CTD13-30 m文庫的3條nifH基因與太平洋和大西洋寡營養(yǎng)海水中的未培養(yǎng)nifH基因親緣關(guān)系較近(圖5)。

    圖4 CTD13-30 m和CTD13-125 m nifH基因文庫稀釋度曲線Fig.4 Rarefaction curves of OTUs of uncultured nifH clones

    41條nifH基因序列屬于Cluster III,該分支包含厭氧的硫還原菌(如δ變形菌)、梭狀芽胞桿菌和紫硫菌(Purple Sulfur Bacteria)等。Cluster III又可以分成2個聚類分支,1支包含CTD13-125 m文庫的8條nifH基因,它們與來自開放性大洋和深海沉積物的未培養(yǎng)nifH基因親緣關(guān)系較近,另一支包含CTD13-30 m文庫的23條nifH基因和CTD 13-125 m文庫的10條nifH基因,它們與來自活性污泥環(huán)境中的未培養(yǎng)nifH基因親緣關(guān)系較近(圖6)。

    Cluster II僅包含3條CTD13-30 m文庫的nifH基因,也可分成2個分支,CTD13-30 m-2和CTD13-30 m-44與來自波羅的海的未培養(yǎng)nifH基因親緣關(guān)系較近,CTD13-30 m-13與密螺旋體的nifH基因親緣關(guān)系較近(圖6)。

    3 討論與結(jié)論

    本文結(jié)果顯示西南印度洋海域固氮細(xì)菌nifH基因具有較為獨特的多樣性和群落結(jié)構(gòu)特征。與熱帶、亞熱帶的太平洋[7]和大西洋[30]海域相比,本文構(gòu)建的克隆文庫中未獲得來自束毛藻的nifH基因,這可能是受溫度限制的結(jié)果。溫度是限制固氮細(xì)菌生長和分布的一個重要因素,例如束毛藻,主要分布在水溫大于20 ℃的海域[31],相比較而言,單細(xì)胞固氮生物對溫度的耐受性更強(qiáng),是低溫環(huán)境下固氮生物的優(yōu)勢類群[12]。調(diào)查站位200 m以淺水層溫度在16~19 ℃之間(圖2),30 m和125 m層的溫度分別為17.9 ℃和16.6 ℃,均低于束毛藻的適宜溫度下限,不利于束毛藻的生長繁殖。

    Group B是本文研究站位固氮生物的主要優(yōu)勢類群,溫度較高的CTD13-30 m文庫中共獲得16條屬于Group B的nifH基因,占其文庫百分比的34.8%,溫度較低的CTD13-125 m文庫獲得9條屬于Group B的nifH基因,占文庫總序列數(shù)的30.0%,前者明顯高于后者。有研究表明,單細(xì)胞固氮藍(lán)細(xì)菌Group B的數(shù)量與所在環(huán)境水溫之間存在顯著的正相關(guān)[13],溫度越高數(shù)量越多。這與本文的結(jié)果一致,說明本文調(diào)查站位中Group B的數(shù)量與溫度可能也存在正相關(guān)關(guān)系。

    固氮生物是一類全球性分布的生物,LANGLOIS et al[30]認(rèn)為熱帶、亞熱帶大西洋和太平洋的單細(xì)胞藍(lán)細(xì)菌和Cluster III的nifH基因具有很高的同源性。BIRD et al[32]發(fā)現(xiàn)北印度洋阿拉伯海的變形菌nifH基因與大西洋和太平洋的變形菌nifH基因相似性很高。本文獲得的Group BnifH基因與大西洋Group B的非培養(yǎng)nifH基因相似性較高(98%),與太平洋Group B的非培養(yǎng)nifH基因相似性較低(93%),這意味著,西南印度洋Group B固氮藍(lán)細(xì)菌與大西洋固氮藍(lán)細(xì)菌具有更近的親緣關(guān)系。來自變形菌分支的nifH基因與大西洋和太平洋的變形菌nifH基因都具有很高的相似性,而Cluster III的nifH基因與大西洋和太平洋同一分支中的nifH基因相似性都很低。本文還獲得了屬于Cluster II分支的3條nifH基因,這一分支很少出現(xiàn)在熱帶、亞熱帶的太平洋和大西洋,在深海沉積物[33]、熱液噴口[34]等海洋環(huán)境中更為常見,受這一環(huán)境特征影響,可能與產(chǎn)甲烷細(xì)菌具有較近的親緣關(guān)系[34]。

    盡管固氮生物廣泛分布在全球海洋當(dāng)中,但不同海域固氮生物的群落結(jié)構(gòu)和多樣性存在顯著差異[18]。在熱帶、亞熱帶北大西洋,以束毛藻為主的固氮藍(lán)細(xì)菌在固氮生物中占絕對優(yōu)勢,是生物固氮的主要來源[12,35];在亞熱帶北太平洋單細(xì)胞藍(lán)細(xì)菌(Group A)則成為固氮生物的優(yōu)勢類群,同時也是生物固氮的主要來源[8,11]。本文的研究結(jié)果表明西南印度洋固氮細(xì)菌群落結(jié)構(gòu)與西南太平洋[18]相似,以藍(lán)細(xì)菌和Cluster III的單細(xì)胞固氮細(xì)菌為主,這可能是由兩者之間相似的環(huán)境特征造成的。盡管還缺乏更為直接的證據(jù),本文結(jié)果仍表明Cluster III可能在西南印度洋的生物固氮和單循環(huán)中扮演重要角色。實際的固氮酶表達(dá)活性和細(xì)菌豐度仍需要進(jìn)一步研究,以證明Cluster III固氮細(xì)菌的固氮能力。

    圖5 nifH基因Cluster I分支的系統(tǒng)發(fā)育樹Fig.5 Neighbor-joining phylogenetic tree of Cluster I nifH gene sequences

    圖6 nifH基因Cluster II和Cluster III分支的系統(tǒng)發(fā)育樹Fig.6 Neighbor-joining phylogenetic tree of Cluster II and Cluster III nifH gene sequences

    綜上所述,西南印度洋固氮細(xì)菌nifH基因與其他海區(qū)相比具有較為獨特的群落結(jié)構(gòu),Group B和Cluster III是固氮細(xì)菌的主要優(yōu)勢類群。從系統(tǒng)發(fā)育上來看,該海域的nifH基因與大西洋更為接近。

    致謝 感謝“大洋一號”全體船員在采樣過程中提供的幫助,感謝國家海洋局第二海洋研究所王淵助理研究員提供溫、鹽數(shù)據(jù)。

    [1] GALLOWAY J, DENTENER F J, CAPONE D G, et al. Nitrogen cycles: Past, present and future[J]. Biogeochemistry,2004,70(2):153-226.

    [2] DEUTSCH C A, GRUBER N P, KEY R M, et al. Denitrification and N2fixation in the Pacific Ocean[J]. Global Biogeochem Cycles,2001,15(2):483-506.

    [3] MICHAELS A F, BATES N R, BUESSELER K O, et al. Carbon-cycle imbalances in the Sargasso Sea[J]. Nature,1994,372(6506):537-540.

    [4] KARL D, MICHAELS A, BERGMAN B, et al. Dinitrogen fixation in the world’s oceans[J]. Biogeochemistry,2002,57(1):47-98.

    [5] CAPONE D G, ZEHR J P, PAERL H W, et al.Trichodesmium: a globally significant marine cyanobacterium[J]. Science,1997,276(5316):1 221-1 229.

    [6] ZEHR J P, CARPENTER E J, VILLAREAL T A. New perspectives on nitrogen-fixing microorganisms in tropical and subtropical oceans[J]. Trends Microbiol,2000,8(2):68-73.

    [7] ZEHR J P, WATERBURY J B, TURNER P J, et al. Unicellular cyanobacteria fix N2in the subtropical North Pacific Ocean[J]. Nature,2001,412(6847):635-638.

    [8] MONTOYA J P, HOLL C M, ZEHR J P, et al. High rates of N2fixation by unicellular diazotrophs in the oligotrophic Pacific Ocean[J]. Nature,2004,430:1 027-1 032.

    [9] ZEHR J P, JENKINS B D, SHORT S M, et al. Nitrogenase gene diversity and microbial community structure: a cross-system comparison[J]. Environmental Microbiology,2003,5(7):539-554.

    [10] FALCON L I, CARPENTER E J, CIPRIANO F, et al. N2fixation by unicellular bacterioplankton from the Atlantic and Pacific Oceans: Phylogeny and in situ rates[J]. Appl Environ Microbiol,2004,70(2):765-770.

    [11] CHURCH M J, BJORKMAN K M, KARL D M. Regional distribution of nitrogen-fixing bacteria in the Pacific Ocean[J]. Limnol Oceanogr,2008,53(1):63-77.

    [12] LANGLOIS R J, HUMMER D, LAROCHE J. Abundances and distribution of dominantnifHphylotypes in the northern Atlantic Ocean[J]. Appl Environ Microbiol,2008,74(6):1 922-1 931.

    [13] MOISANDER P H, BEINART R A, HEWSON I, et al. Unicellular cyanobacterial distributions broaden the oceanic N2fixation domain[J]. Science,2010,327(5972):1 512-1 514.

    [14] ZEHR J P, MONTOYA J P, JENKINS B D, et al. Experiments linking nitrogenase gene expression to nitrogen fixation in the North Pacific subtropical gyre[J]. Limnol Oceanogr,2007,52(1):169-183.

    [15] BURNS J A, ZEHR J P, CAPONE D G. Nitrogen fixing phylotypes of Chesapeake Bay and Neuse River Estuary sediments[J]. Microb Ecol,2002,44(4):336-343.

    [16] DANG Hong-yue, LUAN Xi-wu, ZHAO Jing-yi, et al. Diverse and novelnifHandnifH-like genes sequences in the deep-sea methane seep sediments of the Okhotsk Sea[J]. Appl Environ Microbiol,2009,75(7):2 238-2 245.

    [17] CHIEN Y T, ZINDER S H. Cloning, functional organization, transcript studies, and phylogenetic analysis of the complete nitrogenase structural genes (nifHDK2) and associated genes in the archaeonMethanosarcinabarkeri227[J]. J Bacteriol,1996,178(1):143-148.

    [18] FARNELID H, ANDERSSON A F, BERTILSSON S, et al. Nitrogenase gene amplicons from global marine surface waters are dominated by genes of non-cyanobacteria[J]. Plos One,2011,6(4):e19223.

    [19] POMMIER T, PINHASSI J, HAGSTROM A. Biogeographic analysis of ribosomal RNA clusters from marine bacterioplankton[J]. Aquat Microb Ecol,2005,41(1):79-89.

    [20] KITAJIMA S, FURUYA K, HASHIHAMA F, et al. Latitudinal distribution of diazotrophs and their nitrogen fixation in the tropical and subtropical western North Pacific[J]. Limnol Oceanogr,2009,54(2):537-547.

    [21] BIRD C, MARTINEZ J M, O’DONNEL A G, et al. Spatial distribution and transcriptional activity of an uncultured clade of planktonic diazotrophic γ-Proteobacteria in the Arabian Sea[J]. Appl Environ Microbiol,2005,71(4):2 079-2 085.

    [22] POULTON A J, STINCHCOMBE M C, QUARTLY G D. High numbers ofTrichodesmiumand diazotrophic diatoms in the southwest Indian Ocean[J]. Geophysical Research Letters,2009,36(15):139-156.

    [23] RAJAKUMAR A, ALAGARSAMY R, KHARE N, et al. Studies on the nutrient distribution in the Southern Ocean waters along the 45°E transect[J]. Indian Journal of Marine Sciences,2008,37(4):424-429.

    [24] READ J F, LUCAS M I, HOLLEY S E, et al. Phytoplankton, nutrients and hydrography in the frontal zone between the Southwest Indian Subtropical gyre and the Southern Ocean[J]. Deep-Sea Research I,2000,47(12):2 341-2 367.

    [25] WAJIH S, NAQVI A. Nitrogen in the marine environment[M]. Second Edition. Elsevier,2008:631-671.

    [26] ZEHR J P, MCREYNOLDS L A. Use of degenerate oligonucleotides for amplification of thenifHgene from the marine cyanobacteriumTrichodesmiumthiebautii[J]. Appl Environ Microbiol,1989,55(10):2 522-2 526.

    [27] ZANI S, MELLON M T, COLLIER J L, et al. Expression ofnifHgenes in natural microbial assemblages in Lake George, New York, detected by reverse transcriptase PCR[J]. Appl Environ Microbiol,2000,66(7):3 119-3 124.

    [28] COLE J R, CHAI B, MARSH T L, et al. The Ribosomal Database Project (RDP-II): Previewing a new autoaligner that allows regular updates and the new prokaryotic taxonomy[J]. Nucleic Acids Research,2003,31(1):442-443.

    [29] GB/T 12763.6-2007 Specifications for oceanographic survey[S]. Beijing: China Standard Press,2007. GB/T 12763.6—2007 海洋調(diào)查規(guī)范[S].北京:中國標(biāo)準(zhǔn)出版社,2007.

    [30] LANGLOIS R J, LAROCHE J, RAAB P A. Diazotrophic diversity and distribution in the tropical and subtropical Atlantic Ocean[J]. Appl Environ Microbiol,2005,71(12):7 910-7 919.

    [31] BREITBARTH E, OSCHLIES A, LAROCHE J. Physiological constraints on the global distribution ofTrichodesmium-effect of temperature on diazotrophy[J]. Biogeosciences Discussions,2006,3(3):779-801.

    [32] BIRD C, MARTINEZ J M, O’DONNELL A G, et al. Spatial distribution and transcriptional activity of an uncultured clade of plankton diazotrophic γ-proteobacteria in the Arabian Sea[J]. Appl Environ Microbiol,2005,71(4):2 079-2 085.

    [33] DANG Hong-yue, YANG Jin-ying, LI Jing, et al. Environment-dependent distribution of the sedimentnifH-harboring microbiota in the Northern South China Sea[J]. Appl Environ Microbiol,2013,79(1):121-132.

    [34] MEHTA M P, BUTTERFIELD D A, BAROSS J A. Phylogenetic diversity of nitrogenase (nifH) genes in deep-sea and hydrothermal vent environments of the Juan de Fuca Ridge[J]. Appl Environ Microbiol,2003,69(2):960-970.

    [35] FOSTER R A, SUBRAMANIAM A, ZEHR J P. Distribution and activity of diazotrophs in the Eastern Equatorial Atlantic[J]. Environ Microb,2009,11(4):741-750.

    Phylogenetic diversity ofnifHgenes in the euphotic zone of the southwest Indian Ocean

    XIE Shang-wei1,2, YANG Jun-yi*1,2, ZHANG Dong-sheng1,2, WANG Chun-sheng1,2

    (1.TheSecondInstituteofOceanography,SOA,Hangzhou310012,China; 2.LaboratoryofMarineEcosystemandBiogeochemistry,SOA,Hangzhou310012,China)

    Our knowledge on diazotrophs is limited in the Indian Ocean, which is considered an important source of biological N2fixation of the global ocean. Phylogenetic diversity ofnifHgenes in the southwest Indian Ocean was studied through cloning and sequencing methods. After sequences were checked, twonifHgene libraries were established. Overall, 76nifHgene sequences are recovered from Station CTD13, including 46nifHgene sequences of 10 OTUs from libraries CTD13-30 m, and 30nifHgene sequences of 8 OTUs from libraries CTD13-125 m, respectively. Results of phylogenetic analysis show that mostnifHgene sequences fall into Cluster I and III, only 3 recoverednifHgene sequences are grouped with Cluster II. In Cluster I, most ofnifHgene sequences cluster with unicellular cyanobacteria Group B, while in Cluster III, mostnifHgene sequences are closely related with uncultured clone from active sludge. NonifHgene sequences clustered withTrichodesmiumare recovered, which may be attributed to low water temperature. In general, diazotroph community in the southwest Indian Ocean is different from those in tropical and subtropical oligotrophic ocean.

    diazotroph;nifHgene; diversity; the southwest Indian Ocean

    10.3969/j.issn.1001-909X.2015.03.008.

    2015-03-24

    2015-05-19

    國家自然科學(xué)基金項目資助(41206104); 國際海域資源調(diào)查與開發(fā)“十二五”課題資助(DY125-11-E-03)

    謝尚微(1990-),女,浙江溫州市人,主要從事海洋生物學(xué)方面的研究。E-mail:410074144@qq.com

    *通訊作者:楊俊毅(1970-),男,研究員,主要從事海洋生物與生態(tài)環(huán)境監(jiān)測技術(shù)方面的研究。E-mail:junyiyang@sio.org.cn

    Q938.1

    A

    1001-909X(2015)03-0054-08

    10.3969/j.issn.1001-909X.2015.03.008

    謝尚微,楊俊毅,張東聲,等. 西南印度洋真光層海水中固氮細(xì)菌多樣性[J]. 海洋學(xué)研究,2015,33(3):54-61,

    XIE Shang-wei, YANG Jun-yi, ZHANG Dong-sheng, et al. Phylogenetic diversity ofnifHgenes in the euphotic zone of the southwest Indian Ocean[J]. Journal of Marine Sciences, 2015,33(3):54-61, doi:10.3969/j.issn.1001-909X.2015.03.008.

    猜你喜歡
    藍(lán)細(xì)菌固氮親緣
    巧用問題串 拆解藍(lán)細(xì)菌
    谷子近緣野生種的親緣關(guān)系及其利用研究
    土壤中的天然化肥廠
    中國醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    生命的開路先鋒藍(lán)細(xì)菌
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    殖民火星,藍(lán)細(xì)菌來幫忙
    小白菜種質(zhì)遺傳多樣性與親緣關(guān)系的SRAP 和SSR分析
    杉木與固氮樹種混交對土壤有機(jī)質(zhì)及氮含量的影響
    土壤中的天然化肥廠
    ——固氮微生物
    最近最新中文字幕免费大全7| 亚洲人与动物交配视频| 国产乱人偷精品视频| 热99久久久久精品小说推荐| 亚洲av欧美aⅴ国产| 成人国产麻豆网| 精品人妻在线不人妻| 夫妻性生交免费视频一级片| 日本免费在线观看一区| 中文欧美无线码| 最近最新中文字幕免费大全7| 国产精品国产av在线观看| a级片在线免费高清观看视频| 高清视频免费观看一区二区| 日本黄大片高清| 国产不卡av网站在线观看| 精品酒店卫生间| 肉色欧美久久久久久久蜜桃| 婷婷色综合大香蕉| 国产欧美日韩综合在线一区二区| 亚洲欧洲国产日韩| 国产日韩一区二区三区精品不卡| 国内精品宾馆在线| 青春草国产在线视频| 国产在线一区二区三区精| 久久国产精品男人的天堂亚洲 | 久久99蜜桃精品久久| 极品人妻少妇av视频| 欧美国产精品一级二级三级| 午夜福利影视在线免费观看| 国产精品久久久久久精品古装| 涩涩av久久男人的天堂| 熟女av电影| 人体艺术视频欧美日本| 国产精品成人在线| 啦啦啦在线观看免费高清www| 亚洲精品,欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 色哟哟·www| 国产av国产精品国产| 女性被躁到高潮视频| 欧美成人精品欧美一级黄| 一个人免费看片子| 欧美激情 高清一区二区三区| 日韩欧美一区视频在线观看| 99精国产麻豆久久婷婷| 国产一区二区三区av在线| 一本—道久久a久久精品蜜桃钙片| 久久ye,这里只有精品| 精品国产国语对白av| 国产 精品1| 18禁观看日本| 亚洲国产精品国产精品| 国产一区有黄有色的免费视频| 高清视频免费观看一区二区| 久久久久久久久久久久大奶| 亚洲成国产人片在线观看| 综合色丁香网| 国产老妇伦熟女老妇高清| 天美传媒精品一区二区| 一级片免费观看大全| 亚洲伊人色综图| 国产老妇伦熟女老妇高清| 热re99久久国产66热| 曰老女人黄片| av在线老鸭窝| 亚洲av电影在线进入| 一区二区三区精品91| av免费在线看不卡| 777米奇影视久久| 青青草视频在线视频观看| 99久久精品国产国产毛片| 婷婷成人精品国产| 亚洲欧美成人精品一区二区| 免费大片黄手机在线观看| 母亲3免费完整高清在线观看 | 成年av动漫网址| 国产熟女午夜一区二区三区| 午夜视频国产福利| 欧美变态另类bdsm刘玥| 成人午夜精彩视频在线观看| av视频免费观看在线观看| videos熟女内射| 欧美精品高潮呻吟av久久| 久久久久久久大尺度免费视频| 日韩制服丝袜自拍偷拍| 高清毛片免费看| 麻豆精品久久久久久蜜桃| a级毛色黄片| 久久久国产欧美日韩av| 日本黄色日本黄色录像| 欧美精品人与动牲交sv欧美| 色哟哟·www| av电影中文网址| 久久久精品区二区三区| 久久久久久久久久成人| 中文欧美无线码| 亚洲成人av在线免费| 少妇高潮的动态图| 国产av一区二区精品久久| 妹子高潮喷水视频| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| 免费日韩欧美在线观看| 91精品伊人久久大香线蕉| 91精品伊人久久大香线蕉| av一本久久久久| 啦啦啦中文免费视频观看日本| 欧美精品一区二区免费开放| 国产片内射在线| 中文字幕精品免费在线观看视频 | av在线app专区| 精品少妇黑人巨大在线播放| 国产男女超爽视频在线观看| 久久这里只有精品19| 国产综合精华液| 欧美精品亚洲一区二区| 久久国产精品男人的天堂亚洲 | 免费久久久久久久精品成人欧美视频 | 亚洲情色 制服丝袜| 欧美亚洲日本最大视频资源| 一级毛片黄色毛片免费观看视频| 国产亚洲欧美精品永久| 999精品在线视频| 女人久久www免费人成看片| av片东京热男人的天堂| 亚洲四区av| 国产精品国产av在线观看| 免费观看a级毛片全部| 国产1区2区3区精品| 国产乱人偷精品视频| 国产免费一级a男人的天堂| 欧美国产精品一级二级三级| 内地一区二区视频在线| 母亲3免费完整高清在线观看 | 国语对白做爰xxxⅹ性视频网站| 成人综合一区亚洲| 欧美日韩视频精品一区| 亚洲成人手机| 日韩在线高清观看一区二区三区| 亚洲色图 男人天堂 中文字幕 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲美女视频黄频| 欧美国产精品va在线观看不卡| 久久精品aⅴ一区二区三区四区 | 又粗又硬又长又爽又黄的视频| 日韩 亚洲 欧美在线| 亚洲精品一二三| 又大又黄又爽视频免费| 国产 精品1| 亚洲久久久国产精品| 亚洲伊人久久精品综合| 少妇熟女欧美另类| 美女脱内裤让男人舔精品视频| 国产精品免费大片| 欧美成人午夜免费资源| 母亲3免费完整高清在线观看 | 亚洲四区av| 18+在线观看网站| 波多野结衣一区麻豆| 大片免费播放器 马上看| 一级爰片在线观看| 91aial.com中文字幕在线观看| 久久 成人 亚洲| 国产色爽女视频免费观看| 一级片免费观看大全| 久久精品久久精品一区二区三区| 成年女人在线观看亚洲视频| 一区二区av电影网| 少妇人妻精品综合一区二区| 一级毛片黄色毛片免费观看视频| xxxhd国产人妻xxx| 在线观看人妻少妇| 91午夜精品亚洲一区二区三区| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 蜜臀久久99精品久久宅男| 日日爽夜夜爽网站| 国产av一区二区精品久久| 日韩伦理黄色片| 一本大道久久a久久精品| 80岁老熟妇乱子伦牲交| 极品人妻少妇av视频| 日韩 亚洲 欧美在线| 久久久久国产网址| 97在线视频观看| 飞空精品影院首页| 天天操日日干夜夜撸| 天堂8中文在线网| 十八禁高潮呻吟视频| 伊人亚洲综合成人网| 91在线精品国自产拍蜜月| 久久久国产欧美日韩av| 中文字幕人妻熟女乱码| 久久人人爽人人片av| 免费高清在线观看日韩| 中文字幕最新亚洲高清| 国产色爽女视频免费观看| 超色免费av| 大片免费播放器 马上看| 精品人妻在线不人妻| 国产欧美日韩综合在线一区二区| 久久久久精品性色| 久久女婷五月综合色啪小说| 国国产精品蜜臀av免费| 国产精品麻豆人妻色哟哟久久| 国产黄频视频在线观看| 午夜老司机福利剧场| 少妇被粗大的猛进出69影院 | 成人二区视频| 亚洲国产精品一区二区三区在线| 久久综合国产亚洲精品| 国产淫语在线视频| 草草在线视频免费看| 51国产日韩欧美| 国产白丝娇喘喷水9色精品| 精品久久国产蜜桃| 国产乱人偷精品视频| 欧美亚洲 丝袜 人妻 在线| 咕卡用的链子| 少妇高潮的动态图| 黑人欧美特级aaaaaa片| 日韩欧美一区视频在线观看| 爱豆传媒免费全集在线观看| 天天躁夜夜躁狠狠躁躁| 在线观看免费日韩欧美大片| 精品一区二区三区视频在线| 国产av码专区亚洲av| 成人漫画全彩无遮挡| av黄色大香蕉| 视频区图区小说| 人妻少妇偷人精品九色| 熟女电影av网| 91精品国产国语对白视频| 欧美变态另类bdsm刘玥| 啦啦啦在线观看免费高清www| 欧美亚洲日本最大视频资源| 国产精品国产三级专区第一集| 亚洲综合精品二区| 国产在线免费精品| 观看美女的网站| 一级爰片在线观看| 亚洲国产av新网站| 亚洲精品成人av观看孕妇| 国产av精品麻豆| 性色avwww在线观看| 久久毛片免费看一区二区三区| 国产麻豆69| 亚洲人与动物交配视频| 成人免费观看视频高清| 伦理电影免费视频| 国产乱人偷精品视频| 两性夫妻黄色片 | 五月开心婷婷网| 欧美国产精品一级二级三级| 午夜影院在线不卡| 七月丁香在线播放| 热re99久久精品国产66热6| 免费观看在线日韩| 亚洲婷婷狠狠爱综合网| 岛国毛片在线播放| 欧美少妇被猛烈插入视频| 女的被弄到高潮叫床怎么办| 亚洲一区二区三区欧美精品| 久久久久精品性色| 在线看a的网站| 在线观看免费高清a一片| 欧美老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 亚洲美女搞黄在线观看| 波多野结衣一区麻豆| 亚洲欧美一区二区三区国产| 人成视频在线观看免费观看| www日本在线高清视频| 国产精品一区二区在线观看99| 制服人妻中文乱码| 久久久欧美国产精品| 一本—道久久a久久精品蜜桃钙片| 晚上一个人看的免费电影| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 最近中文字幕2019免费版| 人人澡人人妻人| 亚洲国产看品久久| 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 美女主播在线视频| 国产精品久久久久久久久免| 91国产中文字幕| 久久久久网色| 免费观看av网站的网址| 亚洲欧美清纯卡通| 国产熟女欧美一区二区| 国产免费一区二区三区四区乱码| 九草在线视频观看| 亚洲精品久久午夜乱码| 精品人妻熟女毛片av久久网站| 夫妻午夜视频| 国产色婷婷99| 中文精品一卡2卡3卡4更新| 高清毛片免费看| 午夜av观看不卡| 一级黄片播放器| 中文天堂在线官网| 午夜激情av网站| 成年美女黄网站色视频大全免费| 99久久精品国产国产毛片| 久久精品人人爽人人爽视色| 91精品伊人久久大香线蕉| 久久久久精品人妻al黑| 成人18禁高潮啪啪吃奶动态图| 亚洲丝袜综合中文字幕| 国产成人精品婷婷| 美女中出高潮动态图| 亚洲综合色惰| 91成人精品电影| 国产麻豆69| 纯流量卡能插随身wifi吗| 91aial.com中文字幕在线观看| 97精品久久久久久久久久精品| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 国产欧美亚洲国产| 久久 成人 亚洲| 秋霞伦理黄片| 精品久久国产蜜桃| 国产av一区二区精品久久| 又粗又硬又长又爽又黄的视频| 久久久国产欧美日韩av| 日韩中文字幕视频在线看片| 国产在线一区二区三区精| 免费在线观看黄色视频的| 美女视频免费永久观看网站| 制服诱惑二区| 欧美成人精品欧美一级黄| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| 午夜日本视频在线| 久久精品夜色国产| 亚洲av国产av综合av卡| 99热这里只有是精品在线观看| 纵有疾风起免费观看全集完整版| 国产免费视频播放在线视频| 国产成人午夜福利电影在线观看| 激情五月婷婷亚洲| 在线观看人妻少妇| 日本午夜av视频| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 黄色 视频免费看| 国产一区二区三区av在线| 久久久精品区二区三区| 一区在线观看完整版| 婷婷色av中文字幕| 日本色播在线视频| 久久久久久人妻| 又大又黄又爽视频免费| 精品亚洲乱码少妇综合久久| 极品人妻少妇av视频| 国产深夜福利视频在线观看| 七月丁香在线播放| 巨乳人妻的诱惑在线观看| av又黄又爽大尺度在线免费看| 国产亚洲精品第一综合不卡 | 日韩精品免费视频一区二区三区 | 韩国精品一区二区三区 | 亚洲第一av免费看| 国产激情久久老熟女| h视频一区二区三区| 久久人人爽人人爽人人片va| 日韩视频在线欧美| 777米奇影视久久| 日韩视频在线欧美| 男女国产视频网站| 高清视频免费观看一区二区| 午夜免费鲁丝| 久久精品国产亚洲av涩爱| xxxhd国产人妻xxx| 国产在线视频一区二区| 热re99久久国产66热| 国产探花极品一区二区| 亚洲精品av麻豆狂野| 伦理电影大哥的女人| 国产av国产精品国产| 精品酒店卫生间| 亚洲成人手机| 免费观看av网站的网址| 国产 精品1| 视频在线观看一区二区三区| 日韩一区二区三区影片| 日本午夜av视频| 狠狠精品人妻久久久久久综合| 日韩不卡一区二区三区视频在线| 久久久国产欧美日韩av| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 男女高潮啪啪啪动态图| 亚洲,欧美精品.| 久久热在线av| 精品酒店卫生间| 大片免费播放器 马上看| 亚洲成人av在线免费| 一本久久精品| 国产成人精品无人区| 九色亚洲精品在线播放| 亚洲美女黄色视频免费看| 欧美成人精品欧美一级黄| 免费播放大片免费观看视频在线观看| 亚洲av在线观看美女高潮| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| 中文精品一卡2卡3卡4更新| 久久这里有精品视频免费| av天堂久久9| 在线精品无人区一区二区三| 热re99久久国产66热| 亚洲成色77777| 免费观看无遮挡的男女| 亚洲,欧美精品.| av黄色大香蕉| 久久99精品国语久久久| 午夜91福利影院| 人人妻人人爽人人添夜夜欢视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲少妇的诱惑av| 中文欧美无线码| www日本在线高清视频| 少妇高潮的动态图| 成人影院久久| 国产日韩欧美亚洲二区| 亚洲国产av新网站| 麻豆精品久久久久久蜜桃| 丝袜脚勾引网站| 欧美日韩综合久久久久久| 最近中文字幕高清免费大全6| 国产一区有黄有色的免费视频| 99热国产这里只有精品6| 国精品久久久久久国模美| 亚洲精品一二三| 色婷婷av一区二区三区视频| 亚洲成人一二三区av| 日本-黄色视频高清免费观看| av片东京热男人的天堂| 51国产日韩欧美| 亚洲国产精品一区二区三区在线| av在线播放精品| 午夜福利视频在线观看免费| 最黄视频免费看| 建设人人有责人人尽责人人享有的| 两个人免费观看高清视频| 如日韩欧美国产精品一区二区三区| 春色校园在线视频观看| 国产伦理片在线播放av一区| 亚洲国产av新网站| 热re99久久精品国产66热6| 只有这里有精品99| 久久久久精品性色| 国产成人精品无人区| 日韩成人av中文字幕在线观看| 国产av一区二区精品久久| 亚洲av福利一区| 亚洲美女搞黄在线观看| 亚洲综合色惰| 在线观看免费视频网站a站| 久久综合国产亚洲精品| 久久久久久久久久人人人人人人| 国产精品一区www在线观看| 日本与韩国留学比较| 中文字幕亚洲精品专区| 最近最新中文字幕免费大全7| 久久这里有精品视频免费| 在线观看美女被高潮喷水网站| 久久久久久久久久久免费av| 男人爽女人下面视频在线观看| 亚洲精品日韩在线中文字幕| 成年av动漫网址| 亚洲av福利一区| 尾随美女入室| 午夜老司机福利剧场| 美女大奶头黄色视频| 国产亚洲一区二区精品| 国产精品秋霞免费鲁丝片| 亚洲欧美色中文字幕在线| 丰满乱子伦码专区| 女人久久www免费人成看片| 亚洲av免费高清在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲美女搞黄在线观看| 午夜久久久在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久人人爽人人爽人人片va| 在线精品无人区一区二区三| 亚洲精品国产av成人精品| 视频在线观看一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 国产毛片在线视频| 最后的刺客免费高清国语| 交换朋友夫妻互换小说| 999精品在线视频| 18禁裸乳无遮挡动漫免费视频| 国产又色又爽无遮挡免| 美女内射精品一级片tv| 男女啪啪激烈高潮av片| 中文字幕制服av| 爱豆传媒免费全集在线观看| 亚洲性久久影院| 一区二区三区乱码不卡18| 91午夜精品亚洲一区二区三区| 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 不卡视频在线观看欧美| 老司机亚洲免费影院| 香蕉国产在线看| 一个人免费看片子| 久久人人爽av亚洲精品天堂| 欧美日韩一区二区视频在线观看视频在线| 少妇熟女欧美另类| 国产成人91sexporn| 热99久久久久精品小说推荐| 久久午夜福利片| 男女边摸边吃奶| 丝瓜视频免费看黄片| 免费av中文字幕在线| videosex国产| 国产黄色视频一区二区在线观看| 看非洲黑人一级黄片| 国产淫语在线视频| 啦啦啦视频在线资源免费观看| 亚洲欧美一区二区三区黑人 | www.av在线官网国产| 大香蕉97超碰在线| 亚洲天堂av无毛| 日本欧美视频一区| 日本av手机在线免费观看| 永久网站在线| 婷婷成人精品国产| 一级片'在线观看视频| 久久久久国产网址| 丝袜脚勾引网站| 视频区图区小说| 久久久精品免费免费高清| 亚洲精品av麻豆狂野| 中文天堂在线官网| 韩国精品一区二区三区 | 成人免费观看视频高清| 51国产日韩欧美| 亚洲伊人久久精品综合| a级片在线免费高清观看视频| 美国免费a级毛片| 中文欧美无线码| 日本欧美视频一区| 高清不卡的av网站| 男的添女的下面高潮视频| 久久人人97超碰香蕉20202| 九色成人免费人妻av| 性高湖久久久久久久久免费观看| 97精品久久久久久久久久精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美成人精品一区二区| 又黄又粗又硬又大视频| 国产精品无大码| 精品国产一区二区三区四区第35| 久久国产亚洲av麻豆专区| 成人毛片a级毛片在线播放| 在线亚洲精品国产二区图片欧美| 多毛熟女@视频| 老司机影院毛片| 性色av一级| 少妇的丰满在线观看| 久久久久久久久久成人| av女优亚洲男人天堂| 亚洲美女视频黄频| 中国美白少妇内射xxxbb| 国产在视频线精品| 波多野结衣一区麻豆| 成年动漫av网址| 18禁在线无遮挡免费观看视频| 视频在线观看一区二区三区| 日本91视频免费播放| 日本爱情动作片www.在线观看| av网站免费在线观看视频| 久久久久久人妻| 中文字幕最新亚洲高清| 天堂中文最新版在线下载| 亚洲精品国产av蜜桃| 欧美另类一区| 99精国产麻豆久久婷婷| 欧美亚洲日本最大视频资源| 国产亚洲午夜精品一区二区久久| 欧美日韩综合久久久久久| 999精品在线视频| 国产亚洲午夜精品一区二区久久| 欧美日韩综合久久久久久| 亚洲精品中文字幕在线视频| 国产激情久久老熟女| 最近中文字幕2019免费版| 国产一级毛片在线| 国产男女内射视频| 欧美日韩综合久久久久久| 黄色 视频免费看| 丰满乱子伦码专区| 狠狠婷婷综合久久久久久88av| 欧美亚洲日本最大视频资源| 免费av中文字幕在线| 肉色欧美久久久久久久蜜桃| 狠狠精品人妻久久久久久综合| av天堂久久9| 亚洲一码二码三码区别大吗| 久久精品国产鲁丝片午夜精品| 黄色怎么调成土黄色| 伊人亚洲综合成人网| 日本91视频免费播放| 久久 成人 亚洲|