Skip to main content

Pyrazinamide resistance-conferring mutations in pncA and the transmission of multidrug resistant TB in Georgia

Abstract

Background

The ongoing epidemic of multidrug-resistant tuberculosis (MDR-TB) in Georgia highlights the need for more effective control strategies. A new regimen to treat MDR-TB that includes pyrazinamide (PZA) is currently being evaluated and PZA resistance status will largely influence the success of current and future treatment strategies. PZA susceptibility testing was not routinely performed at the National Reference Laboratory (NRL) in Tbilisi between 2010 and September 2015. We here provide a first insight into the prevalence of PZA resistant TB in this region.

Methods

Phenotypic susceptibility to PZA was determined in a convenience collection of well-characterised TB patient isolates collected at the NRL in Tbilisi between 2012 and 2013. In addition, the pncA gene was sequenced and whole genome sequencing was performed on two isolates.

Results

Out of 57 isolates tested 33 (57.9%) showed phenotypic drug resistance to PZA and had a single pncA mutation. All of these 33 isolates were MDR-TB strains. pncA mutations were absent in all but one of the 24 PZA susceptible isolate. In total we found 18 polymorphisms in the pncA gene. From the two major MDR-TB clusters represented (94–32 and 100–32), 10 of 15, 67.0% and 13 of 14, 93.0% strains, respectively were PZA resistant. We also identified a member of the potentially highly transmissive clade A strain carrying the characteristic I6L substitution in PncA. Another strain with the same MLVA type as the clade A strain acquired a different mutation in pncA and was genetically more distantly related suggesting that different branches of this particular lineage have been introduced into this region.

Conclusion

In this high MDR-TB setting more than half of the tested MDR-TB isolates were resistant to PZA. As PZA is part of current and planned MDR-TB treatment regimens this is alarming and deserves the attention of health authorities. Based on our typing and sequence analysis results we conclude that PZA resistance is the result of primary transmission as well as acquisition within the patient and recommend prospective genotyping and PZA resistance testing in high MDR-TB settings. This is of utmost importance in order to preserve bacterial susceptibility to PZA to help protect (new) second line drugs in PZA containing regimens.

Peer Review reports

Background

Pyrazinamide (PZA) has been part of the standard first line TB treatment since the 1980s. Adding PZA to the treatment regimen had a significant effect on the duration of treatment needed to cure TB (six instead of nine months) (reviewed in [1, 2]). Its unique activity against metabolically inactive, slow growing bacteria make PZA not only an invaluable component of the current first line treatment regimen, but also attractive for inclusion in novel treatment regimens for susceptible TB [3, 4] as well as for the treatment of drug resistant TB [5]. In high MDR-TB settings worse treatment outcomes are observed if the infecting strain is resistant to PZA, especially when resistance to fluoroquinolones is also present [6, 7]. It is therefore of utmost importance to preserve bacterial susceptibility to PZA in order to ensure the effective treatment of (MDR-)TB.

Drug susceptibility testing (DST) for PZA is not routinely used to guide (first and second line) anti-TB therapy. Due to the technical complexity of DST for PZA the results of phenotypic tests assessing the susceptibility of TB bacteria to PZA are not always considered to accurately reflect the resistance status of the bacteria [8]. The MGIT system is the most standardised and reliable method to perform DST for PZA [1, 9] and is currently the only commercially available phenotypic test to investigate PZA susceptibility. Here, clinical specimens are cultured in the presence of PZA under acidic conditions mimicking the in vivo conditions under which pyrazinoic acid (the activated drug) is active [1]. These issues have stimulated efforts to inform the design of molecular assays for the detection of PZA resistance [10, 11] but these methods are not currently widely applied.

For the reasons explained above the WHO recommends inclusion of PZA in the treatment regimen regardless of the results of PZA susceptibility testing [12]. Consequently, PZA susceptibility testing is often not prioritised. As a result there is often only limited data available on PZA resistance trends and epidemiology.

Resistance to PZA is largely conferred by mutations in the pncA gene, which encodes the pyrazinamidase enzyme that activates the drug [10, 13]. Studies that have assessed PZA resistance typically report about 40–50% in primary MDR-TB [14,15,16,17] and more than 90% in XDR-TB isolates (TBNET) [18]. Alignment of the protein sequence of PZA from various bacterial species revealed three conserved regions [19], moreover three mutational hot spot regions in pncA have been proposed [20]. However recent systematic reviews comparing large data sets do not support extensive clustering of mutations in pncA [10, 11] but rather define large numbers of high confidence mutations that when combined can identify molecular PZA resistance with an accuracy between 89.5 and 98.9% [10]. The observed diversity of pncA mutations found in clinical PZA resistant strains is in strong contrast to the limited spectrum of mutations conferring resistance to isoniazid and rifampicin [21, 22]. Clusters of PZA resistant strains with identical pncA mutations have been reported but so far have been geographically restricted [23,24,25,26,27]. The presence of clustered strains all having the same rifampicin- and isoniazid-resistance conferring mutations but dispersed pncA mutations has stimulated the idea that PZA resistance mutations are acquired post-MDR and that the presence of PZA conferring mutations could represent a bottleneck for transmission [28].

The National Tuberculosis Reference centre in Tbilisi, Georgia, suspended DST for PZA in 2010 (M. Akhalaia/ N. Tukvadze, personal communication) but have resumed DST for PZA in September 2015 (N. Bablishvili, personal communication). PZA remains an important component of the standard first line and second line treatment regimen in this setting. Georgia is recognised as one of the high-burden MDR-TB countries with 16.6% MDR-TB patients reported among all laboratory confirmed pulmonary TB patients in 2014 [29]. Although the percentage of new MDR-TB cases remains constant the percentage of retreatment MDR-TB cases is increasing. The incidence and prevalence of PZA resistance is unknown in this setting. Here we report phenotypic susceptibility and molecular resistance to PZA in a convenience collection of primary MDR-TB isolates from Georgia collected between 2012 and 2013.

Material and methods

Patient samples

Samples were part of a larger study of 399 clinical isolates of newly diagnosed patients with pulmonary TB collected between 2012 and 2013 at the NCTLD in Tbilisi, Georgia [30]. Of the 399 samples 67 were subjected to more detailed genetic investigation. These samples were selected to reflect the high proportion of MDR-TB strains and strains of the Beijing genotype previously found in this collection. Coletsos slopes of all 67 samples were subjected to PZA DST and Sanger sequencing of the pncA gene. Interpretable PZA DST results were obtained for 57 strains. Two strains were subjected to WGS.

The study was exempt from the approval of the Local Ethics Committee of the National Center for Tuberculosis and Lung Diseases (50 Maruashvili St, 0101 Tbilisi, Georgia) and Informed Consent was not required as the patient information used was anonymised before linking to the results of the analysis of the bacterial cultures and could not be linked back to individual patients.

Drug susceptibility testing

Drug susceptibility testing for PZA was performed according to the manufacturer’s instructions in the BD BACTEC MGIT 960 system. Bacteria from an inoculum of at least one day positive and a maximum of 5 day positive, according to the BACTEC MGIT, were exposed to 100 ul/ml PZA in an pH 5.9 environment. A strain was determined resistant if at the moment when the growth control (GC) reaches 400 growth units (GU) the PZA tube had reached at least a 100 growth units.

Sanger sequencing

Of the 57 samples tested by PZA DST 33 samples that were PZA resistant by DST and all 24 strains that were susceptible had their pncA gene amplified and sequenced. The pncA gene was sequenced using pncA_1F: GGC CGC GAT GAC ACC TCT, pncA_1R: GCC GCA GCC AAT TCA GCA GT, pncA_2F: CGA AGC GGC GGA CTA CCA TCA CG and pncA_2R: CCC CAC CTG CGG CTG CGA ACC and the Sanger sequencing technology (BaseClear, Leiden, The Netherlands). Sequences were analysed using BioNumerics (Applied Maths).

Whole genome sequencing

CTAB purified DNA was prepared from isolates 12-17889 and 12-15893 and whole genome sequences were prepared at GATC Biotech (GATC, Konstanz, Germany) on a Illumina HiSeq 2500 device using paired-end reads of 2x150bp and minimum coverage of 400 reads in the core genome. The FASTQ files were mapped to the H37Rv reference genome using bwa-mem software. From the resulting alignments, SAMtools and GATK software suites (default parameter settings) were used to call SNPs and small indels, and the intersection of variants between the methods retained. A full description of the pipeline is described elsewhere [31, 32]). The strain 12–17,889 was shared with N. Casali (Devision of Infectious Diseases, Imperial College London). Sequence analysis as well as phylogenetic analysis of strain 12-17889 was performed as previously described [26].

MIRU-VNTR typing

All 57 samples were typed by 24-locus VNTR [33] either at the RIVM or by Genoscreen (Lille, France). MLVA MtbC 15–9 types were assigned using the MIRU-VNTRplus web application [34] [35] based on the MIRU-VNTRplus database. A cluster was defined as two or more strains having the same MLVA type.

Results

PZA DST results were obtained for 57 strains (Fig. 1). Of these 57.9% (33/57) were resistant and 42.1% (24/57) were susceptible to PZA. Of the tested 57 strains 63.1% (36/57) were MDR-TB and 7.0% (4/57) were XDR-TB by DST (Table 1). These 40 M(X)DR-TB samples represent 54.0% (40/74) of all the M(X)DR-TB isolates identified in a the larger collection of 399 strains [30]. Of the remaining 17 strains nine were monoresistant to isoniazid and seven were pan-susceptible and one strain was streptomycin mono-resistant. Of all 57 strains 24 (24/57, 42.1%) were from patients with a previous TB treatment history and 33 (33/57, 57.9%) were from newly diagnosed patients. Of the 24 retreatment cases 17 (17/24, 70%) were PZA resistant and of the 33 new patients 16 (16/33, 48%) were PZA resistant.

Fig. 1
figure 1

Flowchart of available information. XDR = extensively drug resistant, MDR = multidrug resistant tuberculosis, INH-R = isoniazid monoresistant, STR-R = streptomycin mono resistant. a = [30]

Table 1 Summary of pncA polymorphisms of all 57 strains with interpretable PZA DST results

A mutation in pncA was found in all 33 isolates that were identified as PZA resistant by DST (Table 1). In total 18 different mutations were identified in the 33 PZA resistant strains, of which 17 were non-synonymous mutations and one strain was identified with a 1-bp insertion resulting in a frameshift. No pncA mutations were observed in 23 of 24 PZA susceptible strains. In one of the 24 PZA susceptible isolates a point mutation at amino acid position 16 (ATC > CTG) was identified leading to an I6L substitution. Of the 18 polymorphisms identified in pncA 16 have been described previously in the literature (Fig. 2 and Table 1).

Fig. 2
figure 2

Distribution of pncA mutations identified in PZA resistant strains by sequencing and number of mutants per nucleotide polymorphisms

Among the group of 57 strains with interpretable PZA DST results, four MLVA clusters were identified two of which contained members with different pncA mutations (Fig. 3 and Table 4). The largest cluster, MLVA type 94–32, consisted of 16 strains. Of these, 6/16 (37.5%) were PZA susceptible and had no mutation in pncA; eight strains were phenotypically PZA resistant and four had a codon 212 mutation leading to a CAT212CCT (4/16, 25.0%) and four a CAT212CGT (4/16, 25.0%) change. These eight isolates represented 45.0% (8/17) of MLVA 94–32 strains among all sequenced isolates. Two other mutations were identified within this cluster (TGG356 TTG, ACG425ATG) each in a single isolate (1/16, 6.2%).

Fig. 3
figure 3

Distribution of pncA polymorphisms among the two largest MLVA clusters 94–32 and 100–32. Among all 57 strains with interpretable PZA DST results four clusters (94–32, 100–32, 407–32 and 1065–32) were identified. Clusters 407–32 and 1065–32 each consisted of two members, all carrying unique pncA mutations and are therefore not depicted in this figure. Distribution of strains with mutations in the pncA gene and their respective percentage in the cluster are listed. The legend indicates the position of the mutation in the pncA gene and the nucleotide change. In cluster 100–32 all but mutation A422C are present in individual isolates. No mutations are shared between the two clusters

The MLVA cluster 100–32 was formed by 14 isolates of which 13 were phenotypically resistant to PZA. One isolate was PZA susceptible by DST and did not carry a mutation in pncA. Seven of the 13 PZA resistant isolates (7/13, 53.8%) carried the same pncA mutation (CAG422CCG) and six isolates each (1/13, 7.8%) had different mutations in pncA (Table 1).

The CAG422CCG mutation resulting in a Gln141Pro substitution was observed in three MLVA types, 100–32, 100–35 and 9368–32 with the MLVA type 100–32 having the most representatives: this mutation was found in 50% (7/14) of all sequenced 100–32 strains.

MLVA clusters 407–32 and 1065–32 were both formed by two isolates each (Table 1 and Table 4). The strain carrying the I6L substitution (12–17,889) was one of two XDR-TB strains in the MLVA 1065–32 cluster. Unlike the 12–17,889 isolate the other member of the cluster, isolate 12–15,893, was PZA resistant by DST. PZA resistance was most likely conferred in 12–15,893 by a mutation in the pncA gene at position GTC20GCC leading to a V7A substitution. The I6L mutation, found in the isolate 12–17,889 has been previously reported in a large set of highly clustered isolates, called clade A, from Russia [26] together with four other characteristic genetic markers (rpoB E761D, Peis c-37a, pks15 T46I, Rv3254 P271P) [26]. In order to investigate whether these two isolates where part of the clade A cluster the genome sequences of isolates 12-15893 and 12-17889 were compared to the clade A strains from Russia. All five clade A specific markers were identified in the genome of strain 12–17,889, the isolate carrying the I6L mutation, but not in 12-15893. Genetic comparison revealed 43 SNPs difference between the two isolates (12-15893 and 12-17889) from this study and four SNPs difference between 12-17889 and the closest neighbour within the previously determined phylogeny of clade A (Fig. 2 in [26].

Discussion

The high diversity of PZA resistance-conferring mutations reported in TB patient isolates as well as the absence of confidence in phenotypic testing challenge routine PZA resistance testing. In addition, it is unclear what clinical implications there are for the detection of PZA resistance in vitro. Although, as PZA resistance is associated with worse outcome, information on (the epidemiology of) PZA resistance will undoubtedly become increasingly important. Molecular methods aiming to identify resistance-conferring mutations in pncA seem more reliable in identifying PZA resistance than complex DST methods, as the majority of mutations observed in pncA have been shown to confer phenotypic resistance, albeit in varying levels [10, 11, 36, 37]. Reliable detection of PZA resistance even seems achievable using next generation sequencing with reported sensitivity of between 70.9% to 97% and a specificity of 94% [38, 39].

Sixteen of the 18 pncA polymorphisms identified in our study have been described previously [10] and 10/18 and categorized as very high confidence resistance mutations. In our collection of isolates the correlation between phenotypic testing and genetic analysis was very good; only one isolate that carried a mutation in the pncA gene, strain 12-17889, was susceptible to PZA by conventional DST. This finding was corroborated by a recent study where the pncA-I6L, mutation was also associated with susceptibility to PZA in a panel of clustered isolates [26] and was not categorised as a high confidence resistance mutation [10]. Interestingly, strain 12-17889 was tested phenotypically resistant by an alternative method for PZA DST not using low pH [40].

Three pncA mutations, A212G, A212C and A422C, were observed in multiple isolates. The mutations A212G and A422C were previously categorised as very high confidence resistance mutations ([10], Table 1). These mutations were not only present in closely related (Beijing) strains (MLVA 94–32 and 100–32) but were also identified in non-clustered (Beijing) strains. This may indicate a mild propensity to fix mutations in this particular locus for specific (Beijing) strains, a phenomenon called epistasis that we previously investigated but were unable to confirm for rpoB mutations [22]. The two clusters identified in our dataset (MLVA 94–32 and 100–32) are also known for their epidemic potential outside of Georgia [41,42,43].

PZA resistance among new cases as well as clustering of PZA resistant strains with identical pncA mutations, as we found here (Table 1, Figs. 3 and 4), suggests that PZA resistant strains are transmissible. Contrary to the majority of first- and second line drug resistance mutations, multiple studies show that neither the fitness nor the virulence are significantly reduced by pncA mutations and data from Quebec moreover show that pncA mutants can maintain virulence and transmit [23, 44, 45]. However, as MIRU-VNTR typing has limited power to discriminate within genetically very similar strains only the comparison of the genomes of the respective strains can confirm our observation and rule out the possibility that the mutations in pncA have not been individually acquired.

Fig. 4
figure 4

UPGMA-tree of isolates with a PZA DST profile. Of all 57 isolates with a PZA DST results, 53 isolates with a valid MIRU-VNTR profile were included

Mutations that confer resistance to rifampicin or isoniazid and that result in fit transmissible bacterial strains are significantly more constrained. A very limited number of rare RIF and INH resistance mutations is found in genetically diverse clinical isolates. Consequentially, as genetically similar clinical isolates carry mostly identical MDR mutations (e.g. katG-S315 T and rpoB-S531 L) a proportion of isolates with identical MLVA types and these identical MDR mutations may be the result of independent mutational events rather than transmission of resistant strains. This is less likely for PZA where the available pool of fit mutations is much larger. Thus there might be a tendency to underestimate the de novo acquisition of rifampicin and isoniazid resistance. A recent study where longitudinal genomic comparison was performed supports the repeated acquisition theory for multiple drugs [25].

In addition to the large clusters, there were two smaller clusters in our collection, as identified by MLVA: 407–32 and 1065–32. Both clusters consisted of two members, all of which had acquired different pncA mutations, indicating independent events. One of the 1065–32 strains, 12-17889, carried an I6L mutation, which seems specific for a subgroup within the clade A MDR outbreak strain in the Samara region [26]. The other 1065–32 strain, 12–15,893, acquired resistance to PZA most likely via the V7A pncA mutation. Genomic comparison showed that there were only four SNPs difference between strain 12-17889 and the nearest neighbour of the clade A cluster, whereas there were 43 SNPs difference between strain 12–17,889 and strain 12-15893 indicating that they were much more distantly related.

The high prevalence of PZA resistant strains in Georgia, a highly endemic MDR(+) region, is a cause for concern as PZA resistance in MDR-TB cases is associated with unfavourable treatment outcomes. Bastos et al. showed that anti-TB treatment was more successful if patients infected with PZA susceptible MDR-TB received PZA, this association remained when they adjusted for fluoroquinolone and/or amikacin/kanamycin resistance [7]. Similar results were found in a study by Aung et al. showing that additional PZA resistance in the context of MDR-TB and fluoroquinolone resistance is associated with a worse treatment outcome [6]. In addition, PZA resistance may significantly jeopardise any PZA containing treatment regimens in the future. The sterilizing capacity of PZA as well as the synergistic action with drugs such as bedaquilline [3, 4] make it an attractive drug to include in both first- as well as second line treatment regimens.

This study has limitations. The isolates tested for PZA resistance were not specifically sampled for this study but were part of a previous study. We also did not sequence genes associated with PZA resistance other than pncA.

Conclusions

We report an alarmingly high percentage of PZA resistance in MDR-TB isolates in Georgia. In this study we see evidence of both the transmission of pncA resistant strains as well as the acquisition of PZA resistance in the patient. This needs to be confirmed by whole genome sequencing of the respective strains and should also be extended to a larger study population ideally collected within a larger time period. Our data contribute to the rising pool of evidence showing the high incidence of PZA resistance among MDR-TB isolates. Health authorities and TB control programs should consider prospective genotyping and PZA testing to assure current effective MDR-TB treatment and to inform the design of new MDR-TB treatment trials. Knowledge of PZA resistance in MDR-TB strains is needed to ensure sufficient active drugs are provided to protect against the emergence of additional resistance, particularly against newly introduced drugs.

Abbreviations

CTAB:

cetylmethylammonium bromide

DST:

drug susceptibility testing

INH:

isoniazid

MDR-TB:

multidrug-resistant tuberculosis

MGIT:

mycobacteria growth indicator tube

MIRU-VNTR:

mycobacterial interspersed repetitive unit-variable number of tandem repeats

MLVA:

multiple locus variable number of tandem repeat analysis

Mtb:

Mycobacterium tuberculosis

NRL:

National Reference Laboratory

PZA:

pyrazinamide

SNP:

single nucleotide polymorphism

STR:

streptomycin

TB:

tuberculosis

WGS:

whole genome sequencing

WHO:

World Health Organization

XDR-TB:

extensively drug-resistant tuberculosis

References

  1. Zhang Y, Mitchison D. The curious characteristics of pyrazinamide: a review. Int J Tuberc Lung Dis. 2003;7(1):6–21.

    CAS  PubMed  Google Scholar 

  2. Almeida Da Silva PE, Palomino JC. Molecular basis and mechanisms of drug resistance in mycobacterium tuberculosis: classical and new drugs. J Antimicrob Chemother. 2011;66(7):1417–30.

    Article  CAS  PubMed  Google Scholar 

  3. Diacon AH, Dawson R, von Groote-Bidlingmaier F, Symons G, Venter A, Donald PR, van Niekerk C, Everitt D, Hutchings J, Burger DA, et al. Bactericidal activity of pyrazinamide and clofazimine alone and in combinations with pretomanid and bedaquiline. Am J Respir Crit Care Med. 2015;191(8):943–53.

    Article  CAS  PubMed  Google Scholar 

  4. Diacon AH, Dawson R, von Groote-Bidlingmaier F, Symons G, Venter A, Donald PR, van Niekerk C, Everitt D, Winter H, Becker P, et al. 14-day bactericidal activity of PA-824, bedaquiline, pyrazinamide, and moxifloxacin combinations: a randomised trial. Lancet. 2012;380(9846):986–93.

    Article  CAS  PubMed  Google Scholar 

  5. Nunn AJ, Rusen ID, Van Deun A, Torrea G, Phillips PP, Chiang CY, Squire SB, Madan J, Meredith SK. Evaluation of a standardized treatment regimen of anti-tuberculosis drugs for patients with multi-drug-resistant tuberculosis (STREAM): study protocol for a randomized controlled trial. Trials. 2014;15:353.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Aung KJ, Van Deun A, Declercq E, Sarker MR, Das PK, Hossain MA, Rieder HL. Successful '9-month Bangladesh regimen' for multidrug-resistant tuberculosis among over 500 consecutive patients. Int J Tuberc Lung Dis. 2014;18(10):1180–7.

    Article  CAS  PubMed  Google Scholar 

  7. Bastos ML, Hussain H, Weyer K, Garcia-Garcia L, Leimane V, Leung CC, Narita M, Pena JM, Ponce-de-Leon A, Seung KJ, et al. Treatment outcomes of patients with multidrug-resistant and extensively drug-resistant tuberculosis according to drug susceptibility testing to first- and second-line drugs: an individual patient data meta-analysis. Clin Infect Dis. 2014;59(10):1364–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Hoffner S, Angeby K, Sturegard E, Jonsson B, Johansson A, Sellin M, Werngren J. Proficiency of drug susceptibility testing of mycobacterium tuberculosis against pyrazinamide: the Swedish experience. Int J Tuberc Lung Dis. 2013;17(11):1486–90.

    Article  CAS  PubMed  Google Scholar 

  9. Simons SO, van Ingen J, van der Laan T, Mulder A, Dekhuijzen PN, Boeree MJ, van Soolingen D. Validation of pncA gene sequencing in combination with the mycobacterial growth indicator tube method to test susceptibility of mycobacterium tuberculosis to pyrazinamide. J Clin Microbiol. 2012;50(2):428–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Miotto P, Cabibbe AM, Feuerriegel S, Casali N, Drobniewski F, Rodionova Y, Bakonyte D, Stakenas P, Pimkina E, Augustynowicz-Kopec E, et al. Mycobacterium tuberculosis pyrazinamide resistance determinants: a multicenter study. MBio. 2014;5(5):e01819–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ramirez-Busby SM, Valafar F. Systematic review of mutations in pyrazinamidase associated with pyrazinamide resistance in mycobacterium tuberculosis clinical isolates. Antimicrob Agents Chemother. 2015;59(9):5267–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Guidelines for treatment of tuberculosis. http://www.who.int/tb/publications/2010/9789241547833/en/ . Accessed 21 May 2017.

  13. Singh P, Mishra AK, Malonia SK, Chauhan DS, Sharma VD, Venkatesan K, Katoch VM. The paradox of pyrazinamide: an update on the molecular mechanisms of pyrazinamide resistance in mycobacteria. J Commun Dis. 2006;38(3):288–98.

    PubMed  Google Scholar 

  14. Pankhurst LJ, Del Ojo EC, Votintseva AA, Walker TM, Cole K, Davies J, Fermont JM, Gascoyne-Binzi DM, Kohl TA, Kong C, et al. Rapid, comprehensive, and affordable mycobacterial diagnosis with whole-genome sequencing: a prospective study. Lancet Respir Med. 2016;4(1):49–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Stoffels K, Mathys V, Fauville-Dufaux M, Wintjens R, Bifani P. Systematic analysis of pyrazinamide-resistant spontaneous mutants and clinical isolates of mycobacterium tuberculosis. Antimicrob Agents Chemother. 2012;56(10):5186–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ahmad N, Javaid A, Basit A, Afridi AK, Khan MA, Ahmad I, Sulaiman SA, Khan AH. Management and treatment outcomes of MDR-TB: results from a setting with high rates of drug resistance. Int J Tuberc Lung Dis. 2015;19(9):1109–14. i-ii

    Article  CAS  PubMed  Google Scholar 

  17. Chang KC, Yew WW, Zhang Y. Pyrazinamide susceptibility testing in mycobacterium tuberculosis: a systematic review with meta-analyses. Antimicrob Agents Chemother. 2011;55(10):4499–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Gunther G, van Leth F, Altet N, Dedicoat M, Duarte R, Gualano G, Kunst H, Muylle I, Spinu V, Tiberi S, et al. Beyond multidrug-resistant tuberculosis in Europe: a TBNET study. Int J Tuberc Lung Dis. 2015;19(12):1524–7.

    Article  CAS  PubMed  Google Scholar 

  19. Lemaitre N, Sougakoff W, Truffot-Pernot C, Jarlier V. Characterization of new mutations in pyrazinamide-resistant strains of mycobacterium tuberculosis and identification of conserved regions important for the catalytic activity of the pyrazinamidase PncA. Antimicrob Agents Chemother. 1999;43(7):1761–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Scorpio A, Lindholm-Levy P, Heifets L, Gilman R, Siddiqi S, Cynamon M, Zhang Y. Characterization of pncA mutations in pyrazinamide-resistant mycobacterium tuberculosis. Antimicrob Agents Chemother. 1997;41(3):540–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Van Rie A, Warren R, Mshanga I, Jordaan AM, van der Spuy GD, Richardson M, Simpson J, Gie RP, Enarson DA, Beyers N, et al. Analysis for a limited number of gene codons can predict drug resistance of mycobacterium tuberculosis in a high-incidence community. J Clin Microbiol. 2001;39(2):636–41.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bergval I, Kwok B, Schuitema A, Kremer K, van SD, Klatser P, Anthony R. pre-existing isoniazid resistance, but not the genotype of mycobacterium tuberculosis drives rifampicin resistance codon preference in vitro. PLoS One. 2012;7(1):e29108.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Nguyen D, Brassard P, Menzies D, Thibert L, Warren R, Mostowy S, Behr M. Genomic characterization of an endemic mycobacterium tuberculosis strain: evolutionary and epidemiologic implications. J Clin Microbiol. 2004;42(6):2573–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Thomas TA, Heysell SK, Houpt ER, Moore JL, Keller SJ. Outbreak of pyrazinamide-monoresistant tuberculosis identified using genotype cluster and social media analysis. Int J Tuberc Lung Dis. 2014;18(5):552–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Cohen KA, Abeel T, Manson McGuire A, Desjardins CA, Munsamy V, Shea TP, Walker BJ, Bantubani N, Almeida DV, Alvarado L, et al. Evolution of extensively drug-resistant tuberculosis over four decades: whole genome sequencing and dating analysis of mycobacterium tuberculosis isolates from KwaZulu-Natal. PLoS Med. 2015;12(9):e1001880.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Casali N, Nikolayevskyy V, Balabanova Y, Harris SR, Ignatyeva O, Kontsevaya I, Corander J, Bryant J, Parkhill J, Nejentsev S, et al. Evolution and transmission of drug-resistant tuberculosis in a Russian population. Nat Genet. 2014;46(3):279–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Eldholm V, Monteserin J, Rieux A, Lopez B, Sobkowiak B, Ritacco V, Balloux F. Four decades of transmission of a multidrug-resistant mycobacterium tuberculosis outbreak strain. Nat Commun. 2015;6:7119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. den Hertog AL, Sengstake S, Anthony RM: Pyrazinamide resistance in Mycobacterium tuberculosis fails to bite? Pathog Dis 2015, 73(6):ftv037.

  29. Tuberculosis surveillance and monitoring in Europe 2016. http://ecdc.europa.eu/en/publications/_layouts/forms/Publication_DispForm.aspx?List=4f55ad51-4aed-4d32-b960-af70113dbb90&ID=1452 . Accessed 21 May 2017.

  30. Tukvadze N., Bergval IL., Bablishvili N., Bzekalava N., Schuitema ARJ., de Beer J., de Zwaan R., Alba S., van Soolingen D., Aspindzelashvili R. et al: Evaluation of SNP-based genotyping to monitor tuberculosis control in a high MDR-TB setting. In.: bioRxiv; 2016.

  31. Coll F, McNerney R, Guerra-Assuncao JA, Glynn JR, Perdigao J, Viveiros M, Portugal I, Pain A, Martin N, Clark TG. A robust SNP barcode for typing mycobacterium tuberculosis complex strains. Nat Commun. 2014;5:4812.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Phelan J, Coll F, McNerney R, Ascher DB, Pires DE, Furnham N, Coeck N, Hill-Cawthorne GA, Nair MB, Mallard K, et al. Mycobacterium tuberculosis whole genome sequencing and protein structure modelling provides insights into anti-tuberculosis drug resistance. BMC Med. 2016;14:31.

    Article  PubMed  PubMed Central  Google Scholar 

  33. de Beer JL, Akkerman OW, Schurch AC, Mulder A, van der Werf TS, van der Zanden AG, van IJ, van SD: Optimization of standard in-house 24-locus variable-number tandem-repeat typing for mycobacterium tuberculosis and its direct application to clinical material. J Clin Microbiol 2014, 52(5):1338–1342.

  34. Weniger T, Krawczyk J, Supply P, Niemann S, Harmsen D: MIRU-VNTRplus: a web tool for polyphasic genotyping of Mycobacterium tuberculosis complex bacteria. Nucleic Acids Res 2010, 38(Web Server issue):W326-W331.

  35. Allix-Beguec C, Harmsen D, Weniger T, Supply P, Niemann S. Evaluation and strategy for use of MIRU-VNTRplus, a multifunctional database for online analysis of genotyping data and phylogenetic identification of mycobacterium tuberculosis complex isolates. J Clin Microbiol. 2008;46(8):2692–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Sheen P, Mendez M, Gilman RH, Pena L, Caviedes L, Zimic MJ, Zhang Y, Moore DA, Evans CA. Sputum PCR-single-strand conformational polymorphism test for same-day detection of pyrazinamide resistance in tuberculosis patients. J Clin Microbiol. 2009;47(9):2937–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Whitfield MG, Soeters HM, Warren RM, York T, Sampson SL, Streicher EM, van Helden PD, van Rie A. A global perspective on pyrazinamide resistance: systematic review and meta-analysis. PLoS One. 2015;10(7):e0133869.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Maningi NE, Daum LT, Rodriguez JD, Mphahlele M, Peters RP, Fischer GW, Chambers JP, Fourie PB. Improved detection by next-generation sequencing of pyrazinamide resistance in mycobacterium tuberculosis isolates. J Clin Microbiol. 2015;53(12):3779–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Coll F, McNerney R, Preston MD, Guerra-Assuncao JA, Warry A, Hill-Cawthorne G, Mallard K, Nair M, Miranda A, Alves A, et al. Rapid determination of anti-tuberculosis drug resistance from whole-genome sequences. Genome Med. 2015;7(1):51.

    Article  PubMed  PubMed Central  Google Scholar 

  40. den Hertog AL, Menting S, Pfeltz R, Warns M, Siddiqi SH, Anthony RM. Pyrazinamide is active against mycobacterium tuberculosis cultures at neutral pH and low temperature. Antimicrob Agents Chemother. 2016;60(8):4956–60.

    Article  CAS  Google Scholar 

  41. Devaux I, Kremer K, Heersma H, van Soolingen D. Clusters of multidrug-resistant mycobacterium tuberculosis cases, Europe. Emerg Infect Dis. 2009;15(7):1052–60.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Mokrousov I. Insights into the origin, emergence, and current spread of a successful Russian clone of mycobacterium tuberculosis. Clin Microbiol Rev. 2013;26(2):342–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Merker M, Blin C, Mona S, Duforet-Frebourg N, Lecher S, Willery E, Blum MG, Rusch-Gerdes S, Mokrousov I, Aleksic E, et al. Evolutionary history and global spread of the mycobacterium tuberculosis Beijing lineage. Nat Genet. 2015;47(3):242–9.

    Article  CAS  PubMed  Google Scholar 

  44. Cheng SJ, Thibert L, Sanchez T, Heifets L, Zhang Y. pncA mutations as a major mechanism of pyrazinamide resistance in mycobacterium tuberculosis: spread of a monoresistant strain in Quebec, Canada. Antimicrob Agents Chemother. 2000;44(3):528–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Yee DP, Menzies D, Brassard P. Clinical outcomes of pyrazinamide-monoresistant mycobacterium tuberculosis in Quebec. Int J Tuberc Lung Dis. 2012;16(5):604–9.

    CAS  PubMed  Google Scholar 

  46. Sandgren A, Strong M, Muthukrishnan P, Weiner BK, Church GM, Murray MB. Tuberculosis drug resistance mutation database. PLoS Med. 2009;6(2):e2.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

We thank N. Tukvadze, N. Bablishvili, N. Bzeklava and R. Aspindzelashvili (National TB Reference Laboratory, National Center for Tuberculosis and Lung Diseases, Tbilisi, Georgia) for collection of the strains, N.Casali and F. Drobniewski (Infectious Diseases and Immunity, Imperial College London, UK) for genomic comparison of our pncA–I6L mutant to the clade A cluster and A. Mulder (Mycobacteria Diagnostic Laboratory for Bacteriology and Parasitology Center for Infectious Disease Research, Diagnostics and Perinatal Screening National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands) for excellent technical assistance.

Funding

TGC is funded by the UK Medical Research Council (Grant no. MR/K000551/1, MR/M01360X/1, MR/N010469/1). JP is supported by a BBSRC PhD studentship. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and writing the manuscript.

Availability of data and materials

The sequence files of the two strains 12–17,889 (Accession: SRX2961214) and 12–15,893 (Accession: SRX2961215) have been uploaded to the Sequence Read Archive (http://trace.ncbi.nlm.nih.gov/Traces/sra/), NCBI’s repository for high-throughput sequencing data and are available under the submission number SRA581793 and Identifier: BioProject: PRJNA 392103. All other data supporting the conclusions of this article are included within the article.

Author information

Authors and Affiliations

Authors

Contributions

SS, ILB and RMA designed the study. SS, ILB, ARS, RMA, JLvB, RdZ, DvS, JP, TGC acquired, analysed and interpreted the data. SS, ILB and RMA wrote the first draft and adapted the manuscript. All authors participated in drafting the final version of the manuscript and gave final approval to the version to be published.

Corresponding author

Correspondence to Sarah Sengstake.

Ethics declarations

Ethics approval and consent to participate

The study was exempt from the approval of the Local Ethics Committee of the National Center for Tuberculosis and Lung Diseases (50 Maruashvili St, 0101 Tbilisi, Georgia) and Informed Consent was not required as the patient information used was anonymised before linking to the results of the analysis of the bacterial cultures and could not be linked back to individual patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sengstake, S., Bergval, I.L., Schuitema, A.R. et al. Pyrazinamide resistance-conferring mutations in pncA and the transmission of multidrug resistant TB in Georgia. BMC Infect Dis 17, 491 (2017). https://doi.org/10.1186/s12879-017-2594-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-017-2594-3

Keywords