Molecular and clinical characterization of hypervirulent Klebsiella pneumoniae isolates from individuals with urinary tract infections

Despite being a significant public health concern, hypervirulent Klebsiella pneumoniae (hvKP) has rarely been investigated in urinary tract infections (UTIs). To investigate the molecular and clinical characterization of hvKP in UTIs, we collected K. pneumoniae strains and clinical data from patients with UTIs. HvKP was confirmed by virulence-related genes and the Galleria mellonella model and sequenced by next-generation sequencing. Our data showed that 30/121 isolates were hvKP [17 carbapenem-resistant hypervirulent K. pneumoniae (CR-hvKP), 12 hvKP, and 1 extended-spectrum β-lactamase-producing hvKP]; these had higher resistance to most antimicrobials and were more likely to cause complicated UTIs (cUTIs). Notably, the mucoid phenotype-regulating genes _prmpA and _prmpA2 were truncated in 3 and 19 hvKP, respectively. Eight serotypes were detected and divided into three groups: K64 (n = 17), K1/K2 (n = 6), and others (n = 7). Furthermore, 16/17 K64 hvKP isolates were CR-hvKP but with a lower mortality rate of G. mellonella as the truncated _prmpA/_prmpA2 incurred high fitness cost to the isolates. In addition, all K64 isolates belonged to ST11 with the same cluster, and in two of these strains (KP88 and KP92) bla_KPC-2 gene was successfully transferred to EC600. Genetic environment analysis showed that IS26–tnpR–ISKpn27–bla_KPC−2–ISKpn6 may be the core structure in the horizontal transfer of bla_KPC-2. The highest mortality rate among the infected G. mellonella was observed in the K1/K2 group. In conclusion, hvKP had a higher resistance rate and was more likely to lead to cUTIs. Convergence of hypervirulence and carbapenem resistance in a transmissible ST11 clone of K64 K. pneumoniae was mediated by a plasmid in UTIs. Therefore, surveillance of hvKP in UTIs should be strengthened.

Introduction

Globally, more than 150 million people are affected with urinary tract infections (UTIs) yearly, making it as the most frequent infections (Öztürk and Murt, 2020). To confirm a diagnosis of UTI, a combination of positive urine culture, pyuria, and the relevant signs and symptoms is needed. A characteristic finding of UTIs is the presence of microbial pathogens in the urinary tract; it is divided into symptomatic UTI (with associated signs and symptoms) and asymptomatic bacteriuria (without associated signs and symptoms) (Mody and Juthani-Mehta, 2014). Typically, symptomatic UTI is divided into complicated and uncomplicated UTIs, and in complicated UTIs (cUTIs), patients often have functional or structural abnormalities or systemic diseases (e.g., renal insufficiency, immunodeficiency, or diabetes); furthermore, they are at a greater risk of developing systemic infection symptoms, such as fever and sepsis (Mody and Juthani-Mehta, 2014).

Previous studies revealed the following as some of the pathogens causing UTIs: Escherichia coli, Klebsiella pneumoniae, Enterococcus faecium, and Enterococcus faecalis (Nicolle, 2015). K. pneumoniae was the second most prominent pathogen. Hypervirulent K. pneumoniae (hvKP) is a more easily transmitted, more pathogenic, and a more fatal variant of K. pneumoniae (KP) (Lee et al., 2017; Russo and Marr, 2019). Notably, hvKP can cause infection at multiple locations: endophthalmitis, abdominal disease, thoracic disease, central nervous system disease, and genitourinary tract infections (Russo and Marr, 2019). In recent times, its multidrug-resistant and extensively drug-resistant isolates have also been detected in UTIs (Mataseje et al., 2019; Zafar et al., 2019; Zhu et al., 2020). The finding that hvKP infection is more likely to cause serious consequences and bacteremia than classical K. pneumoniae (cKP) as a causative agent of UTIs is particularly concerning (Kim et al., 2017).

Although hvKP has been reported in UTIs, systematic reports of hvKP in UTIs have been rare (Zafar et al., 2019; Taraghian et al., 2020). Therefore, herein, we obtained hvKP isolates from patients with UTIs and compared the clinical and molecular characteristics of these isolates with those of cKP to provide laboratory basis for the empirical treatment of UTIs.

Methods

Materials and methods

We collected urine samples from patients with symptomatic UTIs from October 2019 to October 2020 at Central South University’s Xiangya Hospital located at Changsha, China. Patients were considered as having symptomatic UTI if they exhibited typical signs and symptoms of a UTI, with bacteriuria [≥10⁵ colony forming units (CFUs)/mL] and pyuria (≥10 white blood cells/high-powered field) (Mody and Juthani-Mehta, 2014). Furthermore, we collected demographic data as well as information pertaining to underlying diseases and clinical manifestations. Freshly voided midstream urine (5–10 mL) was aseptically collected in wide-mouthed, sterile plastic bottles having a tight cap. Then, 0.001 mL of uncentrifuged, uniformly mixed samples were aseptically inoculated onto blood agar plates (Guangzhou Boret Biotechnology Co., Ltd., Guangzhou, China) using a calibrated inoculating loop. After incubating overnight for 24–48 h at 37°C, any significant growth was checked by counting colonies. Bacteriuria was indicated by the presence of colonies yielding a bacterial growth of 10⁵ CFUs per one milliliter of urine (Mody and Juthani-Mehta, 2014). Matrix-assisted laser desorption/ionization–time-of-flight mass spectrometry (MALDI–TOF MS; Bruker Daltonics GmbH, Bremen, Germany) was used to identify all isolates, and the quality control strain used was E. coli ATCC 25922 (National Center for Clinical Laboratories, Beijing, China).

HvKP was screened as carriage of at least two of the _prmpA, _prmpA2, iroB, peg344, and iucA, and confirmed with G. mellonella infection model (Liu et al., 2020).

Antimicrobial susceptibility testing

The minimum inhibitory concentrations (MICs) of the antimicrobial agents used were determined using the broth microdilution test according to the manufacturer instructions. The agents included cefepime (FEP), ceftazidime (CAZ), levofloxacin (LEV), ceftriaxone (CRO), piperacillin/tazobactam (TZP), meropenem (MEM), amikacin (AK), imipenem (IPM), ciprofloxacin (CIP), gentamicin (GEN), aztreonam (ATM), trimethoprim/sulfamethoxazole (SXT), nitrofurantoin (F), ceftazidime/avibactam (CZA), tigecycline (TGC), and polymyxin B (PB), and all of these agents were obtained from Hangzhou Kangtai Biotechnology Co., Ltd.. The guidelines of the Clinical and Laboratory Standards Institute’s (2020) were used to interpret susceptibility breakpoints (CLSI, 2020). The breakpoints indicated by the US Food and Drug Administration determined the MIC of TGC. E. coli ATCC 25922 was used as the quality control strain.

Phenotypic screening of extended-spectrum β-lactamase

The presence of ESBLs was evaluated using CLSI criteria for screening these bacteria using the double-disk diffusion method. An increase of >5 mm in the inhibition zone of the ceftazidime–clavulanic acid disk when compared with the ceftazidime alone disk was interpreted as phenotypic evidence of ESBL production.

Detection of hypermucoviscosity and virulence-associated genes

The hypermucoviscosity of K. pneumoniae was determined by string test as described previously (Li et al., 2020). Briefly, after subculturing K. pneumoniae isolates overnight at 37°C on 5% sheep blood agar, a bacterial loop was used to stretch a “string” of mucous from the bacterial colony. A mucoviscous string longer than 5 mm indicated hypermucoviscosity (Li et al., 2020).

For all K. pneumoniae strains, their genomic DNA was extracted using overnight-cultured strains using the boiling method (Li et al., 2016). For detecting capsular serotype genes (K1, K2, K5, K20, K54, K57, and K64), _prmpA, _prmpA2, iroB, peg344, and iucA, we performed by polymerase chain reaction (PCR) as described previously (Yu et al., 2018; Liu et al., 2020). This was followed by subjecting the positive PCR products to direct Sanger sequencing.

Galleria mellonella infection model

We used the G. mellonella infection model (Tianjin Huiyude Biotech Company, Tianjin, China) for evaluating the virulence of the collected isolates (Zhou et al., 2018). After adjusting the K. pneumoniae overnight cultures to 1 × 10⁸ CFUs/mL with the use of phosphate-buffered saline, the bacterial culture (10 μL) was injected into G. mellonella larvae, followed by incubation of the larvae in the dark for 5 days at 37°C; the survival was monitored during this period. K. pneumoniae ATCC 700603 and K. pneumoniae NTUH-K2044 were the low- and high-virulence control strains used, respectively. All experiments were performed in triplicate.

Biofilm formation assay

To test if KP isolates formed a biofilm, the previously described crystal violet staining method was implemented (Coffey and Anderson, 2014). The absorbance was observed at 570 nm, and for each assay performed in triplicate, data are represented in the form of mean ± standard deviation (SD). K. pneumoniae NTUH-K2044 and ATCC 700606 were the positive and negative controls used, respectively. The optical density (OD) cut-off (ODc) value was calculated according to the parameters defined by Ragupathi et al. (Devanga Ragupathi et al., 2020) as follows: ODc = average OD of the negative control + (3 × SD of the negative control). Using the OD results, the biofilm formation ability was categorized as follows: (1) medium biofilm producer (4 × ODc ≥ OD > 2 × ODc); (2) strong biofilm producer (OD > 4 × ODc); (3) non-biofilm producer (OD ≤ ODc); and (4) weak biofilm producer (2 × ODc ≥ OD > ODc).

Whole-genome sequencing (WGS)

Genomic DNA was isolated using MagAttract HMW DNA Kit (Qiagen, Hilden, Germany) and subjected to high-throughput next-generation sequencing (NGS) on a HiSeq 2000™ platform (Illumina Inc., San Diego, CA, USA) with 2 × 100-bp paired-end reads and to long-read high-throughput sequencing (LRS) on a MinION platform (Oxford Nanopore Technologies, Oxford, UK). The long reads generated by MinION were assembled using Canu v. 1.6 (Koren et al., 2017) and polished with the short reads generated by HiSeq using Pilon v1.22 (Walker et al., 2014) to obtain the whole genome and complete plasmid sequences. The chromosome and plasmid sequences were annotated using the prokaryotic gene prediction tool Prokka (Seemann, 2014). The plasmid type and the resistant genes were identified using the CGE server (https://cge.cbs.dtu.dk). The phylogenetic tree was constructed using the BacWGSTdb software (Institute of Translational Medicine, Zhejiang University, China) based on the previously reported single nucleotide polymorphism strategy (Ruan and Feng, 2016). For MLST, the Oxford scheme was used, and the sequence type (ST) allocation was carried out using the MLST database (https://bigsdb.pasteur.fr/klebsiella/klebsiella.html). Notably, the sequence data have been deposited in NCBI with the accession number PRJNA806839.

Conjugation experiment

For assessing the transferability of carbapenemase-producing genes and identify whether these genes were located on plasmids, conjugation experiment was carried out as previously reported (Zeng et al., 2021). Briefly, the recipient (E. coli EC600) and donor strains (carbapenem-resistant K. pneumoniae, CRKP) were both mixed in a 4:1 ratio in Luria–Bertani broth; then, the mixtures were placed on a membrane, which was followed by incubation at 35°C for 24 h. Transconjugant selection was done on meropenem (1 μg/mL)- and rifampicin (600 μg/mL)-supplemented Mueller–Hinton agar plates. Next, the colonies that grew on the selective medium were identified using MALDI-TOF MS. Transconjugants were identified as carbapenemase-producing strains that exhibited a higher MIC of resistance to carbapenems than EC600, and then, PCR was used to confirm the presence of resistance determinants.

Statistical analyses

In patients with cKP and hvKP strains, the following aspects were compared: clinical characteristics, biofilm formation, and antimicrobial susceptibility. For categorical variables, the Fisher’s exact test or χ² test was used, and for continuous variables, the Student’s t-test was used. For statistical analysis, SPSS 21.0 was used, and P < 0.05 indicated statistical significance.

Results

Clinical characteristics of cKP and hvKP

We isolated a total of 121 K. pneumoniae strains from urine samples of 121 patients with UTIs. The mean age of these patients was 57.53 years (range, 1–96 years). The demographic data of the patients revealed that most patients had underlying diseases, including bacteremia (97, 80.2%), hypertension (49, 40.5%), diabetes mellitus (37, 30.6%), and bladder catheter-related infection (35, 28.9%) (Table 1).

TABLE 1

Table 1 Clinical characteristics and susceptibility of hvKP and cKP isolates.

To evaluate the virulence of these strains, the virulence-associated genes were identified and confirmed using the G. mellonella infection model. Among the 121 K. pneumoniae strains, 30 (24.7%) strains simultaneously carried at least two virulence genes (_prmpA, _prmpA2, iroB, peg344, and iucA), and the mortality rate of G. mellonella infected with these 30 strains was 50%–100%, which was significantly higher than that of those infected with ATCC 700603, suggesting that these strains were hvKP (Figure 1).

FIGURE 1

Figure 1 Evaluation of the virulence of KP isolates. Using a G. mellonella infection model, we investigated the virulence of 1 × 10⁶ CFU of each KP strain (KP1, KP7, KP8, KP14, KP16, KP23, KP28, KP92, and KP110). NTUH-K2044 and ATCC 700603 were used as high- and low-virulence KP controls, respectively. The survival rates of KP7- and KP16-infected G. mellonella were similar with ATCC 700603, while those of KP1, KP8, KP14, KP23, KP28, KP92, and KP110 were similar with NTUH-K2044, indicating that they were hypervirulent.

The age, sex, underlying disease, cUTI incidence, and inflammatory indices (e.g., CRP and PCT) did not significantly differ between patients with hvKP and cKP infections. In contrast to hvKP, cKP isolates were mainly derived from inpatients (70.0% vs. 89.0%, P = 0.0283) and had superior biofilm-forming ability (70.0% vs. 87.9%, P = 0.0447). In addition, carbapenem resistance was significantly more prevalent in hvKP than in cKP (56.7% vs. 23.1%, P = 0.0010), and hvKP showed significantly more resistance than cKP to most antimicrobials, including TZP, FEP, ATM, IPM, MEM, and AK, whereas resistance to CAZ, CRO, CIP, SXT, and F did not significantly differ between the two. The detection rate of extended-spectrum β-lactamase (ESBL)-producing strains was higher in cKP than in hvKP (3.3% vs. 38.5%, P = 0.0003) (Table 1).

Antimicrobial susceptibility, virulence-associated genes, and resistant genes of hvKP

Among the 30 hvKP isolates, 17 were carbapenem-resistant Enterobacterales (CRE) and 1 was an ESBL-producing strain (Table 1). The resistance rates of hvKP to TZP, CAZ, CRO, FEP, ATM, IPM, MEM, AK, CIP, SXT, and F were 56.7%, 56.7%, 66.7%, 60.0%, 63.3%, 56.7%, 56.7%, 50.0%, 70.0%, 36.7%, and 76.7%, respectively. There were 3 TGC-resistant strains, and no isolate was resistant to polymyxins B and ceftazidime/avibactam. In addition, 15 in 17 CR-hvKP isolates carried the bla_KPC-2 gene, and the remaining harbored mutations in the Ompk37-encoding gene, suggesting that the bla_KPC-2 gene is the main reason for their resistance to carbapenems (Figure 2).

FIGURE 2

Figure 2 Molecular characterization of 30 hvKP isolates. NA, not determined. Homology analysis of 30 hvKP isolates was performed using BacWGSTdb software, and all K64 ST11 CR-hvKP strains belonged to the same cluster.

The results of the string test revealed that 13 (five hvKP and eight cKP strains) of 121 KP strains showed hypermucoviscosity. Furthermore, _prmpA and _prmpA2 were identified as important regulators of the formation of the capsular polysaccharide of K. pneumoniae, and 18, 20 were detected in hvKP, respectively. Notably, the _prmpA gene was truncated in three hvKP strains (KP1, KP9, and KP30), and the _prmpA2 gene was truncated in 19 hvKP strains (Figure 2).

Eight serotypes were detected, including K1, K2, K5, K20, K23, K57, K62, and K64, with the rates of 3.3% (1/30), 16.7% (5/30), 3.3% (1/30), 3.3% (1/30), 6.7% (2/30), 6.7% (2/30), 3.3% (1/30), and 56.7% (17/30), respectively. K64 was the predominant serotype in our study. In addition, these serotypes were divided into three groups (K64, K1/K2, and others), and the drug resistance and virulence of different groups of strains were compared. Herein, 16 of 17 K64 isolates were CRE and 15 carried the bla_KPC-2 gene, which is predominantly responsible for conferring carbapenem resistance to K. pneumoniae in China. Moreover, one of six K1/K2 isolates was CRE and did not carry the bla_KPC-2 gene, and seven isolates of other serotypes were not detected as CRE. K64 isolates showed 64.7% resistance to SXT, and the resistance rates of the remaining tested antimicrobials were above 82.4% and showed higher resistance rates to all tested antimicrobials than the other two groups. The distribution of the mucoid phenotype regulating genes _prmpA and _prmpA2 differed among the three groups: _prmpA2 was mainly distributed in the K64 group with all truncated gene sequences; _prmpA was mainly distributed in the K1/K2 group; and the other group had both distributions. The highest mortality rate in infected G. mellonella was observed in the K1/K2 group, followed by the other group and the K64 group (Table 2).

TABLE 2

Table 2 Comparison of drug resistance and virulence among K64, K1/K2 and other serotype hvKP.

Molecular characteristics

We identified the following 10 ST types in the 30 hvKP isolates by MLST: ST11 (n = 17), ST23 (n = 1), ST25 (n = 3), ST36 (n = 1), ST39 (n = 1), ST65 (n = 1), ST412 (n = 2), ST420 (n = 1), ST881 (n = 1), and ST1049 (n = 1); moreover, one novel ST type was identified as well. K64 hvKP isolates were all ST11. K1 isolates were ST23, and K2 isolates were ST25, ST65, and ST881; K23 and K62 isolates were ST39 and ST36, respectively. Homology analysis of 30 hvKP isolates was performed using BacWGSTdb software (Institute of Translational Medicine, Zhejiang University, China), and all K64 ST11 CR-hvKP strains belonged to the same cluster, suggesting clonal transmission. Other serotype isolates were clustered in another clusters (Figure 2).

Conjugation tests

The conjugation tests were performed on bla_KPC-2-positive hvKP strains, and the results showed that two K64 ST11 CR-hvKP strains (KP88 and KP92) successfully transferred the bla_KPC-2 gene into the recipient strain (EC600). In addition, transconjugants (KP88* and KP92*) were concurrently resistant to CAZ, TZP, FEP, ATM, MEM, IPM, and LEV, and KP92* also showed resistance to AK, suggesting that the genes responsible for resistance to these agents were also transferred to the recipient strain. Simultaneously, we also tested virulence genes by PCR methods in KP88* and KP92*; however, these were all negative (Table 3).

TABLE 3

Table 3 Susceptibility of KP88, KP88*, KP92, KP92*, and EC600.

The sequence of the KP92 strain comprises a chromosome and three plasmids (plasmids 1, 2, and 3). Plasmid 1 (size, 144551 bp) is an IncFII plasmid harboring bla_KPC-2, bla_SHV-12, bla_CTX-M-65, and bla_TEM-1B acquired antimicrobial resistance genes. Plasmid 2 (size, 216439 bp) is an IncFIB plasmid that harbors rmpA2 (258 bp, truncated rmpA2), iucABCD, and other virulence genes. Finally, Plasmid 3 (size, 47106 bp) is an incN plasmid harboring qnrs1 resistance gene (Table 4). The genetic environment of the bla_KPC-2 gene is IS26–tnpR–ISKpn27–bla_KPC−2–ISKpn6–korC–ORF–klcA–ISKpn19, which is highly similar to the CP065952.1 sequence identified in patients of Taipei Veterans General Hospital.

TABLE 4

Table 4 Genome characteristics of isolate KP92.

Discussion

The emergence of hypervirulent K. pneumoniae makes clinical anti-infective therapy extremely challenging. HvKP was originally reported in patients with liver abscesses and is currently found to cause several diseases, such as pneumonia and endophthalmitis; however, it has been less frequently reported in UTIs (Russo and Marr, 2019; Zafar et al., 2019; Taraghian et al., 2020). In this study, we screened hvKP isolated from the urine samples of patients diagnosed as having UTIs and investigated the molecular and clinical characteristics of these strains to provide empirical treatment for UTIs.

The detection rate of hvKP in UTIs in our study was 24.7%, indicating a higher incidence than that reported in the study of Taraghian et al. (Taraghian et al., 2020) but lower than that reported in the study of Liu et al. (Liu et al., 2020). In addition, hvKP was predominantly derived from hospitalized patients in this study and not from community-acquired infections as in previous reports (Shankar et al., 2018; Gao et al., 2020). Diabetes has been suspected to be a risk factor for hvKP infection, but patients with hvKP and cKP did not significantly differ in this regard in this study (Zhang et al., 2016). However, this study found that compared to other serotypes, K1/K2 serotypes were associated with a higher incidence of diabetes, suggesting diabetes as a risk factor for K1/K2 hvKP infection. Consistent with our previous results, in this study as well, hvKP was found to have a weaker biofilm-forming ability than cKP (Li et al., 2022). This finding may be attributed to hvKP being more likely to escape the phagocytosis by urinary macrophages to cause bacteremia and serious consequences; the specific mechanism behind this needs to be further investigated (Kim et al., 2017). Furthermore, hvKP was more likely to cause cUTIs; however, systemic infection indicators, such as CRP and PCT, did not significantly differ between cKP and hvKP infections.

The hypervirulence combined with drug-resistant strains is associated with increasing emergence of such strains and poses a major challenge to clinical anti-infective therapy (Lee et al., 2017; Mataseje et al., 2019; Zhu et al., 2020; Zhang et al., 2021). In this study, the detection rate of CRE in hvKP was higher than it was in cKP; this could be attributed to bla_KPC-2, which was identified to be successfully transferred into EC600 by conjugation test, indicating that this hypervirulent strain could acquire resistance genes by horizontal gene transfer (HGT) (Lan et al., 2021). We also confirmed the location of the bla_KPC-2 gene on the IncFII plasmid, which was an important plasmid type for HGT of the bla_KPC-2 gene (Andrade et al., 2011; Shi et al., 2018). The genetic environment of the bla_KPC-2 gene analysis revealed that the gene was located on IS26–tnpR–ISKpn27–bla_KPC−2–ISKpn6–korC–ORF–klcA–ISKpn19, which is highly similar to the CP065952.1 sequence identified in patients of Taipei Veterans General Hospital, suggesting that the IS26–tnpR–ISKpn27–bla_KPC−2–ISKpn6 may be the core structure in the horizontal transfer of bla_KPC-2. In addition, all CR-hvKP belonged to K64-ST11 and had the same cluster with the possibility of clonal transmission. ST11-KPC-2 K pneumoniae is the most prevalent clone in China, and it was also found in CR-hvKP (Zhang et al., 2020). Therefore, we should strengthen its monitoring and take effective control measures to prevent its further dissemination.

Initial studies suggested that hypermucoviscosity (string test > 5 mm) was an important feature of hvKP, and strains with hypermucoviscosity were considered as hvKP. However, subsequent researches have reported that the correlation between hypermucoviscosity and hypervirulence varies greatly (Fang et al., 2004; Yu et al., 2008; Lee et al., 2010; Lin et al., 2011) and that a certain proportion of cKP strains also has hypermucoviscosity (Fang et al., 2004; Yu et al., 2008; Lee et al., 2010; Lin et al., 2011), indicating that hypermucoviscosity cannot be used as a reliable method to identify hvKP. In our study, only five of 30 hvKP strains were positive for hypermucoviscocity, and eight of 91 cKP strains were also hypermucoviscous, further confirming that hypermucoviscocity was not specific to hvKP (Fang et al., 2004; Lee et al., 2010). Capsular polysaccharide is an important reason behind the pathogenicity of K. pneumoniae, and the pathogenic ability of different capsular polysaccharides may differ with specific serotypes. The common serotypes of hvKP are K1, K2, K5, K20, K54, and K57, mainly K1 and K2. Recently, K64 identified by wzc or wzi sequencing was considered as the main serotype of CR-hvKP, which has been seldom reported in the past due to the limitations of the traditional serotype identification methods (Pan et al., 2015). Zhang et al. (Zhang et al., 2020) reported that K64 has become a common type of CR-hvKP in China (80%, 44/55). In our study, eight serotypes, namely K1, K2, K5, K20, K23, K57, K62, and K64 (mainly K64), were detected and divided into three groups (K64, K1/K2, and others), and K64 hvKP isolates had higher resistance rate to commonly used antimicrobials but were associated with a lower mortality rate in infected G. mellonella, which may be attributed to the truncation of_prmpA/_prmpA2 incurring high fitness cost in CR-hvKP (Lin et al., 2020; Shankar et al., 2021). Choi et al. (Choi and Ko, 2015) suggested that impaired virulence and reduced capsular polysaccharide production increases fitness cost in the acquisition of colistin-resistant K. pneumoniae. However, in the present study, the K64 serotype CR-hvKP strains were still associated with high mortality of G. mellonella, which is in line with the results of Shankar et al. (Shankar et al., 2021), and this may be related to the increasing simultaneous occurrence of drug-resistance and hypervirulence in K. pneumoniae strains observed globally. Thus, hvKP should be paid more attention to when surveillance is performed in the areas of high hvKP endemicity. The K1/K2 group had the highest mortality rate in infected G. mellonella, followed by the other group and the K64 group.

In summary, the clinical and molecular characteristics of hvKP isolated from UTIs were investigated in this study. Notably, compared with cKP, hvKP had higher resistance rate and was more likely to cause cUTIs. The incidence of CR-hvKP isolates was increasing among hvKP, particularly K64-ST11 CR-hvKP, which was associated with higher resistance but lower mortality in infected G. mellonella because of the truncation of_prmpA/_prmpA2 incurring high fitness cost in CR-hvKP. In addition, the clonal distribution of CR-hvKP isolates was mediated by a plasmid in UTIs in our study. Therefore, surveillance of hvKP in UTIs should be strengthened.

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: NCBI with the accession number PRJNA806839.

Author contributions

Study design: JL and MZ. Study conduct: FX, CM, ZL, and MT. Data collection: MT, FX, and CM. Data analysis: JL and MZ. Data interpretation: YH and HW. Drafting manuscript: JL. Revising manuscript content: MZ. Approving the final version of the manuscript: JL and MZ. All authors contributed to the article and approved the submitted version.

Funding

The National Natural Science Foundation of China (No. 81702068), Natural Science Foundation of Hunan Province (No. 2020JJ4886, No. 2021JJ40980), the Science Foundation of Hunan Health Commission in Hunan province (No. 202111000066), and Fundamental Research Funds for the Central Universities of Central South University (No. 2021zzts1043) supported this study.

Acknowledgments

We would like to express our gratitude to the staff of the Department of Clinical Laboratory at Central South University’s Xiangya Hospital for their assistance in collecting and identifying the K. pneumoniae isolates.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.

Abbreviations

hvKP, hypervirulent K. pneumoniae; UTI, urinary tract infection; CR-hvKP, carbapenem-resistant hypervirulent K. pneumoniae; cUTI, complicated urinary tract infection; KP, K. pneumoniae; cKP, classical K. pneumoniae; CFU, colony forming unit; MALDI–TOF MS, Matrix-assisted laser desorption/ionization–time-of-flight mass spectrometry; MIC, minimum inhibitory concentration; FEP, cefepime; CAZ, ceftazidime; LEV, levofloxacin; CRO, ceftriaxone; TZP, piperacillin/tazobactam; MEM, meropenem; AK, amikacin; IPM, imipenem; CIP, ciprofloxacin; GEN, gentamicin; ATM, aztreonam; SXT, trimethoprim/sulfamethoxazole; F, nitrofurantoin; CZA, ceftazidime/avibactam; TGC, tigecycline; PB, polymyxin B; PCR, polymerase chain reaction; SD, standard deviation; OD, optical density; NGS, next-generation sequencing; WGS, whole-genome sequencing; ST, sequence type; CRKP, carbapenem-resistant K. pneumoniae; ESBL, extended-spectrum β-lactamase; CRE, carbapenem-resistant Enterobacterales; HGT, horizontal gene transfer; PCT, procalcitonin.

References

Andrade, L. N., Curiao, T., Ferreira, J. C., Longo, J. M., Clímaco, E. C., Martinez, R., et al. (2011). Dissemination of bla_KPC-2 by the spread of Klebsiella pneumoniae clonal complex 258 clones (ST258, ST11, ST437) and plasmids (IncFII, IncN, IncL/M) among enterobacteriaceae species in Brazil. Antimicrob. Agents Chemother. 55, 3579–3583. doi: 10.1128/AAC.01783-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, M. J., Ko, K. S. (2015). Loss of hypermucoviscosity and increased fitness cost in colistin-resistant Klebsiella pneumoniae sequence type 23 strains. Antimicrob. Agents Chemother. 59, 6763–6773. doi: 10.1128/AAC.00952-15

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI. (2020) The performance stands for antimicrobial susceptibility testing, 30th ed CLSI supplement M100. Wayne, PA: Clinical and laboratory standards institute.

Google Scholar

Coffey, B. M., Anderson, G. G. (2014). Biofilm formation in the 96-well microtiter plate. Methods Mol. Biol. 1149, 631–641. doi: 10.1007/978-1-4939-0473-0_48

PubMed Abstract | CrossRef Full Text | Google Scholar

Devanga Ragupathi, N. K., Muthuirulandi Sethuvel, D. P., Triplicane Dwarakanathan, H., Murugan, D., Umashankar, Y., Monk, P. N., et al. (2020). The influence of biofilms on carbapenem susceptibility and patient outcome in device associated K. pneumoniae infections: Insights into phenotype vs genome-wide analysis and correlation. Front. Microbiol. 11, 591679. doi: 10.3389/fmicb.2020.591679

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, C. T., Chuang, Y. P., Shun, C. T., Chang, S. C., Wang, J. T. (2004). A novel virulence gene in Klebsiella pneumoniae strains causing primary liver abscess and septic metastatic complications. J. Exp. Med. 199, 697–705. doi: 10.1084/jem.20030857

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Q., Shen, Z., Qin, J., Liu, Y., Li, M. (2020). Antimicrobial resistance and pathogenicity determination of community-acquired hypervirulent Klebsiella pneumoniae. Microb. Drug Resist. 26, 1195–1200. doi: 10.1089/mdr.2019.0439

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. J., Kim, S. I., Kim, Y. R., Wie, S. H., Lee, H. K., Kim, S. Y., et al. (2017). Virulence factors and clinical patterns of hypermucoviscous Klebsiella pneumoniae isolated from urine. Infect. Dis. (Lond) 49, 178–184. doi: 10.1080/23744235.2016.1244611

PubMed Abstract | CrossRef Full Text | Google Scholar

Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., Phillippy, A. M. (2017). Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Lan, P., Jiang, Y., Zhou, J., Yu, Y. (2021). A global perspective on the convergence of hypervirulence and carbapenem resistance in Klebsiella pneumoniae. J. Glob Antimicrob. Resist. 25, 26–34. doi: 10.1016/j.jgar.2021.02.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C. R., Lee, J. H., Park, K. S., Jeon, J. H., Kim, Y. B., Cha, C. J., et al. (2017). Antimicrobial resistance of hypervirulent Klebsiella pneumoniae: Epidemiology, hypervirulence-associated determinants, and resistance mechanisms. Front. Cell Infect. Microbiol. 7, 483. doi: 10.3389/fcimb.2017.00483

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C. H., Liu, J. W., Su, L. H., Chien, C. C., Li, C. C., Yang, K. D. (2010). Hypermucoviscosity associated with Klebsiella pneumoniae-mediated invasive syndrome: a prospective cross-sectional study in Taiwan. Int. J. Infect. Dis. 14, e688–e692. doi: 10.1016/j.ijid.2010.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Li, Y., Tang, M., Xia, F., Min, C., Hu, Y., et al. (2022). Distribution, characterization, and antibiotic resistance of hypervirulent Klebsiella pneumoniae isolates in a Chinese population with asymptomatic bacteriuria. BMC Microbiol. 22, 29. doi: 10.1186/s12866-021-02413-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y. C., Lu, M. C., Tang, H. L., Liu, H. C., Chen, C. H., Liu, K. S., et al. (2011). Assessment of hypermucoviscosity as a virulence factor for experimental Klebsiella pneumoniae infections: comparative virulence analysis with hypermucoviscosity-negative strain. BMC Microbiol. 11, 50. doi: 10.1186/1471-2180-11-50

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z. W., Zheng, J. X., Bai, B., Xu, G. J., Lin, F. J., Chen, Z., et al. (2020). Characteristics of hypervirulent Klebsiella pneumoniae: Does low expression of rmpA contribute to the absence of hypervirulence? Front. Microbiol. 11, 436. doi: 10.3389/fmicb.2020.00436

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Shi, J., Zhao, Y., Xie, Y., Tang, Y., Jiang, X., et al. (2020). Identification of hypervirulent Klebsiella pneumoniae isolates using the string test in combination with galleria mellonella infectivity. Eur. J. Clin. Microbiol. Infect. Dis. 39, 1673–1679. doi: 10.1007/s10096-020-03890-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Du, P., Xiao, N., Ji, F., Russo, T. A., Guo, J. (2020). Hypervirulent Klebsiella pneumoniae is emerging as an increasingly prevalent K. pneumoniae pathotype responsible for nosocomial and healthcare-associated infections in Beijing, China. Virulence 11, 1215–1224. doi: 10.1080/21505594.2020.1809322

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zou, M., Dou, Q., Hu, Y., Wang, H., Yan, Q., et al. (2016). Characterization of clinical extensively drug-resistant Pseudomonas aeruginosa in the hunan province of China. Ann. Clin. Microbiol. Antimicrob. 15, 35. doi: 10.1186/s12941-016-0148-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Mataseje, L. F., Boyd, D. A., Mulvey, M. R., Longtin, Y. (2019). Two hypervirulent Klebsiella pneumoniae isolates producing a bla_KPC-2 carbapenemase from a Canadian patient. Antimicrob. Agents Chemother. 63, e00517–19. doi: 10.1128/AAC.00517-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Mody, L., Juthani-Mehta, M. (2014). Urinary tract infections in older women: A clinical review. Jama. 311, 844–854. doi: 10.1001/jama.2014.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicolle, L. E. (2015). Asymptomatic bacteriuria and bacterial interference. Microbiol. Spectr. 3. doi: 10.1128/microbiolspec.UTI-0001-2012

CrossRef Full Text | Google Scholar

Öztürk, R., Murt, A. (2020). Epidemiology of urological infections: A global burden. World J. Urol. 38, 2669–2679. doi: 10.1007/s00345-019-03071-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Y. J., Lin, T. L., Lin, Y. T., Su, P. A., Chen, C. T., Hsieh, P. F., et al. (2015). Identification of capsular types in carbapenem-resistant Klebsiella pneumoniae strains by wzc sequencing and implications for capsule depolymerase treatment. Antimicrob. Agents Chemother. 59, 1038–1047. doi: 10.1128/AAC.03560-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, Z., Feng, Y. (2016). BacWGSTdb, a database for genotyping and source tracking bacterial pathogens. Nucleic Acids Res. 44, D682–D687. doi: 10.1093/nar/gkv1004

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, T. A., Marr, C. M. (2019). Hypervirulent Klebsiella pneumoniae. Clin. Microbiol. Rev. 32. doi: 10.1128/CMR.00001-19

CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics. 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Shankar, C., Basu, S., Lal, B., Shanmugam, S., Vasudevan, K., Mathur, P., et al. (2021). Aerobactin seems to be a promising marker compared with unstable RmpA2 for the identification of hypervirulent carbapenem-resistant Klebsiella pneumoniae: In silico and In vitro evidence. Front. Cell Infect. Microbiol. 11, 709681. doi: 10.3389/fcimb.2021.709681

PubMed Abstract | CrossRef Full Text | Google Scholar

Shankar, C., Veeraraghavan, B., Nabarro, L. E. B., Ravi, R., Ragupathi, N. K. D., Rupali, P. (2018). Whole genome analysis of hypervirulent Klebsiella pneumoniae isolates from community and hospital acquired bloodstream infection. BMC Microbiol. 18, 6. doi: 10.1186/s12866-017-1148-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, L., Feng, J., Zhan, Z., Zhao, Y., Zhou, H., Mao, H., et al. (2018). Comparative analysis of bla_(KPC-2)- and rmtB-carrying IncFII-family pKPC-LK30/pHN7A8 hybrid plasmids from Klebsiella pneumoniae CG258 strains disseminated among multiple Chinese hospitals. Infect. Drug Resist. 11, 1783–1793. doi: 10.2147/IDR.S171953

PubMed Abstract | CrossRef Full Text | Google Scholar

Taraghian, A., Nasr Esfahani, B., Moghim, S., Fazeli, H. (2020). Characterization of hypervirulent extended-spectrum β-Lactamase-Producing Klebsiella pneumoniae among urinary tract infections: The first report from Iran. Infect. Drug Resist. 13, 3103–3111. doi: 10.2147/IDR.S264440

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PloS One. 9, e112963. doi: 10.1371/journal.pone.0112963

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, W. L., Ko, W. C., Cheng, K. C., Lee, C. C., Lai, C. C., Chuang, Y. C. (2008). Comparison of prevalence of virulence factors for Klebsiella pneumoniae liver abscesses between isolates with capsular K1/K2 and non-K1/K2 serotypes. Diagn. Microbiol. Infect. Dis. 62, 1–6. doi: 10.1016/j.diagmicrobio.2008.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, F., Lv, J., Niu, S., Du, H., Tang, Y. W., Pitout, J. D. D., et al. (2018). Multiplex PCR analysis for rapid detection of Klebsiella pneumoniae carbapenem-resistant (Sequence type 258 [ST258] and ST11) and hypervirulent (ST23, ST65, ST86, and ST375) strains. J. Clin. Microbiol. 56. doi: 10.1128/JCM.00731-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Zafar, S., Hanif, S., Akhtar, H., Faryal, R. (2019). Emergence of hypervirulent K. pneumoniae causing complicated UTI in kidney stone patients. Microb. Pathog. 135, 103647.

PubMed Abstract | Google Scholar

Zeng, L., Yang, C., Zhang, J., Hu, K., Zou, J., Li, J., et al. (2021). An outbreak of carbapenem-resistant Klebsiella pneumoniae in an intensive care unit of a major teaching hospital in chongqing, China. Front. Cell Infect. Microbiol. 11, 656070. doi: 10.3389/fcimb.2021.656070

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Jin, L., Ouyang, P., Wang, Q., Wang, R., Wang, J., et al. (2020). Evolution of hypervirulence in carbapenem-resistant Klebsiella pneumoniae in China: a multicentre, molecular epidemiological analysis. J. Antimicrob. Chemother. 75, 327–336. doi: 10.1093/jac/dkz446

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Wang, X., Wang, S., Sun, S., Li, H., Chen, H., et al. (2021). Emergence of colistin resistance in carbapenem-resistant hypervirulent Klebsiella pneumoniae under the pressure of tigecycline. Front. Microbiol. 12, 756580. doi: 10.3389/fmicb.2021.756580

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zhao, C., Wang, Q., Wang, X., Chen, H., Li, H., et al. (2016). High prevalence of hypervirulent Klebsiella pneumoniae infection in China: Geographic distribution, clinical characteristics, and antimicrobial resistance. Antimicrob. Agents Chemother. 60, 6115–6120. doi: 10.1128/AAC.01127-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, K., Tang, X., Wang, L., Guo, Z., Xiao, S., Wang, Q., et al. (2018). An emerging clone (ST457) of Acinetobacter baumannii clonal complex 92 with enhanced virulence and increasing endemicity in south China. Clin. Infect. Dis. 67, S179–s188. doi: 10.1093/cid/ciy691

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Sun, C., Chen, H., Zhan, X., Liu, C., Liu, X., et al. (2020). Co-Occurrence of three different plasmids in an extensively drug-resistant hypervirulent Klebsiella pneumoniae isolate causing urinary tract infection. J. Glob Antimicrob. Resist. 23, 203–210. doi: 10.1016/j.jgar.2020.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar