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

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

關鍵詞：暹羅芽孢桿菌；鏈格孢菌；揮發(fā)性有機化合物；轉錄組分析

中圖分類號：S435.72? ? ? ? ? ? ? ? ? ? 文獻標識碼：A? ? ? ? ? ? ? ? ? ? 文章編號：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.）引起的一種重要的葉部真菌性病害，具有潛育期短、流行速度快等特點[1]，是我國煙草生產上威脅最大的病害之一，直接影響煙葉的產量和品質。目前，煙草生產上對赤星病的防治主要以化學防治為主，但是長期大量使用化學藥劑會導致病原菌產生抗藥性，并造成煙葉農藥殘留問題，影響煙葉品質[2-3]。生物防治能有效避免或減少對生物和環(huán)境造成的不利影響，在煙草病蟲害的防控中具有極大的應用前景。

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

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

1? 材料與方法

1.1? 材料

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

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

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

1.2? 方法

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

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

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

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

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

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

2? 結? 果

2.1? 2-甲基丁酸和3-甲基丁酸混合物對鏈格孢菌生長的影響

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

2.2? 測序結果的評估

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

2.3? 基因表達差異

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

2.4? 富集分析

2.4.1? 差異表達基因GO功能富集分析? 對下調的DEGs進行GO富集分析發(fā)現(xiàn)（圖5），細胞組分中細胞結構體、膜組成部分和膜固有成分富集基因數(shù)較多，而生物過程中脂質代謝過程富集最顯著。對上調的DEGs進行GO富集分析發(fā)現(xiàn)細胞氨基酸生物合成過程、核糖核蛋白復合體組裝和有機酸分解代謝過程3個功能最顯著富集。

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

2.5? 與鏈格孢菌生長相關代謝通路的差異表達基因分析

結合COG功能注釋、GO富集以及KEGG富集，從所有的差異表達基因中篩選出與鏈格孢菌生長相關的差異表達基因（p<0.01），下調的基因主要參與甘油磷脂代謝、MAPK信號通路及聚酮類物質的合成，上調的基因主要參與谷胱甘肽代謝、淀粉和蔗糖代謝等過程（表3）。

2.6? qRT-PCR驗證

為了檢驗轉錄組測序結果的準確性，從篩選出的差異表達基因中隨機挑選8個差異顯著的基因進行qRT-PCR驗證，其中包括5個下調表達基因和3個上調表達基因。結果表明，所選基因的表達水平與轉錄組結果趨勢一致（圖8），說明轉錄組測序的結果是可靠的。

3? 討? 論

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

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

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

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

4? 結? 論

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

參考文獻

[1]張中義，冷懷瓊，張志銘，等. 植物病原真菌學[M]. 成都：四川科學技術出版社，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ā)性有機物活性研究及菌劑研制[D]. 武漢：華中農業(yè)大學，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-甲基丁酸對甘薯長喙殼菌的抑制作用研究[J]. 微生物學報，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]王笑連. 楊樹炭疽病菌CgMk1 MAPK及上游受體蛋白的功能分析[D]. 北京：北京林業(yè)大學，2021.

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

[15]張楠，柳志強，吳曼莉，等. 膠孢炭疽菌CgSho1基因的克隆與功能分析[J]. 植物病理學報，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.