《mSphere》：From gene detection to resistome ecology: targeted hybrid capture in AMR surveillance
抗菌素耐药性（Antimicrobial Resistance, AMR)监测需要能够超越基于培养的方法和靶向分子诊断、全面捕获完整耐药组(Resistome)的方法。本研究中，研究人员比较了商用实时聚合酶链式反应(Real-time PCR)检测碳青霉烯酶的方法与靶向杂交捕获下一代测序技术(xHYB, Targeted Hybrid Capture Next-Generation Sequencing)在直肠拭子样本中进行全面抗菌素耐药基因(Antimicrobial Resistance Gene, ARG)分析的效果。xHYB检测板与PCR检测在*bla*KPC和*bla*NDM检测方面显示出高度一致性，但总体一致性因基因靶点不同而有所差异。测序方法表现出更高的分析灵敏度，能够检测到PCR未能识别的额外ARG，特别是OXA样变体。测序读数计数与PCR循环阈值(Cycle Threshold, CT)呈负相关，支持两种方法之间存在半定量关系，而读数计数与总DNA浓度无关联。超越单基因检测，耐药组分析揭示了碳青霉烯酶阳性与阴性样本之间显著的组成差异。阳性样本显示出更广泛的多重耐药谱，包括β-内酰胺酶、氨基糖苷类修饰酶、四环素、磺胺类、大环内酯类、甲氧苄啶耐药基因以及多重耐药外排系统。跨耐药数据库的比较分析显示，基因和肽标记物水平存在变异，而基因家族水平的谱图更为一致。总体而言，靶向xHYB测序能够灵敏检测临床相关ARG，并揭示与碳青霉烯酶携带相关的结构化耐药组模式，支持其整合到AMR监测策略中。
抗菌素耐药性(Antimicrobial Resistance, AMR)已成为21世纪全球最紧迫的公共卫生威胁之一，2019年全球因细菌AMR导致的死亡人数估计达495万，其中127万例直接归因于耐药菌感染。尽管基于培养的传统诊断方法和靶向分子检测在临床中广泛应用，但前者耗时长、可能低估低细菌负荷或难培养微生物的定植情况，后者则受限于预设基因靶标，难以捕捉新型或罕见变异。更重要的是，现有方法多聚焦于单一或少数耐药基因的检测，无法反映微生物群落中耐药组(Resistome)的整体复杂性和生态结构。人类肠道作为多重耐药菌的主要储存库，其耐药基因的分布与传播对感染控制和公共卫生监测至关重要。直肠拭子作为一种便捷采样方式，可直接反映胃肠道微生物群落状态，但基于高通量测序的深度分析仍面临技术挑战：常规鸟枪法宏基因组测序虽然在理论上全面，但对于复杂样本中低丰度ARG的检测灵敏度不足，尤其在低生物量样本中。
在此背景下，研究人员开展了一项比较研究，旨在评估靶向杂交捕获测序技术(xHYB)相对于传统PCR检测在AMR监测中的性能，并系统刻画直肠拭子样本的耐药组特征。研究发现发表于《mSphere》，研究证实靶向xHYB测序不仅能高灵敏度检测临床相关ARG，还能揭示与碳青霉烯酶(Carbapenemase)携带相关的结构化耐药组模式，为整合测序技术到AMR监测体系提供了重要依据。
研究采用的主要关键技术方法包括：样本来源于120份用于常规AMR监测的直肠拭子标本；使用AB Analitica商用实时PCR试剂盒进行碳青霉烯酶基因的筛选与鉴定；采用QIAseq xHYB AMR检测板进行靶向杂交捕获富集，随后通过AVITI LT测序平台进行2×150 bp双端测序；数据分析使用CLC Genomics Workbench v. 24.0中的"Find QIAseq xHYB AMR Markers"流程，耐药基因注释综合比对CARD、NCBI、ResFinder和ARGA四个数据库；统计方法包括Spearman秩相关分析、Wilcoxon秩和检验、置换多元方差分析(PERMANOVA)、主成分分析(PCA)和层次聚类等。
研究结果部分，"队列中碳青霉烯酶基因及分离菌种的分布"显示：120份样本中，PCR检测阳性74份(61.7%)，阴性44份(36.7%)，未确定2份。阳性样本中，*a*+*d*基因组(*bla*
KPC和/或*bla*
OXA-48-like)占67.6%，*b*基因组(金属β-内酰胺酶基因*bla*
VIM、*bla*
NDM和*bla*
IMP)占47.3%，AcOXA(Acinetobacter OXA-like基因)占17.6%。最常分离的菌种为肺炎克雷伯菌(*Klebsiella pneumoniae*)，占阳性样本的63.5%。在鉴定板中，*bla*
KPC确认40份样本(平均CT 24.34±5.33)，*bla*
OXA-48-like基因9份(平均CT 24.86±4.65)，*bla*
NDM 11份(平均CT 24.75±7.00)，*bla*
VIM 13份(平均CT 28.42±6.67)，鲍曼不动杆菌OXA样基因3份(平均CT 25.80±3.10)。商检面板显示*a*+*d*基因高度流行，主要与肺炎克雷伯菌相关，而鉴定分析确认*bla*
KPC为最常见的碳青霉烯酶基因。
"xHYB检测板跨数据库的耐药基因分布"部分显示：通过xHYB检测板评估ARG时，研究人员从四个主要数据库获取注释信息。CARD数据库表现出最高的注释深度，识别出632个不同肽标记物、328个基因、129个基因家族和332个基因ARO注释；NCBI次之，有475个肽标记物、106个基因、39个基因家族和105个基因ARO注释；ResFinder和ARGA的检测能力更为有限，ARGA在所有类别中计数最低。表型ARO注释在所有数据库中均有限(范围5-13)。跨数据库比较显示，基因水平的共享核心较小，大部分基因为数据库特异性；但在基因家族水平，重叠显著增加，主要耐药基因家族(APH、*dfr*、OXA、*tet*、*van*等)在多个数据库中共享。Spearman相关分析证实基因数量与肽标记物丰度在各数据库内呈正相关，其中CARD(*ρ*=0.71)、NCBI(*ρ*=0.75)和ResFinder(*ρ*=0.72)显示强相关，而ARGA相关较弱且未达显著性(*ρ*=0.35，*P*=0.32)。
"xHYB检测板与AB Analitica商用检测法的碳青霉烯酶基因检测比较"部分显示：对118份可评估样本中四种主要碳青霉烯酶基因家族(*bla*
KPC、*bla*
NDM、*bla*
OXA和*bla*
VIM)的分析表明，测序读数丰度与CT值之间存在清晰关系，高读数计数通常对应低CT值。特别是*bla*
OXA家族显示出最宽的读数分布范围和最大的亚型多样性，部分高测序信号样本对应PCR阴性或高CT值，反映了xHYB能检测到PCR未靶向的多种*bla*
OXA变体。诊断性能指标分析显示：*bla*
KPC检测中，NGS与ABA的一致性极佳，阳性百分一致性(PPA)为0.92，阴性百分一致性(NPA)为0.97，Cohen's kappa系数为0.90；*bla*
NDM-1检测中，PPA为0.92，NPA为0.98，kappa系数为0.87；*bla*
OXA-48检测的一致性较低，PPA为1.00但NPA仅0.80，kappa系数为0.390，主要由于NGS检测到21例假阳性；*bla*
VIM检测的kappa系数为0.600，属中等一致性。Spearman相关分析显示，仅*bla*
KPC的CT值与测序读数之间存在统计学显著的中度负相关(*ρ*=-0.468，*n*=41，*P*=0.002)，而*bla*
NDM、*bla*
OXA-48-like和*bla*
VIM的相关性均未达显著性。DNA浓度分析显示其在PCR阳性与阴性组间分布重叠，表明总DNA产量不影响碳青霉烯酶耐药基因的检测。
"碳青霉烯酶阳性与阴性样本的耐药组结构差异"部分显示：基于合并耐药表和单个数据库的基因水平丰度矩阵，PCA揭示了ABA阳性与阴性样本在PC1(16%方差)和PC2(9.1%方差)上的明显分离趋势。PERMANOVA证实两组间整体耐药组组成存在显著差异(*P*=0.001，*R*
2=0.042)，层次聚类以较高的轮廓分数(k2: 0.18)和聚类纯度(0.63)重现了ABA分类。对PCA分离贡献最大的基因主要与多重耐药机制相关，包括外排系统(Mex、Opr和Mux家族)、β-内酰胺酶(*bla*
OXA、*bla*
SHV、ADC和*ampH*)、氨基糖苷类修饰酶(APH和AAC)、四环素(*tet*)、磺胺类(*sul*)、大环内酯类(*mph*)和甲氧苄啶(*dfr*)耐药基因。这一分离模式在独立分析CARD(*P*=0.001，*R*
2=0.040)、NCBI(*P*=0.001，*R*
2=0.045)、ResFinder(*P*=0.001，*R*
2=0.037)和ARGA(*P*=0.001，*R*
2=0.070)数据库时均观察到。基因水平层次聚类显示ABA阳性样本倾向于聚类在一起，具有更广泛和更高的基因丰度模式，而许多ABA阴性样本形成基因表达明显降低的独立聚类。差异丰度分析鉴定出大量在两组间丰度显著不同的基因，表明组间差异是定量和组成性的，而非仅限于特定碳青霉烯酶基因的有无。
讨论部分，研究人员首先强调AMR对人类健康构成重大威胁，全球卫生组织已认可行动计划以维持抗感染治疗的有效性，而监测在制定政策和控制策略中发挥关键作用。现有监测主要基于ARG的流行率，但主动监测常依赖细菌培养(需24-72小时，可能低估定植)或局限于预设基因的靶向分子检测。胃肠道作为MDR细菌的主要储存库，肠道定植与后续感染风险、医疗环境传播及ARG的水平 mobilization相关。在该研究背景下，靶向杂交捕获测序技术为克服传统方法局限提供了新途径。
研究人员指出，xHYB检测板与商用PCR检测(ABA)相比具有高阳性一致性，尤其在*bla*
KPC检测方面表现优异，证明了其在复杂生物基质中识别临床关键碳青霉烯酶的强性能。阴性一致性较低主要源于xHYB阳性/ABA阴性的额外检出，特别是*bla*
OXA-48相关检测，这归因于xHYB能识别多种非*bla*
OXA-48-like的*bla*
OXA-like变体，而PCR检测未能靶向这些变体。值得注意的是，DNA浓度未影响两种方法的性能，且*bla*
KPC的CT值与测序读数之间存在负相关。
该研究明确提出了靶向杂交捕获AMR监测中实用读数阈值的定义问题。基于观察到的读数分布，1-10读数范围内的低丰度检测占相当比例，常与微生物培养阴性结果相关，可能反映背景噪声或生物学无关检测。因此，≥10读数的阈值被采纳为增强可解释性的保守截断值。然而，研究人员强调固定阈值在宏基因组AMR分析中仍是未解决的方法学问题，读数受测序深度、探针杂交效率、微生物负荷和基因拷贝数等多因素影响，未来需通过更大数据集、独立验证队列及基因特异性或概率模型来优化阈值策略。
数据库差异分析显示，基因家族水平的聚类产生更一致的模式，而精确基因变体识别可能因数据库不同而受损。整合多个数据库可提高稳健性，但需谨慎解释以平衡灵敏度和特异性。研究人员同时指出xHYB效率的靶点变异性可能引入捕获偏倚，影响基因丰度的定量估计；该方法对探针设计参考数据库具有内在依赖性，高度分歧序列可能因捕获效率降低而被低估，完全新颖的耐药基因或机制仍需通过非靶向方法识别。
超越单基因检测的核心发现是：ABA阳性状态反映的不仅是碳青霉烯酶决定因子的存在，而是耐药基因库的更广泛重组，对应于多种耐药机制协调累积的明确MDR状态。多重外排系统(Mex、Acr和Mux)、β-内酰胺酶(*bla*
OXA、*bla*
SHV和ADC)、氨基糖苷类修饰酶(APH和AAC)以及对四环素、磺胺类、大环内酯类和甲氧苄啶的耐药决定因子对组间区分贡献最大。这一分离在不同数据库中的可重复性支持了信号的生物学稳健性，表明观察到的差异反映一致的耐药组结构而非数据库特异性假象。
研究人员进一步将人类耐药组置于更广泛的生态和进化框架中讨论。人类耐药组是嵌入微生物组中的复杂动态遗传生态系统，涵盖共生菌和致病菌分布的基因，代表一个古老且生态结构化的ARG储存库。耐药组组成因解剖部位、抗生素暴露、炎症和人体测量因素而异。碳青霉烯酶检测状态与耐药组组成的关联提示，其可作为更广泛MDR基因组背景的间接标志物，ABA阳性相关菌株可能代表耐药机制随时间协调累积的进化稳定MDR谱系。然而，缺乏相关临床和生态数据限制了确定这些耐药组转变和ARG获得事件驱动因素的能力。
在移动遗传元件方面，质粒介导的水平基因转移是ARG传播的主要驱动力。质粒可在选择压力下跨物种和生态位移动，作为ARG mobilization的高效载体。因此，物种特异性关联可能具有误导性，并阻碍新型分子或培养检测策略的开发。疾病相关微生物组常显示扩展的耐药组，尽管分类学变化温和，表明耐药变化可能发生在菌株水平或通过移动ARG的获得而不伴随主要组成改变。这一认识支持更广泛的"泛耐药组"(Pan-resistome)概念：以ARG为中心而非物种为中心的方法更为适当，监测策略应优先追踪高度移动且临床相关的ARG及其载体，而非仅关注致病菌种。
最后，研究人员讨论了将宏基因组测序信号转化为临床可解释诊断结果的固有限难。尽管NGS能够实现耐药决定因子的全面检测，但测序数据的定量性质在定义阳性阈值时引入重要分析考量。固定读数截断值与参考分子检测表现出良好一致性，但基于阈值的分类策略可能需要优化或基因特异性校准。未来工作应聚焦于整合额外参数以改善宏基因组AMR检测的可解释性和诊断性能。
研究结论部分指出：该研究证明靶向xHYB测序能够灵敏检测临床相关ARG，并揭示与碳青霉烯酶携带相关的结构化耐药组模式，支持其整合到AMR监测策略中。研究人员强调，虽然xHYB在识别预设靶标方面具有优势，但其对参考数据库和探针设计的依赖性意味着完全新颖的耐药机制仍需非靶向方法发现。将宏基因组测序整合到临床和公共卫生AMR监测中，有望提供更全面的耐药趋势图景，辅助早期识别高风险耐药模式，最终加强感染控制和抗生素管理策略。ogy) represent one of the most pressing health challenges of the 21st century. In this study, the integration of xHYB sequencing into AMR surveillance strategies demonstrated its capability to capture structured resistome patterns associated with carbapenemase carriage, offering a more comprehensive view of resistance gene diversity than conventional diagnostics and supporting its role in improving surveillance precision.
