关于克雷伯氏菌氧酸杆菌(Klebsiella oxytoca)种内分化以及抗生素抗性/毒力因子基因潜在进化机制的基因组学研究

时间:2026年5月30日
来源:Microbial Pathogenesis

编辑推荐:

王玉童|韩海伟|青云明|魏冠云|张贤中南大学湘雅公共卫生学院职业与环境健康系,中国长沙410013摘要产氧克雷伯菌(Klebsiella oxytoca)因其临床重要性和生态适应性而广为人知,然而对其全基因组特征的系统性研究以及与临床和环境适应相关的基因分布的探究仍然有限。本研究

广告
   X   

王玉童|韩海伟|青云明|魏冠云|张贤
中南大学湘雅公共卫生学院职业与环境健康系,中国长沙410013

摘要

产氧克雷伯菌(Klebsiella oxytoca)因其临床重要性和生态适应性而广为人知,然而对其全基因组特征的系统性研究以及与临床和环境适应相关的基因分布的探究仍然有限。本研究对598株产氧克雷伯菌进行了比较基因组分析,以描绘该物种的基因组结构,并分析其中的抗生素抗性基因(ARGs)、毒力因子基因(VFGs)和移动遗传元件(MGEs)的分布情况。分析结果显示,所有菌株共拥有1239个保守的核心基因,辅助基因组平均包含4373.5个基因,菌株特异性基因平均为11个,表明该物种具有显著的基因多样性。从功能上看,核心基因组主要涉及细胞基本过程,而辅助基因组和菌株特异性基因则与适应性代谢相关。不同菌株间的ARGs和VFGs分布存在差异,平均每个菌株分别含有18个和62个相关基因。值得注意的是,MGEs与ARGs的丰度之间存在显著正相关(Spearman's R = 0.65,P < 2.2e–16),而MGEs与VFGs之间的相关性较弱(Spearman's R = 0.29,P = 2.4e–13)。此外,这些因素的共定位进一步表明MGEs可能在产氧克雷伯菌中介导抗性和毒力基因的传播。总体而言,本研究为理解该物种的基因分化及其抗性和毒力基因的进化动态提供了基础基因组框架。

引言

产氧克雷伯菌作为克雷伯菌属中具有临床意义的成员,是一种机会性病原体,与多种感染相关,并且在临床环境中表现出日益增加的抗菌素耐药性[1]。它的生态适应性不仅限于医院[2]、[3],还在人类肠道[4]、废水[5]、土壤[6]和植物根际[7]中被发现,这凸显了其在临床和环境储存库中持续存在和传播的潜力。其适应性的关键可能在于其基因组结构,该结构包含一个稳定的核心基因组和一个由移动遗传元件(MGEs)构成的高度动态的辅助基因组,这些元件通过水平基因转移(HGT)促进外源DNA的整合[8]、[9]。尽管最近的研究已经开始从更广泛的层面探索该物种的基因组特征[10],但对产氧克雷伯菌菌株全基因组及其功能基因分布的系统性研究仍然不足。建立物种水平的基因型基线对于阐明抗生素耐药性和毒力等表型特征的遗传基础至关重要[11]、[12]、[13]。
产氧克雷伯菌因同时携带多重耐药性和毒力决定因子而日益受到关注,这些因素会复杂化感染管理并带来公共卫生挑战。值得注意的是,这些基因型特征并非随机分布,而是常常共位于MGEs上,尤其是质粒上,质粒是这些因子水平转移的载体[8]、[14]、[15]。为了全面了解这一潜在威胁,有必要超越个别病例报告,通过系统地识别和量化多种产氧克雷伯菌菌株中的抗生素抗性基因(ARGs)、毒力因子基因(VFGs)和MGEs来进行全面的物种范围调查。此外,分析基因组内容以确定共享基因和变异基因的分布情况,有助于揭示物种内的基因分化和适应潜力[16]、[17]、[18]。
高风险产氧克雷伯菌克隆的进化与通过MGEs进行的遗传物质移动和重组密切相关。HGT的速率和方向通常受到MGEs之间复杂相互作用、宿主染色体背景以及外部选择压力的影响[19]、[20]、[21]、[22]。例如,克雷伯菌属中的高毒力和多重耐药性的出现已被证明是由质粒杂交和插入序列介导的染色体整合所导致的[23]。类似的现象也出现在其他细菌类群中,表明这些相互作用在选择压力下对细菌进化具有保守作用[24]。目前的研究尚未完全阐明由核心基因组系统发育和变异基因内容所塑造的基因组环境如何调节携带ARGs和VFGs的MGEs的整合、稳定性和转移效率。因此,确定可识别的MGEs与产氧克雷伯菌谱系中特定ARGs和VFGs组合之间的潜在关联模式具有重要意义。
在本研究中,我们专注于产氧克雷伯菌分离株,以提供一个基于分类学的框架来进行种内比较,认识到这种方法能够捕捉多个菌株的基因组内容。截至2026年4月,已有633个产氧克雷伯菌基因组在NCBI平台上完成测序并公开(见表S1),进一步通过平均核苷酸身份(ANI)进行比较,选择了ANI值≥95%的598个基因组用于物种水平的基因组分析。基于这些测序基因组进行了比较分析,以识别产氧克雷伯菌分离株中的核心基因、辅助基因和菌株特异性基因,并对其基因组结构进行了表征,以评估物种内的基因多样性。随后通过与KEGG和COG数据库的序列比对,对这些核心基因、辅助基因和菌株特异性基因进行了全面的功能注释。此外,还对所有产氧克雷伯菌基因组中的ARGs、VFGs和MGEs进行了广泛的鉴定和量化,并评估了特定类型MGEs与ARGs/VFGs谱型之间的潜在关联。这种基于产氧克雷伯菌基因组的综合方法旨在提供一个基于分类学的框架,以理解物种内的基因分化以及基因组结构、基因内容和水平转移如何共同驱动产氧克雷伯菌的适应性进化。

章节片段

识别产氧克雷伯菌分离株中的核心基因、辅助基因和菌株特异性基因

参考我们之前的研究[25]、[26]、[27],对基因库进行了计数和比较,以识别产氧克雷伯菌分离株的基因多样性。为了分析产氧克雷伯菌的基因组结构,我们从633个公开可用的基因组中选择了598个基于基因组质量指标和ANI值≥95%的高质量全基因组序列(使用FASTANI v1.34 [28]),并使用Roary v3.12.0 [29]进行处理

产氧克雷伯菌菌株中的核心基因、辅助基因和菌株特异性基因

从NCBI数据库下载了633个公开可用的产氧克雷伯菌基因组。基于这些基因组进行了ANI导向的分析,以确保物种水平的一致性,阈值设为≥95%。最终保留了ANI值≥98%的598个基因组用于后续的全基因组分析(图1A和表S2–S3)。
使用OrthoFinder进行同源性推断,并结合Roary进行层次聚类(使用95%的蛋白质身份阈值和核心定义)

讨论

基于598株产氧克雷伯菌的比较基因组分析揭示了该菌作为适应性病原体和环境生物体的基本遗传结构。将基因库划分为相对较小的稳定核心基因组(1239个基因)、较大的灵活辅助基因组(中位数约4373.5个基因)以及大量的菌株特异性基因(中位数约11个基因),为理解物种内的基因多样性提供了定量框架

结论

通过对598个产氧克雷伯菌基因组的分析,本研究发现其核心基因组包含1239个基因,而辅助基因组具有高度变异性(每个菌株平均4373.5个基因)。ARGs和VFGs的分布不均,MGEs的丰度与ARGs之间存在强相关性(R = 0.65,P < 2.2e–16),而与VFGs的相关性较弱(R = 0.29,P = 2.4e–13),表明MGEs是抗生素抗性传播的主要驱动因素,而在传播毒力方面的作用较为有限

CRediT作者贡献声明

韩海伟:撰写——初稿。王玉童:撰写——初稿。青云明:可视化、软件和方法学。魏冠云:数据管理。张贤:撰写——审稿与编辑、监督、资金获取、概念构思

数据可用性

本研究的数据可在文章中找到。

利益冲突声明

作者声明他们没有已知的可能会影响本文工作的财务利益或个人关系。

致谢

本研究得到了中国湖南省自然科学基金(2025JJ50123)、国家自然科学基金(32101368和32200676)以及湖南省研究生科学研究创新项目(CX20250393)的财政支持。

生物通微信公众号
微信
新浪微博


生物通 版权所有