不丹广泛传播的多药耐药结核分枝杆菌:一项回顾性基因组流行病学研究

Extensive Endemic Transmission of Multidrug-Resistant Mycobacterium tuberculosis in Bhutan: A Retrospective Genomic-Epidemiological Study

作者信息Thinley Dorji, Karchung Tshering, Lila Adhikari, Thinley Jamtsho, Pavitra Bhujel, Pema Lhaden, Sonam Wangchuk, Wytamma Wirth, Kristy Horan, Justin T Denholm, Norelle L Sherry, Michelle Sait, Timothy P Stinear, Benjamin P Howden, Patiyan Andersson
PMID41626299
发布时间2026-01-29
DOI10.1093/ofid/ofaf802

实验完整度

中

研究基于回顾性样本的全基因组测序和生物信息学分析,包括耐药突变检测、系统发育分析、聚类分析和统计关联分析,但缺乏功能验证和体内实验。

主要模型

203株结核分枝杆菌临床分离株 谱系2.2.1亚型 MDR-TB簇1 公共谱系2序列数据集

重点核对

SNP阈值(≤12 SNPs)用于定义近期传播簇 样本选择策略(每年约40% DR-TB和2.5% DS-TB) 测序平台为Illumina NextSeq 谱系2.2.1占主导(86.2%) 最大簇包含88% MDR-TB序列

摘要

背景:不丹的多药耐药结核病(MDR-TB)病例比例正在上升。我们进行了首次回顾性基因组流行病学研究,以了解不丹结核分枝杆菌的群体结构、耐药模式及近期传播情况。方法:对2018年至2022年间在澳大利亚墨尔本微生物诊断单位公共卫生实验室收集的不丹耐药结核病(DR-TB)和药物敏感结核病(DS-TB)分离株进行全基因组测序。进行生物信息学分析以识别耐药突变和病例的基因组聚类。结果:每年约对40%的DR-TB和2.5%的DS-TB进行测序。在通过质量控制的203个序列中,126个(62.1%)为MDR-TB,15个(7.4%)为异烟肼耐药结核病。共有4个不同的流行谱系,其中大多数属于谱系2(86.2%)。使用≤12个SNP的阈值,71%的序列形成12个基因组簇;最大的簇包含所有MDR-TB序列的88%,并覆盖了整个研究时期和全国。这些病例高度克隆,平均成对SNP距离为10(范围0-25)。与公开可用的国际序列数据的系统发育分析显示,该MDR-TB簇形成一个独特的分支。结论:与目前关于重复输入的假设相反,不丹MDR-TB的主要负担似乎是由于近期本地传播导致的一个大型地方性簇,主张针对该MDR-TB分支进行有针对性的加强接触者追踪和筛查。这项研究强调了在资源有限的环境中投资结核病基因组学以获得可操作信息来指导政策决策的重要价值。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
不丹MDR-TB的主要负担是源于反复的跨境输入还是本地传播?
核心机制
不丹的MDR-TB主要由单一引入事件后的地方性克隆传播所致,形成大型流行簇,而非多次输入。
主要证据
对203个结核分枝杆菌分离株进行全基因组测序,发现71%的序列形成12个簇,其中最大的簇包含88%的MDR-TB序列,且与邻国序列形成独特分支。
研究意义
研究结果挑战了先前关于反复输入的传统观点,强调加强接触者追踪和筛查对控制本地传播的重要性,并展示了在资源有限地区投资结核病基因组学的价值。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本收集与测序

获取代表不丹DR-TB和DS-TB的基因组数据。

从2018-2022年不丹国家结核病参考实验室的DR-TB和DS-TB培养阳性样本中随机选择一部分,送至澳大利亚进行全基因组测序。

2

生物信息学分析

鉴定耐药突变和系统发育关系。

使用fastp进行质量控制,TBProfiler检测耐药突变和谱系,Bohra管道进行核心比对和系统发育分析。

3

聚类分析

推断近期传播事件和基因组簇。

使用snp-dists计算成对SNP距离,采用单连锁层次聚类,以≤5和≤12 SNPs阈值定义簇。

4

统计分析和流行病学关联

评估与聚类相关的因素。

使用Fisher精确检验和多元logistic回归分析人口学因素和细菌因素与聚类的关联。

5

国际背景分析

将不丹序列与区域和全球序列进行系统发育比较。

将不丹谱系2序列与来自印度、中国、孟加拉国、越南和CRyPTIC研究的公开序列进行系统发育分析。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
7H9液体培养基----
Nextera XT制备试剂盒----
Illumina NextSeq平台Illumina--
GenoType MTBDRplus----
GenoType MTBDRsl----
Xpert MTB/RIF----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本选择
病例定义、随机选择比例(DR-TB约40%,DS-TB约2.5%)
阅读提示:Methods-Sample Selection
DNA提取与测序
DNA提取方法(珠磨法和醋酸钠沉淀)、文库制备试剂盒(Nextera XT)和测序平台(Illumina NextSeq)
阅读提示:Methods-DNA Extraction and Sequencing Methods
生物信息学分析
参考基因组(H37Rv)、比对最小比对百分比(95%)、最小深度(30x)、重复区域masking、系统发育树构建模型(GTR+G4)
阅读提示:Methods-Bioinformatic Analysis
聚类分析
SNP阈值(≤5和≤12)、聚类方法(单连锁层次聚类)
阅读提示:Methods-Cluster Analysis
统计关联
统计检验(Fisher精确检验、多元logistic回归)、变量(年龄、性别、地区、感染部位、病例类型、年份、谱系)
阅读提示:Methods-Statistical Analysis及Table 1