- Research article
- Open Access
- Open Peer Review
Prevalence and molecular characterization of pyrazinamide resistance among multidrug-resistant Mycobacterium tuberculosis isolates from Southern China
BMC Infectious Diseasesvolume 17, Article number: 711 (2017)
Pyrazinamide (PZA) plays a unique role in the treatment for multidrug-resistant tuberculosis (MDR-TB) in both first- and second-line regimens. The aim of this study was to investigate the prevalence and molecular characterization of PZA resistance among MDR-TB isolates collected in Chongqing municipality.
A total of 133 MDR-TB isolates were collected from the smear-positive tuberculosis patients who were registered at local TB dispensaries of Chongqing. PZA susceptibility testing was determined with a Bactec MGIT 960 system. In addition, the genes conferring for PZA resistance were screened by DNA sequencing.
Of these 133 MDR-TB isolates, 83 (62.4%) were determined as PZA-resistant by MGIT 960. In addition, streptomycin- (83.1% vs. 56.0%, P < 0.01), ofloxacin- (51.8% vs. 18.0%, P < 0.01), kanamycin- (22.9% vs. 2.0%, P < 0.01), amikacin- (18.1% vs. 2.0%, P = 0.01), capromycin-resistance (12.0% vs. 2.0%, P = 0.05), were more frequently observed among PZA-resistant isolates compared with PZA-susceptible isolates. Sequence analysis revealed that 73 out of 83 (88.0%) MDR strains harbored a mutation located in the pncA gene, including 55 (75.3%, 55/73) of single nucleotide substitutions and 18 (24.7%, 18/73) of frameshift mutation, while no genetic mutation associated with PZA resistance was found in the rpsA gene. The pncA expression of strains harboring substitution from A to G at position −11 in the promoter region of pncA was significantly lower than that of H37Rv (P < 0.01).
In conclusion, our data have demonstrated that the analysis of the pncA gene rather than rpsA gene provides rapid and accurate information regarding PZA susceptibility for MDR-TB isolates in Chongqing. In addition, loss of pncA expression caused by promoter mutation confers PZA resistance in MDR-TB isolates.
Multidrug resistant tuberculosis (MDR-TB) is a major concern hampering global tuberculosis control efforts [1, 2]. According to new estimates from World Health Organization (WHO), there were around 0.48 million new cases of MDR-TB cases, and approximately 0.19 million deaths from MDR-TB worldwide in 2014 (WHO, 2015). Only behind India, China has the second burden of MDR-TB globally, with 52, 000 prevalent MDR-TB cases annually (WHO, 2015). A recent national survey of drug-resistant tuberculosis in China revealed that 5.7% of new TB cases and 25.6% of previously treated cases were MDR-TB, respectively, which were higher than the global average rates . Given its high rate of treatment failure, the epidemic of MDR-TB constitutes a serious public health problem in China [3, 4].
Pyrazinamide (PZA) is one of cornerstone first-line anti-tuberculosis agents that is also commonly used as essential component in the therapeutic treatment of MDR-TB [5, 6]. As the prodrug, PZA requires conversion into its active form pyrazinoic acid (POA) by the enzyme pyrazinamidase (PZase). PZase is encoded by the 561-nucleotide pncA gene. Loss of PZase activity cause by genetic mutations in the pncA gene is the main mechanism of resistance to PZA in M. tuberculosis . Several recent literatures have reported that several PZA-resistant M. tuberculosis isolates harbored the mutations in the promoter of pncA gene, indicating these alternations may result in PZA resistance by influencing the expression of pncA [7,8,9]. In addition to pncA, another gene named rpsA, which encodes the 30S ribosomal protein S1, has been demonstrated to confer PZA resistance by structural alternation of POA binding . To date, the contribution of the rpsA mutations conferring PZA resistance is controversial, which requires more experimental evidences to elucidate the role of rpsA as a PZA resistance mechanism .
Chongqing is the largest municipality in the southwestern China . Due to the underdeveloped setting, Chongqing is considered as a hotspot for both TB and MDR-TB in China . However, limited data is available for the prevalence and molecular characteristics of PZA resistance in M. tuberculosis isolates, especially MDR-TB. In this study, our main goal was to investigate the prevalence of PZA resistance among MDR-TB isolates collected in Chongqing municipality. We also sought to analyze the mutant profiles of MDR-TB isolates conferring PZA resistance in this area.
All smear-positive tuberculosis patients who were registered at local TB dispensaries between November 2014 and February 2016 were enrolled in this study. Information was obtained from patient’s medical record. Two sputum samples were obtained from each smear-positive patient for culturing in the county-level laboratories. After 4–8 weeks of incubation, cultures with growing colonies on Löwenstein-Jensen (L-J) medium were sent to the Chongqing Tuberculosis Control Institute for further drug susceptibility testing. Re-treated cases were defined as patients having previously received more than one month of anti-TB treatment.
Drug susceptibility testing
Drug susceptibility testing (DST) was performed with the proportional method recommended by WHO [3, 13]. The concentrations of drugs in L-J medium were as follows: isoniazid (INH), 0.2 μg/ml; rifampicin (RIF), 40 μg/ml; ethambutol (EMB), 2 μg/ml; streptomycin (SM), 4 μg/ml; ofloxacin (OFLX), 2 μg/ml; kanamycin (KAN), 30 μg/ml; amikacin (AMK), 30 μg/ml; capromycin (CAP), 40 μg/ml; and protionamide (PTO), 40 μg/ml and p-aminosalicylic acid (PAS) 1 μg/ml . A strain was declared resistant to an antimicrobial agent when the growth rate exceeded 1% compared with the control. In addition, PZA susceptibility testing was determined with a Bactec MGIT 960 system according to the manufacturer’s instructions. The critical concentration of PZA in the liquid medium was 100 μg/ml . The MDR-TB strains were defined as those resistant to both isoniazid and rifampicin. In addition, XDR-TB was defined as MDR-TB resistant to any member of the quinolone family and at least one of the remaining second-line anti-TB injectable drugs .
DNA extraction and amplification
Genomic DNA was extracted from freshly cultured bacteria as previously described . Briefly, the harvested bacteria from the surface of L-J medium were suspended in 500 μl Tris-EDTA (TE) buffer and heated in a 95 °C water bath for 30 min. The crude DNA was used as template for amplification. The fragments of pncA and rpsA were amplified with the following primers: pncA-F 5′-AACAGTTCATCCCGGTTC-3′ and pncA-R 5′-GCGTCATGGACCCTATATC-3′; rpsA-F 5′-CGGAGCAACCCAACAATA-3′ and rpsA-R 5′-GTGGACAGCAACGACTTC-3′, respectively . The 50 μL PCR mixture was prepared as follows: 25 μl 2 × GoldStar MasterMix (CWBio, Beijing, China), 5 μL of DNA template, and 0.2 μM of each primer set. PCR parameters for amplification were 5 min at 94 °C followed by 35 cycles of 94 °C for 1 min, 58 °C for 1 min, 72 °C for 1 min and a final extension of 72 °C for 5 min. PCR products were sent to Tsingke company for sequencing. All sequence data were aligned with pncA and rpsA of reference strain H37Rv (ATCC) using BioEdit (version 188.8.131.52) software.
Quantitative reverse transcription PCR (QRT-PCR)
Bacteria were harvested from L-J medium after inoculation for 4 weeks. The total RNA was isolated following a standard Trizol RNA extraction protocol supplied by Invitrogen (Invitrogen, Life Technologies, USA) . Followed by treatment with DNaseI (Invitrogen, Life Technologies, USA), the reverse transcription was carried out using SuperScriptIII RT kit (Invitrogen, Life Technologies, USA). The relative expression level of pncA gene was detected by QRT-PCR in a 20 μl system containing 10 μL of 2 × UltraSYBR Mixture (CWBio, Beijing, China), 2 μL of cDNA template, and 0.2 μM of each primer set (pncA-QF 5′-GAAGCGGCGGACTACCATC-3′ and pncA-QR 5′-AGTGGCGTGCCGTTCTCG-3′). PolyA was set as the internal control in respective PCR experiments .
Chi square test was used to evaluate the associations among multiple categorical variables, and the statistical results were summarized with odds ratios (ORs) with 95% confidence intervals (CIs). All calculations were performed in SPSS 13.0 (SPSS Inc., USA). Differences with a P value less than 0.05 were declared statistically significant.
Demographic characteristics and drug susceptibility profiles
A total of 133 (10.8%) of 1236 clinical isolates were identified as MDR-TB, including 38 (28.6%) pre-XDR and 17 (12.8%) XDR. Overall, 80 (60.2%) strains were isolated from male patients and 53 (39.8%) from female patients. The average age of the patients was 46.6 years (range 19–76 years). In addition, 36.8% of isolates were from new cases, and 63.2% from re-treated cases. Out of the 133 MDR isolates tested, 97 (72.9%) were resistant to SM, 40 (30.1%) to EMB, 52 (39.1%) to OFLX, 20 (15.0%) to KAN, 16 (12.0%) to AMK, 11 (8.3%) to CAP and 15 (11.3%) to PAS (Table 1).
Factors associated with PZA resistance
As shown in Table 1, 83 (62.4%) isolates were determined as PZA-resistant by MGIT, while the other 50 (37.6%) were susceptible to PZA. In regard to the distribution of MDR-TB cases treatment history, the percentage of re-treated MDR-TB patients in the PZA-resistant group was significantly higher than in the PZA-susceptible group (odd ratio (OR) [95% confidence interval (CI): 3.26[1.55–6.82], P < 0.01). In contrast, our data revealed no statistically significant difference between the PZA-resistant and PZA-susceptible group in gender and age (P > 0.05).
We further analyzed the resistance profiles of other drugs among PZA-resistant and PZA-susceptible M. tuberculosis strains. Statistical analysis revealed that SM- (83.1% vs. 56.0%, P < 0.01), OFLX- (51.8% vs. 18.0%, P < 0.01), KAN- (22.9% vs. 2.0%, P < 0.01), AMK- (18.1% vs. 2.0%, P = 0.01), CAP-resistance (12.0% vs. 2.0%, P = 0.05), pre-XDR (36.1% vs. 16.0%, P < 0.02) and XDR (19.3% vs. 2.0%, P < 0.01) were more frequently observed among PZA-resistant isolates compared with PZA-susceptible isolates, while no difference was identified between EMB-resistant and EMB-susceptible strains (31.3% vs. 28.0%, P = 0.69) (Table 1).
Mutations in the pncA and rpsA gene
A total of 73 out of 83 (88.0%) MDR strains harbored a mutation located in the pncA gene, including 55 (75.3%, 55/73) of single nucleotide substitutions and 18 (24.7%, 18/73) of frameshift mutation. As summarized in Table 2, we observed great mutant diversity in pncA gene, and there were 48 different mutant types conferring PZA resistance among MDR strains in Chongqing. The most prevalent mutation associated with PZA resistance was found in codon 132 of pncA (12.3%, 9/73), resulting in the amino aicd substitution of Gly to Asp. As the second affected codon, five isolates had a mutation in codon 62 (6.8%, 5/73). Notably, the third frequent mutation was the nucleotide substitution from A to G at position −11 in the promoter region of pncA, accounting for 5.4% of mutant isolates. We also found that 4 PZA-susceptible isolates carried a genetic mutation in pncA, including 2 strains in codon 62, one in codon 67 and one in codon 154 (Table 3). In addition, no genetic mutation associated with PZA resistance was found in the rpsA gene in this study.
We analyzed the performance of pncA mutations for predicting PZA susceptibility. When setting the phenotypic PZA susceptibility as a gold standard, we found that detection of mutations in pncA gene exhibited a sensitivity of 88.0% (95% CI, 80.9%–95.0%) and a specificity of 92.0% (95% CI, 84.5%–99.5%) (Table 4).
Loss of pncA expression due to promoter mutation
We further explore whether the substitution at position −11 affected the expression level of pncA. As shown in Fig. 1, compared with reference strain H37Rv, the relative expression level of pncA in MDR87, MDR88, MDR114 and MDR126 carried the substitution at position −11 were 0.25 fold, 0.19 fold, 0.36 fold and 0.22 fold, respectively. Statistical analysis revealed that the pncA expression of strains harboring promoter mutation at position −11 was significantly lower than that of H37Rv (P < 0.01).
PZA plays a unique role in the treatment for MDR-TB in both first- and second-line regimens [19, 20]. The prevalence of PZA among MDR-TB thus is a determining factor for initiation of PZA in the therapy regimens for these refractory patients . Here, our data demonstrated that 62.4% of MDR-TB exhibited resistance against PZA in Chongqing, which was similar to a recent literature from Beijing (57.7%) , while higher than those from Zhejiang (43.1%) , Shanghai (38.5%) , United States (38.0%) , and Thailand (49.0%) . The high prevalence of PZA resistance among MDR-TB patients from our report indicates that Chongqing is a hotspot of PZA resistance in China. In our study, more than 60% MDR patients received previous anti-TB therapy with PZA, which is significantly higher than the average national level (21.8%) . Hence, we speculate that the high proportion of PZA resistance may be contributed to the high rate of re-treated TB patients. The serious issue on PZA resistance highlights the diminished role of PZA in the treatment for MDR-TB in this setting with high MDR-TB burden. Prior to the use of PZA for treatment of MDR-TB cases, it is essential to perform in vitro susceptibility testing against PZA to formulate a suitable regimen .
Another important finding from our observation was that we observed that there were high correlation between PZA resistance and several other drugs’ resistance, including OFLX, second-line injectable drugs, and PAS. Similar to our findings, a recent report from Alame-Emane and colleagues has revealed that PZA resistance in M. tuberculosis arises after RIF and fluoroquinolone (FQ) resistance . Genetic mutations constitute the most important mechanism conferring drug resistance in M. tuberculosis . Exposure to bacterial species to antimicrobial agents, including RIF, FQ and the aminoglycosides, induces the production of oxygen radicals, thereby conferring high frequency mutagenesis [25,26,27]. Considering long duration of anti-TB treatment, we hypothesize that MDR bacteria will harbor more genetic mutations induced by prolonged exposure to these drugs, which may be responsible for the potential cross resistance between PZA and other drugs in our study.
In vitro susceptibility against PZA is essential for proper management of MDR-TB with regimen containing PZA . However, phenotypic DST for PZA is not routinely performed due to the requirement of harshly acidic environment . Molecular method based on detecting the mutations in pncA and rpsA serves as an alternative to predict the PZA susceptibility in M. tuberculosis . In this study, our data demonstrated that genetic alternations in pncA confer 88.0% of PZA resistance among MDR-TB in Chongqing. A number of studies have demonstrated a diverse prevalence of pncA mutation among PZA resistant isolates in different regions, ranging from 45.7% in Brazil , 70.6% in Iran , 75.0% in Thailand , 78.0% in Zhejiang , 84.6% in Southern China , and 94.1% in Sweden . Hence, pncA mutations may differ from one geographic region to another. In addition, we found that pncA mutations exhibited great diversity, and the most frequent mutant type in codon 142 only accounted for approximate 12% of PZA resistant isolates, which was also different from reports from other regions [8, 9]. Given the diversity of pncA mutations within more than 500-bp long segment, DNA sequencing of the entire pncA is more effective for verification of PZA resistance rather than the routine methods by covering the mutant hotspots.
In addition, the third frequent mutation identified in this study was located at position −11 of the pncA promoter region. In line with our observation, numerous literatures have observed this mutant type in PZA-resistant M. tuberculosis isolates [31, 32]. We found that this substitution at position −11 was associated with low level of pncA expression, which was also consistent to the observation from Sheen et al. . The anti-TB activity of PZA depends on the transformation to POA by PZAse, which is encoded by pncA gene. The loss of pncA transcriptional level may result in the relative low PZAse activity, which is further associated with the phenotypic PZA resistance. Our results suggest that the promoter region of pncA is recommended to be included in the sequence analysis of pncA gene.
We acknowledge several limitations of this study. First, the small sample size is a major limitation of our report. And all the strains collected from one region also reduce their representativeness and thereby constrict the generalizability of findings. Further wider sampling will give more credence to this study. Second, although the primary data of this study suggest that loss of pncA expression caused by promoter mutation confers PZA resistance in MDR-TB isolates, we could give no biochemical or transgenic substantiation of this statement. Therefore, there is an urgent need to confirm our findings with more experimental evidences in the future. Third, the high diversity of pncA mutations in MTB-TB isolates from Chongqing underscores previous findings that there is no clear hotspot for pncA mutations , and several novel mutations in pncA gene were found for first time among PZA-resistant isolates. Despite being highly correlated with the loss of PZA susceptibility, further experiments will be carried out to clarify the potential contributions of these mutations to PZA resistance. Nevertheless, our report firstly described the molecular characteristics of PZA resistance among MDR-TB isolates from Southern China, which provides important hints to diagnose PZA resistance and help guide therapy with PZA for MDR-TB patients in this region with high MDR-TB burden.
In conclusion, our data have demonstrated that the analysis of the pncA gene rather than rpsA gene provides rapid and accurate information regarding PZA susceptibility for MDR-TB isolates in Chongqing. In addition, loss of pncA expression caused by promoter mutation confers PZA resistance in MDR-TB isolates. Considering the high degree of diversity of pncA gene mutations, DNA sequencing of the entire pncA is more effective for verification of PZA resistance rather than the routine methods by covering the mutant hotspots. The high prevalence of PZA resistance among MDR highlights the diminished role of PZA in the treatment for MDR-TB in this setting with high TB burden.
Drug susceptibility testing
Gandhi NR, Nunn P, Dheda K, Schaaf HS, Zignol M, van Soolingen D, Jensen P, Bayona J. Multidrug-resistant and extensively drug-resistant tuberculosis: a threat to global control of tuberculosis. Lancet. 2010;375(9728):1830–43.
Zhang Z, Pang Y, Wang Y, Liu C, Zhao Y. Beijing genotype of Mycobacterium tuberculosis is significantly associated with linezolid resistance in multidrug-resistant and extensively drug-resistant tuberculosis in China. Int J Antimicrob Agents. 2014;43(3):231–5.
Zhao Y, Xu S, Wang L, Chin DP, Wang S, Jiang G, Xia H, Zhou Y, Li Q, Ou X, et al. National survey of drug-resistant tuberculosis in China. N Engl J Med. 2012;366(23):2161–70.
Orenstein EW, Basu S, Shah NS, Andrews JR, Friedland GH, Moll AP, Gandhi NR, Galvani AP. Treatment outcomes among patients with multidrug-resistant tuberculosis: systematic review and meta-analysis. Lancet Infect Dis. 2009;9(3):153–61.
Verdugo D, Fallows D, Ahuja S, Schluger N, Kreiswirth B, Mathema B. Epidemiologic Correlates of Pyrazinamide-Resistant Mycobacterium tuberculosis in New York City. Antimicrob Agents Chemother. 2015;59(10):6140–50.
Zhang Y, Mitchison D. The curious characteristics of pyrazinamide: a review. Int J Tuberc Lung Dis. 2003;7(1):6–21.
Scorpio A, Zhang Y. Mutations in pncA, a gene encoding pyrazinamidase/nicotinamidase, cause resistance to the antituberculous drug pyrazinamide in tubercle bacillus. Nat Med. 1996;2(6):662–7.
Xia Q, Zhao LL, Li F, Fan YM, Chen YY, Wu BB, Liu ZW, Pan AZ, Zhu M. Phenotypic and genotypic characterization of pyrazinamide resistance among multidrug-resistant Mycobacterium tuberculosis isolates in Zhejiang. China Antimicrob Agents Chemother. 2015;59(3):1690–5.
Tan Y, Hu Z, Zhang T, Cai X, Kuang H, Liu Y, Chen J, Yang F, Zhang K, Tan S, et al. Role of pncA and rpsA gene sequencing in detection of pyrazinamide resistance in Mycobacterium tuberculosis isolates from southern China. J Clin Microbiol. 2014;52(1):291–7.
Shi W, Zhang X, Jiang X, Yuan H, Lee JS, Barry CE 3rd, Wang H, Zhang W, Zhang Y. Pyrazinamide inhibits trans-translation in Mycobacterium tuberculosis. Science. 2011;333(6049):1630–2.
Zhang D, An J, Wang J, Hu C, Wang Z, Zhang R, Wang Y, Pang Y. Molecular typing and drug susceptibility of Mycobacterium tuberculosis isolates from Chongqing Municipality, China. Infect Genet Evol. 2013;13:310–6.
Zhang D, An J, Wang Y, Pang Y. Genetic diversity of multidrug-resistant tuberculosis in a resource-limited region of China. Int J Infect Dis. 2014;29:7–11.
World Health Organization: Global tuberculosis report 2015. Geneva: World Health Organization; 2015.
World Health Organization. Anti-tuberculosis drug resistance in the world. Geneva: World Health Organization; 2008.
Pang Y, Zhou Y, Zhao B, Liu G, Jiang G, Xia H, Song Y, Shang Y, Wang S, Zhao YL. Spoligotyping and drug resistance analysis of Mycobacterium tuberculosis strains from national survey in China. PLoS One. 2012;7(3):e32976.
Zhang Z, Lu J, Wang Y, Pang Y, Zhao Y. Prevalence and molecular characterization of fluoroquinolone-resistant Mycobacterium tuberculosis isolates in China. Antimicrob Agents Chemother. 2014;58(1):364–9.
Pang Y, Wang Z, Zheng H, Song Y, Wang Y, Zhao Y. Pyrazinamide resistance determined by liquid culture at low pH better correlates with genetic mutations in MDR tuberculosis isolates. J Microbiol Methods. 2015;119:142–4.
Pang Y, Lu J, Wang Y, Song Y, Wang S, Zhao Y. Study of the rifampin monoresistance mechanism in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2013;57(2):893–900.
Kurbatova EV, Cavanaugh JS, Dalton T, Click ES, Cegielski JP. Epidemiology of pyrazinamide-resistant tuberculosis in the United States, 1999-2009. Clin Infect Dis. 2013;57(8):1081–93.
Zhang Y, Yew WW. Mechanisms of drug resistance in Mycobacterium tuberculosis: update 2015. Int J Tuberc Lung Dis. 2015;19(11):1276–89.
Zhang Y, Chiu Chang K, Leung CC, Wai Yew W, Gicquel B, Fallows D, Kaplan G, Chaisson RE, Zhang W. ‘Z(S)-MDR-TB’ versus ‘Z(R)-MDR-TB’: improving treatment of MDR-TB by identifying pyrazinamide susceptibility. Emerg Microbes Infect. 2012;1(7):e5.
Gu Y, Yu X, Jiang G, Wang X, Ma Y, Li Y, Huang H. Pyrazinamide resistance among multidrug-resistant tuberculosis clinical isolates in a national referral center of China and its correlations with pncA, rpsA, and panD gene mutations. Diagn Microbiol Infect Dis. 2016;84(3):207–11.
Xu P, Wu J, Yang C, Luo T, Shen X, Zhang Y, Nsofor CA, Zhu G, Gicquel B, Gao Q. Prevalence and transmission of pyrazinamide resistant Mycobacterium tuberculosis in China. Tuberculosis (Edinb). 2016;98:56–61.
Jonmalung J, Prammananan T, Leechawengwongs M, Chaiprasert A. Surveillance of pyrazinamide susceptibility among multidrug-resistant Mycobacterium tuberculosis isolates from Siriraj Hospital. Thailand BMC Microbiol. 2010;10:223.
Alame-Emane AK, Xu P, Pierre-Audigier C, Cadet-Daniel V, Shen X, Sraouia M, Siawaya JF, Takiff H, Gao Q, Gicquel B. Pyrazinamide resistance in Mycobacterium tuberculosis arises after rifampicin and fluoroquinolone resistance. Int J Tuberc Lung Dis. 2015;19(6):679–84.
Ysern P, Clerch B, Castano M, Gibert I, Barbe J, Llagostera M. Induction of SOS genes in Escherichia coli and mutagenesis in Salmonella typhimurium by fluoroquinolones. Mutagenesis. 1990;5(1):63–6.
Baharoglu Z, Mazel D. Vibrio cholerae triggers SOS and mutagenesis in response to a wide range of antibiotics: a route towards multiresistance. Antimicrob Agents Chemother. 2011;55(5):2438–41.
Bhuju S, Fonseca Lde S, Marsico AG, de Oliveira Vieira GB, Sobral LF, Stehr M, Singh M, Saad MH. Mycobacterium tuberculosis isolates from Rio de Janeiro reveal unusually low correlation between pyrazinamide resistance and mutations in the pncA gene. Infect Genet Evol. 2013;19:1–6.
Doustdar F, Khosravi AD, Farnia P. Mycobacterium tuberculosis genotypic diversity in pyrazinamide-resistant isolates of Iran. Microb Drug Resist. 2009;15(4):251–6.
Jureen P, Werngren J, Toro JC, Hoffner S. Pyrazinamide resistance and pncA gene mutations in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2008;52(5):1852–4.
Sheen P, Lozano K, Gilman RH, Valencia HJ, Loli S, Fuentes P, Grandjean L, Zimic M. pncA gene expression and prediction factors on pyrazinamide resistance in Mycobacterium tuberculosis. Tuberculosis (Edinb). 2013;93(5):515–22.
Portugal I, Barreiro L, Moniz-Pereira J, Brum L. pncA mutations in pyrazinamide-resistant Mycobacterium tuberculosis isolates in Portugal. Antimicrob Agents Chemother. 2004;48(7):2736–8.
We would like to thank the staffs from the National Tuberculosis Reference Laboratory for their technical support.
This study was supported by the Major State Basic Research Development Program of China (2014CB744403).
Availability of data and materials
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Chongqing Tuberculosis Control Institute. Patients enrolled in this study provided the signed Informed Consent form.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.