Na Wang1,2
Lexiu Deng3
Huiying Li4Na Jia5
Xiaocui Peng3
Jianliang Chang3
Jiatong Hao3
Jianhua Tang6
Chunmei Lei1
Bu Wang3*
Jianhua Liu3*
Wei Zhang7*- 1Microbiology Department, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
- 2Infection Management Department, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
- 3Respiratory Department, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
- 4Obstetrics and Gynecology Department, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
- 5Clinical Laboratory, Tianjin Nankai Tianyun Hospital, Tianjin, China
- 6Department of Pharmacy, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
- 7Central Laboratory, The First Affiliated Hospital of Hebei North University, Hebei North University, Zhangjiakou, Hebei, China
Background: Carbapenem-resistant hypervirulent Klebsiella pneumoniae (CR-hvKP) represents a critical public health threat in China, characterized by the convergence of multidrug resistance and hypervirulence. The emergence of ceftazidime/avibactam (CZA) resistance further complicates clinical management. This study aimed to elucidate the molecular epidemiology, risk factors, and resistance mechanisms of CZA resistance in CR-hvKP across China, providing evidence for targeted interventions.
Methods: A single-center molecular epidemiological analysis was conducted on 81 Carbapenem-resistant Klebsiella pneumoniae (CRKP) clinical isolates collected. All isolates underwent whole-genome sequencing for MultiLocus Sequence Typing, capsule typing, and identification of resistance genes (blaKPC-2 and blaNDM-1) and virulence factors (iucA, iroB, rmpA, rmpA2, and peg-344). CZA resistance mechanisms were investigated through broth microdilution minimum inhibitory concentration (MIC) testing and bioinformatics analysis. Galleria mellonella infection models were employed to assess virulence potential. Risk factors were analyzed using multivariate regression of clinical variables. Phylogenetic reconstruction employed single-nucleotide polymorphism-based analysis.
Results: ST11 accounted for 96.15% (50/52) of CR-hvKP isolates, with K64 being the predominant capsule type (92.31%, 48/52). Additionally, 98.77% (80/81) of CRKP carried ≥1 virulence gene; 64.2% (52/81) of isolates with all five virulence genes exhibited lethality. Galleria mellonella revealed that the survival rate of CR-hvKP was lower than that of carbapenem-resistant non-hypervirulent Klebsiella pneumoniae (p<0.05). Antibiotic usage time (odds ratio [OR]=1.076, 95% confidence interval [CI]: 1.026–1.138), carbapenem antibiotic (OR = 0.117, 95% CI: 0.02266–0.4602), and malignant tumors (OR = 65.1, 95% CI: 7.078–1798) predicted CR-hvKP infection. Transferable blaKPC-2 on IncFII/IncR plasmids conferred CZA resistance (MIC>128 mg/L) without compromising carbapenem resistance, facilitated by a unique genetic context (TnpR_Tn3-ISKpn27-blaKPC-2-ISKpn6).
Conclusion: China faces a rapid dissemination of ST11 CR-hvKP clones carrying diversified CZA resistance mechanisms. The convergence of hypervirulence and resistance in ST11 lineages—accelerated by invasive procedures and international transmission—demands enhanced genomic surveillance. CZA resistance arises through multiple pathways, necessitating combination therapies and stewardship programs limiting prolonged CZA use. Our findings underscore an urgent need for rapid diagnostics targeting emergent resistance determinants and infection control measures to contain high-risk clones.
Graphical Abstract.
1 Introduction
Carbapenem-resistant Klebsiella pneumoniae (CRKP) has become a major global public health threat, and the prevalence of carbapenem-resistant hypervirulent Klebsiella pneumoniae (CR-hvKP) has sharply increased in China, leading to high morbidity and mortality rates due to the lack of effective treatment options (Zhang et al., 2018; Yang et al., 2022). Klebsiella pneumoniae Carbapenemase (KPC) belongs to class A serine carbapenemases and is the main mechanism of resistance in Enterobacteriaceae. The proportion of blaKPC-2 producing CRKP in China exceeds 70% (Yang et al., 2022). In China, ST11 CRKP has become the major clonal type (Hu et al., 2024), whereas in the United States and Europe, ST258/512 is predominant (Tryfinopoulou et al., 2023). Owing to ceftazidime/avibactam (CZA) having good in vitro sensitivity and safety against serine carbapenemases, it has become the last line of defense against CRKP infections. However, CRKP can generate blaKPC mutations via various mechanisms, causing resistance to CZA (Pu et al., 2023a; Shi et al., 2024). Research has found that the diversity of the CRKP genome is mainly due to horizontal transfer, including plasmids, phages, integration, and binding elements. The spread of blaKPC-2 is typically mediated by two mobile elements, Tn4401 and NTEKPC (Yang et al., 2021).
Globally, CR-hvKP is increasingly reported, and most clinical cases occur in Asia, particularly in China, where studies have shown that the prevalence of CR-hvKP has increased from 28.2% in 2016 to 45.7% in 2020 (Hu et al., 2024), causing serious nosocomial infection in the intensive care unit ward (Gu et al., 2018). CR-hvKP virulence factors include capsules, siderophores, virulence plasmids, and other virulence genes, which help differentiate hypervirulent strains (Pu et al., 2023b). Various complex mechanisms have led to the rapid spread of strains in the clinic (Pu et al., 2023b). Previous studies have mainly focused on resistant or virulent strains of hvKP, with little attention paid to the correlations among virulence genes, resistance genes, and antimicrobial susceptibility of CR-hvKP. Therefore, this study aimed to summarize the molecular epidemiological characteristics of CR-hvKP isolated from this region and discuss the evolution of virulence and resistance genes in CR-hvKP and their relationship with clinical phenotypes.
2 Materials and methods
2.1 Isolate and antimicrobial susceptibility testing
In this study, 81 non-duplicate CRKP strains were collected from the First Affiliated Hospital of Hebei North University between 2021 and 2022 for analysis (Figure 1). The PubMLST database was downloaded for comparative analysis of Klebsiella pneumoniae containing blaKPC-2 and/or blaNDM-1 in November 2024 (https://bigsdb.pasteur.fr/klebsiella/). This study employed the microbroth dilution method recommended by the Clinical and Laboratory Standards Institute to perform in vitro AST of 81 CRKP strains. Sensitivity, intermediate resistance, and resistance were determined according to the Clinical and Laboratory Standards Institute-M100 ED33 guidelines. A minimum inhibitory concentration (MIC) of ≥4 mg/L for imipenem or meropenem against KP was defined as CRKP.
Figure 1. Design and experimental flowchart of this study.
2.2 Polymerase chain reaction
Polymerase chain reaction experiments were conducted on the 81 collected CRKP strains using the following reaction system: The total reaction volume was 25 μL, which included 12.5, 8.5, 1, 1, and 2 μL of 2x Taq Plus Master Mix II (Dye Plus), deionized water, forward primer, reverse primer, and template DNA, respectively. The reaction program was set as follows: initial denaturation at 95°C for 3 min, denaturation at 95°C for 15 s, annealing at 60°C for 20 s, and extension at 72°C for 60 s, with a total of 30 cycles. After the reaction, agarose gel electrophoresis was performed to record the experimental results.
2.3 Whole genome sequencing and annotation
Genomic DNA was extracted using the sodium dodecyl sulfate method, and the harvested DNA was analyzed through agarose gel electrophoresis. Quantitative analysis of the DNA was performed using a Qubit® 2.0 fluorometer (Thermo Scientific). Sequencing libraries were prepared using the NEBNext Ultra DNA Library Prep Kit (Illumina, USA), and index codes were added to classify the sequences for each sample. Whole-genome sequencing was performed on an Illumina NovaSeq PE150 platform using a 2 × 150 bp paired-end strategy.
The genome was annotated using Prokka V1.14.6, and resistance, virulence genes, and plasmid replicons were predicted using Abricate V1.0.1. Capsular serotypes were predicted using PathogenWatch (https://pathogen.watch). Multilocus sequence typing (MLST) of all strains was conducted using the Pasteur database (https://bigsdb.pasteur.fr/klebsiella/), and newly identified sequence typing (ST) was submitted to the MLST database administrator for approval and were assigned ST numbers. A minimum spanning tree based on the allelic difference between isolates of the seven housekeeping genes was constructed using PHYLOViZ (Ribeiro-Gonçalves et al., 2016).
2.4 Phylogenetic analysis
Single-nucleotide polymorphisms were extracted using Snippy v4.6.0 (https://github.com/tseemann/snippy) to generate a core genome alignment. This core genome alignment was employed to construct a maximum likelihood phylogenetic tree using FastTree V2.1.11-2, with KP subsp. pneumoniae HS11286 (GCA_000240185.2) as the reference genome. The resulting phylogenetic tree was visualized using iTOL (https://itol.embl.de/) (Letunic and Bork, 2024).
2.5 Galleria mellonella killing assay
Bacterial suspensions were adjusted to a 0.5 McFarland standard (approximately 1.5 × 108 CFU/mL) and then serially diluted in phosphate-buffered saline (PBS) to obtain target concentrations. A 10 µL volume of each dilution was injected into the hemocoel of Galleria mellonella larvae via the last right proleg, resulting in final inocula of approximately 1.5 × 105 and 1.5 × 104 CFU/larva (corresponding to 10-fold and 100-fold dilutions of the standard suspension, respectively). The control group received an equal volume of sterile PBS. After injection, larvae were incubated at 37°C, and survival was monitored every 2 hours for 48 hours. The mortality rate was calculated at the endpoint of the observation period.
2.6 Antibiotic-resistant plasmids and mobile genetic elements containing the blaKPC-2 gene
The assembled contigs of six CR-hvKP strains resistant to cefotaxime/avibactam were input into the VRprofile2 pipeline (https://tool2-mml.sjtu.edu.cn/VRprofile/home.php) to divide the contigs into chromosome or plasmid fragments and identify whether they carry resistance genes. To further identify plasmids carrying resistance genes, we compared these contigs with the NCBI database (https://www.ncbi.nlm.nih.gov/) using BLASTn to search for reference sequences. For the plasmid resistant to cefotaxime/avibactam, we used pC76 KPC (NZ_CP080299.1) with a coverage of 93.82% and an identity of 99.37% as a reference sequence to study the structure of the plasmid. MAUVE (Darling et al., 2004) was used to identify all contigs located in drug-resistant plasmids by comparing the contigs of the strain with the reference plasmid pC76 KPC (NZ_CP080299.1). Plasmid maps were presented using BRIG (Alikhan et al., 2011). To determine whether contigs carrying the blaKPC-2 gene all contain insertion and repeat sequences, we submitted them to Isfinder (Siguier et al., 2006) (https://www-is.biotoul.fr/index.php). Easyfig (Sullivan et al., 2011) was used to visualize the structure of the blaKPC-2 gene.
2.7 Statistical methods
Categorical variables were compared between the CR-hvKP and carbapenem-resistant non-hypervirulent Klebsiella pneumoniae (CR-non-hvKP) groups using the Chi-square test, and presented as numbers (percentages). The Log-rank (Mantel-Cox) test was used to analyze survival differences in the Galleria mellonella infection model. Differences in AST results among groups were analyzed using the Wilcoxon rank-sum and Kruskal–Wallis tests, as appropriate. Spearman’s correlation analysis was conducted to explore the relationships between plasmids and virulence genes, plasmids and resistance genes, resistance genes and clinical data, and virulence genes and clinical data.
Univariate logistic regression analyses were performed using IBM SPSS Statistics for Windows, version 27.0 (IBM Corp., Armonk, N.Y., USA) to identify potential risk factors for CR-hvKP infection. Variables with a p-value < 0.05 in the univariate analysis were included in the multivariate model. The multivariate logistic regression analysis was conducted using GraphPad Prism (version 8). The goodness-of-fit of the final multivariate model was assessed using the Hosmer–Lemeshow test, and its discriminatory power was evaluated through the area under the receiver operating characteristic curve. A two-tailed p-value < 0.05 was considered statistically significant.
Statistical analyses were performed using IBM SPSS Statistics for Windows, version 27.0 (IBM Corp., Armonk, N.Y., USA) and GraphPad Prism 8 (GraphPad Software, LaJolla, CA, USA). Visualization was performed using Xiantao Academic (https://www.xiantaozi.com/), the Wekemo Bioincloud (https://www.bioincloud.tech) (Gao et al., 2024), and Origin 2021 (OriginLab Co., MA, USA). Plots were generated using R software (v.4.2.2) with the ggplot2 (v.3.3.6) (Villanueva and Chen, 2019) package through Hiplot Pro (https://hiplot.com.cn/).
2.8 Ethical statements
The study was approved by the Ethics Committee of the First Affiliated Hospital of Hebei North University (ethical approval No. K2019147), which waived the requirement of written informed consent from patients. All strains are part of the routine laboratory procedures of the hospital and do not involve any human genetic resources. This study was conducted in accordance with the principles outlined in the Declaration of Helsinki.
3 Results
3.1 Molecular epidemiology
The global prevalence of CRKP is highlighted in Figure 2A, with China, the United States, and Brazil being the most affected countries. Within China, the highest frequency of CRKP isolates was reported in Beijing, followed by Shanghai and Hunan (Figure 2B). The molecular epidemiology at our hospital reflected this national trend (Figure 3A), with ST11 being the predominant type (69.14%), followed by ST15 (11.11%) among the 12 STs detected (Figure 3B). Capsular (K) antigen serotyping revealed 16 different types, with K64 being the most common (60.49%, 49/81). Similarly, lipopolysaccharide (O) antigen serotyping identified O1/O2v1 as the predominant serotype (66.67%, 54/81). These strains were primarily isolated from sputum (65.43%, 53/81), followed by urine (16.05%, 13/81) and secretions (4.94%, 4/81) (Figure 3C), which is consistent with provincial surveillance data (Wang et al., 2023, 2022). Phylogenetic analysis showed that evolutionary clustering was associated with ST type rather than specimen source. The primary clusters identified were ST11-K64-O1/O2v1, ST11-K47-OL101, and ST15-K19-O1/O2v2. Notably, ST11 strains diverged into three distinct branches based on serotype (ST11-K64-O1/O2v1, ST11-K47-OL101, and ST11-K25-other), suggesting divergent evolution within this lineage (Figure 4). The overall distributions of ST, K-type, O-type, and sample sources among the 81 CRKP strains are summarized in Figure 5.
Figure 2. Global and Chinese distribution of CRKP infection by country. (A) Global distribution of CRKP infection. (B) Chinese distribution of CRKP infection. CRKP, carbapenem-resistant Klebsiella pneumonia.
Figure 3. ST-type minimum-spanning tree. (A) Genetic relationship of carbapenem-resistant Klebsiella pneumoniae isolates of China. (B) Genetic relationship of carbapenem-resistant Klebsiella pneumonia isolates of this study; (C) Genetic relationship between MLST and K-type of carbapenem-resistant Klebsiella pneumonia isolates of this study. The size of the dots is proportional to the number of strains. Numbers on the connecting line represent genetic distance. MLST, MultiLocus Sequence Typing.
Figure 4. Evolutionary relationship diagram of 81 CRKP strains. CRKP: carbapenem-resistant Klebsiella pneumonia.
Figure 5. Overall distribution of 81 CRKP strains. (A) ST distribution; (B) K-type distribution; (C) O-type distribution; (D) Sample distribution CRKP: carbapenem-resistant Klebsiella pneumonia; ST, sequence typing.
3.2 Multivariate logistic regression analysis of independent risk factors for CR-hvKP infection
For the convenience of analysis, we categorized 81 CRKP strains into the following two groups: the high virulence (rmpA, rmpA2, iucA, iroB, and peg-344) (Russo et al., 2018) and low virulence groups. We collected clinical data from 81 patients with CRKP, including patient Basic Data, Department, Underlying diseases, Infection type, Invasive procedures and devices, Antimicrobial exposure, and Outcomes. After analysis, significant differences (p<0.05) were found between the two groups of patients regarding Antimicrobial usage time, Pulmonary disease, Malignant tumors, and Carbapenem antimicrobial exposure and Outcomes. Logistic multiple regression analysis was conducted on the items with significant differences mentioned above, and it was found that Antibiotic usage time, Carbapenem antibiotic exposure, and Malignant tumors were independent risk factors for CR-hvKP infection (Tables 1, 2).
Table 1. Clinical characteristics of CR-hvKP and CR-non-hvKP.
Table 2. Multivariate logistic regression analysis of independent risk factors for CR-hvKP infection.
3.3 Antimicrobial susceptibility tests
We conducted a new antibiotics combination β-lactam/β-lactase inhibitor, tigecycline and polymyxin susceptibility test on 81 CRKP strains, including Aztreonam/Avibactam, Ceftaroline/Avibactam, CZA, Imipenem/Avibactam, Meropenem/Avibactam, and Ertapenem/Avibactam, and classified them into high (4 μg/mL) and low (8 μg/mL) concentration inhibitor groups. Except for Ceftaroline, significant differences (p<0.05) were found in the drug sensitivity results of the other five antibiotics after adding low-concentration avibactam inhibitors (Figure 6C). The MIC values of CZA, Meropenem/Avibactam, and Ertapenem/Avibactam increased as the inhibitor concentration rose (p<0.05) (Figures 6B, E, F). However, no significant change was observed in MIC values after adding high-concentration inhibitors, including Aztreonam/Avibactam and Imipenem/Avibactam (Figures 6A, D). The resistance rate of polymyxin (79.01%) was lower than that of tigecycline (96.30%), and the resistance rate and MIC range values of the CR-hvKP group were generally higher than those of the CR-non-hvKP group (Tables 3, 4).
Figure 6. Comparison of MIC values of CRKP strains against novel antibiotics β-lactam/β-lactase inhibitors. (A) MIC values of Aztreonam, Aztreonam/Avibactam-4, and Aztreonam/Avibactam-8; (B) MIC values of Ceftazidim, CZA-4, and CZA-8; (C) MIC values of Ceftaroline, Ceftaroline/Avibactam-4, and Ceftaroline/Avibactam-8; (D) MIC values of Imipenem, Imipenem/Avibactam-4, and Imipenem/Avibactam-8; (E) MIC values of Meropenem, Meropenem/Avibactam-4, and Meropenem/Avibactam-8; (F) MIC values of Ertapenem, Ertapenem/Avibactam-4, and Ertapenem/Avibactam-8. ns represents not significant (p ≥ 0.05), **represents p < 0.01, and ***represents P < 0.001. CRKP: carbapenem-resistant Klebsiella pneumoniae; MIC: minimum inhibitory concentration; ST, sequence typing; Azt, Aztreonam; Ceftazi, Ceftazidime; Ceftaro,Ceftaroline; Imi, Imipenem; Mero, Meropenem; Eta, Ertapenem; Avi, Avibactam; CZA, Ceftazidime-avibactam.
Table 3. Drug susceptibility characteristics of novel enzyme inhibitor antibiotics.
Table 4. Resistance rates of CR-hvKP and CR-non-hvKP.
3.4 Antimicrobial resistance genes, virulence genes, and plasmids of CR-hvKP
All 52 CR-hvKP strains carried the carbapenemase gene blaKPC-2. Carbapenem resistance genes included blaKPC-2 and blaOXA-1. Chloramphenicol resistance genes were found only in the ST11-K47-OL101 strain. The blaCTX-M-15 gene was detected only in ST307-K102-O1/O2v2, while SHV-182 was found in both ST307-K102-O1/O2v2 and ST1-K64-O1/O2v1. All strains carry blaKPC-2, ompk35, and ompk36 genes (Figure 7, Figure 8A). Carbapenem resistance genes blaOXA-1 and blaCTX-M-15 were strongly correlated with IncFII_1_pKP91, while tetracycline resistance genes were strongly associated with IncFII(pCRY) (Figure 9B). The resistance genes of macrolide antibiotics were positively correlated with CRP (Figure 9D).
Figure 7. Evolutionary relationship diagram of 52 CR-hvKP strains. CR-hvKP: carbapenem-resistant hypervirulent Klebsiella pneumonia.
Figure 8. The characteristic distribution of resistance and virulence genes. (A) Distribution of resistance genes; (B) Distribution of virulence genes.
Figure 9. Correlations among virulence genes, resistance genes, plasmid replicons, and clinical data of the CR-hvKP strain. (A) Spearman correlation between virulence genes and plasmid replicon. (B) Spearman correlation between resistance genes and plasmid replicon. (C) Spearman correlation between virulence genes and clinical data. (D) Spearman correlation between resistance genes and clinical data. "*" represents the size of the P value, ns represents p≥0.05, *represents p<0.05, **represents p<0.01, and ***represents p<0.001. CR-hvKP: carbapenem-resistant hypervirulent Klebsiella pneumonia.
Virulence gene analysis revealed that 28 CR-hvKP strains carried all tested virulence genes. The lowest detection rates were observed for iroB (61.54%) (Figure 8B). Siderophore virulence genes showed high carrying rates in all strains. Salmonella was detected in 61.54% of the strains, which was lower than that for aerobactin (92.31%), yersiniabactin (98.08%), and enterobactin (100%). The yersiniabactin gene was not detected in ST307-K102-O1/O2v2. No strains tested positive for peg-344 (metabolite transporter) virulence genes. Among the 52 CR-hvKP strains, 10 plasmid replicon types were identified (Figure 7). Aerobactin and yersiniabactin virulence genes were strongly correlated with the IncR_1 (Figure 9A). The virulence genes of Salmochelin were positively correlated with CRP (Figure 9C).
The experiment on the Galleria mellonella showed that the survival rate of CR-non-hvKP was higher than that of CR- hvKP. The survival rate of CR-hvKP-47 was higher than that of CR-hvKP-64, indicating that the virulence of CR-hvKP-64 was higher than that of CR-hvKP-47 (Figure 10).
Figure 10. The Galleria mellonella infection model. (A) Survival rates of Galleria mellonella infected (1:10) with CR-hvKP and CR-non-hvKP; (B) Survival rates of Galleria mellonella infected (1:100) with CR-hvKP and CR-non-hvKP; (C) Survival rates of Galleria mellonella infected (1:10) with CR-hvKP-K47 and CR-hvKP-K64; (D) Survival rates of Galleria mellonella infected (1:100) with CR-hvKP-K47 and CR-hvKP-K64 CR-hvKP: carbapenem-resistant hypervirulent Klebsiella pneumonia CR-non-hvKP: carbapenem-resistant non-hypervirulent Klebsiella pneumonia. ** indicates p < 0.01, and *** indicates p < 0.001.
3.5 Antibiotic resistance plasmids and mobile genetic elements
Blast analysis showed that these drug-resistant plasmids were highly similar to the reference plasmid pC76 KPC (NZ_CP080299.1) with 99.37% identity and 93.25% coverage (Figure 11). This plasmid is located on IncFII/IncR. pKpn30, pKpn45, and pKpn46 were similar in structure and belonged to the plasmid replicons of IncFIB(AP001918)_1, IncFIC(FII)_1, and IncFII(pHN7A8)_1. The structures of pKpn32, pKpn38, and pKpn70 were similar, belonging to repB_KLEB_VIR, IncHI1B(pNDM-MAR)_1, and IncFII(pHN7A8)_1 plasmid replicons. After analysis, the blaKPC-2 genes of the six CR-hvKP strains were all located on plasmids and carried blaTEM-1B, blaSHV-12, and rmtB_1 resistance genes. Research on the genetic environment of blaKPC-2 has shown that its horizontal transfer is mediated by TnpR_Tn3-ISKpn27-blaKPC-2-ISKpn6. Various transposons and insertion sequences around the resistance genes predicted the possibility of horizontal transfer of the resistance genes. ISKpn27 contained IRL and IRR, while ISKpn6 only had IRR (Figure 12).
Figure 11. Comparative genomic analysis of antibiotic-resistant plasmids in CR-hvKP strains. From the inner circle to the outer circle, use this as follows: circle 1, GC Content; circle 2, GC Skew; circle 3, reference plasmid pC76 KPC(NZ_CP080299.1); circle 4, pKpn30; circle 5, pKpn45; circle 6, pKpn46; circle 7, pKpn32; circle 8, pKpn38; circle 9, pKpn70; and circle 10, Annotated mobile genetic elements and carbapenem resistance genes. CR-hvKP: carbapenem-resistant hypervirulent Klebsiella pneumonia.
Figure 12. Linear alignment of gene environment in blaKPC-2.
4 Discussion
The global prevalence of CR-hvKP is increasing (Han et al., 2022). Surveillance data from China indicate that the carbapenem resistance rate in K. pneumoniae has risen approximately ninefold, from 2.95% in 2005 to 25.4% in 2023, making it the second most common clinical isolate after Escherichia coli (https://www.chinets.com/Data/AntibioticDrugFast). A similar upward trend was observed in Europe, where resistance increased from 7.1% in 2017 to 11.7% in 2021 (https://www.ecdc.europa.eu/). CR-hvKP combines high pathogenicity with multidrug resistance, frequently resulting in poor clinical outcomes. Its capacity for rapid spread and outbreak formation poses significant challenges for infection control. Known risk factors for CR-hvKP infection include underlying chronic comorbidities, diabetes, age <65 years, and mechanical ventilation (Liang et al., 2022; Li et al., 2023).
In our hospital, the predominant sequence type was ST11, which aligns with findings from previous studies in China (Pu et al., 2023a). This stands in contrast to the epidemic ST258 clone prevalent in Europe and the United States (David et al., 2019; Ernst et al., 2020). The global epidemiology of CRKP is diverse, with ST101 reported as dominant in Turkey (Ibik et al., 2023), ST23 and ST25 linked to hypervirulence in Madagascar (Rakotondrasoa et al., 2020), and ST25 identified in Argentinean hypermucoviscous strains (Vargas et al., 2019). In Greece, ST39 carrying blaKPC-2 is a high-risk clone causing bloodstream infections (Tryfinopoulou et al., 2023). Geographically, the Class B carbapenemase NDM is endemic in regions including Egypt, India, Pakistan, and Serbia (Wu et al., 2019), while the Class D enzyme OXA-48 frequently cooperates with other mechanisms to enhance resistance. Notably, the ST15 clone carrying blaOXA-232 has spread globally from a likely origin in the United States, with significant prevalence in China (Wu et al., 2023). The emergence of CRKP strains producing multiple carbapenemases is a growing concern (Gao et al., 2020). We identified a high-risk ST307 clone lacking the siderophore genes yersiniabactin, aerobactin, and salmochelin. The presence of ST307 is concerning, as it has been reported that this sequence type can develop elevated MICs to CZA, leading to clinical resistance (Hernández-García et al., 2022). Furthermore, surveillance data indicate that in some regions, ST307 has replaced established high-risk clones such as ST512 and ST258 (Peirano et al., 2020). These findings from other settings necessitate heightened vigilance and monitoring for the emergence of this clone in our region.
K. pneumoniae expresses over 79 K-capsule types. In China, K64 (50.4%) and K47 (25.9%) are the most common among CRKP, with K64 increasing in eastern and central China while K47 declines, now primarily found in the north and northeast (Hu et al., 2024). Capsule-deficient strains exhibit impaired transmission and phagocytosis resistance but may display enhanced antibiotic resistance, causing persistent urinary tract infections. Conversely, hypercapsular strains show increased resistance to phagocytosis, dissemination, and mortality (Hu et al., 2024; Ernst et al., 2020). The evolving serotype prevalence underscores the urgent need for CRKP vaccines. Our finding that CR-hvKP-K64 is more virulent than CR-hvKP-K47 is consistent with that of the report by Jia et al (Jia et al., 2021).
Avibactam is a non-β-lactam β-lactamase inhibitor active against serine β-lactamases but not metallo-β-lactamases. Due to its efficacy and safety against serine carbapenemases, CZA has become a first-line treatment for CRKP (Tumbarello et al., 2021). However, CZA resistance is emerging, driven by blaKPC mutations, including blaKPC-135 and blaKPC-112, frequently selected by antibiotic pressure (Shi et al., 2024; Shen et al., 2022). The spread of blaKPC involves multiple mechanisms, including mobile genetic elements, plasmid transfer, and clonal spread. In our study, CZA resistance was primarily mediated by a blaKPC-2-bearing plasmid, pC76 KPC (NZ_CP080299.1), a non-classical type first cataloged in the RefSeq database in the United States (Tatusova et al., 2014). The blaKPC gene is typically mobilized by elements, including the Tn4401 transposon (common in the USA) and the NTEKPC-I/II elements (common in China and Brazil) (Naas et al., 2012). Our analysis of the blaKPC-2 genetic context revealed flanking insertion sequences ISKpn27 and ISKpn6. The presence of incomplete reverse repeats adjacent to ISKpn27 may facilitate gene mobilization. Different blaKPC-2 subtypes are associated with distinct mobile elements and plasmids, highlighting the complexity of their dissemination.
Beyond CZA resistance, we observed that polymyxin resistance rates were lower than those for tigecycline. However, hypervirulent ST11-K64 can rapidly develop resistance during tigecycline or polymyxin therapy (Jin et al., 2021). The emergence of the tmexCD-toprJ gene complex further threatens tigecycline efficacy (Yao et al., 2024). Concurrently, increasing reports of polymyxin-resistant, hypervirulent K. pneumoniae pose a major public health challenge (Liu et al., 2022). Given these multifaceted resistance threats, novel therapeutic approaches, including non-ribosomal tobramycin-cyclohexane conjugates, could enhance the efficacy of β-lactam/β-lactamase inhibitor combinations (Idowu et al., 2019), offering a promising direction for future research.
This study had some limitations that should be considered. First, the single-center design may mean that the observed clonal composition is specific to our local epidemiology and antibiotic practices. Second, while providing detailed insights, the sample size may not capture the full diversity of CR-hvKP in the region. These factors could limit the generalizability of our results. To overcome these constraints, future research will focus on employing long-read sequencing for complete genomic analysis and expanding into a multi-center study to verify the national relevance of our findings.
5 Conclusion
In conclusion, our study underscores the pivotal role of the ST11-K64 clone in converging CR-hvKP in China, linked clinically to prolonged antibiotic exposure and malignancies. We delineated a novel plasmid-borne genetic context (TnpR_Tn3-ISKpn27-blaKPC-2-ISKpn6) responsible for high-level CZA resistance. These findings highlight an urgent need for enhanced surveillance and stringent infection control to contain this public health threat. Future efforts should focus on developing rapid diagnostics for these resistance mechanisms and exploring novel combination therapies.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, BioProject number:PRJNA1200658.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
NW: Writing – review & editing, Methodology, Writing – original draft, Software, Investigation, Visualization. LD: Methodology, Investigation, Writing – review & editing. HL: Visualization, Data curation, Writing – review & editing. NJ: Investigation, Writing – review & editing, Software. XP: Validation, Investigation, Writing – review & editing. JC: Software, Investigation, Writing – review & editing. JH: Software, Writing – review & editing, Validation. JT: Writing – review & editing, Data curation. CL: Formal analysis, Writing – review & editing. BW: Funding acquisition, Writing – review & editing. JL: Project administration, Supervision, Writing – review & editing. WZ: Methodology, Writing – original draft, Resources, Writing – review & editing, Conceptualization, Project administration, Supervision.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was supported by the Scientific Research Fund of Hebei Health Commission (No.20231461).
Acknowledgments
We thank all colleagues and classmates in the laboratory for their hard work during the experiment. We also thank colleagues in the Microbiology Department of the First Affiliated Hospital of Hebei North University for their technical support. We would like to express our gratitude to Minghua Zhan, Baoliang Li, Teng Wang, Shuting Liu, and Caiqing Li for their contributions to the manuscript. The author thanks Editage (www.editage.cn) for English language editing.
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alikhan, N. F., Petty, N. K., Ben Zakour, N. L., and Beatson, S. A. (2011). BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12, 402. doi: 10.1186/1471-2164-12-402, PMID: 21824423
Darling, A. C., Mau, B., Blattner, F. R., and Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704, PMID: 15231754
David, S., Reuter, S., Harris, S. R., Glasner, C., Feltwell, T., Argimon, S., et al. (2019). Epidemic of carbapenem-resistant Klebsiella pneumoniae in Europe is driven by nosocomial spread. Nat. Microbiol. 4, 1919–1929. doi: 10.1038/s41564-019-0492-8, PMID: 31358985
Ernst, C. M., Braxton, J. R., Rodriguez-Osorio, C. A., Zagieboylo, A. P., Li, L., Pironti, A., et al. (2020). Adaptive evolution of virulence and persistence in carbapenem-resistant Klebsiella pneumoniae. Nat. Med. 26, 705–711. doi: 10.1038/s41591-020-0825-4, PMID: 32284589
Gao, H., Liu, Y., Wang, R., Wang, Q., Jin, L., and Wang, H. (2020). The transferability and evolution of NDM-1 and KPC-2 co-producing Klebsiella pneumoniae from clinical settings. EBioMedicine. 51, 102599. doi: 10.1016/j.ebiom.2019.102599, PMID: 31911273
Gao, Y., Zhang, G., Jiang, S., and Liu, Y. X. (2024). Wekemo Bioincloud: A user-friendly platform for meta-omics data analyses. Imeta. 3, e175. doi: 10.1002/imt2.175, PMID: 38868508
Gu, D., Dong, N., Zheng, Z., Lin, D., Huang, M., Wang, L., et al. (2018). A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: a molecular epidemiological study. Lancet Infect. Dis. 18, 37–46. doi: 10.1016/s1473-3099(17)30489-9, PMID: 28864030
Han, Y. L., Wen, X. H., Zhao, W., Cao, X. S., Wen, J. X., Wang, J. R., et al. (2022). Epidemiological characteristics and molecular evolution mechanisms of carbapenem-resistant hypervirulent Klebsiella pneumoniae. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.1003783, PMID: 36188002
Hernández-García, M., Castillo-Polo, J. A., Cordero, D. G., Pérez-Viso, B., García-Castillo, M., Saez de la Fuente, J., et al. (2022). Impact of ceftazidime-avibactam treatment in the emergence of novel KPC variants in the ST307-klebsiella pneumoniae high-risk clone and consequences for their routine detection. J. Clin. Microbiol. 60, e0224521. doi: 10.1128/jcm.02245-21, PMID: 35107303
Hu, F., Pan, Y., Li, H., Han, R., Liu, X., Ma, R., et al. (2024). Carbapenem-resistant Klebsiella pneumoniae capsular types, antibiotic resistance and virulence factors in China: a longitudinal, multi-centre study. Nat. Microbiol. 9, 814–829. doi: 10.1038/s41564-024-01612-1, PMID: 38424289
Ibik, Y. E., Ejder, N., Sevim, E., Rakici, E., Tanriverdi, E. S., and Copur Cicek, A. (2023). Evaluating molecular epidemiology of carbapenem non-susceptible Klebsiella pneumoniae isolates with MLST, MALDI-TOF MS, PFGE. Ann. Clin. Microbiol. Antimicrob. 22, 93. doi: 10.1186/s12941-023-00640-9, PMID: 37891628
Idowu, T., Ammeter, D., Arthur, G., Zhanel, G. G., and Schweizer, F. (2019). Potentiation of β-lactam antibiotics and β-lactam/β-lactamase inhibitor combinations against MDR and XDR Pseudomonas aeruginosa using non-ribosomal tobramycin-cyclam conjugates. J. Antimicrob. Chemother. 74, 2640–2648. doi: 10.1093/jac/dkz228, PMID: 31139830
Jia, X., Li, C., Chen, F., Li, X., Jia, P., Zhu, Y., et al. (2021). Genomic epidemiology study of Klebsiella pneumoniae causing bloodstream infections in China. Clin. Transl. Med. 11, e624. doi: 10.1002/ctm2.624, PMID: 34841702
Jin, X., Chen, Q., Shen, F., Jiang, Y., Wu, X., Hua, X., et al. (2021). Resistance evolution of hypervirulent carbapenem-resistant Klebsiella pneumoniae ST11 during treatment with tigecycline and polymyxin. Emerg. Microbes Infect. 10, 1129–1136. doi: 10.1080/22221751.2021.1937327, PMID: 34074225
Letunic, I. and Bork, P. (2024). Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 52, W78–w82. doi: 10.1093/nar/gkae268, PMID: 38613393
Li, L., Li, S., Wei, X., Lu, Z., Qin, X., and Li, M. (2023). Infection with Carbapenem-resistant Hypervirulent Klebsiella Pneumoniae: clinical, virulence and molecular epidemiological characteristics. Antimicrob. Resist. Infect. Control. 12, 124. doi: 10.1186/s13756-023-01331-y, PMID: 37953357
Liang, S., Cao, H., Ying, F., and Zhang, C. (2022). Report of a fatal purulent pericarditis case caused by ST11-K64 carbapenem-resistant hypervirulent klebsiella pneumoniae. Infect. Drug Resist. 15, 4749–4757. doi: 10.2147/idr.S379654, PMID: 36034175
Liu, X., Wu, Y., Zhu, Y., Jia, P., Li, X., Jia, X., et al. (2022). Emergence of colistin-resistant hypervirulent Klebsiella pneumoniae (CoR-HvKp) in China. Emerg. Microbes Infect. 11, 648–661. doi: 10.1080/22221751.2022.2036078, PMID: 35086435
Naas, T., Cuzon, G., Truong, H. V., and Nordmann, P. (2012). Role of ISKpn7 and deletions in blaKPC gene expression. Antimicrob. Agents Chemother. 56, 4753–4759. doi: 10.1128/aac.00334-12, PMID: 22733068
Peirano, G., Chen, L., Kreiswirth, B. N., and Pitout, J. D. D. (2020). Emerging antimicrobial-resistant high-risk klebsiella pneumoniae clones ST307 and ST147. Antimicrob. Agents Chemother. 64, e01148-20. doi: 10.1128/aac.01148-20, PMID: 32747358
Pu, D., Zhao, J., Chang, K., Zhuo, X., and Cao, B. (2023a). Superbugs" with hypervirulence and carbapenem resistance in Klebsiella pneumoniae: the rise of such emerging nosocomial pathogens in China. Sci. Bull. (Beijing). 68, 2658–2670. doi: 10.1016/j.scib.2023.09.040, PMID: 37821268
Pu, D., Zhao, J., Lu, B., Zhang, Y., Wu, Y., Li, Z., et al. (2023b). Within-host resistance evolution of a fatal ST11 hypervirulent carbapenem-resistant Klebsiella pneumoniae. Int. J. Antimicrob. Agents. 61, 106747. doi: 10.1016/j.ijantimicag.2023.106747, PMID: 36758779
Rakotondrasoa, A., Passet, V., Herindrainy, P., Garin, B., Kermorvant-Duchemin, E., Delarocque-Astagneau, E., et al. (2020). Characterization of Klebsiella pneumoniae isolates from a mother-child cohort in Madagascar. J. Antimicrob. Chemother. 75, 1736–1746. doi: 10.1093/jac/dkaa107, PMID: 32303060
Ribeiro-Gonçalves, B., Francisco, A. P., Vaz, C., Ramirez, M., and Carriço, J. A. (2016). PHYLOViZ Online: web-based tool for visualization, phylogenetic inference, analysis and sharing of minimum spanning trees. Nucleic Acids Res. 44, W246–W251. doi: 10.1093/nar/gkw359, PMID: 27131357
Russo, T. A., Olson, R., Fang, C. T., Stoesser, N., Miller, M., MacDonald, U., et al. (2018). Identification of Biomarkers for Differentiation of Hypervirulent Klebsiella pneumoniae from Classical K. pneumoniae. J. Clin. Microbiol. 56, e00776-18. doi: 10.1128/jcm.00776-18, PMID: 29925642
Shen, S., Tang, C., Ding, L., Han, R., Yin, D., Yang, W., et al. (2022). Identification of KPC-112 from an ST15 klebsiella pneumoniae strain conferring resistance to ceftazidime-avibactam. mSphere. 7, e0048722. doi: 10.1128/msphere.00487-22, PMID: 36374086
Shi, Q., Shen, S., Tang, C., Ding, L., Guo, Y., Yang, Y., et al. (2024). Molecular mechanisms responsible KPC-135-mediated resistance to ceftazidime-avibactam in ST11-K47 hypervirulent Klebsiella pneumoniae. Emerg. Microbes Infect. 13, 2361007. doi: 10.1080/22221751.2024.2361007, PMID: 38801099
Siguier, P., Perochon, J., Lestrade, L., Mahillon, J., and Chandler, M. (2006). ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 34, D32–D36. doi: 10.1093/nar/gkj014, PMID: 16381877
Sullivan, M. J., Petty, N. K., and Beatson, S. A. (2011). Easyfig: a genome comparison visualizer. Bioinformatics. 27, 1009–1010. doi: 10.1093/bioinformatics/btr039, PMID: 21278367
Tatusova, T., Ciufo, S., Fedorov, B., O'Neill, K., and Tolstoy, I. (2014). RefSeq microbial genomes database: new representation and annotation strategy. Nucleic Acids Res. 42, D553–D559. doi: 10.1093/nar/gkt1274, PMID: 24316578
Tryfinopoulou, K., Linkevicius, M., Pappa, O., Alm, E., Karadimas, K., Svartström, O., et al. (2023). Emergence and persistent spread of carbapenemase-producing Klebsiella pneumoniae high-risk clones in Greek hospitals 2013 to 2022. Euro Surveill. 28, 2300571. doi: 10.2807/1560-7917.Es.2023.28.47.2300571, PMID: 37997662
Tumbarello, M., Raffaelli, F., Giannella, M., Mantengoli, E., Mularoni, A., Venditti, M., et al. (2021). Ceftazidime-Avibactam Use for Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae Infections: A Retrospective Observational Multicenter Study. Clin. Infect. Dis. 73, 1664–1676. doi: 10.1093/cid/ciab176, PMID: 33618353
Vargas, J. M., Moreno Mochi, M. P., Nuñez, J. M., Cáceres, M., Mochi, S., Del Campo Moreno, R., et al. (2019). Virulence factors and clinical patterns of multiple-clone hypermucoviscous KPC-2 producing K. pneumoniae. Heliyon. 5, e01829. doi: 10.1016/j.heliyon.2019.e01829, PMID: 31286076
Villanueva, R. A. M. and Chen, Z. J. (2019). ggplot2: elegant graphics for data analysis (Boca Raton, FL, USA: Taylor & Francis).
Wang, N., Zhan, M., Liu, J., Wang, Y., Hou, Y., Li, C., et al. (2022). Prevalence of carbapenem-resistant klebsiella pneumoniae infection in a northern province in China: clinical characteristics, drug resistance, and geographic distribution. Infect. Drug Resist. 15, 569–579. doi: 10.2147/idr.S347343, PMID: 35228807
Wang, N., Zhan, M., Wang, T., Liu, J., Li, C., Li, B., et al. (2023). Long term characteristics of clinical distribution and resistance trends of carbapenem-resistant and extended-spectrum β-lactamase klebsiella pneumoniae infections: 2014-2022. Infect. Drug Resist. 16, 1279–1295. doi: 10.2147/idr.S401807, PMID: 36910515
Wu, W., Feng, Y., Tang, G., Qiao, F., McNally, A., and Zong, Z. (2019). NDM metallo-β-lactamases and their bacterial producers in health care settings. Clin. Microbiol. Rev. 32, e00115-18. doi: 10.1128/cmr.00115-18, PMID: 30700432
Wu, Y., Jiang, T., He, X., Shao, J., Wu, C., Mao, W., et al. (2023). Global Phylogeography and Genomic Epidemiology of Carbapenem-Resistant bla(OXA-232)-Carrying Klebsiella pneumoniae Sequence Type 15 Lineage. Emerg. Infect. Dis. 29, 2246–2256. doi: 10.3201/eid2911.230463, PMID: 37877525
Yang, X., Dong, N., Chan, E. W., Zhang, R., and Chen, S. (2021). Carbapenem resistance-encoding and virulence-encoding conjugative plasmids in klebsiella pneumoniae. Trends Microbiol. 29, 65–83. doi: 10.1016/j.tim.2020.04.012, PMID: 32448764
Yang, X., Sun, Q., Li, J., Jiang, Y., Li, Y., Lin, J., et al. (2022). Molecular epidemiology of carbapenem-resistant hypervirulent Klebsiella pneumoniae in China. Emerg. Microbes Infect. 11, 841–849. doi: 10.1080/22221751.2022.2049458, PMID: 35236251
Yao, H., Zhang, T., Peng, K., Peng, J., Liu, X., Xia, Z., et al. (2024). Conjugative plasmids facilitate the transmission of tmexCD2-toprJ2 among carbapenem-resistant Klebsiella pneumoniae. Sci. Total Environ. 906, 167373. doi: 10.1016/j.scitotenv.2023.167373, PMID: 37758131
Keywords: carbapenem-resistant hypervirulent Klebsiella pneumoniae, ceftazidime/avibactam resistance, ST11 clone, molecular epidemiology, China
Citation: Wang N, Deng L, Li H, Jia N, Peng X, Chang J, Hao J, Tang J, Lei C, Wang B, Liu J and Zhang W (2025) Molecular epidemiology of carbapenem-resistant hypervirulent Klebsiella pneumoniae: risk factors and resistance mechanism of ceftazidime/avibactam in China. Front. Cell. Infect. Microbiol. 15:1698033. doi: 10.3389/fcimb.2025.1698033
Received: 03 September 2025; Accepted: 26 November 2025; Revised: 24 November 2025;
Published: 19 December 2025.
Edited by:
Mariusz Dyląg, University of Warsaw, PolandReviewed by:
Miguel A. Ares, Mexican Social Security Institute, MexicoDongjie Chen, Fujian Provincial Hospital, China
Copyright © 2025 Wang, Deng, Li, Jia, Peng, Chang, Hao, Tang, Lei, Wang, Liu and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Zhang, MTUzNjkzMTgzMThAMTYzLmNvbQ==; Jianhua Liu, amlhbmh1YWxpdTE5NzlAMTYzLmNvbQ==; Bu Wang, Mjg5OTc0NzI2QHFxLmNvbQ==