Abstract
In this study, we investigated the molecular characteristics of extended-spectrum β-lactamase (ESBL)-producing Klebsiella pneumoniae isolates that were recovered from an outbreak in a Korean hospital. A new multilocus sequence typing (MLST) scheme for K. pneumoniae based on five housekeeping genes was developed and was evaluated for 43 ESBL-producing isolates from an outbreak as well as 38 surveillance isolates from Korea and also a reference strain. Overall, a total of 37 sequence types (STs) and six clonal complexes (CCs) were identified among the 82 K. pneumoniae isolates. The result of MLST analysis was concordant with that of pulsedfield gel electrophoresis. Most of the outbreak isolates belonged to a certain clone (ST2), and they produced SHV-1 and CTX-M14 enzymes, which was a different feature from that of the K. pneumoniae isolates from other Korean hospitals (ST20 and SHV-12). We also found a different distribution of CCs between ESBL-producing and -nonproducing K. pneumoniae isolates. The MLST method we developed in this study could provide unambiguous and well-resolved data for the epidemiologic study of K. pneumoniae. The outbreak isolates showed different molecular characteristics from the other K. pneumoniae isolates from other Korean hospitals.
Klebsiella pneumoniae is a relevant opportunistic pathogen that accounts for up to 10% of all nosocomial infections (1, 2). It is a frequent causal agent of urinary tract infections, septicemia, and pneumonia in immunocompromised patients, and it is also an important pathogen for community-acquired infections. Its importance in clinical settings is increasing due to its ability to produce extended-spectrum β-lactamase (ESBL) (3). It has been reported that ESBL-producing organisms tend to be multidrug-resistant and they have an increased risk of treatment failure (3, 4). Thus, it is indispensable to detect these ESBL-producing isolates and to control their emergence and spread.
Several molecular methods have been used for characterizing K. pneumoniae isolates such as randomly amplified polymorphic DNA analysis, pulsed-field gel electrophoresis, and amplified fragment length polymorphism (5-7). However, these methods are band-based and so the results may be ambiguous, but generally these methods are well suited for local epidemiologic investigation. Yet, it is very difficult to compare the results generated at different laboratories for conducting a global epidemiologic analysis. To overcome the limitations of the band-based typing methods, multilocus sequence typing (MLST) has been developed, which is a nucleotide sequence-based method for bacterial or fungal typing (8-10). It is highly discriminative and it is easy to standardize, store, and exchange information electronically. It has been applied successfully for the epidemiologic characterization of a variety of clinically important bacterial pathogens, such as Streptococcus pneumoniae, Staphylococcus aureus, Enterococcus faecium, Neisseria meningitidis, and Campylobacter jejuni (9-13). Very recently, a MLST method for K. pneumoniae by using seven housekeeping genes, such as rpoB, gapA, mdh, pgi, phoE, infB, and tonB, has also been developed (14).
We developed another MLST scheme for K. pneumoniae isolates based on five housekeeping genes, and we characterized the ESBL-producing K. pneumoniae isolates from an outbreak in a Korean hospital.
From August 2004 to May 2005, a total of 66 ESBL-producing K. pneumoniae were isolated from patients in an intensive care unit (ICU) of the Kangbuk Samsung Hospital (Seoul, Korea), which is a tertiary care hospital with 645 beds. Of these, 43 isolates were preserved and included in this study (Table 1). Twenty-nine isolates were from sputum, and 11 isolates were from urine; the others were isolated from blood or bile. Thirty-eight K. pneumoniae isolates, which were isolated at another three Korean tertiary hospitals as a part of a surveillance study that was conducted in 1999 and 2004, were also investigated in this study for comparisons (Table 2). Thirty-one isolates were collected from the Samsung Medical Center (Seoul, Korea), and others were from the Dong-A University Hospital (Busan, Korea) and the Chungbuk National University Hospital (Cheongju, Korea). Of these 38 surveillance isolates, 15 produced ESBL and the others did not. For screening for ESBL, the GNS-650 card (Vitek system, bio-Merieux, Hazelwood, MO, U.S.A.), which is based on the use of cefotaxime and ceftazidime with or without clavulanic acid, was used. For a reference, K. pneumoniae ATCC 13883 was also analyzed. Thus, a total of 82 K. pneumoniae isolates were analyzed.
We selected five housekeeping genes to establish the MLST scheme; rpoB (RNA polymerase β-subunit), gyrA (DNA gyrase subunit A), gapA (glyceraldehydes 3-phosphate dehydrogenase A), groEL (GroEL protein), and gyrB (DNA gyrase subunit B). The primers we used for amplification of the gene fragments are shown in Table 3. PCR amplifications were carried out under the following conditions: 35 cycles of denaturation at 94℃ for 30 sec, annealing at 50-55℃ for 30 sec, and extension at 72℃ for 1 min; this was all preceded by a 5 min denaturation at 94℃ and then it was all followed by a 5 min final extension at 72℃. The PCR products were purified using a PCR purification kit (Bioneer, Daejeon, Korea), according to the manufacturer's recommendations, and then the products were sequenced on an ABI3700 DNA sequencer (Applied Biosystems, Foster City, CA, U.S.A.).
The raw sequences were concatenated and edited by using the EditSeq and MegAlign programs (DNASTAR, Madison, WI, U.S.A.). For each locus, distinct allele sequences were assigned as an arbitrary allele number. Each isolate was characterized by its allelic profile, which was represented as a series of 5 integers corresponding to the alleles at each of the loci, in the order of rpoB, gyrA, gapA, groEL, and gyrB. The sequence type (ST) was designated for each unique allelic profile. Clusters of related STs sharing three or more identical alleles out of the five were defined as clonal complexes (CCs), which were determined by the eBURST program (http://eburst.mlst.net/loci.asp) (15). The dN/dS ratios, with dN being the number of nonsynonymous substitutions per nonsynonymous site and dS being the number of synonymous substitutions per synonymous site (16), were determined by using the START program (http://outbreak.ceid.ox.ac.uk/software.htm). The index of association (IA) was also calculated with the START program to test to what extent the K. pneumoniae population was clonal (17). IA is a measure of the variance in the number of allelic mismatches relative to that expected under the hypothesis of panmixia (random association), where the mean number of pairwise allelic mismatches is a measure of genetic distance. When the IA value does not deviate significantly from zero, the alleles are in linkage disequilibrium, i.e., there is frequent recombination or non-clonal population structure.
The 43 ESBL-producing K. pneumoniae isolates recovered from an outbreak at the Kangbuk Samsung Hospital were characterized by PFGE as described elsewhere (18, 19). Agarose plugs containing genomic DNA were digested with SmaI (Gibco, BRL, Gaithersburg, MD, U.S.A.), according to the manufacturer's recommendations. The DNA restriction fragments were separated with using a CHEF-Mapper apparatus (Bio-Rad Lab.), at 6 V/cm for 22 hr. The PFGE patterns were interpreted with using the published criteria (20). A divergence in more than three bands was interpreted as indicative of different PFGE types.
The bla genes related to ESBL enzymes were assayed by PCR with the corresponding primers for the TEM, SHV, CTX-M, and OXA ESBL types, as described previously (21). To specify the subtype of the bla genes, the amplified PCR products were sequenced on both strands as was done in the MLST analysis. The amino acid sequences were deduced from the nucleotide sequences using the MegAlign program (DNASTAR), and they were compared with the database of the website (http://www.lahey.org/Studies/).
Sequence data were obtained for the five selected loci from all 82 K. pneumoniae isolates included in this study. The nucleotide sequences at each locus ranged from 271 bp (gapA) to 752 bp (gyrA). The number of different alleles ranged from 9 (gapA) to 20 (gyrA). The proportion of sites at which nucleotide variation was observed within each locus that was examined ranged from 3.93% (rpoB) to 8.49% (groEL) (Table 4). By calculating the dN/dS ratio, the degree of selection operating on each locus was evaluated. No nonsynonymous substitutions were found in the gapA and gyrB loci, so the dN/dS ratios could not obtained for them. The dN/dS ratios for the other loci were significantly less than 1 (Table 4), indicating there was no strong positive selective pressure on the selected genes. Thus, these five loci could be appropriate for the K. pneumoniae MLST scheme. The index of association (IA) was used to test for linkage disequilibrium between the alleles at the five loci. The IA value was 1.951 for all the K. pneumoniae isolates. This value was reduced to 0.524 when including only one representative isolate for each sequence type (STs), and this reflects the moderate linkage equilibrium or clonal population structure of the K. pneumoniae isolates included in this study.
Based on the nucleotide variations of the five genetic loci, a total of 37 different STs could be identified among 82 K. pneumoniae isolates (Table 1, 2). The majority of these (29 out 37 STs, 78.4%) were represented by single isolate. Among the STs shared by multiple isolates, the most frequently encountered were ST2 (25 isolates, 30.5%), ST20 (8 isolates, 9.8%), and ST15 (5 isolates, 6.1%). Other STs (ST9, -11, -13, -14, and -21) represented two to four isolates. These STs were grouped, using the eBURST program, into CCs that consisted of all the isolates of the STs that shared three or more identical alleles (Fig. 1). As a result, six CCs were identified among the 82 K. pneumoniae isolates. CC20 included six STs (ST17, -18, -20, -21, -26, and -32), of which ST20 was assumed as a founder. Fourteen isolates (17.1%) were included in the CC20. CC2 included four STs (ST2, -13, -15, and -37), of which ST2 was assumed as a founder. Thirty-four isolates (41.5%) belonged to the CC2, but most of them were isolated at the Kangbuk Samsung Hospital during an outbreak. CC36 also included 7 isolates (8.5%) of four STs (ST4, -9, -31, and -36), of which ST36 was assumed as a founder. The other three CCs each consisted of two STs (Table 1, 2, Fig. 1).
PFGE was performed for 43 outbreak isolates from Kangbuk Samung Hospital. Among the 43 outbreak isolates, seven PFGE types (A to G) were identified (Fig. 2 and Table 1). The PFGE patterns were well concordant with the STs in the MLST. Two subtypes, AA and DD, could be differentiated from A and D, respectively. All 31 isolates of the PFGE type A belonged to CC2, which were differentiated into ST2 or ST15. ST2 and ST15 were single-locus variants of each other, differing at rpoB. One isolate (KS412-6) of PFGE type AA belonged to ST13, a single-locus variant of ST2 with a different allele in gyrA. PFGE type B was correlated with ST36 of CC36. Isolates of the PFGE type D and DD showed ST11 and ST14, respectively, both of which belonged to CC11-14. Both isolates of the PFGE type G belonged to CC20, but they could be differentiated into ST17 and ST18, which are different at gyrA. The other three PFGE types, C, E, and F, represented single isolates.
Of the 43 ESBL-producing K. pneumoniae isolates recovered from the Kangbuk Samsung Hospital, all but one produced SHV enzymes (29 isolates of SHV-1, 5 of SHV-5, 8 of SHV-11, and one each of SHV-12 and SHV-14). All but four isolates produced CTX-M-14 and 10 isolates were positive to blaTEM-1 (Table 1). One isolate (KS412-4) was negative on PCR assays for the four bla genes tested in this study.
All 33 outbreak isolates of CC2 and PFGE type A produced SHV and CTX-M enzymes. All the isolates produced CTXM-14, and all but three isolates produced SHV-1. The other three isolates produced SHV-5. Nine K. pneumoniae isolates produced TEM-1. Except for one isolate (KS502-5) of CC3-16 and PFGE type F, all the TEM-1-producing isolates belonged to CC2 and they showed PFGE type A. Four isolates of CC36 and PFGE type B produced SHV-11 and CTX-M-14. Of the three isolates of CC11-14, two of ST11 and the PFGE type D produced SHV-11 and CTX-M-14, while one of ST13 and PFGE type DD produced SHV-1 and CTX-M-14. Two isolates of CC20 and PFGE type G produced both SHV-11 and CTX-M-14, or only SHV-11.
Of the 38 K. pneumoniae isolates collected from three Korean hospitals as a part of a surveillance program, 15 isolates produced ESBL. All 15 ESBL-producing isolates produced SHV enzymes. Among these surveillance isolates, SHV-12 was the most frequent (10 isolates), and this was in contrast to the outbreak isolates, of which SHV-1 was the most common. Five isolates produced CTX-M-14, and 8 isolates were positive to blaTEM-1 (Table 2). Of the 15 ESBL-producing isolates, nine belonged to CC20, six isolates belonged to ST20, two belonged to ST21, and one belonged to ST26. There might be no obvious correlation between CC and the bla genes. Twenty-three K. pneumoniae isolates did not produce ESBL. These ESBL-negative isolates were more diverse than were the ESBL-producers. Only three isolates belonged to CC20 (ST20 and ST32), and three isolates belonged to CC2 (ST13 and ST37). Most of the isolates formed singletons in MLST analysis.
The objectives of this study were to demonstrate development of a typing method for K. pneumoniae and to investigate the characteristics of the ESBL-producing isolates recovered during an outbreak within a Korean hospital. Molecular typing is prerequisite for elucidating the epidemiology and population structure of bacterial pathogens. Out of several molecular typing methods, MLST is becoming popular owing to its several advantages over the other typing methods, which have been discussed repeatedly elsewhere; a high level of discrimination, it is unambiguous and reproducible, and its scalable nature resulting from use of nucleotide sequences, it has electronic portability via the Internet, it can easily analyze the generated data, and it has a wide applicability (9, 10, 22, 23). Very recently, the MLST scheme for K. pneumoniae was reported on and this was based on seven housekeeping genes (14). Diancourt et al. (14) characterized 67 K. pneumoniae isolates that were collected from several European countries with using their MLST scheme and they compared the results with those from ribotyping.
In this study, we presented another MLST scheme for K. pneumoniae based on five housekeeping genes. The appropriate nucleotide variations of each gene fragment (3.93% to 8.49%) indicated that these genes were applicable for the MLST scheme. The low of dN/dS ratio (0 to 0.1068), which indicated the neutral feature of the nucleotide substitution, also suggested that the selected housekeeping genes are probably suitable for population genetic studies of K. pneumoniae. Although our MLST scheme could not be compared with that of Diancourt et al. (14) because of the different loci and isolates, it showed results concordant with the PFGE. That is, the isolates with a distinct PFGE pattern belonged to a different ST. Moreover, some PFGE types could be differentiated into two STs, for example, the PFGE type A was differentiated into ST2 and ST15. Thus, our MLST procedure presented here is probably useful for a short-term epidemiologic study as well as for a global epidemiologic study, although only five loci were included in the scheme.
Using the MLST scheme, we investigated the molecular characteristics of ESBL-producing K. pneumoniae isolates from an outbreak within a Korean hospital. At the Kangbuk Samsung Hospital, at most three or less ESBL-producing K. pneumoniae had been isolated per month by July 2004. From August 2004 to May 2005, we observed an outbreak of ESBL-producing K. pneumoniae isolates primarily at the ICU. During this period, a total of 66 ESBL-producing K. pneumoniae were isolated at up to 14 isolates a month. Although we did not preserve the isolates found in August 2004, we could include the 43 ESBL-producing K. pneumoniae isolates from September 2004 to May 2005 in this study.
The results of MLST indicated that most of the isolates collected during the outbreak period belonged to the same clone, ST2 (allelic profile, 2-9-2-12-3). Particularly, all the isolates except two belonged to ST2 early during the outbreak period (September and November, 2004). After November 2004, ST2 was still the most prevalent, but its proportion decreased. In addition, five isolates emerged with its single-locus variant, ST15. ST15 (7-9-2-12-3) differed from ST2 by a single nucleotide of rpoB, which was probably due to a point mutation, and so it must have diverged from ST2 (24, 25).
Most of these outbreak isolates used in this study produced SHV and CTX-M enzymes as the ESBL enzymes (Table 1). A previous study performed in other Korean hospitals also reported that SHV and CTX-M were the most common ESBL enzymes (21). SHV-12 was the most frequently encountered SHV enzyme in the previous studies (21, 24-26), which was also confirmed in this study by the investigation of ESBL-producing K. pneumoniae isolates from Korean hospitals other than the Kangbuk Samsung Hospital (Tabel 2). However, SHV-11 was the most common among the outbreak K. pneumoniae isolates. This indicates that the outbreak K. pneumoniae isolates that produce ESBL are different from the prevailing isolates among the other Korean hospitals. However, we could not conclude whether this clone had previously existed in the Kangbuk Samsung Hospital and had spread during the outbreak period or whether it was transferred from another hospital or the community because we did not investigate the characteristics of the ESBL-producing K. pneumoniae isolates prior to the outbreak.
One of the interesting findings in this study was probably the different distribution of CCs between the ESBL-producing isolates and the non-producing isolates. Of the 12 ESBLproducing K. pneumoniae isolates from the Samsung Medical Center in 2004, eight isolates (66.7%) belonged to CC20. However, only three out of 23 ESBL-nonproducing isolates (13.0%) that were collected from the same hospital during the same period belonged to CC20 (Table 2). Although this result was based on the small number of isolates and thus it might be premature to draw a conclusion, it may suggest that the ESBL genes prefer to incorporate into a particular clone (s) as in methicillin-resistant S. aureus and vancomycinresistant E. faecium (12, 27, 28). Further investigation by including more isolates from diverse localities will be required to prove such a hypothesis.
In summary, we developed a new MLST scheme for K. pneumoniae based on five housekeeping genes. We investigated the molecular characteristics of ESBL-producing K. pneumoniae isolates from an outbreak within a Korean hospital with using the developed MLST method. Most of the outbreak isolates belonged to ST2 and they produced SHV-1 and CTXM-14, which was different from the features of the K. pneumoniae isolates from other Korean hospitals.
ACKNOWLEDGMENTS
Drs. Hyuck Lee (Dong-A University Hospital) and Jun Seong Son (Chungbuk National University Hospital) provided K. pneumoniae isolates.
References
1. Podschun R, Ullmann U. Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors. Clin Microbiol Rev. 1998. 11:589–603.
2. Spencer RC. Predominant pathogens found in the European Prevalence of Infection in Intensive Care Study. Eur J Clin Microbiol Infect Dis. 1996. 15:281–285.
3. Bradford PA. Extended-spectrum β-lactamases in the 21st century: characterization, epidemiology, and detection of this important resistance threat. Clin Microbiol Rev. 2001. 14:933–951.
4. Pitout JD, Nordmann P, Laupland KB, Poirel L. Emergence of Enterobacteriaceae producing extended-spectrum beta-lactamases (ESBLs) in the community. J Antimicrob Chemother. 2005. 56:52–59.
5. de Almeida VA, Pessoa-Silva CL, Sampaio JL, Gontijo Filho PP, Teixeira LM, Moreira BM. Genetic relatedness among extendedspectrum β-lactamase-producing Klebsiella pneumoniae outbreak isolates associated with colonization and invasive disease in a neonatal intensive care unit. Microb Drug Resist. 2005. 11:21–25.
6. Li CR, Li Y, Zhang PA. Dissemination and spread of CTX-M extended-spectrum beta-lactamases among clinical isolates of Klebsiella pneumoniae in central China. Int J Antimicrob Agents. 2003. 22:521–525.
7. van der Zee A, Steer N, Thijssen E, Nelson J, van't Veen A, Buiting A. Use of multienzyme multiplex PCR amplified fragment length polymorphism typing in analysis of outbreaks of multiresistant Klebsiella pneumoniae in an intensive care unit. J Clin Microbiol. 2003. 41:798–802.
8. Aanensen DM, Spratt BG. The multilocus sequence typing network: mlst.net. Nucleic Acids Res. 2005. 33:W728–W733.
9. Enright MC, Spratt BG. A multilocus sequence typing scheme for Streptococcus pneumoniae: identification of clones associated with serious invasive disease. Microbiology. 1998. 144:3049–3060.
10. Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, Feavers IM, Achtman M, Spratt BG. Multilocus sequence typing: a portable approach to identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998. 95:3140–3145.
11. Enright MC, Day NP, Davies CE, Peacock SJ, Spratt BG. Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J Clin Microbiol. 2000. 38:1008–1015.
12. Homan WL, Tribe D, Poznanski S, Li M, Hogg G, Spalburg E, van Embden JD, Willems RJ. Multilocus sequence typing scheme for Enterococcus faecium. J Clin Microbiol. 2002. 40:1963–1971.
13. Dingle KE, Colles FM, Wareing DR, Ure R, Fox AJ, Bolton FE, Bootsma HJ, Willems RJ, Urwin R, Maiden MC. Multilocus sequence typing system for Campylobacter jejuni. J Clin Microbiol. 2001. 39:14–23.
14. Diancourt L, Passet V, Verhoef J, Grimont PA, Brisse S. Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates. J Clin Microbiol. 2005. 43:4178–4182.
15. Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing. J Bacteriol. 2004. 186:1518–1530.
16. Ota T, Nei M. Variance and covariances of the numbers of synonymous and nonsynonymous nucleotide substitutions per site. Mol Biol Evol. 1994. 11:613–619.
17. Maynard Smith J, Smith NH, O'Rourke M, Spratt BG. How clonal are bacteria? Proc Natl Acad Sci USA. 1993. 90:4384–4388.
18. Ko KS, Baek JY, Lee JY, Oh WS, Peck KR, Lee N, Lee WG, Lee K, Song JH. Molecular characterization of vancomycin-resistant Enterococcus faecium isolated from Korea. J Clin Microbiol. 2005. 43:2303–2306.
19. Oh WS, Suh JY, Song JH, Ko KS, Jung SI, Peck KR, Lee NY, Yang Y, Chongthaleong A, Chiu CH, Kamarulzaman A, Parasakthi N, Lalitha MK, Perera J, Yee TT, Kumarasinghe G, Carlos CC. Fluoroquinolone resistance in clinical isolates of Streptococcus pneumoniae from Asian countries: ANSORP study. Microb Drug Resist. 2004. 10:37–42.
20. Tenover FC, Arbeit RD, Goering RV, Mickelsen PA, Murray BE, Persing DH, Swaminathan B. Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J Clin Microbiol. 1995. 33:2233–2239.
21. Kim J, Lim YM, Rheem I, Lee Y, Lee JC, Seol SY, Lee YC, Cho DT. CTX-M and SHV-12 β-lactamases are the most common extended-spectrum enzymes in clinical isolates of Escherichia coli and Klebsiella pneumoniae collected from 3 university hospitals within Korea. FEMS Microbiol Lett. 2005. 245:93–98.
22. Cooper JE, Feil EJ. Multilocus sequence typing - what is resolved? Trends Microbiol. 2004. 12:373–377.
23. Urwin R, Maiden MC. Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 2003. 11:479–487.
24. Feil EJ, Smith JM, Enright MC, Spratt BG. Estimating recombinational parameters in Streptococcus pneumoniae from multilocus sequence typing data. Genetics. 2000. 154:1439–1450.
25. Feil EJ, Cooper JE, Grundmann H, Robinson DA, Enright MC, Berendt T, Peacock SJ, Smith JM, Murphy M, Spratt BG, Moore CE, Day NP. How clonal is Staphylococcus aureus? J Bacteriol. 2003. 185:3307–3316.
26. Jeong SH, Bae IK, Lee JH, Sohn SG, Kang GH, Jeon GJ, Kim YH, Jeong BC, Lee SH. Molecular characterization of extended-spectrum beta-lactamases produced by clinical isolates of Klebsiella pneumoniae and Escherichia coli from a Korean nationwide survey. J Clin Microbiol. 2004. 42:2902–2906.
27. Katayama Y, Robinson DA, Enright MC, Chambers HF. Genetic background affects stability of mecA in Staphylococcus aureus. J Clin Microbiol. 2005. 43:2380–2383.
28. Oliveira DC, Tomasz A, de Lencastre H. Secrets of success of a human pathogen: molecular evolution of pandemic clones of methicillin-resistant Staphylococcus aureus. Lancet Infect Dis. 2002. 2:180–189.
29. Mollet C, Drancourt M, Raoult D. rpoB sequence analysis as a novel basis for bacterial identification. Mol Microbiol. 1997. 26:1005–1011.