Genomic Analysis of KPC-2-Producing Klebsiella pneumoniae ST11 Isolates at the Respiratory Department of a Tertiary Care Hospital in Beijing, China

Background: Carbapenem-resistant Klebsiella pneumoniae (CRKP) is an important pathogen causing hospital-associated outbreaks worldwide. The spread of K. pneumoniae carbapenemase-2 (KPC-2)-producing CRKP is primarily associated with sequence type (ST) 11.

Methods: A total of 152 KPC-2-producing K. pneumoniae ST11 isolates were collected from the respiratory department of a tertiary care hospital in Beijing, China between 2009 and 2018. The genome sequencing of these isolates was performed on the HiSeq X Ten sequencer. Multilocus sequence typing (MLST), capsular type, plasmid replicon types and resistance genes were identified. Fifteen isolates were selected for the subsequent single-molecule real-time (SMRT) sequencing on the PacBio RS II. Alignment of the complete sequences of the plasmids carrying bla_KPC–2 and/or virulence genes was performed by using BRIG and Easyfig.

Results: From 2012 to 2018, the detection rate of the bla_KPC–2-carrying CRKP rose rapidly from 3.3 to 28.1%. KPC-2-producing K. pneumoniae ST11 isolates were dominant in CRKP, which emerged in 2012 and caused several outbreaks. Most isolates exhibited multidrug-resistant to commonly used antibiotics, while all the isolates remained susceptible to tigecycline and polymyxin B. The single nucleotide polymorphism (SNP) analysis showed that all these 152 KPC-2-producing K. pneumoniae ST11 isolates could be divided into three genetically distinct clades (A, B, and C) and eleven subclades (A1–A9 and B1–B2). The majority belonged to clade A with KL47 serotype (n = 117, 77.0%), while KL64 and KL16 were identified in clades B and C, respectively. The bla_KPC–2-carrying plasmids exhibited diverse types, namely, IncFII (pHN7A8)/IncR(6/15), IncFII (pHN7A8)/Inc_pA1763–KPC (5/15), IncFII (pHN7A8) (1/15), IncR (1/15), and Inc_pA1763–KPC (1/15). The genetic environment of bla_KPC–2 showed nine IS26-based composite transposons, which had a basic core structure ISKpn27-bla_KPC–2-ΔISKpn6. About 27.6% (42/152) isolates co-carried 2 to 4 virulence marker genes (namely, peg344, iucA, iroB, rmpA, and rmpA2) for hvKp strains. At least three isolates were identified to harbor virulence gene-carrying plasmids.

Conclusion: KPC-2-producing K. pneumoniae ST11 was highly heterogeneous in our hospital. Transmission of these strains was mainly mediated by twelve high-risk clones. The bla_KPC–2-carrying plasmids and genetic environment of bla_KPC–2 genes exhibited active evolution in K. pneumoniae ST11. More attention should be paid to the tendency of KPC-2-ST11 to acquire hypervirulent plasmids.

Introduction

The increasing prevalence of carbapenem-resistant Klebsiella pneumoniae (CRKP) poses a continuously increasing threat to public health worldwide (Munoz-Price et al., 2013). The production of carbapenemases, namely, K. pneumoniae carbapenemase (KPC), NDM, and OXA-48, is the most important mechanism among carbapenem-resistant K. pneumoniae (Lee et al., 2016). Among these carbapenemases, KPC is the most prevalent enzyme since it was reported in 1996 (Nordmann and Poirel, 2019). In China, about 60% of CRKP strains produce the KPC and the KPC-2 is the major variant of KPC family (Hu et al., 2020). Transmission of the KPC-2-coding gene involves multiple mechanisms ranging from clonal spread to horizontal transfer mediated by plasmids and other transposons (Naas et al., 2008).

The spread of KPC-producing K. pneumoniae is primarily associated with different sequence types (STs). ST258 is most common in North America, Latin America, and Europe, while ST11 is most common in Asia. ST11 is a single-locus variant (tonB) of ST258 and both of them belong to clone complex CC258. In China, the mass dissemination of KPC-2-producing CRKP has been primarily restricted to clone ST11 (Yang et al., 2013; Cheng et al., 2016; Shen et al., 2016; Wang et al., 2016; Gu et al., 2019; Li et al., 2019). The genomic analysis has revealed that ST11 strains were highly heterogeneous and can be grouped into two or three genetic lineages according to the capsule types, while KL47 and KL64 are the main types, which have caused outbreaks in many districts (Dong et al., 2018b; Zhao et al., 2020).

The pandemic spread of bla_KPC–2 among ST11 CRKP in China is mainly associated with horizontal transfer mediated by incompatibility group F (IncF) plasmids. The IncFII replicon often coexists with one or more other replicons such as IncFIA (Shen et al., 2016), IncFIB (Wang et al., 2016), and IncR (Yang et al., 2013). The IncpA1763-KPC replicon is a new type of replicon and was designated in 2019 (Qu et al., 2019). Plasmids with multireplicon have a better chance to replicate in a broader range of hosts by encoding different sets of replication initiation proteins (Li et al., 2019). There were also some new types of plasmids carrying bla_KPC–2 discovered.

Genetic structures surrounding the bla_KPC–2 gene are mostly associated with a Tn3-based transposon, Tn4401, with bla_KPC–2 sandwiched between insertion sequences, namely, ISKpn6 and ISKpn7 (Naas et al., 2008). The genetic environments of bla_KPC–2 from China can be assigned into three main categories: Tn4401, the Tn1722-based unit transposon Tn6296, and the IS26-based truncated Tn6296 (ΔTn6296). Both the Tn6296 and ΔTn6296 contain the ISKpn27-bla_KPC–2-ISKpn6 core structure (Zhang et al., 2017). The ΔTn6296 has one or two IS26 elements at one or both ends, while different lengths of truncated Tn6296 can be observed from different plasmids.

Due to the rapid evolution of ST11, carbapenem-resistant hypervirulent ST11 has also been consistently reported (Dong et al., 2018b; Zhao et al., 2020), which has led to a growing therapeutic dilemma (Chen et al., 2018; Gu et al., 2018).

In our hospital, a nosocomial outbreak of CRKP was first reported in 2013 and KPC-2-producing CRKP ST11 was the main pathogen of the transmission (Yang et al., 2013). To further characterize the nosocomial transmission and genomic characteristics of this high-risk clone, all the KPC-2-producing CRKP ST11 isolated from the respiratory department were collected from 2010 to 2018 and their genomic characteristics were analyzed.

Materials and Methods

Bacterial Strains

All the non-duplicated K. pneumoniae isolates were consecutively collected from various clinical specimens between 2010 and 2018 at the respiratory department, including four general wards and one intensive care unit (ICU) in our hospital. All the clinical isolates were identified using VITEK 2 GN cards or VITEK MS (bioMérieux SA, Marcy-l’Etoile, France). The isolates that exhibited non-susceptibility to carbapenems were screened by PCR amplification for the common carbapenemase genes, namely, bla_KPC, bla_NDM, bla_OXA–48, and bla_IMP, as described previously (Poirel et al., 2011). Isolates harboring bla_KPC genes were further analyzed.

Antimicrobial Susceptibility Test

The susceptibility of strains was tested by the VITEK 2 or KB method at first. The minimal inhibition concentrations (MICs) of cefotaxime (CTX), ceftazidime (CAZ), piperacillin-tazobactam (TZP), imipenem (IMP), meropenem (MEM), ertapenem (ETP), amikacin (AK), ciprofloxacin (CIP), tigecycline (TGC), sulfamethoxazole/trimethoprim (SXT), and polymyxin B (PB) were measured with broth microdilution method by using the Biofosun^® Gram-negative panels (Biofosun Biotechnology Corporation Ltd., Shanghai, China) for 15 strains of complete genome sequencing. Results were interpreted according to the interpretive standards of the Clinical and Laboratory Standards Institute (M100, 30th edition). Susceptibility to tigecycline and polymyxin B was interpreted according to the 2017 European Committee on Antimicrobial Susceptibility Testing (EUCAST) breakpoints (Supplementary Table 1).

Whole-Genome Sequencing

Genome DNA from overnight cultures was extracted using the DNeasy^® UltraClean^® Microbial Kit (QIAGEN GmbH, 40,724 Hilden, Germany). Genome sequencing was performed with a paired-end library with an average insert size of 350 bp on the HiSeq X Ten sequencer (Illumina, California, CA, United States). The draft genome was assembled into a scaffold. Quality assessment was performed with Fastqc (Version 0.11.8) and all reads score above Q30 were used with follow-up analysis. After removing the adapter, barcode, and trimming of the raw reads, sequences were assembled using SOAP de novo (SOAP version 2.21) with default settings. N50, N90, coverage rate, and scaffold number were used to identify de novo characters.

Data Analysis

Klebsiella multilocus sequence typing¹ databases were performed to identify ST using data from the assembled contigs. Whole-genome sequencing (WGS) data for all the isolates were analyzed using ResFinder and PlasmidFinder tools provided by the Center for Genomic Epidemiology² to identify genes encoding antimicrobial resistance and plasmid replicons using an identity threshold of 95% and a selected minimum length of 60%. Determination of capsular type and wzi alleles was conducted by using Kleborate.³ The presence of the virulence genes or gene clusters was screened by using the reference databank of the Pasteur Institute, Kleborate, and the National Center for Biotechnology Information (NCBI) blast. We screened for genes associated with increased pathogenicity, namely, peg344, iucA (aerobactin siderophore gene), iroB (salmochelin siderophore biosynthesis gene), rmpA, and rmpA2 (regulators of the mucoid phenotype via increased capsule production), which could be markers for hypervirulent strains and demonstrated >0.95% diagnostic accuracy for identifying hypervirulent K. pneumoniae (Russo et al., 2018; Supplementary Table 2).

Phylogenic Analysis

The CSI Phylogeny version 1.4⁴ was used to call single nucleotide polymorphism (SNP) and alignment. The complete genome sequence of ST11 K. pneumoniae strain HS11286 (CP003200.1) was used as a reference. The tree file was visualized by iTOL⁵ and annotated information was edited by iTOL editor v1_1.

Complete Genome Sequencing Analysis

One or two strains from each subclade were selected for subsequent single-molecule real-time (SMRT) sequencing using a sheared DNA library with an average size of 10 kb on the PacBio RS II (Pacific Biosciences, Menlo Park, CA, United States), and a paired-end library with an average insert size of 350 bp on an Illumina NovaSeq platform (Illumina, California, CA, United States). Preliminary assembly was performed using the Hierarchical Genome Assembly Process 3 (HGAP3) within SMRT link version 5.0.1. The paired-end short Illumina reads were used to correct the polymer errors for long PacBio reads. Reads were mapped into assembled genome sequences using BWA version 0.5.9. Pilon version 1.13 was subsequently used to polish genome sequences using the obtained alignments. The assemblies yielded a circular chromosome and several plasmids for each strain. Genome sequences were annotated with the Rapid Annotation using Subsystem Technology (RAST) pipeline (Overbeek et al., 2014), screening for insertion sequences in ISfinder (Siguier et al., 2006) and BLAST searches. Alignment of bla_KPC–2-carrying plasmid complete sequences of each strain used BRIG and Easyfig.

GenBank Accession Number

All the genome sequences used in this study were submitted to GenBank with BioProject accession numbers PRJNA782392 and PRJNA835045.

Results

Clinical Characteristics

A total of 1,043 non-duplicated K. pneumoniae isolates were collected from various clinical specimens of the respiratory departments between 2010 and 2018, of which about 20.6% (n = 215) were carbapenem-insensitive CRKPs. The majority of CRKP carried bla_KPC–2 gene (n = 170, 79.1%), bla_OXA–48 (n = 39, 18.1%), and bla_IMP (n = 3, 1.4%). All the 173 bla_KPC–2-carrying strains were isolated from various clinical sources from 170 unique patients. Almost half (49.7%) of CRKP was isolated from the respiratory intensive care unit (RICU). The most common sources were sputum (74.6%), urine (14.5%), and blood (5.2%). From 2012 to 2018, the detection rate of the bla_KPC–2-carrying CRKP rose rapidly from 3.3 to 28.1%. Most isolates exhibited multidrug-resistant to commonly used antibiotics, while all the isolates remained susceptible to tigecycline and polymyxin B (Supplementary Table 1).

Multilocus Sequence Typing Analysis

Among the 170 bla_KPC–2-carrying CRKP, six STs were identified. Approximately, 89.4% (n = 152) of them belonged to ST11, followed by ST15 (n = 9), ST395 (n = 5), ST392 (n = 2), ST163 (n = 1), and ST13 (n = 1). CRKP ST11 was first discovered in January 2012 and quickly became the main epidemic clone in our hospital (Figure 1).

FIGURE 1

Figure 1. The number of Klebsiella pneumoniae carbapenemase-2 (KPC-2) producing carbapenem-resistant K. pneumoniae (CRKP) strains isolated each year.

Core-Genome Phylogenetic Analyses of Klebsiella pneumoniae Carbapenemase-2 Producing Carbapenem-resistant Klebsiella Pneumoniae Sequence Type 11

A maximum-likelihood phylogenetic tree was constructed based on core genome single nucleotide polymorphisms (cgSNPs) of the 152 KPC-2-producing CRKP ST11. Three genetically distinct populations, clade A (n = 122), clade B (n = 22), and clade C (n = 8), were observed (Figure 2 and Supplementary Table 2). The isolates of clade A and clade B were further divided into nine and two subclades, respectively. Clade A emerged in 2012 and continued to spread among the different departments until 2018. Clade B emerged in 2016 and spread through 2018. Clade C was associated with an outbreak that occurred between 2014 and 2016 (Figure 3).

FIGURE 2

Figure 2. Core genome single nucleotide polymorphism (SNP) phylogeny of 152 KPC-2-producing sequence type 11 (ST11) isolates. Dotted circles and ovals indicate the three major clades (clades A, B, and C). Letters depict subclades described in the text.

FIGURE 3

Figure 3. The distribution of KPC-2-producing CRKP ST11. ICU, intensive care unit; D1 to D4, four wards of the respiratory department. Clinical isolates with different types are indicated by various symbols (▲: A1; ▼: A2; ◆: A3; ✱: A4; ◼: A5; ●: A6; ❤: A7; ★: A8; ✿: A9; ○: B1; □: B2; and #: C).

Capsular Polysaccharide Synthesis Gene Analysis

In this study, the majority of clade A belonged to KL47 (n = 117) and the rest strains were identified as KL25 (n = 5). Isolates of clade B and clade C belonged to KL64 and KL16, respectively (Supplementary Table 2).

Analysis of bla_KPC–2-Carrying Plasmids

The complete genome sequences of fifteen strains were analyzed. All the bla_KPC–2 genes were located on the plasmids. The characteristics of bla_KPC–2-carrying plasmids are given in Table 1. There were two copies of bla_KPC–2 in tandem on the plasmid in pKPC-IR12079 and pKPC-IR5077.

TABLE 1

Table 1. Characteristics of bla_KPC–2-carrying plasmids.

Genetic Environment of bla_KPC–2

The genetic environment of bla_KPC–2 showed nine distinct variants compared with the Tn6296 structure (accession no. FJ628167). All the nine forms were IS26-based composite transposons that had a basic core structure (ISKpn27-bla_KPC–2-ΔISKpn6). These variants presented different lengths and were truncated by IS26 at one or both ends of Tn6296 in different positions. The bla_KPC–2 was commonly present as a single copy on a plasmid, while two copies of bla_KPC–2 were located on a single plasmid of pKPC-IR12079 and pKPC-IR5077, respectively. Meanwhile, two copies of the backbone structure (IS26-ΔtnpR-ISKpn27-bla_KPC–2-ΔISKpn6-klcA-korC-ΔrepB-IS26) were arranged in tandem with pKPC-IR12079, while they shared one IS26 at the connection in pKPC-IR5077 (Figure 4).

FIGURE 4

Figure 4. Linear comparison of bla_KPC–2 genetic surroundings between the 15 plasmids. Genes are denoted by arrows. Resistance genes, mobile elements, and other features are shown based on functional classification. Two copies of the IS26-ΔtnpR-ISKpn27-bla_KPC–2-ISKpn6-klcA-korC-ΔrepB-IS26 units are arranged in tandem in pKPC-IR12079, while in pKPC-IR5077 sharing one IS26 at the connection.

Analysis of Virulence Genes-Carrying Plasmids

About 27.6% (42/152) isolates, including all the subclade A6 isolates (34/34), 6 of 9 subclade B2 isolates, and 2 of 8 subclade A1 isolates, co-carried 2 to 4 virulence marker genes (namely, peg344, iucA, iroB, rmpA, and rmpA2) for hvKp strains. Three isolates from these subclades: IR12061 (B2), IR5077 (A1), and IR12079 (A6) were analyzed to identify the characteristics of these virulence-carrying plasmids (Table 2). The plasmids pVir-IR12061 and pVir-IR5077 showed high nucleotide sequence identity of more than 99% over the 75 and 65% query coverage with pK2044. Both the plasmids got a resistant fragment based on pK2044 (Figure 5). The pVir-IR12079 from subclade A6 harbored 2 virulence genes and showed 34% query coverage with pK2044.

TABLE 2

Table 2. Characteristics of virulence plasmids.

FIGURE 5

Figure 5. Linear comparison of virulence plasmids pVir-IR12061_1, pVir-IR5077_3, and pK2044. Virulence genes are shown by dark arrows. Other genes are denoted by gray arrows. Regions of synteny between adjacent schematics are indicated by the gray-shaded areas.

Discussion

In this study, we retrospectively investigated the characteristics of KPC-2-producing CRKP ST11 isolated from the respiratory departments over 9 years. We chose the respiratory departments consisting of four ordinary departments and an ICU, which were one of the most common CRPK outbreaking areas in our hospital. The molecular characteristics of the CRKP in these departments could to some extent elucidate the relevant mechanisms related to the emergence, outbreak, and persistence of CRKP in the whole hospital.

Consistent with other reports in China (Zhang et al., 2017; Dong et al., 2018b), KPC-2-producing CRKP ST11 was first discovered in January 2012 and then quickly became the main epidemic clone (Figure 1). The genetic basis of the predominant prevalence of ST11 in China was not well understood. Some studies revealed that ST11 clone may possess features such as a type IV secretion system in ICEKp258.1 and chromosomal recombination, which may potentially contribute to the dissemination of ST11 (Chin et al., 2018). So, it is important to identify the intrinsic factors that caused the rapid dissemination of ST11 and focus on CRKP infectious control and surveillance. Unlike ST258 associated both with KPC-2 and KPC-3(Marsh et al., 2019), ST11 was strictly associated with KPC-2 production. Further study was needed to understand the mechanism of the association such as the adaptability of plasmids encoding KPC-2 in ST11.

In concordance with previous reports, all the KPC-2-producing ST11 isolates were resistant to multiple antibiotics, namely, aminoglycosides, fluoroquinolones, and cephalosporins (Wang et al., 2016). Fortunately, all the isolates in this study were susceptible to tigecycline and colistin, suggesting a choice for clinical treatment.

Phylogenetic analysis demonstrated that KPC-2-producing CRKP ST11 could be divided into three clades and eleven subclades, which suggested that the strains were highly heterogeneous. The spread of these 12 high-risk clones caused several persistent outbreaks in the whole respiratory departments during the 9 years (Figures 2, 3). A hospital in Beijing analyzed the CRKP outbreak in the ICU within 14 months and found that there were three clones (Li et al., 2020) and analysis of ST11 strains in a hospital in Shanghai found twelve clones from 2009 to 2016 (Liu et al., 2018). Even during an outbreak, there might be one or two clones (Chen et al., 2019). All these results indicated that ST11 KPC-2-producing CRKP evolved separately in distinct clonal groups.

CPS of CRKP is a hotspot of genome recombination and plays an important role in the evolution of bacterial resistance, which is possibly related to the spread and epidemic of CRKP. In this study, basically three noticeably divergent clades were coherent with the results of their KL type. Thus, KL typing may be an important factor for the subgroup of isolates belonging to the same ST type. CRKP KL47 had been widely reported in China in recent years, mainly in the ST11 type (Dong et al., 2018b; Hu et al., 2020). CRKP KL47 is not a highly virulent strain, but in recent years, many reports have found that KL47-type strains can also acquire virulence plasmids and become highly virulent strains with high lethality (Gu et al., 2018). In this study, it was also found that the CRKP KL47-A6 strain carried rmpA2 and iucA genes. In some studies, the CRKP KL64 strain was found to be highly virulent and the isolation rate had gradually increased in recent years (Zhao et al., 2020; Zhou et al., 2020). In this study, although the CRKP KL64 strain did not account for a large proportion, there were also highly virulent strains, which possessed all the virulence genes. It is necessary to continuously monitor the epidemic situation for these clones. Among these clades, we noticed that clade C isolates (KL16) caused a small outbreak in ICU between 2015 and 2016 and then disappeared in the departments. Most of these isolates co-possessed bla_KPC–2 and bla_OXA–48. We analyzed these isolates in another study and found that those plasmids also showed a divergent evolution (Wang et al., 2021).

Compared with the background of bla_KPC–2 gene in Tn6296, this study exhibited a considerable diversity of structures surrounding bla_KPC–2 due to the IS26 insertion; however, they had a similar core structure, ISKpn27-bla_KPC–2-ISKpn6, as reported in many other studies (Zhang et al., 2017; Dong et al., 2018b). Form IV was found on seven different plasmids, namely, dual replicon IncFII (pHN7A8)/IncR, IncFII (pHN7A8)/Inc_pA1763–KPC, and IncR plasmids, which belonged to the distinct clones and was identical to ΔTn6296-3 in p12181-KPC (Feng et al., 2017). These results suggested that horizontal transfer of Tn6296 might be the major mechanism responsible for the emergence and rapid transmission of bla_KPC–2. There is no clear relationship between the genetic environment and the plasmid types or the clonal profile of the isolates. After transfer to a plasmid, Tn6296 evolved continuously. There were two copies of ΔTn6296 arranged in tandem in pKPC-IR12079 and pKPC-IR5077, which was similar to the three copies of the IS26-ISKpn27-bla_KPC–2-ISKpn6-IS26 unit in a study by Feng et al. (2019). We also found that most of the units were truncated by the IS26 insertion at different sites and there were seven distinct forms ofΔTn6296 (Figure 4). Previous studies also had shown that IS26 is involved in the dissemination and amplification of bla_KPC (He et al., 2015) and might generate tandem repeats via homologous recombination.

In this study, bla_KPC–2-bearing plasmids were structurally divergent, as various types and sizes were detectable in the strains. Among these, 12 out of 14 plasmids from different subclones harbored IncFII_pHN7A8 replicon (Table 1). A recent report indicated that the IncFII-like plasmid may promote the spread of the bla_KPC gene in K. pneumoniae ST11 in China (Fu et al., 2019). IncFII plasmids can be divided into multiple subgroups, namely, IncFII_Y, IncFII_K, IncFII_pHN7A8, and IncFII_p0716–KPC (Qu et al., 2019). They are usually large (>100 kb), but with low copy numbers and often carry an additional replicon type to initiate replication (Ogbolu et al., 2013). In this study, dual replicons such as IncFII_pHN7A8 replicon/IncR and IncFII_pHN7A8/Inc_pA1763–KPC had also been reported in other studies in China (Feng et al., 2019; Li et al., 2022). Inc_pA1763–KPC replicon often coexists with an auxiliary replicon with the major replicons (such as IncFII and IncR) to compose multireplicon plasmids, represented by plasmids p0716-KPC (Feng et al., 2017). As the sole replicon in the plasmid, the Inc_pA1763–KPC replicon is also capable of maintaining the genetic stability of the plasmid (Qu et al., 2019; Chen et al., 2021). In this study, 5 of 14 plasmids harbored dual replicons IncFII_pHN7A8/Inc_pA1763–KPC, while one plasmid harbored sole Inc_pA1763–KPC replicon. Although some plasmids harbored the same replicon, the bla_KPC-harboring plasmids in different subclades showed different sizes and carried different kinds of resistant genes, suggesting the diversity of the bla_KPC-harboring plasmids.

The hypervirulent KPC-2-producing K. pneumoniae ST11 was first identified in 2018 (Gu et al., 2018). The genetic mechanism by which classic ST11 CRKP strains acquired hypervirulence relies on the horizontal acquisition of a pLVPK-like virulence plasmid. In this study, all the hypervirulent isolates obtained a plasmid containing both the virulence genes and resistant genes. Multidrug-resistant (MDR)-virulent plasmids were rare, but an increasing number of studies had been reported in the past few years (Bialek-Davenet et al., 2014; Dong et al., 2018a; Wyres et al., 2020). The emergence and spread of MDR-virulent plasmids as a result of ongoing recombination was a significant potential health threat. Findings in this study provided evidence for active plasmid evolution in K. pneumoniae and suggested an urgent need for the surveillance of multidrug-resistant and hypervirulent K. pneumoniae.

Conclusion

In conclusion, using WGS-based analysis, we described a population snapshot of CRKP ST11 collected from the respiratory department between 2010 and 2018 in our hospital. Different clones harbored different plasmids, while diverse genetic environments of bla_KPC–2 gene suggested different sources of acquisition. This finding demonstrated that hypervirulent and multidrug-resistant K. pneumoniae strains had already appeared several years ago.

This study has some limitations. First, a collection of about 152 strains across a long period in the respiratory departments may not reflect the whole evolution of the KPC-2-producing ST11 CRKP in a hospital. Second, the complete sequence of one plasmid from each subclade may not reflect the characteristics of all the plasmids from the subclade. Third, the strains are mainly isolated from sputum and urine, which are not easy to check whether these strains cause infection or colonization. Finally, we did not phenotypically evaluate the virulence potential of K. pneumoniae clones encoding virulence traits.

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/, PRJNA782392; https://www.ncbi.nlm.nih.gov/genbank/, PRJNA835045.

Author Contributions

DS and JY designed the study. LG, LW, QZ, LY, KY, and YM did the phenotypic and genotypic analyses. LG and JY drafted the manuscript. All authors have read and approved the final version of the manuscript.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.929826/full#supplementary-material

Footnotes

1. ^ http://bigsdb.pasteur.fr/klebsiella/
2. ^ http://www.genomicepidemiology.org/services/
3. ^ https://github.com/katholt/Kleborate
4. ^ https://cge.food.dtu.dk/services/CSIPhylogeny/
5. ^ https://itol.embl.de

References

Bialek-Davenet, S., Criscuolo, A., Ailloud, F., Passet, V., Jones, L., Delannoy-Vieillard, A. S., et al. (2014). Genomic definition of hypervirulent and multidrug-resistant Klebsiella pneumoniae clonal groups. Emerg. Infect. Dis. 20, 1812–1820. doi: 10.3201/eid2011.140206

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Zhang, Y., Yu, S. L., Zhou, Y., Yang, S. Y., Jin, J. L., et al. (2019). Tracking carbapenem-producing Klebsiella pneumoniae Outbreak in an intensive care unit by whole genome sequencing. Front. Cell. Infect. Microbiol. 9:281. doi: 10.3389/fcimb.2019.00281

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. M., Guo, M. K., Ke, S. C., Lin, Y. P., Li, C. R., Vy Nguyen, H. T., et al. (2018). Emergence and nosocomial spread of ST11 carbapenem-resistant Klebsiella pneumoniae co-producing OXA-48 and KPC-2 in a regional hospital in Taiwan. J. Med. Microbiol. 67, 957–964. doi: 10.1099/jmm.0.000771

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Zhan, Z., Jiang, X., Qing, Y., Yin, Z., Mei, L., et al. (2021). Comparative genomic analyses of Inc(pA1763-KPC) plasmids. J. Basic Microbiol. 61, 219–229. doi: 10.1002/jobm.202000668

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, L., Cao, X. L., Zhang, Z. F., Ning, M. Z., Xu, X. J., Zhou, W., et al. (2016). Clonal dissemination of KPC-2 producing Klebsiella pneumoniae ST11 clone with high prevalence of oqxAB and rmtB in a tertiary hospital in China: results from a 3-year period. Ann. Clin. Microbiol. Antimicrob. 15:1. doi: 10.1186/s12941-015-0109-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chin, C. Y., Tipton, K. A., Farokhyfar, M., Burd, E. M., Weiss, D. S., and Rather, P. N. (2018). A high-frequency phenotypic switch links bacterial virulence and environmental survival in Acinetobacter baumannii. Nat. Microbiol. 3, 563–569. doi: 10.1038/s41564-018-0151-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, N., Zhang, R., Liu, L., Li, R., Lin, D., Chan, E. W., et al. (2018b). Genome analysis of clinical multilocus sequence Type 11 Klebsiella pneumoniae from China. Microb. Genom. 4:e000149. doi: 10.1099/mgen.0.000149

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, N., Lin, D., Zhang, R., Chan, E. W., and Chen, S. (2018a). Carriage of blaKPC-2 by a virulence plasmid in hypervirulent Klebsiella pneumoniae. J. Antimicrob. Chemother. 73, 3317–3321.

Google Scholar

Feng, J., Yin, Z., Zhao, Q., Zhang, D., Jiang, X., Wu, W., et al. (2017). Genomic characterization of novel IncFII-type multidrug resistant plasmids p0716-KPC and p12181-KPC from Klebsiella pneumoniae. Sci. Rep. 7:5830.

Google Scholar

Feng, Y., Liu, L., McNally, A., and Zong, Z. (2019). Coexistence of three blaKPC-2 genes on an IncF/IncR plasmid in ST11 Klebsiella pneumoniae. J. Glob. Antimicrob. Resist. 17, 90–93.

Google Scholar

Fu, P., Tang, Y., Li, G., Yu, L., Wang, Y., and Jiang, X. (2019). Pandemic spread of bla((KPC-2)) among Klebsiella pneumoniae ST11 in China is associated with horizontal transfer mediated by IncFII-like plasmids. Int. J. Antimicrob. Agents 54, 117–124.

Google Scholar

Gu, B., Bi, R., Cao, X., Qian, H., Hu, R., and Ma, P. (2019). Clonal dissemination of KPC-2-producing Klebsiella pneumoniae ST11 and ST48 clone among multiple departments in a tertiary teaching hospital in Jiangsu Province, China. Ann. Transl. Med. 7:716. doi: 10.21037/atm.2019.12.01

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

He, S., Hickman, A. B., Varani, A. M., Siguier, P., Chandler, M., Dekker, J. P., et al. (2015). Insertion Sequence IS26 reorganizes plasmids in clinically isolated multidrug-resistant bacteria by replicative transposition. mBio 6:e00762-15. doi: 10.1128/mBio.00762-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Y., Liu, C., Shen, Z., Zhou, H., Cao, J., Chen, S., et al. (2020). Prevalence, risk factors and molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae in patients from Zhejiang, China, 2008-2018. Emerg. Microbes Infect. 9, 1771–1779.

Google Scholar

Lee, C. R., Lee, J. H., Park, K. S., Kim, Y. B., Jeong, B. C., and Lee, S. H. (2016). Global dissemination of Carbapenemase-Producing Klebsiella pneumoniae: epidemiology, genetic context, treatment options, and detection methods. Front. Microbiol. 7:895. doi: 10.3389/fmicb.2016.00895

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Jiang, X., Yang, T., Ju, Y., Yin, Z., Yue, L., et al. (2022). Genomic epidemiology of carbapenemase-producing Klebsiella pneumoniae in china. Genomics Proteomics Bioinformatics. [Epub ahead of print]. doi: 10.1016/j.gpb.2022.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D. X., Zhai, Y., Zhang, Z., Guo, Y. T., Wang, Z. W., He, Z. L., et al. (2020). Genetic factors related to the widespread dissemination of ST11 extensively drug-resistant carbapenemase-producing Klebsiella pneumoniae strains within hospital. Chin. Med. J. 133, 2573–2585. doi: 10.1097/CM9.0000000000001101

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Shen, H., Zhu, C., and Yu, Y. (2019). Carbapenem-Resistant Klebsiella pneumoniae Infections among ICU Admission Patients in Central China: prevalence and prediction model. Biomed Res. Int. 2019:9767313. doi: 10.1155/2019/9767313

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Yu, J., Chen, F., Yu, J., Simner, P., Tamma, P., et al. (2018). Emergence and establishment of KPC-2-producing ST11 Klebsiella pneumoniae in a general hospital in Shanghai, China. Eur. J. Clin. Microbiol. Infect. Dis. 37, 293–299. doi: 10.1007/s10096-017-3131-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Marsh, J. W., Mustapha, M. M., Griffith, M. P., Evans, D. R., Ezeonwuka, C., Pasculle, A. W., et al. (2019). Evolution of Outbreak-Causing Carbapenem-Resistant Klebsiella pneumoniae ST258 at a Tertiary Care Hospital over 8 Years. mBio 10:e01945-19. doi: 10.1128/mBio.01945-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Munoz-Price, L. S., Poirel, L., Bonomo, R. A., Schwaber, M. J., Daikos, G. L., Cormican, M., et al. (2013). Clinical epidemiology of the global expansion of Klebsiella pneumoniae carbapenemases. Lancet Infect. Dis. 13, 785–796.

Google Scholar

Naas, T., Cuzon, G., Villegas, M. V., Lartigue, M. F., Quinn, J. P., and Nordmann, P. (2008). Genetic structures at the origin of acquisition of the beta-lactamase bla KPC gene. Antimicrob. Agents Chemother. 52, 1257–1263. doi: 10.1128/AAC.01451-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Nordmann, P., and Poirel, L. (2019). Epidemiology and diagnostics of Carbapenem Resistance in Gram-negative Bacteria. Clin. Infect. Dis. 69, S521–S521. doi: 10.1093/cid/ciz824

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogbolu, D. O., Daini, O. A., Ogunledun, A., Terry Alli, O. A., and Webber, M. A. (2013). Dissemination of IncF plasmids carrying beta-lactamase genes in Gram-negative bacteria from Nigerian hospitals. J. Infect. Dev. Ctries. 7, 382–390. doi: 10.3855/jidc.2613

PubMed Abstract | CrossRef Full Text | Google Scholar

Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 42, D206–D214. doi: 10.1093/nar/gkt1226

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirel, L., Walsh, T. R., Cuvillier, V., and Nordmann, P. (2011). Multiplex PCR for detection of acquired carbapenemase genes. Diagn. Microbiol. Infect. Dis. 70, 119–123. doi: 10.1016/j.diagmicrobio.2010.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, D., Shen, Y., Hu, L., Jiang, X., Yin, Z., Gao, B., et al. (2019). Comparative analysis of KPC-2-encoding chimera plasmids with multi-replicon IncR:Inc(pA1763-KPC):IncN1 or IncFII(pHN7A8):Inc(pA1763-KPC):IncN1. Infect. Drug Resist. 12, 285–296. doi: 10.2147/IDR.S189168

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Shen, P., Zhang, Y., Tang, Y., Liang, W., and Jiang, X. (2016). Molecular dissection of blaKPC-2-bearing plasmids evolving in Klebsiella pneumoniae isolated at one teaching hospital in Shanghai, China. FEMS Microbiol. Lett. 363:fnw142. doi: 10.1093/femsle/fnw142

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Guo, L., Ye, K., and Yang, J. (2021). Genetic characteristics of OXA-48-producing Enterobacterales from China. J. Glob. Antimicrob. Resist. 26, 285–291. doi: 10.1016/j.jgar.2021.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L. H., Wei, D. D., Wan, L. G., Yu, Y., Deng, Q., and Liu, Y. (2016). Diversity of the genetic environment of the blaKPC-2 Gene among Klebsiella pneumoniae clinical isolates in a Chinese Hospital. Microb. Drug Resist. 22, 15–21. doi: 10.1089/mdr.2014.0281

PubMed Abstract | CrossRef Full Text | Google Scholar

Wyres, K. L., Nguyen, T., Lam, M., Judd, L. M., van Vinh Chau, N., Dance, D. A. B., et al. (2020). Genomic surveillance for hypervirulence and multi-drug resistance in invasive Klebsiella pneumoniae from South and Southeast Asia. Genome Med. 12:11. doi: 10.1186/s13073-019-0706-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Ye, L., Guo, L., Zhao, Q., Chen, R., Luo, Y., et al. (2013). A nosocomial outbreak of KPC-2-producing Klebsiella pneumoniae in a Chinese hospital: dissemination of ST11 and emergence of ST37, ST392 and ST395. Clin. Microbiol. Infect. 19, E509–E515. doi: 10.1111/1469-0691.12275

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Liu, L., Zhou, H., Chan, E. W., Li, J., Fang, Y., et al. (2017). Nationwide Surveillance of Clinical Carbapenem-resistant Enterobacteriaceae (CRE) Strains in China. EBioMedicine 19, 98–106. doi: 10.1016/j.ebiom.2017.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Liu, C., Liu, Y., Zhang, Y., Xiong, Z., Fan, Y., et al. (2020). Genomic characteristics of clinically important ST11 Klebsiella pneumoniae strains worldwide. J. Glob. Antimicrob. Resist. 22, 519–526.

Google Scholar

Zhou, K., Xiao, T., David, S., Wang, Q., Zhou, Y., Guo, L., et al. (2020). Novel Subclone of Carbapenem-Resistant Klebsiella pneumoniae sequence Type 11 with enhanced virulence and transmissibility, China. Emerg. Infect. Dis. 26, 289–297. doi: 10.3201/eid2602.190594

PubMed Abstract | CrossRef Full Text | Google Scholar