《Journal of Fungi》：Comparative Genomics and Gene Cluster Family Analysis Reveal Conserved and Isolate-Specific Biosynthetic Potential in Medicinal Plant-Associated Endophytic Diaporthe
编辑推荐：
与药用植物相关的内生真菌是特化代谢产物（specialized metabolites）的有前景储库；然而，其生物合成潜力背后的基因组多样性仍未完全表征。在此，研究人员对四株与药用植物相关的内生Diaporthe分离株——D. inconspicua LGMF
与药用植物相关的内生真菌是特化代谢产物（specialized metabolites）的有前景储库；然而，其生物合成潜力背后的基因组多样性仍未完全表征。在此，研究人员对四株与药用植物相关的内生Diaporthe分离株——D. inconspicua LGMF1612、D. cf. heveae LGMF1633、D. vochysiae LGMF1583和Diaporthe sp. HANT25——开展了比较基因组学和基因组挖掘分析。基因组大小为52.45–60.11 Mb，预测蛋白数为16,033–17,434。比较同源分析鉴定出16,196个蛋白簇，其中11,230个（69.3%）在四株分离株中保守，表明存在大量共享蛋白组库以及较小的可变组分。功能注释显示膜转运蛋白、真菌转录调控因子、细胞色素P450s（cytochrome P450s）、氧化还原和解毒酶以及与特化代谢相关的蛋白结构域广泛存在。利用antiSMASH进行基因组挖掘，鉴定出481个生物合成基因簇（biosynthetic gene clusters, BGCs），每株分离株含112–129个BGCs。I型聚酮合酶（type I polyketide synthase, T1PKS）簇占主导，同时存在大量非核糖体肽合成酶（nonribosomal peptide synthetase, NRPS）、NRPS-like、萜烯及杂合/多类群基因簇库。BiG-SLiCE将预测BGCs归为348个基因簇家族（gene cluster families, GCFs），其中仅15个被归类为核心GCFs，其成员在所有四株分离株中均有代表，揭示尽管蛋白组背景大体保守，但生物合成库存在广泛变异。与MIBiG（Minimum Information about a Biosynthetic Gene cluster）相关预测整合后，鉴定出与已表征的squalestatin S1、choline、scytalone/T3HN、cercosporin和nectriapyrone生物合成途径相似的核心GCFs。其他已表征途径相关特征显示部分或分离株限制分布，包括D. inconspicua中α-acorenol和phyllostictine相关BGCs，D. vochysiae中chaetolivacine、chrysogine、patulin和AKML相关BGCs，以及Diaporthe sp. HANT25中enniatin、azaphilone、curvupallide-B、aspyridone和wortmanamide相关BGCs。重要的是，相当一部分GCFs与已表征MIBiG簇缺乏显著对应，表明存在大量化学未表征的生物合成潜力。总体而言，这些发现揭示药用植物相关Diaporthe具有保守基因组基础以及特化代谢产物生物合成库的广泛多样化，并鉴定出保守及分离株限制的BGC家族作为未来代谢组学和功能表征的候选。
研究背景：药用植物相关内生真菌是结构多样的特化代谢产物（specialized metabolites）的重要储库，Diaporthe属兼具内生、腐生和致病等生态类型，且以产生丰富次级代谢产物著称。丝状真菌的特化代谢产物多由生物合成基因簇（biosynthetic gene clusters, BGCs）编码，常规培养条件下许多途径沉默或弱表达，因而基因组挖掘可揭示常规代谢谱难以显示的隐晦生物合成潜力。已有Diaporthe基因组研究提示其次级代谢基因变异较大，但药用植物相关内生Diaporthe分离株之间BGC及基因簇家族（gene cluster families, GCFs）水平的保守与分离株特异性仍不清楚。为此，研究人员比较四株药用植物相关内生Diaporthe分离株的基因组、蛋白组、功能注释和生物合成潜力，旨在区分共享与分离株限制的BGC家族，为后续代谢组学和功能研究提供候选。
方法概括：四株内生Diaporthe分离株分别来自Maytenus ilicifolia（D. inconspicua LGMF1612）、Vochysia divergens（D. cf. heveae LGMF1633、D. vochysiae LGMF1583）和Hydnocarpus anthelminthicus（Diaporthe sp. HANT25）；基因组和蛋白序列取自JGI MycoCosm与NCBI GenBank。采用ITS、TUB2、HIS3、TEF1串联序列在MEGA11构建最大似然树；OrthoVenn3/OrthoMCL进行同源蛋白聚类；eggNOG-mapper、Pfam和Gene Ontology（GO）注释；antiSMASH 8.0.4真菌模式预测BGCs，KnownClusterBlast/MIBiG 4.0以≥45%参考簇基因覆盖率为操作阈值；BiG-SLiCE 1.0以距离0.4聚类GCFs；clinker比较选定BGCs。
研究结果：
3.1 四株Diaporthe分离株的基因组特征：通过基因组比较和多位点系统发育分析，四株基因组大小为52.45–60.11 Mb，GC含量51.68%–52.73%，预测蛋白16,033–17,434个；D. inconspicua LGMF1612、D. cf. heveae LGMF1633和Diaporthe sp. HANT25位于同一主要谱系，其中D. cf. heveae LGMF1633与Diaporthe sp. HANT25形成100% bootstrap支持分支，D. vochysiae LGMF1583与D. ampelina聚类。比较同源分析获得16,196个同源簇，11,230个（69.3%）为四株共享核心簇，2,574个和2,240个分别由三株和两株共享，仅152个（0.9%）为分离株特异；D. vochysiae LGMF1583分离株特异簇最多（85个）。结论是四株具有大体共享的蛋白组背景，同时存在较小的部分共享、分离株特异和未聚类组分。
3.2 与生态适应和特化代谢相关的功能特征：通过Pfam和GO比较分析，四株蛋白结构域谱总体相似，膜转运相关结构域丰富，MFS_1（PF07690）在452–456个蛋白中检出，p450（PF00067）在366–390个蛋白中检出，Zn_clus（PF00172）、Fungal_trans（PF04082）等转录调控相关结构域以及adh_short、FAD_binding、GMC_oxred_N、Aldo_ket_red、过氧化物酶、硫氧还蛋白和谷胱甘肽S-转移酶等氧化还原/解毒相关结构域广泛存在；PKS、NRPS/载体蛋白及甲基转移酶、乙酰转移酶等tailoring相关结构域亦分布于四株。结论是转运、转录调控、氧化还原、解毒、细胞色素P450和特化代谢相关功能在四株中广泛存在。
3.3 生物合成基因簇的分布与多样性：通过antiSMASH基因组挖掘，四株共预测481个BGCs，每株112–129个，Diaporthe sp. HANT25最多（129个），D. vochysiae LGMF1583最少（112个）。T1PKS相关区域在四株中均最丰富（53–58个），NRPS为17–25个，NRPS-like为18–22个，萜烯相关为24–40个且D. cf. heveae LGMF1633最多，indole相关为6–9个；杂合/多类群BGCs分别为24、28、21和23个。结论是四株BGC类别组成总体相似，但萜烯和NRPS等特定类别的相对比例存在差异。
3.4 共享生物合成途径及其跨Diaporthe分离株的GCF组织：通过KnownClusterBlast/MIBiG与BiG-SLiCE整合分析，348个GCFs中15个为核心GCFs；GCF84、GCF109、GCF197、GCF214和GCF241分别在四株中与squalestatin S1、choline、scytalone/T3HN、cercosporin和nectriapyrone相关MIBiG途径相似。另有一些已表征途径在四株中均检出但分布于多个GCFs，如(?)-mellein相关BGCs分布于GCF208、GCF237和GCF319，alternariol相关BGCs分布于GCF200、GCF203、GCF295和GCF193；三株共享的已表征相关GCFs包括GCF22（clavaric acid）、GCF278（alternapyrone）和GCF303（duclauxin，且为部分表征相关）；两株共享的19个GCFs中包括ochrindole A、pimara-8(14),15-diene、huperzine A、betaenones A–C、stipitatic acid和zearalenone等相关家族。结论是四株存在共享的核心已表征途径相关潜力，但BGC在GCF层面的组织存在广泛变异。
3.5 分离株特异性已表征途径相关BGCs：通过分离株限制性GCF与MIBiG关联分析，共有75个已表征途径相关GCFs仅存在于单一分离株，其中D. inconspicua LGMF1612为17个、D. cf. heveae LGMF1633为19个、D. vochysiae LGMF1583为22个、Diaporthe sp. HANT25为17个，另有1个部分表征相关GCF仅见于Diaporthe sp. HANT25。D. inconspicua LGMF1612含α-acorenol相关GCF87、armillyl orsellinate/8α-hydroxy-6-protoilludene相关GCF106、phyllostictine A/B相关GCF341和agnestin A/B相关GCF289；D. cf. heveae LGMF1633含copalyl diphosphate相关GCF74；D. vochysiae LGMF1583含chaetolivacine相关GCF292、chrysogine相关GCF163、patulin相关GCF259和AKML相关GCF227；Diaporthe sp. HANT25含enniatin相关GCF85、azaphilone相关GCF328、curvupallide-B相关GCF310、aspyridone相关GCF261和wortmanamide相关GCF217。结论是四株预测的已表征途径相关生物合成能力存在明显分离株差异。
3.6 未表征基因簇家族的分布与分离株特异性：通过GCF分类，348个GCFs中91个（26.1%）为已表征途径相关GCFs，15个（4.3%）为部分表征相关GCFs，242个（69.5%）为化学未表征GCFs。未表征GCFs包括7个四株核心家族、11个三株共享家族、37个两株共享家族和187个分离株限制家族；7个核心未表征GCFs均无保留的MIBiG关联，但含可识别的PKS、NRPS/NRPS-like或萜烯核心酶；15个部分表征相关GCFs包括4个核心、1个三株共享、9个两株共享和1个分离株限制家族。结论是四株编码大量未表征生物合成潜力，其中分离株限制未表征GCFs数量最多。
讨论部分总结：讨论指出Diaporthe是特化代谢生物合成特别丰富的真菌类群，四株112–129个BGCs的数量远超单一培养条件代谢谱可推断的潜力；尽管同源蛋白簇高度保守，GCF成员关系变异更大，说明特化代谢途径是相对可变的基因组组分。核心GCFs提示四株共享生物合成基础，但GCF成员关系不代表完整序列、基因顺序或严格直系同源。分离株特异性BGCs提示预测生物合成库的化学分化；大量未表征GCFs，尤其分离株限制家族，代表未开发生物合成潜力。该研究结果与药物和农业天然产物发现相关，但BGC预测只表明生物合成潜力，不等同于途径表达、产物生成或生物活性；基因组组装和基因预测不完整、未纳入腐生和致病Diaporthe比较、分离株限制BGCs不能解释为内生特异等均为限制。后续需结合OSMAC（One Strain Many Compounds）培养、共培养、表观遗传调控、异源表达、LC-HRMS/MS代谢组学、分子网络、转录组和遗传破坏或异源表达进行验证。
研究结论翻译：对四株药用植物相关内生Diaporthe分离株的比较基因组分析揭示，其具有保守蛋白组背景并伴随特化代谢生物合成潜力的显著多样化。尽管69.3%的同源蛋白簇在四株中共享，481个预测BGCs被归为348个GCFs，其中仅15个在四株基因组中均有代表。T1PKS、NRPS、NRPS-like、萜烯和杂合生物合成系统占主导，加之相关tailoring酶的多样性，突出其广泛的预测生物合成能力。核心GCFs和分离株限制生物合成特征进一步表明共享及分离株依赖的生物合成潜力，而相当比例的化学未表征GCFs提示基因组引导天然产物发现的机会。重要的是，所鉴定的BGCs和GCFs代表计算预测的候选生物合成位点，而非实验验证的途径，其表达、代谢产物和生物活性仍有待确定。未来研究应基于共享或分离株限制分布、与已表征MIBiG途径的差异、独特tailoring酶组成及不寻常核心生物合成架构优先选择候选BGCs；靶向基因表达分析和基于LC-MS/MS的代谢组学，随后对选定簇进行遗传破坏或异源表达，对于确定其生物合成功能和产物至关重要。
