Analysis of Drug Resistance, Virulence, and Phylogeny of Klebsiella pneumoniae in Patients with Different Infections

authors:

avatar Jianhua Liu # 1 , avatar Rui Li # 1 , avatar Shukun Qian 1 , avatar Jianing Xi 2 , *

Department of Clinical Laboratory, Beijing Rehabilitation Hospital, Capital Medical University, Beijing, China
Department of Respiratory Rehabilitation, Beijing Rehabilitation Hospital, Capital Medical University, Beijing, China
# These authors are contributed equally as the first author.

how to cite: Liu J, Li R, Qian S, Xi J. Analysis of Drug Resistance, Virulence, and Phylogeny of Klebsiella pneumoniae in Patients with Different Infections. Jundishapur J Microbiol. 2023;16(5):e133082. https://doi.org/10.5812/jjm-133082.

Abstract

Background:

Klebsiella pneumoniae is an important pathogen among nosocomial infections, which can cause urinary and respiratory system infections, surgical site infections, and sepsis. Recently, carbapenem-resistant K. pneumoniae showed an upward trend with the wide use of clinical carbapenem antibiotics. However, there are few studies on the relationship between drug resistance, virulence, and phylogeny of K. pneumoniae in patients with different infections.

Objectives:

We investigated the drug resistance, virulence, and phylogeny of K. pneumoniae in patients with different infections.

Methods:

Seventy infected patients were selected as subjects. The extended-spectrum β-lactamases (ESBLs) color screening plane and blood plate medium were used for culture. Identification and drug resistance analysis was carried out by VITEK-2 compact automatic bacterial analyzer. Multilocus sequence typing (MLST) and capsule genotyping analysis were also performed. The Xpert CarbaR cartridge detected the carbapenem resistance genes. Comprehensive Antibiotic Research Database (CARD) and Virulence Factor Database (VFDB) predicted the drug-resistant genes and virulence factor genes, respectively. The phylogenetic tree was constructed, and the correlation was analyzed.

Results:

A total of 43 K. pneumoniae strains were cultured. We found that all K. pneumoniae strains exhibited different multiple drug resistance. MLST analysis indicated that ST11 was the main ST (60.61%). Analysis of the carbapenem-resistance genes showed that all isolates harbored the blaKPC-2 gene and some others blaOXA. The prediction result of capsular blood genotyping and virulence factor genes indicated that K47, K64, and K25 were the main types of capsular blood, and the top three detection rates of virulence genes were fimH (97.67%), mrkD (94.19%), and entB (84.88%). All isolates were clustered into one branch based on the virulence factor genes in the phylogenetic tree, and the strains of the same ST type or capsular blood type showed a closer relationship. Correlation analysis manifested that the drug resistance of K. pneumoniae in infected patients was positively correlated with virulence and phylogeny (r = 0.682, P = 0.000).

Conclusions:

There were complicated differences in the multidrug resistance to K. pneumoniae, resulting in high independent gene-positive rates of strains and a strong correlation with phylogeny, which can provide a reference for the selection of clinical antimicrobial drugs.

1. Background

Klebsiella pneumoniae is an important pathogen causing nosocomial infection, which can cause urinary and respiratory systems, surgical site infections, and sepsis (1). In Singapore, the mortality rates attributed to K. pneumoniae ranged from 20 to 26% (2). In China, K. pneumoniae accounted for 11.9% of pathogens isolated from cases of ventilator-associated pneumonia (VAP) and pneumonia acquired in the intensive care unit (ICU) (3). Recently, carbapenem-resistant K. pneumoniae showed an upward trend with the wide use of clinical carbapenem antibiotics (accounting for 74.36% of ICU cases) (4). The prevalence of carbapenem-resistant K. pneumoniae ranged from 0.9 to 23.6% in different provinces in China (5).

Carbapenemase is a kind of β-lactam enzyme that can hydrolyze carbapenem antibiotics, such as ertapenem, imipenem, and meropenem, and usually can hydrolyze β-lactam antibiotics, such as penicillin, β-lactam enzyme inhibitor complex preparation, cephalosporins, etc. According to the result of nationwide surveillance in China, K. pneumoniae carbapenemase (KPC) and New Delhi metallo-β-lactamase (NDM) were the main enzymes that induced drug resistance in K. pneumoniae (6). Also, KPC-2 was the main subtype of KPC in China, encoded by transferable plasmids and independent of the loss of membrane pore proteins, which are highly contagious and prone to outbreaks (7). Previous studies have shown that infection with carbapenem-resistant K. pneumoniae can increase the mortality of patients (8). It can cause multiple-site infection, disseminated abscess, and invasive infection, such as meningitis, liver abscess, osteomyelitis, etc. Also, it can cause secondary migration and dissemination, resulting in multiple-organ infections and sequelae (9). In addition to the increased drug resistance, high virulence K. pneumoniae with high contagion was more frequently reported than before.

The whole-genome analysis found specific pathogenic factors in highly virulent K. pneumoniae (i.e., the clamping membrane synthesis factor rmpA and iron carrier of K. pneumoniae in a virulence plasmid pLVKP) (10). Klebsiella pneumoniae carrying pLVKP had high toxicity and viscosity and quickly caused invasive infection and liver abscess (11). Patients infected with highly virulent K. pneumoniae with high blood migration infectivity can cause multiple-organ or tissue infections, affecting patients' life quality (12). Therefore, strengthening the study of drug resistance and virulence factor of K. pneumoniae is of great significance to improve the prognosis of patients.

Phylogenetic analysis has been widely used to study the evolutionary relationship between bacterium DNA and protein sequences. It made it possible to analyze the differences and evolutionary relationship between drug-resistant genes and virulence factor genes between different genotypes of K. pneumoniae, which will contribute to the reasonable selection of antibacterial drugs and clinical diagnosis and treatment. However, there are few studies on the relationship between drug resistance, virulence, and phylogeny of K. pneumoniae in patients with different infections (13, 14).

2. Objectives

In this study, the drug resistance, virulence, and phylogeny of K. pneumoniae in patients with different infections were studied.

3. Methods

3.1. Sample Collection

Infected patients (n = 70, marked KP01, KP02, KP03 … KP70) admitted to the hospital from January 2020 to May 2021 were selected as subjects, including 42 males and 28 females. The age ranged from 28 to 94 years, with an average of 71.51 ± 6.98 years. The course of the disease ranged from 1 to 8 days, with an average of 4.12 ± 0.79 days. There were 2 cases in the orthopedics ward, 6 cases in the respiratory medicine ward, 2 cases in the rehabilitation ward, 7 cases in the general geriatric ward, 10 cases in the geriatric intensive care unit, 2 cases in the urology ward, 5 cases in the general surgery ward, 16 cases in the neurological rehabilitation center, 2 cases in the neurology ward, 2 cases in the renal endocrinology ward, 4 cases in the digestive ward, 9 cases in the cardiovascular medicine ward, and 3 cases in the traditional Chinese medicine ward. Specimen type: One case of secretion, one case of cerebrospinal fluid, 36 cases of urine, 22 cases of sputum, and 10 cases of blood.

3.2. Inclusion and Exclusion Criteria

Inclusion criteria: (1) all infected patients admitted to the hospital initially diagnosed after admission, (2) those who completed the sample collection and drug resistance analysis, and (3) complete pathogen culture data and follow-up data available in the hospital.

Exclusion criteria: (1) patients with a mental disorder, cognitive dysfunction, or coagulation dysfunction; (2) patients with infectious diseases difficult to control, chemotherapy, or biological immunotherapy; and (3) death of patients or withdrawal from the experiment.

3.3. Sampling Instruments and Equipment

VITEK-2 Compact automatic bacterial analyzer was purchased from BioMerieux, Inc., France. Real-time fluorescence PCR systems were purchased from Bid-Rad, Inc., USA. Bacterial total RNA lysate was purchased from Sigma, Inc., USA. A bacterial genomic DNA extraction kit was purchased from Biomiga, Inc., USA. The reverse recording kit was purchased from Thermo, Inc., USA. A real-time fluorescent PCR kit was purchased from Toyobo, Inc., Japan.

3.4. Strain Collection and Antimicrobial Susceptibility Testing

Isolation and culture of pathogenic bacteria were completed in all 70 infected patients in the hospital. Secretion, cerebrospinal fluid, urine, sputum, and blood of patients were collected after admission, and extended-spectrum β-lactamase (ESBL) screening plate and blood plate culture were used to culture and isolate K. pneumoniae. Confirmation and anti-microbial susceptibility of isolated strains were conducted by VITEK-2 compact automatic bacterial analyzer. Identification of ESBL production was performed using the double disk confirmation method. Quality control strains included Escherichia coli (ATCC25922), Pseudomonas aeruginosa (ATCC27853), and Klebsiella pneumoniae (ATCC25923). The following 15 drugs were tested: amikacin, amoxicillin/clavulanic, ampicillin, ampicillin/sulbactam, piperacillin/tazobactam, aztreonam, cefepime, cefotaxime, ceftazidime, chloramphenicol, imipenem, meropenem, levofloxacin, sulfamethoxazole, and tetracycline. All quality control strains were purchased from the United States. Resistance to commonly used antibacterial agents was analyzed, and all operating procedures and judgment results were carried out according to Clinical and Laboratory Standards Institute (CLSI) standards (15, 16).

3.5. Xpert® CarbaR Testing

Organisms were assessed in pure culture for the blaNDM, blaIMP, blaKPC, blaVIM, and blaOXA-48 carbapenem-resistant genes using the Xpert CarbaR cartridge (Cepheid, Sunnyvale, CA) as instructed. Carbapenem-resistant organisms were grown on sheep blood agar plates (Hardy Diagnostics) to ensure the resistance gene retention before testing, using a 10 µg Journal Pre-proof 8 meropenem disk located in the inoculum center. Three to five colonies obtained from the zone of inhibition inner edge underwent suspending in Mueller-Hinton broth to 0.5 McFarland Standard, and the suspension (10 µL) was mixed with 5 mL of sample reagent before testing with the Xpert CarbaR cartridge (17).

3.6. Whole-Genome Sequencing

Bacterial DNA of K. pneumoniae strains was extracted using the QIAamp cador Pathogen DNA Mini kit (QIAGEN; Hilden, Germany) according to the instructions, and then the quantity and purity of the DNA were determined using NanoDrop 8000 spectrophotometer (Thermo Fisher Scientific, Waltham). The libraries were created using the VAHTSTM Universal DNA Library Prep kit for Illumina, and the genomes were sequenced using the Illumina Novaseq platform system to obtain reads as recommended by the manufacturer. The sequence of the primers was provided in Appendix 1 in Supplementary File. The library sizes had peaks centered around 300 bp. The paired-end reads were filtered using Fastp (v0.20.0) (18), removing adapter sequences and bases with a quality score (≤ 20) before assembly. Sequence reads were assembled into contigs using SPAdes (v3.41.0) (19), with contigs of ≥ 500 bp being selected,

3.7. Data Analysis

Based on the whole-genome sequencing (WGS) data, multilocus sequence typing (MLST) analysis was performed using the Public Database of Multilocus Sequence Typing (PubMLST, https://pubmlst.org/) (20), and the antibiotic-resistant genes were screened using RGI software (version 4.0.3) and the Comprehensive Antibiotic Research Database (CARD, https://card.mcmaster.ca/) (21). Amino acid sequence data was blasted in the Virulence Factor Database (VFDB, http://www.mgc.ac.cn/VFs/) (22) to detect virulence factors. Capsule types were identified based on the published wzc sequences in the previous report (23). The phylogenetic tree was constructed by IQ-TREE using the maximum likelihood algorithm (24).

4. Results

4.1. Strain Collection and Culture

Pathogenic bacteria were obtained from all kinds of specimens of 70 infected patients in the hospital, but samples KP04, KP07, KP36, KP50, KP55, and KP64 were contaminated or failed to be cultured, and subsequent analysis was not included. Finally, 43 K. pneumoniae strains were cultured and identified from 70 patients with different infections.

4.2. Antimicrobial Susceptibility Testing of K. pneumoniae in Patients with Infection in Hospital

The antimicrobial susceptibility of 43 K. pneumoniae strains was performed, and the results were judged "susceptible," "intermediate," or "resistant" based on the CLSI standards. The antimicrobial susceptibility results of the cultured pathogens are shown in Table 1 and Appendix 2 in Supplementary File. All K. pneumoniae strains were utterly resistant to 11 drugs, except amikacin, chloramphenicol, compound sulfamethoxazole, and tetracycline, accounting for 67.44%, 23.25%, 51.16%, and 48.84%, respectively, exhibiting different multidrug resistance patterns.

Table 1.

Analysis of Drug Resistance of Klebsiella pneumoniae in Patients with Infection in the Hospital (43 Strains)

Types of Antibacterial AgentsDrug-Resistant StrainsDrug Resistance Rate
Amikacin2967.44
Amoxicillin/clarvic acid43100.00
Ampicillin43100.00
Ampicillin/sulbactam43100.00
Aztreonam43100.00
Cefepime43100.00
Cefotaxime43100.00
Ceftazidime43100.00
Chloramphenicol1023.25
Imipenem43100.00
Levofloxacin43100.00
Piperacillin/tazobactam43100.00
Compound sulfamethoxazole2251.16
Tetracycline2148.84
Meropenem43100.00

4.3. MLST Analysis

The most frequently represented sequence typings (STs) were ST11 (40), ST15 (2), and the following ST types, each represented by one isolate: ST14, ST23, ST36, ST65, ST70, ST86, ST105, ST111, ST198, ST299, ST479, ST540, ST1333, ST3688, and ST4415, indicating the diversity of K. pneumoniae infection (Figure 1). UN means that MLST typing was unknown, and the unknown reason was that some housekeeping genes could not be identified or there was no corresponding typing in the database. The drug resistance of K. pneumoniae was predicted according to the sequence alignment results conducted by VFDB (Appendix 3 in Supplementary File). All strains showed complete resistance to 16 of 25 tested drugs, and only four drugs (lincosamide, oxazolidinone, pleuromutilin, and streptogramin) were effective against most strains. Also, the characteristics of drug resistance were familiar with ST typing. The phylogenetic tree was constructed based on the drug-resistant genes, and it could be seen from the evolutionary tree that strains assigned to the same ST types had familiar drug resistance.

Phylogenetic tree construction of Klebsiella pneumoniae with ST typing and schematic diagram of antibiotic. Each isolate represented ST11 (the most frequently represented ST) was marked as green, each isolate represented ST15 (the second frequently represented ST) was marked as light blue, each isolate represented unknown ST was marked purple, and the STs only represented one isolate and marked as light red.
Phylogenetic tree construction of Klebsiella pneumoniae with ST typing and schematic diagram of antibiotic. Each isolate represented ST11 (the most frequently represented ST) was marked as green, each isolate represented ST15 (the second frequently represented ST) was marked as light blue, each isolate represented unknown ST was marked purple, and the STs only represented one isolate and marked as light red.

4.4. Xpert CarbaR Test

The results of the Xpert CarbaR test are shown in Appendix 4 in Supplementary File. All strains were positive for the K. pneumoniae carbapenemase gene (blaKPC), and eight of them (KP14, KP16, KP19, KP20, KP26, KP28, KP40, and KP44) were positive for blaOXA, while the IMP, VIM, and NDM genes were not detected. Furthermore, the prediction of drug-resistant genes was further conducted based on their sequence, and the results indicated that the blaKPC belonged to the blaKPC-2 subtype, but the blaOXA was found to belong to two blaOXA subtypes (blaOXA-48 and blaOXA-1).

The detection of blaIMP, blaVIM, blaNDM, and blaKPC genes obtained by the Xpert CarbaR test was consistent with the gene prediction results, while the detection and prediction result of the blaOXA gene was inconsistent. Also, the results of antimicrobial susceptibility testing (Appendix 2 in Supplementary File) and the Xpert CarbaR test (Appendix 4 in Supplementary File) were compared with drug resistance prediction analysis (Appendix 3 in Supplementary File). The comparison results showed that 43 K. pneumoniae strains were resistant to five types of drugs, including penam, monobactam, cephalosporin, carbapenem, and fluoroquinolone, which is in complete agreement with the predictions of resistant genes. Aminoglycoside, phenicol, sulfonamide, and tetracycline showed no consistent results (Table 2).

Table 2.

Consistency of Drug Resistance Test and Drug Resistance Prediction Analysis

Types of Antibacterial AgentsSampleRatio, %
Aminoglycoside2967.44
Carbapenem43100.00
Cephalosporin43100.00
Fluoroquinolone43100.00
Monobactam43100.00
Penam43100.00
Phenicol1534.88
Sulfonamide3069.77
Tetracycline2455.81

4.5. Capsule Genotyping and Virulence Factor Prediction

The capsule genotyping results (Appendix 5 in Supplementary File) showed that K47, K64, and K25 were the top three types, accounting for 42.18%, 21.88%, and 6.25%, respectively, and were mainly derived from ST11 type isolate. Virulence factor prediction (Figure 2 and Appendix 6 in Supplementary File) showed that the top three virulence factors were fimH, mrkD, and entB, with detection rates of 97.67%, 94.19%, and 84.88%, respectively. Moreover, the composition of virulence factors was correlated with MLST typing, indicating that similar virulence factors were derived from similar ST types (Figure 2).

Virulence factor identification of 43 Klebsiella pneumoniae isolates by VFDB.
Virulence factor identification of 43 Klebsiella pneumoniae isolates by VFDB.

4.6. Phylogenetic Tree Construction and Correlation Analysis of Klebsiella pneumoniae in Patients with Infection in the Hospital

Gene sequencing was performed on representative strains, and homology retrieval was completed. The results showed that all strains were clustered into one branch in the phylogenetic tree, as shown in Figure 3A. The strains assigned to the same ST or capsule type had a closer genetic relationship. K47, K64, and K25 capsule-type strains were clustered into smaller branches, indicating that they may mutate from a common ancestor strain. Furthermore, the phylogenetic tree was constructed based on the K47, K64, and K25 capsule-type strains (Figure 3B), and the results showed that K47 and K25 capsule-type strains had a closer relationship than K64. Correlation analysis manifested that the drug resistance of K. pneumoniae in infected patients was positively correlated with virulence and phylogeny (r = 0.682, P = 0.000).

Phylogenetic tree based on the virulence factor sequence from whole-genome sequencing (WGS) data of all isolates (A) and main capsule genotypes (K47, K64, and K25) strains (B). The horizontal coordinate represented the genetic distance.
Phylogenetic tree based on the virulence factor sequence from whole-genome sequencing (WGS) data of all isolates (A) and main capsule genotypes (K47, K64, and K25) strains (B). The horizontal coordinate represented the genetic distance.

5. Discussion

Klebsiella pneumoniae is a common Gram-negative bacterium in the clinic, ranking as the second most common bacterium in bacterial infection (25). In recent years, the report of high drug resistance of K. pneumoniae increased with the clinical abuse of antibiotics. Beyrouthy et al. revealed that high drug-resistant K. pneumoniae could be isolated from the intestinal and respiratory tract (26). Shi et al. reported two cases of enhanced ceftazidime-avibactam resistance of K. pneumoniae due to the mutation of the drug-resistant gene (27). In this study, 43 K. pneumoniae strains were cultured from 70 infected patients in the hospital. The drug resistance tests obtained by the Xpert CarbaR test and predicted by the CARD database of cultured pathogens all indicated the high-level multidrug resistant ability of all the isolates. The possible reason is related to the clinical abuse of antibiotics. A further genetic study revealed that all K. pneumoniae isolates harboring blaKPC-2 and ST11 were the main STs accounting for 60.61% of isolates. These findings were consistent with previous reports (6, 28, 29). However, the detection and sequence prediction of blaOXA was inconsistent, which needs more research on other levels (transcriptome, proteome, etc.).

Klebsiella pneumoniae with high viscosity in invasive infections can effectively resist the killing effect of complement and phagocytes, reducing sensitivity to serum killing and exhibiting stronger virulence in animal experiments (30). Mucinous strains with clinical manifestations have been identified as highly virulent (31). The rmpA gene encodes the functional protein rmpA, which regulates capsular polysaccharide synthesis and is closely related to the mucus phenotype. In this study, we analyzed the capsular blood type and drug-resistance virulence genes of K. pneumoniae. The results showed that the K47, K64, and K25 types were the main capsular blood types, with a total detection rate of 70.31%. The detection rates of virulence genes were also relatively high, with fimH, mrkD, and entB having detection rates of 97.67%, 94.19%, and 84.88%, respectively. The fimH gene encoded the minor subunit FimH and was responsible for the adhesive properties of Type 1 and 3 fimbriae. Type 1 fimbriae and type 3 fimbriae were the mediators of K. pneumoniae that have been identified as important pathogenicity factors during infections (32).

Type 1 fimbriae bond D-mannosylated glycoproteins and expressed in 90% of both clinical and environmental K. pneumoniae isolates. Type 3 fimbriae were located at the tips of extracellular matrix proteins, such as type IV and V collagens, with the adhesion mrkD, which is crucial for biofilm production. Jagnow and Clegg reported that Type 1 and type 3 fimbriae were all expressed during biofilm formation on catheters and may lead to urinary tract infections. They also promote the colonization of endotracheal tubes and aggravate lung infection (33). Iron plays a significant role in bacterial growth and reproduction, and the aerobactin system is a strong factor in iron uptake. The entB gene encoded enterobactin and the enterobactin expression played an important role in the iron uptake system in K. pneumoniae (34). However, the specific mechanism behind the K47, K64, and K25 dominance as the main capsular type remains unclear, which may be related to the KMV phenotype under these genes.

5.1. Conclusions

According to the clustering result from drug resistance gene sequencing and virulence factor sequencing, we found isolates with the same ST type in closer branches that had similar drug resistance and virulence factor, indicating that drug resistance may have a correlation with virulence factor. Correlation analysis manifested that the drug resistance of K. pneumoniae in infected patients was positively correlated with virulence and phylogeny. Therefore, the drug resistance, virulence, and phylogenetic tree of K. pneumoniae should be considered in the selection of sensitive antibiotics in clinical practice to improve the rationality of clinical medication (26).

References

  • 1.

    Halaby T, Kucukkose E, Janssen AB, Rogers MR, Doorduijn DJ, van der Zanden AG, et al. Genomic characterization of colistin heteroresistance in Klebsiella pneumoniae during a nosocomial outbreak. Antimicrob Agents Chemother. 2016;60(11):6837-43. [PubMed ID: 27600049]. [PubMed Central ID: PMC5075093]. https://doi.org/10.1128/AAC.01344-16.

  • 2.

    Chew KL, Lin RTP, Teo JWP. Klebsiella pneumoniae in Singapore: Hypervirulent infections and the carbapenemase threat. Front Cell Infect Microbiol. 2017;7:515. [PubMed ID: 29312894]. [PubMed Central ID: PMC5732907]. https://doi.org/10.3389/fcimb.2017.00515.

  • 3.

    Wang G, Zhao G, Chao X, Xie L, Wang H. The characteristic of virulence, biofilm and antibiotic resistance of Klebsiella pneumoniae. Int J Environ Res Public Health. 2020;17(17). [PubMed ID: 32872324]. [PubMed Central ID: PMC7503635]. https://doi.org/10.3390/ijerph17176278.

  • 4.

    Bonnin RA, Bernabeu S, Jaureguy F, Naas T, Dortet L. MCR-8 mediated colistin resistance in a carbapenem-resistant Klebsiella pneumoniae isolated from a repatriated patient from Morocco. Int J Antimicrob Agents. 2020;55(4):105920. [PubMed ID: 32062001]. https://doi.org/10.1016/j.ijantimicag.2020.105920.

  • 5.

    Yang X, Sun Q, Li J, Jiang Y, Li Y, Lin J, et al. Molecular epidemiology of carbapenem-resistant hypervirulent Klebsiella pneumoniae in China. Emerg Microbes Infect. 2022;11(1):841-9. [PubMed ID: 35236251]. [PubMed Central ID: PMC8942559]. https://doi.org/10.1080/22221751.2022.2049458.

  • 6.

    Zhang R, Liu L, Zhou H, Chan EW, Li J, Fang Y, et al. Nationwide surveillance of clinical carbapenem-resistant enterobacteriaceae (CRE) strains in China. EBioMedicine. 2017;19:98-106. [PubMed ID: 28479289]. [PubMed Central ID: PMC5440625]. https://doi.org/10.1016/j.ebiom.2017.04.032.

  • 7.

    Han R, Shi Q, Wu S, Yin D, Peng M, Dong D, et al. Dissemination of carbapenemases (KPC, NDM, OXA-48, IMP, and VIM) among carbapenem-resistant enterobacteriaceae isolated from adult and children patients in China. Front Cell Infect Microbiol. 2020;10:314. [PubMed ID: 32719751]. [PubMed Central ID: PMC7347961]. https://doi.org/10.3389/fcimb.2020.00314.

  • 8.

    Chen Z, Qin Z, Jiang L, Tang Z. [Research progress on the mechanism of two-component systems in regulating carbapenem resistance of Klebsiella pneumoniae]. Zhonghua Wei Zhong Bing Ji Jiu Yi Xue. 2021;33(6):761-4. [PubMed ID: 34296703]. https://doi.org/10.3760/cma.j.cn121430-20210402-00540.

  • 9.

    Du F, Huang X, Wei D, Mei Y, Liu P, Liu Y, et al. Distribution of CRISPR-Cas system in Klebsiella pneumoniae infected by blood flow and its relationship with virulence genes and drug resistance. Chin J Antibiotics. 2019:586-590. Chinese.

  • 10.

    Liu J, Du SX, Zhang JN, Liu SH, Zhou YY, Wang XR. Spreading of extended-spectrum beta-lactamase-producing Escherichia coli ST131 and Klebsiella pneumoniae ST11 in patients with pneumonia: a molecular epidemiological study. Chin Med J (Engl). 2019;132(16):1894-902. [PubMed ID: 31408445]. [PubMed Central ID: PMC6708689]. https://doi.org/10.1097/CM9.0000000000000368.

  • 11.

    Walker KA, Miller VL. The intersection of capsule gene expression, hypermucoviscosity and hypervirulence in Klebsiella pneumoniae. Curr Opin Microbiol. 2020;54:95-102. [PubMed ID: 32062153]. [PubMed Central ID: PMC8121214]. https://doi.org/10.1016/j.mib.2020.01.006.

  • 12.

    Fu Y, Dong A, Zhou H, Wang Y, Huang J, Zhang L, et al. Molecular epidemiological characteristics and drug resistance analysis of membrane porin ompK36 variant and virulence in Klebsiella pneumoniae. Chin J Nosocomiol. 2020:2603-2609. Chinese.

  • 13.

    Liao W, Zhu L, Du F, Long D, Zhang W, Liu Y. [Distribution of CRISPR-Cas systems and its potential relationship with virulence and antibiotic resistance in clinical hypervirulent Klebsiella pneumoniae strains]. Chin J Infect Chemother. 2020;20(2):163-8. https://doi.org/10.16718/j.1009-7708.2020.02.010.

  • 14.

    Chen M, Li Y, Zhang J, Chen S, Zhang Y. Molecular genetic characteristics of highly virulence Klebsiella pneumoniae in patients with respiratory system infection. Chin J Nosocomiol. 2020:6-9. Chinese.

  • 15.

    Humphries R, Bobenchik AM, Hindler JA, Schuetz AN. Overview of changes to the clinical and laboratory standards institute performance standards for antimicrobial susceptibility testing, M100, 31st edition. J Clin Microbiol. 2021;59(12). e0021321. [PubMed ID: 34550809]. [PubMed Central ID: PMC8601225]. https://doi.org/10.1128/JCM.00213-21.

  • 16.

    Weinstein MP, Lewis J2. The clinical and laboratory standards institute subcommittee on antimicrobial susceptibility testing: Background, organization, functions, and processes. J Clin Microbiol. 2020;58(3). [PubMed ID: 31915289]. [PubMed Central ID: PMC7041576]. https://doi.org/10.1128/JCM.01864-19.

  • 17.

    Kim DK, Kim HS, Pinto N, Jeon J, D'Souza R, Kim MS, et al. Xpert CARBA-R assay for the detection of carbapenemase-producing organisms in intensive care unit patients of a korean tertiary care hospital. Ann Lab Med. 2016;36(2):162-5. [PubMed ID: 26709264]. [PubMed Central ID: PMC4713850]. https://doi.org/10.3343/alm.2016.36.2.162.

  • 18.

    Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884-90. [PubMed ID: 30423086]. [PubMed Central ID: PMC6129281]. https://doi.org/10.1093/bioinformatics/bty560.

  • 19.

    Prjibelski A, Antipov D, Meleshko D, Lapidus A, Korobeynikov A. Using SPAdes De Novo Assembler. Curr Protoc Bioinformatics. 2020;70(1). e102. [PubMed ID: 32559359]. https://doi.org/10.1002/cpbi.102.

  • 20.

    Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;3:124. [PubMed ID: 30345391]. [PubMed Central ID: PMC6192448]. https://doi.org/10.12688/wellcomeopenres.14826.1.

  • 21.

    Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45(D1):D566-73. [PubMed ID: 27789705]. [PubMed Central ID: PMC5210516]. https://doi.org/10.1093/nar/gkw1004.

  • 22.

    Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis--10 years on. Nucleic Acids Res. 2016;44(D1):D694-7. [PubMed ID: 26578559]. [PubMed Central ID: PMC4702877]. https://doi.org/10.1093/nar/gkv1239.

  • 23.

    Pan YJ, Lin TL, Chen YH, Hsu CR, Hsieh PF, Wu MC, et al. Capsular types of Klebsiella pneumoniae revisited by wzc sequencing. PLoS One. 2013;8(12). e80670. [PubMed ID: 24349011]. [PubMed Central ID: PMC3857182]. https://doi.org/10.1371/journal.pone.0080670.

  • 24.

    Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530-4. [PubMed ID: 32011700]. [PubMed Central ID: PMC7182206]. https://doi.org/10.1093/molbev/msaa015.

  • 25.

    Shi Q, Yin D, Han R, Guo Y, Zheng Y, Wu S, et al. Emergence and recovery of ceftazidime-avibactam resistance in blaKPC-33-Harboring Klebsiella pneumoniae sequence type 11 isolates in China. Clin Infect Dis. 2020;71(Suppl 4):S436-9. [PubMed ID: 33367577]. https://doi.org/10.1093/cid/ciaa1521.

  • 26.

    Beyrouthy R, Dalmasso G, Birer A, Robin F, Bonnet R. Carbapenem resistance conferred by OXA-48 in K2-ST86 hypervirulent Klebsiella pneumoniae, France. Emerg Infect Dis. 2020;26(7):1529-33. [PubMed ID: 32568057]. [PubMed Central ID: PMC7323517]. https://doi.org/10.3201/eid2607.191490.

  • 27.

    Shi Q, Han R, Guo Y, Yang Y, Wu S, Ding L, et al. Multiple novel ceftazidime-avibactam-resistant variants of bla(KPC-2)-positive Klebsiella pneumoniae in two patients. Microbiol Spectr. 2022;10(3). e0171421. [PubMed ID: 35588280]. [PubMed Central ID: PMC9241591]. https://doi.org/10.1128/spectrum.01714-21.

  • 28.

    Giddins MJ, Macesic N, Annavajhala MK, Stump S, Khan S, McConville TH, et al. Successive emergence of ceftazidime-avibactam resistance through distinct genomic adaptations in bla(KPC-2)-harboring Klebsiella pneumoniae sequence type 307 isolates. Antimicrob Agents Chemother. 2018;62(3). [PubMed ID: 29263067]. [PubMed Central ID: PMC5826117]. https://doi.org/10.1128/AAC.02101-17.

  • 29.

    Meng X, Yang J, Duan J, Liu S, Huang X, Wen X, et al. Assessing molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae (CR-KP) with MLST and MALDI-TOF in central China. Sci Rep. 2019;9(1):2271. [PubMed ID: 30783127]. [PubMed Central ID: PMC6381170]. https://doi.org/10.1038/s41598-018-38295-8.

  • 30.

    Tan S, Zhang X, Lvu W, Li L, Zhang L, Bi Q. Differences and changes of antimicrobial resistance of Klebsiella pneumoniae isolated from patients with healthcare-associated infection and community-associated infection. Chin J Infect Control. 2020;19(11). Chinese.

  • 31.

    Chen Y, Marimuthu K, Teo J, Venkatachalam I, Cherng BPZ, De Wang L, et al. Acquisition of plasmid with carbapenem-resistance gene bla(KPC2) in hypervirulent Klebsiella pneumoniae, Singapore. Emerg Infect Dis. 2020;26(3):549-59. [PubMed ID: 32091354]. [PubMed Central ID: PMC7045839]. https://doi.org/10.3201/eid2603.191230.

  • 32.

    Paczosa MK, Mecsas J. Klebsiella pneumoniae: Going on the offense with a strong defense. Microbiol Mol Biol Rev. 2016;80(3):629-61. [PubMed ID: 27307579]. [PubMed Central ID: PMC4981674]. https://doi.org/10.1128/MMBR.00078-15.

  • 33.

    Jagnow J, Clegg S. Klebsiella pneumoniae MrkD-mediated biofilm formation on extracellular matrix- and collagen-coated surfaces. Microbiology (Reading). 2003;149(Pt 9):2397-405. [PubMed ID: 12949165]. https://doi.org/10.1099/mic.0.26434-0.

  • 34.

    El Fertas-Aissani R, Messai Y, Alouache S, Bakour R. Virulence profiles and antibiotic susceptibility patterns of Klebsiella pneumoniae strains isolated from different clinical specimens. Pathol Biol (Paris). 2013;61(5):209-16. [PubMed ID: 23218835]. https://doi.org/10.1016/j.patbio.2012.10.004.