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

    舟山群島海域沉積物厭氧氨氧化細(xì)菌多樣性

    2015-01-19 07:06:22張東聲劉鎮(zhèn)盛張海峰王小谷王春生
    生態(tài)學(xué)報(bào) 2015年19期
    關(guān)鍵詞:親緣文庫站位

    張東聲, 劉鎮(zhèn)盛, 張海峰, 王小谷, 王春生,2

    1 國家海洋局海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點(diǎn)實(shí)驗(yàn)室, 杭州 310012 2 衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 杭州 310012

    舟山群島海域沉積物厭氧氨氧化細(xì)菌多樣性

    張東聲1,*, 劉鎮(zhèn)盛1, 張海峰1, 王小谷1, 王春生1,2

    1 國家海洋局海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點(diǎn)實(shí)驗(yàn)室, 杭州 310012 2 衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 杭州 310012

    通過構(gòu)建16S rRNA基因文庫和克隆測序研究了舟山海域沉積物中厭氧氨氧化細(xì)菌(AAOB)的多樣性。從5個(gè)克隆文庫中共獲得297條16S rRNA基因序列,包含16個(gè)操作分類單元(OTUs)。離岸距離較近的3個(gè)站具有相似的群落結(jié)構(gòu),且與離岸較遠(yuǎn)的2個(gè)站具有明顯差異。系統(tǒng)發(fā)育結(jié)果顯示,Scalindua屬是該海域AAOB的優(yōu)勢類群,95.3%的序列與Scalindua屬AAOB具有較近的親緣關(guān)系;1條序列與Kuenenia屬具有較近的親緣關(guān)系;此外還有15條序列與已知的AAOB相似性較低。相關(guān)性分析表明沉積物有機(jī)碳含量與多樣性指數(shù)具有顯著正相關(guān),可能是該海域AAOB多樣性變化的重要影響因子。

    厭氧氨氧化細(xì)菌(AAOB); 舟山群島; 海洋沉積物; 多樣性

    微生物的硝化作用能夠?qū)h(huán)境中的銨鹽和亞硝酸鹽轉(zhuǎn)化成硝酸鹽[1],因而對保持河口等生境的水質(zhì)和生態(tài)健康具有重要作用。然而,近年來河口生態(tài)系統(tǒng)的季節(jié)性缺氧現(xiàn)象越發(fā)嚴(yán)重,其出現(xiàn)的頻率、范圍、持續(xù)時(shí)間和強(qiáng)度都呈現(xiàn)持續(xù)上升的趨勢[2]。低氧容易導(dǎo)致好氧的硝化作用中斷,使具有生物毒性的亞硝酸鹽大量積聚,進(jìn)而引發(fā)一系列生態(tài)和環(huán)境問題[3]。

    厭氧氨氧化(Anaerobic ammonium oxidation, anammox)是由一類微生物主導(dǎo),在缺氧條件下以亞硝酸鹽為底物進(jìn)行氨氧化并產(chǎn)生氮?dú)獾倪^程[4]。當(dāng)硝化作用受到阻礙時(shí),厭氧氨氧化能夠有效地降低亞硝酸鹽的濃度,有助于保持水質(zhì),維持生態(tài)系統(tǒng)的健康。厭氧氨氧化作用最初在廢水處理裝置中發(fā)現(xiàn)[5],主要受浮霉菌門的一類化能自養(yǎng)細(xì)菌調(diào)控,這類細(xì)菌很難從自然環(huán)境中分離純化[6],目前已有描述的包括5個(gè)屬(Brocadia,Kuenenia,Scalindua,Anammoxoglobus和Jettenia)[7-9]。厭氧氨氧化細(xì)菌(AAOB)在自然環(huán)境中的分布非常廣泛,在淡水沉積物[10-11]、海水沉積物[12]和低氧海水[13]都有發(fā)現(xiàn),厭氧氨氧化作用是水生環(huán)境中維持氮平衡的一個(gè)重要途徑。

    舟山群島位于長江口外海域,受海洋和長江徑流的影響,屬于咸淡水生態(tài)系統(tǒng)的交錯(cuò)群落,具有獨(dú)特的氮元素生物地球化學(xué)循環(huán)特征,是厭氧氨氧化作用的重要發(fā)生場所。近年來長江口外夏季低氧現(xiàn)象愈加嚴(yán)重,在夏季整個(gè)舟山漁場幾乎都被低氧區(qū)所覆蓋[14-15]。低氧現(xiàn)象嚴(yán)重危害該海域的漁業(yè)資源和生態(tài)系統(tǒng)健康、破壞該海域的海洋經(jīng)濟(jì)[2]。Dang等人曾對該海域沉積物中的氨氧化古菌進(jìn)行過報(bào)道[16],但目前還未有關(guān)于舟山群島海域低氧區(qū)AAOB的研究。本文采用16S rRNA基因文庫構(gòu)建和克隆測序法,研究夏季舟山海域沉積物中的AAOB多樣性,旨在初步了解其種類組成和分布情況,為揭示海洋沉積環(huán)境中AAOB對低氧的響應(yīng)機(jī)制提供科學(xué)依據(jù)。

    1 材料和方法

    1.1 樣品采集

    2012年6月在舟山以東海域采集沉積物表層泥樣和底層水樣,采樣站位信息見表1。沉積物泥樣使用抓斗采泥器采集,取少量沉積物樣品裝入預(yù)滅菌過的微生物采樣杯,-20℃冷凍保存,帶回實(shí)驗(yàn)室分析;底層水樣使用Niskin采水器采集,水樣經(jīng)GF/F濾膜過濾,于-20℃冷凍保存。

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

    水樣溫度和鹽度使用HydroLab多參數(shù)分析儀現(xiàn)場測定,水樣營養(yǎng)鹽和溶解氧濃度按照《海洋調(diào)查規(guī)范》[17]使用分光光度儀測定,沉積物有機(jī)碳、總碳和總氮含量參照《海洋監(jiān)測規(guī)范》[18]使用元素分析儀測定。

    表1 采樣站位信息

    1.3 DNA的提取和擴(kuò)增

    取少量(1—2 g)沉積物樣品,使用Fast DNA@SPIN Kit for soil(MP,美國)提取沉積物總DNA。使用AAOB的16S rRNA基因特異引物[19]進(jìn)行DNA擴(kuò)增:Amx368F,5′-TTCGCAATGCCCGAAAGG-3′和Amx820R,5′-AAAACCCCTCTACTTAGTGCCC-3′, 50 μl PCR反應(yīng)體系包含:5 μL 10×PCR buffer、200nmol/L dNTPs、前后引物各0.25 μmol/L、1 U Taq酶和1 μL DNA提取物。PCR反應(yīng)條件:94 ℃ 4 min;94 ℃ 30 s,56 ℃ 30 s,72 ℃ 60 s,30個(gè)循環(huán);72 ℃ 7 min。1%瓊脂糖凝膠電泳檢測PCR擴(kuò)增結(jié)果,并用QIAquick Gel Extraction Kit(QIAGEN,美國)純化、回收DNA目的擴(kuò)增產(chǎn)物。

    1.4 TA克隆與測序

    PCR回收產(chǎn)物與pMD20-T載體(TaKaRa)連接,轉(zhuǎn)化到感受態(tài)細(xì)胞E.coliDH5α(TaKaRa),涂布于LB平板(含Amp、IPTG和X-Gal),37℃培養(yǎng)過夜,篩選陽性克隆構(gòu)建基因克隆文庫。將所獲的陽性克隆送生工生物工程(上海)測序。

    1.5 AAOB 16S rRNA基因序列分析

    測序獲得的16S rRNA 基因序列在NCBI 數(shù)據(jù)庫中進(jìn)行BLAST 比對(http://www.ncbi.him.nih.gov),下載相似性最高的序列作為參比序列,應(yīng)用Clustal X進(jìn)行匹配比對。用DOTUR軟件包將相似性≥97%的序列定義為一個(gè)分類單元(OTU)[20],并構(gòu)建稀釋度曲線。用MEGA4 軟件構(gòu)建系統(tǒng)發(fā)育樹[21]。研究獲得序列已提交GenBank,注冊號(hào)為KF029766-KF030062。

    1.6 多樣性分析

    用DOTUR軟件分析各個(gè)克隆文庫的生物多樣性指標(biāo)(Shannon和Chao 1),根據(jù)OTUs分析結(jié)果計(jì)算各個(gè)克隆文庫的覆蓋度,計(jì)算公式如下[22]:

    Good =[1-(n/N)]×100

    式中,n代表單克隆OUT的數(shù)量,N代表文庫中克隆總數(shù)量。AAOB群落的生態(tài)分布特征使用Fast UniFrac進(jìn)行分析[23]。

    2 結(jié)果

    2.1 環(huán)境特征

    研究海域5個(gè)采樣站位根據(jù)其環(huán)境特征可分為河口和海洋兩類(表2)。EZ3-1站距離長江口最近,受長江徑流影響底層海水溫度較高、鹽度較低,硝酸鹽和硅酸鹽濃度是研究站位中最高,沉積物中碳、氮含量也是研究站位中最高。其它4個(gè)站位環(huán)境特征較為均勻,鹽度較高、營養(yǎng)鹽濃度較低,屬典型的海洋環(huán)境特征,但EZ1-3和EZ3-3兩個(gè)站營養(yǎng)鹽濃度和沉積物碳、氮含量略高于EZ1-5和EZ3-5,表明長江徑流等對距離較近的站位仍有一定的影響。研究海域底層水的溶解氧濃度較低,其中EZ3-3和EZ3-5站位于長江口外低氧區(qū)的核心區(qū)內(nèi)[14],其溶解氧低于其它3個(gè)站(<4 mg/L)。

    表2 各采樣站位環(huán)境參數(shù)

    2.2 克隆文庫與多樣性分析

    圖1 厭氧氨氧化細(xì)菌16S rRNA基因稀釋度曲線Fig.1 Rarefaction curves of AAOB 16S rRNA libraries

    針對5個(gè)站沉積物樣品分別構(gòu)建克隆文庫,總共獲得297條AAOB 16S rRNA基因序列。根據(jù)基因序列相似性≥97%歸為一個(gè)OTU,共獲得16個(gè)OTUs,5個(gè)文庫共有的OTUs有5個(gè),文庫各自的OUT數(shù)在5—8個(gè)之間(表3)。稀釋度曲線顯示5個(gè)基因文庫用于測序的克隆數(shù)量達(dá)到或接近飽和(圖1)。5個(gè)克隆文庫的覆蓋度較高,均大于80%(表3),表明本研究所構(gòu)建的文庫基本涵蓋研究海域沉積物中的AAOB多樣性。各個(gè)文庫的香濃指數(shù)和Chao 1指數(shù)之間有一定差異(表3),結(jié)果顯示EZ1-5、EZ3-1和EZ3-3文庫的多樣性較高,EZ1-3和EZ3-5文庫的多樣性較低,而EZ1-3和EZ1-5文庫的物種數(shù)量較高,EZ3-1、EZ3-3和EZ3-5文庫的物種數(shù)量較低。

    表3 厭氧氨氧化細(xì)菌基因文庫的多樣性指數(shù)

    圖2 基于16S rRNA基因的厭氧氨氧化細(xì)菌群落PCoA分析 Fig.2 The Fast UniFrac weighted PCoA analysis of the anammox communities using 16S rRNA gene sequences

    使用FastUniFrac PCoA分析了研究海域5個(gè)站位AAOB群落組成的空間分布特征。圖2顯示,PCoA的前兩個(gè)主坐標(biāo)軸(PC1和PC2)能夠解釋5個(gè)采樣站位之間AAOB群落差異的70.25%,根據(jù)AAOB的群落結(jié)構(gòu)可以把5個(gè)站分成3類,EZ1-3、EZ3-1和EZ3-3三個(gè)站位具有相似的AAOB群落結(jié)構(gòu),EZ1-5和EZ3-5兩個(gè)站AAOB的群落結(jié)構(gòu)與上述3個(gè)站具有明顯的差別。

    2.3 系統(tǒng)發(fā)育分析

    對從舟山海域沉積物樣品中構(gòu)建的5個(gè)AAOB的16S rRNA基因文庫進(jìn)行了系統(tǒng)發(fā)育分析,結(jié)果顯示,本研究獲得的絕大部份序列(283條)屬于Scalindua屬(圖3),占總序列數(shù)的95.3%。由圖3可以看出,來自Scalindua屬的序列主要聚類于Brodae和Wagneri兩個(gè)分支。Brodae分支包含4個(gè)OUT,245條序列,占Scalindua屬序列數(shù)的86.6%,與CandidatusScalinduabrodae的親緣關(guān)系較近。其中,OTU1包含47條序列,與CandidatusS.brodae的相似性在96.0%—97.7%之間,主要與來自河口沉積物[24]、深海熱液噴口沉積物[25]的AAOB非培克隆具有較近的親緣關(guān)系;OTU2包含53條序列,與CandidatusS.brodae的相似性在96.8%—97.9%之間,與來自河口沉積物中[24]的AAOB非培克隆具有較近的親緣關(guān)系;OTU3包含26條序列,與CandidatusS.brodae(AY254883)的相似性最高(>98%),與來自秘魯上升流和納米比亞上升流區(qū)的氧最小層(OMZ)中[26-27]的AAOB非培克隆具有較近的親緣關(guān)系;OTU4包含119條序列,與CandidatusS.brodae的相似性在95.1%—95.6%之間,主要與來自長江口、南海和黃海的AAOB非培克隆[28-30]具有較近的親緣關(guān)系。Wagneri分支包含4個(gè)OTU,共36條序列,占Scalindua屬序列總數(shù)的12.7%,與CandidatusScalinduawagneri的親緣關(guān)系較近。其中OTU10的10條序列與CandidatusS.wagneri的相似性較高(>97%),其他3個(gè)OTU的26條序列與CandidatusS.wagneri的相似性在94.5%—96.6%之間,與來自潮間帶沉積物、河口沉積物、和近岸低氧沉積物等環(huán)境中的AAOB非培克隆[24,31,32]具有較近的親緣關(guān)系。此外,2條來自EZ3-1文庫的序列EZ3-1-8(OTU5)和EZ3-1-49(OTU6)也聚類于Scalindua分支內(nèi),但與已知CandidatusScalinduabrodae和CandidatusScalinduawagneri的相似性均較低,分別為91%和93%左右,它們可能代表了該屬未知的AAOB菌株。

    圖3 AAOB 16S rRNA基因系統(tǒng)發(fā)育樹Fig.3 Neighbor-joining phylogenetic tree of anammox bacteria-related 16S rRNA gene fragments from the study areaClong names include the sample name and the number of times a sequence among all of the sequenced clones of samples; Bootstrap values represent 1000 replications and only values above 50% are shown; Branch lengths correspond to sequence differences as indicated by the scale bar; Numbers in parentheses refer to the number of clones were assigned to an OUT

    EZ3-5文庫中有一條序列EZ3-5-46(OTU11)聚類于Kuenenia分支,與Kuenenia屬的其代表種CandidatusKueneniastuttgartiensis的相似性為94.6%,與來自河口、陸源淡水和深海熱液噴口的AAOB[24,33]具有較近的親緣關(guān)系。此外,還有5個(gè)OTU共15條序列與已知AAOB相似性較低,這些序列可分為2個(gè)分支,Unknown I分支包含3個(gè)OTU,9條序列,與來自極地海洋沉積物中的非培克隆[34]具有較近的親緣關(guān)系;Unknown II分支包含2個(gè)OUT,6條序列,與來自紅樹林沉積物中的非培養(yǎng)細(xì)菌具有較近的親緣關(guān)系(圖3)。

    3 討論

    厭氧氨氧化細(xì)菌(AAOB)廣泛分布在河流、湖泊、海洋等不同的環(huán)境中,群落組成在不同環(huán)境中具有明顯的差異性。Scalindua屬是海洋環(huán)境中AAOB的優(yōu)勢類群[35],它們具有較高的鹽耐受性,在多種海洋環(huán)境中,特別是上升流區(qū)水柱的OMZ和一些次低氧水體中[35-36]都有報(bào)道。本文研究結(jié)果表明舟山群島海域的AAOB是典型的海洋環(huán)境群落,Scalindua屬是研究海域沉積物中AAOB的優(yōu)勢類群,來自于該屬的序列占總序列數(shù)的95.3%。本文研究海域還獲得了1條與Kuenenia屬AAOB 親緣關(guān)系較近的16S rRNA基因序列。該屬主要分布在淡水或陸源生境中[36-37],在河口等咸淡水交匯的低鹽環(huán)境中也有分布[24,38]。但本文在距離陸地較遠(yuǎn)的EZ3-5站獲得了Kuenenia屬序列,該站位主要受外海水團(tuán)控制,鹽度較高(34.32),營養(yǎng)鹽濃度較低,而在鹽度相對較低的EZ3-1站未發(fā)現(xiàn)Kuenenia屬序列。近期Byrne等在北大西洋中脊的深海熱液區(qū)的貽貝和煙囪體中也獲得了與Kuenenia屬親緣關(guān)系較近的序列[33],表明Kuenenia屬的部分AAOB對鹽度具有較強(qiáng)的耐受能力。本文和Byrne等的研究結(jié)果表明來自Kuenenia屬AAOB的分布范圍可能超出我們原來的認(rèn)識(shí),一方面這可能受到采樣和研究不足的限制,另一方面需要更深水平的二代測序研究證實(shí)。此外,除了已知的AAOB類群,本研究還發(fā)現(xiàn)了一些與來自海洋環(huán)境的浮霉菌門非培養(yǎng)克隆相似性較低(<85%)的序列(Unkown I和Unknown II)。有研究表明,海洋環(huán)境中可能存在除浮霉菌門以外的其它厭氧氨氧化微生物,如β-和γ-變形菌中的Nitrosospira、Nitrosomonas和Nitrosococcus[39-40],甚至在某些海洋環(huán)境中存在具有厭氧氨氧化能力的古菌[41]。因此,這些序列可能代表了某些未知門類的厭氧氨氧化微生物。

    圖4 有機(jī)碳與多樣性指數(shù)相關(guān)性 Fig.4 Correlation between organic carbon and Shannon index of the five sampling sites

    AAOB多樣性和群落結(jié)構(gòu)的分布特征與環(huán)境因子有著密不可分的關(guān)系。Dale等認(rèn)為鹽度是影響河口生態(tài)系統(tǒng)中AAOB群落結(jié)構(gòu)的重要因子[31],但本文采樣站位之間鹽度差別不大,除EZ3-1略低外,其它站位的鹽度十分接近。本文中EZ3-1、EZ3-3和EZ1-3三個(gè)站的群落結(jié)構(gòu)較為接近,而與EZ1-5和EZ3-5站不同(圖2),表明離岸距離和水深可能是影響表層沉積物AAOB群落組成的重要因素[30]。另一方面,各站位的多樣性指數(shù)卻呈現(xiàn)出不同的分布特征,EZ3-1、EZ3-3和EZ1-5三個(gè)站的多樣性指數(shù)高于EZ1-3和EZ3-5站(表2)。相關(guān)分析顯示沉積物有機(jī)碳含量與多樣性指數(shù)具有顯著正相關(guān)關(guān)系(圖4)。Hou等研究認(rèn)為沉積物有機(jī)碳含量是河口生態(tài)系統(tǒng)中影響AAOB多樣性的重要環(huán)境因子[29],高有機(jī)碳含量能夠促進(jìn)硝酸鹽還原作用,進(jìn)而提高亞硝酸鹽濃度,亞硝酸鹽是厭氧氨氧化過程中的電子受體,高濃度的亞硝酸鹽有利于促進(jìn)厭氧條件下的氨氧化作用[42]。Hu等人在椒江口海域的研究表明AAOB多樣性指數(shù)與沉積物有機(jī)碳含量呈正比[24],Li等人在南海北部也有類似結(jié)果[30]。

    [1] Ward B B. Nitrification and the marine nitrogen cycle // Kirchman D L. Microbial Ecology of the Oceans New York: John Wiley and Sons, 2000: 427-453.

    [2] 陳春輝, 王春生, 許學(xué)偉, 劉鎮(zhèn)盛. 河口缺氧生物效應(yīng)研究進(jìn)展. 生態(tài)學(xué)報(bào), 2009, 29(5): 2595-2602.

    [3] Alonso A, Camargo J A. Toxicity of nitrite to three species of freshwater invertebrates. Environmental Toxicology, 2006, 21(1): 90-94.

    [4] Jetten M S M, Logemann S, Muyzer G, Robertson L A, Vries S, van Loosdrecht M C M, Kuenen J G. Novel principles in the microbial conversion of nitrogen compounds. Antonie Van Leeuwenhoek, 1997, 71(1/2): 75-93.

    [5] Mulder A, van de Graff A A, Robertson L A, Kuenen J G. Anaerobic ammonium oxidation discovered in a denitrifying fluidized bed reactor. FEMS Microbiology Ecology, 1995, 16(3): 177-183.

    [6] Strous M, Jetten M S M. Anaerobic oxidation of methane and ammonium. Annual Review of Microbiology, 2004, 58: 99-117.

    [7] Kartal B, Rattray J, van Niftrik LA, van de vossenberg J, Schmid M C, Webb R I, Schouten S, Fuerst J A, Damste J S, Jetten M S M, Strous M.Candidatus"Anammoxglobuspropionicus" a new propionate oxidizing species of anaerobic ammonium oxidizing bacteria. Systematic and Applied Microbiology, 2007, 30(1): 39-49.

    [8] Schmid M, Walsh K, Webb R, Rijpstra W I, van de Pas-Schoonen K, Verbruggen M J, Hill T, Moffett B, Fuerst J, Schouten S, Damste J S, Harris J, Shaw P, Jetten M S M, Strous M.Candidatus"ScalinduaBrodae", sp. nov.,Candidatus"ScalinduaWagneri", sp. nov., two new species of anaerobic ammonium oxidizing bacteria. Systematic and Applied Microbiology, 2003, 26(4): 529-538.

    [9] Schmid M C, Maas B, Dapena A, van de Pas-Schoonen K, van de Vossenberg J, Kartal B, van Niftrik L, Schmidt I, Cirpus I, Kuenen J G, Wagner M, Damste J S, Kuypers M, Revsbech N P, Mendez R, Jetten M S M, Strous M. Biomarkers for the in situ detection of anaerobic ammonium oxidizing (anammox) bacteria. Applied and Environmental Microbiology, 2005, 71(4): 1677-1684.

    [10] Hu B L, Shen L D, Zheng P, Hu A H, Chen T T, Cai C, Liu S, Lou L P. Distribution and diversity of anaerobic ammonium-oxidizing bacteria in the sediments of the Qiantang River. Environmental Microbiology Reports, 2012, 4(5): 540-547.

    [11] Zhu GB, Wang S Y, Wang W D, Wang Y, Zhou L L, Jiang B, Huub J M, Camp O, Risgaard-Petersen N, Schwark L, Peng Y Z, Hefting M M, Jetten M S M, Yin C Q. Hotspots of anaerobic ammonium oxidation at land freshwater interfaces. Nature Geoscience, 2013, 6(2): 103-107.

    [12] Dang H Y, Chen R P, Wang L, Guo L Z, Chen P P, Tang Z W. Environmental factors shape sediment anammox bacterial communities in hyper nitrified Jiaozhou Bay, China. Applied and Environmental Microbiology, 2010, 76(21): 7036-7047.

    [13] Lam P, Jensen M M, Lavik G, McGinnis D F, Müller B, Schubert C J, Amann R, Thamdrup B, Kuypers M M M. Linking crenarchaeal and bacterial nitrification to anammox in the Black Sea, Proceedings of the National Academy of the Sciences of the United States of America, 2007, 104(17): 7104-7109.

    [14] Li D J, Zhang J, Huang D J, Wu Y, Liang J. Oxygen depletion off the Changjiang (Yangtze River) Estuary. Science in China (Series D), 2002, 45(12): 1137-1146.

    [15] Chen C C, Gong G C, Shiah F K. Hypoxia in the East China Sea: One of the largest coastal low-oxygen areas in the world. Marine Environmental Research, 2007, 64(4): 399-408.

    [16] Dang H Y, Zhang X X, Sun J, Li T G, Zhang Z N, Yang G P. Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea. Microbiology, 2008, 154(7): 2084-2095.

    [17] GB/T 12763. 4-2007, 海洋調(diào)查規(guī)范第4部分: 海水化學(xué)要素調(diào)查. 北京: 中國標(biāo)準(zhǔn)出版社, 2008.

    [18] GB/T 17378. 5-2007, 海洋監(jiān)測規(guī)范第5部分: 沉積物分析. 北京: 中國標(biāo)準(zhǔn)出版社, 2008.

    [19] Amano T, Yoshinaga I, Okada K, Yamagishi T, Ueda S, Obuchi A, Sako Y, Suwa Y. Detection of anammox activity and diversity of anammox bacteria-related 16S rRNA genes in coastal marine sediment in Japan. Microbes and Environment, 2007, 22(3): 232-242.

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

    [21] Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4. 0. Molecular Biology and Evolution, 2007, 24(8): 1596-1599.

    [22] Good I J. The population frequencies of species and the estimation of population parameters. Biometrika, 1953, 40(3/4): 237-264.

    [23] Sagaram U S, De Angelis K M, Trivedi P, Andersen G L, Lu S E, Wang N. Bacterial diversity analysis of Huanglongbing pathogen-infected citrus, using PhyloChip arrays and 16S rRNA gene clone library sequencing. Applied and Environmental Microbiology, 2009, 75(6): 1566-1574.

    [24] Hu B L, Shen L D, Du P, Zheng P, Xu X Y, Zeng J N. The influence of intense chemical pollution on the community composition, diversity and abundance of anammox bacteria in the Jiaojiang Estuary (China). PLoS ONE, 2012, 7(3): e33826.

    [25] Hirsch M D, Long Z T, Song B. Anammox bacterial diversity in various aquatic ecosystems based on the detection of hydrazine oxidase genes (hzoA/hzoB). Microbial Ecology, 2011, 61(2): 264-276.

    [26] Woebken D, Fuchs BM, Kuypers MMM, Amann R. Potential interactions of particles-associated anammox bacteria with bacterial and archaeal partners in the Nambian upwelling system. Applied and Environmental Microbiology, 2007, 73(14): 4648-4657.

    [27] Woebken D, Lam P, Kuypers M M M, Naqvi S W A, Kartal B, Strous M, Jetten M S M, Fuchs B M, Amann R. A microdiversity study of anammox bacteria reveals a novelCandidatusScalindua phylotype in marine oxygen minimum zones. Environmental Microbiology, 2008, 10(11): 3106-3119.

    [28] Hong J K and Cho J C. High level of bacterial diversity and novel taxa in continental shelf sediment. Journal of Microbiology and Biotechnology, 2012, 22(6): 771-779.

    [29] Hou L J, Zheng Y L, Liu M, Gong J, Zhang X L, Yin G Y, You L. Anaerobic ammonium oxidation (anammox) bacterial diversity, abundance, and activity in marsh sediments of the Yangtze Estuary. Journal of Geophysical Research: Biogeosciences, 2013, 118(3): 1237-1246.

    [30] Li M, Hong Y G, Cao H L, Gu J D. Community structures and distribution of anaerobic ammonium oxidizing andnirS-encoding nitrite-reducing bacteria in surface sediments of the South China Sea. Microbial Ecology, 2013, 66(2): 281-296.

    [31] Dale O R, Tobias C R, Song B. Biogeographic distribution of diverse anaerobic ammonium oxidizing (anammox) bacteria in Cape Fear River Estuary. Environmental Microbiology, 2009, 11(5): 1194-1207.

    [32] Prokepenko M G, Hirst M B, De Brabandere L, Lawrence D J P, Berelson W M, Granger J, Chang B X, Dawson S, Crane III E J, Chong L, Thamdrup B, Townsend-Small A, Sigman D M. Nitrogen losses in anoxic marine sediments driven byThioploca-anammox bacterial consortia. Nature, 2013, 500(7461): 194-198.

    [33] Byrne N, Strous M, Crepeau V, Kartal B, Birrien J L, Schmid M, Lesongeur F, Schouten S, Jaeschke A, Jetten M, Prieur D, Godfroy A. Presence and activity of anaerobic ammonium-oxidizing bacteria at deep-sea hydrothermal vents. International Symposium on Microbial Ecology Journal, 2009, 3(1): 117-123.

    [34] Tait K, Laverock B, Widdicombe S. Response of an Arctic sediment nitrogen cycling community to increased CO2. Estuaries and Coasts, 2014, 37(3): 724-735.

    [35] Schmid M C, Risgaard-Petersen N, van de Vossenberg J, Kuypers M M, Lavik G, Petersen J, Hulth S, Thamdrup B, Canfield D, Dalsgaard T, Rysgarrd S, Sejr M K, Strous M, den Camp H J, Jetten M S. Anaerobic ammonium-oxidizing bacteria in marine environments: widespread occurrence but low diversity. Environmental Microbiology, 2007, 9(6): 1476-1484.

    [36] Schubert C J, Durisch-Kaiser E, Wehrli B, Thamdrup B, Lam P, Kuypers M M M. Anaerobic ammonium oxidation in a tropical freshwater system (LakeTanganyika). Environmental Microbiology, 2006, 8(10): 1857-1863.

    [37] Penton C R, Devol A H, Tiedje J M. Molecular evidence for the broad distribution of anaerobic ammonium-oxidizing bacteria in freshwater and marine sediments. Applied and Environmental Microbiology, 2006, 72(10): 6829-6832.

    [38] Jetten M S M, Sliekers A O, Kuypers M, Dalsgaard T, van Niftrik L, Cirpus I, van de Pas-Schoonen K, Lavik G, Thamdrup B, Le Paslier D, Op den Camp H J, Hulth S, Nielsen L P, Abma W, Third K, Engstrom P, Kuenen J G, Jorgensen B B, canfield D E, Sinninghe-Damste J S, Revsbech N P, Fuerst J, Weissenbach J, Wagner M, Schmidt I, Schmid M, Strous M. Anaerobic ammonium oxidation by marine and freshwater Planctomycete-like bacteria. Applied Microbiology and Biotechnology, 2003, 63(2): 107-114.

    [39] Freitag T E, Prosser J I. Community structure of ammonia-oxidizing bacteria within anoxic marine sediments. Applied and Environmental Microbiology, 2003, 69(3): 1359-1371.

    [40] O′Mullan G D, Ward B B. Relationship of temporal and spatial variability of ammonia-oxidizing bacteria to nitrification rates in Monterey Bay, California. Applied and Environmental Microbiology, 2005, 71(2): 697-705.

    [41] Wuchter C, Abbas B, Coolen M J, Herfort L, vanBleijswijk J, Timmers P, Strous M, Teira E, Herndl G J, Middelburg J J, Schouten S, Sinninghe Damste J S. Archaeal nitrification in the ocean. The Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(3): 12317-12322.

    [42] Lam P, Lavik G, Jensen M M, van de Vossenberg J, Schmid M C, Woebken D, Gutierrez D, Amann R, Jetten M S M, Kuypers M M M. Revising the nitrogen cycle in the Peruvian oxygen minimum zone. The Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(12): 4752-4757.

    Diversity of anaerobic ammonium oxidizing bacteria in marine sediments from the Zhoushan Islands

    ZHANG Dongsheng1,*, LIU Zhensheng1, ZHANG Haifeng1, WANG Xiaogu1, WANG Chunsheng1,2

    1LaboratoryofMarineEcosystemandBiogeochemistry,StateOceanicAdministration,Hangzhou310012,China2StateKeyLaboratoryofSatelliteOceanEnvironmentDynamics,Hangzhou310012,China

    Anaerobic ammonium oxidation (anammox) is an important process regulating the balance of marine nitrogen and ecosystem health, particularly under anoxic conditions. The Zhoushan Islands are located east of the Changjiang river estuary, and collect a high load of anthropogenic nitrogen, which leads to severe eutrophication and seasonal hypoxia. Therefore, bacteria that mediate the anammox process are of major interest in this area. Although the importance of anammox-mediating bacteria is known, few studies on these bacteria have been conducted in the East China Sea. To the best of our knowledge, this study is the first to report the diversity, community composition, and distribution of anammox bacteria in the Zhoushan Islands. Field surveys were conducted in June 2012; triplicate surface sediment samples were collected at each site and stored in sterile plastic bags at -80℃ for subsequent DNA extraction and molecular analysis. Total genomic DNA was extracted using the Fast DNA? SPIN Kit for soil. Environmental DNA extracted from sediment samples was used as the template for PCR amplification of anammox 16S rRNA genes using primers Amx368f—Amx820r. The purified fragments were cloned and sequenced for phylogenetic and statistical analyses. In total, 297 sequences belonging to 16 operational taxonomic units (OTUs) were obtained from five 16S rRNA gene libraries. The biodiversity of anammox bacteria was examined using rarefaction analysis of the 16S rRNA genes, the Chao1 estimator, and Shannon index calculations. EZ3-1, EZ3-3, and EZ1-5 exhibited higher diversity than EZ1-3 and EZ3-5. A significant positive correlation between Shannon index and organic carbon content indicate that sediment organic carbon content plays an important role in modulating anammox bacterial diversity in the Zhoushan Island area. Weighted UniFrac PCoA analysis of the 16S rRNA genes demonstrated spatial heterogeneity in the community composition of anammox bacteria; the anammox bacteria in the study area could be divided into three distinct groups. EZ3-1, EZ3-3, and EZ1-3 exhibited similar community composition, while EZ1-5 and EZ3-5 clustered separately. The composition might be affected by distance from land and water depth. Phylogenetic analysis indicated that anammox bacterial communities were dominated by the genusScalindua(283 of 297 sequences). TheScalinduacluster comprised of 245Scalinduasequences with 95.4%—98.5% sequence similarity toCandidatusScalinduabrodae, and 36Scalinduasequences with 94.5%—97.6% sequence identity to the 16S rRNA gene ofCandidatusScalinduawagneri. In addition, twoScalinduasequences that grouped in theScalinduacluster were distantly related to knownScalinduaspecies, indicating that they might represent unidentified species ofScalindua. One sequence recovered from the EZ3-5 library was closely related to genusKuenenia, which is traditionally considered an anammox bacterial genus present in freshwater or low-salinity environments. Our results suggest that members of the genusKueneniamay possess the ability to survive in high-salinity marine environments. Additionally, two clusters of unknown sequences (unknown cluster I and II) were not classified under any known anammox bacterial genus, but were most closely related to 16S rRNA gene sequences recovered from arctic sediment and mangrove sediment, respectively.

    anaerobic ammonium oxidation bacteria (AAOB); Zhoushan islands; marine sediments; diversity

    浙江省自然科學(xué)基金(Y5110185, Y5110171); 國家自然科學(xué)基金青年基金項(xiàng)目(41206104, 41203085)

    2014-02-21; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2014-12-04

    10.5846/stxb201402210303

    *通訊作者Corresponding author.E-mail: dszhang7@163.com

    張東聲, 劉鎮(zhèn)盛, 張海峰, 王小谷, 王春生.舟山群島海域沉積物厭氧氨氧化細(xì)菌多樣性.生態(tài)學(xué)報(bào),2015,35(19):6250-6258.

    Zhang D S, Liu Z S, Zhang H F, Wang X G, Wang C S.Diversity of anaerobic ammonium oxidizing bacteria in marine sediments from the Zhoushan Islands.Acta Ecologica Sinica,2015,35(19):6250-6258.

    猜你喜歡
    親緣文庫站位
    谷子近緣野生種的親緣關(guān)系及其利用研究
    提高政治站位 對標(biāo)國內(nèi)一流
    中國醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    建黨百年說“站位”
    水上消防(2021年3期)2021-08-21 03:12:00
    專家文庫
    優(yōu)秀傳統(tǒng)文化啟蒙文庫
    幽默大師(2020年10期)2020-11-10 09:07:22
    關(guān)于推薦《當(dāng)代詩壇百家文庫》入選詩家的啟事
    中華詩詞(2019年1期)2019-11-14 23:33:56
    提升站位講政治 創(chuàng)新?lián)?dāng)爭出彩
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    專家文庫
    男女之事视频高清在线观看| 亚洲欧洲精品一区二区精品久久久| 免费在线观看日本一区| 亚洲成人免费电影在线观看| 美女免费视频网站| 叶爱在线成人免费视频播放| 一区二区三区精品91| 99热这里只有精品一区 | 变态另类丝袜制服| 香蕉国产在线看| 国产亚洲欧美精品永久| 日本熟妇午夜| av在线播放免费不卡| av在线天堂中文字幕| av福利片在线| 精品国内亚洲2022精品成人| 一区二区三区激情视频| 99精品欧美一区二区三区四区| 国产精品 国内视频| 真人一进一出gif抽搐免费| 久久伊人香网站| 欧美精品啪啪一区二区三区| 亚洲美女黄片视频| 久久久久久大精品| 三级毛片av免费| 精品第一国产精品| 一区二区三区激情视频| 哪里可以看免费的av片| 可以在线观看的亚洲视频| 两个人看的免费小视频| 不卡一级毛片| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 国产高清videossex| 女性被躁到高潮视频| 亚洲真实伦在线观看| 成人免费观看视频高清| 欧美日韩福利视频一区二区| 九色国产91popny在线| 女同久久另类99精品国产91| 国产高清激情床上av| 日日干狠狠操夜夜爽| 国产亚洲欧美98| 精品卡一卡二卡四卡免费| 欧美一区二区精品小视频在线| 男女床上黄色一级片免费看| www.www免费av| 亚洲国产精品久久男人天堂| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美日韩无卡精品| 久久久精品欧美日韩精品| 99久久国产精品久久久| 久久国产亚洲av麻豆专区| 999精品在线视频| 国产亚洲精品第一综合不卡| 国产黄a三级三级三级人| 久久久国产成人免费| 夜夜躁狠狠躁天天躁| 曰老女人黄片| 国产熟女xx| 男人操女人黄网站| 首页视频小说图片口味搜索| 波多野结衣巨乳人妻| 一二三四在线观看免费中文在| 12—13女人毛片做爰片一| 999精品在线视频| 国产av在哪里看| 亚洲无线在线观看| 午夜激情福利司机影院| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 国产亚洲精品久久久久5区| 成年免费大片在线观看| 老熟妇乱子伦视频在线观看| 亚洲国产看品久久| 啪啪无遮挡十八禁网站| www.精华液| 欧美不卡视频在线免费观看 | АⅤ资源中文在线天堂| 无人区码免费观看不卡| av视频在线观看入口| 可以免费在线观看a视频的电影网站| 久久久久免费精品人妻一区二区 | 久久性视频一级片| 午夜免费鲁丝| 午夜视频精品福利| 亚洲狠狠婷婷综合久久图片| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 国产在线观看jvid| 亚洲五月天丁香| 日本三级黄在线观看| 亚洲av五月六月丁香网| 欧美激情 高清一区二区三区| a在线观看视频网站| 中文在线观看免费www的网站 | 亚洲成av片中文字幕在线观看| 成年人黄色毛片网站| 高潮久久久久久久久久久不卡| 久久久久国产精品人妻aⅴ院| 亚洲精品国产区一区二| 精品欧美国产一区二区三| 国产色视频综合| 日韩大码丰满熟妇| 午夜福利视频1000在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品爽爽va在线观看网站 | 老汉色av国产亚洲站长工具| 免费女性裸体啪啪无遮挡网站| 精品久久久久久久末码| 欧美人与性动交α欧美精品济南到| 两人在一起打扑克的视频| 老司机午夜福利在线观看视频| 午夜久久久久精精品| 黄色毛片三级朝国网站| 超碰成人久久| 亚洲精品国产区一区二| 精华霜和精华液先用哪个| 99在线人妻在线中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 长腿黑丝高跟| 日本精品一区二区三区蜜桃| 一进一出好大好爽视频| 亚洲久久久国产精品| 这个男人来自地球电影免费观看| 国产日本99.免费观看| 国产v大片淫在线免费观看| 热99re8久久精品国产| 桃色一区二区三区在线观看| 女警被强在线播放| 黑人欧美特级aaaaaa片| 日韩欧美 国产精品| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 日本免费一区二区三区高清不卡| 色综合站精品国产| 美女扒开内裤让男人捅视频| 欧美成人午夜精品| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 欧美日韩精品网址| 99在线人妻在线中文字幕| 99久久精品国产亚洲精品| 一本大道久久a久久精品| 一区二区日韩欧美中文字幕| 久久性视频一级片| 亚洲va日本ⅴa欧美va伊人久久| 草草在线视频免费看| 午夜福利视频1000在线观看| 亚洲av片天天在线观看| 少妇被粗大的猛进出69影院| 麻豆国产av国片精品| netflix在线观看网站| 亚洲第一av免费看| 国产熟女午夜一区二区三区| 91在线观看av| 亚洲第一青青草原| www国产在线视频色| 黄色成人免费大全| 国产欧美日韩精品亚洲av| 国产日本99.免费观看| 欧美日韩精品网址| 日韩视频一区二区在线观看| 国产精品 欧美亚洲| 久久婷婷成人综合色麻豆| 国产精品av久久久久免费| 国产成人啪精品午夜网站| 国产黄片美女视频| 中文字幕人成人乱码亚洲影| 中国美女看黄片| 国产高清videossex| 日本成人三级电影网站| 欧美激情极品国产一区二区三区| 亚洲第一欧美日韩一区二区三区| 成人国产综合亚洲| 亚洲,欧美精品.| 给我免费播放毛片高清在线观看| 国产精品野战在线观看| 久久狼人影院| 热99re8久久精品国产| 香蕉久久夜色| 黄色女人牲交| 欧美日韩福利视频一区二区| 欧美日韩精品网址| 亚洲精品国产精品久久久不卡| 欧美日韩一级在线毛片| 在线观看舔阴道视频| 成人三级黄色视频| 在线观看午夜福利视频| 亚洲av片天天在线观看| 亚洲精品久久成人aⅴ小说| 麻豆国产av国片精品| 精品卡一卡二卡四卡免费| 日韩成人在线观看一区二区三区| 色播亚洲综合网| 91大片在线观看| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸| 97超级碰碰碰精品色视频在线观看| 国产免费男女视频| 国产精品久久久人人做人人爽| 久久久久国内视频| 脱女人内裤的视频| 精品国产一区二区三区四区第35| 久久午夜亚洲精品久久| 一本大道久久a久久精品| 久热这里只有精品99| 久久九九热精品免费| 一级毛片精品| 亚洲精品久久成人aⅴ小说| 久久人妻福利社区极品人妻图片| 欧美日韩亚洲国产一区二区在线观看| 老汉色∧v一级毛片| 琪琪午夜伦伦电影理论片6080| 夜夜夜夜夜久久久久| 欧美激情久久久久久爽电影| 午夜福利一区二区在线看| 99久久国产精品久久久| 午夜激情av网站| 成人国产综合亚洲| 欧美+亚洲+日韩+国产| 黑丝袜美女国产一区| 国产高清videossex| 国产99白浆流出| 亚洲av电影在线进入| 久久伊人香网站| 成人av一区二区三区在线看| 免费女性裸体啪啪无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| 色哟哟哟哟哟哟| 成人亚洲精品一区在线观看| 久久久久国产一级毛片高清牌| 色哟哟哟哟哟哟| 91av网站免费观看| 99热这里只有精品一区 | 一进一出抽搐动态| 熟妇人妻久久中文字幕3abv| 国产蜜桃级精品一区二区三区| 给我免费播放毛片高清在线观看| 中文字幕av电影在线播放| 欧美日本亚洲视频在线播放| 亚洲黑人精品在线| 久久久久免费精品人妻一区二区 | www.www免费av| 午夜免费鲁丝| www日本黄色视频网| 亚洲人成电影免费在线| 熟妇人妻久久中文字幕3abv| 国产av一区二区精品久久| 1024手机看黄色片| www日本在线高清视频| 国产不卡一卡二| 禁无遮挡网站| 欧美 亚洲 国产 日韩一| 久久99热这里只有精品18| 久久午夜亚洲精品久久| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 精品欧美一区二区三区在线| 一边摸一边做爽爽视频免费| 欧美性猛交╳xxx乱大交人| 国产激情久久老熟女| 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 少妇的丰满在线观看| 女同久久另类99精品国产91| 嫩草影院精品99| 欧美中文综合在线视频| 日韩欧美国产一区二区入口| 欧美激情高清一区二区三区| 99精品在免费线老司机午夜| www日本黄色视频网| 精品欧美国产一区二区三| 亚洲专区字幕在线| av天堂在线播放| 色播亚洲综合网| 人人澡人人妻人| 久久久国产欧美日韩av| 久久性视频一级片| 亚洲午夜精品一区,二区,三区| 免费在线观看完整版高清| 国产精品 国内视频| 熟女少妇亚洲综合色aaa.| 久久国产精品人妻蜜桃| 美国免费a级毛片| 午夜视频精品福利| 午夜福利欧美成人| 亚洲无线在线观看| 人人澡人人妻人| 天天躁狠狠躁夜夜躁狠狠躁| 高清毛片免费观看视频网站| 黄色a级毛片大全视频| 色在线成人网| 香蕉国产在线看| 麻豆成人av在线观看| 久久久久久久午夜电影| 免费看日本二区| 久久 成人 亚洲| 777久久人妻少妇嫩草av网站| 97碰自拍视频| 免费在线观看影片大全网站| 老汉色av国产亚洲站长工具| 免费人成视频x8x8入口观看| 两人在一起打扑克的视频| 熟女电影av网| 亚洲精品中文字幕一二三四区| 美女扒开内裤让男人捅视频| 亚洲 国产 在线| 一本一本综合久久| 亚洲av五月六月丁香网| 嫩草影院精品99| 欧美久久黑人一区二区| 丝袜在线中文字幕| 丁香欧美五月| av有码第一页| 男人操女人黄网站| 中文字幕最新亚洲高清| 国产三级在线视频| 视频在线观看一区二区三区| 一区二区三区国产精品乱码| 国产私拍福利视频在线观看| 国产亚洲欧美精品永久| 国产精品影院久久| 久久性视频一级片| 亚洲一区二区三区色噜噜| 黑丝袜美女国产一区| 真人做人爱边吃奶动态| 亚洲精品美女久久久久99蜜臀| 人妻久久中文字幕网| 一区二区三区高清视频在线| 嫁个100分男人电影在线观看| 国产亚洲精品一区二区www| 黄片小视频在线播放| 亚洲免费av在线视频| 亚洲中文字幕日韩| 极品教师在线免费播放| 色播在线永久视频| 日本精品一区二区三区蜜桃| svipshipincom国产片| 亚洲成人精品中文字幕电影| 操出白浆在线播放| 日本黄色视频三级网站网址| 国产精品 国内视频| 日本一本二区三区精品| 亚洲av片天天在线观看| 欧美激情久久久久久爽电影| 亚洲美女黄片视频| 国产av一区二区精品久久| av视频在线观看入口| 人妻久久中文字幕网| 国产高清有码在线观看视频 | 日韩欧美在线二视频| 欧美+亚洲+日韩+国产| 日韩大码丰满熟妇| 国产精品98久久久久久宅男小说| 日本免费a在线| 日日夜夜操网爽| 精品乱码久久久久久99久播| 18禁裸乳无遮挡免费网站照片 | www国产在线视频色| 国产亚洲精品一区二区www| 国产午夜福利久久久久久| 国产真人三级小视频在线观看| 久久精品影院6| 在线观看免费日韩欧美大片| 日本三级黄在线观看| 日本黄色视频三级网站网址| 国产黄a三级三级三级人| 国内精品久久久久久久电影| 国产伦在线观看视频一区| 亚洲国产看品久久| 99国产精品一区二区蜜桃av| 欧美色欧美亚洲另类二区| 人妻久久中文字幕网| 国产爱豆传媒在线观看 | 国产成人系列免费观看| 亚洲一码二码三码区别大吗| 夜夜夜夜夜久久久久| 岛国在线观看网站| avwww免费| а√天堂www在线а√下载| 成人18禁在线播放| 欧美三级亚洲精品| 少妇被粗大的猛进出69影院| 俺也久久电影网| 中文字幕久久专区| 黄色女人牲交| 亚洲国产精品久久男人天堂| 国产又爽黄色视频| 可以免费在线观看a视频的电影网站| 伊人久久大香线蕉亚洲五| 亚洲精品久久国产高清桃花| 久久人妻av系列| 麻豆成人av在线观看| 国产黄a三级三级三级人| 无人区码免费观看不卡| 欧美一级a爱片免费观看看 | e午夜精品久久久久久久| 欧美日本视频| 两个人免费观看高清视频| 成年免费大片在线观看| 日本五十路高清| 亚洲成人久久爱视频| 国产午夜福利久久久久久| 免费搜索国产男女视频| 国产黄片美女视频| 两个人视频免费观看高清| 国产伦一二天堂av在线观看| 成在线人永久免费视频| 成人欧美大片| 亚洲久久久国产精品| 91大片在线观看| 日本 av在线| 日韩大码丰满熟妇| 精品不卡国产一区二区三区| 麻豆av在线久日| 757午夜福利合集在线观看| 免费观看人在逋| 成人精品一区二区免费| 国产黄色小视频在线观看| 91在线观看av| 午夜福利免费观看在线| 国产精品野战在线观看| 女生性感内裤真人,穿戴方法视频| 人人妻人人看人人澡| 黄色视频不卡| 脱女人内裤的视频| 中文字幕人妻丝袜一区二区| 午夜久久久久精精品| 亚洲成av片中文字幕在线观看| 观看免费一级毛片| 欧美激情 高清一区二区三区| 国产成人啪精品午夜网站| 亚洲成人国产一区在线观看| 国产精品永久免费网站| 久久精品国产亚洲av香蕉五月| 欧美一级毛片孕妇| 18禁黄网站禁片免费观看直播| 久久伊人香网站| 国内揄拍国产精品人妻在线 | xxx96com| 宅男免费午夜| 91大片在线观看| 哪里可以看免费的av片| 国产1区2区3区精品| 国产av不卡久久| 精品久久久久久久久久免费视频| 国产不卡一卡二| 久久精品成人免费网站| 亚洲国产高清在线一区二区三 | 亚洲在线自拍视频| 国产v大片淫在线免费观看| 午夜日韩欧美国产| 免费高清视频大片| 高清在线国产一区| 精品国产美女av久久久久小说| 观看免费一级毛片| 国产精品,欧美在线| 最近最新中文字幕大全电影3 | 国产成人av教育| 一区二区三区高清视频在线| 青草久久国产| 亚洲人成网站在线播放欧美日韩| 久久久久九九精品影院| 日本黄色视频三级网站网址| 欧美成人一区二区免费高清观看 | 在线av久久热| 丁香欧美五月| 人妻久久中文字幕网| www日本黄色视频网| 天堂影院成人在线观看| 午夜免费观看网址| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇中文字幕五十中出| 一本精品99久久精品77| 亚洲熟妇熟女久久| 看片在线看免费视频| 美女高潮到喷水免费观看| 丰满人妻熟妇乱又伦精品不卡| 在线观看一区二区三区| 久久伊人香网站| 香蕉av资源在线| 日日摸夜夜添夜夜添小说| 国产精品爽爽va在线观看网站 | 国产高清激情床上av| 成熟少妇高潮喷水视频| 日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 国产一区二区激情短视频| 我的亚洲天堂| 99热只有精品国产| 91麻豆av在线| 亚洲国产毛片av蜜桃av| 中文字幕人妻熟女乱码| 婷婷精品国产亚洲av| 亚洲狠狠婷婷综合久久图片| 老司机福利观看| 亚洲成人久久性| 欧美乱妇无乱码| 亚洲av片天天在线观看| 国产av一区二区精品久久| 午夜激情福利司机影院| 国产精品1区2区在线观看.| 91成年电影在线观看| 日本一本二区三区精品| 在线视频色国产色| 日日干狠狠操夜夜爽| 亚洲欧洲精品一区二区精品久久久| 亚洲第一av免费看| 真人一进一出gif抽搐免费| 757午夜福利合集在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久久久久久久 | 亚洲成人久久爱视频| 国产又爽黄色视频| 日本三级黄在线观看| av欧美777| 国产99久久九九免费精品| 十八禁网站免费在线| 黑人操中国人逼视频| 欧美日韩乱码在线| 精品国产超薄肉色丝袜足j| 99在线视频只有这里精品首页| 神马国产精品三级电影在线观看 | 欧美 亚洲 国产 日韩一| 国产精品久久久久久人妻精品电影| 久久久久久免费高清国产稀缺| 国产麻豆成人av免费视频| 日日干狠狠操夜夜爽| 欧美日韩精品网址| 精品欧美国产一区二区三| 国产在线观看jvid| 亚洲天堂国产精品一区在线| av在线天堂中文字幕| 亚洲av五月六月丁香网| 久久性视频一级片| 国产精品影院久久| 欧美大码av| 国产精品日韩av在线免费观看| 美女国产高潮福利片在线看| 91国产中文字幕| 少妇粗大呻吟视频| 法律面前人人平等表现在哪些方面| 久久久久久久久久黄片| 超碰成人久久| 嫁个100分男人电影在线观看| 国产野战对白在线观看| 丰满人妻熟妇乱又伦精品不卡| 露出奶头的视频| 成人国产综合亚洲| 女警被强在线播放| 久久狼人影院| 一级毛片高清免费大全| 午夜福利在线观看吧| 日韩免费av在线播放| 日韩一卡2卡3卡4卡2021年| 日韩欧美一区视频在线观看| 男男h啪啪无遮挡| 69av精品久久久久久| 黄网站色视频无遮挡免费观看| 性欧美人与动物交配| 婷婷六月久久综合丁香| 99久久99久久久精品蜜桃| 日日干狠狠操夜夜爽| 国产精品综合久久久久久久免费| 精品欧美一区二区三区在线| 欧美黑人精品巨大| 色播在线永久视频| 国产成人精品无人区| 999久久久国产精品视频| 首页视频小说图片口味搜索| 天堂动漫精品| 香蕉久久夜色| 亚洲五月天丁香| 女人高潮潮喷娇喘18禁视频| svipshipincom国产片| 亚洲性夜色夜夜综合| 手机成人av网站| 黄色片一级片一级黄色片| 国产真人三级小视频在线观看| svipshipincom国产片| 琪琪午夜伦伦电影理论片6080| 人人妻人人看人人澡| 99riav亚洲国产免费| 久久久久久九九精品二区国产 | 麻豆一二三区av精品| 中文字幕另类日韩欧美亚洲嫩草| 999精品在线视频| 久久草成人影院| 草草在线视频免费看| 人妻丰满熟妇av一区二区三区| 久久中文字幕人妻熟女| 丝袜在线中文字幕| 欧美一级毛片孕妇| 成人亚洲精品一区在线观看| 精品福利观看| 免费在线观看成人毛片| 欧美黑人精品巨大| 一本综合久久免费| 大香蕉久久成人网| 在线观看一区二区三区| 成人一区二区视频在线观看| 视频在线观看一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲在线自拍视频| 叶爱在线成人免费视频播放| 日韩欧美 国产精品| 校园春色视频在线观看| 国产精品精品国产色婷婷| 久久香蕉精品热| 亚洲av电影不卡..在线观看| 日韩三级视频一区二区三区| 精品不卡国产一区二区三区| 日韩一卡2卡3卡4卡2021年|