Background
Tuberculosis (TB) is a global public health issue, getting worse due to emergence of resistance. Pyrazinamide (PZA) is first-line antimicrobial drugs used against non-replicated Mycobacterium tuberculosis (MTB). Data is scarce about whole genome sequencing of PZA resistance (PZA-R) in Khyber Pakhtunkhwa (KP) province of high burden country, Pakistan. In the current study we aimed to find the most common mutations in PZA-R MTB isolates in association with other candidate genes in a whole genome sequence (WGS). Samples were collected from TB suspects and drug susceptibility testing (DST) was performed according according to the WHO standards. The resistant samples were subjected for whole genome sequencing (WGS). The sequence data was through MTBseq and Total Genotyping Solution for Mycobacterium tuberculosis (TGS-TB). Metabolic model was analyzed, using RAST server.
Results
Among the three whole genome sequences, (NCBI BioProject Accession: PRJNA629298, PRJNA629388) 1997, 1162, and 2053 mutations including indel, was detected. Diverse variability has been detected in the membrane proteins PE and PPE, modulating the host immune response. Nine mutations in coding and promotor region have been detected in pncA with one novel (T-4C) variants. Mutations in the other drug candidate genes, KatG, rpoB have also been detected.
Conclusion
The metabolic model shows a distinct property. Diversity of variants has been detected in majority of MTB essential genes, functions from cell growth to cell signaling. The current study provides useful information, associated with geographic specific strains for biomarkers development and better management of drug resistance isolates.