• <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.

    日韩欧美一区二区三区在线观看| 欧美成人一区二区免费高清观看| 欧美色欧美亚洲另类二区| 欧美在线黄色| 日韩中文字幕欧美一区二区| 91九色精品人成在线观看| 成人无遮挡网站| 欧美性猛交╳xxx乱大交人| 日本与韩国留学比较| 韩国av一区二区三区四区| 亚洲人成伊人成综合网2020| 欧美xxxx黑人xx丫x性爽| 国产一区二区亚洲精品在线观看| 五月伊人婷婷丁香| 每晚都被弄得嗷嗷叫到高潮| 在线观看舔阴道视频| 人妻丰满熟妇av一区二区三区| 亚洲人与动物交配视频| 色综合站精品国产| 国产三级中文精品| 色5月婷婷丁香| 国产成人啪精品午夜网站| 日韩欧美精品免费久久 | 嫩草影院新地址| 精品午夜福利在线看| 无遮挡黄片免费观看| 身体一侧抽搐| 亚洲三级黄色毛片| 欧美成人免费av一区二区三区| 最后的刺客免费高清国语| 免费看美女性在线毛片视频| 丰满乱子伦码专区| 97超视频在线观看视频| 在线观看一区二区三区| 久久婷婷人人爽人人干人人爱| 久久欧美精品欧美久久欧美| 国产精品久久久久久久电影| 又黄又爽又免费观看的视频| 国产熟女xx| 丁香六月欧美| 夜夜躁狠狠躁天天躁| 中文字幕免费在线视频6| 久久久国产成人免费| 精品久久久久久久末码| 岛国在线免费视频观看| 久久久久久大精品| 香蕉av资源在线| 色综合亚洲欧美另类图片| 亚洲精品亚洲一区二区| 亚洲经典国产精华液单 | avwww免费| 欧美三级亚洲精品| 亚洲一区二区三区不卡视频| 精品人妻一区二区三区麻豆 | 久久99热6这里只有精品| 欧美在线一区亚洲| 日韩欧美免费精品| 亚洲av日韩精品久久久久久密| 热99在线观看视频| 婷婷精品国产亚洲av在线| 两性午夜刺激爽爽歪歪视频在线观看| 国内精品久久久久精免费| 国产大屁股一区二区在线视频| 美女cb高潮喷水在线观看| 男女床上黄色一级片免费看| 91av网一区二区| 国产免费一级a男人的天堂| 男女那种视频在线观看| 国产精品久久久久久久久免 | 欧美成人免费av一区二区三区| 18美女黄网站色大片免费观看| 天美传媒精品一区二区| 亚洲精品色激情综合| 亚洲av成人不卡在线观看播放网| 伦理电影大哥的女人| 麻豆成人av在线观看| 久久香蕉精品热| 国产黄片美女视频| 一进一出抽搐gif免费好疼| 午夜福利在线在线| 久久国产乱子伦精品免费另类| 亚洲美女视频黄频| 久久精品综合一区二区三区| 黄色女人牲交| 国产三级在线视频| 亚洲一区二区三区色噜噜| 午夜福利成人在线免费观看| 1024手机看黄色片| 一边摸一边抽搐一进一小说| 日本免费a在线| 高清日韩中文字幕在线| 成人精品一区二区免费| 国产精品亚洲av一区麻豆| 99国产综合亚洲精品| 黄色女人牲交| 少妇高潮的动态图| 国产伦精品一区二区三区四那| 欧美另类亚洲清纯唯美| 亚洲国产色片| 亚洲五月天丁香| 欧美激情在线99| 亚洲五月天丁香| 69av精品久久久久久| 日韩欧美免费精品| 舔av片在线| 亚洲七黄色美女视频| 51午夜福利影视在线观看| 亚洲av熟女| 嫩草影院新地址| 日韩欧美精品免费久久 | 精品一区二区三区视频在线| 国产亚洲av嫩草精品影院| 亚洲aⅴ乱码一区二区在线播放| 亚洲性夜色夜夜综合| 成人特级av手机在线观看| 深夜精品福利| 亚洲国产欧美人成| 色哟哟·www| 国产精品日韩av在线免费观看| 久久这里只有精品中国| 亚洲片人在线观看| 亚洲精品日韩av片在线观看| 亚洲乱码一区二区免费版| 国产av在哪里看| 亚洲av.av天堂| 免费大片18禁| 国产成人av教育| 久久草成人影院| 国内精品久久久久久久电影| 国产成人aa在线观看| 夜夜夜夜夜久久久久| 丰满人妻熟妇乱又伦精品不卡| 最近最新中文字幕大全电影3| 一级av片app| 久久亚洲真实| 欧美日韩乱码在线| 国产高清有码在线观看视频| 国产一区二区三区在线臀色熟女| 偷拍熟女少妇极品色| 久久天躁狠狠躁夜夜2o2o| 成年人黄色毛片网站| 久久亚洲真实| 一级毛片久久久久久久久女| 少妇裸体淫交视频免费看高清| 亚洲av免费在线观看| 91午夜精品亚洲一区二区三区 | 国产单亲对白刺激| 国产精品久久久久久久久免 | 免费av观看视频| 国产免费男女视频| 成人国产综合亚洲| 一区二区三区激情视频| 国产高潮美女av| 人人妻,人人澡人人爽秒播| 日韩中字成人| 男人舔女人下体高潮全视频| 一a级毛片在线观看| 久久婷婷人人爽人人干人人爱| 又黄又爽又免费观看的视频| 麻豆一二三区av精品| 国内精品久久久久精免费| 国产色爽女视频免费观看| 久久久久免费精品人妻一区二区| 免费无遮挡裸体视频| 精品人妻视频免费看| 国产精品av视频在线免费观看| 亚洲激情在线av| 十八禁网站免费在线| 99国产精品一区二区蜜桃av| 久久久久国内视频| 非洲黑人性xxxx精品又粗又长| 别揉我奶头 嗯啊视频| 国内久久婷婷六月综合欲色啪| 国产久久久一区二区三区| 国产69精品久久久久777片| 蜜桃亚洲精品一区二区三区| 亚洲天堂国产精品一区在线| 脱女人内裤的视频| 国产免费男女视频| 波野结衣二区三区在线| 国产精品综合久久久久久久免费| 亚洲国产高清在线一区二区三| 亚洲色图av天堂| 国产精品电影一区二区三区| 国产午夜精品论理片| 国产毛片a区久久久久| 亚洲人成网站在线播放欧美日韩| 亚洲精品一区av在线观看| 日本 av在线| 亚洲av美国av| 乱人视频在线观看| 99精品久久久久人妻精品| 久99久视频精品免费| 美女 人体艺术 gogo| 欧美另类亚洲清纯唯美| 精品人妻一区二区三区麻豆 | 精品国内亚洲2022精品成人| 国产精品野战在线观看| 91九色精品人成在线观看| 我的老师免费观看完整版| 久久精品人妻少妇| 日韩欧美三级三区| 国产亚洲欧美在线一区二区| 舔av片在线| 国产精品久久久久久人妻精品电影| av女优亚洲男人天堂| 欧美另类亚洲清纯唯美| 午夜久久久久精精品| 国产精品野战在线观看| 天天一区二区日本电影三级| 国产熟女xx| 伦理电影大哥的女人| 亚洲专区中文字幕在线| 丁香欧美五月| 色噜噜av男人的天堂激情| 亚洲av.av天堂| 久久久久久久午夜电影| 亚洲狠狠婷婷综合久久图片| 国产黄a三级三级三级人| 久久久久久大精品| 一区二区三区四区激情视频 | 美女高潮喷水抽搐中文字幕| 99国产综合亚洲精品| 日韩欧美精品免费久久 | 午夜两性在线视频| av福利片在线观看| 国产爱豆传媒在线观看| 欧美在线黄色| АⅤ资源中文在线天堂| 91在线精品国自产拍蜜月| 亚洲欧美清纯卡通| 亚洲人与动物交配视频| 欧美日韩亚洲国产一区二区在线观看| av黄色大香蕉| 看黄色毛片网站| 国产毛片a区久久久久| 1000部很黄的大片| 国产三级在线视频| 久久精品影院6| 久久久久久久精品吃奶| 国产人妻一区二区三区在| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 搡老熟女国产l中国老女人| 亚洲第一区二区三区不卡| 日本与韩国留学比较| 国产主播在线观看一区二区| 久久精品国产自在天天线| 国产欧美日韩一区二区三| 久久久久久大精品| 免费观看人在逋| 毛片女人毛片| 久久久久性生活片| 久久久久精品国产欧美久久久| 人人妻人人看人人澡| 麻豆久久精品国产亚洲av| www日本黄色视频网| 美女高潮喷水抽搐中文字幕| 深夜精品福利| 最后的刺客免费高清国语| 欧美xxxx性猛交bbbb| 亚洲欧美精品综合久久99| 精品一区二区三区视频在线观看免费| 精品人妻视频免费看| 婷婷色综合大香蕉| 精品午夜福利在线看| 一夜夜www| 十八禁人妻一区二区| 一进一出抽搐动态| 欧美丝袜亚洲另类 | 国产高清视频在线观看网站| 我要看日韩黄色一级片| 99精品久久久久人妻精品| 校园春色视频在线观看| 波多野结衣巨乳人妻| 中文资源天堂在线| 亚洲中文字幕日韩| 亚洲激情在线av| 亚洲内射少妇av| 亚洲av第一区精品v没综合| 一个人免费在线观看电影| 精品熟女少妇八av免费久了| 国产成人a区在线观看| 又粗又爽又猛毛片免费看| 超碰av人人做人人爽久久| 极品教师在线视频| 在线观看66精品国产| 色吧在线观看| 亚洲欧美日韩高清专用| 国产 一区 欧美 日韩| 国产成+人综合+亚洲专区| 亚洲乱码一区二区免费版| 成人特级av手机在线观看| 国产精华一区二区三区| 非洲黑人性xxxx精品又粗又长| 精品久久久久久成人av| av天堂中文字幕网| 女生性感内裤真人,穿戴方法视频| 精品午夜福利视频在线观看一区| 国内精品一区二区在线观看| 国产精品一及| 亚洲五月天丁香| 男女那种视频在线观看| 午夜福利在线在线| 在线观看美女被高潮喷水网站 | 久久精品国产99精品国产亚洲性色| 99精品在免费线老司机午夜| 男人和女人高潮做爰伦理| 国产高清有码在线观看视频| 日本黄色片子视频| 欧美+日韩+精品| 一区福利在线观看| 少妇裸体淫交视频免费看高清| 成熟少妇高潮喷水视频| 九色国产91popny在线| 国产精品影院久久| 午夜老司机福利剧场| 亚洲精品乱码久久久v下载方式| 久久久久国产精品人妻aⅴ院| 我的女老师完整版在线观看| 色播亚洲综合网| 99热6这里只有精品| 9191精品国产免费久久| 久久国产乱子免费精品| 一边摸一边抽搐一进一小说| 国产色婷婷99| 国产亚洲精品综合一区在线观看| 欧美+亚洲+日韩+国产| 男女之事视频高清在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美三级亚洲精品| 嫩草影院新地址| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久久久亚洲av毛片大全| 亚洲av熟女| 中文字幕高清在线视频| 欧美bdsm另类| 别揉我奶头~嗯~啊~动态视频| 757午夜福利合集在线观看| www.999成人在线观看| 女人被狂操c到高潮| 国产精品久久久久久久电影| 亚洲av日韩精品久久久久久密| 美女cb高潮喷水在线观看| 久久久精品欧美日韩精品| 欧美3d第一页| 成年女人看的毛片在线观看| 亚洲黑人精品在线| a级一级毛片免费在线观看| 亚洲美女视频黄频| 非洲黑人性xxxx精品又粗又长| 欧美精品啪啪一区二区三区| 九色成人免费人妻av| 精品日产1卡2卡| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区精品| 中文资源天堂在线| 日本五十路高清| 嫩草影院入口| 韩国av一区二区三区四区| avwww免费| www.www免费av| 亚洲不卡免费看| 成人av一区二区三区在线看| 国产精品三级大全| 三级毛片av免费| 国产老妇女一区| 91久久精品国产一区二区成人| 国产午夜精品论理片| 有码 亚洲区| 国产精品永久免费网站| 国产野战对白在线观看| 亚洲人成电影免费在线| 中文亚洲av片在线观看爽| 精品免费久久久久久久清纯| 99精品在免费线老司机午夜| 国产精品久久久久久精品电影| 久久久国产成人精品二区| 欧美+亚洲+日韩+国产| 亚洲欧美日韩高清专用| 亚洲色图av天堂| 美女xxoo啪啪120秒动态图 | 日韩欧美国产一区二区入口| 美女cb高潮喷水在线观看| 欧美不卡视频在线免费观看| 90打野战视频偷拍视频| 九色国产91popny在线| 一区福利在线观看| 一级av片app| 国产一区二区三区视频了| 亚洲精华国产精华精| 在线观看一区二区三区| 男女之事视频高清在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲经典国产精华液单 | 日本三级黄在线观看| 无遮挡黄片免费观看| 很黄的视频免费| ponron亚洲| 亚洲av免费高清在线观看| 国产精品98久久久久久宅男小说| 18禁在线播放成人免费| 看黄色毛片网站| АⅤ资源中文在线天堂| 国产综合懂色| 少妇被粗大猛烈的视频| 在线看三级毛片| 国产高清视频在线播放一区| 日韩欧美免费精品| 午夜视频国产福利| 色吧在线观看| 日韩欧美在线二视频| 午夜福利欧美成人| 夜夜夜夜夜久久久久| 麻豆久久精品国产亚洲av| 高清毛片免费观看视频网站| 激情在线观看视频在线高清| 舔av片在线| 久久久久久久久中文| 十八禁网站免费在线| 级片在线观看| 人人妻人人澡欧美一区二区| 成年人黄色毛片网站| av福利片在线观看| 国产黄色小视频在线观看| 欧美日本视频| 人人妻人人澡欧美一区二区| 欧美日韩福利视频一区二区| 黄色一级大片看看| 久久人人爽人人爽人人片va | 又黄又爽又刺激的免费视频.| 波多野结衣高清无吗| 午夜精品久久久久久毛片777| 日韩精品青青久久久久久| 99久久精品国产亚洲精品| 亚洲国产精品999在线| 成熟少妇高潮喷水视频| 亚洲,欧美精品.| 欧美潮喷喷水| 免费看美女性在线毛片视频| 一个人看的www免费观看视频| 男人的好看免费观看在线视频| 色哟哟哟哟哟哟| 午夜激情欧美在线| 精品久久久久久久久亚洲 | 欧美最新免费一区二区三区 | 欧美最新免费一区二区三区 | av天堂中文字幕网| 精品一区二区三区视频在线| av中文乱码字幕在线| 一级av片app| 黄色女人牲交| 窝窝影院91人妻| 成年女人永久免费观看视频| 欧美绝顶高潮抽搐喷水| 国模一区二区三区四区视频| 1000部很黄的大片| 国产视频内射| 黄色女人牲交| 国产av麻豆久久久久久久| 国产精品女同一区二区软件 | 如何舔出高潮| 中文亚洲av片在线观看爽| 成人性生交大片免费视频hd| 欧美日本亚洲视频在线播放| 色尼玛亚洲综合影院| 可以在线观看的亚洲视频| 波野结衣二区三区在线| 久久久国产成人免费| 国产精品影院久久| 老司机午夜福利在线观看视频| 国产精品av视频在线免费观看| 观看免费一级毛片| 国产私拍福利视频在线观看| 高清在线国产一区| 深夜精品福利| 一本综合久久免费| 国产成人影院久久av| 亚洲 国产 在线| 蜜桃亚洲精品一区二区三区| 亚洲av美国av| 欧美日韩国产亚洲二区| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 搡老熟女国产l中国老女人| 最近最新免费中文字幕在线| 国产爱豆传媒在线观看| 亚洲av电影不卡..在线观看| 长腿黑丝高跟| 欧美+日韩+精品| www.www免费av| 成人国产一区最新在线观看| 校园春色视频在线观看| 变态另类丝袜制服| 久久精品国产亚洲av涩爱 | 日本免费一区二区三区高清不卡| 国产单亲对白刺激| 深夜精品福利| 特级一级黄色大片| 国产亚洲精品久久久com| 亚洲av成人不卡在线观看播放网| 国产久久久一区二区三区| 长腿黑丝高跟| 亚洲国产精品合色在线| 国产淫片久久久久久久久 | 午夜激情欧美在线| 婷婷六月久久综合丁香| av专区在线播放| 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看| 一级黄片播放器| 99久久精品热视频| 色综合婷婷激情| 老鸭窝网址在线观看| 日韩免费av在线播放| 一卡2卡三卡四卡精品乱码亚洲| 欧美精品国产亚洲| 老司机深夜福利视频在线观看| 亚洲五月天丁香| 好男人电影高清在线观看| 中文字幕免费在线视频6| 少妇的逼好多水| 日韩精品中文字幕看吧| 亚洲人成电影免费在线| 久久精品夜夜夜夜夜久久蜜豆| 窝窝影院91人妻| 欧美激情在线99| 日韩欧美 国产精品| 亚洲第一欧美日韩一区二区三区| 一二三四社区在线视频社区8| 免费在线观看亚洲国产| 男女做爰动态图高潮gif福利片| 午夜视频国产福利| 一区二区三区免费毛片| 国产精品影院久久| 国产综合懂色| 国产伦人伦偷精品视频| 亚洲国产高清在线一区二区三| 亚洲成a人片在线一区二区| 久久中文看片网| 中文字幕久久专区| 最近视频中文字幕2019在线8| 日本成人三级电影网站| 夜夜夜夜夜久久久久| 97超视频在线观看视频| 日韩国内少妇激情av| 最近在线观看免费完整版| 国产不卡一卡二| 看片在线看免费视频| 色精品久久人妻99蜜桃| 人人妻人人澡欧美一区二区| 不卡一级毛片| 91麻豆av在线| 亚洲精品粉嫩美女一区| 亚洲精品一卡2卡三卡4卡5卡| 三级毛片av免费| 国产欧美日韩一区二区三| 好看av亚洲va欧美ⅴa在| 日韩欧美精品免费久久 | 又爽又黄a免费视频| 舔av片在线| 免费观看av网站的网址| 亚洲精品中文字幕在线视频 | 22中文网久久字幕| 久久久久久久久久久丰满| 久久影院123| 在线天堂最新版资源| 免费看日本二区| 日本猛色少妇xxxxx猛交久久| 亚洲经典国产精华液单| 国产视频内射| 国产男人的电影天堂91| 99久久九九国产精品国产免费| 国产 精品1| 九九久久精品国产亚洲av麻豆| 欧美老熟妇乱子伦牲交| 男人和女人高潮做爰伦理| 精品国产乱码久久久久久小说| 80岁老熟妇乱子伦牲交| 插阴视频在线观看视频| 日本午夜av视频| 秋霞伦理黄片| 最后的刺客免费高清国语| 深夜a级毛片| 日韩一区二区视频免费看| 欧美日韩国产mv在线观看视频 | 国产黄片视频在线免费观看| 交换朋友夫妻互换小说| 高清毛片免费看| 成人漫画全彩无遮挡| 身体一侧抽搐| 国产91av在线免费观看| 我要看日韩黄色一级片| 三级国产精品欧美在线观看| av线在线观看网站| 成人漫画全彩无遮挡| 国产免费一级a男人的天堂| 久久精品国产a三级三级三级| 免费大片18禁| 午夜福利高清视频| 国产乱人偷精品视频| 校园人妻丝袜中文字幕| 久久精品人妻少妇| 观看免费一级毛片| 丝袜美腿在线中文| 2021天堂中文幕一二区在线观| 联通29元200g的流量卡| 97热精品久久久久久| 免费观看a级毛片全部| 小蜜桃在线观看免费完整版高清| 国产亚洲av嫩草精品影院| 少妇人妻 视频| 久久久久久久久久成人|