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

    暹羅芽孢桿菌揮發(fā)性物質(zhì)拮抗下煙草赤星病菌轉(zhuǎn)錄組分析

    2023-06-11 20:58:57李藝池梁金昌韋承建王東坤王衛(wèi)民陳勇華宋光龍胡希好程德杰任廣偉王曉強(qiáng)
    中國(guó)煙草科學(xué) 2023年2期

    李藝池 梁金昌 韋承建 王東坤 王衛(wèi)民 陳勇華 宋光龍 胡希好 程德杰 任廣偉 王曉強(qiáng)

    摘? 要:暹羅芽孢桿菌(Bacillus siamensis)LZ88是一株對(duì)多種病原真菌具有較強(qiáng)拮抗能力的生防細(xì)菌,產(chǎn)生的兩種揮發(fā)性有機(jī)化合物(VOCs)2-甲基丁酸和3-甲基丁酸對(duì)鏈格孢菌具有顯著的抑制作用。為進(jìn)一步研究其拮抗機(jī)制,本研究通過(guò)2-甲基丁酸和3-甲基丁酸混合物與鏈格孢菌的平板對(duì)扣試驗(yàn)觀察其拮抗作用,采用RNA-Seq技術(shù)對(duì)揮發(fā)性混合物拮抗下鏈格孢菌的轉(zhuǎn)錄組進(jìn)行研究,并對(duì)部分差異表達(dá)顯著的基因進(jìn)行了qRT-PCR驗(yàn)證。結(jié)果表明,2種揮發(fā)性有機(jī)化合物混合處理能顯著抑制鏈格孢菌的生長(zhǎng),菌絲表現(xiàn)出皺縮、塌陷的現(xiàn)象;通過(guò)差異表達(dá)分析共獲得3461個(gè)差異表達(dá)基因(DEGs),其中1714個(gè)上調(diào)表達(dá),1747個(gè)下調(diào)表達(dá)。對(duì)獲得的DEGs進(jìn)行功能注釋及富集分析發(fā)現(xiàn),揮發(fā)性有機(jī)化合物主要是通過(guò)影響鏈格孢菌的甘油磷脂代謝、MAPK信號(hào)通路以及聚酮類化合物的合成來(lái)抑制鏈格孢菌的生長(zhǎng)發(fā)育。本研究揭示了2種揮發(fā)性化合物對(duì)鏈格孢菌的抑菌機(jī)制,為生防菌代謝產(chǎn)物的拮抗機(jī)制分析提供了新的思路。

    關(guān)鍵詞:暹羅芽孢桿菌;鏈格孢菌;揮發(fā)性有機(jī)化合物;轉(zhuǎn)錄組分析

    中圖分類號(hào):S435.72? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼:A? ? ? ? ? ? ? ? ? ? 文章編號(hào):1007-5119(2023)02-0027-08

    Abstract: Bacillus siamensis LZ88 is a biocontrol bacterium with strong antagonistic ability against a variety of pathogenic fungi. Previous study indicated that volatile organic compounds (VOCs), 2-methylbutyric acid and 3-methylbutyric acid, produced by B. siamensis LZ88 showed significant inhibitory effects on Alternaria alternata. In the present study, the mechanism of 2-methylbutyric acid and 3-methylbutyric acid on A. alternate were further explored. Additionally, RNA-Seq was used to study the transcriptome of two volatile mixtures on A. alternata, and some significantly differentially expressed genes were verified by qRT-PCR. The results showed that the mixed treatment of the two VOCs could significantly inhibit the growth of A. alternata, and the hyphae were shrinkage and collapse. Through differential expression analysis, a total of 3461 differentially expressed genes (DEGs) were obtained, of which 1714 genes were up-regulated and 1747 genes were down-regulated. Functional annotation and enrichment analysis of the obtained DEGs revealed that various biological functions of A. alternata were affected by VOCs, including glycerophospholipid metabolism, MAPK signaling pathway and the synthesis of polyketides. The results give new insight into the inhibition mechanism of biocontrol bacteria metabolites products to A. alternata.

    Keywords: Bacillus siamensis; Alternaria alternata; volatile organic compounds; transcriptome analysis

    煙草赤星病是在煙葉成熟后期由半知菌亞門鏈格孢屬真菌(Alternaria spp.)引起的一種重要的葉部真菌性病害,具有潛育期短、流行速度快等特點(diǎn)[1],是我國(guó)煙草生產(chǎn)上威脅最大的病害之一,直接影響煙葉的產(chǎn)量和品質(zhì)。目前,煙草生產(chǎn)上對(duì)赤星病的防治主要以化學(xué)防治為主,但是長(zhǎng)期大量使用化學(xué)藥劑會(huì)導(dǎo)致病原菌產(chǎn)生抗藥性,并造成煙葉農(nóng)藥殘留問(wèn)題,影響煙葉品質(zhì)[2-3]。生物防治能有效避免或減少對(duì)生物和環(huán)境造成的不利影響,在煙草病蟲(chóng)害的防控中具有極大的應(yīng)用前景。

    芽孢桿菌是一種常見(jiàn)的生防細(xì)菌,并被廣泛應(yīng)用于植物病害的防控。芽孢桿菌可以產(chǎn)生多種抗菌物質(zhì),包括非揮發(fā)性化合物和揮發(fā)性化合物。其中,揮發(fā)性有機(jī)化合物(VOCs)主要包括酯類、烷烴類、醇類、萜類等[4],可以抑制病原菌的生長(zhǎng),或誘導(dǎo)植物產(chǎn)生系統(tǒng)抗性[5-6]。目前,利用芽孢桿菌揮發(fā)性化合物抑菌病原菌的研究主要以描述病原菌的菌絲、孢子變化和細(xì)胞形態(tài)為主。例如,特基拉芽孢桿菌XK29產(chǎn)生的揮發(fā)物2-甲基丁酸顯著抑制甘薯長(zhǎng)喙殼菌菌絲生長(zhǎng)和孢子萌發(fā),降低其產(chǎn)孢能力,導(dǎo)致菌絲折疊彎曲并形成不連續(xù)的空腔等[7],但鮮有關(guān)于揮發(fā)性化合物抑制病原真菌生長(zhǎng)的分子機(jī)制的報(bào)道。

    本課題組前期研究發(fā)現(xiàn)暹羅芽孢桿菌(Bacillus siamensis)LZ88對(duì)多種病原真菌都具有顯著的拮抗活性,通過(guò)對(duì)其揮發(fā)性有機(jī)化合物進(jìn)行研究,發(fā)現(xiàn)2-甲基丁酸和3-甲基丁酸具有抗真菌活性,能夠抑制鏈格孢菌的生長(zhǎng)[8],但其拮抗機(jī)制尚不清楚。因此,本研究對(duì)2-甲基丁酸和3-甲基丁酸處理后的鏈格孢菌進(jìn)行轉(zhuǎn)錄組分析,探究?jī)煞N揮發(fā)性化合物對(duì)鏈格孢菌的抑菌機(jī)制,以期為利用暹羅芽孢桿菌揮發(fā)性化合物防控鏈格孢菌提供理論依據(jù),并為相關(guān)病害的防控提供新思路。

    1? 材料與方法

    1.1? 材料

    供試病原菌:鏈格孢菌(A. alternata)由本實(shí)驗(yàn)室從煙草赤星病發(fā)病煙葉上分離鑒定并保存。將鏈格孢菌在PDA培養(yǎng)基上28 ℃活化培養(yǎng)。

    供試培養(yǎng)基:馬鈴薯葡萄糖瓊脂培養(yǎng)基(PDA)和馬鈴薯葡萄糖肉湯培養(yǎng)基(PDB)購(gòu)自海博生物技術(shù)有限公司。

    供試試劑:2-甲基丁酸(純度:98%;CAS:116-53-0)和3-甲基丁酸(純度≥99.0%;CAS:503-74-2)購(gòu)自北京索萊寶生物科技有限公司;Trizol試劑購(gòu)自Invitrogen公司;Yeasen Hifair?Ⅲ 1st Strand cDNA Synthesis SuperMix for qPCR(gDNA digester plus)購(gòu)自翌圣生物科技有限公司;SparkJade 2×SYBR Green qPCR Mix(With ROX)購(gòu)自思科捷生物技術(shù)有限公司。

    1.2? 方法

    1.2.1? 2-甲基丁酸和3-甲基丁酸混合物對(duì)鏈格孢菌生長(zhǎng)的影響? 根據(jù)前期的試驗(yàn)結(jié)果[8],提前將2-甲基丁酸和3-甲基丁酸純品用無(wú)水乙醇分別配置成83.10和104.19 mg/mL的溶液。處理組分別吸取10 μL配置好的2-甲基丁酸和3-甲基丁酸溶液添加至直徑6 mm的濾紙片上,放在上層PDA培養(yǎng)基中心,取直徑約5 mm的新鮮鏈格孢菌菌餅接種在下層PDA培養(yǎng)基中心,將兩個(gè)培養(yǎng)皿對(duì)扣培養(yǎng)。對(duì)照組只在下層PDA培養(yǎng)基中心接種鏈格孢菌菌餅。處理組和對(duì)照組各重復(fù)3次,于28 ℃培養(yǎng)7 d后,測(cè)定抑菌效果。收集各組菌絲體,用掃描電子顯微鏡(SEM)拍照觀察不同處理組菌絲形態(tài)。

    1.2.2? 轉(zhuǎn)錄組樣品準(zhǔn)備? 在150 mL無(wú)菌PDB液體培養(yǎng)基中接入4個(gè)鏈格孢菌菌餅,180 r/min、28 ℃振蕩培養(yǎng)2 d后,處理組加入2-甲基丁酸(終濃度0.066 mg/mL)和3-甲基丁酸(終濃度0.058 mg/mL)各15 μL,對(duì)照組不做任何處理。對(duì)照組和處理組各3個(gè)重復(fù),180 r/min、28 ℃繼續(xù)振蕩培養(yǎng)4 d,用無(wú)菌濾膜過(guò)濾出菌絲,再用無(wú)菌水清洗3次,抽濾去掉水分后用液氮速凍,置于–80 ℃保存,用于提取鏈格孢菌的總RNA。

    1.2.3? RNA提取與測(cè)序文庫(kù)構(gòu)建? 鏈格孢菌總RNA的提取使用Trizol法。分別提取處理組和對(duì)照組的RNA,并通過(guò)Nanodrop 2000檢測(cè)總RNA的濃度和純度,瓊脂糖凝膠電泳檢測(cè)RNA的完整性,Agilent 2100測(cè)定RIN值。將總RNA送至上海美吉生物醫(yī)藥科技有限公司進(jìn)行測(cè)序(Illumina HiSeq 2500測(cè)序平臺(tái))。轉(zhuǎn)錄組數(shù)據(jù)上傳至NCBI(https://www.ncbi.nlm.nih.gov/)數(shù)據(jù)庫(kù),登錄號(hào)為PRJNA900390。

    1.2.4? 測(cè)序數(shù)據(jù)處理? 對(duì)過(guò)濾得到的高質(zhì)量clean reads進(jìn)行拼接組裝,獲得unigenes庫(kù)。通過(guò)Blast將篩選到的unigene序列與NR、Swiss-prot、Pfam、GO、COG和KEGG等公共數(shù)據(jù)庫(kù)進(jìn)行比對(duì),獲得注釋信息。利用FPKM(Fragments Per Kilobases per

    Millionreads)公式計(jì)算unigene的表達(dá)量,DESeq2軟件篩選差異表達(dá)基因(Differentially expressed genes,DEGs),以差異倍數(shù)[|log2(Fold change)|]>1和錯(cuò)誤發(fā)現(xiàn)率(false discovery rate, FDR)<0.05作為篩選的標(biāo)準(zhǔn)。通過(guò)Goatools和KOBAS對(duì)篩選出的差異表達(dá)基因分別進(jìn)行GO功能富集分析和KEGG通路分析。

    1.2.5? qRT-PCR驗(yàn)證? 將鏈格孢菌的總RNA利用Yeasen HifairⅢ 1st Strand cDNA Synthesis SuperMix for qPCR(gDNA digester plus)反轉(zhuǎn)錄成cDNA。以GAPDH為內(nèi)參基因,選取8 個(gè)與鏈格孢菌生長(zhǎng)有關(guān)的DEGs進(jìn)行qRT-PCR反應(yīng),以驗(yàn)證轉(zhuǎn)錄組測(cè)序結(jié)果的準(zhǔn)確性。引物序列見(jiàn)表1,反應(yīng)體系參考SparkJade 2×SYBR Green qPCR Mix(With ROX)試劑盒,利用ABI SetpOne? Real-Time PCR System進(jìn)行PCR擴(kuò)增,具體程序如下:94 ℃ 2 min;94 ℃ 5 s,60 ℃ 30 s,共40個(gè)循環(huán)。結(jié)果根據(jù)2?ΔΔCT法進(jìn)行數(shù)據(jù)分析。

    2? 結(jié)? 果

    2.1? 2-甲基丁酸和3-甲基丁酸混合物對(duì)鏈格孢菌生長(zhǎng)的影響

    如圖1所示,通過(guò)SEM觀察到對(duì)照組菌絲茂密,粗細(xì)較均勻,表面光滑,處理組的菌絲體損傷嚴(yán)重,菌絲粗細(xì)不均勻,出現(xiàn)明顯的皺縮、扭曲等現(xiàn)象。此外,與對(duì)照組相比,處理組樣品中孢子變形和塌陷,未附著分生孢子的菌絲出現(xiàn)破裂、收縮和扭曲,表明2-甲基丁酸和3-甲基丁酸的混合物抑制了鏈格孢菌絲的生長(zhǎng)。

    2.2? 測(cè)序結(jié)果的評(píng)估

    根據(jù)VOCs對(duì)鏈格孢菌生長(zhǎng)的影響效果,選取對(duì)照組和處理組各3個(gè)樣本進(jìn)行轉(zhuǎn)錄組測(cè)序。如表2所示,獲得的raw reads序列數(shù)量在23 477 010~55 234 710之間,過(guò)濾掉低質(zhì)量的序列reads后,獲得的clean reads序列的數(shù)量在23 476 806~55 234 282之間;滿足Q30標(biāo)準(zhǔn)的堿基數(shù)均達(dá)到95%以上,堿基錯(cuò)誤率均值約為0.0232%,表明測(cè)序質(zhì)量良好,可以為后續(xù)的分析提供可靠的原始數(shù)據(jù)。根據(jù)PCA分析結(jié)果(圖2),PC1解釋78.13%的變量,對(duì)照組和處理組分別聚成一簇,對(duì)照組和處理組間具有明顯差異。

    2.3? 基因表達(dá)差異

    與對(duì)照組相比,經(jīng)VOCs處理后的鏈格孢菌中有3461個(gè)差異表達(dá)基因(DEGs),其中上調(diào)表達(dá)的基因有1714個(gè),下調(diào)表達(dá)的基因有1747個(gè)(圖3)。對(duì)差異表達(dá)基因進(jìn)行COG(Cluster of Orthologous Groups of proteins)注釋,共有3255個(gè)DEGs注釋到COG數(shù)據(jù)庫(kù),包含1635個(gè)上調(diào)表達(dá)基因,1620個(gè)下調(diào)表達(dá)基因。在上調(diào)表達(dá)基因中翻譯、核糖體結(jié)構(gòu)與生物合成(J),翻譯后修飾、蛋白質(zhì)折疊、分子伴侶(O),氨基酸運(yùn)輸與代謝(E),細(xì)胞內(nèi)運(yùn)輸、分泌與囊泡運(yùn)輸(U)這4類最多,占到78.72%;而細(xì)胞內(nèi)運(yùn)輸、分泌與囊泡運(yùn)輸(U),碳水化合物運(yùn)輸與代謝(G),翻譯后修飾、蛋白質(zhì)折疊、分子伴侶(O)以及轉(zhuǎn)錄(K)這4類在下調(diào)基因中所占比例較高,達(dá)到85.31%(圖4)。

    2.4? 富集分析

    2.4.1? 差異表達(dá)基因GO功能富集分析? 對(duì)下調(diào)的DEGs進(jìn)行GO富集分析發(fā)現(xiàn)(圖5),細(xì)胞組分中細(xì)胞結(jié)構(gòu)體、膜組成部分和膜固有成分富集基因數(shù)較多,而生物過(guò)程中脂質(zhì)代謝過(guò)程富集最顯著。對(duì)上調(diào)的DEGs進(jìn)行GO富集分析發(fā)現(xiàn)細(xì)胞氨基酸生物合成過(guò)程、核糖核蛋白復(fù)合體組裝和有機(jī)酸分解代謝過(guò)程3個(gè)功能最顯著富集。

    2.4.2? 差異表達(dá)基因的KEGG富集分析? 對(duì)1714個(gè)上調(diào)和1747個(gè)下調(diào)的DEGs分別進(jìn)行KEGG富集分析(圖6-7),發(fā)現(xiàn)有937個(gè)上調(diào)的DEGs定位到112條KEGG通路中。顯著富集的通路有核糖體,纈氨酸、亮氨酸和異亮氨酸降解,RNA轉(zhuǎn)運(yùn),其中富集基因數(shù)目最多的是核糖體通路。另外有391個(gè)下調(diào)的DEGs定位到97條KEGG通路中,顯著富集的有甘油磷脂代謝、醚脂代謝、MAPK信號(hào)通路,富集基因數(shù)目最多的是甘油磷脂代謝通路。

    2.5? 與鏈格孢菌生長(zhǎng)相關(guān)代謝通路的差異表達(dá)基因分析

    結(jié)合COG功能注釋、GO富集以及KEGG富集,從所有的差異表達(dá)基因中篩選出與鏈格孢菌生長(zhǎng)相關(guān)的差異表達(dá)基因(p<0.01),下調(diào)的基因主要參與甘油磷脂代謝、MAPK信號(hào)通路及聚酮類物質(zhì)的合成,上調(diào)的基因主要參與谷胱甘肽代謝、淀粉和蔗糖代謝等過(guò)程(表3)。

    2.6? qRT-PCR驗(yàn)證

    為了檢驗(yàn)轉(zhuǎn)錄組測(cè)序結(jié)果的準(zhǔn)確性,從篩選出的差異表達(dá)基因中隨機(jī)挑選8個(gè)差異顯著的基因進(jìn)行qRT-PCR驗(yàn)證,其中包括5個(gè)下調(diào)表達(dá)基因和3個(gè)上調(diào)表達(dá)基因。結(jié)果表明,所選基因的表達(dá)水平與轉(zhuǎn)錄組結(jié)果趨勢(shì)一致(圖8),說(shuō)明轉(zhuǎn)錄組測(cè)序的結(jié)果是可靠的。

    3? 討? 論

    微生物在其生長(zhǎng)代謝過(guò)程中會(huì)產(chǎn)生多種非揮發(fā)性和揮發(fā)性次級(jí)代謝產(chǎn)物[9]。在自然環(huán)境狀態(tài)下,揮發(fā)性有機(jī)化合物是以多種化合物混合的形式產(chǎn)生。因此,本研究通過(guò)模擬自然環(huán)境狀態(tài),將暹羅芽孢桿菌LZ88[8]產(chǎn)生的揮發(fā)性有機(jī)化合物2-甲基丁酸和3-甲基丁酸混合后處理鏈格孢菌,利用轉(zhuǎn)錄組學(xué)進(jìn)一步分析了2-甲基丁酸和3-甲基丁酸混合物對(duì)鏈格孢菌的抑制作用機(jī)制,發(fā)現(xiàn)VOCs通過(guò)調(diào)控鏈格孢菌的甘油磷脂代謝途徑、MAPK信號(hào)通路以及聚酮類物質(zhì)的合成等抑制菌絲的生長(zhǎng)繁殖。

    甘油磷脂是真核生物中生物膜的重要組成部分,參與細(xì)胞膜對(duì)蛋白質(zhì)的識(shí)別和信號(hào)傳導(dǎo),與菌絲生長(zhǎng)發(fā)育密切相關(guān)。線粒體膜由磷脂雙分子層組成,參與調(diào)節(jié)線粒體的形態(tài)及線粒體功能發(fā)揮的過(guò)程[10],為生命提供能量[11]。在本研究中甘油磷脂代謝途徑中編碼磷脂酰絲氨酸合酶(AFT_927247)、磷脂酰絲氨酸脫羧酶(AFT_189401)、磷脂酰乙醇胺N-甲基轉(zhuǎn)移酶(AFT_789367)、FabD/溶血磷脂酶類蛋白(AFT_1009916)、溶血磷脂酶類蛋白(AFT_1061421)的基因顯著下調(diào),線粒體內(nèi)外膜的磷脂生物合成受阻。研究發(fā)現(xiàn)釀酒酵母中,PE(磷脂酰乙醇胺)和PS(磷脂酰絲氨酸)的生物合成喪失會(huì)導(dǎo)致線粒體缺陷,從而使小突變體的形成增加[12];白色念珠菌磷脂酰絲氨酸(PS)合酶突變體cho1Δ/Δ缺乏PS,磷脂酰乙醇胺降低,在一定程度上影響了線粒體呼吸功能,突變體表現(xiàn)出絲狀生長(zhǎng)缺陷[13]。因此,本研究中VOCs可能通過(guò)抑制磷脂的合成,使線粒體內(nèi)外膜不能行使正常的功能,進(jìn)而影響鏈格孢菌菌絲的正常生長(zhǎng)。

    MAPK信號(hào)通路是調(diào)控許多病原真菌生長(zhǎng)、發(fā)育和致病性胞外信號(hào)轉(zhuǎn)導(dǎo)的重要途徑[14]。在本研究中MAPK信號(hào)通路中細(xì)胞壁完整性途徑(cell wall intergrity, CWI)的基因Bck1(AFT_927542)和Fks2(AFT_1005247)顯著下調(diào),表明VOCs處理干擾了鏈格孢菌細(xì)胞壁完整性途徑中的級(jí)聯(lián)組分,阻礙了細(xì)胞壁的合成,對(duì)菌絲細(xì)胞壁完整性造成影響。此外,MAPK信號(hào)通路中高滲透壓途徑(High osmolarity pathway,HOG)相關(guān)的Sln1(AFT_946700)、Ssk1(AFT_186722)和Sko1顯著下調(diào),表明VOCs處理影響菌絲對(duì)外界壓力做出適應(yīng)性反應(yīng)。這與張楠等[15]在膠孢炭疽菌中發(fā)現(xiàn)的HOG途徑上游感受器Sho1同源基因CgSho1缺失,造成突變菌株?duì)I養(yǎng)生長(zhǎng)緩慢、產(chǎn)孢量下降、菌絲稀疏等特征是一致的。因此,鏈格孢菌經(jīng)VOCs處理后可能通過(guò)抑制MAPK信號(hào)通路從而影響其菌絲細(xì)胞壁的完整性并對(duì)壓力做出適應(yīng)性反應(yīng),造成菌絲生長(zhǎng)的缺陷。

    真菌聚酮類化合物(polyketides,PKs)主要由真菌聚酮合酶合成,被認(rèn)為是病原體毒素和致病因子的來(lái)源[16]。例如,來(lái)自于曲霉菌的聚酮類化合物aflatoxin B1是侵染花生、玉米、水稻、大豆等多種作物的毒力因子[17];鐮刀屬真菌產(chǎn)生的zearalenone是侵染玉米的重要毒力因子,主要通過(guò)聚酮化合物途徑合成[18]。鏈格孢菌中含有大量的聚酮合酶基因,包括PKSA(AFT_1057721)、PKSB(AFT_176278)、PKSC(AFT_1012953)、PKSD(AFT_1096114)、PKSH(AFT_976935)等,參與合成芳香族和復(fù)雜的還原型聚酮類化合物,包括alternaric acid、alternariol、macrosporin、zinnolide、 zinniol等[19-21]。本研究中,聚酮合酶基因顯著下調(diào),說(shuō)明VOCs處理后其表達(dá)受到抑制,聚酮類物質(zhì)的生物合成受阻,從而降低鏈格孢菌的毒性,保護(hù)煙株免受鏈格孢菌的脅迫。

    4? 結(jié)? 論

    本研究表明,暹羅芽孢桿菌產(chǎn)生的揮發(fā)性有機(jī)化合物2-甲基丁酸和3-甲基丁酸可以顯著抑制鏈格孢菌絲的生長(zhǎng);與對(duì)照組相比,經(jīng)VOCs處理后的鏈格孢菌中有3461個(gè)差異表達(dá)基因;結(jié)合COG功能注釋、GO富集以及KEGG富集分析發(fā)現(xiàn),鏈格孢菌生長(zhǎng)受阻與甘油酯代謝途徑、MAPK信號(hào)通路及聚酮類物質(zhì)的合成有關(guān),但具體機(jī)制還需進(jìn)一步深入研究。

    參考文獻(xiàn)

    [1]張中義,冷懷瓊,張志銘,等. 植物病原真菌學(xué)[M]. 成都:四川科學(xué)技術(shù)出版社,1988.

    ZHANG Z Y, LENG H Q, ZHANG Z M, et al. Plant pathogenic mycology[M]. Chengdu: Sichuan Science and Technology Press, 1988.

    [2]HOU Y, MA X, WAN W, et al. Comparative genomics of pathogens causing brown spot disease of tobacco: Alternaria longipes and Alternaria alternata[J]. Plos One, 2016, 11(5): e0155258.

    [3]JING C, ZHAO J, HAN X, et al. Essential oil of Syringa oblata Lindl. as a potential biocontrol agent against tobacco brown spot caused by Alternaria alternata[J]. Crop Protection, 2018, 104: 41-46.

    [4]孫力. 貝萊斯芽孢桿菌CanL-30揮發(fā)性有機(jī)物活性研究及菌劑研制[D]. 武漢:華中農(nóng)業(yè)大學(xué),2020.

    SUN L. Bioactivities of the volatile organic compounds of Bacillus velezensis CanL-30 and formulation of a wettable powder[D]. Wuhan: Huazhong Agricultural University, 2020.

    [5]ZAMIOUDIS C, KORTELAND J, VAN PELT J A, et al. Rhizobacterial volatiles and photosynthesis-related signals coordinate MYB72 expression in Arabidopsis roots during onset of induced systemic resistance and iron-deficiency responses[J]. Plant Journal, 2012, 84(2): 309-322.

    [6]GAO Z, ZHANG B, LIU H, et al. Identification of endophytic Bacillus velezensis ZSY-1 strain and antifungal activity of its volatile compounds against Alternaria solani and Botrytis cinerea[J]. Biological Control, 2017, 105: 27-39.

    [7]李騰杰,梁孫妍,郭健衡,等. 特基拉芽孢桿菌XK29揮發(fā)物2-甲基丁酸對(duì)甘薯長(zhǎng)喙殼菌的抑制作用研究[J]. 微生物學(xué)報(bào),2022,62(12):5018-5028.

    LI T J, LIANG S Y, GUO J H, et al. Inhibitory effects of volatile 2-methylbutyric acid produced by Bacillus tequilensis XK29 on Ceratocystis fimbriata[J]. Acta Microbiologica Sinica, 2022, 62(12): 5018-5028.

    [8]WANG D K, LI Y C, YUAN Y, et al. Identification of non-volatile and volatile organic compounds produced by Bacillus siamensis LZ88 and their antifungal activity against Alternaria alternata[J]. Biological Control, 2022, 169: 1049-9644.

    [9]MORATH S, HUNG R, BENNETT J W. Fungal volatile organic compounds: A review with emphasison their biotechnological potential[J]. Fungal Biology Reviews, 2012, 26(2-3): 73-83.

    [10]BAUSEWEIN T, MILLS D J, LANGER J D, et al. Cryo-EM structure of the TOM core complex from Neurospora Crassa[J]. Cell, 2017, 170(4): 693-700.

    [11]GUO R, ZONG S, WU M, et al. Architecture of human mitochondrial respiratory megacomplex Ⅰ2Ⅲ2Ⅳ2[J]. Cell, 2017, 170(6): 1247-1257.

    [12]BIRNER R, BURGERMEISTER M, SCHNEITER R, et al. Roles of phosphatidylethanolamine and of its several biosynthetic pathways in Saccharomyces cerevisiae[J]. Molecular Biology of the Cell, 2001, 12(4): 997-1007.

    [13]CHEN Y L, MONTEDONICO A E, KAUFFMAN S, et al. Phosphatidylserine synthase and phosphatidylserine decarboxylase are essential for cell wall integrity and virulence in Candida albicans[J]. Molecular Microbiology, 2010, 75(5): 1112-1132.

    [14]王笑連. 楊樹(shù)炭疽病菌CgMk1 MAPK及上游受體蛋白的功能分析[D]. 北京:北京林業(yè)大學(xué),2021.

    WANG X L. Functional analysis of CgMk1 MAPK and the upstream receptors in Colletotrichum gloeosporioides[D]. Beijing: Beijing Forestry University, 2021.

    [15]張楠,柳志強(qiáng),吳曼莉,等. 膠孢炭疽菌CgSho1基因的克隆與功能分析[J]. 植物病理學(xué)報(bào),2017,47(1):40-49.

    ZHANG N, LIU Z Q, WU M L, et al. Gene cloning and functional analysis of CgSho1 in Colletotrichum gloeosporioides[J]. Acta Phytopathologica Sinica, 2017, 47(1): 40-49.

    [16]MINAMI A, UGAI T, OZAKI T, et al. Predicting the chemical space of fungal polyketides by phylogeny-based bioinformatics analysis of polyketide synthase-nonribosomal peptide synthetase and its modification enzymes[J]. Scientific Reports, 2020, 10(1): 13556.

    [17]RUSHING B R, SELIM M I. Aflatoxin B1: A review on metabolism, toxicity, occurrence in food, occupational exposure, and detoxification methods[J]. Food and Chemical Toxicology, 2019, 124: 81-100.

    [18]RAI A, DAS M, TRIPATHI A. Occurrence and toxicity of a fusarium mycotoxin, zearalenone[J]. Critical Reviews in Food Science and Nutrition, 2020, 60(16): 2710-2729.

    [19]COURTIAL J, HELESBEUX J J, OUDART H, et al. Characterization of NRPS and PKS genes involved in the biosynthesis of SMs in Alternaria dauci including the phytotoxic polyketide aldaulactone[J]. Scientific Reports, 2022, 12(1): 8155.

    [20]FUJII I. Functional analysis of fungal polyketide biosynthesis genes[J]. The Journal of Antibiotics, 2010, 63: 207-218.

    [21]UGAI T, MINAMI A, GOMI K, et al. Genome mining approach for harnessing the cryptic gene cluster in Alternaria solani: production of PKS-NRPS hybrid metabolite, didymellamide B[J]. Tetrahedron Letters, 2016, 57: 2793-2796.

    国产1区2区3区精品| 90打野战视频偷拍视频| 日韩有码中文字幕| 脱女人内裤的视频| 纯流量卡能插随身wifi吗| 考比视频在线观看| 亚洲精品久久久久久婷婷小说| 欧美日韩亚洲国产一区二区在线观看 | 精品一区二区三卡| 一个人免费在线观看的高清视频 | 欧美性长视频在线观看| 美女午夜性视频免费| 精品国产一区二区三区久久久樱花| av福利片在线| 91精品伊人久久大香线蕉| 欧美亚洲日本最大视频资源| 一二三四社区在线视频社区8| 亚洲av国产av综合av卡| 制服人妻中文乱码| 不卡一级毛片| 19禁男女啪啪无遮挡网站| 亚洲一区二区三区欧美精品| 国产成人啪精品午夜网站| 91九色精品人成在线观看| 18禁观看日本| 国产精品自产拍在线观看55亚洲 | 91国产中文字幕| 国产男人的电影天堂91| 精品福利观看| av视频免费观看在线观看| 大型av网站在线播放| 少妇的丰满在线观看| 丰满少妇做爰视频| 午夜福利视频在线观看免费| 国产成人精品无人区| 少妇的丰满在线观看| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区| 熟女少妇亚洲综合色aaa.| 99国产精品一区二区三区| 丁香六月天网| 精品国产一区二区三区久久久樱花| 精品国产一区二区久久| 丝袜美足系列| 国产精品久久久av美女十八| 亚洲国产欧美一区二区综合| 青春草亚洲视频在线观看| 国产精品免费视频内射| 免费在线观看完整版高清| 超色免费av| 国产高清国产精品国产三级| 97精品久久久久久久久久精品| 亚洲国产看品久久| 伊人亚洲综合成人网| www.精华液| 美女国产高潮福利片在线看| a级片在线免费高清观看视频| 亚洲av国产av综合av卡| 制服人妻中文乱码| 悠悠久久av| 亚洲精品美女久久久久99蜜臀| 青草久久国产| 亚洲情色 制服丝袜| 另类精品久久| 人人妻人人澡人人爽人人夜夜| 99久久国产精品久久久| 精品久久久精品久久久| avwww免费| 中亚洲国语对白在线视频| 欧美国产精品一级二级三级| 国产精品一区二区精品视频观看| 99re6热这里在线精品视频| 两性夫妻黄色片| 亚洲情色 制服丝袜| 久久精品国产a三级三级三级| 人人妻人人澡人人爽人人夜夜| 精品久久久久久电影网| 日韩免费高清中文字幕av| 日本91视频免费播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产中文字幕在线视频| 天堂俺去俺来也www色官网| 欧美激情 高清一区二区三区| 欧美亚洲 丝袜 人妻 在线| 成人18禁高潮啪啪吃奶动态图| 亚洲精品在线美女| 捣出白浆h1v1| 18禁裸乳无遮挡动漫免费视频| 欧美国产精品va在线观看不卡| 天天影视国产精品| 久久天躁狠狠躁夜夜2o2o| 青春草亚洲视频在线观看| 999久久久精品免费观看国产| 十八禁网站免费在线| 欧美人与性动交α欧美精品济南到| av国产精品久久久久影院| 自线自在国产av| 一个人免费在线观看的高清视频 | 中文字幕制服av| 成人亚洲精品一区在线观看| 亚洲国产欧美在线一区| 亚洲精品第二区| 亚洲第一青青草原| 午夜老司机福利片| 欧美精品啪啪一区二区三区 | 中文字幕色久视频| 精品少妇一区二区三区视频日本电影| 十八禁网站免费在线| 日韩有码中文字幕| www日本在线高清视频| 日韩精品免费视频一区二区三区| 嫁个100分男人电影在线观看| 黄色a级毛片大全视频| 三上悠亚av全集在线观看| 性高湖久久久久久久久免费观看| 性少妇av在线| 国产精品1区2区在线观看. | 男人舔女人的私密视频| 亚洲av日韩在线播放| 法律面前人人平等表现在哪些方面 | 午夜福利免费观看在线| 亚洲成国产人片在线观看| a 毛片基地| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 亚洲欧美一区二区三区久久| 国产精品国产av在线观看| 国产精品成人在线| 无限看片的www在线观看| 久久精品熟女亚洲av麻豆精品| a在线观看视频网站| 亚洲精品久久午夜乱码| av在线app专区| 久久99一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 动漫黄色视频在线观看| 黑丝袜美女国产一区| 午夜激情久久久久久久| 老汉色∧v一级毛片| 91成年电影在线观看| 日韩中文字幕视频在线看片| 精品国产超薄肉色丝袜足j| av有码第一页| 各种免费的搞黄视频| 国产精品一区二区精品视频观看| 性高湖久久久久久久久免费观看| 大片电影免费在线观看免费| 丰满饥渴人妻一区二区三| www.av在线官网国产| av片东京热男人的天堂| 国产精品国产三级国产专区5o| 亚洲精品美女久久久久99蜜臀| 国产一区二区在线观看av| 老司机影院成人| kizo精华| 国产精品亚洲av一区麻豆| 水蜜桃什么品种好| 十八禁网站网址无遮挡| 欧美一级毛片孕妇| 美女视频免费永久观看网站| 国产福利在线免费观看视频| 亚洲精品av麻豆狂野| 久久久国产欧美日韩av| 亚洲五月婷婷丁香| 美女福利国产在线| 日韩欧美一区二区三区在线观看 | 亚洲久久久国产精品| 免费不卡黄色视频| 国产精品一二三区在线看| 国产精品一区二区在线观看99| 91麻豆av在线| 男人舔女人的私密视频| 久久精品成人免费网站| av电影中文网址| 国产又色又爽无遮挡免| 久久这里只有精品19| 欧美激情久久久久久爽电影 | 亚洲美女黄色视频免费看| netflix在线观看网站| av欧美777| 69精品国产乱码久久久| 精品福利观看| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 国产深夜福利视频在线观看| 一区二区三区精品91| 色播在线永久视频| 人妻一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 手机成人av网站| 成年人免费黄色播放视频| 首页视频小说图片口味搜索| 午夜视频精品福利| 99国产精品一区二区蜜桃av | 久久ye,这里只有精品| 韩国精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 国产免费一区二区三区四区乱码| 黄片小视频在线播放| 操美女的视频在线观看| 悠悠久久av| 亚洲av日韩在线播放| 免费在线观看日本一区| 国产淫语在线视频| 最新在线观看一区二区三区| 飞空精品影院首页| 老鸭窝网址在线观看| 窝窝影院91人妻| 日韩 欧美 亚洲 中文字幕| 国产熟女午夜一区二区三区| 国产精品一区二区精品视频观看| 日韩欧美免费精品| 黄片播放在线免费| 精品国产一区二区三区四区第35| 欧美日韩福利视频一区二区| 国产1区2区3区精品| 国产淫语在线视频| 女人精品久久久久毛片| 国产成人欧美| 亚洲精华国产精华精| 男男h啪啪无遮挡| 精品一区二区三卡| 日韩制服骚丝袜av| 亚洲熟女毛片儿| 丝袜美腿诱惑在线| 一区二区三区四区激情视频| 久久中文看片网| e午夜精品久久久久久久| 天堂俺去俺来也www色官网| 国内毛片毛片毛片毛片毛片| 亚洲精品久久午夜乱码| 久久国产精品男人的天堂亚洲| 亚洲熟女精品中文字幕| 国产精品久久久久久精品电影小说| 亚洲国产日韩一区二区| √禁漫天堂资源中文www| 好男人电影高清在线观看| 欧美精品人与动牲交sv欧美| 老司机深夜福利视频在线观看 | 黑人操中国人逼视频| 青春草视频在线免费观看| 人人妻人人澡人人看| 9热在线视频观看99| 日韩,欧美,国产一区二区三区| 两人在一起打扑克的视频| 国产日韩欧美视频二区| 久久人妻熟女aⅴ| 亚洲情色 制服丝袜| 国产又爽黄色视频| 欧美黄色淫秽网站| 天天躁夜夜躁狠狠躁躁| 最新的欧美精品一区二区| 我要看黄色一级片免费的| 亚洲精品一二三| 欧美黄色淫秽网站| 亚洲精品一区蜜桃| 在线观看免费午夜福利视频| 婷婷色av中文字幕| 汤姆久久久久久久影院中文字幕| 9热在线视频观看99| 老司机影院毛片| 精品第一国产精品| 午夜福利一区二区在线看| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 欧美日韩福利视频一区二区| 丰满饥渴人妻一区二区三| 69av精品久久久久久 | 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 一边摸一边抽搐一进一出视频| 精品久久久久久久毛片微露脸 | 日韩欧美一区视频在线观看| 各种免费的搞黄视频| 狠狠精品人妻久久久久久综合| 午夜福利在线观看吧| 一本综合久久免费| 国产一区二区三区在线臀色熟女 | 肉色欧美久久久久久久蜜桃| 大片免费播放器 马上看| 亚洲久久久国产精品| 国产亚洲午夜精品一区二区久久| 男人添女人高潮全过程视频| 少妇的丰满在线观看| 好男人电影高清在线观看| 久久精品亚洲av国产电影网| 亚洲成人免费电影在线观看| 亚洲黑人精品在线| 99国产精品一区二区三区| 久久精品人人爽人人爽视色| 久久久久精品国产欧美久久久 | 日本精品一区二区三区蜜桃| 国产极品粉嫩免费观看在线| 欧美 日韩 精品 国产| 亚洲av国产av综合av卡| 亚洲专区国产一区二区| av超薄肉色丝袜交足视频| 午夜福利影视在线免费观看| 俄罗斯特黄特色一大片| 丰满少妇做爰视频| 人人妻人人澡人人爽人人夜夜| 淫妇啪啪啪对白视频 | 老司机福利观看| 精品亚洲成a人片在线观看| 操出白浆在线播放| 亚洲国产成人一精品久久久| 久9热在线精品视频| 亚洲欧美激情在线| 亚洲av国产av综合av卡| 黑人巨大精品欧美一区二区mp4| 在线看a的网站| 久久精品国产综合久久久| 国产在线视频一区二区| 久久青草综合色| 久久中文看片网| 欧美少妇被猛烈插入视频| 国产日韩欧美在线精品| 777米奇影视久久| videos熟女内射| 最黄视频免费看| 亚洲成人免费电影在线观看| av超薄肉色丝袜交足视频| 老司机亚洲免费影院| 免费高清在线观看日韩| 国产成人精品久久二区二区免费| 91成年电影在线观看| 亚洲精品av麻豆狂野| 亚洲久久久国产精品| 国产日韩一区二区三区精品不卡| 国产伦理片在线播放av一区| 91老司机精品| 亚洲精品美女久久av网站| 淫妇啪啪啪对白视频 | 亚洲精品国产一区二区精华液| 波多野结衣av一区二区av| 亚洲国产av新网站| 午夜免费成人在线视频| 国产97色在线日韩免费| 少妇精品久久久久久久| 日韩一卡2卡3卡4卡2021年| 亚洲全国av大片| 久久香蕉激情| 久久人人97超碰香蕉20202| 日韩有码中文字幕| 丝瓜视频免费看黄片| 欧美日韩福利视频一区二区| www.自偷自拍.com| 色婷婷av一区二区三区视频| 精品亚洲乱码少妇综合久久| 操美女的视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av男天堂| 久久久国产欧美日韩av| 婷婷丁香在线五月| 午夜福利乱码中文字幕| 亚洲精品国产色婷婷电影| 一区福利在线观看| 2018国产大陆天天弄谢| 少妇被粗大的猛进出69影院| 欧美精品av麻豆av| 亚洲,欧美精品.| 国产人伦9x9x在线观看| 美女高潮喷水抽搐中文字幕| 黄色怎么调成土黄色| 亚洲精品美女久久久久99蜜臀| 黄色怎么调成土黄色| 老汉色av国产亚洲站长工具| 少妇猛男粗大的猛烈进出视频| 欧美激情极品国产一区二区三区| 亚洲成人手机| 香蕉丝袜av| 又黄又粗又硬又大视频| cao死你这个sao货| 亚洲精品成人av观看孕妇| 久久久久国内视频| 丰满饥渴人妻一区二区三| 欧美精品一区二区免费开放| 久久久久久人人人人人| 69精品国产乱码久久久| av线在线观看网站| e午夜精品久久久久久久| 丁香六月天网| 午夜激情久久久久久久| 午夜影院在线不卡| 午夜成年电影在线免费观看| a级片在线免费高清观看视频| 亚洲国产精品一区三区| 老汉色av国产亚洲站长工具| 国产成人精品久久二区二区免费| 中文字幕人妻丝袜制服| 午夜激情av网站| 久久精品成人免费网站| 亚洲精品中文字幕在线视频| 久9热在线精品视频| 在线亚洲精品国产二区图片欧美| 精品乱码久久久久久99久播| 国精品久久久久久国模美| 国产精品.久久久| 99久久综合免费| 久久国产精品大桥未久av| netflix在线观看网站| 欧美日韩亚洲综合一区二区三区_| 人妻人人澡人人爽人人| 欧美激情高清一区二区三区| 午夜福利视频精品| 99国产精品免费福利视频| 国精品久久久久久国模美| 久久中文看片网| 91麻豆av在线| 男女无遮挡免费网站观看| 欧美国产精品一级二级三级| 国产精品二区激情视频| 男男h啪啪无遮挡| 俄罗斯特黄特色一大片| 亚洲一区中文字幕在线| a在线观看视频网站| 亚洲精品国产区一区二| 热99re8久久精品国产| 久久精品国产a三级三级三级| 女人精品久久久久毛片| 久久亚洲精品不卡| 淫妇啪啪啪对白视频 | 亚洲国产精品一区三区| 美女高潮到喷水免费观看| 王馨瑶露胸无遮挡在线观看| 满18在线观看网站| 欧美成人午夜精品| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲少妇的诱惑av| 国产激情久久老熟女| 亚洲成人国产一区在线观看| 日日爽夜夜爽网站| 国产伦人伦偷精品视频| 亚洲欧美日韩另类电影网站| 亚洲一卡2卡3卡4卡5卡精品中文| 新久久久久国产一级毛片| 国产男人的电影天堂91| 欧美激情 高清一区二区三区| 亚洲 欧美一区二区三区| 老司机影院毛片| 免费日韩欧美在线观看| 亚洲av国产av综合av卡| 伊人久久大香线蕉亚洲五| 精品国产一区二区久久| 国产亚洲精品第一综合不卡| 免费av中文字幕在线| 男女床上黄色一级片免费看| 丝袜美足系列| 亚洲七黄色美女视频| 侵犯人妻中文字幕一二三四区| 制服诱惑二区| 国产日韩欧美亚洲二区| 欧美少妇被猛烈插入视频| 99国产精品免费福利视频| 99久久综合免费| 亚洲国产精品一区二区三区在线| 亚洲中文字幕日韩| 一区二区av电影网| 欧美在线黄色| 天堂8中文在线网| 成人亚洲精品一区在线观看| 色老头精品视频在线观看| 久久久久国内视频| 91字幕亚洲| 91麻豆av在线| 深夜精品福利| 精品久久久精品久久久| 亚洲一区二区三区欧美精品| 18禁观看日本| 欧美日韩亚洲国产一区二区在线观看 | 日本av手机在线免费观看| 岛国毛片在线播放| 老司机亚洲免费影院| 免费高清在线观看视频在线观看| 自线自在国产av| 国产男人的电影天堂91| 一本一本久久a久久精品综合妖精| 中文字幕高清在线视频| 啪啪无遮挡十八禁网站| 久久久久久久大尺度免费视频| 99热网站在线观看| 热99久久久久精品小说推荐| 国产成人a∨麻豆精品| 黄频高清免费视频| 一级黄色大片毛片| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 国产激情久久老熟女| 深夜精品福利| 日韩免费高清中文字幕av| 国产精品久久久久久人妻精品电影 | 超碰97精品在线观看| 一区二区av电影网| 精品人妻在线不人妻| 国精品久久久久久国模美| 中文字幕另类日韩欧美亚洲嫩草| 一级黄色大片毛片| 亚洲国产av新网站| 在线永久观看黄色视频| 91精品国产国语对白视频| 69av精品久久久久久 | 亚洲精品国产精品久久久不卡| 欧美 亚洲 国产 日韩一| 欧美日韩成人在线一区二区| 国产日韩欧美亚洲二区| 中国国产av一级| 国产男女超爽视频在线观看| 亚洲五月色婷婷综合| 精品国产乱码久久久久久小说| 亚洲第一av免费看| xxxhd国产人妻xxx| 国产福利在线免费观看视频| 飞空精品影院首页| 2018国产大陆天天弄谢| 岛国毛片在线播放| 欧美日韩av久久| 日韩有码中文字幕| 成年动漫av网址| av国产精品久久久久影院| 亚洲七黄色美女视频| 婷婷成人精品国产| 精品久久蜜臀av无| www.999成人在线观看| 超碰成人久久| 国产欧美日韩一区二区精品| 欧美97在线视频| 啦啦啦中文免费视频观看日本| 狠狠狠狠99中文字幕| 丝袜在线中文字幕| 2018国产大陆天天弄谢| 国产精品九九99| 丝袜脚勾引网站| 女警被强在线播放| 午夜福利一区二区在线看| 亚洲中文日韩欧美视频| 免费看十八禁软件| 涩涩av久久男人的天堂| 亚洲情色 制服丝袜| 肉色欧美久久久久久久蜜桃| 亚洲精品中文字幕在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 日本黄色日本黄色录像| 美女国产高潮福利片在线看| 国产国语露脸激情在线看| 欧美黑人精品巨大| 亚洲av片天天在线观看| 这个男人来自地球电影免费观看| 精品亚洲成a人片在线观看| 丰满迷人的少妇在线观看| 免费观看人在逋| 不卡一级毛片| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕精品免费在线观看视频| 成年人黄色毛片网站| 亚洲第一av免费看| 窝窝影院91人妻| 亚洲国产欧美一区二区综合| 亚洲精品一区蜜桃| 国产一区二区 视频在线| 免费观看人在逋| 国产一级毛片在线| 日本猛色少妇xxxxx猛交久久| 欧美日韩中文字幕国产精品一区二区三区 | 午夜91福利影院| 大香蕉久久成人网| 久久国产精品男人的天堂亚洲| 一区二区三区四区激情视频| 电影成人av| 午夜免费成人在线视频| 久久久久久久久免费视频了| 狠狠狠狠99中文字幕| 老司机靠b影院| 欧美中文综合在线视频| 亚洲性夜色夜夜综合| 欧美日韩黄片免| 午夜激情久久久久久久| videosex国产| 侵犯人妻中文字幕一二三四区| 国产区一区二久久| 十八禁网站网址无遮挡| 老司机在亚洲福利影院| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品999| av不卡在线播放| 免费在线观看黄色视频的| 80岁老熟妇乱子伦牲交| 女人爽到高潮嗷嗷叫在线视频| 如日韩欧美国产精品一区二区三区| 亚洲国产看品久久| 丝袜美足系列| 久久性视频一级片| 久久狼人影院| 黑人操中国人逼视频| 久久九九热精品免费| 日韩精品免费视频一区二区三区| 大片免费播放器 马上看| 亚洲五月婷婷丁香| 精品一区二区三卡| 大片免费播放器 马上看| 久久九九热精品免费| 成年动漫av网址| 天天躁日日躁夜夜躁夜夜| 一区在线观看完整版| 亚洲精品久久成人aⅴ小说| 丰满少妇做爰视频| 国产熟女午夜一区二区三区| 日韩欧美一区二区三区在线观看 | 人妻久久中文字幕网| 女人久久www免费人成看片| 亚洲成人免费电影在线观看| 人妻久久中文字幕网| 国产精品自产拍在线观看55亚洲 | 国产精品国产三级国产专区5o| 自拍欧美九色日韩亚洲蝌蚪91|