背景：慢性伤口感染因其多重感染（polymicrobial）性质及日益加重的抗菌素耐药性（Antimicrobial Resistance, AMR）负担，构成重要的临床与公共卫生挑战。常规培养法诊断常无法全面捕获感染相关的微生物多样性及耐药潜能。方法：收集巴基斯坦某三级医院病程超过15天愈合的慢性伤口标本，按解剖部位分为小腿（ll-H1）、大腿（ul-H2）、足部（ft-H3）、胸部（ct-H4）及留置导管（ca-H5）五组，采用鸟枪法宏基因组测序（Shotgun Metagenomic Sequencing）联合常规培养方法表征微生物群落、抗菌素耐药基因（Antibiotic Resistance Genes, ARGs）及毒力决定簇（Virulence Determinants），并进行分类学和功能谱分析以评估不同伤口亚组的微生物多样性与耐药模式。结果：宏基因组分析显示变形菌门（Proteobacteria）、拟杆菌门（Bacteroidetes）和放线菌门（Actinobacteria）占优势；频繁检出临床相关病原菌，包括木糖氧化无色杆菌（Achromobacter xylosoxidans）、金黄色葡萄球菌（Staphylococcus aureus）和铜绿假单胞菌（Pseudomonas aeruginosa），以及较少报道的菌种如异样无色杆菌（Achromobacter insolitus）和嗜麦芽窄食单胞菌（Stenotrophomonas maltophilia）；鉴定到多重ARGs簇及生物膜相关毒力因子，提示显著的多药耐药（Multidrug Resistance, MDR）潜能；部位特异性分析显示铜绿假单胞菌在ul-H2中占优（约32%），产酸肠杆菌（Enterobacter hormaechei）在ft-H3中最丰富（约40%）；培养法主要恢复常见好氧病原菌，而宏基因组检测到额外机会致病菌及难培养类群，凸显常规诊断局限性；耐药基因组分析鉴定到赋予β-内酰胺类、氨基糖苷类、氟喹诺酮类、四环素类及大环内酯类耐药的ARGs。结论：巴基斯坦慢性伤口感染携带具实质AMR与毒力潜能的多样多重微生物群落；鸟枪法宏基因组较培养法能更全面地表征各伤口部位的附加病原体与耐药决定簇；上述发现支持将宏基因组诊断整合入临床决策，以加强资源受限医疗环境中的抗菌药物管理（Antimicrobial Stewardship）及感染控制策略。
本文对发表于《Journal of Infection and Public Health》的研究论文《慢性伤口多重感染中的抗菌素耐药性与毒力特征：基于鸟枪法宏基因组（Shotgun Metagenomics）的视角》进行解读总结。
【研究背景】
抗菌素耐药性（Antimicrobial Resistance, AMR）是全球公共卫生紧急事件，每年导致约495万例死亡，中低收入国家负担最重。慢性伤口感染因持久性微生物定植、生物膜形成及多重感染（Polymicrobial Infection）特征，临床治疗困难。传统培养法仅能分离约1%的环境微生物，无法检测苛求菌、活的不可培养（Viable but Non-Culturable, VBNC）状态菌及多数厌氧或慢生长菌，且不能全面揭示耐药基因组（Resistome）与毒力因子（Virulence Factors, VFs）。巴基斯坦地区多重耐药（Multidrug Resistance, MDR）与广泛耐药（Extensively Drug-Resistant, XDR）细菌感染率上升，但针对慢性伤口感染的全面宏基因组研究仍较缺乏。因此研究人员旨在借助鸟枪法宏基因组测序（Shotgun Metagenomic Sequencing）联合培养法，系统解析巴基斯坦三级医院术后慢性伤口感染的微生物组成、ARGs及VFs，并探讨解剖部位特异性差异。
【研究目的与总体结论及意义】
研究人员采集50例愈合延迟超15天的术后慢性伤口标本，按解剖部位合并为小腿（ll-H1）、大腿（ul-H2）、足部（ft-H3）、胸部（ct-H4）、留置导管（ca-H5）五个混合样本进行Illumina MiSeq鸟枪法宏基因组测序，同时平行做需氧培养及16S rRNA鉴定。研究证实伤口微生态以变形菌门（Proteobacteria）、拟杆菌门（Bacteroidetes）、放线菌门（Actinobacteria）为主，检出多种临床重要及少见机会致病菌；宏基因组比培养法多检出木糖氧化无色杆菌（Achromobacter xylosoxidans）、异样无色杆菌（Achromobacter insolitus）、嗜麦芽窄食单胞菌（Stenotrophomonas maltophilia）、普通变形杆菌（Proteus mirabilis）、洋葱伯克霍尔德菌相关种（Bordetella trematum）等。Resistome分析发现130种独特ARGs，涵盖β-内酰胺酶（含blaNDM-1等碳青霉烯酶）、氨基糖苷修饰酶（aac(6′)、aph(3′)-IIIa等）、氟喹诺酮靶位突变基因（gyrA、parC）、四环素（tetA等）、磺胺类（sul1/sul2）、大环内酯（erm、msr(E)/mph(E)）等，部分ARGs共定位于同一Contig提示可移动遗传元件介导水平基因转移（Horizontal Gene Transfer, HGT）。VF分析鉴定到铜绿假单胞菌藻酸盐生物合成（algD、algG等）、IV型菌毛（pilG/H/I）、III型分泌系统（Type III Secretion System, T3SS：popB、popD等）、金黄色葡萄球菌胞间黏附/荚膜（ica/cap8）及鞭毛（fliE、flgC-G等）相关基因，表明强生物膜形成与免疫逃逸能力。部位分析示ll-H1与ul-H2 ARG负荷最高（MDR热点），ft-H3与ct-H4 ARG多样性较低；ca-H5多态微生物负荷最高并伴消毒品耐受基因。研究表明鸟枪法宏基因组能弥补培养局限，为MDR伤口感染精准诊疗、抗菌药物管理及"全健康（One Health）"框架下AMR监测提供实证依据。
【主要关键技术方法】
研究人员自2022年11月至12月于巴基斯坦拉合尔本娜齐尔·布托医院骨科及普外科病房采集50例术后慢性伤口组织拭子/去痂组织，按解剖部位各10例混池为ll-H1（小腿/腓肠）、ul-H2（大腿/股部）、ft-H3（足）、ct-H4（胸）、ca-H5（留置导管周围）共5个混合样本。同步行常规需氧培养（血琼脂、巧克力琼脂、营养琼脂、甘露醇高盐琼脂、胰蛋白大豆琼脂Tryptic Soy Agar, TSA）及分离株16S rRNA Sanger测序鉴定。混合样本用DNeasy PowerSoil Kit提取总DNA，Covaris打断建库，Illumina MiSeq 2×150 bp双端测序。原始读段经FastQC质控、SOAPnuke去接头及低质量过滤、KneadData比对去除人源序列（GRCh38.p14）、BBMap去duplicate。Clean reads用Kraken2+Bracken基于MiniKraken数据库做物种分类与相对丰度估算；用MEGAHIT（k-mer=29~141）拼接Contig（>200 bp），PROKKA注释，CARD（Comprehensive Antibiotic Resistance Database）与ResFinder注释ARGs，VFDB注释毒力因子；α多样性（Chao1、Shannon）及β多样性（Jaccard、Bray–Curtis）经phyloseq/vegan包在R中计算并做主坐标分析（Principal Coordinate Analysis, PCoA）；ARG共现及同一Contig定位用自制脚本分析，统计用SPSS（p<0.05为显著）。
【研究结果】
3.1 微生物组（Microbiome）分析
通过分类学 profiling 发现所有样本变形菌门（Proteobacteria）平均相对丰度最高，其次为拟杆菌门（Bacteroidetes）和放线菌门（Actinobacteria）。ll-H1以沙雷氏菌属（Serratia sp. SSNIH 32%）、木糖氧化无色杆菌（Achromobacter xylosoxidans 15%）、鲍曼不动杆菌（Acinetobacter baumannii 20%）为主；ul-H2铜绿假单胞菌（Pseudomonas aeruginosa ~32%）占优；ft-H3产酸肠杆菌（Enterobacter hormaechei ~40%）与鲍曼不动杆菌（32%）为主；ct-H4产酸肠杆菌（29%）与鲍曼不动杆菌（23%）为主；ca-H5多态性最强，铜绿假单胞菌（22%）、洋葱伯克霍尔德菌群（Bordetella spp. 19%）、气味类香味菌（Myroides odoratimimus 16%）、普通变形杆菌（12%）、肺炎克雷伯菌（Klebsiella pneumoniae 11%）共存。培养法仅获得葡萄球菌属（Staphylococcus aureus、Staphylococcus cohnii等）、产酸肠杆菌及少量污染菌，与宏基因组共有检出仅为金黄色葡萄球菌，证明培养法低估革兰阴性及苛求菌。α多样性显示ca-H5 Chao1估计值最高（~2000），ll-H1 Shannon指数最高（3.4）；PCoA显示ft-H3与ca-H5群落结构明显分离，ll-H1与ul-H2相近，证实解剖部位塑造特异微生态。
3.2 伤口感染中抗菌素耐药基因（ARGs）流行情况
CARD注释共识别130个独特ARGs。各部位ARG相对丰度排序：ca-H5(33.7%) > ul-H2(27.9%) > ct-H4(20.9%) > ll-H1(10.5%) > ft-H3(7.0%)。β-内酰胺类（含ESBL blaCTX-M、TEM、SHV及碳青霉烯酶blaNDM、blaOXA）、氟喹诺酮类（gyrA/parC突变型及qnr家族）、氨基糖苷类（aac(6′)-Ib-cr、aph(3″)-Ib、aph(6)-Id、rmtB）、四环素（tetA/B）、磺胺（sul1/sul2）、大环内酯（erm、msr/mph）广泛分布；ca-H5磺胺与消毒剂抗性基因富集明显。层次聚类示ll-H1与ul-H2因共享β-内酰胺/氨基糖苷/氟喹诺酮 ARGs 聚为一支，ft-H3与ct-H4因较低ARG多样性聚为另一支。
3.3 ARGs共现与关联（Cooccurrence and Association of ARGs）
Contig组装发现多位点ARGs物理连锁：ll-H1同Contig携带多个blaNDM家族（blaNDM-1/2/3/18/23/24）；ul-H2同Contig携aph(6)-Id–aph(3″)-Ib及大环内酯外排基因msr(E)-mph(E)；ft-H3同Contig携aac(6′)-Ib-cr与blaOXA-1（β-内酰胺+氨基糖苷双重耐药）；ct-H4单Contig跨7类ARGs（含aadA、blaTEM-1B、rmtB）；ca-H5 qnrD1与磺胺/四环素耐药基因共定位。表明质粒、转座子等可移动遗传元件驱动多药耐药共选择（Co-selection），ll-H1与ul-H2为MDR基因热点区。
3.4 差异丰度ARGs与细菌的相关系数
Pearson相关性显示ul-H2与ca-H5 ARG谱正相关最强（r=0.76），ll-H1与ul-H2次之（r=0.55），ft-H3与ct-H4呈中等正相关（r=0.43），提示相似选择性压或交叉污染可能。
3.5 毒力因子（Virulence Factor, VF）丰度
VFDB比对示铜绿假单胞菌藻酸盐合成基因簇（algU/K/D/G/Q/8/Z等）、IV型菌毛（pilG/H/I）、T3SS效应分子（exsC/popB/popD/pscJ等）、鞭毛组装基因（fliE/fliG/flgC-G）、金黄色葡萄球菌荚膜（cap8）及胞间黏附（icaADBC）基因广泛存在。部位VF hits数：ul-H2(238, 35.68%) > ct-H4(202, 30.28%) > ft-H3(143, 21.44%) > ll-H1(61, 9.15%) > ca-H5(23, 3.45%)，提示上腿伤口毒力潜能较高。
3.6 与临床元数据的相关性
年龄与感染部位（r=0.44）、糖尿病（r=0.37）正相关，与免疫抑制负相关（r=-0.46）；免疫抑制与感染持续时间弱正相关（r=0.23），符合慢创与基础病/免疫状态关联特征。
【讨论与结论翻译总结】
讨论指出变形菌门主导符合慢性伤口"微生物特征（microbial signature）"，培养盲区由宏基因组弥补；检出的多种葡萄球菌、摩根氏菌（Morganella morganii）、脆弱拟杆菌（Bacteroides fragilis）均与生物膜相关难治感染有关；导管相关样本中aac(6′)及sul基因富集反映器械相关选择性压；ARG共定位Contig证实HGT在院内AMR播散中作用。解剖部位差异影响Resistome构成，小腿/大腿（ll-H1/ul-H2）高ARG负荷可能与反复干预及抗生素暴露相关。
结论（翻译）：本研究通过对五级解剖部位慢性/术后伤口样本的宏基因组与培养联用评估，揭示高度多样且具临床意义的ARGs与毒力决定簇储库；各部位优势菌不同——ll-H1以沙雷氏菌属（32%）及无色杆菌属（合计35%）为主，ul-H2铜绿假单胞菌（32%）占优，ft-H3产酸肠杆菌（40%），ct-H4产酸肠杆菌（29%），ca-H5多态性最强（铜绿假单胞菌44%、Bordetella spp.19%、Myroides odoratimimus 16%等）。鸟枪法宏基因组实现未培养类群深度剖析，广谱ARGs（氨基糖苷修饰酶aac(6′)-Ib、磺胺sul1/sul2、碳青霉烯酶、头孢菌素酶、四环素及大环内酯耐药基因）与生物膜-VF簇共存提示感染迁延与治疗难度。16S rRNA验证到多重生物膜形成葡萄球菌及MDR肠杆菌科。研究强调需改进诊断、制定靶向治疗并建立结构化抗菌药物管理。院内需强化术后护理与感染控制，并将宏基因组监测纳入常规，以遏制MDR病原播散。未来需扩大样本量开展个体水平宏基因组研究验证上述发现。
