Whole genome sequencing of drug resistant Mycobacterium tuberculosis isolates from a high burden tuberculosis region of North West Pakistan

Scientific Reports - Tập 9 Số 1
Abdul Jabbar1, Jody Phelan2, Paola Flórez de Sessions3, Taj Ali Khan4, Hazır Rahman5, Sadiq Noor Khan1, Daire Cantillon6, Letícia Muraro Wildner6, Sajid Ali7, Susana Campino2, Simon J. Waddell6, Taane G. Clark8
1Department of Medical Lab Technology, University of Haripur, Haripur, Pakistan
2Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London WC1E 7HT, UK
3Genome Institute of Singapore, 60 Biopolis St, Biopolis, Singapore
4Department of Microbiology, Kohat University of Science and Technology, Kohat, Pakistan
5Department of Microbiology, Abdul Wali Khan University, Mardan, Pakistan
6Department of Global Health and Infection, Brighton and Sussex Medical School, University of Sussex, Brighton, BN1 9PX, UK
7Provincial Tuberculosis Reference Laboratory, Hayatabad Medical Complex Peshawar, Khyber Pakhtunkhwa, Pakistan
8Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK

Tóm tắt

Abstract

Tuberculosis (TB), caused by Mycobacterium tuberculosis bacteria, is a leading infectious cause of mortality worldwide, including in Pakistan. Drug resistant M. tuberculosis is an emerging threat for TB control, making it important to detect the underlying genetic mutations, and thereby inform treatment decision making and prevent transmission. Whole genome sequencing has emerged as the new diagnostic to reliably predict drug resistance within a clinically relevant time frame, and its deployment will have the greatest impact on TB control in highly endemic regions. To evaluate the mutations leading to drug resistance and to assess for evidence of the transmission of resistant strains, 81 M. tuberculosis samples from Khyber Pakhtunkhwa province (North West Pakistan) were subjected to whole genome sequencing and standard drug susceptibility testing for eleven anti-TB drugs. We found the majority of M. tuberculosis isolates were the CAS/Delhi strain-type (lineage 3; n = 57; 70.4%) and multi-drug resistant (MDR; n = 62; 76.5%). The most frequent resistance mutations were observed in the katG and rpoB genes, conferring resistance to isoniazid and rifampicin respectively. Mutations were also observed in genes conferring resistance to other first and second-line drugs, including in pncA (pyrazinamide), embB (ethambutol), gyrA (fluoroquinolones), rrs (aminoglycosides), rpsL, rrs and giB (streptomycin) loci. Whilst the majority of mutations have been reported in global datasets, we describe unreported putative resistance markers in katG, ethA (ethionamide), gyrA and gyrB (fluoroquinolones), and pncA. Analysis of the mutations revealed that acquisition of rifampicin resistance often preceded isoniazid in our isolates. We also observed a high proportion (17.6%) of pre-MDR isolates with fluoroquinolone resistance markers, potentially due to unregulated anti-TB drug use. Our isolates were compared to previously sequenced strains from Pakistan in a combined phylogenetic tree analysis. The presence of lineage 2 was only observed in our isolates. Using a cut-off of less than ten genome-wide mutation differences between isolates, a transmission analysis revealed 18 M. tuberculosis isolates clustering within eight networks, thereby providing evidence of drug-resistant TB transmission in the Khyber Pakhtunkhwa province. Overall, we have demonstrated that drug-resistant TB isolates are circulating and transmitted in North West Pakistan. Further, we have shown the usefulness of whole genome sequencing as a diagnostic tool for characterizing M. tuberculosis isolates, which will assist future epidemiological studies and disease control activities in Pakistan.

Từ khóa


Tài liệu tham khảo

Global tuberculosis report 2018 (2018).

Tahseen, S. et al. Use of Xpert® MTB/RIF assay in the first national anti-tuberculosis drug resistance survey in Pakistan. Int. J. Tuberc. Lung Dis. 20, 448–455 (2016).

Onozaki, I. et al. National tuberculosis prevalence surveys in Asia, 1990–2012: an overview of results and lessons learned. Trop. Med. Int. Heal. 20, 1128–1145 (2015).

Falzon, D. et al. World Health Organization treatment guidelines for drug-resistant tuberculosis, 2016 update. Eur. Respir. J. 49, 1602308 (2017).

Mishra, G. Current updates in tuberculosis. npj Prim. Care Respir. Med. 27, 1 (2017).

Ayaz, A. et al. Characterizing Mycobacterium tuberculosis isolates from Karachi, Pakistan: drug resistance and genotypes. Int. J. Infect. Dis. 16, (e303–e309 (2012).

Hasan, R. et al. Extensively Drug-Resistant Tuberculosis, Pakistan. Emerg. Infect. Dis. 16, 1473 (2010).

Gandhi, N. R. et al. Multidrug-resistant and extensively drug-resistant tuberculosis: a threat to global control of tuberculosis. Lancet 375, 1830–1843 (2010).

Coll, F. et al. Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis. Nat. Genet. 50, 307–316 (2018).

Coll, F. et al. Rapid determination of anti-tuberculosis drug resistance from whole-genome sequences. Genome Med. 7, 51 (2015).

Perdigão, J. et al. Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting. BMC Genomics 15, 991 (2014).

Kanji, A. et al. Single nucleotide polymorphisms in efflux pumps genes in extensively drug resistant Mycobacterium tuberculosis isolates from Pakistan. Tuberculosis 107, 20–30 (2017).

Witney, A. A. et al. Clinical application of whole-genome sequencing to inform treatment for multidrug-resistant tuberculosis cases. J. Clin. Microbiol. 53, 1473–83 (2015).

Phelan, J. E. et al. Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs. Genome Med. 11, 41 (2019).

Coll, F. et al. A robust SNP barcode for typing Mycobacterium tuberculosis complex strains. Nat. Commun. 5, 4812 (2014).

Gagneux, S. Host-pathogen coevolution in human tuberculosis. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 367, 850–9 (2012).

Benavente, E. D. et al. PhyTB: Phylogenetic tree visualisation and sample positioning for M. tuberculosis. BMC Bioinformatics 16, 155 (2015).

Guerra-Assunção, J. et al. Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area. Elife 4 (2015).

Ali, A. et al. Whole Genome Sequencing Based Characterization of Extensively Drug-Resistant Mycobacterium tuberculosis Isolates from Pakistan. PLoS One 10, e0117771 (2015).

Tanveer, M. et al. Genotyping and drug resistance patterns of M. tuberculosis strains in Pakistan. BMC Infect. Dis. 8, 171 (2008).

Watanabe Pinhata, J. M. et al. Use of an immunochromatographic assay for rapid identification of Mycobacterium tuberculosis complex clinical isolates in routine diagnosis. J. Med. Microbiol. 67, 683–686 (2018).

Belisle, J. T., Mahaffey, S. B. & Hill, P. J. Isolation of Mycobacterium Species Genomic DNA. In Methods in molecular biology (Clifton, N.J.) 465, 1–12 (2009).

Woods, G. L. et al. Susceptibility testing of mycobacteria, nocardiae, and other aerobic actinomycetes. M24-A2 31 ( 5 ) (2011).

Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. (2013).

Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993 (2011).

Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).

Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256–W259 (2019).

Coll, F. et al. SpolPred: rapid and accurate prediction of Mycobacterium tuberculosis spoligotypes from short genomic sequences. Bioinformatics 28, 2991–2993 (2012).

Plinke, C., Walter, K., Aly, S., Ehlers, S. & Niemann, S. Mycobacterium tuberculosis embB codon 306 mutations confer moderately increased resistance to ethambutol in vitro and in vivo. Antimicrob. Agents Chemother. 55, 2891–6 (2011).

Malinga, L., Brand, J., Olorunju, S., Stoltz, A. & van der Walt, M. Molecular analysis of genetic mutations among cross-resistant second-line injectable drugs reveals a new resistant mutation in Mycobacterium tuberculosis. Diagn. Microbiol. Infect. Dis. 85, 433–437 (2016).

Falzon, D. et al. Resistance to fluoroquinolones and second-line injectable drugs: impact on multidrug-resistant TB outcomes. Eur. Respir. J. 42, 156–68 (2013).

Hasan, Z. et al. Spoligotyping of Mycobacterium tuberculosis isolates from Pakistan reveals predominance of Central Asian Strain 1 and Beijing isolates. J. Clin. Microbiol. 44, 1763–8 (2006).

Singh, U. B. et al. Predominant tuberculosis spoligotypes, Delhi, India. Emerg. Infect. Dis. 10, 1138–42 (2004).

Manson, A. L. et al. Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance. Nat. Genet. https://doi.org/10.1038/ng.3767 (2017).

Thông tin
Thông tin xuất bản

Scientific Reports

Tập 9 Số 1

Thông tin tác giả