Genetic Characterization of Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates in a Tertiary Hospital in Greece, 2018–2022
by
, , and
Charalampos Zarras
1,2,
Theodoros Karampatakis
1,*,
Styliani Pappa
1,
Elias Iosifidis
3,
Eleni Vagdatli
2,
Emmanuel Roilides
3 and
Anna Papa
1 1
Department of Microbiology, Medical Faculty, School of Health Sciences, Aristotle University of Thessaloniki, 541 24 Thessaloniki, Greece
2
Microbiology Department, Hippokration General Hospital, 546 42 Thessaloniki, Greece
3
Infectious Disease Unit, 3rd Department of Pediatrics, Medical Faculty, School of Health Sciences, Hippokration General Hospital, 546 42 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Antibiotics 2023, 12(6), 976; https://doi.org/10.3390/antibiotics12060976
Submission received: 4 May 2023
/
Revised: 23 May 2023
/
Accepted: 25 May 2023
/
Published: 28 May 2023
(This article belongs to the Special Issue Molecular Characterization of Gram-Negative Bacteria: Antimicrobial Resistance, Virulence and Epidemiology)
Abstract
:Background: Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a serious public health issue. The study aimed to identify the antimicrobial resistance and accessory genes, the clonal relatedness, and the evolutionary dynamics of selected CRKP isolates recovered in an adult and pediatric intensive care unit of a tertiary hospital in Greece. Methods: Twenty-four CRKP isolates recovered during 2018–2022 were included in the study. Next-generation sequencing was performed using the Ion Torrent PGM Platform. The identification of the plasmid content, MLST, and antimicrobial resistance genes, as well as the comparison of multiple genome alignments and the identification of core genome single-nucleotide polymorphism sites, were performed using various bioinformatics software. Results: The isolates belonged to eight sequence types: 11, 15, 30, 35, 39, 307, 323, and 512. A variety of carbapenemases (KPC, VIM, NDM, and OXA-48) and resistance genes were detected. CRKP strains shared visually common genomic regions with the reference strain (NTUH-K2044). ST15, ST323, ST39, and ST11 CRKP isolates presented on average 17, 6, 16, and 866 recombined SNPs, respectively. All isolates belonging to ST15, ST323, and ST39 were classified into distinct phylogenetic branches, while ST11 isolates were assigned to a two-subclade branch. For large CRKP sets, the phylogeny seems to change approximately every seven SNPs. Conclusions: The current study provides insight into the genetic characterization of CRKP isolates in the ICUs of a tertiary hospital. Our results indicate clonal dispersion of ST15, ST323, and ST39 and highly diverged ST11 isolates.
1. Introduction
Klebsiella pneumoniae is recognized as an important Gram-negative opportunistic pathogen that causes community- and hospital-acquired infections globally [1]. It can cause ventilator-associated pneumonia among patients in intensive care units (ICUs), bloodstream infections, and urinary tract infections [2,3,4]. The emergence of multidrug-resistant (MDR) K. pneumoniae has become a major public health problem, causing a significant increase in morbidity and mortality worldwide [5]. The spread of carbapenem-resistant K. pneumoniae (CRKP) has generally been considered an increasingly serious issue for clinical practice due to the limitation of therapeutic options [6,7,8]. Greece is endemic for CRKP [9,10,11,12,13]; data from the European Antimicrobial Resistance Surveillance Network committee for the year 2020 indicate that the prevalence of CRKP in Greece was 66.3%, the highest among the European countries [14].
The rapid dissemination of CRKP is closely related to the antimicrobial resistance genes carried by plasmids and transferable genetic elements [15,16,17,18]. CRKP strains accumulate antimicrobial resistance genes due to inadequate implementation of infection control measures in healthcare settings and the irrational use of antimicrobials, resulting in the emergence of MDR, extremely drug-resistant, and pan-drug-resistant phenotypes [19]. The main mechanism of resistance to carbapenems is through the production of enzymes called carbapenemases, with K. pneumoniae carbapenemase (KPC) being the most prevalent, followed by Verona integron-encoded metallo-β-lactamase (VIM), New Delhi metallo-β-lactamase (NDM), and OXA-48 carbapenemases [20,21,22]. The contribution of outer membrane protein (Omp) deficiencies, such as OmpK35, OmpK36, and OmpK37, is considered a secondary mechanism in the emergence of CRKP [23,24,25]. The pathogenic potential of these strains is determined by virulence factors such as capsular polysaccharide synthesis, fimbriae, pili, outer membrane proteins, and iron acquisition systems [26].
The application of next-generation sequencing (NGS) has become a powerful tool for obtaining whole genome sequences (WGS), resulting in genetic characterization of the strains and a better understanding of the genomic diversity, providing fast and in-depth information about bacterial pathogenicity [27]. WGS analysis is recommended for genotyping and determining relatedness between clinical K. pneumoniae strains, as well as for the surveillance of virulence genes [28]. The aim of the present study was to identify the antimicrobial resistance and accessory genes of selected CRKP isolates recovered in an adult intensive care unit (ICU) and a pediatric ICU (PICU) of a tertiary hospital in Greece, to evaluate their clonal relatedness through core genome single-nucleotide polymorphism (cgSNP) analysis, and to calculate the evolutionary dynamics of this bacterial population.
2. Results
2.1. Antimicrobial Susceptibility Testing
All CRKP isolates were resistant to β-lactam antibiotics (including penicillins, cephalosporins, monobactams, and carbapenems). Nineteen isolates (79.2%) were resistant to at least one aminoglycoside, 23 isolates (95.8%) showed resistance to quinolones, eight isolates (33.3%) were resistant to fosfomycin and tigecycline, and 21 (87.5%) to cotrimoxazole. Four out of five isolates that carried the class A β-lactamase blaKPC showed susceptibility to ceftazidime-avibactam (CAZ-AVI). The phenotypic resistance rates to antimicrobials of CRKP isolates are shown in the heatmap in Figure 1.
2.2. MLST, Antimicrobial Resistance Genes, Plasmids
The CRKP isolates belonged to eight different STs: ST11 (7, 29.1%), ST39 (7, 29.1%), ST323 (3, 12.5%), ST15 (7, 29.1%), ST307 (1, 4.2%), ST30 (1, 4.2%), ST512 (1, 4.2%), and ST35 (1, 4.2%).
The resistance genes and the plasmid content per isolate are detailed in Table 1. The prevalence of genes conferring resistance to beta-lactams was 100%, to aminoglycosides 79.2%, to quinolones 95.8%, to fosfomycin and tigecycline 33.3%, and to cotrimoxazole 87.5%. The most prevalent plasmid was IncF, which is associated with the spread of several extended-spectrum β-lactamase (ESBL) genes or carbapenemases and was detected in all CRKP strains; additional plasmids were found in fifteen strains, such as Col and IncC.
2.3. Virulence Factors and Efflux and Regulator Systems
The virulence factors detected included fimbrial genes, as 23/24 isolates possessed the mrk cluster (type 3 fimbriae); siderophores, as yersiniabactin genes (ybt, irp, and fyuA) were present in 17 isolates; and genes encoding aerobactin siderophore (iut) were identified in five isolates (three of them also carried yersiniabactin genes). In addition, the kfu gene (Klebsiella ferric iron uptake, which is a regulator of the iron transport system) was detected in four isolates.
Efflux and regulator system genes encoding the AcrAB efflux pump (acr, mar, sox, rob, ram, sdiA, fis, and envR) or the OqxAB efflux pump (oqx and rar) were present in all isolates (Table 1).
2.4. Genomic Comparison among CRKP Strains
Circular genome maps were generated for the four STs with more than one isolate (ST11, ST39, ST323, and ST15). It was found that isolates belonging to the same ST shared many identical regions with the reference strain. CRKP isolates C4112 and C746 of the ST11, D2452, and D6184 of the ST39, Z866 of the ST323, and B2562 of the ST15 displayed visually the highest similarity with the reference strain throughout their whole genome (Figure 2).
2.5. CgSNP-Based Phylogenetic Analysis
The length of the multiple core genome alignments of all CRKP isolates was 1,892,845 base pairs. The core genome length corresponds to about 36.1% of the reference genome length. Nearly 1.9 million columns in the core genome alignment (almost 98.5%) showed identical nucleotides for all isolates, meaning that 1.5% of the genome is polymorphic.
A phylogenetic tree including the cgSNP-based phylogenetic analysis of all isolates, along with the reference strain, is seen in Figure 3.
When performing all possible pairwise analyses among ST15 isolates, an average of 17 recombined SNPs (range 5–26) were detected; the SNP diversity along the core genome was very low, that is, 0–3 cgSNPs per kilobase. ST323 CRKP strains presented an average of six recombined SNPs (range 5–6), and the SNP diversity was extremely low, that is, 0-2 cgSNPs per kilobase. ST39 CRKP strains displayed an average of 16 recombined SNPs (range 4–43), with the SNP diversity being low, that is, 0–5 cgSNPs per kilobase. Finally, ST11 CRKP strains showed an average of 866 recombined SNPs (range 1–1652), with the SNP diversity being clearly higher, that is, 0–66 cgSNPs per kilobase.
For pairs of CRKP isolates, for instance, C4112 and C746 with relatively low diversity (7.5 × 10−4) (Figure 4G), the impact of recombination is almost directly evident in the pattern of SNP density along the genome. While the SNP density is very low along most of the genome, that is, 0–2 SNPs per kilobase, there are a few segments where the SNP density is much higher (Figure 4A) and comparable to the typical SNP density between randomly selected highly diverged CRKP isolates, that is, 30–65 SNPs per kilobase (Figure 4C). For some pairs of isolates, for example, A18940 and D2856, with higher genetic diversity (1.5 × 10−3) (Figure 4G), the frequency of these high-spiking recombined regions increases (Figure 4B), until most of the genome is covered by such regions when comparing highly divergent isolates (5.5 × 10−3), for instance, A5051 and C2482 (Figure 4C,G). The distribution of SNP densities was described by a majority of clonally inherited regions with mostly up to three SNPs per kilobase and a long tail of recombined regions with up to 50 or 65 SNPs per kilobase (Figure 4D–F).
Figure 5 highlights the distribution of the lengths of tree-compatible stretches. This distribution had a mode at n = 2, and stretches were normally around 7–12 consecutive SNPs and very rarely longer than 15–20 consecutive SNPs (Figure 5A). Similarly, tree-compatible segments were typically just a few hundreds of base pairs long (around 1000–1500 base pairs) and very rarely more than 2000 base pairs (Figure 5B).
Figure 6 underscores the ratio C/S for random subsets of all 24 isolates. When comprising all 24 isolates, a ratio C/S = 0.15 was obtained for the 5% homoplasy-corrected full alignment. For small subsets of isolates, the C/S ratio displayed significant fluctuations. For instance, for subsets of n = 7 isolates, the C/S ratio presented a mean of 0.094 ± 0.046. However, as the number of isolates in the subset increased, the ratios converged to the value C/S = 0.15, and for large subsets of isolates, there was little variation in this ratio. Consequently, for alignments of large sets of isolates, the phylogeny should change at least approximately every seven SNPs (Figure 6). In addition, on average, each randomly chosen position on the core genome has been overwritten at least T = 4344 times (5% homoplasies removed).
3. Materials and Methods
3.1. Clinical Isolates-Setting
The present study included 24 CRKP isolates collected from June 2018 to July 2022 from 18 patients hospitalized in the ICU (8 females and 10 males) and six patients hospitalized in the PICU (4 females and 2 males) of a tertiary university-affiliated general hospital in Greece. The median age of the adult patients was 57.5 years (range 30–88), while the median age of the children was 6.2 years (range 0.25–17). The selection of the isolates was based on the type of carbapenemase(s) previously detected in order to include all possible combinations, as they were pre-screened as previously described [11]. Eight isolates were collected from colonization sites (rectal swabs), while 16 isolates were taken from various types of clinical samples (Table 2).
3.2. Ethics Approval
The study was approved by the Ethics Committee of Aristotle’s University Medical Faculty (no. of approval 5.160/18/12/2019).
3.3. Microbiological Methods-Antimicrobial Susceptibility Testing
CRKP isolates were identified through the automated system VITEK 2 (bioMérieux, Marcy-l’Étoile, France) using the Gram-negative identification card (GN ID). Antimicrobial susceptibility testing (AST) of isolates was performed using the AST318 and XN10 cards, while the use of minimum inhibitory concentration (MIC) test strips (Liofilchem srl, Roseto, Italy) was applied to confirm susceptibility to ceftazidime-avibactam (CAZ-AVI). Identification and susceptibility cards were interpreted according to the manufacturer’s instructions. The Clinical and Laboratory Standards Institute (CLSI) susceptibility breakpoints were applied for the interpretation of results [29]. The breakpoints approved by the US Food and Drug Administration (FDA) were used to interpret the MICs of tigecycline [30].
All CRKP isolates were previously screened for carbapenemase production using a lateral flow immunochromatographic assay (LFIA) (NG-Test CARBA 5, NG Biotech, Paris, France). The results were confirmed by multiplex PCR in a single reaction following a modified protocol that was previously described [11,31].
3.4. DNA Extraction and Whole Genome Sequencing
Bacterial DNA was extracted using the DNA extraction kit (Qiagen, Hilden, Germany). The DNA concentration (sample starting concentration between 10–100 ng/μL) was measured using the Qubit double-strand DNA HS assay kit (Q32851, Life Technologies Corporation, Grand Island, NY, USA).
NGS was performed using the Ion Torrent PGM Platform (Life Technologies Corporation). All procedures regarding purification, ligation, barcoding, library preparation, emulsion PCR, and enrichment were performed following the manufacturer’s instructions. PCR products were loaded on an Ion-316™ chip. The Ion PGM Hi-Q View (200) chemistry was applied using the Ion PGM Hi-Q view sequencing kit (A25592).
3.5. Assembly Assessment and Genome Annotation
The software Geneious Prime version 2021.2.1 and BLAST Ring Image Generator (BRIG) were used to produce genome assemblies and annotation data [32,33]. K. pneumoniae NTUH-K2044 strain (GenBank accession number NC_012731) was used as a reference sequence, which is an invasive and hypervirulent K. pneumoniae strain with the hypervirulent plasmid pK2044 [34].
3.6. MLST and Detection of Antimicrobial Resistance Genes and Plasmids
The sequence types (ST), plasmid types, and antimicrobial resistance genes were identified through the Center for Genomic Epidemiology website using the MLST web server and the related online databases, PlasmidFinder version 2.0 and Resfinder version 4.1, respectively [35,36,37]. The Comprehensive Antibiotic Resistance Database was also used, with the selection criteria for hits set to perfect (100%) and strict (>95%) identity to the reference sequence [38]. The detection of genes related to virulence capsule, efflux, and regulator systems was performed according to the protocols available from the Pasteur Institute (https://bigsdb.pasteur.fr/cgi-bin/bigsdb/bigsdb.pl?db=pubmlst_klebsiella_seqdef (accessed on 20 October 2022).
3.7. Genomic Comparison–Core Genome Single-Nucleotide Polymorphism (cgSNP)-Based Phylogenetic Analysis
The BRIG software was used to visually compare multiple genome alignments using a default minimum threshold of 50% [33]. The reference sequence alignment-based phylogeny builder version 1.12 (REALPHY; https://realphy.unibas.ch/realphy/ (accessed on 15 November 2022) was used to identify relevant SNP sites for core genome phylogenetic analysis [39]. Default input parameters were applied for cgSNP identification. The K. pneumoniae NTUH-K2044 sequence was used as a reference for the multiple alignments and the visual comparison of the whole genomes [34]. The produced alignment was used for the creation of an unrooted tree on the entire core genome through the approximate maximum likelihood (ML) method using PhyML [40]. Pairwise analysis between isolates was also performed using REALPHY, as previously described by Dixit et al. [41]. The calculation of SNPs between isolates was conducted through REALPHY, using a Poisson mixture model plus a negative binomial for the recombined regions. REALPHY was also used to assess to what extent mutually consistent SNPs are clustered along the alignment; the lengths of segments along the alignment that are consistent with a single phylogeny were calculated [39]. In addition, the lower bound for the ratio between the total number of phylogeny changes (C) and substitution events (S) (C/S) that occur along the core genome alignment was calculated. The same ratio was calculated for any subset of CRKP isolates by removing 5% of potentially homoplastic sites [42]. Moreover, a rough estimate for the average number of times T that a randomly chosen position in the core genome alignment has been overwritten by recombination in its history was calculated, that is, the time since the genetic ancestors of the position in the alignment diverged from a common ancestor.
4. Discussion
In the current study, the WGSs of 24 CRKP isolates were analyzed. It was found that they belonged to eight different STs, with ST11 and ST39 being the most prevalent. Ten out of 24 isolates (41.6%) carried blaNDM-1, and two of them (8.3%) transferred concurrently blaNDM-1 and blaOXA-48. Six out of 10 NDM-producers (60.0%) were assigned to ST11, which is common in NDM-positive CRKP strains and is globally distributed, especially in Asia and particularly in China, but also in Europe [43,44]. In Greece, the first outbreak of an NDM-1-producing K. pneumoniae ST11 clonal strain occurred during 2011–2012 [44], while the first report of a K. pneumoniae ST11 clinical isolate co-producing NDM-1 and OXA-48 carbapenemases was reported in 2016 [22]. NDM-1-producing ST11 CRKP strains have been isolated during the same period in other Balkan countries [45]. However, three out of 10 NDM producers (30.0%) were assigned to ST15; this has also been previously described [46].
ST39 was also the most prevalent type, encompassing seven isolates (29.1%); two of them carried blaKPC-2, while one of them (14.3%) harbored blaKPC-33 (a KPC-2 variant). CRKP strains belonging to ST39 have been detected only in sporadic cases so far, like ST39 carrying blaKPC-2 [47] or blaVIM-1 [48]. ST39 isolates transferring KPC-33 have been recently detected in Greece [8]. The present study also included the recently reported four ST39 isolates carrying both blaKPC-2 and blaVIM-1 [21]. Recently, Kuzina et al. highlighted the concurrent isolation of blaNDM-1, blaKPC-2, and blaOXA-48 in ST39 CRKP isolates [49]. Three CRKP isolates (12.5%) were assigned to ST323; sporadic cases of ST323 CRKP have been previously reported in Greece and other European countries [50,51]. Another study from Tunisia reports the isolation of clonally related ST323 CRKP isolates [52]. It is of interest that one isolate, A1746/22, carried the VEB-25 gene (article submitted for publication); VEB-25 is a variant of VEB-1 and confers resistance to CAZ-AVI [53,54].
Twenty of the 24 isolates harbored at least two genes encoding aminoglycoside-modifying enzymes (AMEs), and six of them carried one 16S rRNA methylase gene (armA, rmtB). In addition, 12 of these isolates carried at least one gene of the three AME classes (acetyltransferases-aac, phosphotransferase-aph, nucleotransferases-ant). A disagreement between phenotype and genotype has been detected in five isolates (20.8%) (detection of AME-encoding resistance genes but phenotypically sensitive to both aminoglycosides tested). A disagreement of approximately 10% between genotype and phenotype for aminoglycosides for the VITEK 2 automated system has been reported previously [55].
Fluoroquinolone resistance was observed in seven isolates (29.2%), which carried the plasmid-encoded qnr gene, as previously described for CRKP isolates [56]. The multidrug efflux pump oqxAB gene was detected in all strains [57], while the aac(6′)-Ib-cr gene (which encodes an aminoglycoside acetyltransferase that modifies not only aminoglycosides but also fluoroquinolones) has been identified in 16 isolates, as previously described [56].
For fosfomycin resistance, the gene fosA was detected in all isolates, but 16 of them (66.7%) were phenotypically susceptible, evidenced by an increase in the MIC, as previously described [58]. Resistance genes for cotrimoxazole (sul or dfrA) were found in all isolates except one. The Inc-type plasmids detected in all isolates and the Col-type plasmids detected in 15 isolates have been reported in other studies [21,59].
In terms of adhesive structures, almost all isolates possessed the mrk cluster, characteristic of type 3 fimbriae, which promotes bacterial adhesion to both abiotic and biotic surfaces [60]. Regarding siderophores, the majority of isolates carried yersiniabactin genes that are significantly associated with pathogenesis and invasive infections in humans [61], while aerobactin (detected in five isolates) is more specific for hypervirulent K. pneumoniae [62]. The acquisition of a chimeric plasmid including genes for the siderophores yersiniabactin and aerobactin grants cefiderocol activity against CRKP strains [63].
NGS allows for more accurate analyses regarding the clonal relatedness of strains compared to conventional typing methods [64]. Several studies have described the molecular epidemiology of CRKP isolates through NGS using a gene-by-gene comparison or SNP analysis [65,66].
In pairwise comparisons, isolates of ST15, ST323, and ST39 presented an average of 17, 6, and 16 recombined SNPs, respectively, implying the clonal relatedness of these isolates in accordance with their ST classification. Similarly, Miro et al. report the clonal relatedness of OXA-48-producing ST405 CRKP strains further studied with core genome MLST, showing between 6 and 17 cgSNPs [64]. Moreover, Onori et al. reveal the clonal relatedness of sporadic KPC-producing ST258 CRKP strains, presenting an average of 20 or 27 cgSNPs [67]. ST15 and ST39 strains display a higher number of cgSNPs compared to ST323 strains, and this is highlighted as a slightly increased diversity in the corresponding clades of the phylogenetic tree (Figure 3). However, ST15 and ST39 has been isolated over a 33- and 36-month period, respectively, while ST323 have been collected over a 2.5-month period (Table 2). This could indicate a difference in the measured paces of the molecular clock between these clades, as previously proposed [67].
On the contrary, ST11 isolates showed an average of 866 cgSNPs, indicating the polyclonality of these isolates. This indication is further strengthened by the variability in their content of carbapenemase and antimicrobial resistance genes (Table 1) and by their phylogenetic division into two distinct subclades (Figure 3). The high average of cgSNPs occurred between ST11 strains belonging to these two different subclades, resulting in SNP patterns (Figure 4A). Most probably, these high SNP density patterns result from horizontal recombination events, e.g., the transfer of a plasmid or phage, as previously proposed for other bacterial species [42]. These results underscore the higher discriminatory power of cgSNP analysis compared to MLST [64].
In addition, the current study highlighted novel ways in which these cgSNPs could be used to quantify phylogenetic structures and the role of recombination in genome evolution. It was shown that there are pairs of closely related isolates and that most of their DNA has been clonally inherited; however, most of the changes are derived from recombination events. When focusing on the core genome alignment, it has been shown that, until the current genetic state of CRKP strains, each genomic locus has been overwritten thousands of times by recombination in the history of their evolution. Sakoparnig et al. reached the same conclusions when focusing on E. coli isolates, and they disclosed a similar C/S ratio [42]. Moreover, an attempt was made to calculate the probability for a pair of SNPs to be consistent with a common phylogeny as a function of their genome distance [68]. Several other studies have also attempted to quantify recombination in bacteria and assess the relative rate with which different lineages have recombined, focusing on different types of bacteria [42,68]. To the best of our knowledge, our study is the first to accomplish such a task by studying CRKP strains.
One of the basic challenges among WGS-based typing methods is the definition of a cut-off value for the identification of the clonal relatedness between isolates; however, this remains doubtful, as it seems to be not only species-dependent but also population-dependent [64,67,69]. Thus, a threshold range should be evaluated along with clinical epidemiological data.
In conclusion, the current study provides insight into the genetic characterization of CRKP isolates circulating during a 4-year period in the ICUs of a Greek tertiary hospital. Phylogenetic analysis showed that some of them were epidemiologically related. Given that NGS has become more affordable, its use in hospital settings is extremely helpful, as the obtained data could guide infection control and prevention strategies.
Author Contributions
C.Z.: investigation, review, and editing. T.K.: methodology, conceptualization, formal analysis, writing the initial draft. S.P.: investigation. E.I.: provision of study materials, review, and editing. E.V.: investigation. E.R.: provision of study materials, review, and editing. A.P.: conceptualization, methodology, data curation, formal analysis, visualization, reviewing and editing, supervision, funding acquisition. All authors have read and agreed to the published version of the manuscript.
Funding
The current study was supported by the European Union’s Horizon 2020 projects COMPARE (grant number 643476) and VEO (grant number 874735).
Institutional Review Board Statement
The study was approved by the Ethics Committee of Aristotle’s University Medical Faculty (no. of approval 5.160/18/12/2019).
Informed Consent Statement
Not Applicable.
Data Availability Statement
The data of this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession numbers PRJEB42192 and PRJEB46676.
Acknowledgments
We thank Surbhi Malhotra-Kumar and Minh Nguyen-Ngoc from the University of Antwerp in Belgium for their technical support.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Juan, C.H.; Fang, S.Y.; Chou, C.H.; Tsai, T.Y.; Lin, Y.T. Clinical characteristics of patients with pneumonia caused by Klebsiella pneumoniae in Taiwan and prevalence of antimicrobial-resistant and hypervirulent strains: A retrospective study. Antimicrob. Resist. Infect. Control 2020, 9, 4. [Google Scholar] [CrossRef] [PubMed]
- Martin, R.M.; Bachman, M.A. Colonization, Infection, and the Accessory Genome of Klebsiella pneumoniae. Front. Cell. Infect. Microbiol. 2018, 8, 4. [Google Scholar] [CrossRef] [PubMed]
- Shao, C.; Wang, W.; Liu, S.; Zhang, Z.; Jiang, M.; Zhang, F. Molecular Epidemiology and Drug Resistant Mechanism of Carbapenem-Resistant Klebsiella pneumoniae in Elderly Patients With Lower Respiratory Tract Infection. Front. Public Health 2021, 9, 669173. [Google Scholar] [CrossRef]
- Gupta, A.; Bhatti, S.; Leytin, A.; Epelbaum, O. Novel complication of an emerging disease: Invasive Klebsiella pneumoniae liver abscess syndrome as a cause of acute respiratory distress syndrome. Clin. Pract. 2018, 8, 1021. [Google Scholar] [CrossRef] [PubMed]
- Ferri, M.; Ranucci, E.; Romagnoli, P.; Giaccone, V. Antimicrobial resistance: A global emerging threat to public health systems. Crit. Rev. Food Sci. Nutr. 2017, 57, 2857–2876. [Google Scholar] [CrossRef] [PubMed]
- Nelson, K.; Hemarajata, P.; Sun, D.; Rubio-Aparicio, D.; Tsivkovski, R.; Yang, S.; Sebra, R.; Kasarskis, A.; Nguyen, H.; Hanson, B.M.; et al. Resistance to Ceftazidime-Avibactam is Due to Transposition of KPC in a Porin-Deficient Strain of Klebsiella pneumoniae with Increased Efflux Activity. Antimicrob. Agents Chemother. 2017, 61, e00989-17. [Google Scholar] [CrossRef] [PubMed]
- Wang, B.; Pan, F.; Wang, C.; Zhao, W.; Sun, Y.; Zhang, T.; Shi, Y.; Zhang, H. Molecular epidemiology of Carbapenem-resistant Klebsiella pneumoniae in a paediatric hospital in China. Int. J. Infect. Dis. 2020, 93, 311–319. [Google Scholar] [CrossRef]
- Galani, I.; Karaiskos, I.; Angelidis, E.; Papoutsaki, V.; Galani, L.; Souli, M.; Antoniadou, A.; Giamarellou, H. Emergence of ceftazidime-avibactam resistance through distinct genomic adaptations in KPC-2-producing Klebsiella pneumoniae of sequence type 39 during treatment. Eur. J. Clin. Microbiol. Infect. Dis. 2021, 40, 219–224. [Google Scholar] [CrossRef]
- Galani, I.; Nafplioti, K.; Adamou, P.; Karaiskos, I.; Giamarellou, H.; Souli, M. Study Collaborators Nationwide epidemiology of carbapenem resistant Klebsiella pneumoniae isolates from Greek hospitals, with regards to plazomicin and aminoglycoside resistance. BMC Infect. Dis. 2019, 19, 167. [Google Scholar]
- Karampatakis, T.; Tsergouli, K.; Politi, L.; Diamantopoulou, G.; Iosifidis, E.; Antachopoulos, C.; Karyoti, A.; Mouloudi, E.; Tsakris, A.; Roilides, E. Molecular Epidemiology of Endemic Carbapenem-Resistant Gram-Negative Bacteria in an Intensive Care Unit. Microb. Drug Resist. 2019, 25, 712–716. [Google Scholar] [CrossRef]
- Zarras, C.; Pappa, S.; Zarras, K.; Karampatakis, T.; Vagdatli, E.; Mouloudi, E.; Iosifidis, E.; Roilides, E.; Papa, A. Changes in molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae in the intensive care units of a Greek hospital, 2018–2021. Acta Microbiol. Immunol. Hung. 2022, 69, 104–108. [Google Scholar] [CrossRef] [PubMed]
- Karampatakis, T.; Antachopoulos, C.; Iosifidis, E.; Tsakris, A.; Roilides, E. Molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae in Greece. Future Microbiol. 2016, 11, 809–823. [Google Scholar] [CrossRef] [PubMed]
- Protonotariou, E.; Meletis, G.; Pilalas, D.; Mantzana, P.; Tychala, A.; Kotzamanidis, C.; Papadopoulou, D.; Papadopoulos, T.; Polemis, M.; Metallidis, S.; et al. Polyclonal Endemicity of Carbapenemase-Producing Klebsiella pneumoniae in ICUs of a Greek Tertiary Care Hospital. Antibiotics 2022, 11, 149. [Google Scholar] [CrossRef] [PubMed]
- Surveillance Atlas of Infectious Diseases. 2020. Available online: https://atlas.Ecdc.Europa.Eu/public/index.Aspx (accessed on 5 December 2022).
- Hatrongjit, R.; Kerdsin, A.; Akeda, Y.; Hamada, S. Detection of plasmid-mediated colistin-resistant and carbapenem-resistant genes by multiplex PCR. MethodsX 2018, 5, 532–536. [Google Scholar] [CrossRef]
- Mari-Almirall, M.; Ferrando, N.; Fernandez, M.J.; Cosgaya, C.; Vines, J.; Rubio, E.; Cuscó, A.; Muñoz, L.; Pellice, M.; Vergara, A.; et al. Clonal Spread and Intra- and Inter-Species Plasmid Dissemination Associated with Klebsiella pneumoniae Carbapenemase-Producing Enterobacterales during a Hospital Outbreak in Barcelona, Spain. Front. Microbiol. 2021, 12, 781127. [Google Scholar] [CrossRef]
- Yanat, B.; Rodriguez-Martinez, J.M.; Touati, A. Plasmid-mediated quinolone resistance in Enterobacteriaceae: A systematic review with a focus on Mediterranean countries. Eur. J. Clin. Microbiol. Infect. Dis. 2017, 36, 421–435. [Google Scholar] [CrossRef]
- Bosch, T.; Lutgens, S.P.M.; Hermans, M.H.A.; Wever, P.C.; Schneeberger, P.M.; Renders, N.H.M.; Leenders, A.C.A.P.; Kluytmans, J.A.J.W.; Schoffelen, A.; Notermans, D.; et al. Outbreak of NDM-1-Producing Klebsiella pneumoniae in a Dutch Hospital, with Interspecies Transfer of the Resistance Plasmid and Unexpected Occurrence in Unrelated Health Care Centers. J. Clin. Microbiol. 2017, 55, 2380–2390. [Google Scholar] [CrossRef]
- Silva, D.D.C.; Rampelotto, R.F.; Lorenzoni, V.V.; Santos, S.O.D.; Damer, J.; Horner, M.; Horner, R. Phenotypic methods for screening carbapenem-resistant Enterobacteriaceae and assessment of their antimicrobial susceptibility profile. Rev. Soc. Bras. Med. Trop. 2017, 50, 173–178. [Google Scholar] [CrossRef]
- Tzouvelekis, L.S.; Miriagou, V.; Kotsakis, S.D.; Spyridopoulou, K.; Athanasiou, E.; Karagouni, E.; Tzelepi, E.; Daikos, G.L. KPC-producing, multidrug-resistant Klebsiella pneumoniae sequence type 258 as a typical opportunistic pathogen. Antimicrob. Agents Chemother. 2013, 57, 5144–5146. [Google Scholar] [CrossRef]
- Karampatakis, T.; Zarras, C.; Pappa, S.; Vagdatli, E.; Iosifidis, E.; Roilides, E.; Papa, A. Emergence of ST39 carbapenem-resistant Klebsiella pneumoniae producing VIM-1 and KPC-2. Microb. Pathog. 2022, 162, 105373. [Google Scholar] [CrossRef]
- Protonotariou, E.; Meletis, G.; Chatzopoulou, F.; Malousi, A.; Chatzidimitriou, D.; Skoura, L. Emergence of Klebsiella pneumoniae ST11 co-producing NDM-1 and OXA-48 carbapenemases in Greece. J. Glob. Antimicrob. Resist. 2019, 19, 81–82. [Google Scholar] [CrossRef]
- Partridge, S.R.; Ginn, A.N.; Wiklendt, A.M.; Ellem, J.; Wong, J.S.; Ingram, P.; Guy, S.; Garner, S.; Iredell, J.R. Emergence of blaKPC carbapenemase genes in Australia. Int. J. Antimicrob. Agents 2015, 45, 130–136. [Google Scholar] [CrossRef] [PubMed]
- Tang, H.J.; Ku, Y.H.; Lee, M.F.; Chuang, Y.C.; Yu, W.L. In Vitro Activity of Imipenem and Colistin against a Carbapenem-Resistant Klebsiella pneumoniae Isolate Coproducing SHV-31, CMY-2, and DHA-1. BioMed Res. Int. 2015, 2015, 568079. [Google Scholar] [CrossRef] [PubMed]
- Kaczmarek, F.M.; Dib-Hajj, F.; Shang, W.; Gootz, T.D. High-level carbapenem resistance in a Klebsiella pneumoniae clinical isolate is due to the combination of bla(ACT-1) beta-lactamase production, porin OmpK35/36 insertional inactivation, and down-regulation of the phosphate transport porin phoe. Antimicrob. Agents Chemother. 2006, 50, 3396–3406. [Google Scholar] [CrossRef] [PubMed]
- Paczosa, M.K.; Mecsas, J. Klebsiella pneumoniae: Going on the Offense with a Strong Defense. Microbiol. Mol. Biol. Rev. 2016, 80, 629–661. [Google Scholar] [CrossRef]
- Brhelova, E.; Antonova, M.; Pardy, F.; Kocmanova, I.; Mayer, J.; Racil, Z.; Lengerova, M. Investigation of next-generation sequencing data of Klebsiella pneumoniae using web-based tools. J. Med. Microbiol. 2017, 66, 1673–1683. [Google Scholar] [CrossRef]
- Enany, S.; Zakeer, S.; Diab, A.A.; Bakry, U.; Sayed, A.A. Whole genome sequencing of Klebsiella pneumoniae clinical isolates sequence type 627 isolated from Egyptian patients. PLoS ONE 2022, 17, e0265884. [Google Scholar] [CrossRef]
- CLSI standard M02; Performance Standards for Antimicrobial Susceptibility Testing; 32nd Ed. Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2020.
- U.S. FDA. FDA approves new antibacterial drug Avycaz. In FDA News Release; U.S. FDA: Silver Spring, MD, USA, 2015. [Google Scholar]
- Tsakris, A.; Pournaras, S.; Woodford, N.; Palepou, M.F.; Babini, G.S.; Douboyas, J.; Livermore, D.D. Outbreak of infections caused by Pseudomonas aeruginosa producing VIM-1 carbapenemase in Greece. J. Clin. Microbiol. 2000, 38, 1290–1292. [Google Scholar] [CrossRef]
- Geneious Prime 2021.2.1. Available online: https://www.geneious.com (accessed on 1 December 2022).
- Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 2011, 12, 402. [Google Scholar] [CrossRef]
- Wu, K.M.; Li, L.H.; Yan, J.J.; Tsao, N.; Liao, T.L.; Tsai, H.C.; Fung, C.P.; Chen, H.J.; Liy, Y.M.; Chen, Y.T.; et al. Genome sequencing and comparative analysis of Klebsiella pneumoniae NTUH-K2044, a strain causing liver abscess and meningitis. J. Bacteriol. 2009, 191, 4492–4501. [Google Scholar] [CrossRef]
- Larsen, M.V.; Cosentino, S.; Rasmussen, S.; Friis, C.; Hasman, H.; Marvig, R.L.; Jelsbak, L.; Sicheritz-Ponten, T.; Ussery, D.W.; Aarestrup, F.M.; et al. Multilocus sequence typing of total-genome-sequenced bacteria. J. Clin. Microbiol. 2012, 50, 1355–1361. [Google Scholar] [CrossRef] [PubMed]
- Carattoli, A.; Zankari, E.; Garcia-Fernandez, A.; Voldby Larsen, M.; Lund, O.; Villa, L.; Møller Aarestrup, F.; Hasman, H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 2014, 58, 3895–3903. [Google Scholar] [CrossRef] [PubMed]
- Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
- Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017, Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef] [PubMed]
- Bertels, F.; Silander, O.K.; Pachkov, M.; Rainey, P.B.; van Nimwegen, E. Automated reconstruction of whole-genome phylogenies from short-sequence reads. Mol. Biol. Evol. 2014, 31, 1077–1088. [Google Scholar] [CrossRef] [PubMed]
- Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed]
- Dixit, P.D.; Pang, T.Y.; Studier, F.W.; Maslov, S. Recombinant transfer in the basic genome of Escherichia coli. Proc. Natl. Acad. Sci. USA 2015, 112, 9070–9075. [Google Scholar] [CrossRef]
- Sakoparnig, T.; Field, C.; van Nimwegen, E. Whole genome phylogenies reflect the distributions of recombination rates for many bacterial species. elife 2021, 10, e65366. [Google Scholar] [CrossRef]
- Zhang, R.; Liu, L.; Zhou, H.; Chan, E.W.; Li, J.; Fang, Y.; Li, Y.; Liao, K.; Chen, S. Nationwide Surveillance of Clinical Carbapenem-resistant Enterobacteriaceae (CRE) Strains in China. eBioMedicine 2017, 19, 98–106. [Google Scholar] [CrossRef]
- Voulgari, E.; Gartzonika, C.; Vrioni, G.; Politi, L.; Priavali, E.; Levidiotou-Stefanou, S.; Tsakris, A. The Balkan region: NDM-1-producing Klebsiella pneumoniae ST11 clonal strain causing outbreaks in Greece. J. Antimicrob. Chemother. 2014, 69, 2091–2097. [Google Scholar] [CrossRef]
- Todorova, B.; Sabtcheva, S.; Ivanov, I.N.; Lesseva, M.; Chalashkanov, T.; Ioneva, M.; Bachvarova, A.; Dobreva, E.; Kantardjiev, T. First clinical cases of NDM-1-producing Klebsiella pneumoniae from two hospitals in Bulgaria. J. Infect. Chemother. 2016, 22, 837–840. [Google Scholar] [CrossRef] [PubMed]
- Berglund, B.; Hoang, N.T.B.; Lundberg, L.; Le, N.K.; Tarnberg, M.; Nilsson, M.; Bornefall, E.; Khu, D.T.K.; Welander, J.; Le, H.T.; et al. Clonal spread of carbapenem-resistant Klebsiella pneumoniae among patients at admission and discharge at a Vietnamese neonatal intensive care unit. Antimicrob. Resist. Infect. Control 2021, 10, 162. [Google Scholar] [CrossRef] [PubMed]
- Liu, H.; Wilksch, J.; Li, B.; Du, J.; Cao, J.; Zhang, X.; Zhou, T. Emergence of ST39 and ST656 extensively drug-resistant Klebsiella pneumoniae isolates in Wenzhou, China. Indian. J. Med. Microbiol. 2017, 35, 145–146. [Google Scholar] [CrossRef]
- Cabanel, N.; Rosinski-Chupin, I.; Chiarelli, A.; Botin, T.; Tato, M.; Canton, R.; Glaser, P. Evolution of VIM-1-Producing Klebsiella pneumoniae Isolates from a Hospital Outbreak Reveals the Genetic Bases of the Loss of the Urease-Positive Identification Character. mSystems 2021, 6, e0024421. [Google Scholar] [CrossRef] [PubMed]
- Kuzina, E.S.; Kislichkina, A.A.; Sizova, A.A.; Skryabin, Y.P.; Novikova, T.S.; Ershova, O.N.; Savin, I.A.; Khokhlova, O.E.; Bogun, A.G.; Fursova, N.K. High-Molecular-Weight Plasmids Carrying Carbapenemase Genes bla(NDM-1), bla(KPC-2), and bla(OXA-48) Coexisting in Clinical Klebsiella pneumoniae Strains of ST39. Microorganisms 2023, 11, 459. [Google Scholar] [CrossRef]
- Giakkoupi, P.; Papagiannitsis, C.C.; Miriagou, V.; Pappa, O.; Polemis, M.; Tryfinopoulou, K.; Tzouvelekis, L.S.; Vatopoulos, A.C. An update of the evolving epidemic of blaKPC-2-carrying Klebsiella pneumoniae in Greece (2009-10). J. Antimihcrob. Chemother. 2011, 66, 1510–1513. [Google Scholar] [CrossRef]
- Piekarska, K.; Zacharczuk, K.; Wolkowicz, T.; Wolaniuk, N.; Rzeczkowska, M.; Gierczynski, R. Emergence of Enterobacteriaceae co-producing CTX-M-15, ArmA and PMQR in Poland. Adv. Clin. Exp. Med. 2019, 28, 249–254. [Google Scholar] [CrossRef]
- Mansour, W.; Grami, R.; Ben Haj Khalifa, A.; Dahmen, S.; Chatre, P.; Haenni, M.; Aouni, M.; Madec, J.Y. Dissemination of multidrug-resistant blaCTX-M-15/IncFIIk plasmids in Klebsiella pneumoniae isolates from hospital- and community-acquired human infections in Tunisia. Diagn. Microbiol. Infect. Dis. 2015, 83, 298–304. [Google Scholar] [CrossRef]
- Galani, I.; Karaiskos, I.; Souli, M.; Papoutsaki, V.; Galani, L.; Gkoufa, A.; Antoniadou, A.; Giamarellou, H. Outbreak of KPC-2-producing Klebsiella pneumoniae endowed with ceftazidime-avibactam resistance mediated through a VEB-1-mutant (VEB-25), Greece, September to October 2019. Eurosurveillance 2020, 25, 2000028. [Google Scholar] [CrossRef]
- Voulgari, E.; Kotsakis, S.D.; Giannopoulou, P.; Perivolioti, E.; Tzouvelekis, L.S.; Miriagou, V. Detection in two hospitals of transferable ceftazidime-avibactam resistance in Klebsiella pneumoniae due to a novel VEB beta-lactamase variant with a Lys234Arg substitution, Greece, 2019. Eurosurveillance 2020, 25, 1900766. [Google Scholar] [CrossRef]
- Livermore, D.M.; Struelens, M.; Amorim, J.; Baquero, F.; Bille, J.; Canton, R.; Henning, S.; Gatermann, S.; Marchese, A.; Mittermayer, H.; et al. Multicentre evaluation of the VITEK 2 Advanced Expert System for interpretive reading of antimicrobial resistance tests. J. Antimicrob. Chemother. 2002, 49, 289–300. [Google Scholar] [CrossRef] [PubMed]
- Garcia-Fulgueiras, V.; Magallanes, C.; Reyes, V.; Cayota, C.; Galiana, A.; Vieytes, M.; Vignoli, R.; Márquez, C. In Vivo High Plasticity of Multi-Drug Resistant ST258 Klebsiella pneumoniae. Microb. Drug Resist. 2021, 27, 1126–1130. [Google Scholar] [CrossRef] [PubMed]
- Albarri, O.; AlMatar, M.; Ocal, M.M.; Koksal, F. Overexpression of Efflux Pumps AcrAB and OqxAB Contributes to Ciprofloxacin Resistance in Clinical Isolates of K. pneumonia. Curr. Protein Pept. Sci. 2022, 23, 356–368. [Google Scholar] [CrossRef] [PubMed]
- Diez-Aguilar, M.; Canton, R. New microbiological aspects of fosfomycin. Rev. Esp. Quimioter. 2019, 32 (Suppl. S1), 8–18. [Google Scholar]
- Sugita, K.; Aoki, K.; Komori, K.; Nagasawa, T.; Ishii, Y.; Iwata, S.; Tateda, K. Molecular Analysis of bla(KPC-2)-Harboring Plasmids: Tn4401a Interplasmid Transposition and Tn4401a-Carrying ColRNAI Plasmid Mobilization from Klebsiella pneumoniae to Citrobacter europaeus and Morganella morganii in a Single Patient. mSphere 2021, 6, e0085021. [Google Scholar] [CrossRef]
- Di Martino, P.; Cafferini, N.; Joly, B.; Darfeuille-Michaud, A. Klebsiella pneumoniae type 3 pili facilitate adherence and biofilm formation on abiotic surfaces. Res. Microbiol. 2003, 154, 9–16. [Google Scholar] [CrossRef]
- Lam, M.M.C.; Wick, R.R.; Wyres, K.L.; Gorrie, C.L.; Judd, L.M.; Jenney, A.W.J.; Brisse, S.; Holt, K.E. Genetic diversity, mobilisation and spread of the yersiniabactin-encoding mobile element ICEKp in Klebsiella pneumoniae populations. Microb. Genom. 2018, 4, e000196. [Google Scholar] [CrossRef]
- Zhu, J.; Wang, T.; Chen, L.; Du, H. Virulence Factors in Hypervirulent Klebsiella pneumoniae. Front. Microbiol. 2021, 12, 642484. [Google Scholar] [CrossRef]
- Di Pilato, V.; Henrici De Angelis, L.; Aiezza, N.; Baccani, I.; Niccolai, C.; Parisio, E.M.; Giordano, C.; Camarlinghi, G.; Barnini, S.; Forni, S.; et al. Resistome and virulome accretion in an NDM-1-producing ST147 sublineage of Klebsiella pneumoniae associated with an outbreak in Tuscany, Italy: A genotypic and phenotypic characterisation. Lancet Microbe 2022, 3, e224–e234. [Google Scholar] [CrossRef]
- Miro, E.; Rossen, J.W.A.; Chlebowicz, M.A.; Harmsen, D.; Brisse, S.; Passet, V.; Navarro, F.; Friedrich, A.W.; García-Cobos, S. Core/Whole Genome Multilocus Sequence Typing and Core Genome SNP-Based Typing of OXA-48-Producing Klebsiella pneumoniae Clinical Isolates From Spain. Front. Microbiol. 2019, 10, 2961. [Google Scholar] [CrossRef]
- Snitkin, E.S.; Zelazny, A.M.; Thomas, P.J.; Stock, F.; NISC Comparative Sequencing Program Group; Henderson, D.K.; Palmore, T.N.; Segre, J.A. Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing. Sci. Transl. Med. 2012, 4, 148ra116. [Google Scholar] [CrossRef] [PubMed]
- Marsh, J.W.; Krauland, M.G.; Nelson, J.S.; Schlackman, J.L.; Brooks, A.M.; Pasculle, A.W.; Shutt, K.A.; Doi, Y.; Querry, A.M.; Muto, C.A.; et al. Genomic Epidemiology of an Endoscope-Associated Outbreak of Klebsiella pneumoniae Carbapenemase (KPC)-Producing K. pneumoniae. PLoS ONE 2015, 10, e0144310. [Google Scholar] [CrossRef] [PubMed]
- Onori, R.; Gaiarsa, S.; Comandatore, F.; Pongolini, S.; Brisse, S.; Colombo, A.; Cassani, G.; Marone, P.; Grossi, P.; Minoja, G.; et al. Tracking Nosocomial Klebsiella pneumoniae Infections and Outbreaks by Whole-Genome Analysis: Small-Scale Italian Scenario within a Single Hospital. J. Clin. Microbiol. 2015, 53, 2861–2868. [Google Scholar] [CrossRef]
- Arnold, B.J.; Gutmann, M.U.; Grad, Y.H.; Sheppard, S.K.; Corander, J.; Lipsitch, M.; Hanage, W.P. Weak Epistasis May Drive Adaptation in Recombining Bacteria. Genetics 2018, 208, 1247–1260. [Google Scholar] [CrossRef] [PubMed]
- Dallman, T.; Ashton, P.; Schafer, U.; Jironkin, A.; Painset, A.; Shaaban, S.; Hartman, H.; Myers, R.; Underwood, A.; Jenkins, C.; et al. SnapperDB: A database solution for routine sequencing analysis of bacterial isolates. Bioinformatics 2018, 34, 3028–3029. [Google Scholar] [CrossRef]
Figure 1.
MICs (μg/mL) heatmap of the 24 clinical isolates of carbapenem-resistant K. pneumoniae (CRKP) according to CLSI.
Figure 1.
MICs (μg/mL) heatmap of the 24 clinical isolates of carbapenem-resistant K. pneumoniae (CRKP) according to CLSI.
Figure 2.
Circular genome map of (A) seven ST11 CRKP strains, (B) the seven ST39 CRKP strains, (C) the three ST323 CRKP strains, and (D) the three ST15 CRKP strains. The genome maps were produced through the BLAST Ring Image Generator (BRIG) software. The genomes of the seven ST11 CRKP isolates (A10143, A7213, A9974, C4112, C746, C833, and D2856), the seven ST39 CRKP isolates (B11395, C2482, D1598, D2452, D6184, Z508 and C1909), the three ST323 CRKP strains (Z557, D1463, and Z866) and the three ST15 CRKP strains (A10037, A5051, and B2562) were mapped separately to the reference genome K. pneumoniae NC_012731. From the inner to the external ring, the innermost ring displays the size of the genome in kbp, followed by GC content (black), GC skew (dark green and dark purple), and K. pneumoniae NC_012731 (red). The white vertical gaps represent sequences of the K. pneumoniae NC_012731 reference strain that are absent in the sequences of the ST11, ST39, ST323, and ST15 CRKP strains or present sequences <50% identical to the reference genome.
Figure 2.
Circular genome map of (A) seven ST11 CRKP strains, (B) the seven ST39 CRKP strains, (C) the three ST323 CRKP strains, and (D) the three ST15 CRKP strains. The genome maps were produced through the BLAST Ring Image Generator (BRIG) software. The genomes of the seven ST11 CRKP isolates (A10143, A7213, A9974, C4112, C746, C833, and D2856), the seven ST39 CRKP isolates (B11395, C2482, D1598, D2452, D6184, Z508 and C1909), the three ST323 CRKP strains (Z557, D1463, and Z866) and the three ST15 CRKP strains (A10037, A5051, and B2562) were mapped separately to the reference genome K. pneumoniae NC_012731. From the inner to the external ring, the innermost ring displays the size of the genome in kbp, followed by GC content (black), GC skew (dark green and dark purple), and K. pneumoniae NC_012731 (red). The white vertical gaps represent sequences of the K. pneumoniae NC_012731 reference strain that are absent in the sequences of the ST11, ST39, ST323, and ST15 CRKP strains or present sequences <50% identical to the reference genome.
Figure 3.
The tree built on the entire core genome with PhyML. The branch color indicates the fraction of supporting versus supporting plus clashing SNPs for each branch of the core tree, which supports the initial bipartition of the strains. MLST types: ST11: D2856, C4112, C833, C746, A10143, A9974, A7213 (blue color), ST512: A18940 (brown color), ST30: C251 (brown color), ST307: Z852 (pink color), ST323: Z866, D1463, Z557 (green color), ST35: A1746 (orange color), ST15: B2562, A5051, A10037 (purple color), ST39: D1598, D2452, C2482, Z508, C1909, B11395, D6184 (red color).
Figure 3.
The tree built on the entire core genome with PhyML. The branch color indicates the fraction of supporting versus supporting plus clashing SNPs for each branch of the core tree, which supports the initial bipartition of the strains. MLST types: ST11: D2856, C4112, C833, C746, A10143, A9974, A7213 (blue color), ST512: A18940 (brown color), ST30: C251 (brown color), ST307: Z852 (pink color), ST323: Z866, D1463, Z557 (green color), ST35: A1746 (orange color), ST15: B2562, A5051, A10037 (purple color), ST39: D1598, D2452, C2482, Z508, C1909, B11395, D6184 (red color).
Figure 4.
Pairwise analysis of recombination between carbapenem-resistant K. pneumoniae CRKP strains. (A–C) SNP densities (SNPs per kilobase) along the core genome for three pairs of CRKP strains with various nucleotide diversities. (D–F) Corresponding histograms for the number of SNPs per kilobase (dots) together with fits of the Poisson mixture model for C4112-C746 (green), A18940-D2856 (red), and A5051-C2482 (blue) CRKP strains. The vertical axis is on a logarithmic scale. (G) For each pair of CRKP strains (dots), the fraction of the genome that was inherited clonally is shown as a function of the nucleotide divergence of the pair, shown on a logarithmic scale. The three CRKP pairs that are highlighted in panels (A–F) are depicted as green, red, and blue dots. (H) Fraction of all SNPs that lie in recombined regions as a function of the clonally inherited fraction of the genome. (I) Cumulative distribution of the clonal fractions of CRKP pairs. (J) Cumulative distribution of the lengths of recombined segments for CRKP pairs that are in the mostly clonal regime.
Figure 4.
Pairwise analysis of recombination between carbapenem-resistant K. pneumoniae CRKP strains. (A–C) SNP densities (SNPs per kilobase) along the core genome for three pairs of CRKP strains with various nucleotide diversities. (D–F) Corresponding histograms for the number of SNPs per kilobase (dots) together with fits of the Poisson mixture model for C4112-C746 (green), A18940-D2856 (red), and A5051-C2482 (blue) CRKP strains. The vertical axis is on a logarithmic scale. (G) For each pair of CRKP strains (dots), the fraction of the genome that was inherited clonally is shown as a function of the nucleotide divergence of the pair, shown on a logarithmic scale. The three CRKP pairs that are highlighted in panels (A–F) are depicted as green, red, and blue dots. (H) Fraction of all SNPs that lie in recombined regions as a function of the clonally inherited fraction of the genome. (I) Cumulative distribution of the clonal fractions of CRKP pairs. (J) Cumulative distribution of the lengths of recombined segments for CRKP pairs that are in the mostly clonal regime.
Figure 5.
SNP compatibility along the core genome alignment displays the shortness of tree-compatible segments. (A) Probability distribution of the number of consecutive SNP columns that are consistent with a common phylogeny for the core genome alignment. (B) Probability distribution of the number of consecutive nucleotides consistent with a common phylogeny for the core genome alignment.
Figure 5.
SNP compatibility along the core genome alignment displays the shortness of tree-compatible segments. (A) Probability distribution of the number of consecutive SNP columns that are consistent with a common phylogeny for the core genome alignment. (B) Probability distribution of the number of consecutive nucleotides consistent with a common phylogeny for the core genome alignment.
Figure 6.
Ratio C/S of the minimal number of phylogeny changes C to substitutions S for random subsets of carbapenem-resistant K. pneumoniae (CRKP) strains using the alignment from which 5% of potentially homoplastic positions have been removed. For CRKP strain numbers ranging from n = 4 to n = 24, random subsets of n strains were collected, and the ratios C/S of phylogeny changes to SNPs in the alignment were calculated. The figure shows box-whisker plots that indicate, for each CRKP strain number n, means, and standard deviations for C/S across subsets. The horizontal red line displays C/S = 0.15.
Figure 6.
Ratio C/S of the minimal number of phylogeny changes C to substitutions S for random subsets of carbapenem-resistant K. pneumoniae (CRKP) strains using the alignment from which 5% of potentially homoplastic positions have been removed. For CRKP strain numbers ranging from n = 4 to n = 24, random subsets of n strains were collected, and the ratios C/S of phylogeny changes to SNPs in the alignment were calculated. The figure shows box-whisker plots that indicate, for each CRKP strain number n, means, and standard deviations for C/S across subsets. The horizontal red line displays C/S = 0.15.
Table 1.
Genetic characteristics of the 24 carbapenem-resistant K. pneumoniae isolates of the study.
Table 1.
Genetic characteristics of the 24 carbapenem-resistant K. pneumoniae isolates of the study.
| Isolates | ST | Plasmids | Antibiotic Resistance Genes Affecting | Virulence | Efflux/Regulator | Capsule | ||||
|---|---|---|---|---|---|---|---|---|---|---|
| Other Beta-Lactams | Carba-Penems | Amino-Glycosides | Quinolones | Other | ||||||
| C2482/18 | 39 | ColRNAI IncA/C2 IncFIB IncFIB(pQil) Col (pHAD28) | blaSHV-79 blaTEM-1B | blaKPC-2 blaVIM-1 | aac(3)-Iid aac(6′)-Ib aac(6′)-Im aant(2″)-Ia aph(2″)-Ib aph(3′)-Ia | oqxA, oqxB | fosA, sul1, Sul3, dfrA1 | fyuA, mrkABCDFHIJ ybtASTX | acrB, marAR, soxSR, ramA, rob sdiA, fis, envR, oqxBR, rarA | wzc, |
| C4112/18 | 11 | ColRNAI IncFIA(HI1) IncFIB(K) IncFII(K) | blaCTX-M-15 blaSHV-11 | blaNDM-1 | aac(3)-Iia, aac(6′)-Ib aph(3′)-Ia | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2 dfrA14, catB3, tet(D) | iutA, mrkABCDFHIJ | acrABR marAR soxSR ramAR rob sdiA, fis, envR, oqxR, rarA | wzc, wzi |
| B11395/18 | 39 | IncFIB (AP001918) IncFIB(pQil) IncFII(K) | blaSHV-79 blaTEM-1B | blaKPC-2 | aac(3)-Iid, aadA1, aadA2, aph(3′)-Ia | oqxA, oqxB | fosA, sul2, sul3, dfrA12, cmlA1 | fyuA, mrkABCDFHIJ, ybtEQSTU | acrABR marAR soxSR, ramA, sdiA, fis, envR, oqxA, rarA | wzc, wzi |
| D6184/18 | 39 | ColRNAI IncFIB (AP001918) IncFIB(pQil) | blaSHV-11 blaSHV-40 blaΤΕΜ-1Β | blaKPC-2 | aac(3)-Iid, aac(6′)-Ib, aph(3′)-Ia, aadA1, aadA2 | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2, sul3, tet(M), dfrA12, catA1, atB3, | fyuA, irp1, irp2, mrkABCDFHIJ, ybtAEPQSTUX | acrABR marAR oxSR, ramA, rob, sdiA, fis, envR, qxABR, rarA | wzc, wzi |
| C251/19 | 30 | ColRNAI IncA/C2 IncFIB(K) IncFII(K) IncX3 | blaSHV-12, blaTEM-1A blaOXA-9 | blaKPC-2 blaVIM-1 | aac(6′)-Ib aph(3′)-Ia aadA2 aadA15 | qnrS1, aac(6′)-lb-cr, oqxA, oqxB | fosA, sul1, dfrA1, dfrA12, catA, tet(D) | fyuA, irp1, iutA mrkBCDFIJ, ybtAEPQSTU | acrABR, marAR, soxSR, ramAR, rob, sdiA, envR, oqxABR, rarA | wzc, wzi |
| Ζ557/19 | 323 | IncFIB(K) IncFIB(pQil) IncFII(K) IncFIB(Mar) IncFIB (pKPHS1) IncHI1B | blaSHV-1 blaSHV-99 | blaKPC-2 blaVIM-1 | ant(2″)-Ia, aac(6′)-Il, aadA1 | qnrA1 oqxA, oqxB | fosA, fosA7, sul1, dfrA1 | mrkABCDFIJ | acrAR, marA, soxSR, ramA, rob, sdiA, fis, envR, oqxA, rarA | wzc, wzi |
| C746/19 | 11 | ColRNAI IncFIB(K) IncFIB(pQil) IncFII(K), IncR | blaSHV-1 blaOXA-1 | blaKPC-2 | aac(6′)-lb-cr aph(3′)-Ia aadA2 | aac(6′)-lb-cr, oqxA oqxB | fosA, sul1 | fyuA, irp1, irp2 iutA, mrkABCDFIJ, ybtAEUX | acrAB, marAR, soxSR, ramAR, rob, sdiA, fis, envR, oqxABR, | wzi |
| B2562/19 | 15 | IncFIB(K) IncFII(K) IncFIA(HI1) | blaCTX-M-15 blaTEM-1B blaOXA-1 | blaNDM-1 | aac(3)-Iia aac(6′)-Ib3 aph(3″)-Ib aph(6)-Id | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2, dfrA14, catB3 | fyuA, irp1, kfuABC, mrkABCDFIJ, ybtAPQSTUX | acrABR, marAR, soxSR, ramA, rob, sdiA, fis, envR, rarA | wzi |
| Ζ508/19 | 39 | ColRNAI IncA/C2 IncFIB IncFIB(pQil) | blaSHV-11 blaSHV-40 blaΤΕΜ-1Β blaOXA-116 blaOXA-366 | blaKPC-2 blaVIM-1 | aac(3)-Iid aac(6′)-Ib3 aac(6′)-Im aph(2″)-Ib aph(3′)-Ia aph(3′)-Via aadA24, armA, aac(6’)-Il | qnrS1 aac(6′)-lb-cr, oqxA, oqxB | fosA, sul1, sul3, frA1, tet(M) | fyuA, irp1, irp2, mrkABCDFHIJ, ybtAEPQSTUX | acrBR, marAR, soxSR, ramA, rob, sdiA, fis, envR, oqxABR, rarA | wzc, wzi |
| D1463/19 | 323 | IncFIB(K) IncFIB(Mar) IncFIB (pKPHS1) IncFIB(pQil) IncFII(K) IncHI1B | blaSHV-99 | blaKPC-2 blaVIM-1 | ant(2″)-Ia aac(6′)-Il | oqxA oqxB | fosA, fosA7, dfrA1 | mrkABCDFHJ, | acrAR, marA, soxSR, ramA, rob, sdiA, fis, envR, oqxA, rarA | wzc, wzi |
| D1598/19 | 39 | ColRNAI IncA/C2 IncFIB IncFIB(pQil) | blaSHV-40 blaTEM-1B | blaKPC-2 blaVIM-1 | aac(3)-Iid aac(6′)-Il aac(6′)-Im ant(3)-Ia aph(2″)-Ib aph(3′)-Ia | qnrS1, aac(6′)-lb-cr, oqxA oqxB | fosA, sul1, sul3, dfrA1 | ND | ND | wzi |
| Z852/19 | 307 | IncA/C2 IncFIB(pQil) | blaSHV-28 | blaKPC-2 blaVIM-1 | aac(6′)-Il aadA24 | oqxA oqxB | fosA, sul1, dfrA1, frA14 | fyuA, irp1, irp2 iutA, mrkABCDFHIJ, ybtAEPQSTX | acrABR, marAR, soxSR, ramAR, rob, sdiA, fis, envR, oqxABR, rarA | wzc, wzi |
| Z866/19 | 323 | IncFIB(K) IncFIB(Mar) IncFIB (pKPHS1) IncFIB(pQil) IncFII(K) IncHI1B | blaSHV-99 | blaKPC-2 blaVIM-1 | aadA1 | oqxA oqxB | fosA, fosA7, Sul1 | mrkABCDFHIJ | acrAR, marA, soxSR, ramA, rob, sdiA, fis, envR, oqxAR, rarA | wzc, wzi |
| D2452/19 | 39 | ColRNAI IncA/C2 IncFIB(pQil) | blaSHV-79 blaTEM-1B | blaKPC-2 blaVIM-1 | aac(3)-Iid aac(6′)-Il aac(6′)-Im ant(3)-Ia aph(2″)-Ib aph(3′)-Ia | qnrS1, aac(6′)-lb-cr, oqxA oqxB | fosA, sul1, Sul4, dfrA1 | fyuA, irp1, irp2 mrkABCDFHIJ, ybtAEPQSTX | acrBR, marAR, soxSR, ramA, rob, sdiA, fis, envR, oqxABR, rarA | wzc, wzi |
| C833/21 | 11 | Col(pHAD2) IncFIA(HI1) IncFIB(K) IncFII(K) | blaSHV-182 | blaOXA-48 blaNDM-1 | aph(3′)-Ia aac(6′)-lb-cr | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2, dfrA14, catB3 | mrkABCDFHIJ | acrABR, marA, soxSR, ramAR, rob, sdiA, fis, envR, rarA | wzc, wzi |
| A5051/21 | 15 | IncFIA(HI1) IncFIB(K) IncFII(K) | blaSHV-28 blaSHV-106 | blaNDM-1 | aac(6′)-lb-cr | aac(6′)-lb-cr oqxA | fosA catB3 | fyuA, kfuC, mrkABCDFHIJ ybtQU | acrABR, marA, soxSR, ramAR, rob, sdiA, fis, envR, rarA | wzcwzi |
| D2856/21 | 11 | Col440II IncFIA(HI1) IncFIB(K) IncFII(K) IncR | blaSHV-182 | blaOXA-48 blaNDM-1 | aph(3′)-Ia aac(6′)-lb-cr | oqxA | fosA, sul2, dfrA14, catB3, tet(A) | mrkABCDJ, ybtAPQS | acrABR, marAR, soxSR, ramAR, rob, sdiA, fis, envR, oqxAR, rarA | wzcwzi |
| C1909/21 | 39 | ColRNAI IncFII(K) | blaSHV-40, -79, -85, -89 | blaKPC-33 | aac(6′)-Ib aac(6′)-lb-cr | aac(6′)-lb-cr, oqxA, oqxB | fosA, dfrA12 tet(A) | mrkABCDFHIJ ybtEPQSTU | acrA, marAR, soxSR, ramAR, rob, sdiA, fis, envR, oqxAR, rarA | wzcwzi |
| A18940/21 | 512 | IncFIB(K) IncFII(K) IncN IncX3 | blaSHV-182 blaOXA-10 | blaKPC-2, blaNDM-1 | aac(6′)-Ib-cr, aac(6′)-Ib, aph(3′)-Ia, aadA2, aadA16 | qnrS1, aac(6′)-lb-cr, oqxA, oqxB | fosA, sul1, dfrA12, dfrA27 catA1 | iutA, mrkABCHI | acrR, marAR soxRS, ramA, rob, sdiA, fis, envR, oqxR, rarA | wzi |
| A1746/22 | 35 | IncC, IncR, IncFIA(HI1) IncFIB(K) IncFIB (pKPHS1) IncFIB(pQil) IncFII(K) | blaSHV-33 blaOXA-10 blaTEM-1B blaVEB-25 blaDHA-1 | blaKPC-2 | ant(2″)-Ia, aph(3″)-Ib aph(6)-Id, rmtB,aadA1 | qnrB4 oqxA, oqxB | fosA, sul1 sul2, catA1, cmlA1, tet(A), tet(G) | kfuA, mrkAFHI, ybtEQTX | acrR, marAR, soxRS, ramA rob, sdiA, fis, envR, oqxR rarA | wzc, wzi |
| A7213/22 | 11 | ColRNAI IncC IncFIA(HI1) IncFIB(K) | blaOXA-1 blaOXA-10 blaCTX-M -15 blaTEM-1B blaVEB-1 | blaNDM-1 | aac(6′)-Ib, aac(3)-IIa, aac(6′)-Ib-cr, ant(2″)-Ia, aph(3′)-Ia, rmtB | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2, dfrA14, tet(A) | mrkAC | soxSR, ramAR, sdiA | wzi |
| A9974/22 | 11 | ColRNAI IncC IncFIA(HI1) IncFIB(K) IncFII(K) | blaSHV-159 blaOXA-1 blaOXA-10 blaCTX-M-15 blaTEM-1B blaVEB-1 | blaNDM-1 | aac(6′)-Ib-cr aac(6′)-Ib aac(3)-IIa aph(6)-Id aph(3′)-Ia aph(3″)-Ib ant(2″)-Ia aadA1, rmtB | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2, dfrA14, tet(A), tet(G) | fyuA, mrkABCDFHIJ, ybtAEPQST | acrA, marAR soxR, ramAR rob, sdiA, fis, envR, oqxAB, rarA | wzi |
| A10037/22 | 15 | IncC IncFIA(HI1) IncFIB(K) IncFII(K) | blaSHV-100 blaOXA-1 blaOXA-10 blaCTX-M-15 blaTEM-1B blaVEB-1 | blaNDM-1 | aac(6′)-Ib-cr aac(6′)-Ib aac(3)-IIa aph(3″)-Ib aph(6)-Id ant(2″)-Ia aadA1, rmtB | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul2 dfrA14, cmlA1, tet(A), tet(G) | fyuA, kfuAC, mrkABCFHJ, ybtAEPQSTU | acrAR, marAR soxSR, rob, sdiA fis, envR, oqxAR rarA | wzi |
| A10143/22 | 11 | ColRNAI IncFII(K) | blaOXA-1 blaOXA-10 blaCTX-M-15 blaTEM-1B blaVEB-1 | blaNDM-1 | aac(6′)-Ib-cr aac(6′)-Ib aac(3)-IIa aph(3″)-Ib aph(6)-Id ant(2″)-Ia aadA1, rmtB | aac(6′)-lb-cr, oqxA, oqxB | fosA, sul1, sul2, dfrA12, dfrA14, cmlA1, tet(A), tet(G) | fyuA, mrkH, ybtEPQSX | marA, soxSR, ramR, sdiA, fis, envR | wzi |
Yersiniabactin cluster: ybt, irp, and fyuA genes; Aerobactin cluster: iuc and iut genes; AcrAB efflux pump: acr, mar, sox, rob, ram sdiA, fis, and envR, genes; OqxAB efflux pump: oqx and rar genes.
Table 2.
Clinical data of carbapenem-resistant Klebsiella pneumoniae isolates.
| Strain | Date of Isolation | Sex | Age (Years) | Ward | Biological Sample | Infection/Colonization |
|---|---|---|---|---|---|---|
| C2482/18 | 9 June 2018 | F | 8 | PICU | Rectal swab | Colonization |
| C4112/18 | 7 September 2018 | F | 56 | ICU | Rectal swab | Colonization |
| B11395/18 | 24 October 2018 | F | 1 | PICU | Urine | Infection |
| D6184/18 | 19 November 2018 | F | 8 | PICU | Bronchial aspirate | Infection |
| C251/19 | 17 January 2019 | F | 17 | PICU | Rectal swab | Colonization |
| Z557/19 | 30 January 2019 | M | 57 | ICU | Rectal swab | Colonization |
| C746/19 | 7 February 2019 | M | 73 | ICU | Rectal swab | Colonization |
| B2562/19 | 4 March 2019 | M | 79 | ICU | Urine | Infection |
| Z508/19 | 20 March 2019 | F | 60 | ICU | Rectal swab | Colonization |
| D1463/19 | 22 March 2019 | M | 55 | ICU | Bronchial aspirate | Infection |
| D1598/19 | 1 April 2019 | M | 40 | ICU | Bronchial aspirate | Infection |
| Z852/19 | 12 April 2019 | M | 3 | PICU | Rectal swab | Colonization |
| Z866/19 | 15 April 2019 | M | 49 | ICU | Rectal swab | Colonization |
| D2452/19 | 24 May 2019 | M | 73 | ICU | Bronchial aspirate | Infection |
| C833/21 | 31 March 2021 | F | 57 | ICU | Wound | Infection |
| A5051/21 | 25 April 2021 | F | 69 | ICU | Blood | Infection |
| D2856/21 | 7 July 2021 | F | 54 | ICU | Bronchial aspirate | Infection |
| C1909/21 | 8 July 2021 | M | 85 | ICU | Wound | Infection |
| A18940/21 | 24 December 2021 | F | 58 | ICU | Blood | Infection |
| A1746/22 | 25 February 2022 | M | 0.25 | PICU | Blood | Infection |
| A7213/22 | 3 May 2022 | F | 30 | ICU | Blood | Infection |
| A9974/22 | 16 June 2022 | F | 67 | ICU | Blood | Infection |
| A10037/22 | 16 June 2022 | F | 53 | ICU | Blood | Infection |
| A10143/22 | 19 June 2022 | M | 88 | ICU | Blood | Infection |
ICU: Intensive Care Unit, PICU: Pediatric Intensive Care Unit.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Zarras, C.; Karampatakis, T.; Pappa, S.; Iosifidis, E.; Vagdatli, E.; Roilides, E.; Papa, A. Genetic Characterization of Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates in a Tertiary Hospital in Greece, 2018–2022. Antibiotics 2023, 12, 976. https://doi.org/10.3390/antibiotics12060976
AMA Style
Zarras C, Karampatakis T, Pappa S, Iosifidis E, Vagdatli E, Roilides E, Papa A. Genetic Characterization of Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates in a Tertiary Hospital in Greece, 2018–2022. Antibiotics. 2023; 12(6):976. https://doi.org/10.3390/antibiotics12060976
Chicago/Turabian StyleZarras, Charalampos, Theodoros Karampatakis, Styliani Pappa, Elias Iosifidis, Eleni Vagdatli, Emmanuel Roilides, and Anna Papa. 2023. "Genetic Characterization of Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates in a Tertiary Hospital in Greece, 2018–2022" Antibiotics 12, no. 6: 976. https://doi.org/10.3390/antibiotics12060976
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
