Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2005, p. 2210-2217, Vol. 49, No. 6
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.6.2210-2217.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Ishmael I. Rosas,1
Edward A. Graviss,3
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, Texas3
Received 13 December 2004/ Returned for modification 12 January 2005/ Accepted 8 February 2005
|
|
|---|
|
|
|---|
Molecular techniques (44) are receiving increased scrutiny as alternatives to traditional culture methods of drug susceptibility testing since they can directly detect drug resistance as DNA mutations (31, 47). Denaturing gradient gel electrophoresis (DGGE) is a simple yet powerful technique that is capable of detecting many mutations over large stretches of DNA (10, 44). With DGGE, mutations within an amplicon are recognized as alterations in melting temperature as the fragment migrates through a gradient of denaturants. GC-rich portions at the 5' end of the PCR primers serve as the highest melting domains and as clamps that prevent DNA fragments from denaturing completely. Mutational alterations in a fragment result in bands that migrate differently from the wild type. For organisms like M. tuberculosis, having a single gene copy, heteroduplexing of the DNA to a reference DNA (usually wild type) generates additional homoduplex and heteroduplex bands that facilitate mutation detection. DGGE can detect point mutations, insertions, and deletions and is more sensitive than DNA sequencing in identifying mutations in complex DNA samples (11, 40). DGGE has previously been used in the detection of rpoB mutations associated with RIF resistance in M. tuberculosis (20, 33), including in primary clinical samples from infected patients (33).
In this study, DGGE was used for the first time to detect PZA resistance mutations in pncA. Overlapping PCR products were used to scan the pncA open reading frame and flanking noncoding regions in a large set of M. tuberculosis isolates from the U.S.-Mexico border states of Texas and Tamaulipas, respectively. By using DGGE, DNA alterations were detected in 99% of PZA-resistant isolates, while only 1% of PZA-susceptible isolates harbored DNA alterations. Coupled with comparable data for detecting RIF resistance mutations in rpoB (20, 33), these results suggest that DGGE is a sensitive, fast, and inexpensive method for detecting drug resistance mutations in M. tuberculosis.
|
|
|---|
DGGE and heteroduplex analysis.
Primer sets used to scan for pncA mutations by DGGE are listed in Table 1. DGGE primers were designed using the SG Primer program (43). Primers contained 18 to 20 bp of unique M. tuberculosis sequence to amplify the region of interest. One primer in each pair also included a 40- to 50-bp GC clamp that functions as a highest melting domain. The GC clamps have a melting temperature of
95°C and prevent the strands of a PCR fragment from separating completely. The melting curve of pncA was generated using the Poland algorithm (http://www.biophys.uni-duesseldorf.de/local/POLAND/poland.html) as modified by Steger (37) to simulate the transition curve of double-stranded DNA in 19 mM NaCl (12).
|
View this table: [in a new window] |
TABLE 1. Sequences of primers used in this study
|
DNA sequencing of pncA. Initial sequencing of pncA alleles (see Fig. 2) was performed as described previously (36). Subsequent sequencing was performed as described elsewhere (20) using the SEQ primers described in Table 1. DNA sequencing was performed by Lark Technologies (Houston, TX), an FDA-approved sequencing facility. DNA changes relative to pncA from H37Rv (Rv2043c) were noted. Previous reports of PZA resistance mutations in pncA (3, 4, 6, 8, 9, 14-19, 22, 24-29, 31, 34-36, 38, 39, 46, 47) were used to identify unique alterations for this report.
![]() View larger version (71K): [in a new window] |
FIG. 2. Detection of pncA mutations by PCR fragments A through D. Isolates were scanned for mutations associated with PZA resistance in pncA using DGGE PCR products A, B, C, and D. Mutations were detected only within a single PCR product, and this allowed for the localization of mutations within the pncA open reading frame. Isolate (mutation): 1 (G17S, GGC AGC); 2 (I5-frameshift, ATC ATG TC...); 3 (Y103stop, TAC TAG); 4 (V93L, GTG CTG); 5 (T160TPE, ACA GCG ACA CCT GAA GCG); 6 (V180F, GTC TTC).
|
|
|
|---|
![]() View larger version (17K): [in a new window] |
FIG. 1. PCR products for DGGE and properties of pncA. A. Distribution of altered codons in pncA that were sequenced, including synonymous alterations. Codons that were identified more than once include identical mutations from independent strains as well as different alterations at the same codon. See Table 2 for details of mutations. The regions of pncA that are scanned by DGGE PCR products A through E are shown. B. Melting profile of pncA. The melting transition of the pncA open reading frame is simulated by the modified Poland application. The nucleotide length of pncA (561 bp) is aligned to the codon length of pncA (187 codons) in panel A above.
|
|
View this table: [in a new window] |
TABLE 2. Sequenced pncA mutations
|
|
View this table: [in a new window] |
TABLE 3. Summary of pncA mutation detection by DGGE and DNA sequencing
|
![]() View larger version (82K): [in a new window] |
FIG. 3. Detection of pncA mutations by fragments A and E. Isolates in lanes 1 to 11 were scanned for mutations with PCR fragments A and E as indicated. The presence of multiple bands on the denaturing gel indicates that DNA alterations were detected within the fragment. The wild type (H37Rv), which did not contain any mutations by definition, generated a single band (lane 11). Mutations in lanes 5 to 10 did not generate multiple bands within fragment A but did generate multiple bands within fragment E. Isolates in lanes 1 and 2 are not within the region scanned by fragment E. See Table 2 for details of mutations and altered codons. Mutations were detected only within one PCR fragment except where indicated. Lane (altered amino acid): 1 (L4S); 2 (I5 bp); 3 (D12A); 4 (G17S); 5 (L27P); 6 (L35P); 7 (Y41stop); 8 (T47A); 9 (D49A); 10 (P54L; also detected within fragment B, Fig. 4, lane 6); 11 (wild type, H37Rv).
|
![]() View larger version (82K): [in a new window] |
FIG. 4. Detection of pncA mutations by fragments B to D. Isolates in lanes 1 to 27 were scanned for mutations with PCR products B, C, and D as indicated. The presence of multiple bands on the denaturing gel indicates that DNA alterations were detected within the fragment. See Table 2 for details of mutations and altered codons. Mutations were detected only within one PCR fragment except where indicated. Lane(s) (altered amino acid): 1 (H71Y); 2 (S65S); 3 (H57D); 4 (H71E); 5 (G75G); 6 (P54L; also detected within fragment E); 7 (S67P); 8 (T76P); 9, 18, and 27 (wild type, H37Rv); 10 (Y103H); 11 (Y103stop); 12 (V93L); 13 (K96T); 14 (G97S); 15 (E91stop); 16 (IS6110 insertion at codon 114); 17 (V125F; also detected within fragment D, not shown); 19 (V180F); 20 (T160TPE); 21 (V139A); 22 (T160P); 23 (A134V); 24 (R154G); 25 (M175V); 26 (A134V).
|
TCT) at codon 65, which encodes a Ser residue. Another synonymous polymorphism at codon 75 (GGT
GGC), both encoding Gly residues, was observed in a PZA-resistant isolate that also harbored a D49A mutation. Interestingly, the latter alteration was not detected within fragment A even though it was within the region scanned by the PCR product. Other mutations in this same region were also not detected within the fragment (see below). One isolate contained an IS6110 insertion after codon 114. The insertion was apparent from agarose gels of the PCR-amplified products from this isolate, and it also resulted in an aberrant migration pattern on the denaturing gels (Fig. 4, lane 16). Two isolates harbored small insertions. One contained a 2-bp insertion after the first nucleotide of codon five and resulted in a frameshift mutation. The other was a 6-bp insertion after the second nucleotide of codon 160 that would result in the insertion of two amino acids (Pro and Glu) within the protein. Interestingly, another isolate contained a T160P mutation. The remaining isolates harbored single nucleotide changes. Three mutations resulted in nonsense mutations at codons 41, 91, and 103. Two strains with differing genetic backgrounds contained the same mutation (GCC
GTC), which encoded an A134V alteration, and the denaturation patterns from these isolates were identical (Fig. 4, lanes 23 and 26). One DNA sample yielded both wild-type and mutant sequence at codon 103; the mutant sequence encodes a Tyr-to-His alteration. As reported previously, DGGE can detect mixtures of DNA that might arise from mixed strains or emerging drug resistance (20, 33). Although these mixed cultures are usually confirmed by DNA sequencing, we have recently reported procedures to confirm the presence of DNA mixtures by DGGE (20). The remaining isolates harbored single nucleotide changes that are predicted to result in single amino acid changes within the open reading frame. Of the 33 DNA alterations within pncA reported here, 11 have not been reported previously (3, 4, 6, 8, 9, 14-19, 22, 24-29, 31, 34-36, 38, 39, 46, 47), although other alterations in the same codon have been previously reported for five of these novel alterations. Three of eighty-three isolates were monoresistant to PZA, and since this is a property of Mycobacterium bovis, which is naturally resistant to PZA, these isolates were investigated in more detail. The pncA gene of M. bovis (Mb2069c) encodes an Asp residue at codon 57, and it appears that this residue prevents the M. bovis PZase from activating PZA. M. tuberculosis pncA normally encodes a His residue at codon 57, and the encoded PZase can activate PZA and render the cells susceptible to the drug. An H57D mutation has also been identified in PZA-resistant M. tuberculosis isolates (31, 47, 48). Codon 57 is amplified by primer set B, and all three monoresistant isolates harbored pncA alterations within this PCR product. One isolate displayed a distinct denaturation pattern (Fig. 4, lane 2), and DNA sequencing revealed that its pncA encoded a P54L alteration (Table 2). The other two isolates displayed the same denaturation patterns (Fig. 4, lane 3), and DNA sequencing of pncA from one isolate revealed the H57D alteration. Strain genotyping indicated that both of these isolates had an IS6110 RFLP pattern (1.001) and spoligotype (676773777777600) characteristic of M. bovis. An investigation of the clinical laboratory records confirmed that these two isolates were M. bovis. Confirming earlier reports, PZA monoresistance alone is not a unique characteristic of M. bovis isolates (14). The DGGE assay is capable of detecting pncA polymorphisms that distinguish M. tuberculosis from M. bovis.
Resolving problem isolates. Problems were encountered in the initial analysis of 18 isolates whose PZA susceptibility did not match their DGGE pattern. These problems were resolved by retesting these isolates for PZA susceptibility, by DNA sequencing their pncA alleles, and by further DGGE analysis. Initially, 11 isolates were identified as PZA resistant but did not harbor DGGE-detectable mutations. Susceptibility was revised from resistant to susceptible for six of these isolates, and DNA sequencing confirmed that they lacked alterations in pncA. A single isolate was confirmed as resistant to PZA but lacked any alteration within pncA. Finally, four isolates were confirmed as resistant to PZA, and mutations within pncA were identified by sequencing. The mutations resided in a region of pncA amplified by primer set A and spanned codons 27 to 47 (Fig. 1, lanes 5 to 8 and 10). Only the nonsense mutation at codon 41 resulted in a partial denaturation pattern. An additional isolate harbored a defect at codon 49 that was also not detected (lane 9). This isolate also contained a synonymous alteration at codon 75 that was detected within fragment B. Inspection of the pncA gene revealed a domain of the gene containing a high GC content. While most of pncA has a melting temperature around 81 to 82°C, this region had a melting temperature peak close to 90°C (Fig. 1B). The undetected mutations were associated with this high melting domain. We hypothesized that these mutations were not detected because the high melting temperature of this region prevented the heteroduplexes from melting properly during electrophoresis. To address this problem, an additional primer set was designed. The fragment amplified by primer set E is slightly shorter than fragment A and contains a different combination of clamps to fine tune denaturation (Table 1). The electrophoretic conditions were also modified to detect denaturation products (see Materials and Methods). All of the mutations not detected in fragment A were detected in fragment E (Fig. 1, lanes 5 to 10). Once this was achieved, every sequenced pncA mutation was detected within one of the DGGE amplicons (Fig. 3 and 4). Only one PZA-resistant isolate, described above, lacked a mutation in pncA as assayed by DGGE and DNA sequencing.
Seven PZA-susceptible isolates were identified with DGGE mutations during the initial screening. One isolate remained PZA susceptible upon retesting and harbored the synonymous S65S polymorphism within pncA, as described above. PZA susceptibility results were revised from susceptible to resistant for three isolates, and DNA sequencing indicated that they harbored alterations in pncA indicative of drug resistance. These alterations included amino acid changes L4S, H71Y, and H71D (Table 2). The L4S mutation has been reported previously, and PZA-resistant isolates harboring this amino acid alteration lack PZase activity (26). The two amino acid changes at codon 71 have not been reported previously, but another alteration at this codon, H71R, has been reported in PZA-resistant isolates lacking PZase activity (26-28). Intuitively, isolates with these amino acid alterations have been defined as PZA resistant. Three other isolates remained PZA susceptible upon retesting. DNA sequencing indicated that these isolates lacked alterations in pncA, and repeat DGGE assays confirmed the absence of polymorphisms.
|
|
|---|
Although pncA was not sequenced from all 181 isolates, the concordance between DNA sequencing and DGGE was 100% for the 41 isolates from which pncA was sequenced. Of the 32 isolates with DGGE-detected DNA alterations, all 32 yielded pncA alterations upon sequencing. Furthermore, all of the codon changes in these isolates were localized to the correct region of pncA probed by the PCR fragments. In addition, all nine isolates without DGGE-detectable alterations were shown to lack pncA alterations by DNA sequencing. Similar concordance was observed between DGGE and DNA sequencing for RIF resistance mutations associated with rpoB (20). Altogether, we identified 32 distinct alterations in pncA, 11 of which have not been reported previously (3, 4, 6, 8, 9, 14-19, 22, 24-29, 31, 34-36, 38, 39, 46, 47). DNA alterations were distributed throughout pncA (Table 3), and representatives detected within the four major DGGE PCR products (A to D) were sequenced to establish a distribution of DNA alterations throughout the gene (Fig. 1). This strategy proved useful when testing whether the four primer sets were sufficient for detecting all pncA mutations in these isolates.
Initial analysis identified 18 isolates whose DGGE profiles did not match their PZA susceptibilities, including 11 PZA-resistant isolates without DGGE patterns and 7 PZA-susceptible isolates with DGGE patterns. The PZA susceptibility profiles of half of these isolates were reassigned upon retesting, and DNA sequencing confirmed that these changes were justified. Six isolates without DNA alterations in pncA were classified as PZA susceptible upon retesting, and three isolates with pncA polymorphisms were revised to PZA resistant. The altered residues of the latter isolates were characteristic of PZA resistance. One isolate remained PZA resistant, but defects in pncA could not be detected by DGGE or by sequencing. A small number of PZA-resistant isolates do not harbor mutations in pncA, and these isolates also usually retain PZase activity (31, 47, 48). PZase activity was not tested, and so this correlation remains to be confirmed for this isolate. Similarly, one PZA-susceptible isolate harbored a S67S synonymous polymorphism.
Problems with the initial DGGE assays were apparent for seven isolates. This included three PZA-susceptible isolates that were incorrectly assigned mutations based on the initial DGGE screen but lacked pncA alterations by sequencing or by repeat DGGE assays. In addition, four PZA-resistant isolates harbored pncA mutations that were not detected within the initial four DGGE PCR products. All four DNA alterations were located in a region scanned by fragment A; a fifth mutation was also identified in this region that was not detected. These mutations were clustered in a region of pncA that had an unusually high melting temperature due to a high GC content. The melting temperature of this region peaked around 90°C, which was about 10°C higher than that of the rest of the gene (Fig. 1B). Indeed, the GC clamps attached to the ends of the amplicons have a peak melting temperature of 95°C, which is only slightly higher than this naturally occurring M. tuberculosis domain. The GC clamps are attached to prevent complete separation of the DNA strands during electrophoresis. These results suggested that this region was not melting properly in the denaturing gradient gels, and therefore the mutations were not detected. Similar problems are expected for any molecular technique that probes for DNA alterations in pncA, probably including sequencing. For instance, Mohamed et al. (25) reported problems detecting a 1-bp deletion within codon 24 using denaturing high-performance liquid chromatography. Branch migration inhibition (18) and single-stranded conformation polymorphism (6, 34) have also been used to detect pncA mutations. These three techniques work on the same basic principles as DGGE in that they detect DNA mismatches that represent mutations. In theory, any molecular technique that relies on denaturation, hybridization, or DNA conformation (44) may face problems in detecting pncA alterations in this region. Once problems associated with this high-melting region were understood, an additional DGGE primer set was used to successfully detect mutations in this region. Mutations from five PZA-resistant isolates were localized to this region from a total of 83 PZA-resistant isolates scanned. Based on these results, it is estimated that 6% of PZA-resistant isolates harbor mutations in this region.
Mutations that inactivate the PZase encoded by pncA are distributed throughout the gene, and no major hot spots have been observed. Alterations, such as T47A and H57D, are commonly observed among groups of PZA-resistant isolates, but these mutations do not represent a significant proportion of the total resistant isolates. The pncA gene can be divided into eight regions that are scanned by the five DGGE PCR products, including four regions in which these fragments overlap. Based on the denaturation patterns that are characteristic of each mutation and on the location of a mutation to one of the aforementioned subdomains, it is possible to predict mutations by DGGE (e.g., Fig. 4, lanes 23 and 26). This property could be used in several ways. First, it is possible to detect the same mutation in different isolates based on their identical denaturation pattern. Second, it can be used to identify synonymous DNA alterations that are not associated with drug resistance. These are infrequent in pncA, although a few have been observed (3, 4, 6, 8, 9, 14-19, 22, 24-29, 31, 34-36, 38, 39, 46, 47). By recognizing these polymorphisms by their denaturation fingerprint, it should be possible to improve the specificity of a DGGE clinical assay to predict PZA susceptibility. Third, it could be used to track outbreaks and transmission of PZA-resistant strains. Fourth, polymorphisms can be used to type strains by multilocus sequence typing (41) or comparable techniques. Multilocus sequence typing can be performed using DGGE because of its high throughput and low cost (13). Finally, as shown here, DGGE fingerprints can also be used to help distinguish M. bovis isolates from other members of the M. tuberculosis complex (14).
The results obtained here are similar to those obtained when DGGE was used to probe rpoB for mutations associated with RIF resistance (20, 33). Using two DGGE PCR products, we were able to detect mutations in 99% of RIF-resistant isolates, which were also collected from the U.S.-Mexico border region (20). By contrast, only 2% of RIF-susceptible isolates produced detectable polymorphisms. RIF resistance mutations in rpoB are highly localized to an 81-bp stretch of rpoB between codons 507 and 533, but mutations outside this region exist as well. The dispersed nature of mutations in both genes makes it expensive to identify mutations by DNA sequencing (2, 7, 44). From the analyses of rpoB and pncA polymorphisms from our data sets, we have estimated that mutation detection by DNA sequencing is 15- to 50-fold more expensive than detection by DGGE. The lower range includes material costs only, while the upper range includes equipment, data analysis, and labor costs. Expense, sensitivity, and mutation detection present problems for other molecular assays since multiple probes (such as oligonucleotides) must be used. Although multiple primer sets are used for DGGE, they are capable of scanning large (50 to 500 bp) stretches of DNA and are capable of detecting all alterations in the amplicon.
Finally, DGGE assays revealed several instances of apparent errors in drug susceptibility testing. These included both false-positive and false-negative errors. These observations point out the inherent problems with drug susceptibility testing. This problem is more apparent with PZA since the pH of the culture assay must be around 6, which is near the limits for growth of mycobacteria. In our analysis of PZA and RIF resistance by DGGE, the number of apparent culture testing errors was comparable to the number of DGGE errors for both false-positive and false-negative scoring of susceptibility. For both RIF and PZA, it will be interesting to perform blinded studies comparing our DGGE assays to culture methods to predict drug susceptibility.
These experiments were initiated in the laboratory of Rebecca A. Cox, and we thank her for her support, encouragement, and enthusiasm. We are also grateful to Denise Dunbar for coordinating the pyrazinamide MIC testing and to Brent Calder for help in obtaining the Poland melting curve data.
Present address: Channing Laboratory, Brigham and Women's Hospital, Department of Medicine, Harvard Medical School, Boston, Massachusetts. ![]()
|
|
|---|
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»