Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic Mutations Associated with Isoniazid Resistance in Mycobacterium tuberculosis: A Systematic Review

  • Marva Seifert,

    Affiliation University of California San Diego, School of Medicine, La Jolla, California, United States of America

  • Donald Catanzaro,

    Affiliation University of California San Diego, School of Medicine, La Jolla, California, United States of America

  • Antonino Catanzaro,

    Affiliation University of California San Diego, School of Medicine, La Jolla, California, United States of America

  • Timothy C. Rodwell

    trodwell@ucsd.edu

    Affiliation University of California San Diego, School of Medicine, La Jolla, California, United States of America

Abstract

Background

Tuberculosis (TB) incidence and mortality are declining worldwide; however, poor detection of drug-resistant disease threatens to reverse current progress toward global TB control. Multiple, rapid molecular diagnostic tests have recently been developed to detect genetic mutations in Mycobacterium tuberculosis (Mtb) genes known to confer first-line drug resistance. Their utility, though, depends on the frequency and distribution of the resistance associated mutations in the pathogen population. Mutations associated with rifampicin resistance, one of the two first-line drugs, are well understood and appear to occur in a single gene region in >95% of phenotypically resistant isolates. Mutations associated with isoniazid, the other first-line drug, are more complex and occur in multiple Mtb genes.

Objectives/Methodology

A systematic review of all published studies from January 2000 through August 2013 was conducted to quantify the frequency of the most common mutations associated with isoniazid resistance, to describe the frequency at which these mutations co-occur, and to identify the regional differences in the distribution of these mutations. Mutation data from 118 publications were extracted and analyzed for 11,411 Mtb isolates from 49 countries.

Principal Findings/Conclusions

Globally, 64% of all observed phenotypic isoniazid resistance was associated with the katG315 mutation. The second most frequently observed mutation, inhA-15, was reported among 19% of phenotypically resistant isolates. These two mutations, katG315 and inhA-15, combined with ten of the most commonly occurring mutations in the inhA promoter and the ahpC-oxyR intergenic region explain 84% of global phenotypic isoniazid resistance. Regional variation in the frequency of individual mutations may limit the sensitivity of molecular diagnostic tests. Well-designed systematic surveys and whole genome sequencing are needed to identify mutation frequencies in geographic regions where rapid molecular tests are currently being deployed, providing a context for interpretation of test results and the opportunity for improving the next generation of diagnostics.

Introduction

While tuberculosis (TB) incidence and mortality declined over the past decade, there were still an estimated 8.6 million new cases of TB and 1.3 million deaths attributed to TB worldwide in 2012 [1]. Additionally, progress toward the reduction of TB morbidity and mortality has been largely overshadowed by the increasing prevalence of TB that fails to respond to standard antibiotic treatment. Multidrug resistant TB (MDR-TB), or TB that does not respond to either isoniazid (INH) or rifampicin (RIF), the “backbone” of the current recommended treatment, accounted for approximately 3.6% of all new TB cases and 20.2% of recurring, or previously treated TB cases reported in 2012 [1]. On a local level, incidence of MDR-TB is even more alarming, with individual countries detecting MDR-TB among as many as 28% of new cases and 62% of previously treated cases [2].

Drug resistance arises when mutations, or chromosomal replication errors, occur in genes that encode drug targets or drug metabolism mechanisms and impact the efficacy of anti-tuberculosis therapies [17]. Suboptimal treatment and poor adherence to drug therapies contribute to resistance amplification; while lack of detection and subsequent inappropriate treatment increase the risk of direct transmission of drug-resistant Mtb [5,8,9]. In 2012 it was estimated that less than 25% of MDR-TB cases were detected worldwide [1].

Low MDR-TB detection rates are due in part to the lack of rapid, accurate, and inexpensive diagnostic tests. The most commonly used method for detection of TB disease—sputum smear microscopy—does not detect drug resistance [1]. Standard drug susceptibility testing (DST), based on mycobacterial culture on solid or liquid media, can take weeks to months to produce results and requires significant biosafety resources, making it unfeasible in many low resource settings [1012]. Multiple genotypic DST methods based on the detection of mutations that confer resistance have been developed during the past decade, including DNA sequencing, line probe assays, hybridization on DNA chips, single-strand-conformation-polymorphisms, pyrosequencing, and real-time PCR [1214]. In order for these genotypic methods to accurately identify all resistant Mtb phenotypes, the genes and specific mutation locations associated with resistance must be known and included in the diagnostic test. However, global frequency and geographic distribution of INH resistance-conferring Mtb mutations have not yet been systematically quantified in the pathogen population [15].

Multiple reviews have identified genes that encode drug targets and have summarized the various mechanisms of resistance to both INH and RIF.[6,1620] While greater than 95% of RIF resistance is associated with mutations in an 81 base pair section of the rpoB gene, INH resistance appears more complex and has been associated with multiple genes, most frequently katG and inhA [17,2025]. Current rapid molecular diagnostic tests for INH resistance have focused on the detection of the “canonical” mutations in codon 315 of katG and position-15 in the inhA promoter region. Previous studies have identified highly variable frequencies of these mutations; with katG315 mutations accounting for 42 to 95% and inhA-15 mutations accounting for 6 to 43% of phenotypic INH resistance [6,19,20,26,27].

To date no systematic review has assessed the geographic variability of the most common mutations associated with INH resistance. Nor have the frequencies for multiple or co-occurring mutations been aggregated in order to understand the overall proportion of phenotypic INH resistance explained by the existing canonical mutations. It is critically important to understand the frequency and geographic distribution of mutations associated with INH resistance. A failure to account for these variations limits the local effectiveness of molecular diagnostic tools currently available and constrains the development of improved genotypic diagnostic tests [28]. The aims of our systematic review were to quantify the frequency and co-occurrence of the most common mutations associated with phenotypic INH drug resistance and to describe the regional differences in the distribution of these mutations as reported in the published literature.

Methods

A PubMed search was conducted to identify all English-language peer-reviewed publications that assessed mutations associated with INH resistance. In order to capture publications not aggregated in previous reviews, the search was restricted to studies published between January 1, 2000 and August 31, 2013. PubMed key search terms included: (INH OR isoniazid) AND (resistance or resistant) AND (tuberculosis) AND (mutation OR sequence).

Publication Selection Criteria

Publications were selected for inclusion if they met the following criteria: 1) presented original data; 2) used clinical strains of Mtb (not laboratory strains); 3) described the phenotypic DST method used as reference standard; 4) used DNA sequencing (including pyrosequencing) as a means of characterizing mutations; and 5) included individual level amino acid mutation data. Data from studies evaluating commercial or in-house molecular drug susceptibility tests were included if the results were verified by DNA sequencing.

Our primary goal was to assemble mutation data that was representative of the global distribution of drug resistance phenotypes. Studies were excluded if DNA sequencing was only performed on isolates with discrepancies between phenotypic and genotypic tests, as data from these studies would not be generalizable to the pathogen population. Studies were also excluded if DNA sequencing was performed on non-random isolate subsets or if studies reported sequential sequencing, i.e. sequencing katG on all isolates then sequencing inhA among isolates without katG mutations then sequencing ahpC-oxyR among isolates without mutations in either inhA or katG. An exception was made for studies where identifying sequencing results for the initial gene sequenced was easily discernible and in these instances, the sequencing results for only the first gene sequenced were included.

Data Extraction

The following information was extracted from each publication: primary author, publication year, geographic origin of specimens, year(s) of specimen collection, reference strains, phenotypic drug susceptibility testing method and cutoff, genotypic testing method, gene loci sequenced, primers used for DNA sequencing, and total number of resistant and susceptible isolates tested. Individual isolate mutation information included: location of gene mutation, amino acid and/or nucleotide changes, number of resistant and susceptible isolates with identified mutations, and number of resistant and susceptible isolates tested. Mutations in codon 463 of katG were not extracted as previous studies have identified katG463 as a genetic marker and not a marker of resistance [29,30]. All data were extracted and compiled using MS Excel software (Microsoft, Redmond, WA).

Mapping of Primer Sequences and Mutation Locations

In order to accurately assess which gene fragments had been sequenced for each individual isolate, the exact start and end points of the sequenced sections were determined using one of two methods. If primer sequences were included or referenced in the selected studies, sequence endpoint coordinates were identified by entering primer sequences into the NCBI Basic Local Alignment Search Tool (BLAST) and mapping the primers onto the complete Mtb H37Rv genome (Accession number NC_000962.3) [31]. If primer sequences were not published, sequenced endpoints coordinates were inferred by using the outermost identified mutations as sequence endpoints. If several primers were included and sequenced fragments overlapped, the final dataset included only the outermost/inclusive primer coordinates. Mutations identified in studies by traditional gene and nucleotide or codon nomenclature were converted to NCBI Mtb H37Rv complete genome (Accession number NC_000962.3) coordinates using the conversions described in S1 Table. In short, mutations occurring in the promoter or regulatory regions of the identified genes were mapped by exact nucleotide position and mutations occurring in the coding region of the gene were mapped by the center nucleotide position of the codon. Due to inconsistencies in the nomenclature and reporting of insertions, deletions, and frameshifts, they were excluded from analysis.

Calculation of Single Mutation Frequencies

As described in 2012 by Georghiou et al., mutation frequencies among phenotypically resistant isolates were calculated by aggregating the number of occurrences of a specific mutation among phenotypically resistant isolates and dividing it by the total number of phenotypically resistant isolates that could have identified that mutation, i.e. all phenotypically resistant isolates that had been sequenced at the location of the specific mutation [32]. Mutation frequencies among phenotypically susceptible isolates were calculated using the same method; the number of occurrences of a specific mutation among phenotypically susceptible isolates was divided by the total number of phenotypically susceptible isolates sequenced at the location of the mutation being assessed. Additionally, mutation frequency data were stratified by WHO region based on specimen origin and regional frequencies were calculated as described above [33].

Calculation of Cumulative Mutation Frequencies

Studies that reported co-occurring or multiple mutations in phenotypically resistant isolates were further analyzed to assess cumulative frequencies of the most commonly occurring mutations. Cumulative mutation frequencies were calculated first for individual genes, then across multiple genes simultaneously. Isolates were only included in data subsets if they had been sequenced at each of the mutation locations being evaluated. Individual mutation frequencies were reassessed in each data subset for global comparison, then the cumulative frequency for all mutations being evaluated was calculated by aggregating the number of isolates containing any of the mutations being assessed and dividing it by the total number of isolates in the subset.

Results

Of the 466 publications identified, 348 were excluded as they did not meet inclusion criteria (S1 Fig.). Mutation data, including gene, mutation location, original amino acid and/or nucleotide, and mutated amino acid and/or nucleotide, were extracted for 11,411 Mtb isolates described in the selected publications. Study size ranged from 4 to 502 isolates per publication. Seventy-seven percent (8,786) of the isolates were phenotypically resistant to INH. Geographic origin of the isolates was reported for 94% (10,745) of the isolates and included 49 countries (S4 Table). Primers were extracted (>96% of publications included primer information) for the three most commonly inspected genes and their surrounding regions: furA-katG, fabG1-inhA, and aphC-oxyR.

Mutations associated with INH resistance were reported in 15 unique gene regions, the most frequently being furA-katG (115 publications), fabG1-inhA (72 publications), and ahpC-oxyR (24 publications). Mutations were also reported in efpA, fadE24, iniA, iniB, iniC, kasA, nat, ndh, Rv1772, Rv1592c, Rv0340, and srmR genes, however data extracted for these genes were very limited and were excluded from further analysis.

Single Mutation Frequencies

Global frequencies for all mutations occurring in the katG coding region, inhA coding region, inhA promoter region, and ahpC-oxyR intergenic region were evaluated, however only mutations occurring at a cumulative frequency of >0.1% and identified in at least two publications are summarized in Table 1 (for a complete listing of all mutations see S2 Table). Of the 8,416 phenotypically resistant isolates sequenced at katG315, 5,400 (64.2%) harbored a mutation in katG315 whereas only 2 of 2,462 (<0.1%) phenotypically susceptible isolates harbored a mutation at that location. The remaining mutations identified in katG had frequencies of less than 1% among phenotypically resistant isolates. The second most commonly occurring mutation in INH resistant isolates, at position-15 of the inhA promoter, was identified in 1,189 of 6,192 (19.2%) phenotypically resistant and 7 of 2,063 (0.3%) phenotypically susceptible isolates. A second mutation in inhA promoter region, -8, was identified among 1.3% of phenotypically resistant isolates. Additionally, mutations occurring in the coding region of the inhA gene at codon positions 21, 94, and 194 were identified in approximately 1% of the isolates. Mutations in the ahpC-oxyR intergenic region reported at positions-10, -46, and -6 occurred in just over 1% of phenotypically resistant isolates. However, the mutation at position-46 of ahpC-oxyR intergenic region was also reported in 0.9% of susceptible isolates as well making it a non-specific marker of INH resistance.

thumbnail
Table 1. Global frequencies for selected mutations among phenotypically isoniazid resistant and phenotypically isoniazid sensitive isolates.

https://doi.org/10.1371/journal.pone.0119628.t001

Four-thousand, five-hundred and five isolates had nucleotide information reported in addition to amino acid information for the katG315 mutation (Table 2). The most common nucleotide change for the katG315 codon, AGC to ACC (serine to threonine), occurred among 93.4% of isolates. AGC to AAC (serine to asparagine), the second most common nucleotide change, occurred among 3.6% of isolates.

thumbnail
Table 2. Amino acid and nucleotide mutation frequency of 4,505 katG315 (wild type AGC) isolates.

https://doi.org/10.1371/journal.pone.0119628.t002

Mutation data at the country level was sparse, allowing for meaningful aggregation only at the regional level. Aggregating mutation frequencies by WHO regions revealed geographic differences in the individual frequencies for the two most commonly occurring mutations, katG315 and inhA-15 (Table 3). Among phenotypically INH resistant isolates, South East Asia had the highest frequency (78.4%) and the Western Pacific Region had the lowest frequency (55.5%) of the katG315 mutation. Frequencies also varied widely for inhA-15 mutation with the Americas region reporting the highest frequency (24.6%) and the Eastern Mediterranean region reporting the lowest frequency (13.0%) of this mutation in phenotypically INH resistant isolates.

thumbnail
Table 3. Regional distribution of the individual frequencies for katG315 and inhA-15.

https://doi.org/10.1371/journal.pone.0119628.t003

Cumulative Mutation Frequencies

The cumulative frequencies of multiple or co-occurring mutations associated with INH resistance, were assessed first by individual genes and then in combination across multiple gene regions (Table 4). Eighty-nine publications assessed 6,134 phenotypically resistant isolates for multiple mutations in katG inclusive of codon positions 306 to 316. The mutation frequency for katG315 in this subset of INH resistant isolates was 66.2%. After the inclusion of mutations at codons 311, 309, and 316 as markers of INH resistance, an additional 9 phenotypically resistant isolates were identified, for a total cumulative frequency of 66.3%. For the inhA promoter region, a total of 4,484 isolates from 51 publications were sequenced inclusive of positions -8 to -47. The cumulative frequency of inhA -15 in this subset was 19.0%, however unlike katG, other inhA promoter mutations did not typically co-occur with inhA -15. By including mutations at inhA -8, inhA -47, and inhA -17 as markers of resistance, an additional 72 (1.6%) phenotypically resistant isolates were identified. Seventeen publications, consisting of 1,696 phenotypically resistant isolates assessed mutations in the ahpC-oxyR intergenic region. In this subset no co-occurring mutations in the ahpC-oxyR intergenic region were observed. Mutations at positions -48, -39, -15, -12, -10, -9, and -6, of the ahpC-oxyR intergentic region occurred independently, resulting in a global cumulative frequency of 5.4% for ahpC-oxyR mutations among phenotypically resistant isolates in this subset.

thumbnail
Table 4. Single and cumulative mutation frequencies on katG, inhA, and ahpC-oxyR among data subsets which assessed co-occurring mutations.

https://doi.org/10.1371/journal.pone.0119628.t004

Forty-six studies, consisting of 4,179 phenotypically resistant isolates, assessed mutations in both the katG gene and the inhA promoter region, specifically encompassing katG315, inhA -15, inhA -8, and inhA -17. Frequencies for katG315 and inhA -15 in this subset were 65.2% and 19.0% respectively. KatG315 and inhA -15 mutations co-occurred in 181 isolates, resulting in a combined cumulative frequency for katG315 and inhA -15 of 79.9%. Only fifteen publications, consisting of 1,582 individual isolates, assessed mutations in katG coding region, inhA promoter region, and ahpC-oxyR intergenic region. In this subset, katG315 occurred among 66.8% of phenotypically resistant isolates and inhA -15 occurred among 16.2% of phenotypically resistant isolates. With the addition of inhA -8, inhA -47, and inhA -17, cumulative frequency reached 80.3%; and with the addition of the ahpC -10, ahpC -6, ahpC -39, ahpC -48, ahpC -15, ahpC -12, and ahpC -9 cumulative frequency reached 83.9%.

Discussion

Few studies to date have attempted to assess the prevalence of the most commonly occurring “canonical” mutations associated with phenotypic INH resistance among a globally representative set of isolates, making it difficult to compare the global frequencies calculated in this systematic review to previous global estimates. However, a study conducted by Lin et al. using 127 INH resistant isolates from California, a population which is thought to mirror global MDR-TB diversity due to immigration, identified a global frequency of 61% for katG315, 23% for inhA -15 mutations, and 83% for the cumulative frequency of either mutation, approximating the frequencies of these mutations as quantified in this systematic review [34]. In contrast, Campbell et al., using 212 INH resistant isolates from both WHO and CDC laboratory archives estimated the global frequency of the katG315 mutation to be 85%, inhA -15 to be 17%, and their cumulative frequency 91%; however, isolates used for that study were selected to provide a diverse set of mutation patterns, and therefore may not accurately represent true global frequencies [35]. Finally, a more recent study conducted by Rodwell et al. using 348 INH resistant isolates from four geographically diverse countries estimated the global frequencies of katG315 and inhA -15 to be significantly higher at 86% and 34% respectively. The cumulative frequency of both mutations with the addition of inhA -17 was 96% [36]. However, as with the previously mentioned study, isolates were selected to maximize diversity, potentially affecting the generalizability of the frequency estimates.

Prevalence of drug resistant TB is commonly understood to vary significantly by region, making it reasonable to expect regional variation in the frequencies of specific mutations that drive INH resistance [37,38]. Based on this systematic review, frequencies of specific mutations do in fact appear to vary quite significantly by geographic region. Although the total number of isolates analyzed varied by country and for some countries data was very sparse, we suspect that regional differences in the frequencies of the most common mutations could play a significant role in the performance of molecular-based diagnostic tests in some regions. For example, of the 808 INH resistant isolates sequenced from Japan and Korea, only 43% harbored the katG315 mutation, in contrast katG315 mutations in the post-soviet states of Belarus, Lithuania, Kazakhstan, Latvia, Moldova, and Russia were identified in 94% of the 751 INH resistant isolates analyzed from those countries. This regional variability of the katG315 mutation frequency may contribute to the poor performance of molecular diagnostics for INH resistance detection in Japan [39]. However, due to the lack of large, geographically diverse cross-sectional studies, regional patterns in single mutations should be interpreted with caution.

The frequency patterns of the most common mutations associated with INH resistance appear to differ between individual genes. It is clear that the overwhelming majority (64%) of phenotypic INH resistance among Mtb isolates is associated with a single mutation, katG315. The dominance of this mutation is hypothesized to be the result of a low or zero fitness cost for this mutation, allowing it to propagate without negative selection pressure [40,41]. Mutations other than katG315 in the katG gene appear to occur at low (<1%) frequencies and occur overwhelmingly in conjunction with the katG315 mutation. Patterns of co-occurring mutations in the inhA promoter region appear to differ markedly from co-occurring mutations in the katG gene. Although the inhA -15 mutation is the dominant (19%) mutation in the inhA promoter region, other resistance associated mutations (~1%) in the inhA promoter region appear to occur independently of the inhA -15 mutation and frequently contribute to the detection of INH resistance. In contrast, no single mutation in the ahpC-oxyR intergenic region showed an individual frequency above 1.3%; rather multiple mutations (at positions -48, -39, -15, -12, -10, -9, and -6,) in the ahpC-oxyR intergenic region occurred at a total cumulative frequency of 5.4% among phenotypically resistant isolates.

Equally important to our understanding of the distribution of INH resistance associated mutations is their cumulative frequency across genes. Analysis of the cumulative frequencies of mutations in katG, inhA, and ahpC-oxyR was performed only in a small subset of the data (n = 1,582), as few publications sequenced portions of all three genes of interest. In this subset, 80% of the isolates had at least one mutation reported at the following locations: katG315, inhA -15, inhA -8, inhA -47, or inhA -17. With the addition of mutations occurring at ahpC -48, ahpC -39, ahpC -15, ahpC -12, ahpC -10, ahpC -9, or ahpC -6, the percentage of isolates with at least one reported mutation increased to 84%. This approaches the cumulative frequency previously documented by Lin et al. of 88% [42]. The increased detection capacity achieved by including mutations in the ahpC-oxyR intergenic region appears to contradict previous studies which have found no significant contribution of INH resistance detection from mutations occurring in the ahpC-oxyR intergenic region [18,27,43].

Limitations

Given the nature of a systematic review, we were unable to obtain data from all geographic regions, thereby introducing a theoretical potential for bias in the global cumulative frequencies calculated. However, given that the total number of isolates reviewed exceeded 11,000, the cumulative frequencies most likely approximate true global frequencies of the most common mutations associated with INH resistance. Additionally, there is a potential for publications to selectively sequence isolates that did not show the most common mutations as detected on rapid molecular tests which may have led to artificially lower cumulative frequencies of the most common mutations. To reduce the potential for this reporting bias, we specifically excluded publications that indicated any non-random selection of isolates for sequencing analysis.

Conclusions

The performance of molecular-based diagnostic tests for drug resistant TB is intrinsically linked with not only the regional frequencies of mutations being detected, but also the diversity of resistance-conferring mutations being detected by diagnostic tests. Both of these factors influence the maximum sensitivity which rapid molecular tests can be expected to achieve.

Based on the data we evaluated, approximately 80% of global Mtb isolates with phenotypic resistance to INH appeared to contain mutations in codon 315 of the katG gene or position -15 in the inhA promoter. Taken together with other independently occurring mutations in the inhA promoter and the ahpC-oxyR intergenic region that have been previously documented to confer INH resistance, a minimum of 84% of all global Mtb isolates with INH resistance phenotypes should be detectable with molecular diagnostics based on analysis of these mutations. Additionally, as these mutations almost never (<0.1%) occurred in INH susceptible isolates, they should have specificities in excess of 99% as markers of phenotypic INH resistance. However, the cumulative frequencies of these mutations do appear to differ by region which could lead to variation in the sensitivity of molecular diagnostics if they are based only on these mutations. Further research, especially in regions with low prevalence of katG and inhA canonical mutations, is needed to determine if unknown genes or specific mutations account for the unexplained phenotypic INH resistance, or if a more inclusive combination of mutations in previously identified genes such as the inhA promoter, inhA coding, ahpC-oxyR intergenic region, or katG gene account for the remaining phenotypic resistance. Regional surveys of the most commonly occurring mutations, including mutations reported in the ahpC-oxyR intergentic region should continue to be conducted in geographic regions where rapid molecular tests are being deployed. This would allow for tailoring of molecular tests to specific regions, better interpretation of the molecular tests being used, and improved therapy recommendations.

Supporting Information

S1 Table. Calculations for H37Rv gene coordinate locations.

https://doi.org/10.1371/journal.pone.0119628.s003

(DOCX)

S2 Table. Frequency of reported mutations on furA-katG, fabG1-inhA, and aphC-oxyR gene regions.

https://doi.org/10.1371/journal.pone.0119628.s004

(DOCX)

S3 Table. Drug Susceptibility Testing (DST) Methods Employed in Publications.

https://doi.org/10.1371/journal.pone.0119628.s005

(DOCX)

S4 Table. Details of Studies Included in Review.

https://doi.org/10.1371/journal.pone.0119628.s006

(DOCX)

Acknowledgments

The authors would like to acknowledge Elizabeth Y. Kim for contribution during the initial stages of this study.

Author Contributions

Conceived and designed the experiments: TCR AC. Performed the experiments: MS. Analyzed the data: MS. Contributed reagents/materials/analysis tools: DC. Wrote the paper: MS TCR DC.

References

  1. 1. World Health Organization. Global Tuberculosis Report. Geneva, Switzerland; WHO press. 2013.
  2. 2. World Health Organization. Multidrug and Extensively Drug-Resistant TB 2010 Global Report of Surveillance and Response. Geneva, Switzerland; WHO press. 2010. pmid:21141336
  3. 3. Sougakoff W. Molecular epidemiology of multidrug-resistant strains of Mycobacterium tuberculosis. Clin Microbiol Infect. 2011; 17: 800–805. pmid:21682800
  4. 4. Middlebrook G. Isoniazid-resistance and catalase activity of tubercle bacilli; a preliminary report. Am Rev Tuberc. 1954; 69: 471–472. pmid:13138881
  5. 5. Mitchison DA. The diagnosis and therapy of tuberculosis during the past 100 years. Am J Respir Crit Care Med. 2005; 171: 699–706. pmid:15709048
  6. 6. Zhang Y, Yew WW. Mechanisms of drug resistance in Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2009; 13: 1320–1330. pmid:19861002
  7. 7. World Health Organization. Global Tuberculosis Report. Geneva, Switzerland; WHO Press. 2012.
  8. 8. Tiemersma EW, van der Werf MJ, Borgdorff MW, Williams BG, Nagelkerke NJ. Natural history of tuberculosis: duration and fatality of untreated pulmonary tuberculosis in HIV negative patients: a systematic review. PLoS One. 2011; 6: e17601. pmid:21483732
  9. 9. Hoek KG, Van Rie A, van Helden PD, Warren RM, Victor TC. Detecting drug-resistant tuberculosis: the importance of rapid testing. Mol Diagn Ther. 2011; 15: 189–194. pmid:21913741
  10. 10. Urbanczik R, Rieder HL. Scaling up tuberculosis culture services: a precautionary note. Int J Tuberc Lung Dis. 2009; 13: 799–800. pmid:19555526
  11. 11. World Health Organization. Treatment of tuberculosis: guidlines—4th ed. Geneva Switzerland; WHO Press. 2010.
  12. 12. Neonakis IK, Gitti Z, Krambovitis E, Spandidos DA. Molecular diagnostic tools in mycobacteriology. J Microbiol Methods. 2008; 75: 1–11. pmid:18599139
  13. 13. Alcaide F, Coll P. Advances in rapid diagnosis of tuberculosis disease and anti-tuberculous drug resistance. Enferm Infecc Microbiol Clin. 2011; 29 Suppl 1: 34–40. pmid:21420565
  14. 14. O'Grady J, Maeurer M, Mwaba P, Kapata N, Bates M, Hoelscher M, et al. New and improved diagnostics for detection of drug-resistant pulmonary tuberculosis. Curr Opin Pulm Med. 2011; 17: 134–141. pmid:21415753
  15. 15. Lessells RJ, Cooke GS, Newell ML, Godfrey-Faussett P. Evaluation of tuberculosis diagnostics: establishing an evidence base around the public health impact. J Infect Dis. 2011; 204 Suppl 4: S1187–1195. pmid:21996701
  16. 16. Scior T, Meneses Morales I, Garces Eisele SJ, Domeyer D, Laufer S. Antitubercular isoniazid and drug resistance of Mycobacterium tuberculosis—a review. Arch Pharm (Weinheim). 2002; 335: 511–525. pmid:12596216
  17. 17. Laurenzo D, Mousa SA. Mechanisms of drug resistance in Mycobacterium tuberculosis and current status of rapid molecular diagnostic testing. Acta Trop. 2011; 119: 5–10. pmid:21515239
  18. 18. Riska PF, Jacobs WR Jr., Alland D. Molecular determinants of drug resistance in tuberculosis. Int J Tuberc Lung Dis. 2000; 4: S4–10. pmid:10688142
  19. 19. Rattan A, Kalia A, Ahmad N. Multidrug-resistant Mycobacterium tuberculosis: molecular perspectives. Emerg Infect Dis. 1998; 4: 195–209. pmid:9621190
  20. 20. Rossetti ML, Valim AR, Silva MS, Rodrigues VS. [Resistant tuberculosis: a molecular review]. Rev Saude Publica. 2002; 36: 525–532. pmid:12364929
  21. 21. Zhang Y, Heym B, Allen B, Young D, Cole S. The catalase-peroxidase gene and isoniazid resistance of Mycobacterium tuberculosis. Nature. 1992; 358: 591–593. pmid:1501713
  22. 22. Vilcheze C, Wang F, Arai M, Hazbon MH, Colangeli R, Kremer L, et al. Transfer of a point mutation in Mycobacterium tuberculosis inhA resolves the target of isoniazid. Nat Med. 2006; 12: 1027–1029. pmid:16906155
  23. 23. Banerjee A, Dubnau E, Quemard A, Balasubramanian V, Um KS, Wilson T, et al. inhA, a gene encoding a target for isoniazid and ethionamide in Mycobacterium tuberculosis. Science. 1994; 263: 227–230. pmid:8284673
  24. 24. Telenti A, Imboden P, Marchesi F, Lowrie D, Cole S, Colston MJ, et al. Detection of rifampicin-resistance mutations in Mycobacterium tuberculosis. Lancet. 1993; 341: 647–650. pmid:8095569
  25. 25. Ramaswamy S, Musser JM. Molecular genetic basis of antimicrobial agent resistance in Mycobacterium tuberculosis: 1998 update. Tuber Lung Dis. 1998; 79: 3–29. pmid:10645439
  26. 26. Ahmad S, Mokaddas E. Recent advances in the diagnosis and treatment of multidrug-resistant tuberculosis. Respir Med. 2009; 103: 1777–1790. pmid:19660927
  27. 27. Sreevatsan S, Pan X, Zhang Y, Deretic V, Musser JM. Analysis of the oxyR-ahpC region in isoniazid-resistant and-susceptible Mycobacterium tuberculosis complex organisms recovered from diseased humans and animals in diverse localities. Antimicrob Agents Chemother. 1997; 41: 600–606. pmid:9056000
  28. 28. Abate G, Hoffner SE, Thomsen VO, Miorner H. Characterization of isoniazid-resistant strains of Mycobacterium tuberculosis on the basis of phenotypic properties and mutations in katG. Eur J Clin Microbiol Infect Dis. 2001; 20: 329–333. pmid:11453593
  29. 29. Haas WH, Schilke K, Brand J, Amthor B, Weyer K, Fourie PB, et al. Molecular analysis of katG gene mutations in strains of Mycobacterium tuberculosis complex from Africa. Antimicrob Agents Chemother. 1997; 41: 1601–1603. pmid:9210694
  30. 30. Baker LV, Brown TJ, Maxwell O, Gibson AL, Fang Z, Yates MD, et al. Molecular analysis of isoniazid-resistant Mycobacterium tuberculosis isolates from England and Wales reveals the phylogenetic significance of the ahpC-46A polymorphism. Antimicrob Agents Chemother. 2005; 49: 1455–1464. pmid:15793126
  31. 31. National Library of Medicine website. Basic Local Alignment Search Tool. Available: http://blast.ncbi.nlm.nih.gov/Blast.cgi. Accessed: 2015 Feb 13.
  32. 32. Georghiou SB, Magana M, Garfein RS, Catanzaro DG, Catanzaro A, Rodwell TC. Evaluation of genetic mutations associated with Mycobacterium tuberculosis resistance to amikacin, kanamycin and capreomycin: a systematic review. PLoS One. 2012; 7: e33275. pmid:22479378
  33. 33. World Health Organization website. World Health Organization Regional Offices. Available: http://www.who.int/about/regions/en/. Accessed: 2013 Feb 13.
  34. 34. Lin SY, Probert W, Lo M, Desmond E. Rapid detection of isoniazid and rifampin resistance mutations in Mycobacterium tuberculosis complex from cultures or smear-positive sputa by use of molecular beacons. J Clin Microbiol. 2004; 42: 4204–4208. pmid:15365012
  35. 35. Campbell PJ, Morlock GP, Sikes RD, Dalton TL, Metchock B, Starks AM, et al. Molecular detection of mutations associated with first- and second-line drug resistance compared with conventional drug susceptibility testing of Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2011; 55: 2032–2041. pmid:21300839
  36. 36. Rodwell TC, Valafar F, Douglas J, Qian L, Garfein RS, Chawla A, et al. Predicting extensively drug-resistant Mycobacterium tuberculosis phenotypes with genetic mutations. J Clin Microbiol. 2014; 52: 781–789. pmid:24353002
  37. 37. Jenkins HE, Zignol M, Cohen T. Quantifying the burden and trends of isoniazid resistant tuberculosis, 1994–2009. PLoS One. 2011; 6: e22927. pmid:21829557
  38. 38. Abal AT, Ahmad S, Mokaddas E. Variations in the occurrence of the S315T mutation within the katG gene in isoniazid-resistant clinical Mycobacterium tuberculosis isolates from Kuwait. Microb Drug Resist. 2002; 8: 99–105. pmid:12118524
  39. 39. Chikamatsu K, Mizuno K, Aono A, Yamada H, Sugamoto T, Nishiyama H, et al. [Evaluation of GenoType MTBDRplus for the detection of multi-drug-resistant Mycobacterium tuberculosis strains]. Kekkaku. 2011; 86: 697–702. pmid:21922778
  40. 40. Cohen T, Becerra MC, Murray MB. Isoniazid resistance and the future of drug-resistant tuberculosis. Microb Drug Resist. 2004; 10: 280–285. pmid:15650371
  41. 41. Ando H, Kondo Y, Suetake T, Toyota E, Kato S, Mori T, et al. Identification of katG mutations associated with high-level isoniazid resistance in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2010; 54: 1793–1799. pmid:20211896
  42. 42. Lin SY, Rodwell TC, Victor TC, Rider EC, Pham L, Catanzaro A, et al. Pyrosequencing for rapid detection of extensively drug-resistant Mycobacterium tuberculosis in clinical isolates and clinical specimens. J Clin Microbiol. 2014; 52: 475–482. pmid:24478476
  43. 43. Gagneux S, Burgos MV, DeRiemer K, Encisco A, Munoz S, Hopewell PC, et al. Impact of bacterial genetics on the transmission of isoniazid-resistant Mycobacterium tuberculosis. PLoS Pathog. 2006; 2: e61. pmid:16789833