INTRODUCTION
Klebsiella pneumoniae is an important human opportunistic pathogen that causes a wide variety of community-acquired and nosocomial infections (
1–4). The organism is still one of the leading causes of community-acquired pneumonia (
3) and has been described as the major cause of pyogenic liver abscesses (PLA) in some Asia-Pacific countries (
2,
4). Serotypes K1 and K2 in
K. pneumoniae have been linked to invasive infections (
2–9). Using multilocus sequence typing (MLST), Brisse et al. (
10) identified two virulent K1 clones of
K. pneumoniae, clonal complex (CC) 23 and CC82, which were strongly associated with PLA and respiratory infection, respectively.
K. pneumoniae is also an important reservoir of antibiotic resistance determinants (
1). Several clones of
K. pneumoniae with a high prevalence of carbapenem resistance or extended-spectrum β-lactamases (ESBLs) have been distributed worldwide (
11–16). The epidemic drug-resistant clones, called “high-risk” clones (
11), may play an important role in the dissemination of antimicrobial resistance. The major international high-risk clones identified by MLST belong to sequence type 11 (ST11), ST15, ST147, and ST258 (
11–16). ST11 and ST258 were included in clonal group (CG) 258 (
11,
14), while ST15 and ST147 were grouped in CG15 and CG147 (
11,
13,
14), respectively.
High-level carbapenem resistance in
K. pneumoniae may develop via the loss of the OmpK36 porin coupled with the expression of various β-lactamases (
17–19). The OmpK36 porin in
K. pneumoniae may also play a role in bacterial virulence (
18). The correlation of OmpK36 porin variants with specific STs of
K. pneumoniae was first described by Papagiannitsis et al. (
20). We then found that
K. pneumoniae isolates can be separated into four groups (designated groups A to D) by
ompK36 genotyping (
21). Interestingly, three major high-risk clones (ST11, ST15, and ST147) belonged to group A, while 6 PLA-associated STs (ST23, ST65, ST86, ST163, ST373, and ST375) were classified in group C (
21). The present study aimed to investigate the characteristics of the four groups of
K. pneumoniae isolates based on
ompK36 genotyping.
RESULTS
ompK36 typing.
Among the 226 study isolates, ompK36 types A, B, C, and D were identified by PCR in 61 (27.0%), 28 (12.4%), 100 (44.2%), and 36 (15.9%) isolates, respectively. The study isolates were separated into four groups based on the ompK36 types. Sequence analysis of ompK36 confirmed the result of PCR typing. PCR products were not observed with the 4 group-specific primers in 1 isolate. The comparative sequence analysis showed that the untypeable isolate possessed a rare ompK36 allele not belonging to the four major types (see Fig. S1 in the supplemental material).
MLST and PFGE analyses.
The MLST analysis showed 123 different STs among the 226 isolates (
Table 1). Twenty-six new alleles and 44 novel STs were identified in this study. The three most common STs were ST23 (12.8%, 29 isolates), ST65 (4.9%, 11 isolates), and ST15 (4.4%, 10 isolates); the other STs represented 1 to 7 isolates. Five STs (ST15, ST34, ST107, ST776, and ST1544) comprised isolates with different
ompK36 types. Further PFGE analysis showed that the ST15, ST34, ST107, and ST1544 isolates with different
ompK36 types had <80% similarity (
Fig. 1), suggesting the presence of subgroups in these STs.
The eBURST analysis with the 123 STs identified in this study alone showed 8 single-locus variant (SLV) pairs (ST35-ST1948, ST76-ST490, ST86-ST373, ST107-ST219, ST147-ST1922, ST355-ST1899, ST420-ST1544, and ST1918-ST1928), 7 ST clusters (ST14-ST15-ST709, ST23-ST1905-ST1893, ST37-ST726-ST1906-ST1198-ST884, ST65-ST1927-ST25-ST375, ST268-ST36-ST1914-ST1590, ST292-ST1925-ST1180, and ST1612-ST1915-ST1947), and 82 singletons. ST14 was an SLV of ST15, but they were assigned to different
ompK36 groups. ST14 and ST15 were subgroup founders in the eBURST analysis using the whole MLST database and have been allocated to CG14 and CG15 (
34). Moreover, 3 ST clusters (ST23-ST1905-ST1893, ST268-ST36-ST1914-ST1590, and ST292-ST1925-ST1180) included members with different
ompK36 types. The remaining three ST clusters and all SLV pairs contained STs with the same
ompK36 types (see Fig. S2 in the supplemental material).
Nineteen CGs were found to contain two or more STs detected in this study by the eBURST analysis, including CG15, CG22, CG23, CG34, CG35, CG37, CG45, CG65, CG76, CG86, CG107, CG147, CG268, CG292, CG347, CG420, CG1074, CG1612, and CG1928 (
Table 1). In general, most isolates within a single CG had the same
ompK36 type. All CG37 and CG147 isolates and all CG15 isolates except 1 ST15 and 1 ST277 isolates belonged to
ompK36 group A. All CG1612 isolates and both CG1928 isolates were allocated to group B. All CG23 isolates except 2 ST1905 isolates, all CG34 isolates except 1 ST34 isolate, both CG35 and CG76 isolates, all CG65 isolates except 1 ST685 isolate, and all CG86 isolates were assigned to group C. All CG268 isolates and all CG420 isolate except 1 ST1544 isolate belonged to group D.
Hypermucoviscous phenotype, capsular typing, and virulence gene content.
Among the 226 isolates, hypermucoviscosity was observed in 78 (34.5%) isolates, of which 62 (79.5%) isolates were classified in group C. Three (3.8%), 1 (1.3%), and 12 (15.4%) of the 78 mucoid isolates were in groups A, B, and D, respectively. The hypermucoviscous phenotype was predominant (62.0%) among the group C isolates, and the prevalence of the phenotype was significantly different among the four
ompK36 groups (
Table 2).
Thirty (13.3%) of the 226 isolates were identified as the K1 serotype by multiplex PCRs, 25 (11.1%) as K2, 3 (1.3%) as K5, 12 (5.3%) as K20, 1 (0.4%) as K54, and 7 (3.1%) as K57 (
Table 1). All K1, K2, and K5 isolates, 6 (85.7%) of the 7 K57 isolates, and 1 (8.3%) of the 12 K20 isolates belonged to group C (
Table 2). Group D was predominant among the 12 K20 isolates (83.3%, 10 isolates). The K54 isolate and 1 K57 isolate belonged to group B, and 1 K20 isolate was in group A.
Among the 10 tested virulence genes,
wzy_K1 (
magA),
allS, and
rmpA were significantly different among the four
ompK36 groups, and they were most frequent among group C isolates (
Table 2). Group C isolates had higher mean virulence scores than isolates of the other groups.
Antimicrobial resistance and resistance genes.
ESBLs were detected by the phenotypic and genotypic tests in 33 (14.6%) of the 226 isolates. Group A isolates were predominant among the 33 ESBL producers (66.7%, 22 isolates), and group C and D isolates accounted for 15.1% (n = 5) and 18.2% (n = 6) of all ESBL producers. The SHV-type ESBLs were identified in 26 (11.5%) of the 226 isolates; 17 (65.4%), 4 (15.4%), and 5 (19.2%) of the 26 isolates belonged to groups A, C, and D, respectively. Eight (3.5%) of the 226 isolates had the CTX-M-type ESBLs, including 2 CTX-M-3, 3 CTX-M-14, and 3 CTX-M-15; 6 (75.0%), 1 (12.5%), and 1 (12.5%) of the 8 CTX-M producers belonged to groups A, C, and D, respectively. The DHA-1 AmpC enzyme was detected in 20 (8.9%) of the 226 isolates. Group A was also predominant among the 20 DHA-1 producers (70.0%, 14 isolates); 2 (10.0%) and 4 (20.0%) isolates were in groups B and C, respectively. The most common ST among 44 isolates producing ESBLs and/or AmpC was ST15 (20.5%, 9 isolates), followed by ST11 (9.1%, 4 isolates) and ST709 (6.8%, 3 isolates). The other STs represented 1 to 2 isolates (see Table S1 in the supplemental material). Carbapenemase-producing isolates were not detected among the study isolates.
Overall, qnr genes were detected in 54 (23.9%) of the 226 isolates, including 1 qnrA1-like-positive, 6 qnrB1-like-positive, 18 qnrB4-like-positive, and 37 qnrS1-like-positive isolates. Also, group A isolates were the most common and accounted for 31 (57.4%) of the 54 isolates. Three (5.6%), 11 (20.4%), and 9 (16.7%) isolates were in groups B, C, and D, respectively. Group A was predominant among the qnrB-positive isolates (75.0%, 18 isolates) and the most common among the qnrS-positive isolates (45.9%, 17 isolates). Two isolates had the aac (6′)-Ib-cr gene, including 1 group A and 1 group B isolate.
Statistical analysis showed significant differences in antimicrobial resistance and resistance gene content among the four
ompK36 groups (
Table 3). Group A isolates had the highest nonsusceptibility rates against all tested antibiotics except meropenem and had the highest prevalence rates of SHV-type ESBLs, CTX-M ESBLs, the DHA-1-type AmpC enzyme,
qnrB, and
qnrS.
Plasmid typing.
Eight incompatibility groups of plasmids, IncFIA, I1-Iγ, L/M, N, P, A/C, F, and FIIA, were detected by the PCR-based replicon typing method in 1 (0.4%), 4 (1.8%), 2 (0.9%), 8 (3.5%), 2 (0.9%), 10 (4.4%), 2 (0.9%), and 60 (26.5%) of the 226 isolates, respectively. The four groups of
K. pneumoniae differed only in the carriage of IncFIIA (
Table 3), and the prevalence rate of IncFIIA-type plasmids was highest in group A. It was noted that 22 (84.6%) of the 26
blaSHV-ESBL-positive isolates, 4 (50.0%) of the 8 CTX-M-producers, and 15 (75.0%) of the 20 DHA-1 producers had the IncFIIA-type plasmids, suggesting that the IncFIIA-type plasmids might play an important role in the spread of
blaESBL and
ampC in
K. pneumoniae in this institution.
Classification of bloodstream infections.
Bloodstream infections caused by the isolates of the four
ompK36 groups were analyzed, and statistical analysis showed significant differences in the possible site of infection acquisition between the four groups (
Table 4). Group C isolates were more likely to be community acquired (40.0%) than the isolates of the other groups, while group A isolates tended to be nosocomially acquired (50.8%). Twelve (42.9%) of the 28 group B isolates were responsible for bacteremic biliary tract infections, and the rate was significantly higher than the rates in other groups (8.2 to 14.0%) (
Table 4). Twenty-one currently available STs were identified among the 44 isolates causing bacteremic biliary tract infections; each of the STs represented 1 to 2 isolates (see Table S2 in the supplemental material). Seventeen isolates were responsible for PLA-associated bacteremias, and 16 (94.1%) of them were in group C. ST23 was predominant among the 17 PLA-associated isolates (76.5%, 13 isolates) (see Table S2 in the supplemental material). There were no significant differences in bacteremias related to pneumonia, urinary tract infections, peritonitis, and soft tissue infections among the four groups. The majority of bacteremic PLAs, biliary tract infections, soft tissue infections, peritonitis, and urinary tract infections were community onset, and the vast majority (82.4%) of PLAs were community acquired (see Table S3 in the supplemental material). Most of the bacteremic pneumonias were either nosocomial acquired (50.0%) or health care-associated community onset (38.9%). Fifty-four patients who died within 30 days of the first positive blood culture with
K. pneumoniae and 158 patients known to have survived for more than 30 days were analyzed. The 30-day mortality rates of groups A, B, C, and D were 28.1% (16 of 57 patients), 22.2% (6 of 27 patients), 28.7% (27 of 94 patients), and 14.7% (5 of 34 patients), respectively, and the difference was not statistically significant between the four groups (
P = 0.399).
DISCUSSION
Genetic polymorphism of
ompK36 has been described in
K. pneumoniae (
19,
20), and the
ompK36 variants can be divided into four major groups (
21). The present study demonstrated that the vast majority of clinical
K. pneumoniae isolates can be allocated into four distinct groups based on
ompK36 types; only one isolate with a rare
ompK36 type was not assigned to the four groups by the PCR typing method. Moreover, this study only enrolled isolates from a single institution in Taiwan. It is thus very likely that there are still other rare
ompK36 alleles or groups that cannot be defined by the PCR typing method.
Group A included all isolates of three genetically unrelated epidemic resistant clones (ST11, ST15, and ST147) except 1 ST15 isolate. Among the 17 isolates of the three high-risk clones, 13 (76.5%) isolates were ESBL and/or AmpC producers. The inclusion of major high-risk clones may explain partially the distinct characteristics of group A: high prevalences of antimicrobial resistance, ESBLs, AmpC, Qnr-type quinolone resistance determinants, and the IncFIIA-type plasmids. Moreover, among the four groups, group A was most likely associated with nosocomial infections. The characteristics of this group suggest that group A isolates have been more successfully adapted to health care environments than others. The important epidemic clone ST258 (an SLV of ST11) is not common in Taiwan (
16) and was not detected in our present and previous studies (
21). The
ompK36 type of ST258 and the possibility that clones other than ST11, ST15, and ST147 in this group become high-risk clones deserve further investigation.
The capsular serotypes K1 and K2 and the hypermucoviscous phenotype are important virulence determinants of
K. pneumoniae (
4,
6,
7), and community-onset invasive infections are often caused by isolates with those characteristics (
4,
5,
8,
9). Previous studies have demonstrated that among non-K1/K2
K. pneumoniae isolates, serotype K5, K20, K54, and K57 isolates are common causes of community-onset pneumonia and PLA, and the majority of them had the hypermucoviscous phenotype (
5,
9). It is interesting to find that the majority of mucoid isolates, all K1, K2, and K5 isolates, and most K57 isolates were included in group C. ST23, ST65, ST86, ST373, and ST375 have been reported to be associated with PLA (
4,
7), and all isolates of these virulent clones (29 ST23, 11 ST65, 5 ST86, 7 ST373, and 1 ST375 isolates) were also assigned to this group (see Table S2 in the supplemental data). The inclusion of hypervirulent clones can explain the association of this group with community-onset infection and PLA.
It is also interesting to find that as many as 42.9% of group B isolates were responsible for bacteremic biliary tract infections, and no endemic clone was identified in this group. K20 was the third most common serotype in this study, and 11 K20 isolates were detected in this study. Ten of the 11 K20 isolates, which belonged to four STs, were classified in group D. The rate of bacteremic urinary tract infection-associated isolates was the highest in group D; however, the difference among the four ompK36 groups was not statistically significant. A better understanding the characteristics shared by different clones of K. pneumoniae linked to a specific infection might be helpful to identify microbiological factors contributing to the specific type of infection. One major limitation of this study was that only bloodstream isolates, which might represent more virulent strains, were analyzed. Another limitation was the lack of an adequate number of patients for subset studies to investigate the usefulness of our typing method in clinical practice. More studies using isolates from various sample types are therefore needed to confirm the link between specific types of infection and ompK36 allele groups in K. pneumoniae and to investigate if our typing method can be utilized to predict the progression of primary infection or the prognosis of bacteremia in various clinical subsets.
The OmpK36 porin plays an important role in noncarbapenemase-mediated carbapenem resistance (
17–19). Moreover, previous studies revealed the association of OmpK36 with bacterial virulence (
18,
35). The
ompK36-deficient mutant was shown to confer a lower virulence in a mouse peritonitis model (
18), and the loss of OmpK36 may confer decreased resistance to neutrophil phagocytosis and increased resistance to serum killing (
35). It is therefore likely that the separation of high-risk resistant clones and virulent clones into two distinct groups and similar
ompK36 types among genetically unrelated epidemic clones were the consequence of convergent evolution driven under the selective pressures created by antimicrobial use and host factors.
A prevalent multidrug-resistant clone may either derive from a preexisting successful, but antibiotic-susceptible, strain with subsequently acquired resistance determinants or come from an initially minor strain with resistance as the key driver for its prevalence (
11). ST11, ST15, ST37, and ST147 are important STs among carbapenem-resistant
K. pneumoniae isolates in Taiwan, and they were reported to account for >70% of carbapenem-resistant isolates by Chiu et al. (
16). In this study, the four STs only accounted for 24 (10.6%) of all study isolates (
Table 1), suggesting that these high-risk resistant clones had been minor strains rather than evolving from preexisting successful strains. Active continuous surveillance followed by appropriate infection control measures is needed to prevent the further spread of these clones and dissemination of resistance determinants, and the
ompK36 typing methods may be useful for screening of the high-risk clones.