《BMC Genomics》：The interplay of specialized metabolites, antibiotic resistance, and virulence in Enterococcus faecium: in silico analysis of bacterial genomes
编辑推荐：
背景：粪肠球菌（Enterococcus faecium）是全球范围内重要的革兰阳性医院感染病原体，与高发病率和高死亡率密切相关。识别其抗生素耐药性（Antimicrobial Resistance, AMR）基因、毒力因子及编码特化代谢物的生物合成基因簇（B
背景：粪肠球菌（Enterococcus faecium）是全球范围内重要的革兰阳性医院感染病原体，与高发病率和高死亡率密切相关。识别其抗生素耐药性（Antimicrobial Resistance, AMR）基因、毒力因子及编码特化代谢物的生物合成基因簇（Biosynthetic Gene Clusters, BGCs），是开发新型治疗干预措施的重要前期工作。本研究旨在利用antiSMASH、GECCO和Clinker等多种计算工具，探索从基因组分类数据库（Genome Taxonomy Database, GTDB）获取的81株E. faecium基因组的BGCs、AMR与毒力属性之间的关系，并评估这些特征的出现频率、共现情况及关联模式。结果：在检测到的334个BGCs中，近50%被鉴定为RiPP样基因簇，多数预测编码细菌素相关产物；部分环状内酯自诱导剂基因簇（Cyclic Lactone Autoinducer BGCs, CAL-BGCs）表现出基因组重排或独特序列，提示群体感应潜力或信号系统的多样性。此外，还鉴定出若干初级代谢基因簇（Metabolic Gene Clusters, MGCs），包括精氨酸脱亚胺酶系统、PFOR-II途径和没食子酸代谢通路，提示其在能量适应和肠道定植中的作用。BGCs和MGCs均呈现较高的保守性。毒力基因acm、scm、fss3、sgrA和ecbA在不同菌株中存在差异，表明宿主细胞结合能力和生物膜形成能力存在差异。AMR谱显示大多数基因组具有极高的耐药性，共检测到43种AMR基因，包括万古霉素耐药簇（vanA型）、aac(6')-Ii和msrC等；约70%的基因组同时对5类及以上抗生素耐药。流行病学整合分析显示，序列型（Sequence Type, ST）在地理来源和临床分离源之间存在明显结构分化，AMR负担呈谱系依赖性变化，而毒力和BGC谱相对保守。结论：本研究揭示了E. faecium的代谢与基因组可塑性，尤其是其动态变化的AMR谱反映了高致病性和环境适应性。结果强调同步探索特化代谢物、耐药性与毒力特征的重要性，有助于深入理解该病原体的生存机制，并为发现潜在的新治疗靶点提供依据。然而，这些计算结果的进一步实验验证将加深对其生物学特性的认识，并助力应对AMR细菌的挑战。
研究背景与意义
粪肠球菌（Enterococcus faecium）是ESKAPE病原体之一，因其在医院感染中的高流行率、多重耐药特性及高致病性，成为全球公共卫生关注的重点。随着广谱抗生素的过度使用和侵入性医疗操作的增加，耐万古霉素粪肠球菌（Vancomycin-Resistant Enterococcus faecium, VREfm）已成为医疗环境中最具威胁的耐药菌之一。目前，针对E. faecium的特化代谢物（Specialized Metabolites, SMs）与其抗生素耐药性（AMR）、毒力因子的协同作用机制尚缺乏系统性研究。因此，整合基因组学手段全面解析这些特征的共存与互作关系，对揭示其致病机制和制定防控策略具有重要意义。该研究由研究人员完成并发表于《BMC Genomics》。
主要技术方法
研究人员从GTDB数据库中筛选了81株高质量E. faecium完整基因组（CheckM完整性≥99.6%，无污染且无模糊碱基）。采用antiSMASH 7.0和GECCO v0.9.10两种工具预测BGCs，并利用Clinker v0.0.31进行基因簇比对与可视化，同时借助PRISM 4预测SM化学结构。初级代谢基因簇（MGCs）通过gutSMASH v1.0鉴定。AMR基因通过CARD数据库与ABRicate双重筛查，毒力基因则通过VFDB数据库注释。分子分型采用PubMLST平台确定ST型别。基因共现网络分析基于phi系数和相关统计方法构建，并通过R语言包进行数据可视化与统计分析。
研究结果
BGCs检测结果：antiSMASH共检出334个BGCs，包括RiPP样（166个）、RRE包含型（2个）、III型聚酮合酶（81个）和环状内酯自诱导剂（85个）四类；GECCO检出83个BGCs，主要为RiPP类和糖类BGCs，部分未知BGCs经结构域分析推测可能为新型细菌素样基因簇。CAL-BGCs在基因组结构上高度相似，仅少数存在重排或插入缺失；RiPP样BGCs的核心结构高度保守，但侧翼区存在变异。
初级代谢基因簇：共鉴定出四个主要通路，分别为精氨酸转化为碳酸氢盐途径（精氨酸脱亚胺酶系统）、PFOR-II途径、丙酮酸转化为乙酸-甲酸途径，以及没食子酸代谢通路。
毒力因子基因：检测到acm、scm、fss3、sgrA和ecbA五种毒力基因，其中acm在所有基因组中普遍存在。主要克隆谱系为ST736、ST17、ST18和ST80，各占10–20%。
AMR基因：共检出43种AMR基因，每基因组携带2–20个不等。96.3%的菌株携带超过3类AMR基因，约70%携带至少5类抗生素耐药基因，呈现广泛耐药表型。81.5%的菌株为vanA阳性，氨基糖苷类、大环内酯类、四环素类等耐药基因型几乎普遍存在。
相关性及网络分析：基因共现网络显示，vanA、vanB和vanM耐药操纵子内部基因完全共现，并与多种其他AMR基因紧密关联；CAL-BGCs与efmA呈完全互斥分布；lsaE与vanA操纵子呈显著负相关。跨类别关联以AMR主导为特征，AMR与毒力关联有限，AMR与特化代谢物呈负关联。
基因组特征与流行病学整合：AMR负担在直肠/粪便分离株和来自美国的分离株中最高。ST型分布与地理区域、分离源及疾病类型显著相关，ST736多见于血流感染，ST17和ST18多见于直肠标本。AMR负担在ST17和ST18中最高，ST736和ST80中等，ST78较低。毒力基因和BGC含量在各ST间无明显差异。个别耐药基因（如tet(L)、dfrF）和毒力基因（如ecbA、scm）呈现显著的ST特异性分布。
讨论与结论
研究人员指出，E. faecium的BGCs以RiPP样和CAL-BGCs为主，部分未知BGCs可能编码新型细菌素或信号调控模块。初级代谢通路的保守性提示其作为稳定抗菌靶点的潜力。毒力基因的差异分布反映了菌株在宿主黏附与定植能力的异质性。AMR基因的广泛分布及vanA型VRE的高比例凸显了临床治疗的严峻挑战。基因共现网络揭示了非随机的遗传组织模式，耐药操纵子的完整共现与某些代谢基因簇的互斥分布反映了生态位分化和选择压力的影响。流行病学分析显示，高风险克隆（如ST736、ST17、ST18）在全球范围内传播，并与特定耐药及毒力特征紧密相关。研究人员强调，整合代谢、耐药与毒力数据的策略可为E. faecium的功能研究和靶向治疗开发奠定坚实基础。
结论翻译：本研究通过对E. faecium的综合基因组分析，揭示了其生物合成基因簇（BGCs）、抗生素耐药性（AMR）基因、毒力因子及初级代谢通路的多样性。研究发现，E. faecium拥有一组保守的核心BGCs（如RiPP样和环状内酯自诱导剂BGCs），部分基因组携带新颖或分化的基因簇，可能在群体感应调控或其他信号通路中发挥谱系特异作用。相对保守的BGCs和MGCs可作为新的治疗靶点。vanA型VRE的主导地位及广泛的多重耐药性凸显了这些分离株在临床和社区环境中的重大公共卫生威胁。高风险毒力谱提示其具有高致病性、强侵袭性及环境适应能力。基因共现网络分析揭示了非随机的遗传关联模式，正、负相关性反映了塑造基因组组成与适应的环境驱动力。普遍存在的黏附因子、高比例的侵袭酶及广泛的抗生素耐药性共同构成高毒力病原体特征。这些发现强调，整合代谢、耐药与毒力数据是开发针对性防治策略的关键基础。
