Volume 16, Number 9—September 2010
Determinants of Multidrug-Resistant Tuberculosis Clusters, California, USA, 2004–2007
Laboratory and epidemiologic evidence suggests that pathogen-specific factors may affect multidrug-resistant (MDR) tuberculosis (TB) transmission and pathogenesis. To identify demographic and clinical characteristics of MDR TB case clustering and to estimate the effect of specific isoniazid resistance–conferring mutations and strain lineage on genotypic clustering, we conducted a population-based cohort study of all MDR TB cases reported in California from January 1, 2004, through December 31, 2007. Of 8,899 incident culture-positive cases for which drug susceptibility information was available, 141 (2%) were MDR. Of 123 (87%) strains with genotype data, 25 (20%) were aggregated in 8 clusters; 113 (92%) of all MDR TB cases and 21 (84%) of clustered MDR TB cases occurred among foreign-born patients. In multivariate analysis, the katG S315T mutation (odds ratio 11.2, 95% confidence interval 2.2–∞; p = 0.004), but not strain lineage, was independently associated with case clustering.
In 2007, >500,000 cases of multidrug-resistant (MDR) tuberculosis (TB), defined as resistance to at least isoniazid and rifampin, occurred worldwide (1). Although demographic and clinical risk factors for transmission and pathogenesis of both drug-susceptible and drug-resistant Mycobacterium tuberculosis have been well described (2,3), little is known about the microbial factors that influence the generation of secondary MDR TB cases (4,5).
Community- and population-based molecular epidemiologic studies of isoniazid-monoresistant M. tuberculosis (6–8) have shown that specific resistance-conferring mutations are associated with variable degrees of genotypic clustering, a measure of strain pathogenicity that incorporates host factors, transmissibility of the organism, and capacity of the organism to cause active disease. For example, isoniazid-monoresistant strains with a serine-to-threonine substitution at position 315 (S315T) are more often associated with secondary cases than are strains without the S315T mutation (6,7), likely because of reduced or absent catalase–peroxidase production (9). However, the effects of specific isoniazid resistance–conferring mutations on genotypic clustering in multidrug resistance are less well characterized. The studies reported to date have been limited by inadequate genotypic discrimination (10,11) and/or nonrepresentative sampling of cases (10,12–14).
California reports the highest annual number of TB cases (15), more than one fourth of all MDR TB cases (16), and the highest immigration rates in the United States (17). We conducted a population-based cohort study of all incident MDR TB cases in California during a 4-year period (January 2004–December 2007) to 1) describe demographic and clinical characteristics of clustering and 2) estimate the effect of specific isoniazid resistance–conferring mutations and strain lineage on genotypic clustering of MDR M. tuberculosis.
We analyzed culture-positive cases of MDR TB reported to the California TB registry from January 1, 2004, through December 31, 2007. California state law (Health and Safety Code Title 17 §2505) requires reporting of all verified cases of TB, submission of all M. tuberculosis isolates to local public health laboratories, and submission of all MDR M. tuberculosis isolates to the California Department of Public Health Microbial Diseases Laboratory. Testing for first- and second-line drug susceptibilities was performed at local laboratories or at the Microbial Diseases Laboratory by using BACTEC 460 (Becton Dickinson Diagnostic Instruments, Sparks, MD, USA), MGIT 960 (Becton Dickinson), or the agar proportion method. Some isolates were forwarded to the Centers for Disease Control and Prevention (CDC; Atlanta, GA, USA) for additional second-line drug testing or for confirmation of drug resistance. Patients were included in the study if MDR M. tuberculosis was identified on >1 isolate. Demographic and clinical information for all patients with MDR TB was abstracted from state TB surveillance forms (Report of a Verified Case of Tuberculosis), which have high validity (18). All protocols were approved by the California Committee for the Protection of Human Subjects and University of California, San Francisco, Committee for the Protection of Human Subjects.
Characterization of Mutations Associated with Isoniazid and Rifampin Resistance
For each isolate, genomic DNA was extracted from solid media (Lowenstein-Jensen slants, Middlebrook 7H10 or 7H11 agar), liquid media (BACTEC 12 B or MGIT [Becton Dickinson]), or smear-positive sputum sediments. A real-time PCR assay with 6 molecular beacon probes was performed by using an iQ5 iCycler instrument (Bio-Rad, Hercules, CA, USA) to screen for mutations associated with isoniazid and rifampin resistance (19). Two molecular beacons that targeted katG (codon 311–317) and the inhA promoter were used to detect isoniazid resistance–conferring mutations, and 4 molecular beacons that targeted the core region of rpoB were used to detect rifampin resistance–conferring mutations. Isolates with mutations in katG detected by the wild-type probe were further tested with another molecular beacon that specifically targeted katG S315T (AGC–ACC). When molecular beacon analysis did not show katG S315T or -c15t inhA promoter mutations, the entire furA-katG locus (H37Rv: 2153626–2156657, 3,031 bp) was sequenced as described (6). Sequence data were generated by using ABI BigDye v3.1 dye terminator sequencing chemistry and the ABI PRISM 3730xl capillary DNA analyzer (Applied Biosystems, Foster City, CA, USA) at the Genomic Core Facility, University of California, San Francisco (www.genomics.ucsf.edu/Sequencing/index.aspx), and were analyzed with ClustalW (www.ebi.ac.uk/Tools/clustalw/index.html).
Genotyping and Lineage Determination
Spacer oligonucleotide typing (spoligotyping) and mycobacterial interspersed repetitive unit (MIRU) typing were performed in accordance with the Centers for Disease Control and Prevention Universal Genotyping Program procedures (20). Spoligotyping was performed by using Luminex-based methods to detect 43 known spacer sequences in the direct repeat locus (21). MIRU typing was performed by using the protocol described by Cowan et al. (22). A capillary sequencer, CEQ 8000 (Beckman, Fullerton, CA, USA), was used to analyze the number of repeated sequences at each of the 12 loci. IS6110-based restriction fragment length polymorphism (RFLP) genotyping was performed following standardized methods (23). RFLP patterns were compared by using Bioimage Whole Band Analyzer software version 4.2.1 (Bioimage Corp., Ann Arbor, MI, USA) (24). RFLP patterns with <20 identical bands, or >20 bands and differing by no more than a single band, were considered matched. IS6110 RFLP band assignment was edited by 2 independent readers, and the cluster assignment was confirmed visually.
The phylogeographic lineage of strains was determined from spoligotyping results. Spoligotype families H, LAM, and T, X, S were considered to be of Euro-American lineage; Beijing of East-Asian lineage; EAI of Indo-Oceanic lineage; and CAS of East African-Indian lineage (5).
Cases were defined as clustered if >2 isolates from cases reported during the study period shared the same MIRU and spoligotype, matched IS6110 RFLP, and had specific drug resistance–conferring mutations for isoniazid and rifampin. Clustering was assumed to represent both transmission of M. tuberculosis and progression to active disease, leading to secondary case generation within the period of the study. Cases not in a cluster were considered to be the result of reactivation of latent infection. Patients with the earliest case report date within a cluster were regarded as index cases.
The proportion of clustered MDR TB cases was analyzed as the number of clustered cases divided by the total number of culture-positive cases that occurred during the study period. Because of limited sample size, the independent effects of katG S315T and phylogeographic lineage on clustering were estimated by using exact logistic regression methods. Refugee status and sputum smear positivity were included in the model as relevant host risk factors. Refugee resettlement during the study period could bias our results in that overrepresentation of ethnic groups or geographic locations with a high prevalence of particular strain-specific factors (phylogeographic lineage or drug-resistance mutations) could confound the association under study. To examine this influence and potential clustering of TB cases within households or communities related to refugee resettlement, sensitivity analysis was conducted by reestimation of study results after excluding 1) patients immigrating from refugee settings within the past 3 years and 2) the single largest patient cluster, which accounted for 40% of all clustered cases. In a separate analysis, standard logistic regression was used to estimate the effect of katG S315T and phylogeographic lineage on sputum smear positivity, again controlling for refugee status.
Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to measure associations of interest. Categorical data (e.g., sex, foreign vs. US birth, homelessness) were compared by using Fisher exact tests. The Wilcoxon rank-sum test was used to determine differences in the distribution of continuous variables (e.g., age, time from entry into the United States to TB diagnosis). Interaction between the independent variables was assessed separately for each factor by stratification and statistical testing (Breslow-Day with the Tarone correction  and Zelen test ), with p<0.2 assumed to indicate the presence of interaction. All p values were 2-sided with α = 0.05 as the significance level. Analyses were performed by using Stata 10 (StataCorp., College Station, TX, USA) and StatExact 8 (CYTEL Software Corp., Cambridge, MA, USA).
During the study period, 11,395 cases of TB were reported in California, of which 9,037 (79%) had positive cultures. Of these, 8,899 (98%) had isoniazid and rifampin susceptibility results available. Of the 141 (2%) incident MDR TB cases, 123 (87%) had isolates available for MIRU, spoligotyping, and IS6110 RFLP analysis. Isolates unavailable for genotyping (n = 18) were more often from Los Angeles County; other demographic and clinical characteristics were similar to those of analyzed cases.
Twenty-five MDR M. tuberculosis isolates were aggregated in 8 clusters (1 cluster of 10 cases, 1 cluster of 3 cases, and 6 clusters of 2 cases) for an overall cluster proportion of 20% (25/123). Excluding the 8 index cases, 14% (17/123) of all cases were considered to have resulted from recent transmission and rapid progression to disease (secondary cases). Within clusters, a median of 3 months elapsed between successive secondary cases (range 0–20 months). Of the 123 total cases, 113 (92%) occurred among foreign-born patients, with more than half (56%) occurring among immigrants from Mexico, Philippines, the People’s Republic of China, or Vietnam (Table 1). Seven of 8 index cases occurred either in Mexican immigrants (3/8) or recent refugees from Thailand, Lao People’s Democratic Republic, or India (4/8). Median time from patient arrival in the United States to TB diagnosis was approximately twice as long for clustered as for nonclustered cases (4.3 years vs. 2.4 years) and 3 times as long for cases in Mexican-born persons as for cases in persons from other countries (7.3 years vs. 2.3 years).
Younger persons were more likely than older persons to harbor strains involved in MDR TB clusters (Table 1). HIV infection was unusual in this patient population; only 3 (4%) of 75 patients with known HIV status were HIV-infected. Twenty-eight percent (35/123) of patients reported a history of active TB; this proportion did not vary between clustered and nonclustered cases (p = 0.75). Eight patients had documented previous treatment in California. Time to culture negativity (2.4 vs. 2.8 months, p = 0.95), treatment failure, or death did not differ between clustered and nonclustered cases (p = 0.57) among 105 (85%) of 123 cases for which data were available.
Seventy-five percent of MDR M. tuberculosis strains harbored the isoniazid resistance–conferring mutation katG S315T, including all (100%) clustered strains (Table 2). When we controlled for strain lineage and refugee status, katG S315T was inversely associated (OR 0.28, 95% CI 0.09–0.89, p = 0.03) with sputum smear positivity (Table 3). Most rifampin resistance–conferring mutations were found between positions 529 and 534 of the rpoB gene, likely indicative of the serine to leucine substitution at position 531 (S531L) of the rpoB gene (27). The association between this mutation and clustering did not reach statistical significance (OR 2.2, 95% CI 0.8–7.4; p = 0.16).
MDR M. tuberculosis isolates were distributed among East-Asian (47%), Euro-American (30%), Indo-Oceanic (20%), and East African–Indian (3%) phylogeographic lineages. Lineage could not be established for 11 cases. On univariate analysis, East-Asian strain lineage was associated with clustering (Table 4), but not with adverse outcome (death or treatment failure). Indo-Oceanic strains produced no secondary cases.
Clustering was independently associated with katG S315T (OR 11.2, 95% CI 2.2-∞; p = 0.004) and refugee status (OR 6.0, 95% CI 1.2–36.2; p = 0.03) in exact multivariate analyses in which strain lineage and sputum smear positivity were controlled for (Table 5). The estimated association between katG S315T and case clustering did not appreciably change in sensitivity analyses that excluded either all recently arrived refugees or the single largest patient cluster.
In this 4-year population-based molecular epidemiologic study, transmission followed by secondary case generation contributed to ≈14%, or 1 of every 7, MDR TB cases in California. Clustered cases occurred more often among younger persons and persons who had emigrated from Mexico and refugee settings in Southeast Asia. In addition, pathogen-specific factors were associated with clustering of MDR TB cases, independent of traditional clinical and demographic risk factors.
The proportion of MDR TB cases attributed to transmission in California was higher than that reported by other investigators in most low incidence settings (28–32). The largest clusters in our study resulted from MDR TB outbreaks in California after resettlement of Hmong refugees in 2005–2006 (33) and resettlement of Tibetan refugees in 2001–2006. The associations between pathogen-specific factors and case clustering could be due to regional differences in strain prevalence and preferential migration of persons with specific strains to California. We attempted to control for these factors by using highly stringent criteria to define clustered cases and by adjusting for refugee status in our multivariate model. However, without detailed contact information and contact tracing, we cannot be certain of the extent to which transmission or progression to active disease occurred within or outside California. Given that most clustered cases occurred among persons residing in the United States for >3 years and that US-born persons were involved in 3 of 8 clusters, at least some proportion of MDR TB transmission is likely to have occurred in California. This observation suggests that although most MDR TB cases in the United States are related to the migration of persons already infected with drug-resistant M. tuberculosis, a small but notable proportion may be due to ongoing transmission.
Heterogeneity in the reproductive success of drug-resistant M. tuberculosis is now well established (34,35). In this study, the katG S315T mutation was the only isoniazid resistance–conferring mutation found among clustered MDR TB cases. The high prevalence of katG S315T among MDR strains (12,14,36–38) and the association of this mutation with increased secondary case generation among isoniazid-monoresistant strains (6,7) have been documented. We report that the katG S315T isoniazid resistance–conferring mutation retains an independent effect on clustering of MDR TB cases despite the presence of mutations that confer resistance to additional drugs. In particular, secondary case generation did not significantly vary according to site of rpoB mutations that confer rifampin resistance, which supports the hypothesis that these are no-cost mutations or that compensatory mutations commonly exist (39).
The katG S315T mutation is thought to preserve fitness through the relative preservation of catalase-peroxidase production (9), although whether this mutation is associated with different clinical phenotypes is unknown. In this study, we noted an inverse association between presence of the katG S315T mutation and sputum smear positivity. In addition, the katG S315T mutation was associated with case clustering, independent of sputum smear status. These findings suggest that the katG S315T mutation may preserve the ability of M. tuberculosis to transmit and cause secondary cases through mechanisms unrelated to conventional indices of disease severity, such as the presence of abundant acid-fast bacilli in sputum.
Our findings are clinically useful for at least 2 reasons. First, MDR TB in California is occurring predominantly among patients who were not born in the United States, with some cases from recent transmission and rapid progression to disease. Our study suggests that in California, younger persons and persons who have emigrated from Mexico and from refugee settings may be at higher risk for transmitting MDR M. tuberculosis. Likewise, our findings reinforce the need for giving priority to screening and prevention activities in immigrant communities and US investment in international TB control (40). Second, if our results are verified in other settings, TB-control programs should consider pathogen-specific factors such as isoniazid resistance–conferring mutations when planning the intensity of contact investigation and secondary case-finding activities.
This study has several limitations. First, our estimates of case clustering are imprecise because of the limited number of MDR TB cases observed during the study period. However, these estimates are the best currently available, given that the data make up the largest population-based TB registry in the United States. Second, although our definition of genotypic clustering was highly rigorous, the lack of detailed epidemiologic information precluded confirmation of transmission within California. Third, because our study did not include pan-susceptible or isoniazid-monoresistant strains, we cannot comment directly on MDR M. tuberculosis pathogenicity relative to these groups. Lastly, our analyses implicitly assume independence of outcome events, and household or community-level factors potentially associated with clustering were not available. Future studies should be designed so that statistical methods can be used that are able to accommodate the possible effects of within-household clustering.
We found a substantial proportion of MDR TB cases and case clustering in California among non–US-born persons, and the katG S315T mutation was independently associated with clustering. Validation of these findings in larger cohorts and in different population settings may have crucial public health consequences.
Dr Metcalfe is a senior fellow in the Division of Pulmonary and Critical Care Medicine at University of California San Francisco and a doctoral candidate in epidemiology at University of California Berkeley. His research interests include using prediction modeling to identify at-risk populations and to evaluate novel diagnostics.
We thank Cynthia Sanders and Ed Desmond for supplying the M. tuberculosis isolates and James Watt and Julia Hill for their editorial assistance.
The study was funded by grant numbers T32 HL007185 (J.Z.M.), K23 HL094141 (A.C.), and R01 AI034238 (P.C.H., M.K.M.) from the National Institutes of Health.
- World Health Organization. Multidrug and extensively drug-resistant TB (M/XDR-TB): 2010 global report on surveillance and response [cited 2010 Jun 21]. http://www.who.int/tb/features_archive/world_tb_day_2010/en/index.html
- Nava-Aguilera E, Andersson N, Harris E, Mitchell S, Hamel C, Shea B, Risk factors associated with recent transmission of tuberculosis: systematic review and meta-analysis. Int J Tuberc Lung Dis. 2009;13:17–26.
- Kliiman K, Altraja A. Predictors of extensively drug-resistant pulmonary tuberculosis. Ann Intern Med. 2009;150:766–75.
- Dye C. Doomsday postponed? Preventing and reversing epidemics of drug-resistant tuberculosis. Nat Rev Microbiol. 2009;7:81–7.
- Gagneux S, Small PM. Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect Dis. 2007;7:328–37.
- Gagneux S, Burgos MV, DeRiemer K, Encisco A, Munoz S, Hopewell PC, Impact of bacterial genetics on the transmission of isoniazid-resistant Mycobacterium tuberculosis. PLoS Pathog. 2006;2:e61.
- van Doorn HR, de Haas PE, Kremer K, Vandenbroucke-Grauls CM, Borgdorff MW, van Soolingen D. Public health impact of isoniazid-resistant Mycobacterium tuberculosis strains with a mutation at amino-acid position 315 of katG: a decade of experience in the Netherlands. Clin Microbiol Infect. 2006;12:769–75.
- van Soolingen D, de Haas PE, van Doorn HR, Kuijper E, Rinder H, Borgdorff MW. Mutations at amino acid position 315 of the katG gene are associated with high-level resistance to isoniazid, other drug resistance, and successful transmission of Mycobacterium tuberculosis in the Netherlands. J Infect Dis. 2000;182:1788–90.
- Pym AS, Saint-Joanis B, Cole ST. Effect of katG mutations on the virulence of Mycobacterium tuberculosis and the implication for transmission in humans. Infect Immun. 2002;70:4955–60.
- Ang CF, Ong CS, Rukmana A, Pham Thi KL, Yap SF, Ngeow YF, An overview of the phenotypic and genotypic characteristics of multidrug-resistant Mycobacterium tuberculosis isolates from four Asian countries. J Med Microbiol. 2008;57:1039–40.
- Tracevska T, Jansone I, Baumanis V, Marga O, Lillebaek T. Prevalence of Beijing genotype in Latvian multidrug-resistant Mycobacterium tuberculosis isolates. Int J Tuberc Lung Dis. 2003;7:1097–103.
- Lipin MY, Stepanshina VN, Shemyakin IG, Shinnick TM. Association of specific mutations in katG, rpoB, rpsL and rrs genes with spoligotypes of multidrug-resistant Mycobacterium tuberculosis isolates in Russia. Clin Microbiol Infect. 2007;13:620–6.
- Marttila HJ, Soini H, Eerola E, Vyshnevskaya E, Vyshnevskiy BI, Otten TF, A Ser315Thr substitution in katG is predominant in genetically heterogeneous multidrug-resistant Mycobacterium tuberculosis isolates originating from the St. Petersburg area in Russia. Antimicrob Agents Chemother. 1998;42:2443–5.
- Valvatne H, Syre H, Kross M, Stavrum R, Ti T, Phyu S, Isoniazid and rifampicin resistance-associated mutations in Mycobacterium tuberculosis isolates from Yangon, Myanmar: implications for rapid molecular testing. J Antimicrob Chemother. 2009;64:694–701.
- Centers for Disease Control and Prevention. Trends in tuberculosis—United States, 2008. MMWR Morb Mortal Wkly Rep. 2009;58:249–53.
- Granich RM, Oh P, Lewis B, Porco TC, Flood J. Multidrug resistance among persons with tuberculosis in California, 1994–2003. JAMA. 2005;293:2732–9.
- Passel JS, Cohn DV. A portrait of unauthorized immigrants in the United States. Washington: Pew Hispanic Center; April 2009; pp. 2, 13.
- Sprinson JE, Lawton ES, Porco TC, Flood JM, Westenhouse JL. Assessing the validity of tuberculosis surveillance data in California. BMC Public Health. 2006;6:217.
- 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–8.
- National TB Controllers Association/CDC Advisory Group on Tuberculosis Genotyping. Guide to the application of genotyping to tuberculosis prevention and control. Atlanta (GA): US Department of Health and Human Services, Centers for Disease Control and Prevention; 2004.
- Cowan LS, Diem L, Brake MC, Crawford JT. Transfer of a Mycobacterium tuberculosis genotyping method, spoligotyping, from a reverse line-blot hybridization, membrane-based assay to the Luminex multianalyte profiling system. J Clin Microbiol. 2004;42:474–7.
- Cowan LS, Mosher L, Diem L, Massey JP, Crawford JT. Variable-number tandem repeat typing of Mycobacterium tuberculosis isolates with low copy numbers of IS6110 by using mycobacterial interspersed repetitive units. J Clin Microbiol. 2002;40:1592–602.
- van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B, Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993;31:406–9.
- Woelffer GB, Bradford WZ, Paz A, Small PM. A computer-assisted molecular epidemiologic approach to confronting the reemergence of tuberculosis. Am J Med Sci. 1996;311:17–22.
- Tarone RE. On heterogeneity tests based on efficient scores. Biometrika. 1985;72:91–5.
- Zelen M. The analysis of several 2 × 2 contingency tables. Biometrika. 1971;58:129–37.
- Qian L, Abe C, Lin TP, Yu MC, Cho SN, Wang S, rpoB genotypes of Mycobacterium tuberculosis Beijing family isolates from East Asian countries. J Clin Microbiol. 2002;40:1091–4.
- Marttila HJ, Mäkinen J, Marjamäki M, Ruutu P, Soini H. Molecular genetics of drug-resistant Mycobacterium tuberculosis isolates in Finland, 1995–2004. Int J Tuberc Lung Dis. 2008;12:338–43.
- Nitta AT, Knowles LS, Kim J, Lehnkering EL, Borenstein LA, Davidson PT, Limited transmission of multidrug-resistant tuberculosis despite a high proportion of infectious cases in Los Angeles County, California. Am J Respir Crit Care Med. 2002;165:812–7.
- Sun YJ, Lee AS, Wong SY, Heersma H, Kremer K, van Soolingen D, Genotype and phenotype relationships and transmission analysis of drug-resistant tuberculosis in Singapore. Int J Tuberc Lung Dis. 2007;11:436–42.
- Burgos M, DeRiemer K, Small PM, Hopewell PC, Daley CL. Effect of drug resistance on the generation of secondary cases of tuberculosis. J Infect Dis. 2003;188:1878–84.
- Vazquez-Gallardo R, Anibarro L, Fernandez-Villar A, Diaz-Cabanela D, Cruz-Ferro E, Perez del Molino ML, Multidrug-resistant tuberculosis in a low-incidence region shows a high rate of transmission. Int J Tuberc Lung Dis. 2007;11:429–35.
- Oeltmann JE, Varma JK, Ortega L, Liu Y, O’Rourke T, Cano M, Multidrug-resistant tuberculosis outbreak among US-bound Hmong refugees, Thailand, 2005. Emerg Infect Dis. 2008;14:1715–21.
- Comas I, Gagneux S. The past and future of tuberculosis research. PLoS Pathog. 2009;5:e1000600.
- Nicol MP, Wilkinson RJ. The clinical consequences of strain diversity in Mycobacterium tuberculosis. Trans R Soc Trop Med Hyg. 2008;102:955–65.
- Hazbon MH, Brimacombe M, Bobadilla del Valle M, Cavatore M, Guerrero MI, Varma-Basil M, Population genetics study of isoniazid resistance mutations and evolution of multidrug-resistant Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2006;50:2640–9.
- Afanas’ev MV, Ikryannikova LN, Il’ina EN, Sidorenko SV, Kuz’min AV, Larionova EE, Molecular characteristics of rifampicin- and isoniazid-resistant Mycobacterium tuberculosis isolates from the Russian Federation. J Antimicrob Chemother. 2007;59:1057–64.
- Hillemann D, Kubica T, Agzamova R, Venera B, Rusch-Gerdes S, Niemann S. Rifampicin and isoniazid resistance mutations in Mycobacterium tuberculosis strains isolated from patients in Kazakhstan. Int J Tuberc Lung Dis. 2005;9:1161–7.
- Gagneux S, Long CD, Small PM, Van T, Schoolnik GK, Bohannan BJ. The competitive cost of antibiotic resistance in Mycobacterium tuberculosis. Science. 2006;312:1944–6.
- Schwartzman K, Oxlade O, Barr RG, Grimard F, Acosta I, Baez J, Domestic returns from investment in the control of tuberculosis in other countries. N Engl J Med. 2005;353:1008–20.