Pyrazinamide (PZA) is an important first-line drug in the treatment of tuberculosis (TB) and of significant interest to the HIV-infected community due to the prevalence of TB-HIV coinfection in some regions of the world. The mechanism of resistance to PZA is unlike that of any other anti-TB drug. The gene pncA, encoding pyrazinamidase (PZase), is associated with resistance to PZA. However, because single mutations in PZase have a low prevalence, the individual sensitivities are low. Hundreds of distinct mutations in the enzyme have been associated with resistance, while some only appear in susceptible isolates. This makes interpretation of molecular testing difficult and often leads to the simplification that any PZase mutation causes resistance. This systematic review reports a comprehensive global list of mutations observed in PZase and its promoter region in clinical strains, their phenotypic association, their global frequencies and diversity, the method of phenotypic determination, their MIC values when given, and the method of MIC determination and assesses the strength of the association between mutations and phenotypic resistance to PZA. In this systematic review, we report global statistics for 641 mutations in 171 (of 187) codons from 2,760 resistant strains and 96 mutations from 3,329 susceptible strains reported in 61 studies. For diagnostics, individual mutations (or any subset) were not sufficiently sensitive. Assuming similar error profiles of the 5 phenotyping platforms included in this study, the entire enzyme and its promoter provide a combined estimated sensitivity of 83%. This review highlights the need for identification of an alternative mechanism(s) of resistance, at least for the unexplained 17% of cases.

We report the discovery and confirmation of 23 novel mutations with previously undocumented role in isoniazid (INH) drug resistance, in catalase-peroxidase (katG) gene of Mycobacterium tuberculosis (Mtb) isolates. With these mutations, a synonymous mutation in fabG1g609a, and two canonical mutations, we were able to explain 98% of the phenotypic resistance observed in 366 clinical Mtb isolates collected from four high tuberculosis (TB)-burden countries: India, Moldova, Philippines, and South Africa. We conducted overlapping targeted and whole-genome sequencing for variant discovery in all clinical isolates with a variety of INH-resistant phenotypes. Our analysis showed that just two canonical mutations (katG 315AGC-ACC and inhA promoter-15C-T) identified 89.5% of resistance phenotypes in our collection. Inclusion of the 23 novel mutations reported here, and the previously documented point mutation in fabG1, increased the sensitivity of these mutations as markers of INH resistance to 98%. Only six (2%) of the 332 resistant isolates in our collection did not harbor one or more of these mutations. The third most prevalent substitution, at inhA promoter position -8, present in 39 resistant isolates, was of no diagnostic significance since it always co-occurred with katG 315. 79% of our isolates harboring novel mutations belong to genetic group 1 indicating a higher tendency for this group to go down an uncommon evolutionary path and evade molecular diagnostics. The results of this study contribute to our understanding of the mechanisms of INH resistance in Mtb isolates that lack the canonical mutations and could improve the sensitivity of next generation molecular diagnostics.

Pyrazinamide (PZA) is an important first-line drug in all existing and new tuberculosis (TB) treatmentregimens. PZA-resistance in M. tuberculosis is increasing, especially among M/XDR cases. Noted issues with PZA Drug Susceptibility Testing (DST) have driven the search for alternative tests. This study provides a comprehensive assessment of PZA molecular diagnostics in M/XDR TB cases. A set of 296, mostly XDR, clinical M. tuberculosis isolates from four countries were subjected to DST for eight drugs, confirmatory Wayne's assay, and whole-genome sequencing. Three genes implicated in PZA resistance, pncA, rpsA, and panD were investigated. Assuming all non-synonymous mutations cause resistance, we report 90% sensitivity and 65% specificity for a pncA-based molecular test. The addition of rpsA and panD potentially provides 2% increase in sensitivity. Molecular heterogeneity in pncA was associated with resistance and should be evaluated as a diagnostic tool. Mutations near the N-terminus and C-terminus of PZase were associated with East-Asian and Euro-American lineages, respectively. Finally, Euro-American isolates are most likely to have a wild-type PZase and escape molecular detection. Overall, the 8-10% resistance without markers may point to alternative mechanisms of resistance. Confirmatory mutagenesis may improve the disconcertingly low specificity but reduce sensitivity since not all mutations may cause resistance. Pyrazinamide (PZA) is an important first-line drug recommended by the World Health Organization (WHO) for treatment of tuberculosis (TB)1 . Although most patients diagnosed with TB are successfully treated with first-line drugs, rifampicin (RIF), isoniazid (INH), ethambutol, and PZA, the incidence of multidrug resistant TB (MDR-TB) (resistant to first-line drugs) is increasing every year 1 . Along with this rate, the incidence of PZA resistance (PZA R ) is also increasing since over 60% of MDR-TB cases globally are also PZA R (rates vary regionally) 2 . In this manuscript we denote isolates that are resistant or susceptible to PZA by growth-based drug susceptibility testing (DST) as PZA R or PZA S , respectively. PZA is unique among anti-TB drugs in that, while it demonstrates powerful in vivo sterilizing activity, it exhibits no activity against actively growing Mycobacterium tuberculosis bacilli under normal culture conditions at neutral pH 3 . The drug appears to preferentially act against non-replicating persisters with low metabolic activity at acid pH in vitro or in vivo 4 . Importantly, most new drug regimens proposed for treatment of drug-resistant TB (DR-TB) show improved outcomes when combined with PZA 5-9 . Although growth-based PZA susceptibility testing is recommended by the WHO 10 , the method yields inconsistent results in some cases and suffers from a relatively high false resistance rate 11 .

This study assembles DNA adenine methylomes for 93 Mycobacterium tuberculosis complex (MTBC) isolates from seven lineages paired with fully-annotated, finished, de novo assembled genomes. Integrative analysis yielded four key results. First, methyltransferase allele-methylome mapping corrected methyltransferase variant effects previously obscured by reference-based variant calling. Second, heterogeneity analysis of partially active methyltransferase alleles revealed that intracellular stochastic methylation generates a mosaic of methylomes within isogenic cultures, which we formalize as ‘intercellular mosaic methylation’ (IMM). Mutation-driven IMM was nearly ubiquitous in the globally prominent Beijing sublineage. Third, promoter methylation is widespread and associated with differential expression in the ΔhsdM transcriptome, suggesting promoter HsdM-methylation directly influences transcription. Finally, comparative and functional analyses identified 351 sites hypervariable across isolates and numerous putative regulatory interactions. This multi-omic integration revealed features of methylomic variability in clinical isolates and provides a rational basis for hypothesizing the functions of DNA adenine methylation in MTBC physiology and adaptive evolution.

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.