《Journal of Fungi》：Comparative Phylogenomics Reveals Conserved Protoilludene Synthase Neighbourhoods and Unevenly Distributed Biosynthetic Gene Cluster Families in Armillaria
编辑推荐：
真菌生物合成基因簇（biosynthetic gene clusters, BGCs）为天然产物发现提供了基因组切入点，但仅凭预测区域既不能确立通路活性，也不能鉴定代谢产物。蜜环菌属（Armillaria）是一个有用的模型，因为该属化学丰富、与森林相关的物种涵
真菌生物合成基因簇（biosynthetic gene clusters, BGCs）为天然产物发现提供了基因组切入点，但仅凭预测区域既不能确立通路活性，也不能鉴定代谢产物。蜜环菌属（Armillaria）是一个有用的模型，因为该属化学丰富、与森林相关的物种涵盖病原性和腐生性生活方式。研究人员采用角色感知设计分析了九个组装基因组：七个高完整性独立分类群（使用基准通用单拷贝同源基因（benchmarking universal single-copy orthologs, BUSCO）评估）、一个A. mellea BUSCO重复敏感性样本，以及一个作为A. ostoyae种内菌株对照的A. solidipes样本。该分析解析出一个富含萜类的基因库，其中包含一个共享但对边缘稳健的基因簇家族（gene cluster family, GCF）组分，以及分布不均的候选家族。靶向邻域分析恢复出所有组装中的protoilludene合酶（PRO1）直系同源基因（氨基酸一致性92.2%–100.0%）；核心邻域在七个主要分类群中均保留，而侧翼的修饰相关基因则有所变化。值得注意的是，A. gallica的PRO1直系同源基因位于最近的antiSMASH区域之外23,614 bp处。因此，全区域预测可能掩盖担子菌中以通路为中心的保守性。所得框架优先考虑可重复的候选区域，同时将化学同一性、表达和生态功能留给直接验证。
真菌特化代谢物介导竞争、宿主互作和环境胁迫响应，许多合成所需的酶、转运蛋白和调控因子由生物合成基因簇（biosynthetic gene cluster, BGC）编码。基因组挖掘可揭示常规培养实验不明显的生物合成潜力，但预测区段并不代表表达或功能。蜜环菌属（Armillaria）包含病原菌和腐生菌，化学丰富，产生结构多样的倍半萜芳基酯（包括melleolide类），是研究担子菌BGC进化的理想模型。然而，全属范围内以通路为中心的邻域保守性尚未系统解决。原始BGC数量比较有限，需要结合系统基因组框架、明确样本角色和敏感性分析。
研究人员采用角色感知设计分析了9个来自NCBI GenBank的蜜环菌基因组组装，包括7个独立分类群作为核心组，A. mellea作为BUSCO重复敏感性样本，A. solidipes作为A. ostoyae的种内菌株对照。研究整合基因组完整性、单拷贝标记系统发育、antiSMASH预测、BiG-SCAPE GCF推断、边界敏感性、成对GCF相似性、KCB证据和靶向邻域分析。结果显示7个主要分类群中富含萜类的BGC库包含共享且边缘稳健的GCF组分，以及分布不均的候选家族；PRO1直系同源基因在所有组装中保守，核心邻域保留而侧翼修饰基因变化；A. gallica的PRO1位于最近antiSMASH区域外23,614 bp。该研究发表在《Journal of Fungi》，为优先筛选蜜环菌特化代谢位点建立了可重复且尺度感知的基础。
研究用到的主要关键方法包括：从NCBI GenBank获取9个组装；使用BUSCO v5.8.0评估蛋白质模式完整性；用BRAKER3进行基因注释并用AGAT选择最长异构体；用antiSMASH v8.0.4预测BGC区域；用BiG-SCAPE v2.0.3在0.20、0.30、0.40截断值下进行GCF聚类；用MAFFT和IQ-TREE基于661个单拷贝BUSCO标记构建最大似然系统发育树；用DIAMOND搜索PRO1直系同源基因，并用clinker分析邻域同源性。样本队列来源为NCBI GenBank（7个独立分类群加2个对照组装）。
研究结果如下：
3.1 高注释完整性支持比较BGC谱分析：蛋白质模式BUSCO完整性为96.7%–98.7%，基因模型数量为11,382–15,940，表明后续比较具有可靠基础。
3.2 全基因组谱分析揭示萜类主导的BGC景观：483个BGC区域覆盖12种产物组合；萜类单一产物占比49.9%（241/483），含杂交组合时276个区域带萜类标签；NRPS-like次之（98/483，20.3%）。
3.3 共享生物合成核心嵌入不均匀GCF占位：在0.30截断值下，424个区域形成75个GCF，59个单例；74个主组GCF中18个为核心、22个广泛共享、27个限制、7个私有。
3.4 截断值、边界和菌株对照保持主要占位模式：0.20和0.40截断值分别产生78和70个GCF，占位稳定；移除12个contig-edge区域仅改变3个GCF；A. mellea仅新增1个mellea-only GCF；A. ostoyae与A. solidipes共享34/42和43个GCF，并共享730个单拷贝BUSCO标记。
3.5 GCF相似性仅部分遵循系统基因组邻近性：基于661个标记的系统发育树所有内部节点支持率均为100%超快自举和≥99.9% SH-aLRT；成对GCF Jaccard相似性为0.407–0.769，最高为A. altimontana与A. cepistipes，最低为A. cepistipes与A. fumosa。
3.6 大多数GCF缺乏与已知参考簇的匹配：KCB命中53/483区域，覆盖12/75 GCF；63个GCF无KCB命中；18个核心GCF中15个无KCB支持；所有GCF至少含一个proto-core相关CDS。SubClusterBlast无命中，因安装的参考亚簇均为细菌来源，反映数据库覆盖限制。
3.7 正交证据优先考虑边缘稳定的萜类候选：18个核心GCF中16个为边缘稳健；3个核心GCF有KCB支持，15个无；GCF46（A. ostoyae和A. solidipes）和GCF49（A. gallica、A. ostoyae、A. solidipes）被选为序列级候选，但仅是计算靶标，不确证代谢产物身份或功能。
3.8 Protoilludene合酶邻域的保守性与结构变异：所有9个组装均检出PRO1直系同源基因，氨基酸一致性92.2–100%；7个主要分类群中PRO1、脱氢酶/还原酶样和GMC氧化还原酶样直系群保守，甲基转移酶样直系群出现于6个分类群；P450样基因在7个邻域中均有代表，侧翼补体变化。A. fumosa与A. mexicana虽为姐妹类群，但PRO1邻域结构不同。antiSMASH和BiG-SCAPE将8个区域包含的直系同源基因分到不同GCF，而A. gallica的PRO1距最近antiSMASH区域23,614 bp，说明区域边界调用可能掩盖通路中心保守性。
讨论部分指出，稳定的共享组分嵌入更大的不均匀分布库中，萜类区域占优与蜜环菌倍半萜芳基酯化学一致。靶向邻域比较表明，同一通路中心邻域可能被分裂到多个GCF，且边界调用可能掩盖保守性。角色感知的对照设计降低了重复和菌株混淆造成的过度解释。KCB提供独立证据轴，但参考数据库在担子菌中覆盖不均（MIBiG v3.0中担子菌条目仅23条），SubClusterBlast无真菌参考。局限性包括预测对组装连续性和注释敏感、未运行ClusterCompare、相似性分析不等同于功能。最后提出验证流程：先基于独立属性排序GCF，再通过表达分析、定向突变或异源重组和LC-MS/MS确证。
结论部分翻译如下：对九个蜜环菌组装的角色感知比较鉴定出483个预测BGC区域，并将424个聚集区域组织成75个GCF。采样基因库富含萜类，同时包含边缘稳健的共享组分和较大的不均匀分布组分。靶向分析揭示了一个跨主要分类群保守的PRO1中心邻域，尽管侧翼基因可变、多个全区域GCF归属以及A. gallica中一个antiSMASH边界遗漏。参考簇相似性仅覆盖12个GCF，而边缘稳定的无命中萜类GCF为实验表征提供了聚焦候选。这些结果为优先筛选蜜环菌特化代谢位点建立了系统发育受控且尺度感知的基础，同时将通路活性、化学产物和生态角色留给直接验证。
