Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, June 2009, p. 2679-2683, Vol. 53, No. 6
0066-4804/09/$08.00+0 doi:10.1128/AAC.00063-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Menzies School of Health Research, Institute of Advanced Studies, Charles Darwin University, Darwin, Northern Territory, Australia
Received 15 January 2009/ Returned for modification 15 February 2009/ Accepted 27 March 2009
|
|
|---|
|
|
|---|
The increasing prevalence of bacteria producing extended-spectrum β-lactamases (ESBLs) has become a significant problem facing health care across the world. Many ESBLs are derived from plasmid-encoded SHV or TEM β-lactamases (12-14). The conversion of a non-ESBL SHV enzyme into an ESBL is nearly always associated with a G238S (GGC
AGC) substitution (numbering according to that of Ambler et al. [1]), while a further extension of the spectrum of activity is mediated by an E240K (GAG
AAG) substitution (12). Although there are many other substitutions reported, mutations of codons 238 and 240 are by far the most significant for ESBL activity of the blaSHV family (http://www.lahey.org/Studies/).
An unusual aspect of the biology of the SHV-encoding gene, blaSHV, is that it is present on the chromosome of most Klebsiella pneumoniae strains (2, 3). However, it is also disseminated on plasmids, and most SHV ESBLs are plasmid encoded. It appears that there have been two recent mobilizations of non-ESBL blaSHV from the K. pneumoniae chromosome onto plasmids, the mobilized genes being blaSHV-1 and blaSHV-11, which differ only at codon 35 (7). The G238S mutant-encoding derivatives of these are blaSHV-2 and blaSHV-2a, respectively, while the G238S and E240K double-mutant-encoding derivatives are blaSHV-5 and blaSHV-12, respectively. These six blaSHV variants are all abundant.
The existence of both chromosome and plasmid-borne blaSHV means that many K. pneumoniae strains harbor mixtures of blaSHV alleles, and the ratios between alleles can vary widely (6, 8, 9). It has also been determined that the presence of a minority ESBL-encoding allele confers an ESBL-positive phenotype (9). This provides a diagnostic challenge.
An allele-specific real-time PCR (sometimes known as kinetic PCR) method for interrogating the blaSHV codon 238 and 240 SNPs has previously been reported (8, 9). This has proven effective, but it requires four separate PCRs. Also, it is difficult to calibrate for interlab comparison because of the effects of small batch-to-batch variations in primer concentration and quality (unpublished data). In contrast, the allele discrimination in an HRM-based assay takes place at the end of the PCR, so the format is inherently more robust.
Here, we demonstrate a single-tube HRM-based assay for codon 238 and 240 mutations in the blaSHV gene and its comparison with allele-specific real-time PCR. To support this assay, data analysis methods were developed for the inference of confidence limits regarding whether or not an analyte contains a mutant allele. This procedure is similar in some respects to the probe-based melting analyses of Chia et al. and Randegger and Hachler (4, 15), but it is inherently simpler and provides more information.
Two isolate collections of K. pneumoniae were used in this study (Table 1). Isolates A1 to L1 are clinical isolates from the Princess Alexandra Hospital (PAH) in Brisbane, Australia, have previously been characterized, and contain various mixtures of blaSHV-11, blaSHV-2a, and blaSHV-12 (8, 9, 11, 16). Isolates 14 to 121 are derived from the SENTRY collection (10). All strains were cultured in Luria-Bertani (LB) broth and stored in cryovials with 12% glycerol at –80°C. The isolates that are asterisked in Table 1 are derivatives of the clinical isolates that have been subjected to selection for resistance to 128 µg/ml cefotaxime (8).
|
View this table: [in a new window] |
TABLE 1. blaSHV allele distributions and results from allele-specific real-time PCR and HRM assays
|
Allele-specific PCR was carried out essentially as previously described by Hammond et al. (9). The reactions were performed in 10-µl volumes, containing SensiMix NoRef master mix at a 1x concentration (Quantace), Sybr green I at a 1x concentration, 5 pmol of common primer (Shv238reverse [5'-CGGCGTATCCCGCAGATAA-3'] or Shv240reverse [5'-CCGGCGGGCTGGTTTAT-3']), 5pmol of either the wild-type (WT)-specific primer (Shv238wt [5'-CGCCGATAAGACCGGAGCTG-3'] or Shv240wt [5'-GCGCGCACCCCGCTC-3']) or mutant-specific primer (Shv238mt [5'-CGCCGATAAGACCGGAGCTA-3'] or Shv240mt [5'-GCGCGCACCCCGCTT-3']), and 50 ng of DNA template. Thermocycling parameters were as described by Hammond et al. (9). The
CT values were calculated by subtracting the CT for the WT-specific reaction from the CT for the mutant-specific reaction. It has previously been shown by Hammond and coworkers that ESBL-negative isolates have codon 238
CT values of >2.5 and codon 240
CT values of >5 (8), and these are the cutoffs used for presumptive mutant allele detection by allele-specific real-time PCR.
The HRM assay was carried out using primers SHVmutHRM_F (5'-CGCCGATAAGACCGGAGCT-3') and SHVmutHRM_R (5'-CCGCGCGCACCCCGCT-3'). These were designed to anneal very close to the SNPs, and the amplified fragment was only 39 bp. The small size of this fragment maximizes the differences in melting temperature (Tm) conferred by the SNPs and allows for short cycling times. PCR was performed in a 10-µl reaction, containing SensiMix NoRef master mix at a 1x concentration (Quantace), Sybr green I at 1x concentration, 5 pmol of each primer, Q-solution at a 1x concentration (Qiagen), and 1 ng of DNA template. Thermocycling parameters were as follows: 50°C for 2 min; 95°C for 10 min; 40 cycles of 95°C for 5 s, 65°C for 5 s, and 72°C for 10 s; 95°C for 2 min; and 50°C for 30s. The amplification was followed by HRM from 70 to 86°C at 0.05°C increments, remaining at each step for 2 s. In order to calculate confidence limits for conclusions regarding the sequence of the amplified fragment, the numerical data defining the normalized HRM curves were exported using the "export" function of the software supplied with the Corbett 6000 device, and the mean and standard deviation (SD) were defined at each temperature of the melting protocol. From this, the 95% confidence limits were calculated as the mean ± 1.96 x SD. These numbers were used to define a 95% confidence limit area on the difference graph that is obtained by setting the mean of all of the WT HRM data as the baseline. The amplitudes at the nadirs of the HRM difference graphs were calculated with reference to the mean for all WT isolates.
The fully characterized samples for the development of this assay were K. pneumoniae clinical isolates with blaSHV allele ratios that had previously been determined by allele-specific real-time PCR and also by the cloning and analysis of PCR products (5, 10). These were subjected to HRM analysis, and in addition, the allele-specific real-time PCR assays were repeated so as to control for the effect on primer extension efficiency of different primer batches (Table 1 and Fig. 1A and B). The ESBL-positive isolates were clearly discriminated from the ESBL-negative isolates, with the mutant samples having lower Tm values, consistent with the G
A substitutions. The ESBL-positive isolates included J2, which has 80% WT alleles and 20% double-mutant (blaSHV-12) alleles. This yielded a very different HRM curve from the WT isolates. Therefore, this method can detect mutant alleles when they in the minority. The clear identification of a minority allele probably reflects the formation of heteroduplexes which depress the Tm. In the difference graph (Fig. 1B), it can be seen that the double-mutant (blaSHV-12)-containing isolates J2, L1, and I1 can be easily discriminated by eye because the curve nadirs are displaced to the left. However, isolate D1, which contains just 10% double-mutant (blaSHV-12) alleles and a high proportion of single-mutant (blaSHV-2a) alleles, was not discriminated from A1, in which only the WT and single-mutant alleles have been detected. In addition, there was correlation between the difference graph amplitudes and the ratio between the WT and mutant alleles. Because the difference graph amplitudes are derived from melt curves that are always normalized to the same number of arbitrary fluorescence units by the Corbett 6000 software, they can be compared between runs.
![]() View larger version (21K): [in a new window] |
FIG. 1. HRM data. (A) Normalized HRM graphs from the PAH isolates. Isolates containing mutated alleles (A1, D1, J2, L1, and I1) are clearly separated from isolates containing only WT alleles (B2, J3, and K2). (B) Difference graph from the PAH isolates, with the fluorescence of isolate J3 set as the baseline. For each isolate, the number and allelic distribution are stated. (C) Difference graph of the additional 27 SENTRY isolates (18 WT, 7 single mutant, 2 double mutant) isolates. The baseline used was the average fluorescence data of 21 WT (ESBL negative) isolates. The area representing the 95% confidence limit for the WT sequence is delineated with dotted lines. All ESBL-positive isolates yielded difference graph curves outside the 95% confidence limit area. The two curves with nadirs displaced to the left are the two double-mutant isolates, 14 and 18.
|
CT values versus the amplitude at the nadirs of the HRM difference graphs for the "single-mutant" ESBL-positive isolates is 0.94 (P < 0.0001). The value for "double-mutant" isolates is 0.98 (P = 0.0023). This confirms that this assay can provide an indication of allele ratios. The data from the 18 SENTRY WT isolates together with the 3 PAH WT isolates (B2, J3, and K2) were used to obtain 95% confidence limits for the HRM curve corresponding to the WT allele. This not only indicates when the presence of a mutation should be called but also facilitates the portability of this method. The Tm of the WT sequence using our Rotor-Gene 6000 device is 79.00 ± 0.15°C. It is our experience that an individual Corbett 6000 device is highly accurate with respect to relative temperatures during runs, but absolute temperature calibration between different Corbett 6000 devices can differ by up to 0.5°C. The 95% confidence limit curves for the WT sequence can easily be adjusted with reference to a Tm determined from a small number of runs of the procedure against a WT sequence on a particular machine. In conclusion, an HRM-based method for interrogating the clinically significant codon 238 and codon 240 SNPs of the blaSHV gene has been developed. The method requires no probes, makes use of a low-cost Sybr green-based master mix, detects mutant alleles if they are in the minority, and provides an indication of the allele ratio. The method costs <$0.60 per sample in materials, even when carried out in duplicate, and takes approximately 1 h.
We thank John Turnidge and the SENTRY program for the SENTRY isolates and Graeme Nimmo and Queensland Health Pathology Service for the PAH isolates.
Published ahead of print on 6 April 2009. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»