Klebsiella variicola Is a Frequent Cause of Bloodstream Infection in the Stockholm Area, and Associated with Higher Mortality Compared to K. pneumoniae

1Klebsiella pneumoniae is a gram-negative pathogen causing a wide spectrum of both hospital- and community-acquired infections such as urinary tract infection, pneumonia, intraabdominal infection, bloodstream infection (BSI), meningitis and pyogenic liver abscess (PLA) –. The mortality in invasive infection is high, ranging between 17.5 and 23% –.

2Recently, K. pneumoniae has emerged as an increasingly resistant pathogen; it has shown an unprecedented ability to express several intrinsic and acquired mechanisms making the species frequently multidrug-resistant (MDR) to clinically important antimicrobial classes. The clonal dissemination of resistant strains has been the focus of significant attention because of the lack of effective treatment options,. Therefore, the dramatic increase in antibiotic-resistant outbreaks faced in several hospital and geographic locations has a high potential of causing mortality and morbidity –.

3K. pneumoniae displays a large spectrum of virulence factors associated with the infective potential of community-acquired-isolates –. The bacterial capsule is important for its virulence, with 78 capsular types identified,. Among them K1, K2, K5, K20, K54, K57 were associated with invasive disease,,,. Other potential pathogenic factors are the lipopolysaccharide,, iron scavenging systems, fimbrial and non-fimbrial adhesion factors,, hypermucoviscosity,, the plasmid-borne rmpA (regulator of mucoid phenotype A), wcaG; encoding the capsular fucose synthesis, which may enhance the antiphagocytic activities, and allS; encoding the activator of allantoin regulon.

4Klebsiella pneumoniae isolates display metabolic versatility, enabling bacteria of this genus to thrive in a variety of environmental niches,,,.

5Phylogenetic analysis of clinical, carriage and environmental isolates classically identified as K. pneumoniae demonstrated the existence of three main phylogenetic lineages (phylogroups) called KpI, KpII and KpIII as initially demonstrated based on the sequence analysis of gyrA and parC genes, and later by their association with specific families of chromosomal β-lactamase genes,. From a taxonomical standpoint, Klebsiella pneumoniae comprises three subspecies,,,: K. pneumoniae subsp. pneumoniae, K. pneumoniae subsp. ozaenae and K. pneumoniae subsp. rhinoscleromatis. The two latter subspecies are rarely encountered and are associated to specific diseases (rhinoscleroma and ozena, respectively). From a genetic viewpoint, these two subspecies represent homogeneous genotypic groups (MLST clonal complexes; CC) that are nested within the main phylogroup of K. pneumoniae, named KpI. Subsequent taxonomic work has proposed the names Klebsiella variicola, for phylogroup KpIII and K. quasipneumoniae for phylogroup KpII. As a consequence, the name K. pneumoniae should now be used only for strains that belong to phylogroup KpI (K. pneumoniae sensu stricto).

6It is difficult to distinguish K. pneumoniae from K. quasipneumoniae and K. variicola by biochemical tests,. Until now, K. variicola and K. quasipneumoniae strains have therefore been generally misidentified as K. pneumoniae in clinical microbiology laboratories. Currently these three species can be more reliably differentiated by genotyping methods,,. The introduction of mass spectrometry might offer a possibility of rectifying this problem in the future. As a result, the clinical importance of K. variicola and K. quasipneumoniae is currently underestimated, and there are few studies that have reported these species from clinical samples,. K. variicola isolates were frequently isolated from various plants. Previous investigations have shown that approximately 20% of human isolates thought to be K. pneumoniae are in fact K. variicola/KpIII or K. quasipneumoniae/KpII,.

7A number of studies have shown association between pathogenicity of K. pneumoniae and virulence genes, serotypes, sequence types and mucoid phenotypes,–. The purpose of this study was to characterize clinical invasive isolates of K. pneumoniae sensu lato (i.e., in the classical sense and including K. variicola, K. quasipneumoniae and K. pneumoniae sensu stricto) from the Stockholm area and to analyze their population structure. A specific focus was to identify the main phylogroups. Further, our aims were to analyze their clonal diversity; to describe the association of virulence factors and serotypes with phylogroups; to study the association between phylogroups, bacterial traits and severity of disease, with primary endpoint being 30-day mortality.

8All adult (≥18 years old) patients admitted to Karolinska University Hospital, Solna, Sweden, between 2007 and 2009, with growth of K. pneumoniae either in blood (n = 137) or cerebrospinal fluid (CSF) (n = 1) or both (n = 2) were included in this retrospective cohort study. Isolates and patients were detected by searches in the clinical microbiology laboratory information system. Species identification was done with API 20E system (bioMérieux, Marcy l’Etoile, France) or VITEK2 (bioMérieux). Antimicrobial susceptibility testing was performed with the disk diffusion method on Isosensitest agar (Oxoid, Basingstoke, UK) interpreted according to the guidelines of the Swedish Reference Group for Antibiotics (SRGA). Genotypic analysis of isolates producing extended-spectrum β-lactamases and/or carbapenemases was performed with Check-MDR (Checkpoints, Wageningen, The Netherlands). One isolate was excluded due to non-typability with MLST, hence the total amount of included isolates were 139. One blood or CSF isolate from each patient was used for the subsequent analyses.

9Detection of serotypes and virulence genes was performed by PCR,. We sought for six major serotypes that have been reported to be strongly associated with community-acquired invasive disease: K1, K2, K5, K20, K54, and K57. The virulence genes allS, rmpA, and wcaG were also investigated as described previously,. Hypermucoviscous phenotype was determined using the string test, i.e. by pulling colonies with an inoculation loop, using a threshold of 5 mm.

10MLST was carried out as previously described,. Updated details of protocol and primers for amplification and sequencing are given on the web site (http://bigsdb.web.pasteur.fr). Templates were sequenced on both strands with the published primers using the BigDye Terminator Ready Reaction Mix v3.1. Nucleotide sequences were determined by ABI Prism 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA). New allelic variants were repeated and confirmed in triplicate. Alleles and profiles were determined by direct comparison with the K. pneumoniae MLST database. The alignment of sequences was performed with the Mega v5. Neighbor-joining trees were drawn for concatenate sequence of the seven MLST loci for the total population. To identify phylogenetic groups, three reference strains (group KpI: K. pneumoniae ATCC13883: ST3; group KpII: SB59: ST1118; and group KpIII: Kp342: ST146) were included in phylogenetic analyses. The phylogenetic structure was also assessed by drawing separately a Neighbor-joining tree for each allele. Phylogenetic groupings were identified also by NeighborNet analysis based on the concatenated sequences using SplitsTree v4.

11DnaSP v5 was used to calculate the nucleotide diversity of each gene for each phylogroup and for the whole population. To provide a graphical representation of the populations, a minimal spanning tree (MST) was produced with BioNumerics v7.0 (Applied Maths, Sint Maartens-Latem, Belgium). Clonal complexes were defined as groups of isolates sharing six loci of their allelic profiles with at least one other member of the group, using the program eBURST v3.

12The medical records were retrieved for the 139 cases, studying mortality, patient risk factors, hospital- versus community-acquired disease and antibiotic treatment. Charlson comorbidity index were constructed to assess comorbidity. Infection was classified as polymicrobial if at least one more different species was recovered from specimens drawn within 24 h from the K. pneumoniae bloodstream infection. Skin contaminants (i.e., Corynebacterium spp., Bacillus spp., Propionibacterium spp., and coagulase-negative staphylococci [CoNS]) had to be present in at least two blood cultures, while for other species presence in one blood culture was sufficient. Infections were defined as hospital-acquired if the sample was taken>48 h after admittance to hospital or if the patient had been admitted to hospital within the previous 30 days.

13Fisher’s test and Chi-squared test were used as appropriate for categorical variables. The Mann Whitney test was used to compare continuous variables. For all tests a two-sided p-value <0.05 was considered significant. Odds Ratio (OR) and 95% Confidence Intervals (CI) were calculated using logistic regression in the STATA 12.0 software. The multivariate analysis was restricted by the relatively few outcomes (mortality in analyzed groups, n = 23). Since our primary interest was to investigate the effect of phylogeny on mortality risk we decided to build a first model by including age and Charlson comorbidity index to adjust for host factors (these have previously been correlated to severe outcome in patients with severe bacterial disease). Both were adjusted for as continuous variables. In a second model we further included polymicrobial infection (present/absent) and gender. To avoid adjusting for factors that could be on the causal pathway between phylogroup and mortality, we did not adjust for disease manifestation (severity of disease and location). Time to adequate antibiotic therapy could be regarded as being on the causal pathway since it is likely influenced by disease manifestation. However, since it is often included in multivariate models of bacterial strains/clones and mortality we included time to adequate therapy as a continuous variable in a third model for comparison.

14Ethical approval for the study was obtained from the Karolinska Institutet Regional Ethics Committee of Stockholm (recordal 2009/1985-31/4). The committee approved that no written or verbal consent had to be given by the study subjects, as the study only pertained to extracting limited clinical data from patient charts. Patient records were anonymized and de-identified prior to analysis.

15The phylogenetic analysis of the concatenated sequences of the 7 MLST loci (Fig. 1A, Fig. S1) clearly revealed a well-defined structure, with three phylogroups that were each supported by maximal bootstrap values. The obtained phylogeny corroborates previous findings,,. The three clades in the figures correspond to phylogroups KpI, KpII and KpIII (Fig. 1A, B, Fig. S1). KpI was the largest phylogroup consisting of 96 isolates, followed by KpIII (Klebsiella variicola) consisting of 34 isolates, and KpII (K. quasipneumoniae) consisting of 9 isolates. As previously noted, KpI and KpII were branched together while KpIII was branched in an external position (Fig. 1A). The mean genetic distance between clusters was 4.0% (KpI-KpII); 4.4% (KpI-KpIII) and 4.5% (KpII-KpIII). Furthermore the three phylogroups were confirmed by conducting a Neighbor-joining tree analysis with each of the individual alleles (data not shown). The NeighborNet graph (Fig. 1C) also demonstrated a clear demarcation among the three phylogroups and suggested a high level of recombination within each of the groups as well as among them.

16The alignment of the seven housekeeping genes revealed the presence of minor events of insertion/deletion (INDEL) only in the tonB locus. The polymorphism parameters were assessed in the whole population, and in the three phylogroups (Table 1). Considering all isolates, the polymorphic sites ranged from 5.1% (gapA) to 19.6% (tonB), whereas the concatenate sequence revealed 10.8% polymorphic sites (Table 1). The average nucleotide diversity π was 2.1% for all the population. However, the internal level of π within the three phylogroups was much lower and the diversity within KpI (0.37%) and KpIII (0.63%) was lower than for KpII (1.3%). A ratio of Ka/Ks <1 was observed, as typical for housekeeping genes evolving predominantly under purifying selection.

17The population analysis based on the MLST data showed an extensive haplotypic diversity, with 116 distinct STs. The eBURST and MST analyses (Fig. 1B) were both used to depict the clonal relatedness among sequence types (STs). The 139 isolates were subdivided into 12 clonal complexes comprising 48 isolates (Fig. 1B, Table 2). The remaining 91 isolates were divided into 79 singletons (STs that differed by two or more alleles from any other ST).

18Among the 139 isolates (Fig. S2B, Table S1) only 18 strains (12.9%) were serotypable for the investigated capsular serotypes; K1 (1.4%; 2/139), K2 (5.0%; 7/139), K20 (1.4%; 2/139), K54 (2.2%; 3/139); K57 (2.9%; 4/139). K2 was the most frequently occurring serotype and K5 was not detected in any of the isolates. Of the 18 capsular types found in this study, 16 belonged to KpI and 2 to KpIII.

19Regarding virulence genes (Fig. S2C, Table S1), rmpA was positive for 5 strains (3.6%; 5/139; ST25, 380, 592, 593 and 701), wcaG was positive for 12 strains (8.6%; 12/139; ST14, 29, 29, 581, 599, 604, 605, 616, 637, 695, 699 and 700) and allS was positive for 6 strains (4.3%; 6/139; ST104, 604, 605, 613, 620 and 702). A mucoid phenotype (Fig. S2D, Table S1) was detected in 8 isolates (5.8%; 8/139; ST20, 111, 268, 280, 380, 504, 542 and 592) out of which two were also positive for rmpA. With the exception of ST380, which was previously associated with rmpA, there were no other observed correlations between virulence genes and STs or clonal complexes. The 21 isolates with virulence genes were distributed as follows: KpI (n = 13), KpIII (n = 5) and KpII (n = 3).

20Among the 139 infectious episodes, 24 events were associated with a fatal outcome within 30 days translating to a total 30-day mortality of 17.3% (Table S1). Isolates belonging to phylogroup KpIII had the highest 30-day mortality (29.4%, 10/34 patients; Fig. S2A, Table S1). In phylogroup KpI and KpII the 30-day mortality was 13.5% (13/96) and 11.1% (1/9) (Fig. S2A, Table S1), respectively. In further comparisons of the association between phylogroup and mortality phylogroup KpII was excluded due to the small amount of patients. There was a significant difference in 30-day mortality between KpIII and KpI in crude analyses (OR 2.66 (95% CI: 1.04–6.82, p = 0.037) (Table 3). Increasing Charlson comorbidity index was strongly associated with mortality (p = 0.005). Metastatic cancer, which attributes 6 points to the Charlson comorbidity index, was also associated with mortality (p = 0.001). After adjustment for age and Charlson comorbidity index, KpIII remained associated with increased mortality (OR: 3.03 (95% CI: 1.10–8.36). Furthermore, after including age, Charlson comorbidity index, presence of polymicrobial infection and gender, as well as inclusion of time to adequate antibiotic therapy, KpIII was associated with increased mortality [OR: 3.11 (95% CI: 1.10–8.80) and 3.23 (95% CI: 1.12–9.27) respectively]. In a final analysis we excluded 7 patients that were treated at the intensive care unit at the time of onset of the BSI, which did not alter our results (KpIII adjusted OR: 3.53 (95% CI: 1.20–10.39).

21Among the 24 events associated with 30-day mortality only one isolate was serotypable (K54), and strikingly few (three isolates) featured virulence genes. Almost all isolates belonged to new sequence types, with the exception of ST15 (forming a clonal complex with ST14, which is prevalent among carbapenemase-producers),,; ST37 (previously described in bloodstream infection); ST337 (previously described in bloodstream infection in KPC-producing isolate, Colombia),; ST359 (previously described in bloodstream infection, colistin-resistant, South Korea) and ST280 (previously described in bovine mastitis) (http://bigsdb.web.pasteur.fr). Table 4 summarizes patient characteristics. A majority of the patients (62.6%) were male, and the median age was 70 years. Comorbidity measured in Charlson comorbidity index ranged from 0 to 9. Most of the patients, 76/139 (54.7%) had an index of 2–3, 18 (24.5%), 15 (10.8%) and 30 (21.6%) had Charlson comorbidity index of 0–1, 4–5 and>5 respectively. The most common disease was having a malignancy with a prevalence of 81/139 (58.3%). Hematological malignancy was most common (n = 26) followed by urogenital malignancies and colorectal cancer (n = 25 and n = 11, respectively). 68/139 (48.9%) of the infectious episodes were considered hospital-acquired. Polymicrobial infections occurred in almost one fourth of the patients (34/139; 24.5%). The most frequent co-pathogen in the polymicrobial infections, observed in 15/139 (10.8%) patients, was Escherichia coli followed by Staphylococcus aureus n = 5 (3.6%) and Enterococcus spp, n = 5 (3.6%) (E. faecalis n = 2, E. faecium n = 2 and E. casseliflavus n = 1) respectively. There were no significant differences between the phylogenetic groups concerning patient characteristics (Table 4).

22The most common antibiotic treatment, in 58.3% (81/139) of the patients, was cephalosporins (cefuroxime, cefotaxime or ceftazidime) (Table S1). Twenty-six patients (18.7%) received piperacillin-tazobactam as empiric treatment, while 13.7% (19/139) and 6.5% (9/139) started with a carbapenem (imipenem or meropenem) and ciprofloxacin, respectively. Aminoglycosides (gentamicin, amikacin) were added to the first dose of empiric treatment in 23/139 cases (16.5%) and in 4 of these cases this component was the only active against the K. pneumoniae isolate.

23Within 4 hours, adequate antibiotic treatment was given to 77/139 patients (55.4%), 44/139 received adequate treatment within 2 hours (31.7%), and 15/139 (10.8%) within 1 hour (Table S1). After 24 hours 121/139 patients (87.1%) had received appropriate treatment, and after 48 hours most patients, 134/139 (96.4%) had received adequate antibiotic treatment.

24Overall there was a low level of acquired antibiotic resistance (Table 5). The resistance levels against antimicrobials with relevant activities were as follows: trimethoprim/sulfametoxazole 10.8%, ciprofloxacin 8.6%, piperacillin/tazobactam 5.8%, gentamicin 3.6%, ceftazidime 4.3% and cefotaxime 3.6%. The resistance against trimethoprim/sulfametoxazole was significantly higher in KpI than compared to KpIII (P = 0.014). A total of 15.1% of the isolates were resistant to one or several of these antibiotic classes, 18.8% in KpI and 5.9% in KpIII (P = 0.078). These finding are similar to previous reports, in which K. variicola isolates have been found less resistant than KpI and KpII for several classes of antimicrobial agents.

25Five of 139 isolates (3.6%) were considered multidrug-resistant (resistant to ≥three of the following antibiotic classes; trimethoprim/sulfametoxazole, fluoroquinolones, piperacillin/tazobactam, cephalosporins, aminoglycosides or carbapenems). These isolates all produced extended-spectrum β-lactamases (ESBL); two isolates CTX-M group 1 (ST15), one isolate CTX-M- group1 and an SHV-ESBL (ST340) and one isolate an SHV-ESBL (ST17). The fifth isolate (ST383) co-produced a CTX-M group 1 enzyme, CMY II, and the carbapenemase VIM. The patient had previously been hospitalized in Greece. All ESBL-producing isolates belonged to KpI. Among the 4 patients with carbapenem susceptible ESBL-producing isolates, two received adequate empiric antibiotic treatment within 4 hours (carbapenems), in the other two cases there was a delay between 24 and 72 hours until adequate treatment was received (no death occurred). The patient with BSI caused by a VIM-producing strain received adequate antibiotic treatment, colistin, rifampicin and meropenem, after 4 days and had a fatal outcome. However this patient also had significant comorbidities requiring ICU-care, respirator and dialysis before the onset of the bloodstream infection.

26To our knowledge this study is the first merging phylogroup classification, clinical perspective, virulence factors and sequence types on an inpatient consecutive material of bloodstream infections caused by Klebsiella pneumoniae (sensu lato) during a long time period and in a well-defined geographical area. The population structure of K. pneumoniae has been suggested to consist of three phylogroups but this has not been thoroughly explored in specific geographical locations and with respect to clinical outcome and patient characteristics. According to our findings, the Swedish K. pneumoniae population comprises the three main phylogroups (Fig. 1A, B, C). KpI was the most frequent group consisting 96 isolates (69%), followed by KpIII consisting of 34 isolates (24%) and KpII consisting of 9 isolates (6%). Previously, 420 isolates from 26 European hospitals were classified as follows: 345 isolates (82.1%) identified as KpI, 29 (6.9%) as KpII and 46 (11%) as KpIII. Hence, the population in this study was comparable, in terms of phylogroup relative frequencies, to a Europe-wide sampling. Phylogroups KpII and KpIII were given species status as K. quasipneumoniae and K. variicola,. These observations confirm that K. pneumoniae, as defined biochemically and identified in clinical microbiology laboratories, must be regarded as a complex of three species. These three species appear to have wide geographic distribution and to be represented, among isolates initially identified as K. pneumoniae, in similar relative order of frequency in distinct geographical locations,.

27Phylogroup KpI has been shown to be overrepresented among clinical K. pneumoniae isolates in previous studies,. Contrary to previous studies, phylogroup KpIII (K. variicola) isolates were relatively more common in this study. Recently Seki et al. described a case of fatal bloodstream infection caused by K. variicola. Initially, based on the classic and automated identification methods, the isolate was incorrectly identified as K. pneumoniae, but reclassified as K. variicola according to data derived from next-generation sequencing. As it is difficult to distinguish K. variicola from K. pneumoniae by classical methods used in clinical laboratories for species determination, K. variicola has been underreported in the literature. K. variicola was initially shown to be adonitol negative unlike K. pneumoniae,, but this test alone cannot reliably identify K. variicola,. Utilization of 5-keto-D-gluconate was recently found to be specific of K. variicola and may represent a useful feature for identification. Also, the usefulness of mass spectrometry for species identification of K. variicola remains to be determined.

28Here we addressed for the first time to our knowledge, the possible presence of a link between phylogroup and bacterial traits (virulence factor, serotypes) or severity of disease. The overall 30-day mortality was 17.3%, similar to some of the previous studies on BSI caused by K. pneumoniae –. Analyzing phylogroups in relation to mortality, KpIII had the highest 30-day mortality (10/34, 29.4%). To the best to our knowledge, this is the first report indicating a relatively high mortality associated with K. variicola. The difference in mortality between KpI and KpIII was significant (P = 0.037) both in univariate analysis and in multivariate analysis adjusting for age, comorbidity, gender, time to adequate antimicrobial therapy, and polymicrobial infections. The difference in mortality between the clades cannot be explained by differences in clinical characteristics of patients since they were similar (Table 4). Measuring time to adequate antimicrobial treatment and other supportive therapy is complicated in this type of study as severity of disease and presentation of disease often influences choice and timing of empirical treatment. However as resistance rates were low in our study more than half of the patients (55.4%) received adequate antibiotic treatment within 4 hours. After 24 hours most patients (87.1%) had received adequate antibiotic treatment, and there were no differences among the survivors and the deaths in time to adequate antibiotics.

29However there are limitations to consider. The major concern is the small number of events (30-d mortality of 23 patients in the two analyzed phylogroups), causing limitations in the multivariate model. Despite this concern, it is interesting that our data suggest that the underestimated pathogen K. variicola is associated with poor prognosis in BSI. The 30-d mortality is often used as an endpoint in studies measuring mortality. However, when studying BSI and mortality, there is always a risk that the cause of death is not attributable to the BSI, although this risk would be regarded equally large in all phylogroups. It could be debated whether other variables could be adjusted for, but the risk of this strategy is that such variables could be on the causal pathway between phylogroup and mortality.

30Numerous hospital-based studies have notified that several comorbidities such as malignancy, diabetes mellitus, cirrhosis, biliary tract disorders and alcoholism may impair patient defenses and hence contribute to the development of K. pneumoniae infection,,,–. Previous studies on BSI caused by K. pneumoniae have shown that the disease often affects older patients with high comorbidity. Our study corroborates these earlier findings. Most of the patients in our study were male with a median age of 70, 32.4% had Charlson comorbidity index 4 or more, and almost half of the episodes were hospital acquired. Increasing Charlson comorbidity index was strongly associated with mortality (P = 0.005). Advanced age was associated with an increased, but not significant, risk of death. There was a similar amount of hospital-acquired infections among the deaths and survivors.

31In accordance to what was reported previously,, our study demonstrates that isolates belonging to the phylogroup KpI were found to be more resistant than K. variicola isolates. From the present study, we retrieved only five MDR isolates that were ESBL-producers and belonged to KpI. Recently, Valverde et al. determined the phylogenetic structure of a collection of ESBL-producing K. pneumoniae isolates. Expectedly, most of the ESBL-producing isolates belonged to KpI. A minority of isolates was classified in the KpIII group, and almost all of them were CTX-M-10 producers. Lastly, the prevalence of class 1 integrons was screened among these three phylogroups; the class 1 integrons were dominantly associated with KpI but very few were observed in isolates of KpIII. As suggested previously, the high prevalence of KpI isolates among clinical isolates could be a reflection of the higher antimicrobial resistance rates within this group.

32This collection of K. pneumoniae isolates contained few virulence factors compared with reports from previous studies, based on community-acquired BSI and severe disease. Mucoid phenotype (n = 8), virulence genes rmpA, wcaG and allS (n = 21 isolates; two of them had both wcaG and allS), and capsular serotypes K1, K2, K5, K20, K54 and K57 (n = 18) were all relatively rare in the collection. Much effort has been made in the last decades to estimate the link between the capsular types, virulence potential, and the nature of diseases, but mostly with community-acquired isolates, which differ in ST and virulence factors from nosocomial isolates,,,,. Within the phyologroups in our dataset, there was no correlation between the serotypes or virulence genes and CCs or STs.

33Previous reports have expressed concern about numerous bacterial virulence genes which were significantly associated with invasive strains. Herein, we show that rmpA gene was associated only to K2, K20 and K57, and not to K1, K5 and K54 as found in previous reports,. In the same reports wcaG, encoding capsular fucose production, was associated with capsular types K1 and K54,, which was further corroborated by our present results. Although the presence of the wcaG gene has not been well studied, there is data to suggest that strains harbouring this gene are more often observed in severe disease. Interestingly, two fatal strains harboring the wcaG gene were recorded in our study (Table S1); ST616 (K54, wcaG, KpIII) and ST604 (wcaG/allS, KpII). Conversely, the six isolates harboring the allS gene did not belong to any of the investigated capsular types. Noteworthy, this gene was not universally associated to K1 serotype, but considered as a specific marker for the hypervirulent clone (CC23 K1) isolated from patients with pyogenic liver abscesses detected in several countries,,,,,–. In a recent report from France, five cases of fatal BSI were attributed to K. pneumoniae serotype K2, possessing the rmpA gene. The isolates belonged to two clones; ST86 (two isolates) and ST380 (3 isolates). The same genotype (ST380, K2, rmpA) was encountered in our study but was not associated with a fatal outcome. In the 30-day mortality group (24/139 patients) only three isolates featured one or more virulence factors. In this study, where both community- and hospital-acquired BSIs were included, no indication for an existing correlation between above mentioned virulence factors and increased risk of mortality in BSI caused by K. pneumoniae could be found.

34We have analyzed the population structure of invasive clinical isolates of K. pneumonia sensu lato in Sweden, and found K. pneumoniae (KpI), but also K. quasipneumoniae (KpII) and K. variicola (KpIII). The 30-d mortality was more commonly associated with infecting strains belonging to K. variicola, but was not associated with known virulence factors, including some of the most virulent capsular types. More studies with enhanced species identification would be valuable to expand the data on the clinical importance of K. variicola in relation to K. pneumoniae sensu stricto. In general, a high level of comorbidity, equally high among the patients with BSI caused by K. pneumoniae as among the patients with BSI caused by K. variicola, was observed.