《World Journal of Microbiology and Biotechnology》：Integrative phylogenomic and pangenome landscape of Bacillus: insights from 10,000 genomes into taxonomy, functional potential, and biotechnological applications
编辑推荐：
本研究针对芽孢杆菌属内物种界定模糊、功能潜力不明等问题，开展了大规模整合系统发育基因组学与泛基因组学分析。利用10,839个公开基因组，研究人员通过基因组距离度量、网络分析和标签传播算法，对10,276个高质量基因组进行了操作群落划分，揭示了B. cereus sensu lato等类群内的新复合体，并构建了稳健的系统发育树。功能分析鉴定了135个抗真菌基因、15个次级代谢产物簇以及大量抗性和毒力因子，同时发现与生物肥料化相关的基因广泛存在。该研究为深入理解芽孢杆菌的多样性、生态适应及其在农业生物技术中的应用潜力提供了全面的基因组学框架。
在微生物的世界里，芽孢杆菌属是一个无处不在且功能多样的“明星”群体。它们不仅是土壤、水体和空气中的常客，还与植物、动物乃至我们的食物息息相关。从引发食物中毒的“肇事者”蜡样芽孢杆菌，到臭名昭著的炭疽芽孢杆菌，再到作为生物杀虫剂被广泛应用的苏云金芽孢杆菌，这个家族的成员扮演着截然不同的角色。另一方面，以解淀粉芽孢杆菌、贝莱斯芽孢杆菌和枯草芽孢杆菌为代表的类群，则是农业生产中的“益菌”，它们能促进植物生长、防治病害，展现出巨大的生物技术潜力。然而，尽管我们与芽孢杆菌“打交道”已久，但对其内部的“家谱”（系统发育关系）和“技能库”（功能基因）的全貌仍知之甚少。传统的分类方法，如基于16S rRNA基因的鉴定，常常“力不从心”，难以清晰分辨近缘物种，导致分类混乱。随着公共数据库中基因组数据的爆炸式增长，一个关键的科学问题摆在我们面前：如何整合这海量的基因组信息，绘制出一幅能够准确反映芽孢杆菌属多样性、进化关系、功能特性及其潜在应用价值的全景图？这正是《World Journal of Microbiology and Biotechnology》上发表的这项研究所要回答的核心问题。
为了回答这些问题，研究人员进行了一项迄今为止规模最大的芽孢杆菌属整合分析。他们从NCBI RefSeq数据库中收集了10,839个公开的芽孢杆菌基因组，并经过严格的质量控制筛选出10,276个高质量基因组用于后续分析。研究主要采用了以下关键技术方法：首先，利用基因组距离度量工具Mash计算基因组间相似性，并结合网络分析和标签传播算法对基因组进行聚类，划分出不同的“操作群落”，以实现基于全基因组的物种界定。其次，利用平均核苷酸一致性工具对复杂的群落进行精细划分。接着，从每个群落中选取代表性基因组，通过识别单拷贝直系同源基因、多序列比对和构建最大似然系统发育树，来验证基于ANI的群落划分结果。然后，对包含至少10个基因组的群落，使用Panaroo进行泛基因组分析，计算基因饱和度系数和基因组流动性，以评估其泛基因组开放程度。最后，通过序列比对和生物信息学工具，在全基因组范围内系统性地鉴定了与抗真菌活性、抗生素抗性、毒力以及生物肥料化相关的基因和次级代谢产物生物合成基因簇。
基因组筛选与分类界定
研究人员从公共数据库获取了10,839个芽孢杆菌基因组，经过质量控制后保留了10,625个。他们识别了67个具有有效发表名称的模式菌株。通过结合Mash距离、网络分析和标签传播算法，将10,276个基因组初步划分为78个操作群落。对于高度关联的群落，进一步使用FastANI进行精细划分，最终得到103个基因组学上一致（ANI ≥ 95%）的群落。基于ANI的网络分析（图1）清晰展示了这些群落之间的连接关系，特别是揭示了B. cereus sensu lato类群内部及与其他类群间的新复合体。系统发育树（图3）的构建有力地支持了基于ANI的群落划分，并清晰地展示了主要的操作类群。
物种分配的关联性比较
研究比较了NCBI RefSeq、基因组分类数据库和本研究所用标签传播算法三种分类方案对物种的基因组分配情况。结果显示，尽管存在一些显著差异，但三种方案在物种水平的基因组数量分配上具有高度显著的正相关性，表明基于网络的分类方法能有效捕获与参考数据库一致的分类学关系。
泛基因组结构与基因流动性
研究对包含至少10个基因组的48个群落进行了泛基因组分析。结果显示，芽孢杆菌属不同物种的泛基因组开放程度差异很大。泛基因组的饱和度系数与基因流动性之间存在显著的负相关关系，即开放的泛基因组往往伴随着更高的基因含量变异性。例如，B. cereus、B. thuringiensis和B. mycoides等物种表现出开放的泛基因组和高基因流动性，这可能与它们作为机会致病菌的生态广适性和高水平基因转移有关。而B. licheniformis和B. paralicheniformis则显示出相对闭合的泛基因组特征。
抗真菌策略与生物合成潜力
通过生物信息学挖掘，研究在全部103个群落中鉴定出与三大抗真菌机制相关的135个基因。第一，细胞壁降解机制，涉及几丁质酶、葡聚糖酶和壳聚糖酶，其中GH18家族是核心的碳水化合物活性酶。第二，生长抑制机制，涉及由非核糖体肽合成酶、聚酮合酶及其杂合系统编码的多种抗菌肽和次级代谢产物，如表面活性素、丰原素、伊枯草菌素、杆菌霉素等。利用antiSMASH工具预测了15个与抗真菌活性相关的次级代谢产物生物合成基因簇，其中丰原素分布最广。第三，营养竞争机制，主要是铁载体bacillibactin的合成。分析发现，贝莱斯芽孢杆菌、解淀粉芽孢杆菌、枯草芽孢杆菌等植物有益菌类群拥有最广泛的抗真菌基因库，涵盖了全部三种机制。
抗生素抗性、毒力与生物肥料化基因
研究还系统鉴定了抗生素抗性决定簇、毒力因子和生物肥料化相关基因。共鉴定出135个抗生素抗性决定簇和115个毒力因子，它们主要富集在蜡样芽孢杆菌类群中。同时，与磷、氮、铁载体、钾和硫代谢相关的生物肥料化基因在芽孢杆菌属中广泛保守，凸显了该属在可持续农业中的应用潜力。
本研究通过对万余个芽孢杆菌基因组进行大规模整合分析，成功绘制了该属高分辨率的系统发育和泛基因组图谱。研究不仅厘清了其复杂的分类学关系，揭示了新的物种复合体，还阐明了泛基因组开放性与基因流动性之间的普遍规律。功能分析表明，芽孢杆菌属，特别是那些植物促生类群，是抗真菌水解酶、次级代谢产物和营养竞争性状的丰富资源库。尽管蜡样芽孢杆菌类群中富集了抗性和毒力因子，但整个属在生物肥料化相关基因上表现出广泛的保守性。这些发现不仅深化了我们对芽孢杆菌生物学多样性、生态适应和进化机制的理解，更重要的是，为未来定向筛选安全、高效的生防菌株和生物肥料，开发新一代农业生物制剂提供了宝贵的基因组学资源和理论指导。该研究建立的分析框架和产生的数据集，将成为后续探索细菌物种形成、功能进化及其与环境互作的坚实基础。
