基于全基因组测序的我国耐多药结核分枝杆菌耐药突变特征分析

目的 利用全基因组测序数据分析中国耐多药结核分枝杆菌的耐药相关基因突变谱及主要突变类型与菌株基因型的相关性。方法 查询并下载NCBI数据库中截至2019年8月公开发表的中国结核分枝杆菌全基因组测序数据，利用全基因组数据预测菌株分子药敏结果，统计不同药物耐药相关基因的突变类型，并分析耐药突变类型与菌株基因型的相关性。结果 根据分子药敏结果从2 019株菌株中鉴定出1 024株耐多药菌株，对常用抗结核药物的耐药相关基因主要突变类型分别为katG S315T（73.2%，异烟肼）、rpoB S450L（63.1%，利福平）、rpsL K43R（70.0%，链霉素）、embB M306V（37.4%，乙胺丁醇）、pncA启动子区T（-11）C（7.9%，吡嗪酰胺）、gyrA A90V（32.3%，氟喹诺酮类）、rrs A1401G（67.7%，二线注射类）、fabG1启动子区C（-15）T（87.0%，乙硫异烟胺）、folC I43T（30.4%，对氨基水杨酸）。其中，katG S315T、rpsL K43R、embB M306V、gyrA D94G在L2系菌株中的频率显著高于L4系菌株，folC I43T仅在L2系菌株中发现；katG S315T在古典北京型菌株中比例较高，而rpsL K43R在现代北京基因型菌株中的比例更高，其差异均有统计学意义（均P<0.05）。结论 本研究提供了基于全基因组测序的我国耐多药结核分枝杆菌对多种常用抗结核药的耐药相关基因主要突变类型，为研发敏感、特异的快速分子检测方法提供了依据；同时也发现多种耐药相关基因主要突变类型与菌株基因型有关