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

    鐵皮石斛不同部位內(nèi)生細(xì)菌群落高通量分析

    2017-12-20 00:53:51陳澤斌王定斌郭麗紅
    河南農(nóng)業(yè)科學(xué) 2017年11期
    關(guān)鍵詞:鐵皮石斛內(nèi)生

    陳澤斌,李 冰,高 熹,王定斌,靳 松,郭麗紅

    (1.昆明學(xué)院 農(nóng)學(xué)院,云南 昆明 650214; 2.中國科學(xué)院 昆明動物研究所,云南 昆明 650223; 3.云南農(nóng)業(yè)大學(xué) 植物保護(hù)學(xué)院,云南 昆明 650201; 4.宣威市阿都鄉(xiāng)農(nóng)業(yè)綜合服務(wù)中心,云南 宣威 655425; 5.云南省高校特色生物資源開發(fā)與利用重點實驗室,云南 昆明 650214)

    鐵皮石斛不同部位內(nèi)生細(xì)菌群落高通量分析

    陳澤斌1,李 冰2*,高 熹3,王定斌4,靳 松1,郭麗紅5

    (1.昆明學(xué)院 農(nóng)學(xué)院,云南 昆明 650214; 2.中國科學(xué)院 昆明動物研究所,云南 昆明 650223; 3.云南農(nóng)業(yè)大學(xué) 植物保護(hù)學(xué)院,云南 昆明 650201; 4.宣威市阿都鄉(xiāng)農(nóng)業(yè)綜合服務(wù)中心,云南 宣威 655425; 5.云南省高校特色生物資源開發(fā)與利用重點實驗室,云南 昆明 650214)

    為揭示鐵皮石斛不同器官中內(nèi)生細(xì)菌的分布規(guī)律,應(yīng)用Illumina測序平臺的MiSeq高通量測序儀對鐵皮石斛根、莖、葉中內(nèi)生細(xì)菌的16S rDNA-V4區(qū)擴(kuò)增子進(jìn)行測序,對測得數(shù)據(jù)進(jìn)行生物信息學(xué)分析,評價物種Alpha多樣性和Beta多樣性。經(jīng)FLASH軟件對序列進(jìn)行拼接,莖、葉、根樣品分別獲得64 851、77 499、52 157條原始序列,經(jīng)Qiime軟件過濾,再剔除宿主的葉綠體和線粒體序列,最終分別得到64 742、77 351、52 064條有效序列,在97%的序列相似性水平上,這些有效序列被Uparse軟件劃分為61、134、112個操作分類單元(OTUs)。序列比對結(jié)果表明,假單胞菌屬(Pseudomonas)和食酸菌屬(Acidovorax)細(xì)菌在根、莖、葉中都是優(yōu)勢菌,瘤胃球菌屬(Ruminococcus)、Faecalibacterium屬、甲基桿菌屬(Methylobacterium)、希瓦氏菌屬(Shewanella)細(xì)菌表現(xiàn)出組織專一性。豐度聚類分析及主成分分析表明,根和葉中內(nèi)生細(xì)菌群落結(jié)構(gòu)最為相似。鐵皮石斛植株內(nèi)蘊含豐富的細(xì)菌資源,其分布受組織結(jié)構(gòu)及成分因素影響表現(xiàn)出器官差異性。

    鐵皮石斛; 內(nèi)生細(xì)菌; 多樣性; 高通量分析

    石斛可分為黃草、金釵、馬鞭等數(shù)十種,鐵皮石斛(Dendrobiumcandidum)為石斛中經(jīng)濟(jì)價值最高的一種,有“中華九大仙草之首”的美稱,分布于安徽西南部、浙江東部、福建西部、廣西西北部、四川、云南東南部海拔1 600 m的山地半陰濕的巖石上[1-2]。以鐵皮石斛的莖入藥,具有增加胃液分泌、助消化、促進(jìn)腸胃蠕動、降血糖、降血壓、增強機體免疫力等作用[3-6]。

    植物內(nèi)生菌是指整個生活周期都存在于宿主植物體內(nèi)的細(xì)菌、真菌和放線菌,目前還沒有發(fā)現(xiàn)不存在內(nèi)生菌的植物[7]。內(nèi)生菌在石斛屬植物促進(jìn)生長[8]、抑制病原菌[9]、影響宿主代謝成分[10-11]方面起著重要的作用,因此有必要對鐵皮石斛內(nèi)生菌多樣性進(jìn)行研究。洪群艷等[12]認(rèn)為,刺盤孢屬(Colletotrichum)為野生石斛莖和葉的優(yōu)勢種群,占總菌株數(shù)的39.5%。此外還發(fā)現(xiàn),石斛內(nèi)生真菌具有組織差異性和專一性。章華偉等[13]發(fā)現(xiàn),鐵皮石斛內(nèi)生細(xì)菌分布在芽孢桿菌屬(Bacillus)、土壤桿菌屬(Agrobacterium)、乳桿菌屬(Lactobacillus)、泛菌屬(Pantoea)、黃桿菌屬(Flavobacterium)、梭菌屬(Clostridium)、葡萄球菌屬(Staphylococcus)共7個屬。嚴(yán)亮等[14]發(fā)現(xiàn),Bacillusoceanisediminis、Bacillusberingensis和Bacilluslicheniformis3種菌同時存在于根部和葉部。

    目前,對于石斛內(nèi)生菌多樣性的研究仍是采用傳統(tǒng)的組織培養(yǎng)方法進(jìn)行,然而由于表面消毒過重、所用培養(yǎng)基營養(yǎng)成分單一、有氧的培養(yǎng)條件等因素,實際上只分離到真實存在微生物總數(shù)的0.1%~10%[15-16],隨著分子生物學(xué)技術(shù)的不斷發(fā)展,為了無偏差地了解內(nèi)生細(xì)菌群落的真實組成,可采用免培養(yǎng)技術(shù)手段直接從DNA層面解析群落中的物種組成。因此,本研究首次采用Illumina高通量測序平臺的 MiSeq測序儀對鐵皮石斛不同部位內(nèi)生細(xì)菌的16S rDNA擴(kuò)增子進(jìn)行二代測序,相較于傳統(tǒng)的組織培養(yǎng)法以及基于一代Sanger測序的非培養(yǎng)方法,測序成本更低,深度更大,檢測靈敏度更高,從而實現(xiàn)對鐵皮石斛不同部位內(nèi)生細(xì)菌種類的準(zhǔn)確調(diào)查,為豐富植物微生態(tài)學(xué)理論及基因工程菌的研究和應(yīng)用奠定基礎(chǔ)。

    1 材料和方法

    1.1 樣品采集

    2016年8月于云南省普洱市景東彝族自治縣鐵皮石斛人工栽培基地,采用3點取樣法挖取20株無明顯病蟲害的1年生栽培型鐵皮石斛,去除腐葉,將根、莖、葉分別放入無菌樣品袋中,混合均勻,低溫保鮮,于24 h內(nèi)進(jìn)行表面消毒處理,根、莖、葉樣品分別命名為SG、SJ、SY。

    1.2 表面消毒

    用無菌水將鐵皮石斛根、莖、葉沖洗干凈,吸干表面水分,再用75%乙醇浸泡30 s,0.1%升汞振蕩浸泡5 min,無菌水沖洗4次,無菌濾紙吸干表面水分,表面消毒效果檢驗參照文獻(xiàn)[17]。

    1.3 基因組DNA的提取

    采用改進(jìn)的CTAB法[18]提取表面徹底消毒的根、莖、葉組織基因組DNA,并用0.5%的瓊脂糖凝膠電泳檢查DNA的濃度和純度。

    1.4 16S rDNA-V4區(qū)的PCR擴(kuò)增

    以根、莖、葉基因組DNA為模板,PCR擴(kuò)增16S rDNA的V4可變區(qū),引物為515F(5′-GTGCCAGCMGCCGCGGTAA-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)[19],根、莖、葉引物序列前插入6個堿基組成的標(biāo)簽序列,以區(qū)分根、莖、葉樣品。PCR反應(yīng)體系參照文獻(xiàn)[19]。采用2.0%的瓊脂糖凝膠電泳檢測PCR產(chǎn)物并切膠回收目的條帶,送北京諾禾致源生物信息科技有限公司,用Illumina高通量測序平臺的 MiSeq測序儀測序。

    1.5 測序數(shù)據(jù)的質(zhì)控

    根據(jù)插入標(biāo)簽序列,將下機數(shù)據(jù)拆分為根、莖、葉數(shù)據(jù),剔除標(biāo)簽序列和引物序列,使用FLASH軟件[20]分別對不同樣品的讀長序列進(jìn)行拼接,Qiime軟件[21]過濾,經(jīng)過以上處理后得到的序列與數(shù)據(jù)庫(Golddatabase,http://drive5.com/uchime/uchime_download.html)進(jìn)行比對,UCHIME軟件[22]檢測嵌合體序列,并去除其中的嵌合體序列,得到最終的有效序列。

    1.6 測序數(shù)據(jù)的分析

    在97%的相似性水平上用Uparse軟件[23](Uparse v7.0.1001,http://drive5.com/uparse/)對有效序列進(jìn)行操作分類單元(OTUs)聚類,選取代表序列用RDP classifier軟件[24](V2.2,http://sourceforge.net/projects/rdp-classifier/)與GreenGene數(shù)據(jù)庫[25](http://greengenes.lbl.gov/cgi-bin/nph-index.cgi)進(jìn)行物種注釋分析,剔除宿主葉綠體和線粒體序列,利用Mothur軟件做稀釋度曲線分析,計算文庫覆蓋率,采用Chao 1指數(shù)及香農(nóng)(Shannon)指數(shù)評估根、莖、葉內(nèi)生細(xì)菌多樣性,計算方法參見文獻(xiàn)[26]。采用R語言分析繪制樣品OTUs的維恩圖,通過主成分分析(PCA)挖掘根、莖、葉內(nèi)生細(xì)菌群落結(jié)構(gòu)的差異。

    2 結(jié)果與分析

    2.1 測序數(shù)據(jù)的質(zhì)控及分析

    莖、葉、根樣品分別測得64 851、77 499、52 157條原始序列,經(jīng)過濾后得到的序列拼接總數(shù)分別為64 748、77 361、52 069條,去除出現(xiàn)頻度為1、無法被聚類到OTUs的低頻序列,再剔除宿主的葉綠體和線粒體序列,最終分別得到64 742、77 351、52 064條有效序列,在97%的序列相似性水平上,這些有效序列分別可劃分為61、134、112個OTUs,相對應(yīng)的文庫覆蓋率分別為99.91%、99.83%、99.78%(表1),物種豐度稀釋曲線見圖1,3個樣品稀釋曲線上升平緩,結(jié)合覆蓋率數(shù)據(jù)認(rèn)為,設(shè)計的測序深度已足夠,可以代表樣品中真實存在的大多數(shù)細(xì)菌種類,但仍有少數(shù)細(xì)菌種類未能檢測到。

    表1 鐵皮石斛樣品OTUs豐度和α多樣性指數(shù)

    圖1 物種豐度稀釋曲線

    2.2 不同部位內(nèi)生細(xì)菌OTUs分布

    根、莖、葉細(xì)菌文庫共含OTUs的數(shù)量為307個,特有的OTUs數(shù)量分別為112、61、134個,由此可看出,葉片內(nèi)生細(xì)菌多樣性最高,其次為根,最低為莖,這與Chao 1指數(shù)、香農(nóng)指數(shù)得出的結(jié)果一致(表1)。從圖2可知,有30個OTUs在根、莖、葉中均有分布,只占文庫OTUs總數(shù)的9.77%,說明不同部位內(nèi)生細(xì)菌的種類組成差異大;根和葉文庫共含OTUs的數(shù)量為246個,其中107個OTUs在根和葉中均有分布,占相應(yīng)OTUs總數(shù)的43.50%;葉和莖文庫共含OTUs的數(shù)量為195個,其中35個OTUs在葉和莖中均有分布,占相應(yīng)OTUs總數(shù)的17.95%;根和莖文庫共含OTUs的數(shù)量為173個,其中30個OTUs在根和莖中均有分布,占相應(yīng)OTUs總數(shù)的17.34%。可以看出,根和葉中內(nèi)生細(xì)菌的種類組成較相似。

    圖2 內(nèi)生細(xì)菌OTUs分布維恩圖

    2.3 不同部位內(nèi)生細(xì)菌多樣性

    選取在屬分類水平上相對豐度排名前10的屬,生成物種相對豐度分布圖(圖3)。在莖樣品中,內(nèi)生細(xì)菌優(yōu)勢屬為假單胞菌屬(Pseudomonas,54.84%),其次是食酸菌屬(Acidovorax,22.58%);在葉樣品中,內(nèi)生細(xì)菌優(yōu)勢屬為假單胞菌屬(Pseudomonas,26.00%),其次是鞘脂單胞菌屬(Sphingomonas,24.00%)和食酸菌屬(Acidovorax,10.00%);在根樣品中,內(nèi)生細(xì)菌優(yōu)勢屬為食酸菌屬(Acidovorax,36.96%),其次是假單胞菌屬(Pseudomonas,26.09%)和鞘脂單胞菌屬(Sphingomonas,10.87%)。可發(fā)現(xiàn)假單胞菌屬和食酸菌屬細(xì)菌在根、莖、葉中都是優(yōu)勢菌。此外還發(fā)現(xiàn),鹽單胞菌屬(Halomonas)、布赫納氏菌屬(Buchnera)在根、莖、葉中少量分布。瘤胃球菌屬(Ruminococcus,8.00%)和Faecalibacterium屬(8.00%)只在葉樣品中檢測到,甲基桿菌屬(Methylobacterium)和希瓦氏菌屬(Shewanella)只在根(2.17%、4.35%)和葉(8.00%、2.00%)樣品中檢測到,甲基暖菌屬(Methylocaldum,6.45%)只存在于莖部。

    圖3 內(nèi)生細(xì)菌屬水平上的相對豐度

    2.4 內(nèi)生細(xì)菌豐度聚類分析

    根據(jù)根、莖、葉在屬水平上的物種注釋及豐度信息,選取豐度排名前10的屬及其在不同部位中的豐度信息繪制熱圖(圖4),并從樣品間差異層面進(jìn)行聚類,發(fā)現(xiàn)樣品SG和SY處于同一分支,最為鄰近,而與樣品SJ距離較遠(yuǎn),說明根和葉中內(nèi)生細(xì)菌群落結(jié)構(gòu)最為相似。

    圖4 內(nèi)生細(xì)菌OTUs豐度聚類分析

    2.5 內(nèi)生細(xì)菌群落主成分分析

    主成分分析表明(圖5),主成分1(PC1)和主成分2(PC2)對樣品差異性的解釋度分別為59.80%和40.19%,合計為99.99%,是差異的主要來源。樣品SG和SY分別位于主成分PC2軸的正負(fù)兩側(cè),但同時都位于主成分PC1軸的正值區(qū)域;樣品SG和SJ不僅分別位于主成分PC2軸的正負(fù)兩側(cè),而且也位于主成分PC1軸的正負(fù)兩側(cè),且間隔較大的距離,說明根與莖間的主成分變異程度比根與葉的顯著,總體看來,內(nèi)生細(xì)菌群落的分布隨部位的不同有顯著變化。

    圖5 內(nèi)生細(xì)菌群落主成分分析

    3 結(jié)論與討論

    隨著MiSeq平臺雙端PE250和PE300測序策略的發(fā)展,讀長不斷增長,使得基于MiSeq測序研究物種多樣性的準(zhǔn)確性大幅度提高。利用MiSeq平臺對16S rDNA的一個或多個高變區(qū)測序,具有測序深度高、利于鑒定低豐度群落物種以及費用低的特點,已成為研究微生物群落多樣性的首選之策。本研究首次采用二代高通量測序技術(shù)對鐵皮石斛內(nèi)生細(xì)菌16S rDNA-V4區(qū)進(jìn)行測序,經(jīng)序列比對發(fā)現(xiàn)在根、莖、葉中常見的優(yōu)勢細(xì)菌種類為變形菌門(Proteobacteria)的假單胞菌屬(Pseudomonas)和食酸菌屬(Acidovorax),此外還檢測到一些沒有報道過的低豐度細(xì)菌類群(所占比例小于10.00%),例如布赫納氏菌屬(Buchnera)、瘤胃球菌屬(Ruminococcus)、甲基桿菌屬(Methylobacterium)、Faecalibacterium屬、希瓦氏菌屬(Shewanella),顯示出了新一代測序技術(shù)的優(yōu)勢。在以往通過傳統(tǒng)組織分離方法對鐵皮石斛內(nèi)生細(xì)菌進(jìn)行分離鑒定的報道中,均發(fā)現(xiàn)優(yōu)勢菌為厚壁菌門(Firmicutes)的芽孢桿菌屬(Bacillus)細(xì)菌[27],與本研究的結(jié)果不一致,這可能與芽孢桿菌屬細(xì)菌在植物體內(nèi)多是以孢子形式存在有關(guān),而提取植物組織DNA目前多采用CTAB、SDS法,無法使其裂解,導(dǎo)致未能獲取到該類群細(xì)菌的基因組DNA。

    根、莖、葉細(xì)菌文庫共含OTUs的數(shù)量為307個,只有30個OTUs為根、莖、葉共有,只占文庫OTUs總數(shù)的9.77%,此外Chao 1指數(shù)、香農(nóng)指數(shù)得出的結(jié)果為葉>根>莖,說明內(nèi)生細(xì)菌的分布隨組織器官的不同而各有特點,具有明顯的器官差異性,這可能是由于不同器官的組織結(jié)構(gòu)不同、含有的化學(xué)成分不同造成的。瘤胃球菌屬和Faecalibacterium屬只在葉樣品中檢測到,甲基桿菌屬和希瓦氏菌屬只在根和葉樣品中檢測到,甲基暖菌屬只存在于莖部,表現(xiàn)出部分細(xì)菌種群的組織專一性,這與Petrini等[28]的研究結(jié)果一致,宿主植物受體蛋白能夠影響內(nèi)生細(xì)菌在對其侵入時的識別,同時內(nèi)生細(xì)菌的定殖也受到宿主植物的調(diào)控,這些因素共同造成了內(nèi)生細(xì)菌分布的專一性。本研究發(fā)現(xiàn),假單胞菌屬和食酸菌屬細(xì)菌在根、莖、葉中都是優(yōu)勢菌,推測這2個屬細(xì)菌的侵入能力強,且在植株體內(nèi)具有較強的競爭力。已有報道表明,假單胞菌屬具有一定的固氮功能,與蘭科植物共生,為其供給氮源。除此之外,假單胞菌屬細(xì)菌還可以降解土壤中有機磷農(nóng)藥殘留物,避免污染環(huán)境[29]。食酸菌屬細(xì)菌曾多次作為植物內(nèi)生細(xì)菌被報道,宿主植株包括水稻、煙草、油菜、薄荷、煙草等[30],但是關(guān)于其功能、代謝的研究目前還未見報道。

    在將16S rDNA-V4區(qū)測得序列與數(shù)據(jù)庫進(jìn)行比對的過程中發(fā)現(xiàn),一部分序列歸屬于宿主植物的葉綠體和線粒體,因此,在進(jìn)行后續(xù)分析之前,剔除了這部分無關(guān)的序列。本研究只選擇了普遍用于水體、土壤、糞便等環(huán)境樣品[31]細(xì)菌多樣性研究的16S rDNA-V4區(qū)作為測序區(qū)域,但從研究結(jié)果看來,選擇V4區(qū)進(jìn)行植物內(nèi)生細(xì)菌研究,效果并不理想,由于葉綠體和線粒體與細(xì)菌在系統(tǒng)發(fā)育上具有高度的同源性[32],產(chǎn)生了宿主污染的現(xiàn)象。下一步應(yīng)嘗試測序16S rDNA 9個可變區(qū)中的其他區(qū)域或者幾個區(qū)域,尋找能夠避開線粒體和葉綠體干擾的理想測序區(qū)域。

    [1] 于萌萌,吳令上,斯金平,等.基于多糖和黃酮類成分的鐵皮石斛優(yōu)良內(nèi)生真菌篩選[J].中國中藥雜志,2016,41(12):2208-2212.

    [2] 屠艷拉,肖芳,張敬澤,等.發(fā)生在鐵皮石斛栽培田中的2種黏菌病[J].浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版),2016,42(2):137-142.

    [3] 安紅強,范靜,梁易,等.鐵皮石斛泛素結(jié)合酶基因DoUBC24的克隆及表達(dá)分析[J].生物技術(shù)通訊,2016,27(5):643-648.

    [4] 屈小媛,武玉祥,張久磊.石斛化學(xué)成分與藥理活性研究進(jìn)展[J].天津農(nóng)業(yè)科學(xué),2015,21(4):19-23.

    [5] 肖春宏,黃飛燕,楊波.鼓槌石斛研究進(jìn)展[J].山西農(nóng)業(yè)科學(xué),2014,42(6):647-649.

    [6] 武榮花,栗燕.石斛資源保護(hù)研究進(jìn)展[J].河南農(nóng)業(yè)科學(xué),2009(12):9-12,17.

    [7] 王明月,陶茜,李克艷,等.鐵皮石斛內(nèi)生細(xì)菌群落結(jié)構(gòu)分析[J].西部林業(yè)科學(xué),2014,43(5):106-111.

    [8] 侯曉強,郭順星.鐵皮石斛促生長內(nèi)生真菌的篩選與鑒定[J].中國中藥雜志,2014,39(17):3232-3237.

    [9] 何勁,雷幫星,宋貞富,等.石斛內(nèi)生細(xì)菌DEB-2對5種辣椒病原真菌的抑制作用[J].植物保護(hù)學(xué)報,2014,41(2):157-162.

    [10] 楊紹周,吳毅歆,邵德林,等.鼓槌石斛內(nèi)生細(xì)菌分離、鑒定及功能分析[J].中國農(nóng)學(xué)通報,2014,30(25):171-176.

    [11] 趙昕梅,遠(yuǎn)凌威,張?zhí)K鋒,等.鐵皮石斛內(nèi)生真菌的分離鑒定及其促宿主生長作用[J].河南農(nóng)業(yè)科學(xué),2012,41(6):101-105.

    [12] 洪群艷,董文漢,徐文婷,等.幾種云南野生石斛內(nèi)生真菌的鑒定及分布[J].生物技術(shù)進(jìn)展,2012,2(3):190-194.

    [13] 章華偉,鐘超群.鐵皮石斛內(nèi)生細(xì)菌的分離和初步鑒定[J].湖北農(nóng)業(yè)科學(xué),2013,52(8):1811-1813.

    [14] 嚴(yán)亮,楊瑞娟,王橋美.云南鐵皮石斛內(nèi)生菌的分離與鑒定[J].云南農(nóng)業(yè)大學(xué)學(xué)報,2015,30(5):760-765.

    [15] Tholozan J L,Cappelier J M,Tissier J P,etal.Physiological characterization of viable-but-nonculturableCampylobacterjejunicells[J].Appl Environ Microbiol,1999,65(3):1110-1116.

    [16] Ammann R R,Ludwig W,Schleifer K H.Phylogenetic identification and in situ detection of individual microbial cells without cultivation[J].Microbiol Rev,1995,59(1):143-169.

    [17] 陳澤斌,代方平,寸林江,等.煙草內(nèi)生細(xì)菌分離方法的優(yōu)化研究[J].中國煙草學(xué)報,2014,20(1):90-95,102.

    [18] 邱服斌.培養(yǎng)方法與非培養(yǎng)方法對人參根內(nèi)生細(xì)菌的研究[D].北京:首都師范大學(xué),2007.

    [19] 陳澤斌,李冰,王定康,等.Illumina MiSeq高通量測序分析核桃內(nèi)生細(xì)菌多樣性[J].江蘇農(nóng)業(yè)學(xué)報,2015,31(5):1129-1133.

    [20] Magoc T,Salzberg S L.FLASH:Fast length adjustment of short reads to improve genome assemblies[J].Bioinformatics,2011,27(21):2957-2963.

    [21] Caporaso J G,Kuczynski J,Stombaugh J,etal.QIIME allows analysis of high-throughput community sequencing data[J].Nature Methods,2010,7(5):335-336.

    [22] Edgar R C,Haas B J,Clemente J C,etal.UCHIME improves sensitivity and speed of chimera detection[J].Bioinformatics,2011,27(16):2194-2200.

    [23] Edgar R C.UPARSE:Highly accurate OTU sequences from microbial amplicon reads[J].Nature Methods,2013,10(10):996-998.

    [24] Wang Q,Garrity G M,Tiedje J M,etal.Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J].Applied and Environmental Microbiology,2007,73(16):5261-5267.

    [25] DeSantis T Z,Hugenholtz P,Larsen N,etal.Greengenes,a chimera-checked 16S rRNA gene database and workbench compatible with ARB[J].Applied and Environmental Microbiology,2006,72(7):5069-5072.

    [26] 張志東,顧美英,王瑋,等.基于高通量測序的輻射污染區(qū)細(xì)菌群落特征分析[J].微生物學(xué)通報,2016,43(6):1218-1226.

    [27] 童文君,張禮,薛慶云,等.不同產(chǎn)地美花石斛內(nèi)生細(xì)菌分離及促生潛力比較[J].植物資源與環(huán)境學(xué)報,2014,23(1):16-23.

    [28] Petrini O,Fisher P J.A comparative study of fungal endophytes in xylem and whole stem ofPinussylvestrisandFagussylvatica[J].Trans Br Mycol Soc,1988,91(2):233-238.

    [29] 周小風(fēng).鐵皮石斛內(nèi)生細(xì)菌分布規(guī)律的研究[D].杭州:浙江理工大學(xué),2014.

    [30] 盧占慧,周如軍,袁月,等.人參內(nèi)生拮抗細(xì)菌分離、鑒定及其對人參菌核病抑菌作用研究[J].中國植保導(dǎo)刊,2016,36(3):5-10.

    [31] Sogin M L,Morrison H G,Huber J A,etal.Microbial diversity in the deep sea and the underexplored “rare biosphere”[J].Proceedings of the National Academy of Sciences,2006,103(32):12115-12120.

    [32] 胡楷,吳慶書.單細(xì)胞生物進(jìn)化研究的進(jìn)步[J].遺傳,2002,24(1):104-110.

    High Throughput Analysis of the Endophytic Bacterial Community in Different Parts of Dendrobium candidum

    CHEN Zebin1,LI Bing2*,GAO Xi3,WANG Dingbin4,JIN Song1,GUO Lihong5

    (1.Agriculture College of Kunming University,Kunming 650214,China;2.Kunming Institute of Zoology,the Chinese Academy of Sciences,Kunming 650223,China;3.Plant Protection College,Yunnan Agricultural University,Kunming 650201,China;4.Xuanwei City Adu Township Comprehensive Agricultural Service Center,Xuanwei 655425,China; 5.Key Laboratory of Special Biological Resource Development and Utilization of Universities in Yunnan Province,Kunming 650214,China)

    In order to disclose the distribution rule of endophytic bacteria in different organs ofDendrobiumcandidum,the MiSeq high-throughput sequencer of Illumina sequencing platform was used to sequence the amplicon of 16S rDNA-V4 of endophytic bacteria in the root,stem and leaves ofDendrobiumcandidum.Further Bioinformatics analysis was conducted to evaluate the Alpha and Beta diversities of species.After sequence assembly using FLASH software,64 851,77 499 and 52 157 original sequences were obtained from the stem,leaf and root samples,respectively.Filtered by Qiime software,the chloroplast and mitochondria sequences were eliminated,and 64 742,77 351 and 52 064 valid sequences were finally obtained.Based on 97% sequence similarity,these valid sequences were divided into 61,134 and 112 operational taxonomic units(OTUs).The sequence alignment result showed that bothPseudomonasandAcidovoraxbacteria were in the root,stem and leaves,theRuminococcu,Faecalibacterium,MethylobacteriumandShewanellabacteria presented tissue specificity.The abundance clustering analysis and main component analysis showed that the endophytic bacteria in root and leaves had the most similar community structure.TheDendrobiumcandidumplant contains rich bacterial resources,which present organ differentiation under the influence of organizational structure and composition factors.

    Dendrobiumcandidum; endophytic bacteria; diversity; high throughput analysis

    S567.239

    A

    1004-3268(2017)11-0098-06

    2017-06-27

    國家自然科學(xué)基金項目(31460491);云南省科技廳基礎(chǔ)研究項目(2013FD040);云南省教育廳科學(xué)研究項目(2014Y390);昆明學(xué)院大學(xué)生科研項目(XJD16081);昆明學(xué)院人才引進(jìn)項目(YJL14005);云南省高校特色生物資源開發(fā)與利用重點實驗室開放基金項目(GXKJ201621);云南省高校優(yōu)勢特色重點學(xué)科(生態(tài)學(xué))建設(shè)項目(05000511311)

    陳澤斌(1985-),男,云南昆明人,副教授,博士,主要從事環(huán)境微生物的教學(xué)及科研工作。

    E-mail:zbchenkmu@163.com

    *

    李 冰(1973-),女,遼寧沈陽人,高級工程師,碩士,主要從事非臨床藥物安全性評價工作。

    E-mail:libing@mail.kiz.ac.cn

    猜你喜歡
    鐵皮石斛內(nèi)生
    飛鼠與石斛
    鐵皮俠的裝備
    鐵皮園
    趣味(語文)(2020年4期)2020-07-27 01:43:14
    植物內(nèi)生菌在植物病害中的生物防治
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    35 種石斛蘭觀賞價值評價
    “黨建+”激活鄉(xiāng)村發(fā)展內(nèi)生動力
    授人以漁 激活脫貧內(nèi)生動力
    金釵石斛化學(xué)成分的研究
    中成藥(2018年5期)2018-06-06 03:11:56
    鐵皮石斛家庭種植技術(shù)探索
    99精品欧美一区二区三区四区| 桃色一区二区三区在线观看| 一进一出抽搐动态| 观看免费一级毛片| 国产黄片美女视频| 人人妻人人澡欧美一区二区| or卡值多少钱| 三级毛片av免费| 久久久精品大字幕| 国产极品精品免费视频能看的| 亚洲一区二区三区色噜噜| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区黑人| 午夜福利成人在线免费观看| 老熟妇乱子伦视频在线观看| 一a级毛片在线观看| 观看美女的网站| 99热6这里只有精品| 最近最新中文字幕大全电影3| 欧美激情久久久久久爽电影| 中文字幕人成人乱码亚洲影| 亚洲七黄色美女视频| 美女高潮的动态| 99精品在免费线老司机午夜| 亚洲成人久久性| 欧美黄色片欧美黄色片| tocl精华| 亚洲天堂国产精品一区在线| 国产一区二区三区视频了| 亚洲,欧美精品.| 亚洲精品久久国产高清桃花| 真实男女啪啪啪动态图| 999久久久精品免费观看国产| 午夜免费观看网址| 国产精品久久视频播放| 国产综合懂色| 日本黄色视频三级网站网址| 一二三四社区在线视频社区8| 亚洲精品久久国产高清桃花| 男女那种视频在线观看| 亚洲成人久久性| 免费高清视频大片| 亚洲欧美日韩高清专用| 亚洲人成网站高清观看| 亚洲熟女毛片儿| 他把我摸到了高潮在线观看| 精品久久久久久久末码| 亚洲欧洲精品一区二区精品久久久| 欧美激情久久久久久爽电影| 香蕉av资源在线| 精品久久久久久久末码| 国产精品一区二区免费欧美| av在线蜜桃| 国产成年人精品一区二区| 国产精品久久久久久人妻精品电影| 午夜福利视频1000在线观看| 欧美色欧美亚洲另类二区| 波多野结衣巨乳人妻| 国内精品一区二区在线观看| 欧美日韩黄片免| 国产野战对白在线观看| 宅男免费午夜| 岛国视频午夜一区免费看| 19禁男女啪啪无遮挡网站| 国产精品日韩av在线免费观看| 久久久精品欧美日韩精品| 给我免费播放毛片高清在线观看| 俺也久久电影网| 久久精品综合一区二区三区| 午夜福利免费观看在线| 久久天堂一区二区三区四区| or卡值多少钱| 成人无遮挡网站| 非洲黑人性xxxx精品又粗又长| 天堂影院成人在线观看| 老熟妇仑乱视频hdxx| 变态另类丝袜制服| 女警被强在线播放| 色播亚洲综合网| 男女视频在线观看网站免费| www.熟女人妻精品国产| 国产视频一区二区在线看| 1024香蕉在线观看| 久99久视频精品免费| 麻豆一二三区av精品| 男人舔女人的私密视频| 国产欧美日韩精品亚洲av| 欧美成人一区二区免费高清观看 | 黑人欧美特级aaaaaa片| 国内久久婷婷六月综合欲色啪| 中文字幕最新亚洲高清| 国产不卡一卡二| 亚洲精品色激情综合| 国产亚洲欧美在线一区二区| 性色av乱码一区二区三区2| 日日干狠狠操夜夜爽| 国产精华一区二区三区| 中出人妻视频一区二区| 国内精品久久久久精免费| 国产成人精品久久二区二区免费| 亚洲专区国产一区二区| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 国产一区二区在线观看日韩 | 97碰自拍视频| 国产精品香港三级国产av潘金莲| 国产精品,欧美在线| 日本免费a在线| 最近最新免费中文字幕在线| 校园春色视频在线观看| 欧美日韩综合久久久久久 | 国内久久婷婷六月综合欲色啪| 在线观看免费视频日本深夜| 老汉色av国产亚洲站长工具| 在线免费观看不下载黄p国产 | 老鸭窝网址在线观看| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 亚洲性夜色夜夜综合| 亚洲欧美精品综合久久99| h日本视频在线播放| 免费电影在线观看免费观看| 非洲黑人性xxxx精品又粗又长| 又黄又粗又硬又大视频| 丰满人妻一区二区三区视频av | 97超级碰碰碰精品色视频在线观看| netflix在线观看网站| 日韩成人在线观看一区二区三区| 国产精品久久久久久人妻精品电影| 制服人妻中文乱码| 我的老师免费观看完整版| 国产亚洲欧美在线一区二区| 校园春色视频在线观看| 国产爱豆传媒在线观看| 99国产精品一区二区蜜桃av| 国产探花在线观看一区二区| 亚洲乱码一区二区免费版| 51午夜福利影视在线观看| 午夜成年电影在线免费观看| 身体一侧抽搐| а√天堂www在线а√下载| 亚洲第一电影网av| 老汉色∧v一级毛片| 午夜福利免费观看在线| 色综合欧美亚洲国产小说| 国语自产精品视频在线第100页| 午夜免费观看网址| 在线免费观看不下载黄p国产 | 久9热在线精品视频| 人人妻人人看人人澡| 欧美av亚洲av综合av国产av| 亚洲国产欧美一区二区综合| 免费一级毛片在线播放高清视频| av天堂在线播放| 亚洲av熟女| 国产 一区 欧美 日韩| 日韩欧美 国产精品| 亚洲av电影在线进入| 99久久综合精品五月天人人| 中文亚洲av片在线观看爽| 哪里可以看免费的av片| 亚洲精品456在线播放app | 国产日本99.免费观看| 婷婷亚洲欧美| 亚洲av日韩精品久久久久久密| 婷婷丁香在线五月| 18禁观看日本| 欧美成人一区二区免费高清观看 | 亚洲欧美一区二区三区黑人| 日日摸夜夜添夜夜添小说| 搡老岳熟女国产| a在线观看视频网站| 99国产精品99久久久久| 亚洲片人在线观看| e午夜精品久久久久久久| 成年版毛片免费区| 九九热线精品视视频播放| 美女扒开内裤让男人捅视频| 国产乱人视频| 国模一区二区三区四区视频 | 两个人视频免费观看高清| 国产 一区 欧美 日韩| 99国产极品粉嫩在线观看| 久久亚洲真实| 好看av亚洲va欧美ⅴa在| 一进一出好大好爽视频| 成人鲁丝片一二三区免费| 午夜两性在线视频| 偷拍熟女少妇极品色| www.熟女人妻精品国产| 日韩成人在线观看一区二区三区| 久久午夜亚洲精品久久| 老熟妇乱子伦视频在线观看| 中文字幕熟女人妻在线| 人妻丰满熟妇av一区二区三区| 欧美日本亚洲视频在线播放| 别揉我奶头~嗯~啊~动态视频| 国产精品99久久久久久久久| 长腿黑丝高跟| 亚洲第一欧美日韩一区二区三区| 日韩免费av在线播放| 免费人成视频x8x8入口观看| 日日夜夜操网爽| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 丁香欧美五月| 国产淫片久久久久久久久 | 国产精品久久电影中文字幕| 国产黄片美女视频| 琪琪午夜伦伦电影理论片6080| 色av中文字幕| 国产av在哪里看| 在线观看午夜福利视频| 免费看光身美女| 国产av不卡久久| 国产v大片淫在线免费观看| 一a级毛片在线观看| 俄罗斯特黄特色一大片| 国产精品女同一区二区软件 | 欧美乱妇无乱码| 亚洲男人的天堂狠狠| 国产成+人综合+亚洲专区| 亚洲 国产 在线| 十八禁网站免费在线| 又大又爽又粗| 国产精品一区二区免费欧美| 老汉色∧v一级毛片| 宅男免费午夜| 精品久久久久久久末码| 欧美另类亚洲清纯唯美| 久久精品亚洲精品国产色婷小说| 免费看a级黄色片| 亚洲国产欧美人成| 亚洲人成电影免费在线| 高潮久久久久久久久久久不卡| 成人特级黄色片久久久久久久| 亚洲精品色激情综合| 国产美女午夜福利| 亚洲 欧美 日韩 在线 免费| 久久天躁狠狠躁夜夜2o2o| 在线国产一区二区在线| 国产欧美日韩精品亚洲av| 国产激情久久老熟女| 男人舔奶头视频| 免费观看的影片在线观看| 超碰成人久久| 偷拍熟女少妇极品色| 啦啦啦观看免费观看视频高清| www日本黄色视频网| 青草久久国产| 亚洲专区国产一区二区| 一个人观看的视频www高清免费观看 | 久久久国产成人精品二区| 身体一侧抽搐| 禁无遮挡网站| 又紧又爽又黄一区二区| 成人午夜高清在线视频| 国产亚洲精品av在线| 97超视频在线观看视频| 黄色丝袜av网址大全| 亚洲真实伦在线观看| 免费观看的影片在线观看| 又黄又爽又免费观看的视频| 免费看美女性在线毛片视频| netflix在线观看网站| 亚洲18禁久久av| 亚洲国产精品久久男人天堂| 亚洲成人精品中文字幕电影| 日韩有码中文字幕| 12—13女人毛片做爰片一| 国产精品久久视频播放| 成人国产综合亚洲| 搡老妇女老女人老熟妇| 欧美日本视频| 嫩草影院精品99| 怎么达到女性高潮| 亚洲国产欧洲综合997久久,| 国产高清三级在线| 18禁国产床啪视频网站| 久久伊人香网站| 午夜影院日韩av| 成人性生交大片免费视频hd| 两个人的视频大全免费| 亚洲色图 男人天堂 中文字幕| 国语自产精品视频在线第100页| 制服丝袜大香蕉在线| 美女被艹到高潮喷水动态| 国产精品免费一区二区三区在线| svipshipincom国产片| 亚洲欧美日韩无卡精品| 三级毛片av免费| 麻豆久久精品国产亚洲av| 这个男人来自地球电影免费观看| 一个人看的www免费观看视频| 亚洲九九香蕉| 一级毛片高清免费大全| 亚洲av免费在线观看| 精品久久久久久成人av| 精华霜和精华液先用哪个| 久久这里只有精品19| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看 | av国产免费在线观看| 国产成人av教育| 两个人的视频大全免费| 午夜精品久久久久久毛片777| 欧美日韩综合久久久久久 | 99国产精品一区二区三区| 欧美一区二区国产精品久久精品| 久久久精品大字幕| 精品国产美女av久久久久小说| 欧美一级毛片孕妇| 欧美在线黄色| 美女扒开内裤让男人捅视频| 国内揄拍国产精品人妻在线| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利18| 老司机午夜十八禁免费视频| 午夜福利在线在线| 日本撒尿小便嘘嘘汇集6| 91字幕亚洲| 成人性生交大片免费视频hd| 精品电影一区二区在线| 国产成人aa在线观看| 在线观看日韩欧美| 午夜福利免费观看在线| 免费在线观看亚洲国产| 日本黄色视频三级网站网址| 1024香蕉在线观看| 国产精品久久视频播放| 国产激情欧美一区二区| 亚洲国产欧美网| 欧美中文日本在线观看视频| 国产人伦9x9x在线观看| 久久精品人妻少妇| 男插女下体视频免费在线播放| 天天一区二区日本电影三级| 久久天堂一区二区三区四区| 99久久精品国产亚洲精品| 亚洲无线观看免费| 久久久久久久久久黄片| aaaaa片日本免费| 亚洲欧美精品综合一区二区三区| 听说在线观看完整版免费高清| xxx96com| 亚洲电影在线观看av| 欧美在线一区亚洲| 日韩大尺度精品在线看网址| 少妇丰满av| 亚洲欧洲精品一区二区精品久久久| 91在线观看av| 欧美又色又爽又黄视频| 男女之事视频高清在线观看| 制服人妻中文乱码| 丁香欧美五月| 在线观看午夜福利视频| 99热精品在线国产| 日韩 欧美 亚洲 中文字幕| 99国产精品一区二区蜜桃av| 99热精品在线国产| 无遮挡黄片免费观看| 中国美女看黄片| 国产熟女xx| 亚洲,欧美精品.| 桃色一区二区三区在线观看| 国产亚洲av高清不卡| 成人特级av手机在线观看| 成人鲁丝片一二三区免费| 国产精品久久久久久久电影 | 香蕉av资源在线| 欧美日韩福利视频一区二区| 亚洲最大成人中文| 国产午夜精品久久久久久| 人妻久久中文字幕网| 色哟哟哟哟哟哟| 欧美激情在线99| 国产熟女xx| 日本免费a在线| 在线十欧美十亚洲十日本专区| 最近最新中文字幕大全免费视频| 国产免费男女视频| 欧美在线一区亚洲| 国产视频一区二区在线看| 一本久久中文字幕| 老汉色∧v一级毛片| 99精品久久久久人妻精品| 精品电影一区二区在线| e午夜精品久久久久久久| 精品久久久久久,| 又爽又黄无遮挡网站| 亚洲成人中文字幕在线播放| 国内久久婷婷六月综合欲色啪| 一进一出好大好爽视频| 日韩大尺度精品在线看网址| 综合色av麻豆| 国产精品亚洲美女久久久| 精品国产乱码久久久久久男人| 好看av亚洲va欧美ⅴa在| 黄频高清免费视频| 熟妇人妻久久中文字幕3abv| 哪里可以看免费的av片| 午夜免费观看网址| 母亲3免费完整高清在线观看| 淫妇啪啪啪对白视频| 国产高潮美女av| 最近最新中文字幕大全免费视频| 国产精品电影一区二区三区| 热99在线观看视频| 欧美色欧美亚洲另类二区| 日韩中文字幕欧美一区二区| 久久中文字幕人妻熟女| 九色国产91popny在线| 身体一侧抽搐| 国产精品 国内视频| 91在线观看av| 亚洲18禁久久av| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 青草久久国产| 色哟哟哟哟哟哟| 国产成人av教育| 午夜影院日韩av| 色av中文字幕| 岛国视频午夜一区免费看| 看免费av毛片| 叶爱在线成人免费视频播放| 国产黄片美女视频| 真人一进一出gif抽搐免费| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av在线| 成人鲁丝片一二三区免费| 最新中文字幕久久久久 | 在线a可以看的网站| 高清毛片免费观看视频网站| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 18禁黄网站禁片免费观看直播| 亚洲精品456在线播放app | 一级毛片高清免费大全| 精品免费久久久久久久清纯| 九九在线视频观看精品| 亚洲人与动物交配视频| 国产真人三级小视频在线观看| 日本三级黄在线观看| 美女黄网站色视频| 真人做人爱边吃奶动态| 欧美成人性av电影在线观看| 亚洲 欧美 日韩 在线 免费| 很黄的视频免费| 美女大奶头视频| 午夜成年电影在线免费观看| АⅤ资源中文在线天堂| 国产熟女xx| 免费搜索国产男女视频| 国产99白浆流出| 亚洲人成伊人成综合网2020| 日本免费a在线| 在线观看舔阴道视频| 欧美在线一区亚洲| av在线蜜桃| e午夜精品久久久久久久| 国产精品香港三级国产av潘金莲| 91老司机精品| av在线天堂中文字幕| 白带黄色成豆腐渣| 久久久国产成人免费| 日本与韩国留学比较| www.精华液| 亚洲av美国av| 国产亚洲精品综合一区在线观看| 久久亚洲真实| 国产精品综合久久久久久久免费| 亚洲精品456在线播放app | 亚洲成av人片在线播放无| 51午夜福利影视在线观看| 久久香蕉精品热| 三级毛片av免费| 国产又色又爽无遮挡免费看| 精品久久久久久久久久免费视频| or卡值多少钱| 老司机在亚洲福利影院| 成人性生交大片免费视频hd| 欧美在线黄色| 午夜精品一区二区三区免费看| 综合色av麻豆| 亚洲欧美精品综合一区二区三区| 可以在线观看的亚洲视频| 日韩成人在线观看一区二区三区| 最近最新中文字幕大全电影3| 成人特级av手机在线观看| 成人国产综合亚洲| 一级作爱视频免费观看| 亚洲国产精品sss在线观看| 制服人妻中文乱码| 黄色视频,在线免费观看| 久久精品亚洲精品国产色婷小说| 午夜亚洲福利在线播放| 757午夜福利合集在线观看| 国产亚洲av嫩草精品影院| 一区二区三区高清视频在线| 狂野欧美白嫩少妇大欣赏| 日韩免费av在线播放| 色综合欧美亚洲国产小说| 五月玫瑰六月丁香| 国产午夜精品论理片| www.www免费av| 日韩精品青青久久久久久| 亚洲国产中文字幕在线视频| АⅤ资源中文在线天堂| 亚洲国产精品合色在线| 真人一进一出gif抽搐免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲在线自拍视频| 精品免费久久久久久久清纯| 日本黄大片高清| 狂野欧美激情性xxxx| 欧美黑人欧美精品刺激| 亚洲一区二区三区不卡视频| 美女cb高潮喷水在线观看 | 最新在线观看一区二区三区| 国产免费av片在线观看野外av| 99国产精品99久久久久| 亚洲精华国产精华精| 久久久久久人人人人人| 亚洲国产中文字幕在线视频| 色综合站精品国产| 三级男女做爰猛烈吃奶摸视频| 99精品久久久久人妻精品| 国产一区二区在线观看日韩 | 日韩国内少妇激情av| 男插女下体视频免费在线播放| 不卡一级毛片| 一级作爱视频免费观看| 亚洲第一欧美日韩一区二区三区| 国产精品综合久久久久久久免费| 国产精品美女特级片免费视频播放器 | 成人三级做爰电影| 在线观看一区二区三区| www.精华液| 男女下面进入的视频免费午夜| av黄色大香蕉| 88av欧美| 美女被艹到高潮喷水动态| 久久久久久久午夜电影| 欧美日本视频| 91在线精品国自产拍蜜月 | 欧美日韩黄片免| 操出白浆在线播放| 精品人妻1区二区| 又爽又黄无遮挡网站| av欧美777| 国产av一区在线观看免费| 免费电影在线观看免费观看| 免费观看精品视频网站| 久久久久久大精品| 国产免费男女视频| 99久久无色码亚洲精品果冻| 全区人妻精品视频| 精品熟女少妇八av免费久了| av天堂在线播放| 欧美另类亚洲清纯唯美| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 听说在线观看完整版免费高清| 不卡av一区二区三区| 一本久久中文字幕| 18禁黄网站禁片免费观看直播| 美女被艹到高潮喷水动态| 中亚洲国语对白在线视频| 97碰自拍视频| 久久国产精品人妻蜜桃| 狠狠狠狠99中文字幕| 12—13女人毛片做爰片一| 国产久久久一区二区三区| 天天添夜夜摸| 亚洲一区二区三区色噜噜| 亚洲男人的天堂狠狠| 久久久精品欧美日韩精品| 成年女人看的毛片在线观看| 亚洲第一欧美日韩一区二区三区| 女生性感内裤真人,穿戴方法视频| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 亚洲精品美女久久av网站| 国产黄片美女视频| 亚洲欧美日韩高清专用| 波多野结衣高清无吗| 91麻豆av在线| 黄片小视频在线播放| 精品欧美国产一区二区三| 又黄又粗又硬又大视频| 亚洲aⅴ乱码一区二区在线播放| 在线观看66精品国产| 欧美xxxx黑人xx丫x性爽| 婷婷精品国产亚洲av在线| 搡老妇女老女人老熟妇| 一个人看的www免费观看视频| 欧美成狂野欧美在线观看| 男人舔奶头视频| 精品久久久久久成人av| 九九久久精品国产亚洲av麻豆 | 中文字幕久久专区| 免费看美女性在线毛片视频| 国产黄片美女视频| 久9热在线精品视频| 亚洲精品在线观看二区| 亚洲国产中文字幕在线视频| 十八禁网站免费在线| 韩国av一区二区三区四区| a级毛片在线看网站| 午夜免费成人在线视频| 熟女少妇亚洲综合色aaa.| 亚洲av美国av|