结核分枝杆菌耐药相关单碱基突变的计算方法比较  

Evalutation of three methods searching resistant-related mutations in mycobacterium tuberculosis genome

作　　者：李志远 刘飞[2] 曹德民 潘元龙[2] 律娜[2] 刘东鑫 贺文从 程城 王亚楠 赵雁林[3] 朱宝利 Zhiyuan Li;Fei Liu;Demin Cao;Yuanlong Pan;Na Lü;Dongxin Liu;Wencong He;Cheng Cheng;Yanan Wang;Yanlin Zhao;Baoli Zhu(University of Chinese Academy of Sciences,Beijing 100049,China;CAS Key Laboratory of Pathogenic Microbiology and Immunology,Institute of Microbiology,Chinese Academy of Sciences,Beijing 100101,China;National Tuberculosis Reference Laboratory,National Center for Tuberculosis Control and Prevention,Chinese Center for Disease Control and Prevention,Beijing 102206,China)

机构地区：[1]中国科学院大学,北京100049 [2]中国科学院微生物研究所,病原微生物与免疫学重点实验室,北京100101 [3]中国疾病预防控制中心结核病预防控制中心,国家结核病参比实验室,北京102206

出　　处：《微生物学报》2021年第11期3653-3666,共14页Acta Microbiologica Sinica

基　　金：国家重点研发计划(2018YFC1603803);中国科学院战略性先导科技专项(XDB29020203);北京市科技计划(Z201100005520041)。

摘　　要：【目的】耐药结核分枝杆菌(drug-resistant Mycobacterium tuberculosis)的产生给结核病(tuberculosis)的治疗带来巨大困难。【方法】使用基于全基因组测序的关联分析探究耐药强相关的单核苷酸多态性(single nucleotide polymorphism,SNP)突变,主要有GEMMA、phyc、plink。为了阐明其中最优的耐药相关SNP计算方法,本研究下载NCBI上已有的1504株结核分枝杆菌数据,并获取它们对于3种常见的一线抗结核治疗药物(isoniazid、rifampicin、ethambutol)的耐药性检验结果。并使用这3种耐药相关SNP计算方法计算与结核分枝杆菌耐药相关的SNP;并评估计算得到的耐药相关SNP在预测耐药表型的敏感性和特异性。【结果】发现通过phyc可以预测到最多的已知耐药相关SNP和最少的耐药无关SNP,而且phyc预测的耐药相关SNP的敏感性和特异性恒定大于52.49%。【结论】phyc在预测结核分枝杆菌耐药相关SNP中结果最准确,但考虑到运行时间和表型数据的更新,GEMMA和plink的结果也应作为参考。[Objective]The emergence of durg-resistant Mycobacterium tuberculosis brought a tough challenge to tuberculosis treatment.[Methods]Considering the lack of homologous recombination in Mycobacterium tuberculosis’s genome,the mutation strongly related to resistance could be efficiently confirmed by genome-wide association analysis.However,many resistance-related mutations had yet been found by existing methods.To calculate resistant-related mutations,researchers commonly used some methods similar to genome-wide association analysis(GWAS),which mainly included genome-wide efficient mixed model association(GEMMA),phylogenetic convergence test(phyc),plink.To find the better method among the three methods when calculating resistant-related SNPs on non-mobile antibiotic resistance gene,the genomes of 1504 M.tuberculosis strains from Hunan province and National Center of Biotechnology Information was obtained with their phenotypes of three first-line antibiotics(isoniazid,rifampicin,ethambutol).The three methods were performed to calculate the association between phenotype and known or novel SNPs related to resistance,and their sensibility and specificity were evaluated by the resistant-related SNPs got by the three methods.[Results]Phyc was able to search more known resistance-related SNPs with the sensibility and specificity higher than 52.49%.[Conclusion]Phyc is the most accurate in predicting the SNP related to drug resistance of Mycobacterium tuberculosis,but considering the update of running time and phenotypic data,the results of GEMMA and plink should also be used as a reference.

关 键 词：结核分枝杆菌 耐药基因 比较基因组学 

分 类 号：R378.911[医药卫生—病原生物学]