Skip Navigation

Genetic Polymorphism in Mycobacterium tuberculosis Isolates from Patients with Chronic Multidrug-Resistant Tuberculosis

  1. Frank A. Post1,6,
  2. Paul A. Willcox2,
  3. Barun Mathema5,
  4. Lafras M. Steyn3,
  5. Karen Shean4,
  6. Srinivas V. Ramaswamy7,
  7. Edward A. Graviss7,
  8. Elena Shashkina5,
  9. Barry N. Kreiswirth5 and
  10. Gilla Kaplan6
  1. 1Infectious Diseases Unit,Cape Town, South Africa
  2. 2Departments of Medicine, University of Cape Town, Cape Town, South Africa
  3. 3Departments of Clinical Laboratory Sciences, University of Cape Town, Cape Town, South Africa
  4. 4Brooklyn Chest Hospital, Cape Town, South Africa
  5. 5TB Center, Public Health Research Institute, Newark, New Jersey
  6. 6Laboratory of Mycobacterial Immunity and Pathogenesis, Public Health Research Institute, Newark, New Jersey
  7. 7Department of Pathology, Baylor College of Medicine, Houston, Texas
  1. Reprints or correspondence: Dr. Gilla Kaplan, Laboratory of Mycobacterial Immunity and Pathogenesis, Public Health Research Institute, International Center for Public Health, 225 Warren St., Newark, NJ 07203-3535 (kaplan{at}phri.org).
  1. Presented in part: Keystone Symposia on Tuberculosis: Integrating Host and Pathogen Biology, Taos, New Mexico, 25-30 January 2003 (abstract 446).

Abstract

Multidrug-resistant tuberculosis (MDR-TB) is a major public health problem because treatment is complicated, cure rates are well below those for drug-susceptible tuberculosis (TB), and patients may remain infectious for months or years, despite receiving the best available therapy. To gain a better understanding of MDR-TB, we characterized serial isolates recovered from 13 human immunodeficiency virus-negative patients with MDRTB, by use of IS6110 restriction fragment-length polymorphism analysis, spacer oligonucleotide genotyping (i.e., “spoligotyping”), and sequencing of rpoB, katG, mabA-inhA (including promoter), pncA, embB, rpsL, rrs, and gyrA. For all 13 patients, chronic MDR-TB was caused by a single strain of Mycobacterium tuberculosis; 8 (62%) of the 13 strains identified as the cause of MDR-TB belonged to the W-Beijing family. The sputumderived isolates of 4 (31%) of the 13 patients had acquired additional drug-resistance mutations during the study. In these 4 patients, heterogeneous populations of bacilli with different resistance mutations, as well as mixtures of drug-susceptible and drug-resistant genotypes, were observed. This genetic heterogeneity may require treatment targeted at both drug-resistant and drug-susceptible phenotypes.

Multidrug-resistant tuberculosis (MDR-TB), which is caused by Mycobacterium tuberculosis isolates that are resistant to, at least, rifampin (RIF) and isoniazid (INH), is a serious public health hazard [1, 2]. Treating MDRTB can be difficult, because loss of use of the 2 most potent anti-TB drugs (i.e., INH and RIF) means that only less effective “first-line” and “second-line” therapy, which is more toxic and less efficacious than INH and RIF, is available [3]. Thus, although, for some patients, MDR-TB can be cured by short-course chemotherapy [46], for other patients, bacillary growth is merely suppressed as long as treatment is continued [7]. Furthermore, 8%–35% of patients with MDR-TB have persistently active disease that is refractory to multidrug therapy [711]. Consequently, in most studies, the cure rates for MDR-TB remain well below those for drugsusceptible TB, and mortality rates may be substantial, even among HIV-negative patients [7]. In addition, patients with MDR-TB that does not respond to treatment are a constant source of transmission of multidrugresistant M. tuberculosis [1216].

In contrast to most bacteria, for M. tuberculosis, acquisition of drug resistance does not occur as a result of horizontal transfer of resistance-bearing genetic elements. Rather, acquisition of drug resistance by M. tuberculosis results from mutations (caused by nucleotide substitutions, insertions, or deletions) in specific resistance-determining regions of the genetic targets (or their promoters) or activating enzymes of anti-TB chemotherapeutic agents [17]. In M. tuberculosis, resistance-conferring mutations occur spontaneously at a very low frequency. Inadequate therapy or subtherapeutic drug levels may provide a selective growth advantage and, thus, may favor the growth of a resistant phenotype that can ultimately predominate in persons in whom the disease was originally caused by drug-susceptible isolates [3]. Moreover, in patients with MDR-TB, selection for additional mutations may be accomplished by adding a single drug to a failing regimen [18].

In the human lung, selection of drug-resistance mutations in M. tuberculosis occurs predominantly within lung cavities for which high bacterial loads, active mycobacterial replication, and reduced exposure to host defense mechanisms have been reported [1921]. Because M. tuberculosis in sputum samples obtained from patients with TB originates from lung cavities, molecular analysis of serially recovered sputum isolates allowed us to study aspects of the genetic evolution of drug resistance in the human host.

PATIENTS AND METHODS

Patients. We studied M. tuberculosis isolates recovered from 13 HIV-negative patients with MDR-TB that was refractory to chemotherapy given for ⩾12 months. These 13 patients comprised the entire South African subset of a multinational, randomized, placebo-controlled study that evaluated the role of aerosolized recombinant human interferon-γ (rhuIFN-γ; Inter-Mune) as an adjunct to chemotherapy for patients with MDRTB [22]. Outcome analysis of this study of the use of rhuIFN-γ as an adjunct to antituberculous therapy is in progress and will be reported elsewhere. Consent was obtained from all patients studied, and the study was approved by the institutional review board of the University of Medicine and Dentistry of New Jersey (Newark) and the ethics committee of the University of Cape Town (Cape Town, South Africa).

All subjects were residents of the Western Cape region of South Africa, an area where TB is endemic (incidence of TB in 1999 was 546 cases/100,000 individuals). Because of a low prevalence of multidrug resistance in this area (0.3% in new patients and 1.7% in previously treated patients), all patients received treatment for drug-susceptible TB at diagnosis. Once our patients had MDR-TB diagnosed, drug regimens were individualized on the basis of the results of phenotypic drugsusceptibility tests. Drug susceptibility was evaluated every 2–3 months. For all patients, treatment regimens were adjusted, on the basis of the results of these evaluations, at ~6-month intervals. We performed a detailed microbiological analysis of the M. tuberculosis isolates recovered from these patients.

Sputum smears. Sputum smears obtained from all 13 patients remained positive for M. tuberculosis throughout the study. For each patient, the number of acid-fast bacilli present in sputum was determined, at entry into the interferon-γ (IFN-γ) study (i.e., at baseline), in 2 samples obtained separately (at least 2 weeks apart). By use of a standardized concentration procedure (incubation with 5% bleach, followed by high-speed centrifugation), smears of ⩾2 cm2 were stained with auramine and were examined using direct fluorescence microscopy with 400X magnification. Sputum smears were recorded as having <4, 4–40, or >40 bacilli per high-power field, and they were given a score of 1, 2, or 3, respectively. For each patient, the mycobacterial load in sputum was calculated as the mean of the loads of 2 sputum specimens that were obtained 1 week apart and before administration of rhuIFN-γ and/or placebo.

Radiography. Chest radiographs were obtained at baseline and were analyzed to determine the extent of TB. Each lung was divided into 6 zones (upper [apex to second rib], middle [second to fourth rib], and lower zones [fourth rib to diaphragm]). Each zone was given a score of 0–2 points, according to the degree of pulmonary involvement: 0 denoted no involvement, 1 denoted that <50% of the lung was diseased, and 2 denoted that >150% of the lung was diseased. Single lung cavities were given a score of 0.25, 0.5, or 1 point if they measured <2 cm, 2–4 cm, or >4 cm, respectively. Multiple lung cavities were given a score of 0.5, 1, or 2 points if they measured <2 cm, 2–4 cm, or >4 cm in diameter, respectively. The disease severity score, as determined by evaluation of chest radiographs, thus ranged from 0 to 24 points [23].

Sputum cultures and drug-susceptibility testing. Sputum samples were processed in an automated mycobacterial culture system (MGIT; Becton Dickinson). Cultures positive for growth were examined by fluorescence microscopy, for the presence of acid-fast bacilli, and by polymerase chain reaction, to confirm the presence of M. tuberculosis [24]. The phenotypic resistance of all isolates was determined at baseline. Resistance to RIF and INH was confirmed by growth of all isolates in liquid culture medium that contained 2 µg/mL RIF or 0.1 µg/mL INH, by use of the radiometric Bactec 460 system (Becton Dickinson). The susceptibility of M. tuberculosis to several other antituberculous agents was determined using the indirect proportion method, by use of oleic acid, albumin, dextrose, catalase (OADC)-enriched solid 7H10 medium that contained critical concentrations of 7.5 µg/mL ethambutol (EMB), 10 µg/mL streptomycin, 5 µg/mL kanamycin, 2 µg/mL ofloxacin, 10 µg/mL ethionamide, or 2 µg/mL thiacetazone. Phenotypic susceptibility testing for pyrazinamide (PZA) was not performed, because the results of this test can be difficult to reproduce and may not correlate well with drug susceptibility in vivo [24, 25].

IS6110 restriction fragment-length polymorphism (RFLP) analysis. DNA RFLP analysis was performed on every M. tuberculosis isolate [26]. In brief, Loewenstein-Jensen slants were inoculated with 0.1mL of MGIT broth and were incubated for 3–5 weeks. Mycobacterial DNA was extracted and was digested with PvuII restriction endonuclease. DNA fragments were separated by SDS-PAGE, were Southern blot transferred, and were hybridized with a labeled IS6110 probe. Bands were visualized by enhanced chemiluminescence and were analyzed with Whole Band Analyzer software (version 3.4; BioImage). Spacer oligonucleotide genotyping (i.e., “spoligotyping”) was performed by amplification of the direct repeat region of each isolate. Isolates were assigned a number according to the Centers for Disease Control and Prevention nomenclature [27]. M. tuberculosis isolates that belonged to theW-Beijing family [28] were identified on the basis of high IS6110 copy numbers, the presence of characteristic IS6110 insertions in the origin of replication region, and the presence of spoligotype S00034 [28, 29].

Genetic sequence analysis. Direct DNA sequencing of resistance-determining regions (rpoB, katG, mabA-inhA [including upstream regions], pncA, embB, rpsL, rrs, and gyrA) was performed on isolates obtained at baseline, during follow-up, and at the end of the study [17]. Whenever the sequence of ⩾1 gene had changed during the study, analysis of all isolates available was performed.

Statistical analysis. The 2-tailed Mann-Whitney U test was used to compare sputum bacterial loads and chest radiographic scores at baseline, and the 2-tailed t test was used to compare treatment durations. P<.05 was considered to be statistically significant.

RESULTS

Patients. Thirteen HIV-negative adults (8 men and 5 women; mean age, 35.8 years; age range, 22–60 years) were studied (table 1). All patients were sputum smear positive for MDR-TB and had received therapy for MDR-TB for a mean of 36 months (range, 12–77 months) at entry. Five patients had previously had TB, none of the patients had extrapulmonary TB, and 2 patients (patients 5 and 9) had diabetes mellitus; no other patients had specific risk factors for TB. None of the patients were related or had known epidemiologic links, and none were geographically clustered. Most patients excreted large numbers of bacilli in sputum (median score, 2.0) (table 1), and most patients had extensive disease noted on chest radiographs (median score, 11.0). Eight patients died during the 56 weeks of the study, most likely as a result of cachexia and/or chronic respiratory failure. Patients who died had more extensive disease ease (median score, 14.5), compared with patients who survived (median score, 9.0) (P=.012), despite these 2 groups of patients having had a similar mean duration of disease (41.7 vs. 47.2 months; P=.64) (table 1). A mean of 7.6 cultures/patient (range, 2–16 cultures/patient) was available for analysis.

Table 1.

Characteristics of patients with multidrug-resistant (MDR) tuberculosis (TB).

At the time that TB was originally diagnosed, all patients were treated with World Health Organization (WHO) category I therapy (i.e., treatment with INH, RIF, PZA, and EMB for 2 months, followed by treatment with either INH and RIF or INH, RIF, and PZA for an additional 4 months) for varying lengths of time (table 2). Once MDR-TB was diagnosed, the patients were switched to treatment regimens tailored to the phenotypic drug-susceptibility profile of their isolates (table 2). At entry to the study, therapy was again adjusted according to phenotypic drug susceptibility, treatment history, and the sideeffect profile. Eight of 13 patients received aerosolized rhuIFN-$#x03b3; as an adjunct to treatment (table 2).

Table 2.

Tuberculosis (TB) treatment received before and during the study.

Phenotypic and genotypic resistance profile of M. tuberculosis. The phenotypic and genotypic drug resistance of each patient is shown in table 3. All isolates displayed phenotypic resistance to INH. Eight isolates had single-nucleotide substitutions in katG, and 4 harbored the C15T mutation in the promoter region of mabA-inhA. In 1 INH-resistant strain (in patient 10), neither katG nor mabA-inhA mutations were identified (table 3). Eleven patients continued receiving INH, despite having documented resistance to this drug. Eight of these 11 patients received thiacetazone, which was only available in a combination that contained INH. For all isolates, resistance to RIF was confirmed by phenotypic and genotypic (rpoB) resistance testing. Phenotypic resistance to EMB was not observed in any of the 13 isolates. However, sequence analysis of embB, the genetic target of EMB, revealed resistance mutations in 62$#x0025; of isolates. The presence of aminoglycoside resistance mutations in rrs was poorly predicted by the results of phenotypic susceptibility tests. Taken together, these results suggest that the isolates from all 13 patients were highly resistant to many of the most potent first-and second-line agents.

Table 3.

Phenotypic and genotypic drug resistance in Mycobacterium tuberculosis isolates, as determined by cultures performed at baseline.

Genetic analysis of serially recovered sputum isolates. Despite living in an area with high rates of transmission of M. tuberculosis, the 13 patients who were studied had no evidence of coinfection with a second M. tuberculosis strain. This was demonstrated by the fact that all M. tuberculosis isolates recovered from the same patient showed identical IS6110 RFLP profiles (figure 1) and spoligotype patterns. It is of interest to note that 8 (62$#x0025;) of 13 isolates recovered from the patients with treatment-refractory MDR-TB belonged to either the WBeijing family (n=7) or the W-Beijing lineage (n=1) [29]. Of note, patients who were infected withW-Beijing isolates had higher sputum bacterial loads (median bacillary load, 2.5) than did patients who were infected with non-W-Beijing isolates (median bacillary load, 1.0) (P=.019).

Figure 1.

IS6110 restriction fragment-length polymorphism profiles of Mycobacterium tuberculosis isolated from sputum samples obtained from 3 patients at baseline (0) and after various lengths of treatment (6, 12, 17, 53, or 56 weeks). STD, molecular weight standard.

Sequence analysis of serially recovered M. tuberculosis sputum isolates. M. tuberculosis drug-resistance polymorphism was explored by sequencing the genes of serially recovered isolates (2–3 isolates/patient). For 9 patients, all drug targets that were examined remained unchanged, from baseline, during the 3–56-week study (table 3). However, for 4 patients (31$#x0025;), additional drug-resistance mutations had been selected during treatment (table 4). For patient 1, polymorphisms were noted in katG, embB, and gyrA. The gyrA gene mutation was found in some, but not all, bacilli at multiple time points. Patient 2 had polymorphisms in pncA, gyrA, and rrs; the pncA gene was mutated in some, but not all, bacilli at 2 time points. Patient 6 had polymorphisms in gyrA, pncA, rpoB, and katG, and the gyrA gene was mutated in some, but not all, bacilli at 1 time point. In patient 10, both gyrA and pncA had polymorphisms, and the pncA gene was mutated in some, but not all, bacilli isolated at 1 time point. For 3 of these patients, different nucleotide substitutions accounted for drug resistance in bacillary subpopulations (gyrA D94G, D94Y, and G94C accounted for resistance to fluoroquinolone in patient 1; gyrA A90V and D94Y accounted for resistance in patient 2; and rpoB D516V and S531W accounted for resistance to rifampin in patient 6; both the pncA V157G substitution and the 131GG insertion accounted for resistance to PZA in patient 6) (table 4). As seen in patients 1 and 6, isolates that harbored the wild-type (wt) katG gene (and, thus, potentially were susceptible to INH) were still recovered several years after the patient had developed MDR-TB (table 4).

Table 4.

Genotypic polymorphism of selected Mycobacterium tuberculosis genes in 4 patients.

DISCUSSION

Our studies show that the selection of drug resistance (a stochastic event) during treatment of TB results in a mixture of subpopulations of bacilli, some of which are drug resistant, whereas others are either drug susceptible or bear different genetic mutations that account for the resistance to the same drug (table 4). In the presence of anti-TB therapy, drug-resistant bacilli would be expected to outgrow drug-susceptible bacilli over time. The observed mixture of susceptible and resistant organisms suggests that the mutations may have been acquired recently. Polyclonal and heteroresistant populations of bacilli in sputum have been reported elsewhere [30-34]. In the present study, we demonstrate the presence of clonal discrete founder strains (on the basis of IS6110 RFLP analysis) that are composed of heterogeneous populations of bacilli that differ in their resistance profile. In a recent ex vivo study of multiple lesions obtained from excised human lungs, we noted that, for 3 of 3 cases, discrete founder strains of M. tuberculosis acquired resistance independently in separate physical locales within the same lung. This finding suggested the occurrence of parallel evolution—that is, the independent emergence of different bacillary genotypes within the same lung [21]. Sputum samples obtained from such a patient would likely contain bacillary subpopulations with a mixture of drug-resistant genotypes, similar to what was observed in the present study.

Phenotypic drug resistance is commonly determined by the proportion method, in which an M. tuberculosis isolate is considered to be resistant to a drug if ⩾1% of colonies grow on agar that contains a critical concentration of that drug [35]. Consequently, the particular drug can be withdrawn or withheld [3, 36]. Our observations suggest that it may not be prudent to withdraw the drug(s) after documentation of drug resistance. Rather, as is often practiced in the field, the existing drug(s) should be maintained to target the bacilli that are still drug susceptible. In addition, new drugs (preferably, ⩾2) should be added to the treatment regimen to target the newly drug-resistant organisms. For example, improved outcomes for patients with MDR-TB who were receiving therapy involving INH together with the other drugs have been reported [13]. In accordance with the findings of previous studies, acquisition of drug-resistance mutations was not associated with changes in the IS6110 profile [21, 37, 38]. During the course of our study, the IS6110 patterns of all 13 M. tuberculosis isolates did not change (7.2 person-years of observation), nor was any super infection detected. This was the case despite the high disease burden and rates of transmission of TB in the Western Cape and despite what appears to be constant exposure to TB in the environment [39]. In most published studies, evolutionary changes in IS6110 patterns occurred in <10% of patients, and it has been suggested that these events might be more frequent soon after infection, when replication is most active [37, 40, 41]. The small number of patients in our study and the long periods of active disease before enrollment may have contributed to the lack of any IS6110 instability in the present study.

Of note, 62% of cases of multidrug-resistant disease in our sample were caused by M. tuberculosis strains that belonged to the W-Beijing family of strains (spoligotype 34) (table 1) [28]. These strains have been reported to comprise ~20% of all M. tuberculosis isolates in the Western Cape [42]. The relatively high frequency of W-Beijing strains seen in our study may simply reflect a sampling bias rather than a significant difference between W-Beijing strains and other strains [43]. However, it may also reflect a specific predilection of this family of strains for developing drug resistance, as has been suggested elsewhere [44C46]. Alternatively, multidrug-resistant W-Beijing strains may have a comparatively higher propensity to spread in the community. Our observation of higher bacillary loads, at baseline, in the sputum of patients infected with W-Beijing strains may support this notion (table 1).

In agreement with previous studies, we observed a good correlation between phenotypic and genotypic susceptibility to INH and RIF [47-49]. The single patient in whom no katG or mabA-inhA mutations were identified may have had mutations in kasA, ahpC [17], or the recently identified INH target gene ndh, which encodes the enzyme NADH dehydrogenase, none of which was investigated in this study [50]. In contrast, the phenotypic and genotypic susceptibility to EMB correlated poorly [51, 52]. Together with nonavailability of susceptibility data for PZA (table 3), overestimation of the number of active drugs used to treat most of the patients in our study may have occurred [53].

The results of the present study, together with our previous analysis of multiple lesions in lung samples obtained from patients with chronic progressive TB [21], suggest that, in designing therapeutic strategies, we need to consider that the sputum used for diagnosis and for drug-susceptibility testing represents a static single sampling of a dynamic process. That is, the organisms found in a single sputum sample may not be representative of the entire population of organisms in the affected lung, as well as in other infected sites. New and improved diagnostic and drug-susceptibility testing methods that can detect heterogeneous subpopulations need to be developed to improve treatment of MDR-TB.

Acknowledgments

We thank Monica Chojnacki for her help in performing mycobacterial cultures, Wendy Luckey for assistance in obtaining sputum samples and chest radiographs, and Dorothy Fallows for help with the manuscript.

Footnotes

  • Financial support: Public Health Service (grants AI22616 and AI54361 to G.K.). F.A.P. was a Fogarty International Fellow (grant AITRP TW00231 to Columbia University [New York, NY]).

  • Received June 20, 2003.
  • Accepted January 5, 2004.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
  46. 46.
  47. 47.
  48. 48.
  49. 49.
  50. 50.
  51. 51.
  52. 52.
  53. 53.
| Table of Contents