《Frontiers in Cellular and Infection Microbiology》：Resistance-gene minority variation and delayed microbiological response in rifampicin-resistant tuberculosis
编辑推荐：
背景：常规全基因组测序（WGS）对结核分枝杆菌耐药性的解读主要关注共有水平突变。由等位基因频率（AF）量化的低频变异可能指示药物压力下的耐药亚群，但其与利福平耐药或耐多药结核病（RR/MDR-TB）治疗应答的关联仍不明确。
方法：研究人员回顾性分析了2021至
背景：常规全基因组测序（WGS）对结核分枝杆菌耐药性的解读主要关注共有水平突变。由等位基因频率（AF）量化的低频变异可能指示药物压力下的耐药亚群，但其与利福平耐药或耐多药结核病（RR/MDR-TB）治疗应答的关联仍不明确。
方法：研究人员回顾性分析了2021至2025年间在中国宁波治疗的RR/MDR-TB患者，整合了临床、表型、治疗、微生物学和WGS数据。低频耐药相关变异（RAVs）被定义为在AF≥3%时匹配世界卫生组织（WHO）V2突变目录条目的低频率变异，其中AF≥5%定义了主要的较高频率子集。按药物进行的分析排除了目标基因中位深度<200×的患者。采用协变量校正的逻辑回归评估与微生物学应答的关联，并次要检验了较高频率和治疗匹配的变异。
结果：研究人员分析了229例具有关联WGS、治疗和微生物学随访数据的RR/MDR-TB患者。固定RAVs由220/229例患者携带（中位数3，四分位距[IQR] 2–5）。在13例患者中检测到AF≥3%的低频RAVs。总体低频RAV负荷在校正后与微生物学应答无明确关联。最清晰的信号涉及治疗匹配的AF≥5%低频RAVs，其携带者的6个月培养阳性率高于非携带者（4/7 [57.1%] vs 40/198 [20.2%]；校正优势比[aOR] 7.21，95%置信区间[CI] 1.38–43.68；p = 0.020；q = 0.034），曲线下面积（AUC）从0.675适度增加至0.695。同时具有高固定RAV负荷和治疗匹配AF≥5%低频RAVs的患者6个月阳性率最高（4/5 [80.0%] vs 15/106 [14.2%]在低负荷非携带者参照组；Fisher精确检验p = 0.003）。
结论：低频变异显示出有限的总体预后价值，但当与固定RAV负荷共同考虑时，较高频率的治疗匹配变异为6个月培养阳性提供了连贯的信号。这些探索性发现支持对选定的低频变异进行谨慎的、基于治疗意识的审查，并且在临床使用前需要进行外部验证。
研究背景与科学问题
利福平耐药结核病（RR-TB）和耐多药结核病（MDR-TB）仍是全球结核病控制的主要障碍。尽管世界卫生组织（WHO）的最新指南强调改进诊断、缩短全口服方案和加强耐药监测的重要性，但即使在基线耐药谱相似的患者中，治疗应答仍然存在显著异质性。痰菌培养阴转虽然不是治愈的完美替代指标，但延迟的痰菌培养阴转与不良治疗结局相关，是一个有用的中间表型。全基因组测序（WGS）提高了结核分枝杆菌耐药性解读的分辨率，但大多数常规WGS解读仍聚焦于共有水平或固定突变，这可能导致一个潜在的盲区：耐药亚群可以存在于共有阈值以下，并在药物压力下变得具有临床意义。低频变异和异质性耐药在生物学上是导致治疗延迟应答的可能因素，但其解读非常困难，因为它们可能代表宿主内微进化、混合感染、新发耐药、谱系多样性或技术假象。因此，核心临床问题不仅是低频变异是否存在，而是它们是否出现在与治疗相关的遗传背景下。为此，研究人员开展了一项回顾性研究，旨在明确低频变异在何种条件下与RR/MDR-TB患者的微生物学应答相关联。
研究方法与技术路线
本研究是一项来自中国宁波市的回顾性队列研究，纳入了2021年至2025年期间通过市级结核病防治网络诊断和管理的RR/MDR-TB患者。最终分析队列包括229例具有关联WGS、治疗和微生物学随访数据的患者。关键技术方法包括：①基于Illumina NovaSeq 6000平台的双端150 bp测序进行全基因组测序（WGS）；②使用FreeBayes和LoFreq分别检测固定单核苷酸多态性（SNPs）和低频变异，并对低频变异设定了严格的质控标准（如最小AF 1%、交替读长深度≥10、链偏倚评分≤20）；③利用TB-Profiler和WHO V2突变目录对耐药相关突变（RAVs）进行注释和分级；④采用逻辑回归模型评估固定耐药基因负荷、低频变异特征与微生物学结局（2个月和6个月痰菌培养状态）之间的关联，并校正年龄、性别、糖尿病、基线涂片状态、治疗史和主要谱系等协变量。
主要研究结果
1. 分析队列与基因组变异图谱
从5232例痰菌培养阳性患者中筛选出328例RR-TB患者，经排除后最终纳入229例。大多数分离株属于谱系2（90.0%），以北京相关亚系2.2.1为主（82.5%）。在经典耐药基因中，共鉴定出846个固定RAVs记录，220例（96.1%）患者至少携带一个固定RAV（中位数3，IQR 2–5）。在13例患者中检测到AF≥3%的低频RAVs，其中11例达到AF≥5%阈值。所有14个治疗匹配的AF≥5%低频RAVs均经过人工核查确认，且均发生在已有耐药证据的分离株中。
2. 调整后的基因组关联与微生物学应答
在校正临床协变量后，固定RAV负荷与延迟培养阴转（定义为第2或第6个月培养阳性）仅显示边缘性名义关联（aOR 1.13，95% CI 1.00–1.27；p = 0.044），未通过多重检验校正。对于低频RAVs，任何AF≥3%的低频RAV与6个月培养阳性呈边缘性关联（aOR 3.24，95% CI 0.94–11.12；p = 0.057），而AF≥5%时关联显著（aOR 5.53，95% CI 1.46–22.97；p = 0.013；q = 0.032）。最清晰的信号来自治疗匹配的AF≥5%低频RAVs与6个月培养阳性的关联（aOR 7.21，95% CI 1.38–43.68；p = 0.020；q = 0.034），并通过Firth惩罚回归确认。此关联对6个月终点具有特异性，对2个月培养阳性或无不利结局无明确信号。
3. 耐药面板内的基因水平信号
探索性基因水平Fisher检验显示，固定pncA RAVs（68/229例患者携带）与2个月培养阳性（OR 2.72，95% CI 1.46–5.10；p = 0.0009；q = 0.010）和延迟培养阴转（OR 2.59，95% CI 1.39–4.83；p = 0.0019；q = 0.021）显著关联，并与6个月培养阳性呈名义趋势（OR 2.19，95% CI 1.11–4.31；p = 0.018；q = 0.199）。其他固定RAV基因在FDR校正后均无显著关联。
4. 低频变异的预测价值
将治疗匹配的AF≥5%低频RAVs加入包含临床协变量和固定RAV负荷的基线模型后，6个月培养阳性的AUC从0.675增至0.695（ΔAUC 0.021；似然比p = 0.020；q = 0.034），增量较小。
5. 耐药背景与敏感性分析
将治疗匹配AF≥5%低频RAVs与固定RAV负荷联合分析发现，两者兼具的患者6个月培养阳性率最高（4/5 [80.0%]），而低负荷非携带者仅为14.2%（15/106）（Fisher精确检验p = 0.003）。使用更广泛的WHO目录基因集进行的敏感性分析增加了低频变异检出率，但未增强治疗匹配的关联。作为特异性对照，非WHO目录中的耐药基因低频变异未显示任何结局关联。
讨论与结论
在本宁波RR/MDR-TB队列中，低频变异的临床信号并非由其总体负荷驱动。当变异频率较高、与治疗匹配并结合固定RAV负荷进行评估时，出现了最清晰的模式：治疗匹配AF≥5%低频RAVs与6个月培养阳性之间的关联通过了FDR校正（aOR 7.21，95% CI 1.38–43.68；q = 0.034），尽管仅基于7例携带者且置信区间较宽，AUC增益也有限（0.021）。这些结果支持对选定的低频变异进行解读性角色，而非将其视为独立的预测标志物。研究人员认为，这种模式可以通过治疗压力来理解，但数据并未证明克隆选择。共识水平的WGS仍然是耐药解读最具可重复性的基础，尤其是对于充分验证的固定突变。低频变异与6个月培养阳性的更强关联在临床上具有合理性，因为持续的6个月培养阳性更接近持续的微生物学无应答，而耐药亚群可能需要持续的药物暴露才能变得临床可见。本研究的局限性包括：单一队列回顾性设计、关键亚组稀疏（7例治疗匹配AF≥5%携带者）、估计值不精确、治疗暴露仅记录药物接收情况而非实际活性、以谱系2为主的队列限制了普适性、以及低频变异检测对分析流程的敏感性。结论：在经典耐药基因中，低频耐药相关变异虽不常见但具有生物学连贯性，其临床相关性取决于治疗背景而非总体负荷。治疗匹配的较高频率低频RAVs与持续性6个月培养阳性相关，且在固定RAV负荷较高的背景下最具信息价值。这一发现支持对选定的低频变异进行谨慎的、基于治疗意识的审查，但在任何基于低频变异的规则被考虑用于临床管理之前，需要严格的质量控制、人工核查和外部验证。
