Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2005, p. 2200-2209, Vol. 49, No. 6
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.6.2200-2209.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Edward A. Graviss,3
Barry N. Kreiswirth,4
Jan Vijg,2 and
Teresa N. Quitugua1
Department of Microbiology and Immunology,1 Department of Physiology, University of Texas Health Science Center at San Antonio, San Antonio, Texas,2 Houston Tuberculosis Initiative, Department of Pathology, Baylor College of Medicine, Houston, Texas,3 Public Health Research Institute, Newark, New Jersey4
Received 18 November 2004/ Returned for modification 22 December 2004/ Accepted 23 February 2005
|
|
|---|
|
|
|---|
![]() View larger version (11K): [in a new window] |
FIG. 1. Regions of rpoB scanned for RIF resistance mutations. The NI, NII, RRDR, CII, and CIII fragments that were amplified to scan for mutations in rpoB are shown. Codon numbers refer to alterations that were previously reported for RIF-resistant isolates (6, 28). Primer sequences and fragment properties are described in Table 1.
|
In this report, DGGE has been used to assay for RIF resistance mutations across 775 bp of rpoB. We screened a large set of M. tuberculosis isolates to build a database of rpoB DNA alterations and to identify a suitable combination of PCR target fragments that detect RIF resistance mutations. We detected mutations in 99% of RIF-resistant isolates within two DGGE target fragments. These results suggest that DGGE is a useful technique for detecting drug resistance mutations in M. tuberculosis.
|
|
|---|
DGGE and heteroduplex analysis.
Primer sets used for the scanning of rpoB mutations by DGGE are listed in Table 1. DGGE primers were designed using the SG Primer program (21). Primers containing 18 to 20 bp of M. tuberculosis sequence were designed to amplify the region of interest. One primer in each set also contained a 40- to 50-bp GC-rich portion (GC clamp) at its 5' end to serve as the highest melting domain (21). The GC clamps, with a melting temperature of
95°C, prevent the strands of a PCR fragment from separating completely. Electrophoretic mobility of the partially melted structure slows dramatically on the denaturing gel at its characteristic melting temperature, which is exquisitely sequence dependent (21, 23). A standard 50-µl PCR consisted of 20 ng of genomic DNA, 10 mM Tris-HCl, pH 8.3, 50 mM KCl, 1.1 mM MgCl2, 0.01% gelatin, 0.2 mM concentrations of each deoxynucleoside triphosphate, a 0.2 to 0.8 µM concentration of each primer, 2% dimethyl sulfoxide, and 2 U of REDTaq DNA polymerase (Sigma Chemicals, St. Louis, MO). Standard PCR amplification consisted of 1 min at 94°C, followed by eight cycles of 94°C for 1 min, 65°C for 1 min, and 72°C for 1 min, followed by 27 cycles of 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min, followed by 10 min at 72°C. Amplified products were analyzed on 2% agarose gels in TBE (100 mM Tris, pH 8.4, 90 mM boric acid, 1 mM EDTA).
|
View this table: [in a new window] |
TABLE 1. Primers used and regions scanned for rpoB mutations
|
DNA sequencing of rpoB. Initial sequencing of the isolates was accomplished using primers described previously (10). Subsequent sequencing was performed using the primers shown in Table 1. DNA fragments were amplified as described above, and the DNAs were purified on QiaQuick PCR columns (QIAGEN, Valencia, CA) and quantified by measuring the A260 on a SpectraMax 190 microtiter plate reader (Molecular Devices, Sunnyvale, CA). DNA sequencing was performed by Lark Technologies (Houston, TX), a Food and Drug Administration-approved sequencing facility. Sequencing primers were the same as the amplifying primers, except that only one primer was used in the sequencing reaction. DNA changes relative to rpoB from H37Rv (Rv0667) were noted. Codon numbers are based on Escherichia coli rpoB (accession no. AAC43085) after a BLASTP alignment (1) to Rv0667. Sequence analysis was performed using the Vector NTI version 8 suite of programs (InforMax, Madison, WI).
An initial set of isolates harboring previously sequenced mutations in rpoB (2, 10-11, 14) was used to establish a collection of distinctive mutations. This reference set of isolates harbored 16 distinct RRDR mutations (including double mutations) and four distinct mutations outside the RRDR (Table 2). rpoB mutations in the Texas-Tamaulipas border isolates were predicted based on their denaturation patterns relative to the fingerprints of previously sequenced mutations run in parallel. Samples with novel denaturation patterns were sequenced to identify the mutation, and each mutation was sequenced from at least one isolate.
|
View this table: [in a new window] |
TABLE 2. Mutations in M. tuberculosis rpoB gene detected by DGGE
|
|
|
|---|
The DGGE denaturation profiles of the mutations detected in the two sets of isolates are shown in Fig. 2. Each mutation generated a characteristic denaturation profile that was dependent on the nucleotide alterations and the position of these alterations within the DNA fragment. A mutation could often be predicted by its characteristic denaturation profile. This empirical approach was used to estimate the frequency of each mutation (Table 2). In addition, new mutations could often be predicted by their novel denaturation profiles, and this was used to select isolates for DNA sequencing. Finally, several mutations generated denaturation patterns that were very similar to one another (Fig. 3, top panel), indicating that it was sometimes difficult to unambiguously assign mutations based on the DGGE profile alone. To address this problem, additional DGGE assays were established.
![]() View larger version (55K): [in a new window] |
FIG. 2. DGGE profile of rpoB RRDR mutations. Thirty-four distinct mutations in the RRDR of rpoB were analyzed on 20 to 90% urea-formamide-8% acrylamide gels. The same isolate (an H526S S531A double mutant) was run in lanes 8, 16, 24, 32, and 40 as a denaturation standard. The following lanes have the given mutation(s): lane 1, 514 to 516 deletion; lane 2, 516 to 517 deletion; lane 3, 515F insertion; lane 4, L511R; lane 5, L511P; lane 6, D516V; lane 7, D516Y; lane 9, H526R (heteroresistant; mixed with wild type); lane 10, H526N; lane 11, H526D; lane 12, H526C; lane 13, H526Y; lane 14, H526S no. 1; lane 15, H526S no. 2; lane 17, F514F; lane 18, S522L; lane 19, H526R; lane 20, S531Y; lane 21, S531L; lane 22, S531W; lane 23, S531F; lane 25, Q513L; lane 26, Q513P; lane 27, R528R; lane 28, R528R and S531F; lane 29, Q510H and H526Y; lane 30, D516E and H526N; lane 31, N519K and R529K; lane 33, L511P and D516Y; lane 34, 513 to 517 H; lane 35, D516G; lane 37, D516E and H526N; lane 38, H526L; and lane 39, L533P. Numbered mutations (e.g., H526S no. 1) represent different nucleotide alterations (see Table 2 for details).
|
![]() View larger version (49K): [in a new window] |
FIG. 3. Mutation determination using different template DNAs for the heteroduplexing reaction. The rpoB RRDR was amplified from nine different isolates. The DNAs were heteroduplexed to either the wild type (top panel), an S531L mutant (middle panel), or an F514F mutant (bottom panel). Identity between the template and the isolate is indicated by a single band on the denaturing gel. The wild type was identified in lane 9 of the top panel, S531L mutations were identified in lanes 6 and 8 of the middle panel, and the F514F synonymous polymorphism was identified in lane 1 of the bottom panel.
|
The RRDR and CII region from 59 isolates were sequenced in order to identify mutations (Table 2). Isolates were chosen primarily by unique denaturation patterns, and this strategy was largely successful at identifying novel mutations. Other mutations were sequenced to confirm the DNA alteration. In addition, this region was sequenced in RIF-resistant isolates without DGGE-detectable mutations. DNA alterations were observed in only 3 of 143 RIF-susceptible isolates. One isolate harbored the synonymous F514F polymorphism, while two isolates harbored a synonymous R528R change. Another isolate harbored the R528R polymorphism in conjunction with an S531F alteration; however, this isolate was resistant to RIF. Eight RIF-resistant isolates lacked RRDR or CII mutations.
Detection of heteroresistance or mixed cultures by DGGE. DGGE detected five cultures that contained mixtures of DNAs within the RRDR (Fig. 4). While a homogeneous mutant DNA generates four bands on denaturing gradient gels after it has been heteroduplexed to a wild-type template (Fig. 4A, lane 6), the mixed samples (Fig. 4A, lanes 2 to 5 and 7) generated more than four bands, and the band intensities were variable from one to another. When DNA from a heteroresistant culture was run on a denaturing gel without the heteroduplexing reaction (Fig. 4B), the DNA still generated multiple bands. In contrast, wild-type and homogeneous mutant DNAs generated only a single band (Fig. 4B, lanes 1 and 6, respectively). This indicates that the heteroresistant cultures contain a mixture of DNAs. The mixed nature of these samples was confirmed. First, genotyping data revealed minor IS6110 bands in two cultures of isolates collected from human immunodeficiency virus-infected individuals, suggesting the presence of more than one strain. Second, DNA sequencing of the five samples revealed the presence of minor peaks at specific nucleotides. Interestingly, while two samples contained a mixture of wild-type and mutant DNAs, three samples contained a mixture of two different mutant DNAs. In some cases, the mutations were more consistently detected by DGGE than by DNA sequencing. Published reports indicate that DGGE is between 4- and 17-fold more sensitive than DNA sequencing in detecting mutations in mixed DNA samples (4, 19).
![]() View larger version (59K): [in a new window] |
FIG. 4. Heteroresistant isolates. (A) Heteroresistant samples display more than four bands of variable intensities when heteroduplexed to H37Rv. Heteroduplex of H37Rv to itself produces only one band (lane 1). In lane 6, an S531L mutation produces a characteristic four-band denaturation pattern. (B) PCR products of rpoB RRDR by DGGE without heteroduplexing reaction. Homogeneous DNAs in lane 1 (H37Rv) and lane 6 (S531L) produce a single denaturation band. Heteroresistant isolates in lanes 2 to 5 and 7 produce multiple bands, indicating that the PCR amplified DNA with a mixture of sequences. DNA mixtures were confirmed by DNA sequencing.
|
![]() View larger version (66K): [in a new window] |
FIG. 5. Detection of DNA alterations outside the RRDR. (A) Detection of DNA alterations in the NI and CIII regions of M. tuberculosis isolates. The NI and CIII target fragments were coamplified from isolates, and the heteroduplex products were resolved on the same denaturing gel. The following lanes show the given mutation(s): lanes 1 to 2 and 4 to 7, no mutation; lane 3, P646L in CIII; lane 8, V146F in NI; and lane 9, V146F in NI and V642A in CIII. (B) Typing CIII DNA alterations by heteroduplex analysis. CIII target fragments were amplified from nine isolates. After heteroduplexing to H37Rv (top panel), most of the isolates display a similar denaturation profile. However, after heteroduplexing to isolate 2 (bottom panel), five denaturation patterns can be distinguished. See Table 3 for complete genotypes and RIF phenotypes of these isolates. The following isolates have the given CIII alteration: isolate 1, V642A; isolate 2, P646S; isolate 3, P622Q; isolates 4 to 5 and 7 to 9, P646L; and isolate 6, E641E. (C) DNA alteration in NII and CII regions of rpoB by DGGE. The NII and CII fragments were amplified from isolates, and the heteroduplex products were resolved on denaturing gels. See Table 3 for complete genotypes and RIF phenotypes of these isolates. The following lanes show the given mutation: lane 1, G438G in NII; lane 2, Q482R in NII; lane 3, P564S in CII; lane 4, E562L in CII; lane 5, I572F in CII; and lane 6, I572L in CII.
|
|
View this table: [in a new window] |
TABLE 3. Genotypes at rpoB of isolates with DNA alterations outside of the RRDRa
|
AGC) (Table 3). To our knowledge, this is the first report of V146F and RRDR double mutations. The V146F mutation was not observed in any of the other RIF-resistant isolates with previously identified RRDR mutations. Finally, no DNA alterations in the NI region were detected in any of the 105 RIF-susceptible isolates. DNA alterations were detected within the CIII target fragment in 9 of the 210 isolates (Tables 2 and 3). A V642A mutation was detected in one RIF-resistant isolate that also harbored a V146F mutation, as described above. A P622Q alteration was observed in one RIF-susceptible isolate. The other seven CIII mutations produced similar denaturation patterns when heteroduplexed to H37Rv (Fig. 5B, top panel). However, further analysis revealed that these seven isolates harbored three different mutations. When the CIII fragment from one of the seven isolates was used as the heteroduplex template, three distinct patterns were observed (Fig. 5B, bottom panel). DNA sequence analysis confirmed that these three DGGE patterns corresponded to three distinct DNA alterations in the CIII region. One RIF-resistant isolate harbored an E641E synonymous polymorphism; this isolate also harbored an H526D mutation in the RRDR. A P646S alteration was observed in a RIF-susceptible isolate. Finally five isolates harbored a P646L alteration. Two were RIF resistant and also harbored RRDR changes (D516V and H526Y). The remaining three isolates were RIF susceptible. These CIII DNA alterations do not appear to have been reported previously (7, 20, 28).
Three RIF-resistant isolates harbored alterations in the CII region (Fig. 5C, lanes 3 to 6). Two isolates harbored I572L mutations. These isolates also contained RRDR mutations (D516G and D516Y). Another isolate harbored an I572F alteration. Interestingly, this isolate did not harbor a defect in the RRDR. Similarly, a reference isolate harboring the I572F defect lacked RRDR alterations. Isolates harboring this mutation were previously reported by three other groups (3, 6, 27); none reported an additional RRDR mutation. Reference isolates harbored mutations in both the RRDR and the CII region. One isolate harbored an L511P and E562L combination, while the other contained S531L and P564S alterations. Finally, CII DNA alterations were not detected from any of the RIF-susceptible isolates.
A single alteration was detected in the NII region in the 210 isolates scanned. A G438G synonymous polymorphism at codon 438 was detected in one RIF-susceptible isolate (Fig. 5C). In addition, a Q482R mutation was detected in a reference sample that also harbored an S531L defect.
Revised RIF susceptibility based on DGGE assays. Thirty-five isolates were retested for RIF susceptibility. Half of these isolates harbored common rpoB mutations that were either previously reported to impart resistance to RIF (7, 17, 20, 28) or synonymous RRDR polymorphisms. The remaining 18 isolates lacked identified RRDR mutations by DGGE. The isolates were tested by the BACTEC460 method at three different RIF concentrations: 0.5, 2, and 8 µg/ml. Isolates with either V146F or I572F mutations grew in 8 µg/ml, confirming that these alterations contribute to the drug resistance phenotype. However, the susceptibility results of 10 isolates were revised upon retesting. One isolate that contained a mixture of DNA in the rpoB RRDR grew in 2.0 µg/ml RIF, one isolate with a mixture of DNA grew in 8 µg/ml RIF, while a third isolate with mixed DNA was unable to grow in 0.5 µg/ml RIF and was reclassified as susceptible. Finally, seven of eight isolates without detectable mutations in the amino-terminal region, RRDR, CII region, or CIII region were reclassified as RIF susceptible. An examination of the clinical laboratory records indicated variable and inconsistent results regarding RIF and rifabutin resistance for most of these isolates. The revised susceptibility was justified, since it was consistent with the absence of DNA alterations in rpoB, which was predicted by DGGE and confirmed by DNA sequencing. Once the mutation and susceptibility data were reconciled, only 1 of 153 RIF-resistant isolates lacked an identified rpoB DNA alteration.
|
|
|---|
The most useful DGGE primer set scanned for mutations in the RRDR of rpoB, and 94.8% of the RIF-resistant isolates contained mutations in this fragment (Table 4). This was similar to the accumulated reported global summation derived primarily from DNA sequencing (7, 17, 20, 28). Every RIF-resistant RRDR mutation was detected from the two sets of isolates used in this report. This can be stated with confidence, since this region has been sequenced from every RIF-resistant isolate lacking a DGGE-detectable RRDR mutation, and no DNA alterations were missed. This confirms the sensitivity of the DGGE assay in detecting mutations (4, 19). Using DGGE in combination with DNA sequencing, we have identified 34 different DNA alterations within the RRDR. These mutations are representative of the RRDR defects reported globally (7, 17, 20, 28), although several double mutations appear to be novel.
|
View this table: [in a new window] |
TABLE 4. DGGE detection of rpoB DNA alterations in M. tuberculosis isolates
|
A third primer set scanned for mutations in the CII region of rpoB, and mutations were detected in three RIF-resistant isolates. Two harbored an I572L alteration in conjunction with different alterations at codon 516 in the RRDR, D516G or D516Y. The third isolate harbored an I572F substitution but did not harbor a mutation in the RRDR. Identification of alterations at codon 572 contributed in a small way to the successful detection of rpoB mutations associated with RIF resistance. While the I572F mutation alone may be sufficient for drug resistance, the I572L alteration has been detected only in the presence of RRDR mutations. Further testing will be required to determine if the latter mutation directly contributes to RIF resistance.
Another DGGE target fragment was used to identify five distinct DNA alterations within the CIII region between codons 622 and 646. Alterations in this region were previously reported for RIF-resistant isolates (6, 28). However, only one of the CIII polymorphisms identified in our study, V642A, was uniquely associated with resistance to RIF. However, it is unclear whether the V642A mutation contributes to drug resistance, since it was observed in an isolate that also contained a V146F mutation. The other CIII polymorphisms were either synonymous changes or were not observed exclusively in resistant isolates. Similarly, a single synonymous alteration at codon 438 in a RIF-susceptible isolate was detected in the NII region. Since the DNA alterations detected within these two regions do not appear to be uniquely associated with drug resistance, these rpoB alterations may have limited usefulness in predicting drug resistance. However, the NII and CIII target fragments may be useful in detecting rpoB DNA alterations in other populations of M. tuberculosis strains with a different distribution of mutations.
Since each mutation generated a distinct and reproducible denaturation profile, we tested whether this property could be used to predict mutations. The successful application of this approach was dependent on three main factors: the DGGE target fragment, the number of different mutations detected within the target fragment, and the heteroduplex template. For the CII region, only four different mutations were observed, and these could be readily distinguished by denaturation profile (Fig. 5C). On the other hand, three of the five CIII DNA alterations generated similar denaturation profiles that could not readily be resolved with a single heteroduplex template (Fig. 5B). Prediction was even more complicated when comparing RRDR mutant profiles from 34 different alleles (Fig. 2 and 3). For both the CIII and RRDR target fragments, mutation prediction was greatly facilitated by using mutant DNA as a heteroduplex template. In both Fig. 3 and Fig. 5B, mutations that displayed similar denaturation profiles with H37Rv as a heteroduplex template could be distinguished when mutant heteroduplex templates were used. Since there is a greater degree of mismatches between two mutant DNAs, they can be expected to have a larger difference in melting temperature than mutant-wild-type heteroduplex templates. It was also possible to define mutations that were identical to the mutant heteroduplex template, since the homologous mutations generated a single band in the gels. Further resolution can be achieved with additional mutant templates (Fig. 3). It is expected that this type of analysis could be utilized to assign mutations using a limited number (e.g., two to four) of heteroduplex templates. Identifying the DNA alteration within the RRDR by denaturation profile is of interest, since it would help identify polymorphisms that are not associated with drug resistance, such as the F514F and R528R polymorphisms. This could be used to increase the specificity of a DGGE clinical assay to predict RIF susceptibility.
Using DGGE, we have detected five separate cultures that contain mixtures of DNAs. These heteroresistant cultures may have arisen in several ways (5). First, they may represent the emergence of drug-resistant cells within a strain (9, 15). Second, they may represent multiple infections of an individual by different M. tuberculosis strains (24). Third, they may represent laboratory contamination of one strain by another (8). Fourth, they may represent a heterogeneous population of bacilli with different resistance mutations (15). Regardless of their origin, it is difficult to detect these types of cultures. In this respect, DGGE proved especially useful and, as we have demonstrated, is capable of providing additional information. We observed instances of wild-type and resistant rpoB as well as mixtures of mutations within the RRDR, although the exact causes of these specific mixed cultures have not been determined.
Only 1 RIF-resistant isolate out of 153 scanned did not have a detectable rpoB DNA alteration as analyzed by DGGE and confirmed by DNA sequencing. By scanning in the RRDR and NI regions, mutations in 98.7% of RIF-resistant isolates could be detected (Table 4), and by adding the CII fragment, 99.3% of the resistant isolates had identifiable mutations. Specificity is also high, since only 2.1% of RIF-susceptible isolates harbored DNA sequence variants that might be interpreted as resistant, producing a "false positive" interpretation. All of the DNA alterations in RIF-susceptible isolates were detected in the RRDR. By learning how to recognize sequence variants in the RRDR, it should be possible to improve the specificity of the DGGE assay. Equally impressive is the ability of DGGE to detect inconsistencies in drug susceptibility testing. Seven isolates without rpoB mutations were revised to RIF susceptible upon reanalysis. This resulted in a 5 to 6% improvement in the sensitivity of DGGE assay.
We have reported similar results using DGGE to probe for pyrazinamide resistance mutations in the pncA gene of M. tuberculosis isolates (12). This suggests that DGGE can be used to detect mutations that impart resistance to other antituberculosis drugs. Scarpellini et al. (18) also analyzed RIF-resistant isolates using DGGE, and they were able to detect 14 different mutations in 81 RIF-resistant M. tuberculosis isolates. They were able to detect RRDR mutations in clinical samples, such as sputum, suggesting that the direct diagnosis of drug susceptibility may be possible using a DGGE assay. The RRDR and NI fragments can be coamplified and can be resolved on the same denaturing gel. Many other steps in the DGGE assay lend themselves to high-throughput screening. It will be interesting to test whether the RRDR and NI target fragments can be used to accurately predict RIF susceptibility in blinded studies, where the sensitivity and specificity of the DGGE assay can be compared directly to more traditional culture assays.
These experiments were initiated in the laboratory of Rebecca A. Cox, and we thank her for her support, encouragement, and enthusiasm. We are grateful to Terry Goen for performing the rifampin MIC testing.
Present address: Channing Laboratory, Brigham and Women's Hospital, Department of Medicine, Harvard Medical School, Boston, Mass. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»