基于新疆维吾尔自治区南疆地区169株结核分枝杆菌探讨基因缺失对贝达喹啉的耐药机制

目的挖掘与结核分枝杆菌(MTB)对贝达喹啉耐药相关的基因缺失片段。方法2017—2020年期间从新疆维吾尔自治区南疆地区收集807株MTB分离株,对菌株进行全基因组测序,采用TB-Profiler 4.1.2检测菌株的家系和预测异烟肼、利福平、氟喹诺酮类等7种药物的敏感性;选择44株准广泛耐药(Pre-XDR)菌株和125株随机菌株,采用Delly 0.9.1软件检测基因缺失情况,采用C++脚本程序检测Rv0678、atpE和pepQ基因突变情况,采用微量肉汤稀释微孔板法检测菌株对贝达喹啉的最小抑菌浓度(MIC)及耐药情况。采用χ^(2)检验和秩和检验等方法分析基因缺失与贝达喹啉耐药的相关性及不同家系菌株对贝达喹啉耐药性及敏感程度的差异。结果44株Pre-XDR株中共鉴定出23株贝达喹啉耐药株(MIC≥0.25μg/mL)和49种基因缺失类型,发生缺失菌株数>5株的缺失类型共有8种,其中2729618 bp~2729833 bp(Rv2434c)缺失组贝达喹啉耐药率和MIC水平低于未缺失组(χ^(2)=6.145,P=0.013和Z=−2.224,P=0.026);2235169 bp~2235233 bp(ctpG)、1710766 bp~1711557 bp(Rv1519-Rv1520)和3952924 bp~3954223 bp(echA19-Rv3517)多基因缺失组贝达喹啉耐药率和MIC水平均高于未缺失组(χ^(2)=4.322,P=0.038和Z=−2.526,P=0.012)。但是在125株验证用菌株(含贝达喹啉耐药株23株)中均未发现差异具有统计学意义(χ^(2)=1.660,P=0.198和Z=−1.152,P=0.249)。将44株菌株和125株菌株合并分析,仅发现2235169 bp~2235233 bp(ctpG)、1710766 bp~1711557 bp(Rv1519-Rv1520)和3952924 bp~3954223 bp(echA19-Rv3517)多基因缺失组贝达喹啉耐药率高于未缺失组(χ^(2)=6.396,P=0.011),该多基因缺失仅发生于L3家系菌株,同时两两比较发现L3家系菌株的贝达喹啉耐药率显著高于L4家系(P<0.017)。结论2235169 bp~2235233 bp(ctpG)、1710766 bp~1711557 bp(Rv1519-Rv1520)和3952924 bp~3954223 bp(echA19-Rv3517)多基因缺失为南疆地区L3家系菌株特有的基因缺失类型(与该地其他家系相比),其与贝达喹啉耐药性的关系不稳健,L3家系菌株可能比L4家系对贝达喹啉耐药的易感性更高。由于样本量较小,需扩大样本量做进一步验证,以上结果为了解贝达喹啉的耐药机制及控制相关耐药菌株的传播提供了线索。