- Research article
- Open Access
Diversity of virulence level phenotype of hypervirulent Klebsiella pneumoniae from different sequence type lineage
BMC Microbiology volume 18, Article number: 94 (2018)
Hypervirulent Klebsiella pneumoniae (hvKP) is emerging around the Asian-Pacific region and it is the major cause of the community-acquired pyogenic liver abscesses. Multidrug-resistant hypervirulent Klebsiella pneumoniae (MDR-hvKP) isolates were reported in France, China and Taiwan. However, the international-ally agreed definition for hvKP and the virulence level of hvKP are not clear.
In this study, 56 hvKP isolates were collected from March 2008 to June 2012 and investigated by string test, capsule serotyping, multilocus sequence typing (MLST), virulence gene detection and serum resistance assay. Among the 56 K. pneumoniae isolates, 64.3% had the hypermucoviscosity phenotype, meanwhile, 64.3% were the K1 serotype and 19.6% were the K2 serotype. Within the K1 serotype, 94.4% were ST23, and within the K2 serotype, ST65, ST86 and ST375 accounted for the same percentage 27.3%. The serum resistance showed statistically normal distribution. According to the 50% lethal dose of Galleria. mellonella infection model, hvKP isolates were divided into high virulence level group and moderate virulence level group. The ability of each method evaluating the virulence level of hvKP was assessed using the area under the receiver operating characteristic curve.
K1 ST23 K. pneumoniae was the most prevalent clone of the hvKP. However, K1 ST23 K. pneumoniae was the dominant clone in the moderate virulence level group. MLST was a relatively reliable evaluation method to discriminate the virulence level of hvKP in our study.
Klebsiella pneumoniae, a well-known enteric gram-negative bacillus, often causes hospital-acquired infections, such as pneumonia, urinary tract and bloodstream infections . During the past two decades, a new hypervirulent variant of K. pneumoniae (hvKP) was reported. The hvKP can cause community associated liver abscesses, as well as possible complications including metastatic meningitis and endophthalmitis . The hvKP was usually associated with hypermucoviscosity and antibiotic-sensitivity.
The hvKP can cause serious infections in ambulatory and healthy hosts and result in metastatic spread of infections . The appearance of colonies grown on an agar plate is hypermucoviscous . According to the genomic data, in addition to iroBCD, iucABCD and rmpA/A2, Ye et al. found 21 new genes to be new virulence factors of K. pneumoniae. Several of these genes had been proposed as virulence factors in other bacteria, such as the gene encoding SAM-dependent methyltransferase and pagO which protected bacteria from phagocytosis .
Moreover, multidrug-resistant hypervirulent Klebsiella pneumoniae (MDR-hvKP) isolates had been reported. A K2 ST86 hvKP strain producing CTX-M-3 was reported in France. A K1 hvKP strain producing CTX-M-15 was reported in Taiwan. One K2 ST65 and five K1 hvKP carbapenem-resistant strains were reported in China [6,7,8,9].
Different terms of hvKP were used by investigators as there was no international-ally agreed definition. Basically, the definition of hvKP was based on the clinical manifestations, such as liver abscesses, metastatic meningitis and endophthalmitis . Considering the life-threatening nature of liver abscess and the diversity of the patients’ prognosis, discovering a predictable phenotype or gene to identify hvKP and the virulence level of hvKP is very critical. In this study, we determined the string test, serotyping, multilocus sequence typing, serum resistance and virulence-associated genes of 56 hvKP isolates. The virulence level of these isolates were verified by the LD50 of Galleria. mellonella infection model. Meanwhile, we compared the discriminating power of each method by the area under the receiver operating characteristic (AUROC) curve, the bigger AUROC was associated with better discriminative power.
G. mellonella haemolymph contains haemocytes, which function in a manner similar to phagocytes in mammals , therefore, the LD50 of G. mellonella larvae is a method used to verify the virulence level of hvKP. The LD50 values of all the isolates were distributed from 3.06 to 5.84 log10 CFU/ml. To analyse the characteristic phenotype and genes associated with virulence level, all the 56 isolates were separated into two groups, according to the density graph (Additional file 1: Figure S1). Accordingly, 27 isolates with an LD50 less than 5.06 log10 CFU/ml were defined as the high-virulence level (HVL), while the other 29 isolates fell into the moderate-virulence level (MVL) group. The results of statistical analyses between HVL and MVL groups were shown in Table 1.
To investigate whether the hypermucoviscosity phenotype can determine the virulence level of hvKP, all the isolates were verified by the string test. Among the 56 hvKP isolates, 36 isolates (64.3%) showed hypermucoviscosity. In the hypermucoviscosity isolates, 100% were rmpA positive, whereas in non-hypermucoviscosity isolates, 95% were rmpA positive (Fisher’s exact test, p = 0.3571). By verifying the virulence level, 70.4% (19/27) had the hypermucoviscosity phenotype in the HVL group and 59.6% (17/29) in the MVL group had the hypermucoviscosity phenotype. The distribution of the hypermucoviscosity phenotype was compared between the two groups, and there was no statistically significant difference (Fisher’s exact test, p = 0.4123).
To determine the relationship of serotype and MLST among hvKP, we serotyped the isolates and found that 64.3% (36/56) of the isolates were the K1 serotype, 19.6% (11/56) were the K2 serotype, and the rest were the non-K1/K2 serotype. While the dominant sequence type was ST23 in the K1 serotype (34/36), and the sequence type in the K2 serotype distributed more dispersedly. Within the K2 serotype, ST65, ST86 and ST375 shared the same percentage 27.3% (3/11), with others included one ST374 isolate and one ST380 isolate (Table 2). In the HVL group, 44.4% (12/27) of the isolates were the K1 serotype, and 29.6% (8/27) were the K2 serotype. In the MVL group, however, 82.8% (24/29) of isolates were the K1 serotype, and 10.3% (3/29) were the K2 serotype (Chi-square, p = 0.0112). Likewise, in the HVL group, 44.4% (12/27) of isolates were ST23, and in the MVL group, 75.9% (22/29) of isolates were ST23 (Fisher’s exact test, p = 0.0277). The distributions of serotype and sequence type between the HVL and MVL groups were statistically significantly different, respectively.
To evaluate the viability of hvKP in foetal bovine serum (FBS), a serum resistance test was performed. The viability of hvKP in FBS reflected its virulence level, the lower the value of serum resistance, the higher the proportion of bacteria surviving. The serum resistance of 56 hvKP isolates was distributed from 42.4 to 59.8%, and showed statistically normal distribution (Shapiro-wilk normality test, p = 0.4354). By verifying the virulence level, the serum resistance of the 56 hvKP isolates showed that the diversity of serum resistance between the two groups was not significantly different (Unpaired t test, p = 0.7779).
To determine the connection between the virulence level of hvKP and virulence-associated genes, magA, rmpA, kfu, fimH, wabG, uge, iroN, iutA, allS and entB were screened. Genes such as kfu, fimH and wabG existed in all the 56 isolates, while entB, iutA, iroN, rmpA were detected in more than 90% of the isolates. For magA, 44.4% (12/27) of the isolates in the HVL group were positive and 82.8% (24/29) in the MVL group were positive, which showed a statistically significant difference between the two groups (Fisher’s exact test, p = 0.0048). Regarding allS, in the HVL group, 51.9% (14/27) of isolates were positive, but in the MVL group, 89.7% (26/29) of isolates were positive, which showed a statistically significant difference between the two groups (Fisher’s exact test, p = 0.0214). There was no statistically significant difference between the two groups in genes like rmpA, uge, iutA, iroN and entB. For example, 92.6% (25/27) of the isolates in the HVL group and 89.7% (26/29) of the isolates in the MVL group were iutA positive (Fisher’s exact test, p > 0.9999).
To evaluate the discriminative power of all the methods above, receiver operating characteristic (ROC) curves were performed. MLST was found to be a relatively reliable evaluation method to determine the virulence level of hvKP (AUROC 0.7542, 95% confidence interval [95% CI] 0.6459–0.8625), followed by serotype (AUROC 0.6948, 95% CI 0.5745–0.8150) and magA gene (AUROC 0.6731, 95% CI 0.5542–0.7819) (Fig. 1a). The AUROC results of all the methods were shown in the Fig. 1b. Statistical analysis showed a significant difference between MLST and magA (p = 0.0132).
According to Fang et al., the string test, which reflects the hypermucoviscous phenotype, is the preferable laboratory-based surrogate marker for hvKP presently available and rmpA is important for hypermucoviscosity phenotype . However, in this study, 95% of non-hypermucoviscosity isolates were rmpA positive, evaluated that rmpA gene is not sufficient for hypermucoviscosity phenotype. Second, the hypermucoviscosity phenotype was not statistically different between the two virulence level groups. Therefore, the string test is not a sensitive method to assess virulence level in hvKP. This is in agreement with a recent review, which suggests that hypermucoviscosity is not an appropriate phenotype of hypervirulence, and the hypermucoviscosity of K. pneumoniae is not equal to hvKP .
The K1 serotype of K. pneumoniae isolates contributed to 64.3% of hvKP, which was higher than the previous report in mainland China (39.2% and 43%) [13, 14]. The prevalence of the K1 serotype was statistically higher than the K2 and non-K1/K2 serotype in this study. Our MLST data matched the Bialek-Davenet et al. report, wherein most isolates of serotype K1 belonged to Clone Group 23 (CG23), and STs associated with serotype K2 were distributed into 3 main CGs: CG86, CG375, and CG380 . In summary, the K1 ST23 isolates were the dominant clone of hvKP in China.
The K1 ST23 isolates occupied a higher proportion of the MVL group than the K2 serotype, which suggested that the K1 ST23 isolates had a moderate virulence level among the hvKP isolates. Compared to the AUROC of serotype, the AUROC of MLST showed a bigger area, which suggested that MLST was a relatively reliable method to evaluate the virulence level of hvKP. Further research will be needed to determine whether K1 ST23 isolates are in moderate virulence level than the other K. pneumoniae serotypes.
Virulence-associated genes were great for evaluating the hvKP, and some of the virulence-associated genes generally existed in the isolates of this study [13, 16]. In this study, virulence factors (magA, and allS) were highly associated with virulence level and were more prevalent in the MVL group. Other virulence factors showed no association with virulence level. Because of the complexity of the genetic background in clinical isolates, virulence-associated genes may not be a good choice to verify the virulence level of hvKP.
The string test, serum resistance assay, virulence-associated genes, serotype, and MLST are determined in hvKP isolates. In summary, the string test and serum resistance assay are not proper methods to assess hvKP. The virulence level of these isolates are diverse. The LD50 of G. mellonella is probably suitable for verifying the virulence level, but it is not simple or convenient. MLST is relatively reliable to evaluate the virulence level of hvKP at present. Clinical infections are complicated and determined by both the bacteria and the host. K1 ST23 hvKP isolates are dominant in the Asian-Pacific region. Our data indicates that the K1 ST23 serotype isolates are in relatively moderate virulence level.
Collection and identification of K. pneumoniae clinical isolates
From March 2008 to June 2012, a total of 56 isolates of K. pneumoniae were isolated from patients with liver abscess in the first affiliated hospital of Zhejiang University. These isolates were defined as hvKP according to Shon and Russo .The K. pneumoniae isolates included in this study were identified by 16S rRNA Sanger sequencing and the primers used were 27F (AGAGTTTGATCCTGGCTCAG) and 1492R (TACGGCTACCTTGTTACGACTT).
Hypermucoviscosity phenotypical identification of K. pneumoniae isolates
The hypermucoviscosity phenotypes of K. pneumoniae isolates were determined by the string test . Briefly, a positive string test is the formation of a viscous string > 5 mm in length when a colony is grown on a blood agar plate at 37 °C overnight and stretched by an inoculation loop.
Serum resistance assay
The serum resistance of each bacteria was determined by the method introduced by Mu et al. . Briefly, four independent cultures per strain were grown overnight in Mueller–Hinton (MH) broth, diluted 1:1000 [MH broth or 20% foetal bovine serum (FBS) MH Broth] and aliquoted into a honeycomb 100-well plate in three replicates. The plate was agitated at 37 °C. The OD600 of each culture was determined every 5 mins for 16 h using the Bio-Screen (BioScreen Testing Services, Inc., Finland). The growth rate was estimated by an R script based on the OD600 curves . The results were expressed as the ratio of the growth rate in 20% FBS MH broth and MH broth.
Detection of serotype, virulence-associated genes and multilocus sequence typing (MLST)
In the project, serotype-specific genes for the K1 and K2 capsular serotypes and virulence-associated genes including magA, rmpA, kfu, fimH, wabG, uge, iroN, iutA, allS and entB were all amplified by polymerase chain reaction (PCR). The primers used were shown in the Additional file 2: Table S1. MLST of K. pneumoniae was performed as previously described (http://bigsdb.pasteur.fr/klebsiella/primers_used.html).
50% lethal dose (LD50) of G.mellonella infection model
G. mellonella larvae were acquired from an insect mass-rearing plant in Shandong and kept in darkness at 4 °C. They were used within 12 days since reception. For LD50 experiments, a series of 10-fold serial dilutions of overnight incubation bacteria were diluted in PBS, of which the CFU was approximately 109 per millilitre. Then serial dilutions of K. pneumoniae in PBS were injected into G. mellonella larvae . Ten larvae were injected at each dilution and incubated at 37 °C in the darkness. The dead larvae were counted for three days after injection. For each isolate, the LD50 was calculated according to the probit model, and the results were expressed as log10 LD50 .
The Shapiro-Wilk normality test was performed to test the normality of the distribution of quantitative variables. For categorical variables, comparisons were performed by Chi-square analysis or Fisher’s exact test. All the quantitative variables that were normally distributed were presented as mean ± SD and compared by t-test. A p value < 0.05 was considered statistically significant. The ability of each method to discriminate the virulence level of K. pneumoniae was assessed by the area under the receiver operating characteristic (AUROC) curve, and the discriminative ability of the AUROC curves were compared by Delong’s test. The statistical software used was Prism5 (GraphPad Software, lnc.) and R 3.3.3 (pROC packages) for Windows.
Lery LM, Frangeul L, Tomas A, Passet V, Almeida AS, Bialek-Davenet S, Barbe V, Bengoechea JA, Sansonetti P, Brisse S, et al. Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor. BMC Biol. 2014;12:41.
Ko WC, Paterson DL, Sagnimeni AJ, Hansen DS, Von Gottberg A, Mohapatra S, Casellas JM, Goossens H, Mulazimoglu L, Trenholme G, et al. Community-acquired Klebsiella pneumoniae bacteremia: global differences in clinical patterns. Emerg Infect Dis. 2002;8(2):160–6.
Pomakova DK, Hsiao CB, Beanan JM, Olson R, MacDonald U, Keynan Y, Russo TA. Clinical and phenotypic differences between classic and hypervirulent Klebsiella pneumoniae: an emerging and under-recognized pathogenic variant. Eur J Clin Microbiol Infect Dis. 2012;31(6):981–9.
Shon AS, Bajwa RP, Russo TA. Hypervirulent (hypermucoviscous) Klebsiella pneumoniae: a new and dangerous breed. Virulence. 2013;4(2):107–18.
Ye M, Tu J, Jiang J, Bi Y, You W, Zhang Y, Ren J, Zhu T, Cao Z, Yu Z, et al. Clinical and genomic analysis of liver abscess-causing Klebsiella pneumoniae identifies new liver abscess-associated virulence genes. Front Cell Infect Microbiol. 2016;6:165.
Surgers L, Boyd A, Girard PM, Arlet G, Decre D. ESBL-producing strain of Hypervirulent Klebsiella pneumoniae K2, France. Emerg Infect Dis. 2016;22(9):1687–8.
Zhang Y, Zeng J, Liu W, Zhao F, Hu Z, Zhao C, Wang Q, Wang X, Chen H, Li H, et al. Emergence of a hypervirulent carbapenem-resistant Klebsiella pneumoniae isolate from clinical infections in China. J Inf Secur. 2015;71(5):553–60.
Zhang R, Lin D, Chan EW, Gu D, Chen GX, Chen S. Emergence of Carbapenem-resistant serotype K1 Hypervirulent Klebsiella pneumoniae strains in China. Antimicrob Agents Chemother. 2016;60(1):709–11.
Su SC, Siu LK, Ma L, Yeh KM, Fung CP, Lin JC, Chang FY. Community-acquired liver abscess caused by serotype K1 Klebsiella pneumoniae with CTX-M-15-type extended-spectrum beta-lactamase. Antimicrob Agents Chemother. 2008;52(2):804–5.
Kavanagh K, Reeves EP. Exploiting the potential of insects for in vivo pathogenicity testing of microbial pathogens. FEMS Microbiol Rev. 2004;28(1):101–12.
Fang CT, Chuang YP, Shun CT, Chang SC, Wang JT. A novel virulence gene in Klebsiella pneumoniae strains causing primary liver abscess and septic metastatic complications. J Exp Med. 2004;199(5):697–705.
Catalan-Najera JC, Garza-Ramos U, Barrios-Camacho H. Hypervirulence and hypermucoviscosity: two different but complementary Klebsiella spp. phenotypes? Virulence. 2017;8(7):1111–23.
Luo Y, Wang Y, Ye L, Yang J. Molecular epidemiology and virulence factors of pyogenic liver abscess causing Klebsiella pneumoniae in China. Clin Microbiol Infect. 2014;20(11):O818–24.
Shen DX, Wang J, Li DD. Klebsiella pneumoniae liver abscesses. Lancet Infect Dis. 2013;13(5):390–1.
Bialek-Davenet S, Criscuolo A, Ailloud F, Passet V, Jones L, Delannoy-Vieillard AS, Garin B, Le Hello S, Arlet G, Nicolas-Chanoine MH, et al. Genomic definition of hypervirulent and multidrug-resistant Klebsiella pneumoniae clonal groups. Emerg Infect Dis. 2014;20(11):1812–20.
Candan ED, Aksoz N. Klebsiella pneumoniae: characteristics of carbapenem resistance and virulence factors. Acta Biochim Pol. 2015;62(4):867–74.
Shon AS, Russo TA. Hypervirulent Klebsiella pneumoniae: the next superbug? Future Microbiol. 2012;7(6):669–71.
Yao B, Xiao X, Wang F, Zhou L, Zhang X, Zhang J. Clinical and molecular characteristics of multi-clone carbapenem-resistant hypervirulent (hypermucoviscous) Klebsiella pneumoniae isolates in a tertiary hospital in Beijing, China. Int J Infect Dis. 2015;37:107–12.
Mu X, Wang N, Li X, Shi K, Zhou Z, Yu Y, Hua X. The effect of Colistin resistance-associated mutations on the fitness of Acinetobacter baumannii. Front Microbiol. 2016;7:1715.
Li X, Liu L, Ji J, Chen Q, Hua X, Jiang Y, Feng Y, Yu Y. Tigecycline resistance in Acinetobacter baumannii mediated by frameshift mutation in plsC, encoding 1-acyl-sn-glycerol-3-phosphate acyltransferase. Eur J Clin Microbiol Infect Dis. 2015;34(3):625–31.
Harding CR, Schroeder GN, Collins JW, Frankel G. Use of Galleria mellonella as a model organism to study legionella pneumophila infection. J Vis Exp. 2013;81:e50964.
Gottig S, Riedel-Christ S, Saleh A, Kempf VA, Hamprecht A. Impact of blaNDM-1 on fitness and pathogenicity of Escherichia coli and Klebsiella pneumoniae. Int J Antimicrob Agents. 2016;47(6):430–5.
This work was supported by the National Natural Science Foundation of China (81672067), the Natural Science Foundation of Zhejiang province, China (LY17H190004) and the Key Research and Development Programme of Zhejiang (2015C03046). These funding had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
It was a retrospective study and which was approved by the ethnic committee of 1st affiliated hospital of Zhejiang University School of medicine and the informed consents were waived.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Figure S1. The density graph of LD50 of 56 hypervirulent Klebsiella pneumoniae isolates. (PDF 481 kb)
Additional file 2:
Table S1. PCR primers used for studying Serotype and Virulence-Associated Genes of Klebsiella pneumoniae isolates. (PDF 481 kb) (DOCX 27 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Shi, Q., Lan, P., Huang, D. et al. Diversity of virulence level phenotype of hypervirulent Klebsiella pneumoniae from different sequence type lineage. BMC Microbiol 18, 94 (2018). https://doi.org/10.1186/s12866-018-1236-2
- Hypervirulent K. pneumoniae
- Galleria. mellonella infection model
- 50% lethal dose
- Multilocus sequence typing
- Virulence level