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

    煙草野火病葉際微生物群落結(jié)構(gòu)與多樣性

    2024-12-31 00:00:00郭濤汪漢成余知和蔡劉體王豐陳興江
    廣西植物 2024年10期
    關(guān)鍵詞:高通量測(cè)序病葉群落結(jié)構(gòu)

    摘 要:" 葉際微生物種群穩(wěn)定與煙株健康密切相關(guān)。為揭示感染野火病煙株葉際微生態(tài)特征,該文采用Illumina高通量測(cè)序技術(shù)對(duì)不同發(fā)病程度煙株葉片的病斑和健康組織進(jìn)行微生物群落結(jié)構(gòu)與多樣性分析。結(jié)果表明:(1)病斑組織與健康組織優(yōu)勢(shì)菌為變形菌門(Proteobacteria)和子囊菌門(Ascomycota),并且其在病斑組織的相對(duì)豐度顯著高于健康組織;優(yōu)勢(shì)屬為假單胞菌屬(Pseudomonas)、泛菌屬(Pantoea)、鏈格孢屬(Alternaria)、小不整球殼屬(Plectosphaerella)和錐蓋傘屬(Conocybe),僅假單胞菌屬在病斑組織的相對(duì)豐度顯著高于健康組織。(2)病斑組織細(xì)菌群落多樣性和豐富度高于健康組織,并且隨病害程度的增加,多樣性指數(shù)和豐富度指數(shù)分別呈先減后增和先增后減的趨勢(shì);病斑組織真菌群落豐富度和多樣性指數(shù)隨病害程度的增加先增后減,輕度和中度病葉病斑組織真菌群落多樣性指數(shù)和豐富度高于健康組織,但重度病葉病斑組織真菌群落多樣性指數(shù)和豐富度低于健康組織。(3)病斑和健康組織細(xì)菌的基因功能信息及相對(duì)豐度基本類似,主要為代謝、遺傳信息處理和環(huán)境信息處理3類;真菌的優(yōu)勢(shì)功能類群為植物病原菌、動(dòng)物病原菌-內(nèi)生真菌-植物病原菌-木質(zhì)腐生真菌、未定義腐生真菌、木質(zhì)腐生真菌和植物病原菌-木質(zhì)腐生真菌。該研究結(jié)果為認(rèn)識(shí)煙草野火病演變規(guī)律及葉際微生物群落特征提供了科學(xué)依據(jù)。

    關(guān)鍵詞: 煙草野火病, 高通量測(cè)序, 病葉, 群落結(jié)構(gòu), 多樣性指數(shù), 豐富度指數(shù)

    中圖分類號(hào):" Q948.12

    文獻(xiàn)標(biāo)識(shí)碼:" A

    文章編號(hào):" 1000-3142(2024)10-1894-11

    收稿日期:" 2024-03-24

    接受日期:" 2024-05-27

    基金項(xiàng)目:" 中國(guó)煙草總公司科技項(xiàng)目[110202101048(LS-08),110202001035(LS-04)]; 國(guó)家自然科學(xué)基金(32160522,31960550); 貴州省科技基金(黔科合基礎(chǔ)-ZK [2021]重點(diǎn)036); 貴州省“百層次”創(chuàng)新型人才(黔科合平臺(tái)人才-GCC [2022]028-1,黔科合平臺(tái)人才-GCC [2023]108); 貴州省科技創(chuàng)新人才團(tuán)隊(duì)項(xiàng)目(黔科合平臺(tái)人才-CXTD [2023]021); 貴州省應(yīng)用技術(shù)研究與開發(fā)資金后補(bǔ)助項(xiàng)目; 中國(guó)煙草總公司貴州省公司科技項(xiàng)目(2024XM06)。

    第一作者: 郭濤(1998—),碩士研究生,研究方向?yàn)樯飳W(xué),(E-mail)2634150274@qq.com。

    *通信作者:" 余知和,博士,教授,主要從事真菌生物學(xué)及植物病害研究,(E-mail)zhiheyu@hotmail.com;王豐,博士,研究員,主要從事煙草栽培與微生物學(xué)研究,(E-mail)584228992@qq.com。

    Microbial community structure and diversity of leaf

    phyllosphere in tobacco plants infected with wildfire disease

    GUO Tao1, 2, WANG Hancheng2, YU Zhihe1*, CAI Liuti2,

    WANG Feng2*, CHEN Xingjiang2

    ( 1. College of Life Sciences, Yangtze University, Jingzhou 434025, Hubei, China; 2. Guizhou Academy of Tobacco Sciences, Guiyang 550081, China )

    Abstract:" The stability of phyllosphere microbial community is closely related to plant health. In order to understand the microecological characteristics of tobacco leaves infected with wildfire disease, community structure and diversity of phyllosphere microbiome in diseased spots and healthy tissues of tobacco leaves with different degrees of disease were analyzed using high-throughput sequencing technology in this paper. The results were as follows: (1) The dominant microorganisms in diseased and healthy tissues were Proteobacteria and Ascomycota, and their relative abundance in diseased tissues was significantly higher than that in healthy tissues. The dominant genera were Pseudomonas, Pantoea, Alternaria, Plectosphaerella and Conocybe, and only the relative abundance of Pseudomonas in diseased tissues was significantly higher than that in healthy tissues. (2) The diversity and richness of bacterial communities in diseased spots were higher than those in healthy tissues, and the diversity index showed a trend of decreasing and then increasing with the increase of the disease degree, and the richness index showed a trend of increasing and then decreasing. The richness and diversity index of fungal communities in diseased spots showed a trend of increasing and then decreasing with the increase of the disease degree, and the diversity index and richness of fungal communities in diseased spot tissues of mildly and moderately diseased leaves were higher than those in healthy tissues, but the diversity index and richness of fungal communities in heavily diseased tissues were lower than those in healthy tissues. (3) Gene functional information and relative abundance of bacteria in diseased spots and healthy tissues were generally similar, primarily were metabolism, genetic information processing and environmental information processing. The dominant functional taxa of fungus were plant pathogen, animal pathogen-endophyte-plant pathogen-wood saprotroph, undefined saprotroph, wood saprotroph and plant pathogen-wood saprotroph. The results provide a reference for the understanding of evolutionary pattern of tobacco wildfire disease and characterization of the phyllosphere microbial community.

    Key words: tobacco wildfire disease, high-throughput sequencing, diseased leaf, community structure, diversity index, richness index

    煙草是我國(guó)重要的經(jīng)濟(jì)作物,在山地現(xiàn)代農(nóng)業(yè)發(fā)展實(shí)踐中具有重大作用。近年來(lái),由丁香假單胞菌煙草致病變種(Pseudomonas syringae pv. tabaci)引起的野火病(tobacco wildfire)已演變成煙草種植的主要細(xì)菌性葉斑類病害,常與角斑病、赤星病、棒孢霉葉斑病等混合發(fā)生,造成煙葉生產(chǎn)的嚴(yán)重經(jīng)濟(jì)損失(張廣民等,2002;陳燾等,2018)。

    植物病害發(fā)生在很大程度上與植物微生態(tài)環(huán)境失衡密不可分,葉際微生物是葉際微生態(tài)環(huán)境中的重要組成成分。研究表明,葉際微生物群落結(jié)構(gòu)的穩(wěn)定程度與植物葉片健康具有較強(qiáng)的相關(guān)性(Lindow amp; Brandl, 2003; Shakir et al., 2021)。隨著高通量測(cè)序技術(shù)的發(fā)展與成熟,越來(lái)越多的學(xué)者應(yīng)用該技術(shù)研究植物病害與微生物群落間的關(guān)系(Li et al., 2022)。煙草作為以收獲葉片為主的重要經(jīng)濟(jì)作物,其真菌性葉部病害發(fā)生與葉際微生物群落間的關(guān)系研究已有諸多報(bào)道。例如,黃宇等(2020)對(duì)煙草白粉病研究表明,健康煙葉(健葉)葉際真菌群落豐富度及多樣性高于病葉,病葉和健葉的葉際優(yōu)勢(shì)真菌均為高氏白粉菌屬(Golovinomyces)、曲霉屬(Aspergillus)和鏈格孢屬(Alternaria)且病斑組織的豐度高于健康組織;煙草赤星病葉際微生物群落以子囊菌門和變形菌門為主,鏈格孢屬和假單胞菌屬分別為真菌和細(xì)菌的優(yōu)勢(shì)屬,但健葉真菌和細(xì)菌群落多樣性及豐富度高于病葉(劉亭亭等,2021;Dai et al., 2022);Sun等(2023)報(bào)道了感染靶斑病的煙株葉片病原菌的豐度隨病級(jí)增加而顯著上升,葉際真菌群落多樣性和豐富度均低于健葉,細(xì)菌群落多樣性高于健葉而豐富度低于健葉。由此可見,由真菌引起的煙草葉部病害其葉際微生物群落的變化與發(fā)病程度存在密切關(guān)聯(lián)。

    煙草野火病是一種細(xì)菌性病害,20世紀(jì)40年代末期在我國(guó)云南煙區(qū)零星發(fā)生。隨著煙草栽培面積逐漸擴(kuò)大,野火病常與其他煙草葉部病害混合發(fā)生,病害流行年份造成重大損失。目前,針對(duì)煙草野火病發(fā)生的影響因素研究主要集中于氣候、溫濕度等環(huán)境因子,而對(duì)于葉際微生物的群落結(jié)構(gòu)及多樣性尚缺乏深入的認(rèn)知。程淺(2020)探討了煙草葉際細(xì)菌群落結(jié)構(gòu)與野火病發(fā)生的關(guān)系,但缺少真菌群落結(jié)構(gòu)的相關(guān)研究。Si等(2023)研究了野火病對(duì)雪茄煙的影響,表征感染葉片與健葉之間的微生物群和化學(xué)因子差異,但未對(duì)不同病害程度的野火病煙葉做出進(jìn)一步分析。煙草野火病不但在田間為害煙葉,而且在煙葉采收后至烘干前仍可繼續(xù)為害,導(dǎo)致葉片病斑繼續(xù)擴(kuò)大變成枯焦?fàn)疃ズ婵緝r(jià)值。因此,亟須闡明其病害流行程度與葉際微生物群落的動(dòng)態(tài)變化。鑒于此,本研究以不同發(fā)病程度的煙草野火病煙葉為研究對(duì)象,采用高通量測(cè)序技術(shù),擬探討以下問題:(1)不同程度發(fā)病煙葉病斑組織與健康組織中微生物群落結(jié)構(gòu)與多樣性差異;(2)不同發(fā)病程度病斑組織微生物群落的變化;(3)不同發(fā)病程度煙葉組織微生物的功能變化。旨在為認(rèn)識(shí)煙草野火病的演變規(guī)律提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 樣品采集

    2019年6月,在種植品種為云煙87的貴州省施秉縣煙葉產(chǎn)區(qū)(107°51′36″—108°28′51″ E、26°49′47″—27°21′15″ N)采集移栽60 d后發(fā)生煙草野火病的煙葉樣品,根據(jù)病斑面積占比分為輕度(5%~10%)、中度(11%~20%)和重度(21%~50%)葉片(圖1),依次記為A、B、C。用無(wú)菌剪刀剪取樣品葉片的病斑組織(b)和健康組織(j),分別混合均勻置于50 mL無(wú)菌離心管中,每個(gè)樣品3管,每管5 g。樣品采集后放入低溫保存箱,并迅速帶回實(shí)驗(yàn)室,置于-80 ℃冰箱保存、備用。

    1.2 基因組DNA提取、PCR擴(kuò)增及高通量測(cè)序

    參照Fast DNA SPIN試劑盒(MP Biomedicals)使用說(shuō)明,提取葉際微生物基因組總DNA,經(jīng)超微量分光光度計(jì)NanoDrop 2000(Thermo Fisher Scientific)檢測(cè)質(zhì)量與濃度后,用無(wú)菌水將DNA稀釋至1 ng·μL-1備用。

    細(xì)菌16S rRNA基因和真菌轉(zhuǎn)錄間隔區(qū)ITS1區(qū)分別采用引物515F/806R(Yuan et al., 2020)和ITS5-1737F/ITS2-2043R(Wu et al., 2019)進(jìn)行擴(kuò)增。PCR擴(kuò)增體系及程序參照Apprill等(2005)。擴(kuò)增產(chǎn)物經(jīng)2%濃度瓊脂糖凝膠電泳檢測(cè)、GeneJET膠回收產(chǎn)物純化試劑盒(Thermo Fisher Scientific)純化后送至北京諾禾致源科技股份有限公司,利用Ion Plus Fragment Library試劑盒(Thermo Fisher Scientific)構(gòu)建文庫(kù)并在Ion S5TMXL平臺(tái)上進(jìn)行測(cè)序。

    1.3 數(shù)據(jù)質(zhì)控及群落結(jié)構(gòu)與多樣性分析

    下機(jī)數(shù)據(jù)經(jīng)Cutadapt(https://cutadapt.readthedocs.io/en/v1.9.1/)軟件剔除、過濾低質(zhì)量序列,以及UCHIME算法去除嵌合體序列,最終得到有效序列(clean reads)。使用UPARSE(http://www.drive5.com/usearch/)軟件在相似度≥97%的條件下對(duì)clean reads進(jìn)行OTU(operational taxonomic units)聚類,篩選每個(gè)OTU的代表性序列,并基于Mothur算法將OTU代表性序列與Unite(7.2)數(shù)據(jù)庫(kù)和Silva132數(shù)據(jù)庫(kù)比對(duì)、注釋物種分類(Edgar amp; Robert, 2013; Urmas et al., 2013)。利用R語(yǔ)言工具繪制物種群落柱狀圖、主成分分析(principal component analysis, PCA)圖和Venn分析圖等。利用QIIME(https://docs.qiime2.org)按樣本最小序列數(shù)進(jìn)行樣本序列抽平后,計(jì)算樣品微生物群落α多樣性指數(shù)(Shannon、Simpson、Chao1等);使用DPS軟件(http://www.dpsw.cn/download.html)進(jìn)行不同組別物種相對(duì)豐度的差異性分析。采用Tax4Fun程序語(yǔ)言及FUNGuild和KEGG數(shù)據(jù)庫(kù)分析真菌、細(xì)菌群落的功能類群(孫美麗等,2023)。

    2 結(jié)果與分析

    2.1 測(cè)序深度評(píng)估及數(shù)據(jù)質(zhì)控

    稀釋性曲線結(jié)果顯示,葉際細(xì)菌的OTU數(shù)在測(cè)序數(shù)據(jù)量為2 300時(shí)曲線趨于平坦(圖2:A),葉際真菌的OTU數(shù)在測(cè)序數(shù)據(jù)量為40 000時(shí)曲線趨于平坦(圖2:B),表明測(cè)序數(shù)據(jù)已包含所測(cè)樣品中絕大多數(shù)的物種信息,測(cè)序結(jié)果合理。

    原始測(cè)序數(shù)據(jù)質(zhì)控后,平均每組樣品的細(xì)菌和真菌分別獲得80 092條和70 147條高質(zhì)量序列(表1)。細(xì)菌與真菌原始測(cè)序數(shù)據(jù)的GenBank Accession Number分別為PRJNA1025832、PRJNA1025370。

    2.2 葉際細(xì)菌與真菌群落結(jié)構(gòu)

    2.2.1 門水平葉際群落結(jié)構(gòu) 健康與病斑組織的細(xì)菌群落包括變形菌門(Proteobacteria)、厚壁菌門(Firmicutes)和擬桿菌門(Bacteroidetes),其中病斑組織變形菌門的相對(duì)豐度高于健康組織(圖3:A)。輕度、中度和重度病葉病斑組織的優(yōu)勢(shì)細(xì)菌門為變形菌門(相對(duì)豐度分別為49.44%、72.39%和61.92%),其次為厚壁菌門(相對(duì)豐度分別為0.01%、1.54%和0.73%)和擬桿菌門(相對(duì)豐度分別為0.06%、0.19%和0.06%)。健康組織優(yōu)勢(shì)細(xì)菌為變形菌門(Aj、Bj和Cj的相對(duì)豐度分別為2.89%、8.76%和3.15%)。輕度病葉的病斑與健康組織的變形菌門相對(duì)豐度無(wú)顯著性差異,而中度和重度病葉病斑組織的變形菌門相對(duì)豐度顯著高于健康組織;厚壁菌門和擬桿菌門在6組樣品間無(wú)顯著性差異(表2)。

    健康和病斑組織的真菌群落主要為子囊菌門(Ascomycota)和擔(dān)子菌門(Basidiomycota)(圖3:B)。其相對(duì)豐度病斑組織高于健康組織,其中重度病葉病斑組織與健康組織間子囊菌門相對(duì)豐度存在顯著性差異。在輕度、中度和重度病葉病斑組織中,子囊菌門的相對(duì)豐度分別為10.25%、25.86%和32.22%,擔(dān)子菌門的相對(duì)豐度分別為4.52%、2.48%和4.35%。健康組織子囊菌門相對(duì)豐度隨病害程度的增加而減少(Aj、Bj和Cj的相對(duì)豐度分別為6.31%、5.14%和0.47%),而擔(dān)子菌門相對(duì)豐度隨病害程度的增加而增加(分別為0.19%、0.30%和0.47%)。

    2.2.2 屬水平葉際群落結(jié)構(gòu) 由圖4:A可知,健康與病斑組織的細(xì)菌群落在屬水平上為假單胞菌屬(Pseudomonas)和泛菌屬(Pantoea)。病斑組織Ab、Bb和Cb假單胞菌屬的相對(duì)豐度(分別為39.47%、59.83%和46.05%)與健康組織Aj、Bj和Cj假單胞菌屬的相對(duì)豐度(分別為0.57%、5.95%和1.10%)存在顯著差異。病斑組織泛菌屬的相對(duì)豐度隨病害程度增加而逐級(jí)上升,與健康組織無(wú)顯著性差異。

    煙葉健康和病斑組織的真菌群落優(yōu)勢(shì)屬包括鏈格孢屬(Alternaria)、小不整球殼屬(Plectosphaerella)和錐蓋傘屬(Conocybe),其病斑組織的相對(duì)豐度高于健康組織,但無(wú)顯著性差異(圖4:B),而莖點(diǎn)霉屬(Phoma)在健康組織的相對(duì)豐度高于病斑組織,并且在輕度病葉中二者存在顯著性差異。

    2.2.3 OTU屬水平葉際群落結(jié)構(gòu) 由圖5可知,病斑與健康組織樣品共有的細(xì)菌和真菌OTU數(shù)目分別為9和190,其中共有的細(xì)菌OTU分別屬于假單胞菌屬、泛菌屬、立克次氏體目Rickettsiales等,共有的真菌OTU分別屬于鏈格孢屬、錐蓋傘屬、小不整球殼屬、莖點(diǎn)霉屬等。病葉Bb和健葉Bj樣品細(xì)菌群落獨(dú)有的OTU數(shù)目較其他煙葉樣品高(圖5:A);病葉Bb樣品真菌群落獨(dú)有的OTU數(shù)目較其他煙葉組織高(圖5:B)。

    2.3 葉際微生物群落α多樣性

    由表3可知,煙葉健康與病斑組織的覆蓋度指數(shù)均大于0.997,表明測(cè)序結(jié)果包含絕大多數(shù)微生物群落。供試6組樣品細(xì)菌群落的Simpson指數(shù)、ACE指數(shù)和覆蓋度指數(shù)均無(wú)顯著性差異,僅輕度病斑組織的Shannon指數(shù)與健康組織存在顯著性差異,中度病斑組織的Chao1指數(shù)顯著高于輕度和重度病葉健康組織。病斑組織細(xì)菌群落物種多樣性與豐富度高于健康組織,Shannon指數(shù)和Simpson指數(shù)隨煙葉病害程度的增加呈先減后增的變化趨勢(shì),而Chao1指數(shù)和ACE指數(shù)則呈現(xiàn)先增后減的趨勢(shì);健康組織的細(xì)菌群落多樣性和豐富度指數(shù)隨病害程度的增加呈先增后減的趨勢(shì)。

    輕度和中度病斑組織真菌群落多樣性和豐富度均高于健康組織,但重度病斑組織真菌群落多樣性和豐富度低于其健康組織。Bb樣品的Shannon指數(shù)與Simpson指數(shù)顯著高于Cb、Aj、Bj樣品;6組樣品組織真菌群落的Chao1指數(shù)、ACE指數(shù)和覆蓋度指數(shù)均不存在顯著性差異。病斑組織的真菌群落多樣性和豐富度指數(shù)隨病害程度的增加呈先增后減的趨勢(shì);健康組織的Shannon指數(shù)和Simpson指數(shù)隨病害程度的增加而增加,Chao1指數(shù)和ACE指數(shù)則隨發(fā)病程度的增加而降低。

    2.4 葉際真菌與細(xì)菌群落功能預(yù)測(cè)

    基于KEGG數(shù)據(jù)庫(kù)對(duì)葉際細(xì)菌群落PICRUSt功能預(yù)測(cè)結(jié)果表明,在第一級(jí)功能層水平上,健康和病斑組織葉際細(xì)菌的主要功能為新陳代謝、遺傳信息處理、環(huán)境信息處理、細(xì)胞轉(zhuǎn)化、人類疾病以及有機(jī)系統(tǒng)6類生物代謝通路(圖6:A)。健康和病斑組織的葉際細(xì)菌基因功能大體相同且相對(duì)豐度無(wú)顯著性差異。

    基于FUNGuild功能預(yù)測(cè),健康與病斑組織葉際真菌生態(tài)功能類群相似。健康和病斑組織葉際真菌優(yōu)勢(shì)類群包括植物病原菌、動(dòng)物病原菌-內(nèi)生真菌-植物病原菌-木質(zhì)腐生真菌、未定義腐生真菌、木質(zhì)腐生真菌和植物病原菌-木質(zhì)腐生真菌。病斑組織的植物病原菌類群相對(duì)豐度高于健康組織,中度病葉病斑組織Bb相對(duì)豐度最高(6.23%),其次為Ab(5.27%)、Cb(2.50%),而健康組織中Bj組相對(duì)豐度最高(0.91%)。動(dòng)物病原菌-內(nèi)生真菌-植物病原菌-木質(zhì)腐生真菌類群的相對(duì)豐度在病斑組織Bb(6.51%)遠(yuǎn)高于其他病斑和健康組織。病斑組織的未定義腐生真菌、木質(zhì)腐生真菌類群相對(duì)豐度均高于健康組織,前者在病斑組織Cb中(4.66%)相對(duì)豐度最高,后者在Bb組(1.04%)相對(duì)豐度最高;健康組織的植物病原菌-木質(zhì)腐生真菌類群相對(duì)豐度高于病斑組織,健康組織Aj中(1.84%)相對(duì)豐度最高(圖6:B)。

    3 討論與結(jié)論

    豐富多樣的葉際微生物與宿主植物的健康生長(zhǎng)密切相關(guān)。葉際微生物能夠幫助宿主植物固氮、促進(jìn)其生長(zhǎng)、增強(qiáng)其抗逆和抗病能力,但其存在的病原微生物也可通過侵染葉片組織并在葉際定殖積累導(dǎo)致植物發(fā)?。ˋbdelfattah et al., 2015;楊寬等,2021)。本研究基于高通量測(cè)序技術(shù),從微觀層面解析了感染野火病煙葉健康與病斑組織葉際微生物的異同,并揭示了不同發(fā)病程度葉片葉際微生物群落結(jié)構(gòu)與多樣性。

    本研究發(fā)現(xiàn),感野火病煙葉的健康與病斑組織葉際微生物菌群結(jié)構(gòu)在物種相對(duì)豐度上存在較大差異。病斑和健康組織的葉際細(xì)菌優(yōu)勢(shì)門和屬分別為變形菌門和假單胞菌屬及泛菌屬等,皆為煙草葉際常見菌群,與Qin等(2019)對(duì)煙草野火病田間防治、孫美麗等(2023)對(duì)靶斑病煙葉葉際優(yōu)勢(shì)細(xì)菌和劉暢等(2020)對(duì)赤星病煙葉葉際優(yōu)勢(shì)細(xì)菌等研究報(bào)道一致。煙草野火病作為細(xì)菌性病害,常與真菌性病害協(xié)同發(fā)生,煙田中野火病與靶斑病、 赤星病等復(fù)合侵染造成煙株發(fā)病的案例屢見不鮮,因此其優(yōu)勢(shì)細(xì)菌類群或存在相似之處。假單胞菌屬作為野火病的病原,其相對(duì)豐度在病斑組織顯著高于健康組織,表明病原菌的大量定殖是導(dǎo)致野火病蔓延發(fā)生的主要原因。雖然泛菌屬在病斑和健康組織中的相對(duì)豐度無(wú)顯著性差異,但其豐度隨病害程度的增加而逐級(jí)上升,說(shuō)明泛菌屬與煙草野火病的發(fā)生存在正相關(guān)關(guān)系,可能在病害發(fā)生時(shí)具有潛在促進(jìn)作用。葉際真菌主要分布于鏈格孢屬、小不整球殼屬和錐蓋傘屬等,這與劉亭亭等(2022)、向立剛等(2020)研究結(jié)果類似。鏈格孢屬真菌的存在,表明煙葉后期存在赤星病與野火病混合發(fā)生的可能, 這從目前煙田病害發(fā)生規(guī)律上可見一斑。此外,鏈格孢屬、小不整球殼屬和錐蓋傘屬在病斑組織的富集豐度明顯高于健康組織,表明三者在病害發(fā)生時(shí)可能與病原菌存在協(xié)同促進(jìn)的作用,有助于病害發(fā)生。值得注意的是,莖點(diǎn)霉屬在健康組織的相對(duì)豐度高于病斑組織,并且在輕度病葉中二者存在顯著性差異,表明其在病害抑制中具有潛在作用。針對(duì)葉際微生物中變化顯著的優(yōu)勢(shì)菌群值得進(jìn)一步功能探索與研究。

    微生物之間高度復(fù)雜的關(guān)系使得局部某些病原菌或有益菌的相對(duì)豐度變化直至引起整個(gè)群落組成的改變,而植物病害發(fā)生伴隨著病原菌的入侵定殖,常常會(huì)改變微生物環(huán)境和群落結(jié)構(gòu),對(duì)微生物群落多樣性產(chǎn)生影響(劉宇星等,2022)。Zhang等(2018)發(fā)現(xiàn)南瓜白粉病病原菌Podosphaere的豐度隨病情加重而增加;Sun等(2023)對(duì)不同病害等級(jí)煙草靶斑病煙葉的研究表明,葉際病原菌的相對(duì)豐度隨著煙葉病級(jí)的增加而顯著上升;陳乾麗等(2019)對(duì)不同霉變程度煙葉霉?fàn)€病的研究結(jié)果發(fā)現(xiàn),病原菌Aspergillus在霉變重的煙葉葉片中的相對(duì)豐度較霉變輕的煙葉更高。綜上表明,病原菌豐度及群落特征常與煙葉發(fā)病程度存在一定相關(guān)性。本研究對(duì)選取的3類不同發(fā)病程度的野火病煙葉進(jìn)行微生物群落結(jié)構(gòu)分析發(fā)現(xiàn),在輕度、中度和重度病葉病斑組織中假單胞菌屬的相對(duì)豐度均大于39%,但在中度病葉病斑組織中其相對(duì)豐度最高,這與前人研究結(jié)果和預(yù)測(cè)有所不同。此外,煙葉健康與病斑組織的微生物多樣性指數(shù)隨著病害程度的增加呈先增后減或逐級(jí)遞減。造成這種現(xiàn)象可能的原因是在野火病發(fā)病初期,葉片組織完整性被破壞,引起環(huán)境中包括病原菌在內(nèi)的各類微生物涌入遭破壞的組織中大量定植,進(jìn)而造成病斑組織微生物群落多樣性和豐富度短暫增加,但隨著病斑面積擴(kuò)大,葉片有限的營(yíng)養(yǎng)物質(zhì)和生存空間不斷降低,微生物間的定殖和競(jìng)爭(zhēng)加劇,最終導(dǎo)致病原菌以及其他微生物數(shù)量降低。Huang等(2021)對(duì)不同病級(jí)的亞隔孢殼屬葉斑病煙葉的研究也有類似發(fā)現(xiàn),發(fā)病煙葉組織雖然在功能預(yù)測(cè)中顯示植物病原菌相對(duì)豐度隨病害增加而成比例上升,但葉際病原菌Didymella的相對(duì)豐度在病害后期反而有下降趨勢(shì)。這恰與本研究假單胞菌屬的相對(duì)豐度在重度病斑組織中表現(xiàn)出降低趨勢(shì)一致。葉際微生物群的結(jié)構(gòu)變化因受到環(huán)境等因素的影響往往難以準(zhǔn)確預(yù)測(cè),但其趨勢(shì)大體可循。一般而言,病原菌和與其相關(guān)聯(lián)的微生物菌群大量定植于葉片組織造成的葉際微生態(tài)平衡被破壞是導(dǎo)致煙葉病害發(fā)生的微觀機(jī)制之一。由此可知,煙草病害的防治可圍繞抑制破壞病原菌和相關(guān)菌群的定植進(jìn)而達(dá)到生物防控的目的。劉天波等(2021)利用拮抗菌防治煙草野火病得到較好效果,并發(fā)現(xiàn)假單胞菌屬、泛菌屬等煙草葉際微生物所占比例發(fā)生顯著變化,芽孢桿菌屬作為有益菌相比對(duì)照增加3.9倍。

    葉際微生物在定植積累過程中常發(fā)揮自身特定的生理功能,與其群落結(jié)構(gòu)及組成密切相關(guān)(高爽等,2016)。本研究中功能預(yù)測(cè)發(fā)現(xiàn),健康和病斑組織的細(xì)菌類群基因功能信息及相對(duì)豐度較為趨同,以代謝、遺傳信息處理和環(huán)境信息處理3種功能為主,該結(jié)果與煙草角斑病、靶斑病的葉際細(xì)菌功能預(yù)測(cè)結(jié)果較為一致(吳小軍等,2023;孫美麗等,2023),其原因可能是野火病菌在以上葉際病害中也常被作為優(yōu)勢(shì)類群,因而功能預(yù)測(cè)存在相似之處。本研究發(fā)現(xiàn),病斑組織細(xì)菌多樣性與豐富度高于健康組織,推測(cè)野火病侵染煙葉后,假單胞菌屬、泛菌屬等菌群在葉際養(yǎng)分的利用上占據(jù)主要地位,在定植過程中獲取營(yíng)養(yǎng)物質(zhì)而大量繁殖,進(jìn)一步刺激和開發(fā)其代謝、遺傳和環(huán)境信息處理等功能,進(jìn)而演變?yōu)楦胁熑~葉際優(yōu)勢(shì)細(xì)菌屬。此外,本研究還發(fā)現(xiàn)在病斑與健康組織葉際真菌群落功能預(yù)測(cè)中,病原菌、內(nèi)生真菌及腐生真菌等占據(jù)主導(dǎo)地位,病斑組織的植物病原菌類群相對(duì)豐度高于健康組織,這與Huang等(2021)報(bào)道一致,進(jìn)一步說(shuō)明了植物病原菌和內(nèi)生真菌、腐生真菌為煙草葉際優(yōu)勢(shì)真菌類群,植物病原菌是導(dǎo)致葉際病害的重要原因。

    參考文獻(xiàn):

    ABDELFATTAH A, DESTRI NMG, CACCIOLA SO, et al., 2015.Metabarcoding analysis of fungal diversity in the phyllosphere and carposphere of olive (Olea europaea)" [J]. PLoS ONE, 10(7): e0131069.

    APPRILL A, MCNALLY S, PARSONS R, et al., 2005. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton" [J]. Aquat Microb Ecol, 75(2): 129-137.

    CHEN QL, LI Z, WANG HC, et al., 2019. Fungal composition and diversity of tobacco phyllosphere from cured tobacco leaves" [J]. Acta Microbiol Sin, 59(12): 2401-2409." [陳乾麗, 李忠, 汪漢成, 等, 2019. 烤后不同霉變程度煙葉際真菌群落組成與多樣性分析 [J]. 微生物學(xué)報(bào), 59(12): 2401-2409.]

    CHENG Q, 2020. Study on the relationship between the diversity of tobacco bacterial community and the occurrence of tobacco wildfire and its regulation effect [D]. Chongqing: Southwest University: 23-24." [程淺, 2020. 煙草葉際細(xì)菌群落多樣性與煙草野火病發(fā)生的關(guān)系及調(diào)控效應(yīng)研究 [D]. 重慶: 西南大學(xué): 23-24.]

    CHEN T, ZHOU W, LI HG, et al., 2018. Research progress in disease occurrence and integrated control of tobacco wild fire" [J]. Genomics Appl Biol, 37(1): 469-476." [陳燾, 周瑋, 李宏光, 等, 2018. 煙草野火病的發(fā)生及綜合防治研究進(jìn)展 [J]. 基因組學(xué)與應(yīng)用生物學(xué), 37(1): 469-476.]

    DAI YF, WU XM, WANG HC, et al., 2022. Spatio-temporal variation in the phyllospheric microbial biodiversity of Alternaria alternata-infected tobacco foliage" [J]. Front Microbiol, 13: 920109.

    EDGAR H, ROBERT C, 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads" [J]. Nat Methods, 10(10): 996-998.

    GAO S, LIU XC, DONG Z, et al., 2016. Advance of phyllosphere microorganisms and their interaction with the outside environment" [J]. Plant Sci J, 34(4): 654-661." [高爽, 劉笑塵, 董錚, 等, 2016. 葉際微生物及其與外界互作的研究進(jìn)展 [J]. 植物科學(xué)學(xué)報(bào), 34(4): 654-661.]

    HUANG Y, WANG HC, CAI LT, et al., 2021. Phyllospheric microbial composition and diversity of the tobacco leaves infected by Didymella segeticola" [J]. Front Microbiol, 12: 699699.

    HUANG Y, WANG HC, CHEN QL, et al., 2020. Community structure and diversity of fungi in leaves of tobacco plant infected with powdery mildew" [J]. Guizhou Agric Sci, 48(5): 54-58." [黃宇, 汪漢成, 陳乾麗, 等, 2020. 感染白粉病煙株葉片真菌的群落結(jié)構(gòu)與多樣性 [J]. 貴州農(nóng)業(yè)科學(xué), 48(5): 54-58.]

    LINDOW SE, BRANDL MT, 2003.Microbiology of the phyllosphere" [J]. Appl Environ Microbiol, 69(4): 1875-1883.

    LI PD, ZHU ZR, ZHANG Y, et al., 2022. The phyllosphere microbiome shifts toward combating melanose pathogen" [J]. Microbiome, 10(1): 56-56.

    LIU C, WANG HC, XIE HL, et al., 2020. Biodiversity analysis of phyllosphere bacterial genus from tobacco leaves infected by brown spot disease" [J]. Tob Sci Technol, 53(2): 8-14." [劉暢, 汪漢成, 謝紅煉, 等, 2020. 感染赤星病煙草葉際細(xì)菌的多樣性分析 [J]. 煙草科技, 53(2): 8-14.]

    LIU TB, TENG K, ZHOU XP, et al., 2021. Effects of antagonistic bacteria on tobacco wild fire and responses of phyllosphere microbiota" [J]. Microbiol Chin, 48(8): 2643-2652." [劉天波, 滕凱, 周向平, 等, 2021. 拮抗菌群對(duì)煙草野火病的防治效果及葉際微生物群落多樣性的影響 [J]. 微生物學(xué)通報(bào), 48(8): 2643-2652.]

    LIU TT, ZENG YT, WANG HC, et al., 2021. Metabolic and microbial community structure analysis in the phyllosphere of tobacco leaves with different maturity during brown spot occurring season" [J]. Chin Tob Sci, 42(6): 22-29." [劉亭亭, 曾隕濤, 汪漢成, 等, 2021. 赤星病發(fā)生期不同成熟度煙葉葉際微生物代謝與群落結(jié)構(gòu) [J]. 中國(guó)煙草科學(xué), 42(6): 22-29.]

    LIU TT, WANG HC, SUN ML, et al., 2022. Effects of Bordeaux mixture on the community structure and metabolic function of tobacco phyllosphere microorganisms" [J]. Chin J Pestic Sci, 24(6): 1446-1455." [劉亭亭, 汪漢成, 孫美麗, 等, 2022. 波爾多液對(duì)煙草葉際微生物群落結(jié)構(gòu)與代謝功能的影響 [J]. 農(nóng)藥學(xué)學(xué)報(bào), 24(6): 1446-1455.]

    LIU YX, DONG CB, SHAO QY, et al., 2022. Advances on phyllosphere microorganisms and their association with plant health" [J]. J Microbiol, 42(2): 88-98." [劉宇星, 董醇波, 邵秋雨, 等, 2022. 葉際微生物與植物健康研究進(jìn)展 [J]. 微生物學(xué)雜志, 42(2): 88-98.]

    QIN C, TAO JM, LIU TB, et al., 2019. Responses of phyllosphere microbiota and plant health to application of two different biocontrol agents" [J]. AMB Express, 9(1): 42.

    SHAKIR S, ZAIDI SS, VRIES DTF, et al., 2021. Plant genetic networks shaping phyllosphere microbial community" [J]. Trends Genet, 37(4): 306-316.

    SI H, CUI B, LIU F, et al., 2023. Microbial community and chemical composition of cigar tobacco (Nicotiana tabacum L.) leaves altered by tobacco wildfire disease [J]. Plant Direct, 7(12): e551.

    SUN ML, SHI CH, HUANG Y, et al., 2023. Effect of disease severity on the structure and diversity of the phyllosphere microbial community in tobacco" [J]. Front Microbiol, 13(4): 1081576.

    SUN ML, SHI CH, XIAO BQ, et al., 2023. Composition and diversity of phyllospheric microbial community in tobacco leaves infected by tobacco target spot disease" [J]. Tob Sci Technol, 56(4): 1-9." [孫美麗, 史彩華, 肖本青, 等, 2023. 煙草靶斑病葉際微生物群落結(jié)構(gòu)與多樣性分析 [J]. 煙草科技, 56(4): 1-9.]

    URMAS K, NILSSON RH, KESSY A, et al., 2013. Towards a unified paradigm for sequence-based identification of fungi" [J]. Mol Ecol, 22(21): 5271-5277.

    WU XJ, WANG HC, SUN ML, et al., 2023. Composition and diversity of phyllosphere microbial community on tobacco leaves infected with angular leaf spot" [J]. Tob Sci Technol, 56(10): 1-10." [吳小軍, 汪漢成, 孫美麗, 等, 2023. 煙草角斑病葉際微生物群落結(jié)構(gòu)與多樣性分析 [J]. 煙草科技, 56(10): 1-10.]

    WU YY, GENG JK, XU S, et al., 2019. Isolation and identification of intestinal fungi in Spodoptera frugiperda larvae of Chongqing area" [J]. J Southwest Univ (Nat Sci Ed), 41(9): 1-8.

    XIANG LG, WANG HC, ZHENG P, et al., 2020. Analysis of endophytic fungi and phyllosphere fungi of flue cured tobacco leaves with brown spot disease" [J]. Acta Tab Sin, 26(4): 93-100.nbsp; [向立剛, 汪漢成, 鄭蘋, 等, 2020. 赤星病烤后煙葉內(nèi)生及葉際真菌分析 [J]. 中國(guó)煙草學(xué)報(bào), 26(4): 93-100.]

    YANG K, WANG HL, YE KH, et al., 2021. Advances in research on phyllosphere microorganisms and their interaction with plants" [J]. J Yunnan Agric Univ (Nat Sci), 36(1): 155-164." [楊寬, 王慧玲, 葉坤浩, 等, 2021. 葉際微生物及與植物互作的研究進(jìn)展 [J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)), 36(1): 155-164.]

    YUAN Z, ZHANG CY, PENG XX, et al., 2020. Intestinal microbiota characteristics of mice treated with Folium senna decoction gavage combined with restraint and tail pinch stress" [J]. 3 Biotech, 10(4): 1-11.

    ZHANG GM, LU JH, KAN GF, et al., 2002. Recent development in the study of tobacco wildfire disease" [J]. Acta Tabac Sin, 8(2): 34-37." [張廣民, 呂軍鴻, 闞光鋒, 等, 2002. 煙草野火病研究概況 [J]. 中國(guó)煙草學(xué)報(bào), 8(2): 34-37.]

    ZHANG Z, LUO LY, TAN XQ, et al., 2018. Pumpkin powdery mildew disease severity influences the fungal diversity of the phyllosphere" [J]. PeerJ, 6(2): e4559.

    (責(zé)任編輯 李 莉 王登惠)

    猜你喜歡
    高通量測(cè)序病葉群落結(jié)構(gòu)
    大興安嶺被松櫟柱銹菌侵染的蒙古櫟葉上真菌種群的分布與多樣性
    川明參輪作對(duì)煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    海南省紅樹林區(qū)浮游動(dòng)物多樣性的初步研究
    海南省紅樹林區(qū)大型底棲動(dòng)物多樣性的初步研究
    柳樹病葉的傅里葉變換紅外光譜研究
    濟(jì)南流域輪蟲群落結(jié)構(gòu)及多樣性的研究
    黃河蘭州市區(qū)段浮游動(dòng)物群落結(jié)構(gòu)調(diào)查及水質(zhì)初步評(píng)價(jià)
    亚洲欧美日韩另类电影网站 | 日本黄大片高清| 日韩一本色道免费dvd| 少妇丰满av| 亚洲国产精品一区三区| 免费观看的影片在线观看| 亚洲av综合色区一区| 午夜视频国产福利| 亚洲精品aⅴ在线观看| 七月丁香在线播放| 黄色视频在线播放观看不卡| tube8黄色片| 久久精品国产亚洲av涩爱| 国产乱人视频| 午夜激情久久久久久久| 国产片特级美女逼逼视频| 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 国产精品av视频在线免费观看| 最近中文字幕高清免费大全6| 人妻系列 视频| 成人午夜精彩视频在线观看| 午夜日本视频在线| 精品久久国产蜜桃| 少妇 在线观看| 免费人成在线观看视频色| 老司机影院毛片| 99热国产这里只有精品6| 免费大片18禁| 亚洲国产成人一精品久久久| 女人久久www免费人成看片| 一区二区三区精品91| 国产男女内射视频| 秋霞在线观看毛片| 欧美日本视频| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂 | 99热国产这里只有精品6| 校园人妻丝袜中文字幕| 亚洲人成网站在线播| 国产一区二区三区综合在线观看 | 九九爱精品视频在线观看| 免费看光身美女| 街头女战士在线观看网站| 国产高清国产精品国产三级 | 久久精品国产自在天天线| 免费大片18禁| 纯流量卡能插随身wifi吗| 国产成人一区二区在线| 国产极品天堂在线| 久久综合国产亚洲精品| 久久久久久久大尺度免费视频| 在线精品无人区一区二区三 | 国产乱人偷精品视频| 亚洲色图av天堂| 精品国产乱码久久久久久小说| 亚洲欧美日韩卡通动漫| 伊人久久国产一区二区| 大话2 男鬼变身卡| 国产精品三级大全| 国产中年淑女户外野战色| 欧美变态另类bdsm刘玥| av播播在线观看一区| 美女xxoo啪啪120秒动态图| 国产视频内射| 国产精品.久久久| 看非洲黑人一级黄片| 亚洲精品久久午夜乱码| 亚洲成人一二三区av| 99re6热这里在线精品视频| 国产精品成人在线| 精品少妇久久久久久888优播| 成人二区视频| 最近中文字幕2019免费版| 菩萨蛮人人尽说江南好唐韦庄| 少妇猛男粗大的猛烈进出视频| 男人爽女人下面视频在线观看| 麻豆成人午夜福利视频| 国产精品一二三区在线看| 在线 av 中文字幕| 如何舔出高潮| 老司机影院毛片| 七月丁香在线播放| 国产成人aa在线观看| 又黄又爽又刺激的免费视频.| 成人高潮视频无遮挡免费网站| 五月开心婷婷网| 欧美日韩国产mv在线观看视频 | 欧美成人午夜免费资源| 国产精品人妻久久久久久| 日本-黄色视频高清免费观看| 国产极品天堂在线| 亚洲av成人精品一二三区| 日韩大片免费观看网站| 久久ye,这里只有精品| 久久午夜福利片| 国产视频内射| 国产成人a区在线观看| 亚洲欧美日韩东京热| 91久久精品电影网| 亚洲国产av新网站| 我要看日韩黄色一级片| 亚洲四区av| 香蕉精品网在线| 久久精品熟女亚洲av麻豆精品| 久久久久国产网址| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 在线观看人妻少妇| 免费高清在线观看视频在线观看| 毛片一级片免费看久久久久| 精品亚洲成a人片在线观看 | 亚洲欧美日韩另类电影网站 | 黑人猛操日本美女一级片| 国产中年淑女户外野战色| 久久久久久久久久成人| 一级毛片 在线播放| 五月天丁香电影| 亚洲av中文字字幕乱码综合| 亚洲精品,欧美精品| 校园人妻丝袜中文字幕| 国产成人a区在线观看| 国产永久视频网站| 亚洲aⅴ乱码一区二区在线播放| 欧美高清性xxxxhd video| 国产精品久久久久久久电影| 日本猛色少妇xxxxx猛交久久| 在线免费观看不下载黄p国产| 久久av网站| 3wmmmm亚洲av在线观看| 国产精品久久久久久精品电影小说 | 男女无遮挡免费网站观看| 精品一区在线观看国产| 成年av动漫网址| 高清午夜精品一区二区三区| 激情 狠狠 欧美| 欧美日韩亚洲高清精品| 两个人的视频大全免费| 国产综合精华液| 国产无遮挡羞羞视频在线观看| 欧美激情极品国产一区二区三区 | 久久久久久久国产电影| 中国美白少妇内射xxxbb| kizo精华| 免费观看性生交大片5| 国产中年淑女户外野战色| 男人添女人高潮全过程视频| 日本黄大片高清| 我的女老师完整版在线观看| 91精品国产国语对白视频| 国产精品久久久久久精品古装| 欧美变态另类bdsm刘玥| 国产爱豆传媒在线观看| 一级a做视频免费观看| 性色avwww在线观看| 亚洲图色成人| 精品酒店卫生间| 中文欧美无线码| 99九九线精品视频在线观看视频| 秋霞在线观看毛片| 少妇的逼水好多| 久久久精品免费免费高清| 久久久久久九九精品二区国产| 一级毛片电影观看| 黄色视频在线播放观看不卡| 久久精品熟女亚洲av麻豆精品| 成人亚洲欧美一区二区av| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 美女高潮的动态| 成人美女网站在线观看视频| 欧美+日韩+精品| av播播在线观看一区| 麻豆国产97在线/欧美| 国产 精品1| 最近最新中文字幕大全电影3| 日韩在线高清观看一区二区三区| 一级毛片黄色毛片免费观看视频| 久久精品久久久久久噜噜老黄| 国产综合精华液| 国产乱来视频区| 97超碰精品成人国产| av网站免费在线观看视频| 国产一区亚洲一区在线观看| 一区二区三区免费毛片| 久久精品国产亚洲av天美| 女人久久www免费人成看片| 精华霜和精华液先用哪个| 天美传媒精品一区二区| 亚洲精品成人av观看孕妇| 国产成人a区在线观看| 久久99精品国语久久久| 男人爽女人下面视频在线观看| 精品少妇黑人巨大在线播放| 精品久久久久久久久av| 久久精品国产a三级三级三级| av网站免费在线观看视频| 久久97久久精品| 亚洲精品国产av蜜桃| 国产亚洲av片在线观看秒播厂| 亚洲av免费高清在线观看| 麻豆乱淫一区二区| 日韩电影二区| 最近中文字幕2019免费版| 精品久久久久久久久亚洲| 综合色丁香网| 亚洲精品成人av观看孕妇| 国产伦精品一区二区三区四那| 午夜福利在线观看免费完整高清在| 一区二区三区精品91| 深爱激情五月婷婷| 久久人人爽人人片av| 人妻 亚洲 视频| 中文字幕av成人在线电影| 精品久久国产蜜桃| 国产成人aa在线观看| 久久精品熟女亚洲av麻豆精品| 国产免费福利视频在线观看| 黄色怎么调成土黄色| tube8黄色片| 国产精品福利在线免费观看| 国产片特级美女逼逼视频| 国产成人一区二区在线| 一边亲一边摸免费视频| 欧美极品一区二区三区四区| 国内少妇人妻偷人精品xxx网站| 91精品伊人久久大香线蕉| 小蜜桃在线观看免费完整版高清| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 国产精品不卡视频一区二区| 欧美精品国产亚洲| 高清在线视频一区二区三区| 亚洲成人中文字幕在线播放| 久久精品夜色国产| av天堂中文字幕网| 亚洲欧洲日产国产| 国产一区有黄有色的免费视频| 国产真实伦视频高清在线观看| 成年人午夜在线观看视频| 日本黄色日本黄色录像| 男女边吃奶边做爰视频| 一级爰片在线观看| 最新中文字幕久久久久| 国产精品一区www在线观看| 亚洲av二区三区四区| 在线观看国产h片| 国产色爽女视频免费观看| 亚洲精品日韩av片在线观看| 国产人妻一区二区三区在| 波野结衣二区三区在线| 久久99热这里只频精品6学生| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 高清毛片免费看| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 日韩成人av中文字幕在线观看| 联通29元200g的流量卡| 国产精品一及| 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 成人二区视频| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 久久久久久久久久人人人人人人| 91精品国产国语对白视频| 亚洲欧美成人精品一区二区| 国产亚洲午夜精品一区二区久久| 免费少妇av软件| 亚洲美女视频黄频| 国产v大片淫在线免费观看| 又爽又黄a免费视频| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 成人免费观看视频高清| 国精品久久久久久国模美| 亚洲av中文字字幕乱码综合| 一边亲一边摸免费视频| 国产免费又黄又爽又色| 十分钟在线观看高清视频www | 少妇熟女欧美另类| 国内少妇人妻偷人精品xxx网站| 午夜老司机福利剧场| 99国产精品免费福利视频| 少妇的逼好多水| 欧美成人精品欧美一级黄| 在线观看av片永久免费下载| 久久久久精品性色| 成人毛片60女人毛片免费| 美女内射精品一级片tv| 国产白丝娇喘喷水9色精品| 激情 狠狠 欧美| 免费人成在线观看视频色| 高清日韩中文字幕在线| 又爽又黄a免费视频| 婷婷色综合大香蕉| 嘟嘟电影网在线观看| 国产黄色免费在线视频| 99热6这里只有精品| 男女免费视频国产| 能在线免费看毛片的网站| 免费看光身美女| 舔av片在线| 黄色欧美视频在线观看| 精品午夜福利在线看| 极品少妇高潮喷水抽搐| 国产伦精品一区二区三区视频9| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕| 久久精品夜色国产| 亚洲精品自拍成人| 久久久久国产精品人妻一区二区| 我的老师免费观看完整版| 97热精品久久久久久| 黄色一级大片看看| 久久国产亚洲av麻豆专区| 国产精品99久久久久久久久| 一级二级三级毛片免费看| 成人一区二区视频在线观看| 少妇被粗大猛烈的视频| 亚洲va在线va天堂va国产| 搡女人真爽免费视频火全软件| 日本爱情动作片www.在线观看| 国产精品99久久99久久久不卡 | 卡戴珊不雅视频在线播放| 精品久久久噜噜| 亚洲欧美精品专区久久| 午夜老司机福利剧场| 18+在线观看网站| 亚洲欧美成人综合另类久久久| 老司机影院成人| 国产成人一区二区在线| 国产成人免费观看mmmm| 国产视频内射| 黑丝袜美女国产一区| 国产视频首页在线观看| 亚洲精品乱码久久久v下载方式| 成人毛片60女人毛片免费| 国产黄片视频在线免费观看| 老女人水多毛片| 国产伦精品一区二区三区视频9| 中文字幕制服av| 一区在线观看完整版| 妹子高潮喷水视频| 亚洲欧美精品专区久久| 国产无遮挡羞羞视频在线观看| 婷婷色麻豆天堂久久| 人妻少妇偷人精品九色| 亚洲一区二区三区欧美精品| 人妻系列 视频| 91午夜精品亚洲一区二区三区| 国产片特级美女逼逼视频| av天堂中文字幕网| 一级二级三级毛片免费看| 亚洲欧美中文字幕日韩二区| av国产精品久久久久影院| 高清毛片免费看| 22中文网久久字幕| 国产成人免费无遮挡视频| 美女中出高潮动态图| 观看av在线不卡| 麻豆成人午夜福利视频| 国产在线免费精品| 麻豆国产97在线/欧美| 狠狠精品人妻久久久久久综合| 91久久精品电影网| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 高清日韩中文字幕在线| 亚洲av综合色区一区| 人妻夜夜爽99麻豆av| 久久国产亚洲av麻豆专区| 成人午夜精彩视频在线观看| 精品久久久久久电影网| 亚洲精品第二区| 免费av不卡在线播放| 青春草亚洲视频在线观看| 免费看不卡的av| 在线观看国产h片| 色5月婷婷丁香| 欧美日韩视频高清一区二区三区二| 伦精品一区二区三区| 观看美女的网站| 亚洲欧美日韩东京热| 亚洲,一卡二卡三卡| 性高湖久久久久久久久免费观看| 99久久综合免费| 波野结衣二区三区在线| 免费观看的影片在线观看| 天堂中文最新版在线下载| 大香蕉97超碰在线| 日韩欧美 国产精品| 欧美日韩国产mv在线观看视频 | 日本免费在线观看一区| 国产毛片在线视频| 99热国产这里只有精品6| 欧美精品亚洲一区二区| 亚洲成人av在线免费| 一级毛片久久久久久久久女| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 欧美日韩综合久久久久久| 综合色丁香网| a级一级毛片免费在线观看| 亚洲综合色惰| 欧美3d第一页| 在线观看av片永久免费下载| 亚洲va在线va天堂va国产| 国产一区二区三区综合在线观看 | 嫩草影院入口| 免费观看在线日韩| 日韩 亚洲 欧美在线| 亚洲激情五月婷婷啪啪| 色视频www国产| 国产 一区精品| 久久久久久久国产电影| 丰满人妻一区二区三区视频av| 国产大屁股一区二区在线视频| xxx大片免费视频| 久久这里有精品视频免费| 噜噜噜噜噜久久久久久91| 久久久久久久久久成人| 国产成人a∨麻豆精品| 国产有黄有色有爽视频| 下体分泌物呈黄色| 欧美亚洲 丝袜 人妻 在线| 热re99久久精品国产66热6| 精品人妻视频免费看| 免费黄网站久久成人精品| 少妇人妻久久综合中文| 中文字幕久久专区| 成人漫画全彩无遮挡| 看免费成人av毛片| 一区二区三区乱码不卡18| 丝袜喷水一区| www.色视频.com| 欧美一级a爱片免费观看看| 欧美国产精品一级二级三级 | 亚洲av中文av极速乱| 高清欧美精品videossex| av天堂中文字幕网| 国产高清有码在线观看视频| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 男人爽女人下面视频在线观看| 波野结衣二区三区在线| 天美传媒精品一区二区| 欧美97在线视频| 色婷婷av一区二区三区视频| 夫妻午夜视频| 亚洲精品中文字幕在线视频 | 最近最新中文字幕免费大全7| 亚洲av在线观看美女高潮| 五月伊人婷婷丁香| 热re99久久精品国产66热6| 国产亚洲5aaaaa淫片| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩东京热| 欧美zozozo另类| 成人免费观看视频高清| 高清不卡的av网站| 国产成人一区二区在线| 91久久精品国产一区二区三区| 成年女人在线观看亚洲视频| 在线观看三级黄色| 久久精品久久久久久久性| 国产高清不卡午夜福利| 色视频在线一区二区三区| 精品亚洲乱码少妇综合久久| 国产av精品麻豆| 天天躁日日操中文字幕| 高清视频免费观看一区二区| 国产亚洲91精品色在线| 伊人久久国产一区二区| 国产 一区精品| a级毛色黄片| 欧美精品一区二区大全| 日本一二三区视频观看| .国产精品久久| av在线播放精品| 内射极品少妇av片p| 舔av片在线| 亚洲精品一二三| 久久影院123| 亚洲精品一二三| 国产高清国产精品国产三级 | 国产伦精品一区二区三区视频9| 极品教师在线视频| 哪个播放器可以免费观看大片| 欧美3d第一页| 男人舔奶头视频| 久久久久精品性色| 身体一侧抽搐| 日本一二三区视频观看| 国产av码专区亚洲av| 亚洲高清免费不卡视频| 日本黄色日本黄色录像| 综合色丁香网| 色5月婷婷丁香| 久久久久人妻精品一区果冻| 成年女人在线观看亚洲视频| 国产亚洲欧美精品永久| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 亚洲图色成人| 久久精品国产自在天天线| 亚洲av综合色区一区| 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的 | 天堂中文最新版在线下载| 国产亚洲精品久久久com| 久久人妻熟女aⅴ| 在现免费观看毛片| 国产免费福利视频在线观看| 五月天丁香电影| 国产女主播在线喷水免费视频网站| 国产片特级美女逼逼视频| 国产高潮美女av| 丝袜脚勾引网站| 边亲边吃奶的免费视频| av女优亚洲男人天堂| 国产女主播在线喷水免费视频网站| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| 热re99久久精品国产66热6| 99久久中文字幕三级久久日本| kizo精华| 一本久久精品| 亚洲真实伦在线观看| 精品亚洲成a人片在线观看 | 国产高清不卡午夜福利| 少妇熟女欧美另类| 高清日韩中文字幕在线| 亚洲性久久影院| 高清欧美精品videossex| 亚洲欧美精品专区久久| 亚洲av成人精品一区久久| 国产综合精华液| 天堂中文最新版在线下载| 干丝袜人妻中文字幕| 国产欧美另类精品又又久久亚洲欧美| 91狼人影院| 久久国内精品自在自线图片| 亚洲熟女精品中文字幕| 国产色婷婷99| 一二三四中文在线观看免费高清| 国产免费福利视频在线观看| 特大巨黑吊av在线直播| 色哟哟·www| 亚洲经典国产精华液单| 51国产日韩欧美| 亚洲欧洲国产日韩| 国产精品人妻久久久久久| 最后的刺客免费高清国语| 国产乱人视频| 97在线人人人人妻| 偷拍熟女少妇极品色| 欧美xxⅹ黑人| 在线亚洲精品国产二区图片欧美 | 26uuu在线亚洲综合色| 男女无遮挡免费网站观看| 一区二区三区四区激情视频| 亚洲精品aⅴ在线观看| 欧美高清性xxxxhd video| 亚洲精品国产色婷婷电影| 熟女电影av网| 色视频www国产| 在线观看三级黄色| 在线精品无人区一区二区三 | 成人毛片60女人毛片免费| 国产精品久久久久久精品古装| 亚洲av福利一区| 内射极品少妇av片p| 亚洲国产色片| 九九久久精品国产亚洲av麻豆| 欧美老熟妇乱子伦牲交| 2021少妇久久久久久久久久久| 欧美日韩精品成人综合77777| 日韩欧美一区视频在线观看 | 高清在线视频一区二区三区| 我要看日韩黄色一级片| 亚洲欧美精品自产自拍| 国产亚洲一区二区精品| 久久人人爽人人爽人人片va| 国产精品一区二区三区四区免费观看| 免费观看的影片在线观看| 成人黄色视频免费在线看| 国产淫语在线视频| 亚洲高清免费不卡视频| 国产精品福利在线免费观看| 青春草视频在线免费观看| 美女内射精品一级片tv| 街头女战士在线观看网站| 久久精品国产亚洲av天美| 精品午夜福利在线看| 久久亚洲国产成人精品v| 91久久精品国产一区二区成人| 大陆偷拍与自拍| 丝袜脚勾引网站| 亚洲美女黄色视频免费看| 欧美日韩亚洲高清精品| 99久久中文字幕三级久久日本| 亚洲av中文字字幕乱码综合| 久久久久久久精品精品| 国产有黄有色有爽视频| 国产美女午夜福利| 国产午夜精品久久久久久一区二区三区| 中文字幕制服av| 国产精品一及| 国产伦精品一区二区三区四那| 校园人妻丝袜中文字幕| 熟女电影av网|