Abstract
Various resistance mechanisms facilitate the emergence and spread of multidrug-resistance (MDR) phenotypes of Escherichia coli and Klebsiella pneumoniae. To elucidate the MDR mechanisms of E. coli and K. pneumoniae in China, we analyzed the antimicrobial susceptibilities of strains isolated from clinical samples in a large tertiary care hospital in Beijing, China, during 2007–2009 and characterized the isolates with a cefotaxime–ciprofloxacin–amikacin (CTX–CIP–AK) resistance pattern. In total, 98 and 52 clinical isolates of E. coli and K. pneumoniae, respectively, with a CTX–CIP–AK resistance pattern were subjected to antimicrobial susceptibility testing and screening of common β-lactamase genes, plasmid-mediated quinolone resistance (PMQR) genes, quinolone resistance-determining region (QRDR) substitutions, and 16S rRNA methylase genes by polymerase chain reaction amplification and DNA sequencing. Pulsed-field gel electrophoresis (PFGE) was used to determine the genetic relatedness of the isolates. Approximately 6.86% and 8.05% of the clinical E. coli and K. pneumoniae isolates, respectively, exhibited MDR phenotypes. The MDR K. pneumoniae isolates exhibited significantly higher ceftazidime resistance than the MDR E. coli isolates (90.4% vs. 76.5%, p=0.0339); a similar result was noted for piperacillin–tazobactam resistance (28.8% vs. 2%, p=0.0001). The common resistance determinants among the MDR E. coli and K. pneumoniae isolates were as follows: CTX-M (88.8% vs. 82.7%), PMQR genes (70.4% vs. 90.4%), gyrA mutations (100% vs. 90.4%), and 16S rRNA methylase genes (93.9% vs. 94.2%). Half (50%) of the MDR E. coli isolates belonged to phylogenetic group D, followed by group A (39.8%). For the E. coli isolates, 94 PFGE patterns and 23 clusters were identified, whereas 51 PFGE patterns and 11 clusters were identified for the K. pneumoniae isolates. Clinical E. coli and K. pneumoniae isolates seem to have a low prevalence of MDR phenotypes in China. The great genetic variation indicates a considerable transmission of common resistance determinants, including a high prevalence of QRDR substitutions in E. coli and K. pneumoniae.
Introduction
A growing body of literature evidences the mechanisms of drug resistance in clinical E. coli and K. pneumoniae isolates in China.16,17,28–30,32,33 Plasmid-mediated quinolone resistance (PMQR), β-lactamase, and 16S rRNA methylase genes have been identified as the genetic bases for bacterial drug resistance. However, little is known about the MDR mechanisms in clinical E. coli and K. pneumoniae strains in China. To address this issue, we analyzed the antimicrobial susceptibilities of E. coli and K. pneumoniae strains isolated from clinical samples in Beijing, China, during 2007–2009 and characterized the isolates with a cefotaxime–ciprofloxacin–amikacin (CTX–CIP–AK) resistance pattern. The prevalence of resistance determinants in these isolates was also investigated.
Materials and Methods
Bacterial isolates
In 2007–2009, 2,346 and 994 individual E. coli and K. pneumoniae strains, respectively, were isolated from clinical samples, including sputum, blood, urine, abscesses, secretions, drainage fluids, bile, wounds, and catheters, in a 4,000-bed tertiary care hospital. All isolates were identified by colony morphology, biochemical testing, and VITEK 2 GN ID cards (bioMérieux, Inc.). Of these, 161 E. coli and 80 K. pneumoniae isolates exhibited a CTX–CIP–AK resistance pattern, as indicated by antimicrobial susceptibility testing and the Clinical and Laboratory Standards Institute (CLSI) standards. 6 In early 2011, 98 E. coli and 52 K. pneumoniae isolates with a CTX–CIP–AK resistance pattern were subjected to phenotypic and genotypic characterization. E. coli ATCC 25922 and ATCC 35218 were used as the quality control strains for antimicrobial susceptibility testing. The sodium azide-resistant E. coli J53 strain was used as the recipient for conjugation testing.
No formal ethical approval was obtained to use the clinical samples, because they were collected during routine bacteriologic analyses in public hospitals, and the data were anonymously analyzed.
Antimicrobial susceptibility testing
Antimicrobial susceptibilities of the clinical isolates were determined by measuring the diameters of the zones of complete inhibition in the disk diffusion method. The minimum inhibitory concentrations (MICs) for the isolates and transconjugants were measured by using the agar dilution method. All susceptibility results were interpreted according to the 2010 CLSI performance standards. 6
Screening of resistance determinants
β-lactamase genes, including blaTEM, blaSHV, blaCTX-M (groups 1, 2, 8, 9, and 25/26), blaPER, blaVEB, blaGES, and plasmid-mediated blaAmpC, and PMQR genes, including qnrA, qnrB, qnrC, qnrD, qnrS, qepA, oqxAB, and aac(6′)-Ib-cr were screened by multiplex polymerase chain reaction (PCR) methods as previously described.3,7,15,19 PCR amplification was performed as previously described to screen 16S rRNA methylase genes, including armA, rmtA, rmtB, rmtC, and rmtD. 8 Further, PCR amplification and DNA sequencing were performed to identify gyrA and parC mutations. 9 Amplicons were sequenced in a DNA analyzer (3730xl; Applied Biosystems, Life Technologies).
Phylogenetic analysis
The phylogeny of the E. coli isolates was determined on the basis of the presence or absence of chuA, yjaA, and TspE4.C2 by using a triplex PCR method. 5
Conjugation testing
Conjugation experiments were carried out in Luria–Bertani (LB) broth. Cultures of donor and recipient cells in the logarithmic phase (0.5 ml of each) were added to 4 ml of fresh LB broth and incubated overnight without shaking. Transconjugants were selected on a China blue agar plate containing AK (30 μg/ml) and sodium azide (170 μg/ml; Sigma-Aldrich Co.).
Pulsed-field gel electrophoresis
DNA fingerprints were obtained from pulsed-field gel electrophoresis (PFGE) profiles of genomic DNA digested with XbaI (New England Biolabs) according to the procedures developed by the US Centers for Disease Control and Prevention PulseNet program. 21 The PFGE patterns were interpreted with BioNumerics software (Applied Maths NV) by using the Dice similarity coefficient. A tree indicating relative genetic similarity was constructed on the basis of the unweighted pair-group method of averages with a position tolerance of 1%. Clusters were defined as DNA patterns sharing ≥70% similarity (C1, C2, C3, etc.). Differences in genetic similarity by <5% were considered as representing subtypes within the main group.
Statistical analysis
Pearson's χ 2 -test (for frequencies of higher than five) or Fisher's exact test (for small contingency tables) was used in the Chinese High Intellectualized Statistical Software version 2010 to determine the significance of the prevalence values.
Results
Prevalence of MDR phenotypes
In 2007–2009, 6.86% (161/2,346) E. coli and 8.05% (80/994) K. pneumoniae isolates exhibited MDR phenotypes. Overall, 5.89%, 7.01%, and 7.65% of the clinical E. coli isolates and 8.46%, 8.12%, and 7.56% of the clinical K. pneumoniae isolates exhibited MDR phenotypes in 2007, 2008, and 2009, respectively.
Antimicrobial susceptibilities
Among the 98 E. coli and 52 K. pneumoniae isolates with MDR phenotypes, most were susceptible to imipenem (98% vs. 96.2%) and meropenem (99% vs. 98.1%). These MDR K. pneumoniae isolates exhibited remarkably high ceftazidime (CAZ) resistance compared with the E. coli (90.4% vs. 76.5%, p=0.0339); the cefoxitin resistance was also high (73% vs. 62.2%, p=0.1704). Further, 52% and 92.3% of these E. coli and K. pneumoniae isolates, respectively, were unsusceptible (intermediate or resistance) to piperacillin–tazobactam.
Phylogenetic typing of MDR E. coli isolates
Most E. coli isolates with MDR phenotypes belonged to phylogenetic group D (50%), followed by group A (39.8%), group B1 (5.1%), and group B2 (5.1%).
Prevalence of resistance determinants
Table 1 lists the common resistance determinants. All isolates were negative for blaPER, blaVEB, blaGES, rmtA, rmtC, rmtD, rmtE, and npmA. In addition, 56.1% (55/98) of the MDR E. coli isolates carried qepA and rmtB simultaneously, and 58.2% (32/55) of them were co-transferred by conjugation. Further, 44.2% (23/52) of the MDR K. pneumoniae isolates carried qnrB and armA simultaneously. All (26/26) of the DHA-producing MDR K. pneumoniae isolates were qnrB positive.
Phenotype detected by the double-disk synergy test.
ESBL, extended-spectrum β-lactamase; QRDR, quinolone resistance-determining region; PMQR, plasmid-mediated quinolone resistance.
Transferability of resistance determinants and antimicrobial susceptibility of transconjugants
In total, 74 (75.5%) and 31 (59.6%) transconjugants were obtained from the MDR E. coli and K. pneumoniae isolates, respectively. Among the MDR E. coli isolates, 32.2% (28/87) of blaCTX-M, 42% (29/69) of the PMQR genes, and 56% (51/91) of the 16S rRNA methylase genes were located on transferable plasmids. Further, 60 of 74 (81.1%) transconjugants harbored at least two classes of these resistance determinants, and 14 of 74 (19.1%) transconjugants carried all three classes. For the MDR K. pneumoniae isolates, 34.9% (15/43) of blaCTX-M, 14.6% (7/48) of blaSHV, 30.8% (8/26) of blaDHA, 34% (16/47) of the PMQR genes, and 38.8% (19/49) of the 16S rRNA methylase genes were located on transferable plasmids. Approximately 67.7% (21/31) of the transconjugants harbored at least two classes of these resistance genes, and 32.2% (10/31) of the transconjugants carried all three classes of resistance determinants.
The transconjugants carrying armA or rmtB exhibited similar AK susceptibility to the armA- and rmtB-negative transconjugants. The CTX susceptibility of the transconjugants was significantly affected by the production of CTX-M. CTX MIC50 and MIC90 of the CTX-M-producing transconjugants were 32 and 256 mg/l, respectively, whereas those of the CTX-M-negative strains were 2 and 64 mg/l, respectively. However, the CAZ susceptibility of the transconjugants was rarely affected by CTX-M production: CAZ MIC50 of the CTX-M-producing and CTX-M-negative transconjugants was 4 and 2 mg/l, respectively, and MIC90 of these transconjugants was 64 and 32 mg/l, respectively. The CIP MIC of the transconjugants was slightly affected by the PMQR genes.
Genetic relatedness
Among the 98 MDR E. coli isolates, 94 PFGE patterns and 23 clusters were identified; more than one isolate was found in four PFGE patterns. Among the 52 MDR K. pneumoniae isolates, 51 PFGE patterns and 11 clusters were identified; more than one isolate was found in one PFGE pattern (Fig. 1).

Dendrogram of patterns generated by pulsed-field gel electrophoresis of multidrug-resistance Klebsiella pneumoniae isolates.
Discussion
In this study, the clinical E. coli isolates (6.86%) exhibited a lower prevalence of MDR phenotypes than did the clinical K. pneumoniae isolates (8.05%). However, from 2007 to 2009, the E. coli isolates presented increased MDR prevalence (5.89%–7.68%), whereas the K. pneumoniae isolates showed decreased MDR prevalence (8.46%–7.79%). Considering that both types of MDR isolates were sourced from the same antibiotic use scenarios, the reason for the differential prevalence needs further analysis.
The great genetic variation of the MDR isolates indicates the considerable transmission of common resistance determinants. Despite a similar prevalence of blaCTX-M between the two species (88.8% vs. 82.7%), the MDR E. coli isolates exhibited lower CAZ resistance than did the MDR K. pneumoniae isolates. This difference could be explained by the high prevalence of the SHV-type enzyme among MDR K. pneumoniae isolates, which exhibit strong hydrolytic activity against CAZ compared with CTX, rather than by the effect of CTX-M, because this enzyme slightly affected the CAZ MIC of the transconjugants. 22 Another explanation is that more than half of the MDR K. pneumoniae isolates produced DHA and CMY, which can efficiently hydrolyze cephamycins,2,13 whereas only a few MDR E. coli isolates produced these enzymes (Table 1). Further, the high prevalence of the plasmid-mediated AmpC-type enzymes contributed to the reduced susceptibility to the β-lactam/β-lactamase inhibitor combinations, such as piperacillin–tazobactam, because of the functional traits of AmpC enzymes that are not inhibited by β-lactamase inhibitors but hydrolyze piperacillin significantly.1,2,13 Almost half of the MDR E. coli isolates were susceptible to piperacillin–tazobactam, because of the low prevalence of AmpC-type β-lactamases, whereas only 7.7% of the MDR K. pneumoniae isolates were susceptible to this agent.
In this study, the absolute fluoroquinolone resistance was the result of the high prevalence of QRDR substitutions, because most MDR isolates possessed at least 1 QRDR mutation in both gyrA and parC. However, the two species exhibited significant differences in the prevalence of these substitutions (Table 1). The low frequency of QRDR substitutions in the MDR K. pneumoniae isolates could partly be explained by the high prevalence of qnr (71.2%), because Qnr reduces the amount of enzyme–DNA complexes, which are the target of fluoroquinolones, and protects QRDR domains from mutation.4,25,26,31 PMQR genes were highly prevalent among the MDR E. coli (70.4%) and K. pneumoniae isolates (90.4%) and played an important role in their fluoroquinolone-resistant phenotype. For example, K. pneumoniae strains E9108, E9193, E9194, E9196, and E9218 did not have either gyrA or parC mutations. The mechanisms for their fluoroquinolone-resistant phenotype might include efflux pumps or membrane impermeability, or the generation of synergistic resistance effects with different PMQR genes. 18 In comparison with previous reports,15,32 our study showed a higher prevalence of PMQR genes, indicating a considerable impact of PMQR mechanisms in QRDR substitutions in MDR E. coli and K. pneumoniae strains. In addition, due to the coexistence of PMQR genes with extended-spectrum β-lactamase (ESBL) or carbapenemase genes and their cross-species/genera transferability, 18 a high prevalence of PMQR mechanisms could facilitate the spread of MDR isolates. In this study, qepA exhibited a strong association with rmtB among the MDR E. coli isolates; a close relationship was observed between qnrB and armA and blaDHA among the MDR K. pneumoniae isolates.
High AK susceptibility has been reported for clinical E. coli and K. pneumoniae isolates.12,28 As a semisynthetic aminoglycoside, AK is rarely inactivated by common aminoglycoside-modifying enzymes, 14 making it a good indicator of the presence of 16S rRNA methylase genes in clinical gram-negative isolates. 8 In this study, armA and rmtB played a major role in AK resistance, because most MDR isolates contained either one or both of these genes. However, 6.1% and 5.8% of the MDR E. coli and K. pneumoniae isolates, respectively, tested negative for the known 16S rRNA methylase genes, indicating that other mechanisms also confer AK resistance in E. coli and K. pneumoniae, such as unknown 16S rRNA methylase variants or several isozymic forms of aminoglycoside phosphotransferase.23,24 This finding was confirmed by conjugation testing, because the transconjugants not carrying armA or rmtB also exhibited high-level AK resistance.
In conclusion, the clinical E. coli and K. pneumoniae isolates exhibited diverse phenotypic characteristics because of the differential prevalence of common resistance determinants. Despite the relatively low prevalence in this study, a larger-scale epidemiological investigation on the prevalent characteristics of clinical MDR E. coli and K. pneumoniae isolates and further research on the mechanisms of resistance transmission are warranted to prevent the spread of MDR strains.
This study has several limitations. First, the retrospective investigation was carried out with limited information of the strains isolated from clinical samples. Inadequate sample collection and handling and, especially, etiological diagnoses based on the subjective analyses of patients affect the ability to perform randomized clinical trials, which may be a major source of selection bias when investigating the prevalence of MDR phenotypes in clinical isolates. Second, only 98 E. coli and 52 K. pneumoniae isolates were further analyzed in this study. The loss of strains because of inadequate preservation and contamination could lead to bias. Third, a detailed analysis is required to understand the mechanisms related to MDR phenotypes, such as the Inc-groups of resistance plasmids, sequence type of MDR isolates, and, especially, variant types of resistance determinants. Fourth, the selection of transconjugants on AK only (and not on other antibiotics) could result in a significant underestimation of the contribution of conjugative plasmids to the resistance phenotype.
Footnotes
Acknowledgment
This research was supported by grant 2008ZX10004-001-C from the Ministry of Science and Technology, China.
Authors' Contributions
J.Y. conducted the molecular genetic studies and drafted the manuscript. Y.Z., L.Y., and Y.L. participated in the strain isolation and antimicrobial susceptibility testing. W.W. participated in the molecular genetic studies. W.Z. participated in the sequencing of resistance determinants. Z.C. participated in the PFGE analysis. L.H. helped design of the study. All authors read and approved the final manuscript.
Disclosure Statement
The authors declare that they have no competing interests.
