多重耐药铜绿假单胞菌（MDR P. aeruginosa）群体结构、耐药组、毒力组与移动组的基因组解析及ARG–MGE关联

《Frontiers in Microbiology》：Genomic insights into population structure, resistome, virulome, and mobilome of multidrug-resistant Pseudomonas aeruginosa

　　

背景：铜绿假单胞菌（Pseudomonas aeruginosa）是引起严重医疗保健相关感染的主要机会性革兰阴性病原体。本研究旨在评估从公共数据库获取的铜绿假单胞菌（P. aeruginosa）基因组中群体结构、抗菌药物耐药（AMR）决定因子、毒力因子（VFs）与移动遗传元件（MGEs）之间的相互作用。方法：研究人员对301株来源于临床、与人类相关的多重耐药（MDR）P. aeruginosa基因组进行了基因组谱分析。综合基因组分析主要集中于使用已建立的生物信息学工具对序列型（STs）、抗菌药物耐药基因（ARGs）、VFs、质粒、插入序列（ISs）、整合子、原噬菌体和泛基因组进行表征。结果：该数据集呈现多样化的群体结构，其中MDR基因组由国际公认的高风险STs主导，尤其是ST621、ST277和ST235。MDR P. aeruginosa基因组携带更高的ARG负荷（主要为获得的氨基糖苷修饰酶、β-内酰胺酶基因、磺胺类耐药基因和消毒剂耐药决定因子）。此外，鉴定出更多数量的ISs、完整1类整合子、质粒重叠群和原噬菌体序列，支持MGEs在基因组进化和AMR传播中的关键作用。泛基因组分析揭示了一个广泛的附属基因组，突显了显著的基因组多样性。相比之下，固有耐药组和毒力组保持高度保守。结论：MDR基因组在保守的固有耐药背景下表现出更高的获得性ARGs负荷。该整合基因组框架为与MDR P. aeruginosa相关的基因组特征提供了见解，并支持对新出现的MDR谱系及其相关移动耐药元件的基因组监测。

研究背景：铜绿假单胞菌（Pseudomonas aeruginosa）是临床重要的机会性革兰阴性病原体，可引起肺炎、血流感染、尿路感染和慢性烧伤创面感染等医疗保健相关感染，在囊性纤维化、免疫受损及留置医疗器械患者中危害尤为突出。该菌基因组较大、GC含量高，具有保守核心基因组和扩张的附属基因组，可通过水平基因转移（HGT）获得耐药与毒力性状。与主要依赖获得性耐药的大多数细菌不同，P. aeruginosa具有广泛的固有耐药组（intrinsic resistome），包括染色体编码外排泵、外膜通透性降低和诱导型β-内酰胺酶等；获得性超广谱β-内酰胺酶、碳青霉烯酶、氨基糖苷修饰酶及消毒剂耐药决定因子进一步导致多重耐药（MDR）、广泛耐药（XDR）甚至泛耐药（PDR）。已有全基因组测序（WGS）研究多聚焦单一层面，同时整合耐药基因型-表型、毒力基因谱和移动组（mobilome）的综合性分析仍相对有限。因此，研究人员对301株人源临床MDR P. aeruginosa基因组开展综合基因组分析，评估群体结构、泛基因组、固有与获得性耐药组、毒力组和移动组景观，以理解MDR相关基因组特征及其适应机制。

研究开展与意义：研究人员基于公共数据库来源的301株MDR P. aeruginosa基因组，整合群体结构、泛基因组、耐药组、毒力组和移动遗传元件（MGEs）分析。结论表明，MDR基因组在高度保守的固有耐药背景下携带较高的获得性ARG负荷，并伴随丰富的插入序列（ISs）、完整1类整合子、质粒重叠群和原噬菌体；高风险ST621、ST277和ST235占主导；附属基因组庞大而核心毒力组高度保守。该研究为MDR P. aeruginosa的基因组监测、高风险谱系识别、移动耐药元件追踪及感染控制提供了基因组框架。

关键方法（不超过250字）：样本队列来源于NCBI Pathogen Detection数据库（访问日期2026年2月18日），经分类、宿主、采集时间和抗菌药物敏感性试验（AST）信息筛选，最终纳入301株MDR P. aeruginosa基因组。研究人员以QUAST、BUSCO、CheckM、FastANI和PathogenWatch进行组装质量、污染、物种一致性和表型判定；以PubMLST方案进行序列型（STs）预测并构建最小生成树；以Prokka和COGclassifier进行结构注释与功能分类；以ABRicate结合CARD和VFDB注释ARGs与VFs；以ISEScan、IntegronFinder、MOB-suite和Phigaro注释ISs、整合子、质粒和原噬菌体；以Flankophile分析ARG–MGE基因组上下文；以Panaroo、MAFFT和IQ-TREE进行泛基因组及核心基因组系统发育分析。

研究结果：
3.1 基因组数据集检索与质量评估：研究人员从数据库筛选并通过contig数、Ns/100 kbp、BUSCO完整性、CheckM完整性/污染和ANI等质量评估，保留544株高质量基因组用于后续分析。
3.2 抗生素耐药表型与群体结构：通过AST数据分类，544株中301株（55.3%）为MDR，15株为非MDR，228株为不确定；因磷霉素AST缺失，未可靠判定XDR/PDR。301株MDR基因组来自8个国家、2004—2023年，MLST鉴定出62个STs，以ST621（198株，65.8%）、ST277（23株，7.6%）、ST235（14株，4.7%）和ST973（3株，1.0%）为主。
3.3 结构基因组架构：基因组注释平均含6,460个基因、6,271个CDSs等；COG分类显示原噬菌体和转座子相关基因平均74个/基因组，提示MDR群体MGE负担较高。
3.4 耐药组：共鉴定16,705个ARGs、167个独特基因，平均56个/基因组。3.4.1 固有耐药组：84个基因条目、12,328次检出，aph(3^′)-IIb、arnA和fosA见于全部基因组，Mex/Opr外排系统广泛存在，bla_PDC–374检出283株（94.0%），bla_OXA–914检出198株（65.8%），显示高度保守。3.4.2 获得性耐药组：73个条目、1,673次检出，以bcr-1、sul1、qacEΔ1、aac(6^′)-Ib7和aadA常见；bla_OXA–2检出200株（66.4%），获得性碳青霉烯酶频率较低。3.4.3 调控元件：armR、basS、cprR、cprS、parR、parS、cpxR、soxR和rsmA等见于299—301株。3.4.4 耐药基因的谱系相关分布：ST621与非ST621间固有外排系统分布接近，但bla_OXA–914仅见于全部198株ST621；aac(6^′)-Ib7、aadA、bla_OXA–2、qacEΔ1和sul1在ST621中多呈高频率，非ST621中分布差异明显。
3.5 毒力组：VFs总体高度保守。3.5.1 定植、黏附与生物膜相关决定因子：IV型菌毛、Fap系统、藻酸盐合成和群体感应基因广泛存在。3.5.2 运动、分泌与宿主互作因子：鞭毛基因、III型/VI型分泌系统、HSI和Xcp系统高度保守，exoS、exoT、exoU和exoY分布不一，lasA、lasB、aprA、plcH和toxA广泛检出。3.5.3 免疫逃逸与宿主适应决定因子：LPS和鼠李糖脂相关基因高度普遍。3.5.4 营养获取与代谢毒力决定因子：pyochelin、pyoverdine和pyocyanin合成相关基因总体普遍，但部分pyocyanin基因存在变异。
3.6 移动组：3.6.1 插入序列：共6,332拷贝、平均21拷贝/基因组，IS3家族最多（46.7%），IS5、IS110、IS21和ISL3次之。3.6.2 整合子：251/301株（83.4%）携带整合子，共291个，其中245个为完整整合子，attC、intI和attI1常见。3.6.3 质粒：279/301株（92.7%）含质粒相关contigs，共563个，rep_cluster_322和rep_cluster_2238常见；96个（17.1%）预测为接合型。3.6.4 原噬菌体：全部301株检出，共1,444个区域、平均4.80个/基因组，Siphoviridae和Myoviridae为主，转座性原噬菌体133个。
3.7 ARG与MGE关联：3.7.1 主要获得性ARGs的基因组上下文与共定位：bcr-1常与bpeR相邻，qacEΔ1呈ant1–bla–qacEΔ1–folP复现排列，sul1常与ant1相邻并可伴IS6100，aac(6^′)-Ib7常邻近bla，aadA呈aadA–bla排列，bla_OXA–2常与ant1、aacA4和folP共现。3.7.2 整合子相关耐药决定因子及基因盒组成：aac(6^′)-Ib、qacE、ant(3^″)-Ia/aadA1和bla_OXA–2等基因盒常见，提示1类整合子参与多药耐药基因盒累积。3.7.3 质粒相关抗菌药物耐药基因：362/563个质粒contigs（64.3%）含至少1个耐药决定因子，常见cpxR、qacEΔ1、sul1、aac(6^′)-Ib7、bla_OXA–2和aadA；67个接合型、1个可移动型，部分接合性质粒携带多类耐药决定因子。
3.8 泛基因组：共14,741个基因簇，核心基因组5,124个（34.8%），附属基因组9,617个（65.2%），包括soft-core、shell和cloud基因，显示显著基因组多样性和附属基因变异。

讨论总结：研究人员认为，MDR P. aeruginosa群体以ST621、ST277和ST235等高危克隆为主；耐药组表现为保守固有耐药框架与通过HGT获得的氨基糖苷修饰酶、β-内酰胺酶及1类整合子相关ARGs的叠加；毒力组高度保守，提示成功MDR谱系在获得耐药性的同时保留致病潜力；ISs、完整整合子、质粒和原噬菌体等MGEs促进基因组进化和AMR传播；获得性ARGs常位于复现的耐药相关基因组上下文中，部分位于具接合潜力的质粒contigs上；泛基因组显示大附属基因组和保守核心。研究局限包括未评估染色体突变相关耐药机制、防御系统、HGT预测和重金属耐药等，且ST621及美国和巴西来源占比高，基因组预测缺乏功能验证或转录组分析。

研究结论（翻译）：本研究对301个MDR P. aeruginosa基因组进行了综合基因组表征，整合群体结构、泛基因组结构、耐药组、毒力组和移动遗传元件。研究结果表明获得性抗菌药物耐药决定因子负担较大，并伴随包括质粒、整合子、插入序列和原噬菌体在内的多种移动遗传元件。ARG–MGE关联分析进一步在质粒相关contigs和复现的基因组上下文中鉴定出耐药决定因子，包括涉及主要获得性ARGs的保守排列，突显移动和基因组上下文对AMR持续和传播的潜在贡献。数据集中识别出国际公认的高风险克隆，进一步强调成功耐药谱系的循环。总体而言，本研究为MDR P. aeruginosa的基因组景观和AMR移动相关特征提供了见解，并为基因组监测和分子流行病学建立了有价值资源。持续基于WGS的监测，尤其整合ARG、MGE和基因组上下文信息，对于监测耐药出现和传播并支持感染控制和抗菌药物管理策略具有重要意义。