Previous Article | Next Article ![]()
Journal of Clinical Microbiology, March 2004, p. 1016-1023, Vol. 42, No. 3
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.3.1016-1023.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Pathology and Microbiology, University of Nebraska Medical Center, Omaha, Nebraska 68198-6495,1 Division of AIDS, STD, and TB Laboratory Research, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia 303332
Received 8 April 2003/ Returned for modification 19 August 2003/ Accepted 5 October 2003
|
|
|---|
|
|
|---|
The genetic basis for PZA resistance involves mutation within the pncA gene, which encodes Pzase activity (20, 28). Although cases of PZA-resistant M. tuberculosis isolates with no pncA mutations have been reported, mutations of pncA and its putative promoter remain the major mechanism of PZA resistance (15, 20). Over 40 different mutations in either the pncA structural gene or its putative promoter associated with PZA resistance in M. tuberculosis have been described. The changes are either mutations that involve substitution of nucleotides or mutations in the form of nucleotide insertions or deletions (15, 20, 27). In contrast, the natural resistance to PZA demonstrated by M. bovis strains is uniformly due to a unique single-point mutation (C169G) in pncA. This mutation involves the replacement of histidine (CAC) with aspartic acid (GAC), leading to the production of inactive enzyme (26, 28).
Susceptibility testing to detect PZA resistance has recently received increased attention for a number of reasons, including the important role of PZA in shortening the time course for treatment of tuberculosis and the recent recognition of PZA-monoresistant strains of M. tuberculosis (9), the increasing frequency of tuberculosis infections following intravesical instillation of the naturally PZA-resistant M. bovis BCG strain for the treatment of superficial bladder cancer (1, 17, 19), and the increasing incidence of zoonotic tuberculosis in developing countries due to naturally PZA-resistant M. bovis (6, 16, 24).
Conventional mycobacterial susceptibility testing for PZA is dependent on growth of the organism in the presence of the drug. This technique is both time-consuming and potentially unreliable due to the poor growth of M. tuberculosis in the highly acidic medium required for PZA activity (7, 12). Automated testing systems, such as the BACTEC 460TB and BACTEC MGIT 960 systems, are more sensitive than conventional testing but require from 8 to 12 days to determine antibacterial susceptibility and have the potential for cross-contamination (12, 14, 31).
Genotypic assays for the detection of drug resistance have been applied to both cultured isolates and direct patient specimens. These include amplification techniques, DNA sequence analysis, PCR-single-strand conformation polymorphism electrophoresis, and structure-specific cleavage and DNA probe detection assays, all of which are capable of detecting mutations associated with drug resistance (8, 22, 30).
Temperature-mediated heteroduplex analysis (TMHA) using denaturing high-performance liquid chromatography (DHPLC) was originally applied to the detection of specific gene polymorphisms (21). The technology was recently applied to the detection of mutations associated with antituberculosis drug resistance (5). The technique utilized differential retention of homoduplex and heteroduplex DNAs under partial denaturing conditions for the identification of mutations in rpoB, katG, rspL, embB, and pncA that are responsible for rifampin, isoniazid, streptomycin, ethambutol, and PZA resistance, respectively. Additionally, a separate genetic element (oxyR) was utilized to differentiate between M. tuberculosis and M. bovis. Although the study demonstrated the feasibility of this approach for detecting resistance to multiple antimicrobial agents, detection of mutations in pncA was found to be problematic. The difficulty of detecting pncA mutations was attributed to the diverse natures of the mutations and their distribution throughout the gene and its putative promoter. It was proposed that the potential for highly stable DNA helices due to increased GC content within specific regions of the pncA gene represented a major technical challenge for TMHA methodology (5).
To overcome these difficulties, the analysis conditions of the TMHA assay were reengineered, and a second probe was added. In combination, these changes allowed the rapid identification of pncA mutations associated with PZA resistance and the ability to distinguish between the two closely related species of the complex, M. bovis and M. tuberculosis, using the same genetic target.
|
|
|---|
![]() View larger version (46K): [in a new window] |
FIG. 1. Alignment of pncA genes of wild-type M. tuberculosis and M. bovis and their putative promoters showing the positions of mutations in the 13 different mutant strains used in the study; mutant 1, G233A; mutant 2, C297G; mutant 3, del G71; mutant 4, A410G; mutant 5, T11C; mutant 6, T-7C; mutant 7, A29C; mutant 8, A139G; mutant 9, T398A; mutant 10, T515C; mutant 11, A152C; mutant 12, C185G; and mutant 13, C458A. *, unique mutation of M. bovis (C169G) that conveys natural PZA resistance. The arrow at the top indicates the start of the open reading frame.
|
Primer selection and PCR conditions. Specific primers were designed using Oligo version 6.4 software (Molecular Biology Insight, Inc., Cascade, Colo.) to generate a 638-bp amplicon that included the entire pncA gene and its putative promoter. The sequence of the forward primer, AW-A3 (5'-GTCATGGACCCTATATCTGTGGCTGCCGCGTCG-3'), began at bp -77 upstream of the open reading frame, and that of the reverse primer, AW-A6 (5'-TCAGGAGCTGCAAACCAACTCGACGCTGG-3'), began at the stop codon (bp 561). The PCR assay was performed using 5 µl of template DNA (10 ng/µl) in a total reaction volume of 50 µl, including PCR buffer (20 mM Tris-HCl [pH 8.4], 50 mM KCl); 0.1 mM (each) dATP, dGTP, dTTP, and dCTP; 1.5 mM MgCl2; 0.3 µM (each) primer; and 1.5 U of PlatinumTaq High-Fidelity DNA polymerase (Gibco BRL-Life Technologies, Gaithersburg, Md.). Amplification was performed on a Stratagene (La Jolla, Calif.) Robocycler model 96 thermocycler starting with an initial denaturation step at 95°C for 10 min, followed by 35 cycles, each cycle consisting of a denaturation step at 95°C for 1 min, an annealing step at 64°C for 1 min, and an extension step at 72°C for 1 min. An additional extension step at 72°C for 7 min was performed after the last cycle. The amplicons were stored at 4°C until they were used.
Cloning of PCR products and sequence analysis.
PCR products from selected PZA-resistant M. tuberculosis isolates were cloned directly following amplification using the standard protocol of the Original TA Cloning kit (Invitrogen, San Diego, Calif.). Purified plasmids from selected colonies were screened for the correct insert by digestion with the endonuclease EcoRI (New England Biolabs, Beverly, Mass.) and analyzed by gel electrophoresis for the presence of an
600-bp product. Selected plasmids were sequenced at the Epply Molecular Biology Core Laboratory (UNMC, Omaha) using the universal M13 forward and reverse sequencing primers. The sequences were analyzed for the presence of mutations of interest by alignment against the wild-type M. tuberculosis sequence using MacVector sequence analysis software version 6.5 (Oxford Molecular Group, Inc., Campbell, Calif.).
TMHA assay. The TMHA assay was performed using the commercial WAVE-DHPLC System (Transgenomics Inc. Omaha, Neb.). Since the hydrophobic matrix (polystyrene-divinylbenzene copolymer beads) of the WAVE-DNASep cartridge is electrostatically neutral and does not readily react with DNA, an ion-pairing reagent, triethylammonium acetate (buffer A), was used to adsorb DNA to the cartridge according to the manufacturer's protocol. An elution buffer composed of 0.1 M triethylammonium acetate in 25% acetonitrile (buffer B) was used to elute DNA based on size and/or sequence composition. Once eluted, the DNA was detected spectrophotometrically by UV absorption at 260 nm. The DNA molecules were analyzed for integrity by using nondenaturing conditions at a column temperature of 50°C, and they were analyzed for mutation detection by using partially denaturing conditions at a column temperature range of 52 to 70°C (21).
PCR product integrity analysis. The PCR products of all isolates were analyzed for purity, specificity, and DNA concentration using the universal DNA-sizing gradient concentration program and a column temperature of 50°C with DHPLC. The phiX174 DNA ladder was used as the sizing marker. The sizing capability of the WAVE system allowed analysis of purity, and only those amplicons shown to generate a single uniform peak of the correct size were used for subsequent analysis.
DNA hybridization. DNAs from the reference strains, M. tuberculosis H37Rv (ATCC 25618) and M. bovis ATCC 19210, were used for individual hybridization with each of the test isolates. In a total volume of 50 µl, equimolar ratios of test and reference DNA molecules were mixed together in the presence of polymerization inactivation buffer (5.0 mM EDTA, 60.0 mM NaCl, and 10.0 mM Tris [pH 8.0]). The mixture was heated to 95°C for 4 min and then left at room temperature for gradual cooling to 35°C over 45 min. For heteroduplex analysis, both homoduplex and heteroduplex molecules were generated by hybridization of the PCR product for each of the tested isolates with each of the reference DNA probes.
Heteroduplex analysis. Following hybridization, mixtures of test isolates and reference probes were analyzed for pncA mutations using the partially denatured mode of the DHPLC. A variety of gradient concentrations were examined with different starting concentrations of buffer B at different rates of increase (slope), and a range of column temperatures from 64.8 to 66.8°C was evaluated. A modified gradient concentration program (Table 1) and a column temperature of 65.8°C were chosen for all subsequent mutation detection studies. A set of three mixtures of wild-type reference DNAs (both M. tuberculosis and M. bovis) and reference probes was included with each run of the test isolates. Each of the test isolates was analyzed at least three times on three successive days using three different PCR products from each template to test the reproducibility of the chromatographic patterns. The chromatographic patterns of test isolates were compared with those of reference isolates, and interpretations were made according to the proposed protocol (Table 2). Accordingly, any test isolate which generated a single-peak pattern with the M. tuberculosis reference probe and a double-peak pattern with the M. bovis reference probe was identified as wild-type M. tuberculosis, whereas any test isolate which generated a double-peak pattern with the M. tuberculosis reference probe and a single-peak pattern with the M. bovis reference probe was identified as M. bovis or BCG. Isolates that produced a double-peak pattern with both reference probes were identified as mutant strains of M. tuberculosis (PZA resistant). A double-peak pattern was defined as a negative deflection following a peak that created a visible trough between adjacent peaks. For each of the double-peaked chromatographic patterns, the distance between the peaks was recorded.
|
View this table: [in a new window] |
TABLE 1. Modified gradient buffer concentrations for mutation detection within GC-rich pncA gene in comparison with universal gradient for mutation detection on WAVE system
|
|
View this table: [in a new window] |
TABLE 2. Proposed protocol for identification of test isolates using two different reference probes for hybridization in TMHA-DHPLC assay
|
|
|
|---|
Following optimization of the system, duplexes formed between PCR products of the tested isolates and each of the two reference probes were analyzed using the partially denatured mode of the system at the optimal buffer concentration gradient (Table 1) and column temperature (65.8°C).
Chromatographic patterns produced by the wild-type PZA-susceptible isolates of M. tuberculosis demonstrated single-peak patterns when mixed with the M. tuberculosis reference probe and double-peak patterns when mixed with the M. bovis reference probe, as predicted (Fig. 2A). In contrast, M. bovis isolates produced double-peak patterns when mixed with the M. tuberculosis reference probe and single-peak patterns when mixed with the M. bovis reference probe (Fig. 2B).
![]() View larger version (27K): [in a new window] |
FIG. 2. TMHA of pncA gene PCR products from reference control and test wild-type isolates using M. tuberculosis reference probe (A) and M. bovis reference probe (B). The chromatographic patterns a and b in each panel depict the wild-type reference control isolates of M. tuberculosis and M. bovis, respectively, with the reference probes. Chromatographic patterns 1, 3, and 5 are representative wild-type M. tuberculosis test isolates, and patterns 2, 4, and 6 are representative M. bovis test isolates.
|
![]() View larger version (24K): [in a new window] |
FIG. 3. TMHA of pncA gene PCR products from reference control and test mutant isolates using M. tuberculosis reference probe (panels A) and M. bovis reference probe (panels B). Chromatographic patterns a and b in each panel depict the wild-type reference control isolates of M. tuberculosis and M. bovis, respectively, with the reference probes. Chromatographic patterns 1 to 13 depict the 13 test mutant isolates with each of the reference probes. All mutant isolates demonstrated the predicted double-peak patterns with both probes, with the exception of mutant 3 and mutant 9 (circled).
|
![]() View larger version (19K): [in a new window] |
FIG. 4. (A) TMHA of pncA gene PCR products of mutant isolates 3 and 9 with the M. tuberculosis reference probe. The chromatographs show the difference in shape between the patterns obtained by mutant isolates 3 (Mut.3) and 9 (Mut.9) and that of wild-type M. tuberculosis (WT.). (B) TMHA of pncA gene PCR products of mutant isolates 3 and 9 with the M. bovis reference probe. The differences in retention time between the double-peak patterns of mutant isolates 3 and 9 and that of wild-type M. tuberculosis are illustrated.
|
|
|
|---|
In order to identify the optimal assay conditions, an extended range of column temperatures and various gradient concentrations were studied. This resulted in a modification of the universal gradient concentration recommended by the manufacturer for mutation detection. The modification process included shortening the run time from 18 to <10 min by starting the gradient at a higher elution buffer concentration (61 rather than 40). This change was made based on the predicted retention time of analyzed duplexes according to size. In addition, the slope of elution buffer during the run was reduced from 2 to 1.2% per min. The modification process also included evaluation of a range of column temperatures, starting from the 64.8°C recommended by the system software and ranging up to 66.8°C in 0.1°C increments. The optimal column temperature was determined to be 65.8°C, since all higher and lower temperatures failed to induce the production of the predicted chromatographic patterns. These modifications improved the correlation among the predicted chromatographic patterns based on the theoretical helical structure of heteroduplexes of GC-rich sequences and the observed patterns. The essential outcome of these changes was that the previously cryptic mutations within the GC-rich sequence of pncA could be revealed.
The observed chromatographic patterns following TMHA of the wild-type isolates of M. tuberculosis and M. bovis (Fig. 2) were consistent with the predicted patterns on which the study was based and allowed the differentiation of the two closely related members of the M. tuberculosis complex.
Given the diversity of pncA mutations that convey PZA resistance, it was important to test mutations from within all regions of the coding sequence, as well as the promoter element. To test the clinical applicability of our assay, 13 different PZA-resistant mutant strains of M. tuberculosis were evaluated. Eleven of these mutant isolates generated the predicted chromatographic pattern when mixed with both reference probes, i.e., a double-peak pattern with clear demonstration of an intervening trough between the peaks. Two mutant M. tuberculosis isolates (mutant 3 and mutant 9) did not produce the standard double-peak pattern when mixed with the M. tuberculosis reference probe. These patterns were found to be highly reproducible when mutant isolates 3 and 9 were tested repeatedly. Review of the sequence showed that mutant isolates 3 and 9 had mutations in two different regions of pncA with high GC content. This was consistent with the original suggestion by Cooksey et al. that the difficulty in detecting pncA mutations was due to the presence of GC-rich sequences adjacent to the mutated nucleotides (5). The influence of the GC-rich region on the chromatographic pattern generated by mutations within such sequences was subsequently confirmed by analyzing two additional mutant isolates within GC-rich regions (C401T and G511A). Under the same optimized conditions, these mutants produced patterns similar to those of mutant isolate 9 (data not shown). We concluded that single-point mutations within or near GC-rich regions of pncA were unable to disrupt the helical structure of the heteroduplex DNA under the given conditions, rendering them indistinguishable from the homoduplex DNA. It was subsequently demonstrated that mutations within GC-rich regions could be uncovered through an optimal combination of both column temperature and gradient buffer concentration. Although the entire M. tuberculosis genome is GC rich (65% overall), certain regions have higher GC contents than others, which is expected to affect the melting points of the DNA molecules within those regions. Depending on where the boundaries are drawn for the calculation, the melting temperature is considerably higher; for example, the region surrounding the mutation of isolate 3 has a GC content of 82%.
Production of chromatographic peaks using TMHA-DHPLC (WAVE) technology is a function of temperature and the interaction between the DNA duplex and the cartridge matrix under given buffer gradients. It has been reported that the DNASep cartridge, under nondenaturing conditions, resolves the DNA fragment independently of sequence composition (10). However, in our laboratory, shouldered peaks have been observed with certain GC-rich sequences, even under nondenaturing conditions (D. R. Bastola, P. C. Iwen, and S. H. Hinrichs, unpublished data). We demonstrated that specific sequences with predicted secondary structure generated by these GC-rich sequences were responsible for these shouldered peaks. At higher temperatures and under the optimal gradient concentration used in the present study, the chromatographic patterns generated from mutant isolate mixtures that contain both homoduplex and heteroduplex populations were expected to contain double peaks, or at least shouldered peaks, that were distinguishable from those of wild-type isolates that contain only homoduplex populations.
Another important difference between the chromatographs produced by mutant isolates 3 and 9 and those produced by wild-type M. tuberculosis isolates was apparent when both were analyzed with the M. bovis reference probe. Mutants 3 and 9 produced chromatographic patterns with two peaks that were separated by a greater distance than that of wild-type isolates (Fig. 4B). This increase in peak separation was also seen in all other mutant isolates when they were mixed with the M. bovis probe. The generation of widely separated peaks was a function of an earlier elution time for the heteroduplex formed by the mutant DNA in comparison with the heteroduplex formed by the wild-type M. tuberculosis DNA. One explanation for this observation is that the mutant heteroduplexes have greater secondary structure than the wild-type heteroduplexes. This is due to the presence of two base pair mismatches in the mutant heteroduplex, one in the mutant DNA and one in the M. bovis reference probe, compared to the wild-type heteroduplexes, which have only a single base pair mismatch that is present in the M. bovis reference probe. The greater secondary structure in the mutant-isolate heteroduplexes is believed to result in its elution earlier than the wild-type heteroduplexes.
When the observed patterns from both reference probes were considered together, mutants 3 and 9 could be distinguished from wild-type M. tuberculosis isolates, a characterization that could not be made if only one probe was utilized in the analysis.
Under optimal conditions, mixed DNA populations (homoduplex and heteroduplex) produce two separate double-peak patterns. The first double peak represents the two different heteroduplex DNA populations, and the second double peak represents the two different homoduplex populations. Mixed DNA populations can produce two, three, or four peaks, depending on the type and position of the studied isolate mutation. The reaction conditions used for this study were optimized to recover all types of mutations. In isolate 7, the run temperature, gradient concentration, and gradient slope resulted in three peaks, while in other isolates, the same conditions differentiate only the larger double peaks of homoduplex and heteroduplex combinations.
Although rare, processing of a mixed culture of wild-type M. tuberculosis and M. bovis may generate a double-peak pattern with both reference probes. This pattern could be confused with that of a PZA-resistant mutant M. tuberculosis isolate. However, differentiation between the two cases is still possible based on the distance between the double peaks generated by the unknown isolate when mixed with the M. bovis reference probe (Fig. 4B). The distance between the two peaks of the wild-type M. tuberculosis strain (in the mixed culture) with the M. bovis reference probe is shorter than that of a PZA-resistant mutant M. tuberculosis isolate.
Inclusion of interpretive software to automate the interpretation of chromatographic patterns, especially those that need careful analysis, is under development and is expected to greatly assist the use of this new technology within the clinical setting.
Demonstration of the specificity of the present assay was also important, since cross-contamination with nontuberculous Mycobacterium species is a well-known problem in other standard culture-based automated assays (14, 31). Specificity was achieved through the use of specific primers that selectively amplify the pncA target only from the M. tuberculosis complex and not from nontuberculous mycobacteria (data not shown). The simultaneous screening for PZA resistance and identification of M. tuberculosis complex members was generally accomplished within 24 h of obtaining an isolate. Since PCR can be applied to direct patient specimens such as bronchial wash fluid (30), even faster analysis is feasible and should be further investigated.
Detection of clinically important drug-resistant subpopulations among predominantly susceptible bacilli represents a challenge for most genotypic assays used for drug resistance screening. We are investigating the sensitivity of our assay for the detection of PZA-resistant subpopulations in heterogeneous primary clinical isolates.
The ability to detect mutations within GC-rich sequences, essential for identification of PZA resistance, and the simultaneous ability to distinguish between the closely related Mycobacterium species M. tuberculosis and M. bovis significantly expand the utility of TMHA-DHPLC methodology for clinical applications.
This work was supported in part by an equipment grant from Transgenomic Inc. Amr M. Mohamed was supported by a research fellowship from the Egyptian government.
|
|
|---|
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»