1. Introduction
Klebsiella pneumoniae frequently cause community-acquired and nosocomial infections such as pneumonia, urinary tract infection, liver abscesses and bloodstream infections [
1]. In recent years, antimicrobial resistance in
K. pneumoniae has become problematic [
2,
3]. In particular, resistance to carbapenems is frequently associated with resistance to multiple classes of other antibiotics which leads to limited possibilities for clinical management as alternative treatment options are limited, and lead to higher rates of treatment failures [
4]. In such cases last-resort antibiotics such as colistin [
5,
6] may be used in association with other antibiotics. However, the high worldwide prevalence of carbapenem-resistant, multi-drug resistant (MDR)
K. pneumoniae has fueled and increased the use of colistin over the last years, accelerating the emergence of isolates resistant to this compound [
3]. Besides full colistin resistance (CR), colistin heteroresistance (CHR) has also been increasingly reported over the last years [
7,
8]. Heteroresistance (HR) is defined as a phenotype in which a bacterial isolate contains subpopulations of cells that show a substantial reduction in antibiotic susceptibility compared with the main population (minimum inhibitory concentration (MIC) increase of at least 8-fold), allowing these cells to grow in the presence of the antibiotic [
9]. Detection of these subpopulations can be challenging and raise concern as their frequency may rise during antibiotic exposure and possibly lead to treatment failure [
9,
10,
11]. However, the relevance of CHR in causing reduced clinical effectiveness and negatively affect the treatment outcome remains controversial [
7]. Recent studies have highlighted the high prevalence of CHR but these were often limited by small sample size (mostly relying on local data) and limited investigation of the characteristics of the CHR
K. pneumoniae strains. In this study, we took advantage of the unique availability of a large collection of clinical isolates from various intensive care units (ICUs) across Europe to gain more insight into the frequency of occurrence of CHR among colistin-sensitive (CS)
K. pneumoniae and to investigate its possible association with epidemiological and clinical characteristics such as country of origin, the different therapeutic intervention strategies to which the patients were exposed during their stay in the ICU and sequence type (ST type). We also investigated CHR
K. pneumoniae strains for their possible association with, capsular polysaccharide antigen types (K-antigen type), lipopolysaccharide antigen types (O-antigen type) and colistin MIC value.
4. Discussion
In general, research around HR has been proven to be challenging especially due to a lack standardization of identification methods as well as of precise classification criteria. In the past it has therefore also been difficult to come to reliable estimates on the prevalence of HR as the sample sizes of the performed studies have been generally small and, due to a lack of standardization, both in definitions of HR but also in assays, it has been difficult to compile and compare different studies [
1,
9]. There are various ways through which HR can be screened, including methods that are used routinely for example for MIC determination, but they have been deemed to be not performant enough [
9,
18,
19].
A strength of this study was its large sample size compared to previous studies on CHR using the PAP assay [
20,
21,
22,
23]. Within this study two different classification schemes for CHR were used. Classification 1 was considered to be more stringent due to the higher frequency requirement and the number of isolates fulfilling this more stringent classification CHR was thus limited (n=25). The majority of isolates categorized as CHR fulfilled the second less stringent classification (n=108) which was considered to be less stringent due to a lower frequency requirement though this classification included an additional requirement with regard to the concentration at which growth should occur. It has to be acknowledged that there were also numerous isolates that could not classified as CHR by any of these two classifications but did show growth on PAP assay plates >2 mg/L colistin. Most often, this was growth at a frequency <1x10
-7 (often only one colony). We cannot rule out that these colonies were spontaneous mutants and were thus not caused by the isolate being CHR or were due to an inoculum effect. They were therefore not classified as CHR.
This study was also unique in the associations investigated, which to our knowledge, have not been studied previously in
K. pneumoniae. A study on CHR in
A. baumannii reported the absence of association between CHR and MIC-value and clonal complexes [
22], however, these results are not necessarily applicable to CHR in
K. pneumoniae. In this study, we could not find any association between CHR and country, intervention strategy, K-antigen type or O-antigen type or colistin MIC value. For the ST type an association was found, however, only when classification 2 was used. In contrast, an association was found between CR and country, ST type, O-antigen type and K-type. It is interesting to note that, based on an annual report from the European Centre for Disease Prevention and Control (ECDC) on antimicrobial consumption, Spain, which was found to have a higher number of CR isolates, also had a relatively higher consumption of polymyxins in hospitals compared to other European countries at least in the last period of the trial in which these isolates were collected (2016-2017) [
24]. Also during 2017, colistin consumption in Spain remained relatively much higher compared to other countries included here [
24]. Unfortunately, no data was available for a larger part of the duration of the trial (2013-2016).
For the ST-, O-antigen and K-antigen type, we often saw a relationship between the different molecular indicators and the ST types found to have a statistically significant difference in the amount of CS and CR isolates e.g. a specific K-antigen type was only present in combination with a specific ST type and were both found to have a higher proportion of CR isolates. For CHR, interestingly, no such pattern was found. CHR ST types found to have an association with CR in this study are especially interesting since more than 90% of them are known to be associated with multidrug resistance [
25,
26,
27]. In this study, MDR rates were high in CS, CHR and CR isolates (>75%) and around one third of the isolates determined to be CHR were CP-Kpn. Especially CHR in MDR CP-Kpn is especially important to report as in most sites having a historically high proportion of CP-Kpn, colistin is empirically prescribed and knowledge on the predilection towards development of CHR might facilitate the use of tailored antibiotic combinations instead.
In contrast to CS and CHR, the proportion of CR CP-Kpn isolates was relatively high (±30% vs. ±70%). Since CP-Kpn are more likely to be exposed to colistin, it is to be expected that CR rates will increase over time in this subset of isolates which is also reflected in this study with the far majority of CP-Kpn being CR. Interestingly, rates of CHR among CP-Kpn were similar to CS among CP-Kpn, and CHR rates also did not differ significantly between the baseline, CHX, SOD and SDD. This also suggests that colistin exposure through SOD and SDD had no association with selection of isolates with a CHR phenotype.
Finally, a closer look at the mechanism of colistin resistance for the resistant subpopulation of two CHR isolates showed that mutations in colistin resistance-conferring genes assessed in this study could not always be identified for each isolate. A recent study on CHR in wild-type
K. pneumoniae isolates also reported that for 28% of mutants sequenced, no genetic modification was found in the panel of genes assessed [
28]. Disruption of
mgrB by insertion sequences was the most commonly found genetic modification in the resistant subpopulations. Additionally, a nonsense mutation and complete deletion of
mgrB was found once and twice, respectively. A mutation outside of
mgrB was only found once. The same study in wild-type
K. pneumoniae isolates also reported a high number of
mgrB genetic modifications (54%) [
28]. A complete deletion was only found in 4% of the mutants while disruption by insertion sequences and other amino acid alterations were found in 28% and 22% of mutants respectively [
28]. Luo et al. also reported in their study that there was a high rate of
mgrB insertional mutations and no mutations in
pmrAB or
phoPQ and stated that this was consistent with previous findings which showed that those genes had a significantly lower mutation rate compared to
mgrB. However, they also stated that the high amount of
mgrB disruptions may be related to the high prevalence of ST11 since this ST type showed a significantly higher rate of
mgrB disruptions compared to other ST types in their study [
29].
Though this study has helped in expanding the knowledge on CHR in a clinical setting, there are also some limitations. Firstly, only limited sequencing data was available for CHR isolates. Future studies are needed to further assess the diversity of mechanisms of CHR and whether these mechanisms are ST-specific. Secondly, we cannot exclude confounding factors such as the usage of more colistin in some local settings, outbreaks with specific STs, and the prevalence of CP Kpn which may vary between sites. Thirdly, we only studied the CS population for CHR. However, CHR can also exist as a (sub)proportion of CR isolates. Additionally, we only focused on CHR as a phenomenon in which there is a minor subpopulation with a MIC above the breakpoint in a major population with a MIC below the breakpoint. However, HR may also occur in entirely susceptible populations [
30].
Given the large number of isolates screened, this study is a step forward in elucidating the prevalence and burden of CHR in common ST lineages of K. pneumoniae. Our data prompt for development of more robust and simple diagnostics to enable implementation of HR detection on a larger scale, and for more structured studies to quantify the actual impact of CHR on treatment failures in patients receiving colistin.
Figure 1.
Classification schemes for CHR. Figure describes the growth requirements, more specifically the frequency of growth required, for an isolate to be determined CHR. In case there is no frequency requirement at a specific concentration of colistin but only the requirement that there is visible growth on the plate, this has been indicated as “Growth”. For classification 2 the requirements to be fulfilled depend on the MIC for colistin of the isolates. NR = no requirement, NA = not applicable. 1 Plates containing 16 mg/L of colistin were only included for isolates with a MIC of 2 mg/L. 2 Frequency of ≥1x10-6 only required for either 4 mg/L or 8 mg/L, not for both though allowed Created with BioRender.com.
Figure 1.
Classification schemes for CHR. Figure describes the growth requirements, more specifically the frequency of growth required, for an isolate to be determined CHR. In case there is no frequency requirement at a specific concentration of colistin but only the requirement that there is visible growth on the plate, this has been indicated as “Growth”. For classification 2 the requirements to be fulfilled depend on the MIC for colistin of the isolates. NR = no requirement, NA = not applicable. 1 Plates containing 16 mg/L of colistin were only included for isolates with a MIC of 2 mg/L. 2 Frequency of ≥1x10-6 only required for either 4 mg/L or 8 mg/L, not for both though allowed Created with BioRender.com.
Figure 2.
PAP assay results of the control strains. Graph represents the log10 CFU/ml per concentration of colistin used in the CAMHA plates of the PAP assay. (A) Individual results for each strain for each run, graph illustrates the intra-run variation for the different control strains. (B) Average result for each control strain, graph illustrates the overall result of the control strains.
Figure 2.
PAP assay results of the control strains. Graph represents the log10 CFU/ml per concentration of colistin used in the CAMHA plates of the PAP assay. (A) Individual results for each strain for each run, graph illustrates the intra-run variation for the different control strains. (B) Average result for each control strain, graph illustrates the overall result of the control strains.
Figure 3.
Distribution of isolates across ST types. Graphs show the number of isolates per ST type as well as the number of CR/HR and CS/non-CHR isolates. Only ST types having a statistically significant association with CR and/or CHR are shown. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CR per ST type, (B) CHR per ST type, (C) CR per ST type per country, Slovenia (n=9) and United Kingdom (n=23) were not taken into further consideration due to the low number of isolates, (D) CHR per ST type per country.
Figure 3.
Distribution of isolates across ST types. Graphs show the number of isolates per ST type as well as the number of CR/HR and CS/non-CHR isolates. Only ST types having a statistically significant association with CR and/or CHR are shown. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CR per ST type, (B) CHR per ST type, (C) CR per ST type per country, Slovenia (n=9) and United Kingdom (n=23) were not taken into further consideration due to the low number of isolates, (D) CHR per ST type per country.
Figure 4.
Distribution of isolates across countries. Graphs show the number of isolates per country as well as the number of CR/HR and CS/non-CHR isolates. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CR per country, (B) CHR per country.
Figure 4.
Distribution of isolates across countries. Graphs show the number of isolates per country as well as the number of CR/HR and CS/non-CHR isolates. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CR per country, (B) CHR per country.
Figure 5.
Distribution of isolates across MIC values. Graphs show the number of isolates per MIC as well as the number of CHR and non-CHR isolates. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CHR per MIC value, (B) CHR per MIC value per country, (C) CHR per MIC value per ST type. Only ST types having a statistically significant association with CR and/or CHR are shown.
Figure 5.
Distribution of isolates across MIC values. Graphs show the number of isolates per MIC as well as the number of CHR and non-CHR isolates. Of note C1+C2 represents the total amount of isolates fulfilling C2 whilst C2 alone represents isolates only fulfilling C2. (A) CHR per MIC value, (B) CHR per MIC value per country, (C) CHR per MIC value per ST type. Only ST types having a statistically significant association with CR and/or CHR are shown.
Table 1.
Table contains a detailed breakdown of the PAP assay results (n=288) including the reasons why isolates did not fulfil the definition of CHR and the number of isolates per observed result. Of note, for some isolates that did not fulfil the definition of CHR but did show growth >2 mg/L there were multiple reasons why they were not classified as CHR. Sub-reasons are listed in order of importance. Each isolate was only included once and if it fulfilled multiple sub-reasons it was only included in the sub-reason considered most important.
Table 1.
Table contains a detailed breakdown of the PAP assay results (n=288) including the reasons why isolates did not fulfil the definition of CHR and the number of isolates per observed result. Of note, for some isolates that did not fulfil the definition of CHR but did show growth >2 mg/L there were multiple reasons why they were not classified as CHR. Sub-reasons are listed in order of importance. Each isolate was only included once and if it fulfilled multiple sub-reasons it was only included in the sub-reason considered most important.
Observed results |
No. of isolates |
Fulfilling classification 1 |
Growth on the plates containing 4 and/or 8 mg/L of colistin with a frequency of at least 1•10-6
|
25 |
Fulfilling only classification 2 |
Growth on at least all the plates containing colistin at a concentration up to and including 8-fold the MIC of the isolate at a frequency of minimally 1•10-7, minimum concentration at which there should be growth was 4 mg/L |
83 |
Not fulfilling either classification but growth >2 mg/L |
Growth on 4 and/or 8 mg/L plate but frequency <1•10-7
|
45 |
For MIC 0.0625-0.5 mg/L: growth with frequency ≥1•10-7 on 8 mg/L plate but <1•10-7 on 4 mg/L plate |
8 |
For MIC 1-2 mg/L: growth with frequency ≥1•10-7 on 4 and/or 8 mg/L plate but <1•10-7 on plates ≥8-fold the MIC |
30 |
Growth with frequency ≥1•10-7 on plates ≥8-fold the MIC but frequency of 1•10-7 not reached on all plates <8-fold the MIC |
19 |
No growth at 4 and 8 mg/L |
|
71 |
No growth at 1, 2, 4 and 8 mg/L |
|
7 |
Table 2.
Summary of mutations found in CR and CHR isolates. Table contains information on the mutations found in the resistant subpopulation of two confirmed CHR isolates (AN1505CP2 and IT0244CP) at different concentrations of the PAP assay plates on three separate assays. Additionally, table contains information on three CR isolates with the same ST type as IT0244CP (ST409). For AN1505CP2 there were no CR isolates with the same ST type (ST323).
Table 2.
Summary of mutations found in CR and CHR isolates. Table contains information on the mutations found in the resistant subpopulation of two confirmed CHR isolates (AN1505CP2 and IT0244CP) at different concentrations of the PAP assay plates on three separate assays. Additionally, table contains information on three CR isolates with the same ST type as IT0244CP (ST409). For AN1505CP2 there were no CR isolates with the same ST type (ST323).
Isolate ID |
MIC(mg/L) |
ST type |
PAP assay plate conc. (mg/L) |
Mutations in mgrB
|
Other mutations |
IT0307CP(CR) |
128 |
ST409 |
|
IS1R of IS1 family interruption at nt 107 |
|
IT0636C(CR) |
128 |
ST409 |
|
ISKpn34 of IS3 family interruption at nt 46 |
|
IT0915C(CR) |
64 |
ST409 |
|
IS903B of IS5 family interruption at nt 34 |
|
IT0244CP(CHR1st PAP) |
0.5 |
ST409 |
2 |
ISKpn34 of IS3 family interruption at promoter |
|
8 |
IS903B of IS5 family interruption at nt 117 |
|
16 |
IS1S of IS1 family interruption at promoter |
|
IT0244CP(CHR3rd PAP) |
0.5 |
ST409 |
2 |
IS1X2 of IS1 family interruption at nt 123 |
|
AN1505CP2(CHR1st PAP) |
1 |
ST323 |
4 |
Deleted |
|
8 |
IS903B of IS5 family interruption at nt 70 |
|
AN1505CP2(CHR2nd PAP) |
1 |
ST323 |
8 |
|
pmrB: T157P |
AN1505CP2(CHR3rd PAP) |
1 |
ST323 |
8 |
Q30X |
|
16 |
Deleted |
|