铜绿假单胞菌是重要的机会致病菌，碳青霉烯耐药株已被世界卫生组织列为“高优先级”病原体，提示其对临床抗感染治疗与院内感染防控构成持续压力。碳青霉烯耐药通常由外膜通透性下降、外排泵增强或β-内酰胺酶（如AmpC）过表达等机制导致，但近年来经水平基因转移获得碳青霉烯酶的情况呈增加趋势，使耐药在不同地区与不同谱系间传播的风险上升。在产碳青霉烯酶铜绿假单胞菌中，blaVIM、blaIMP、blaKPC、blaNDM等基因较常见，常嵌入Ⅰ类整合子并位于染色体基因组岛内，部分也可由质粒或Tn125样转座子等结构携带，从而促进全球扩散；并且已有证据表明特定高风险克隆（如ST-235以及ST-111、ST-244、ST-357、ST-308等）与碳青霉烯酶获得密切相关。然而，与这些经典高风险克隆不同，ST-1047在既往研究中仅零星检出，却呈现blaIMP、blaVIM、blaNDM与blaDIM等多类碳青霉烯酶的区域性出现，监测亦提示其在南亚早期循环后，2022年冲突背景下又在欧洲多国的难民或撤离伤员中被发现，提示社会事件驱动的人群流动可能放大其跨国传播。尽管该谱系兼具“多碳青霉烯酶汇聚”与“欧洲出现增多”的特征，仍缺乏群体水平的系统解析来追溯其起源、演化与扩散路径，有必要通过群体系统发育与移动遗传元件结构证据，厘清“哪些分支在扩张、耐药基因装在什么载体上、扩散链条如何形成”，以支持早期预警与精准防控。

基于此，本研究收集全球141株ST-1047 基因组/分离株（覆盖四大洲、15 国）。图1显示该谱系近年样本增多，且以blaIMP-1为主，同时并存blaVIM、blaNDM、blaDIM等多类碳青霉烯酶，提示其具备“多酶型耐药基因汇聚并出现扩张”的高风险特征。因而后续分析将围绕两点展开：（1）识别是否存在主导传播的扩张亚克隆，以及（2）阐明这些碳青霉烯酶基因的载体与获得方式。针对此目标，作者首先基于BEAST2构建时间标定系统发育树，纳入83株ST-1047并以8株近缘ST-377作为外群，结果显示两谱系约相差425个SNP，最近共同祖先可追溯至约1875年。进一步地，图2在系统树上清晰识别出两条主导分支：SC2（n=59）囊括全部blaIMP-1/10携带株，SC1为blaVIM-11携带株的单系群，整体格局支持“单次水平获得后发生克隆扩张”。

为验证猜想，作者进一步转入“载体与结构层面”的证据链构建，依次解析关键亚克隆及其余碳青霉烯酶基因的获得方式。首先，针对SC1，图3的序列比对显示blaVIM-11并非零散获得，而是由一段约20.4 kb的Tn501-like转座子整合至染色体而来；同时，携带exoU的基因组岛（GI-exoU）在SC1中经历两次删除导致exoU缺失，并保留残留片段作为结构证据，提示SC1在起源阶段同时发生“耐药获得”与“毒力岛重塑”。随后，针对SC2，作者在图4中比较PAGI-97B-like基因组岛的不同结构型，发现其整体高度同源且仅有少量indel，blaIMP-1位于I类整合子In1595内，并由xerD相关重组介导整合入染色体，支持SC2的blaIMP获得符合“单次整合后克隆扩张”模式；同时，SC2内个别blaIMP阴性样本可由PAGI-97B-like结构变体解释，即在仍保留In1595/基因组岛骨架的背景下发生blaIMP-1的切除或缺失（对应原文Type B/C），从而形成“岛在而基因缺失”的表型。最后，为补齐谱系内其他碳青霉烯酶的来源机制，图5进一步表明blaDIM-1由约66 kb染色体基因组岛携带（含In1592），在549457/549448中插入gcvP，而不同国家来源的blaNDM-1周边序列缺乏同源性，提示其在ST-1047内为独立、反复的水平获得事件，从而揭示ST-1047既能形成优势扩张亚克隆（SC1/SC2），也具备持续捕获多类碳青霉烯酶的能力。

综上，本研究从全球尺度系统刻画了铜绿假单胞菌ST-1047的耐药演化与传播图景，核心结论是该谱系已具备“高风险克隆”的关键属性，其威胁并不止于携带碳青霉烯酶，而在于能够通过染色体整合型移动遗传元件稳定获得并维持耐药决定子，从而形成以blaVIM-11（SC1）和blaIMP-1/10（SC2）为标志的单系扩张分支，同时又可在不同地区反复、独立捕获blaDIM-1与blaNDM-1等碳青霉烯酶基因，体现出持续吸纳耐药基因的能力；在此基础上，作者将时间标定系统发育证据与样本来源背景相结合，提出社会事件驱动的人群流动可能放大该谱系跨国传播的公共卫生意义，从而为“冲突-医疗体系受损-院感扩散-跨境输出”的风险链条提供了可追溯的分子框架。本研究的创新性主要体现在：以群体水平数据整合与时间标定系统树为骨架，结合长读长闭环基因组对移动遗传元件进行结构解析，明确了SC1的Tn501-like介导blaVIM-11整合与GI-exoU删除、以及SC2的PAGI-97B-like（In1595）介导blaIMP染色体整合及其变体演化，从“载体-插入位点-谱系扩张”三位一体地解释了耐药扩张机制，从而超越了仅报告耐药基因出现的描述性研究。需要指出的局限性在于：样本在时间与地区上存在监测和采样偏倚，且部分基因组仅为短读长草图组装、缺少原始reads或长读长验证，使得某些耐药基因（尤其blaNDM）的载体完整结构与传播关系仍有不确定性；同时，“冲突背景促进传播”更多是基于时间系统树与流行病学线索的关联推断，缺乏跨国、前瞻性的临床与转运数据闭环支撑。未来研究可聚焦三点：扩大并均衡全球代表性采样并系统引入长读长测序以精确解析移动元件结构演化，将基因组监测与院感、患者流动及暴露信息联动以重建可检验的传播路径，并通过功能实验评估关键结构变异对适应性、毒力与传播力的影响，从而把结构关联提升为机制证据，为精准预警与防控提供依据。

图1  全球ST-1047型铜绿假单胞菌141株菌株/基因组的时间与地理分布概览。

图2  时间标定系统发育分析揭示ST-1047型铜绿假单胞菌的群体结构。

图3  SC1起源相关的遗传结构事件（blaVIM-11获得与exoU缺失）的序列比对证据。（A）（blaVIM-11获得）：以无插入参考序列（AP017302.1）为对照，成对比对显示SC1通过一段约20.4 kb的Tn501-like转座子整合获得blaVIM-11。（B）（exoU缺失）：以携带exoU的参考基因组（549457）为对照，比对显示SC1的exoU缺失涉及两次删除事件，且在同源区域仍保留该基因组岛的残留片段（remnants）。

图4  SC2中PAGI-97B-like基因组岛的结构变异及其携带blaIMP-1的分子证据。（A）对SC2分离株中不同PAGI-97B-like基因组岛变体进行成对序列比对，显示基因组岛主体区域高度同源（可达100%），差异主要表现为少量插入或缺失（indel）。（B）（关键功能区放大）：对基因组岛局部进行比对与标注，指出该岛通过XerD家族重组酶介导整合；blaIMP-1位于I类整合子In1595内。

图5  ST-1047中blaDIM-1与blaNDM-1的独立、反复获得事件（载体结构证据）。（A）blaDIM-1由一段约66 kb染色体基因组岛携带并插入gcvP基因；该岛包含携带blaDIM-1的I类整合子In1592，且与PAGI-18共享16.6 kb、>99%的连续同源区，并在5′端均被IS6100侧翼包夹。（B）（blaNDM-1）：对比巴基斯坦（PA20/PA70/PA78）与坦桑尼亚（FFDD_1011）样本中含blaNDM-1的contig，除blaNDM-1本体外几乎无同源性，提示其来自不同移动元件背景，属于独立获得。

原文链接：https://doi.org/10.1128/mbio.02020-25