Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2021.711577MicrobiologyOriginal ResearchWhole Genome Sequencing of Pediatric Klebsiella pneumoniae Strains Reveals Important Insights Into Their Virulence-Associated TraitsFlores-ValdezMauricio12†AresMiguel A.13†Rosales-ReyesRoberto4TorresJavier1GirónJorge A.5WeimerBart C.6Mendez-TenorioAlfonso2*De la CruzMiguel A.1*1Unidad de Investigación Médica en Enfermedades Infecciosas y Parasitarias, Hospital de Pediatría, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social, Mexico City, Mexico2Laboratorio de Biotecnología y Bioinformática Genómica, Escuela Nacional De Ciencias Biológicas (ENCB), Instituto Politécnico Nacional, Mexico City, Mexico3Departamento de Microbiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Mexico City, Mexico4Unidad de Medicina Experimental, Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico5Centro de Detección Biomolecular, Benemérita Universidad Autónoma de Puebla, Puebla, Mexico6Department of Population Health and Reproduction, School of Veterinary Medicine, 100K Pathogen Genome Project, University of California, Davis, Davis, CA, United States
Edited by: Leonard Peruski, Centers for Disease Control and Prevention (CDC), United States
Reviewed by: Jason Sahl, Northern Arizona University, United States; Claire Hennequin, Université Clermont Auvergne, France; Gabriela Jorge da Silva, University of Coimbra, Portugal
*Correspondence: Miguel A. De la Cruz, miguel_angel_81@live.comAlfonso Mendez-Tenorio, amendezt@ipn.mx
†These authors have contributed equally to this work
This article was submitted to Infectious Diseases, a section of the journal Frontiers in Microbiology
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Klebsiella pneumoniae is recognized as a common cause of nosocomial infections and outbreaks causing pneumonia, septicemia, and urinary tract infections. This opportunistic bacterium shows an increasing acquisition of antibiotic-resistance genes, which complicates treatment of infections. Hence, fast reliable strain typing methods are paramount for the study of this opportunistic pathogen’s multi-drug resistance genetic profiles. In this study, thirty-eight strains of K. pneumoniae isolated from the blood of pediatric patients were characterized by whole-genome sequencing and genomic clustering methods. Genes encoding β-lactamase were found in all the bacterial isolates, among which the blaSHV variant was the most prevalent (53%). Moreover, genes encoding virulence factors such as fimbriae, capsule, outer membrane proteins, T4SS and siderophores were investigated. Additionally, a multi-locus sequence typing (MLST) analysis revealed 24 distinct sequence types identified within the isolates, among which the most frequently represented were ST76 (16%) and ST70 (11%). Based on LPS structure, serotypes O1 and O3 were the most prevalent, accounting for approximately 63% of all infections. The virulence capsular types K10, K136, and K2 were present in 16, 13, and 8% of the isolates, respectively. Phylogenomic analysis based on virtual genome fingerprints correlated with the MLST data. The phylogenomic reconstruction also denoted association between strains with a higher abundance of virulence genes and virulent serotypes compared to strains that do not possess these traits. This study highlights the value of whole-genomic sequencing in the surveillance of virulence attributes among clinical K. pneumoniae strains.
Klebsiella pneumoniaewhole-genome sequenceMexicoblood samplesphylogenomicFIS/IMSS/PROT/G16/1617Instituto Mexicano del Seguro Social10.13039/501100004881Introduction
Klebsiella pneumoniae is a ubiquitous Gram-negative encapsulated bacterium that is commonly found in the mucosal surfaces of mammals and the environment. In humans, these bacteria can colonize the gastrointestinal and nasopharynx tracts, from where they gain entry to blood circulation and other tissues causing infection (Piperaki et al., 2017). K. pneumoniae is well-known as one of the major causes of community- and hospital-acquired Gram-negative bloodstream infections (Meatherall et al., 2009). In immunocompromised patients it can cause biliary infection, peritonitis, meningitis, septicemia, and liver abscesses (Ares et al., 2013). Also, K. pneumoniae nosocomial isolates often display high rates of antimicrobial multi-drug resistance (Chong et al., 2018). This opportunistic microorganism causes infection through several virulence factors such as surface antigens, fimbriae, capsule, outer membrane proteins, and siderophores; all of which allow the bacteria to enter and multiply within the host (Podschun and Ullmann, 1998; Ares et al., 2016, 2019). The genus Klebsiella expresses two types of surface antigens: a lipopolysaccharide (O antigen) and a capsular polysaccharide (K antigen). The K. pneumoniae capsule is a complex structure that measures up to 160 nm wide (Clements et al., 2008). The number of genes that code for the capsule of K. pneumoniae vary from 16 to 20 genes according to the capsular serotype. The capsular polysaccharide synthesis (cps) locus comprises genes that are organized into three consecutive transcriptional units that give rise to a large number of proteins involved in biosynthesis and export of capsular components. The capsule is the most thoroughly studied virulence factor of K. pneumoniae. The CPS confers resistance against the bactericidal activity of antimicrobial peptides, complement, and phagocytes (Paczosa and Mecsas, 2016). In addition, the capsule also protects the bacterium from opsonization and suppresses the early inflammatory response (Ko, 2017; Ares et al., 2019). The LPS is essential for the bacterial cell envelope structure, protects against antibiotics, and promotes adhesion. The LPS consists of three domains: the lipid A, an oligosaccharide core and the O antigen. The polysaccharide chains of the O antigen facilitate the initial adhesion process and confer resistance to serum bactericidal activity (Podschun and Ullmann, 1998; Stojkovic et al., 2017). Thus, both capsule and LPS antigens significantly contribute to the pathogenicity of the bacteria (Evrard et al., 2010). Seventy-eight K antigens and eight O antigens have been described in K. pneumoniae strains. The structural variability of these antigens has been used to classify strains into several serotypes, which in turn have elucidated differences between the degree of virulence among different strains (Follador et al., 2016). Strains carrying the K1 and K2 capsular serotypes, which mainly cause liver abscess, are known to be hypermucoviscous (Chung et al., 2012) and the serotype O1 are the most prevalent in clinical infections (Hsieh et al., 2014). Host cell adherence and biofilm formation of K. pneumoniae isolates are mediated by type 1, type 3, and ECP fimbriae whose major fimbrial subunits are encoded by the fimA, mrkA, and ecpA genes, respectively (Alcántar-Curiel et al., 2013). OmpA is one of the major outer membrane proteins of Gram-negative bacteria, and it is highly conserved among the Enterobacteriaceae. K. pneumoniae’s OmpA contributes to the attenuation of the airway epithelial cell-mediated inflammatory response and the absence of OmpA makes the microorganism more susceptible to antimicrobial peptides (Li et al., 2014). Additionally, K. pneumoniae produces two major porins OmpK35 and OmpK36, which are homologs of OmpF and OmpC of Escherichia coli, respectively, functioning as channels for export and import of hydrophilic molecules (e.g., nutrients and cephalosporins/carbapenems). The loss of outer membrane porins in K. pneumoniae leads to increased resistance to cephalosporins/carbapenems and reduced virulence (March et al., 2013; Li et al., 2014). The production of siderophores confers to bacteria the ability to sequester and uptake iron from infected tissues. K. pneumoniae produces different siderophores among which the most relevant are enterobactin, yersiniabactin, aerobactin, and salmochelin. In this study, we characterized the most relevant phenotypic and genomic virulence traits occurring among a set of 38 clinical K. pneumoniae strains isolated from a pediatric population in Mexico City. The outcome of this study highlights the importance whole genome sequence analysis in the surveillance of virulence attributes of clinical K. pneumoniae strains.
Materials and MethodsClinical Isolates
Klebsiella pneumoniae strains were isolated (September 2010 to August 2013) from individuals diagnosed with sepsis from the Laboratory of Clinical Microbiology, Pediatric Hospital, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social (IMSS), Mexico City, Mexico (Ares et al., 2013). All pediatric patients included in this study were hospitalized for more than 48 h. A clinical examination was conducted by the medical staff in order to rule out sepsis acquired outside the hospital. Microbiological methods were performed as described (Ares et al., 2013). To isolate and identify the microbial species, positive blood cultures were selected. An aliquot of 50 μL was inoculated onto MacConkey agar plates and incubated in aerobic conditions for 24–48 h at 37°C. After incubation, colonies grown on agar plates were suspended in 3 mL of sterile saline solution (0.45% NaCl, pH 7.0) to reach a McFarland turbidity of 0.5–0.63. Bacterial identification was carried out using the VITEK-2 system (bioMérieux) and genomic DNA was extracted using the DNeasy Mini Kit (Qiagen, Hilden, Germany). The clinical data are shown in Table 1.
Clinical data of K. pneumoniae strains.
Strain
Clinical service
Age
Gender
Date
Diagnosis
2k
NICU
Newborn
M
September-10
Sepsis
4k
PREESCHOOL CHILDREN
3 years
M
September-10
Fibroblastoma
5k
CHILDREN/ADOLESCENT
7 years
M
October-10
Sepsis
6k
PITU
10 years
F
November-10
Sepsis
7k
INFANTS
10 months
M
November-10
Sepsis
8k
NICU
Newborn
M
Dic-10
Sepsis
9k
PREESCHOOL CHILDREN
2 years
F
February-11
Sepsis
10k
NICU
Newborn
M
March-11
Sepsis
11k
CHILDREN/ADOLESCENT
8 years
M
June-11
Osteosarcoma
12k
PREESCHOOL CHILDREN
3 years
M
May-11
Sepsis
13k
INFANTS
3 months
M
May-11
Sepsis
15k
NICU
Newborn
F
July-11
Sepsis
16k
NICU
Newborn
F
September-11
Sepsis
17k
PREESCHOOL CHILDREN
4 years
M
October-11
Pneumonia
18k
INFANTS
6 months
M
November-11
Sepsis
19k
INFANTS
5 months
M
December-11
Sepsis
20k
CHILDREN/ADOLESCENT
8 years
F
January-12
Sepsis
21k
INFANTS
4 months
M
January-12
Sepsis
22k
NICU
Newborn
F
February-12
Sepsis
23k
NICU
Newborn
F
February-12
Sepsis
24k
INFANTS
5 months
M
May-12
Sepsis
25k
NICU
Newborn
F
June-12
Sepsis
27k
NICU
Newborn
F
July-12
Sepsis
28k
NICU
Newborn
M
August-12
Sepsis
30k
NICU
Newborn
M
September-12
Sepsis
31k
NICU
Newborn
M
September-12
Sepsis
34k
NICU
Newborn
F
October-12
Sepsis
35k
NICU
Newborn
F
October-12
Sepsis
36k
PREESCHOOL CHILDREN
5 years
F
November-12
Sepsis
37k
PREESCHOOL CHILDREN
4 years
M
November-12
Sepsis
40k
INFANTS
3 months
M
November-12
Sepsis
42k
PREESCHOOL CHILDREN
2 years
M
Dicr-12
Sepsis
43k
NICU
Newborn
M
Dic-12
Sepsis
45k
NICU
Newborn
F
January-13
Sepsis
46k
CHILDREN/ADOLESCENT
7 years
M
February-13
Sepsis
48k
NICU
Newborn
F
May-13
Sepsis
49k
NICU
Newborn
M
July-13
Sepsis
50k
NICU
Newborn
M
August-13
Sepsis
Antimicrobial Susceptibility Testing
Antimicrobial susceptibility testing for K. pneumoniae strains was carried out using Gram-negative bacilli cards (VITEK-2, AST- GN04) as was previously reported (Ares et al., 2013).
Whole-Genome Sequencing
Genomic DNA was obtained from 50 random strains of the clinical K. pneumoniae strain collection and subjected to whole genome sequencing with the HiSeq 2000 platform (Illumina) producing 150 bp paired end reads at the laboratory of Dr. Bart Weimer, University of California, Davis. De novo genome assembly was accomplished using SPAdes v.3.13.0 (Bankevich et al., 2012) using a range of k-mer sizes of 21, 33, 55, 77-mers for the assembly. The assembled draft genomes were annotated by Prokka v.1.13 (Seemann, 2014). Finally, 38 K. pneumoniae genomes showed good data quality and were deposited at NCBI. Five Klebsiella quasipneumoniae (29k, 33k, 39k, 41k, and 44k) and three Klebsiella variicola (3k, 26k, and 32k) strains were added to the bioproject. In addition, strain 1k, which is the genome of wild-type K. pneumoniae used in our group, was also included (Ares et al., 2016).
Data Access
Sequence data were submitted to the NCBI’s Gen Bank under BioProject Accession Number PRJNA605025. Each of the strains is stored with the Accession Numbers SAMN14008666-SAMN14008712 and the genomic sequences can be accessed with the Accession Numbers JAAGZT000000000-JAAHBN000000000.
Phylogenomic Analysis
Genome sequences of both K. pneumoniae isolates and NTUH K-2044 reference strain (Wu et al., 2009) were provided in FASTA format to calculate a phylogenomic tree based on Virtual Genome Fingerprints (VGF) using the VAMPhyRE software1. This method consists of carrying out a virtual hybridization analysis in which the recognition sites of a collection of 15,264 13-mer probe sequences previously designed to maximize their diversity, are searched in both strands of each genome and allowing a single-base mismatch during hybridization. The collection of hybridization sites obtained by this technique is known as the Virtual Genomic Fingerprint (VGF) and consists of a sample of sites regularly distributed throughout the genome. The VGFs are compared pairwise to estimate the fraction of shared homologous sites. A site comparison optimization analysis is previously performed, including a symmetric extension to both sites of genome alignment at each hybridization site and applying threshold match value that guarantees that only those sites not shared by chance are counted. From the fraction of homologous sites shared between each pair of genomes, a measure of their genomic distance, equivalent to the number of nucleotide substitutions between each pair, is estimated according to a method described by Nei and Li (1979). The distance is used in the MEGAX program to calculate phylogenetic trees using the Neighbor-Joining method (Kumar et al., 2018). Phylogenomic trees were edited with iTol (Interactive Tree of Life). Likewise, all assembled genomes in the dataset study were annotated using the Prokka prokaryotic annotation pipeline (Seemann, 2014) and the pan-genome was calculated with Roary: The Pan Genome Pipeline (Page et al., 2015) negative bacilli cards (VITEK-2, AST- GN04), as was previously reported (Ares et al., 2013).
Multilocus Sequence Typing
Multilocus sequence typing (MLST) analysis was performed using seven housekeeping genes (gapA, infB, mdh, pgi, phoE, rpoB, and tonB) (Diancourt et al., 2005). The allele sequences and sequence types (STs) were determined by the BIGS database of the Institute Pasteur2 (Jolley and Maiden, 2010). Only exact matches of the allelic profile were reported.
Antibiotic Resistance
The bacterial resistome was obtained by assessing the occurrence of antimicrobial resistance genes and plasmid replicons using both ResFinder 3.1 (Zankari et al., 2012) and PlasmidFinder 2.0 (Carattoli et al., 2014) from the Center for Genomic Epidemiology3. The analysis was performed with an identity threshold of 98% and a minimum match length of 80%. The families of the antibiotics assessed were fosfomycin, quinolones, phenicoles, aminoglycosides, sulphonamides, trimethoprim, tetracycline, and β-lactam.
Virulence Factors
Known virulence genes were searched on the K. pneumoniae genomes obtained above with BLASTN (Altschul et al., 1990; Camacho et al., 2009). For type 1 and type 3 fimbriae the clusters fim (fimB, fimE, fimA, fimI, fimC, fimD, fimF, fimG, fimH, and fimK) and mrk (mrkA, mrkB, mrkC, mrkD, and mrkF) were used, respectively. For ECP, the ecpA, ecpB, ecpC, ecpD, and ecpE gene cluster was analyzed. Likewise, the genes that code for outer membrane proteins (ompA, ompK35, and ompK36), the cluster virB (virB1, virB2, virB3-B4, virB5, virB6, virB7, virB8, virB9, virB10, and virB11) of the T4SS, and the hypermucoviscous associated genes magA and rmpA were screened. Additionally, fes, iroN, fyuA, and iutA genes, which correspond to enterobactin, salmochelin, yersiniabactin, and aerobactin siderophores respectively, were also searched. All the genes used for querying were obtained from the NTUH K-2044 reference strain. All analyses were carried out with an identity threshold of 96% and a minimum match length of 80%.
Serotyping (K and O Antigens)
The genes wzi and wzc were searched in the draft genomes of K. pneumoniae using BLASTN (Altschul et al., 1990; Camacho et al., 2009), a database was assembled utilizing the corresponding locus of K serotypes published by Wyres et al. (2016). In the same manner the locus wb was searched and compared with the data published by Fang et al. (2016) to identify O serotypes. All analyses were carried out with an identity threshold of 96% and a minimum length of 80%. Additionally, the software Kaptive4 (Wick et al., 2018) was used to determine the K and O serotypes. Results that matched in both methodologies were reported as positive for the searched serotype.
Glucuronic Acid Analysis
Capsular polysaccharides were extracted and quantified as previously described (Ares et al., 2016). The glucuronic acid concentration in each bacterium was expressed in micrograms/109 CFU.
Phagocytosis of Bacteria by THP-1 Macrophages
The phagocytosis assay was performed using THP-1 (ATCC TIB-202) human monocytes as previously described (Ares et al., 2016). Phagocytosis resistance was expressed in intracellular colony-forming-units per mL (CFUs/mL).
Antimicrobial Activity of Polymyxin B
All experiments were performed with an initial inoculum of ∼106 CFU/mL in 20 mL of LB medium in 50 mL pyrogen-free and sterile polypropylene tubes. Polymyxin B activity was examined at 4 mg/mL. Serial strains in 20 mL were collected at 4 h for viable counting on LB agar plates after 24 h of incubation at 37°C as described Nang et al., 2019.
ResultsWhole-Genome Sequencing of K. pneumoniae Strains
Genomic sequences of 38 bacterial strains positively identified as K. pneumoniae were obtained and they were named with a number followed by a letter “k” (e.g., 1k; Table 1). The draft genomes had an average assembly length of 5.55 Mb ± 0.14; a GC content of 57.23% ± 0.16%; N50 of 169,445 ± 41,243; and CDS of 5,135 ± 144 (assembly data and accession are shown in Table 2). The K. pneumoniae isolates shared a core genome consisting of 1,192 genes (present in ≥99% isolates and with 95% identity). The genomes analysis showed that with exception of strain 46k, all K. pneumoniae presented at least one plasmid sequence (Table 3).
Virtual Genomic Fingerprinting (VGF) Correlates With MLST
Twenty-four different sequence types (ST) were identified among the K. pneumoniae isolates in this study. The most prevalent STs were ST-76 in six strains followed by ST-70 with four occurrences (Table 4). Interestingly, two K. pneumoniae strains, 46k and 43k, presented two new sequence types ST-4872 and ST-4873, respectively. Our analysis showed that ST-76 and ST-70 are the most prevalent ST in the strains collection and that new STs such as ST-4872 and ST-4873, appeared in the bacterial population of K. pneumoniae. The phylogenetic analysis was carried out based on the virtual genomic fingerprint technique (VAMPhyRE) using 13 k-mer probes and the Nearest-Neighbor method in MEGAX with default settings. The final phylogenomic representation (Figure 1) denotes consistent correlation with MLST findings.
Capsular polysaccharide synthesis (CPS), LPS, and MLST in K. pneumoniae strains.
Strain
LPS
CPS
MLST
2k
O1
K28
ST-4839
4k
O2
K57
ST-471
5k
O3
K123
ST-628
6k
O1
K2
ST-25
7k
O1
K2
ST-25
8k
O1
K136
ST-70
9k
O1
K136
ST-70
10k
O1
K62
ST-348
11k
O1
K136
ST-70
12k
OL101
K49
ST-1401
13k
O4
K58
ST-405
15k
O4
K17
ST-322
16k
O1
K2
ST-380
17k
O4
K17
ST-322
18k
O2
K33
ST-441
19k
O2
K23
ST-280
20k
O4
K58
ST-405
21k
O1
K33
ST-337
22k
O2
K62
ST-45
23k
O2
K24
ST-45
24k
O1
K17
ST-870
25k
O1
K136
ST-70
27k
O2
K136
ST-399
28k
O1
K28
ST-20
30k
O3
K10
ST-76
31k
O1
K28
ST-20
34k
O3
K10
ST-76
35k
O3
K10
ST-76
36k
O1
K30
ST-29
37k
O1
K30
ST-29
40k
O2
K27
ST-1440
42k
O3
K10
ST-76
43k
O3
K24
ST-4873*
45k
O4
K15
ST-37
46k
O3
K39
ST-4872*
48k
O3
K123
ST-834
49k
O3
K10
ST-76
50k
O3
K10
ST-76
*means new MLST found in this study.
Phylogenomic reconstruction of Klebsiella pneumoniae strains by VGF. NTUH-K2044 strain was added as reference. Yellow and orange colors denote the clustering of strains by MLST.
Resistome of K. pneumoniae Strains
The results for the search of resistance genes in the K. pneumoniae strains studied is presented in Table 5. The prevalence of genes encoding resistance to fosfomycin, quinolones, phenicols, aminoglycosides, sulphonamides, trimethoprim, tetracycline was 100, 97, 5, 79, 68, 39, and 45%, respectively. Moreover, β-lactamase genes (Table 5) such as blaSHV, blaOXA, blaTEM, and blaCTX–M were present in 53, 53, 5, and 26%, of the strains, respectively. No blaOKP, blaLEN genes were found in this bacterial collection. All the strains carried at least two different types of resistance genes. Antimicrobial susceptibility testing of K. pneumoniae strains (Table 6) showed a similar profile correspondent to the resistance genes encoded in the bacterial genome.
Antibiotic resistance genes in K. pneumoniae strains.
Strain
Fosfomycin
Quinolones
Phenicols
Aminoglycosides
Sulphonamides
Trimethoprim
Tetracycline
b-lactam
2k
fosA
oqxA, qnrB
aac(3), aadA1, strB
sul1
dfrA
tet(A)
blaSHV–28
4k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
5k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaSHV–1, blaTEM–1B
6k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
7k
fosA
oqxA
aac(3)
blaCTX–M–15, blaTEM–1B
8k
fosA
oqxA
aac(6′)-33, aadA1
sul1
blaOXA–2
9k
fosA
oqxA, qnrB
catB3
aac(3),strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
10k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
11k
fosA
oqxA
aac(6′), aadA1
sul1
blaOXA–2
12k
fosA
oqxA
blaSHV–12
13k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
15k
fosA
oqxA
blaSHV–11
16k
fosA
oqxA
blaSHV–1
17k
fosA
oqxA
blaSHV–11
18k
fosA
blaSHV–38
19k
fosA
oqxA,
catB3
strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
20k
fosA
oqxA, qnrB
catB3
aac(3), strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaOXA–1, blaTEM–1B
21k
fosA
oqxA
strA, strB
sul2
dfrA
tet(A)
blaSHV–11
22k
fosA
oqxA, qnrB
strA, strB
sul2
dfrA
tet(A)
blaCTX–M–15, blaSHV–1, blaTEM–1B
23k
fosA
oqxA
tet(A)
blaSHV–1
24k
fosA
oqxA
strA, strB
tet(A)
blaSHV–1
25k
fosA
oqxA
aac(6′), aadA1
sul1
blaOXA–2
27k
fosA
oqxA
catA1
aac(3), aadA1
sul1
dfrA
blaSHV–28, blaTEM–1B
28k
fosA
oqxA
strA, strB
blaSHV–11
30k
fosA
oqxA
catA1
aac(6′), aadA1
sul1
blaOXA–2, blaTEM–1B
31k
fosA
oqxA
strA, strB
blaSHV–11
34k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
blaOXA–2, blaSHV–129, blaTEM–1B
35k
fosA
oqxA
catA1
aac(6′), aadA1
sul1
blaOXA–2, blaTEM–1B
36k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
dfrA
tet(A)
blaOXA–2, blaTEM–1B
37k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
dfrA
tet(A)
blaOXA–2, blaTEM–1B
40k
fosA
oqxA
blaSHV–11
42k
fosA
oqxA
catA1
aac(6′), aadA1, aadA1
sul1
blaOXA–2, blaTEM–1B
43k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
dfrA
tet(A)
blaOXA–2, blaTEM–1B
45k
fosA
oqxA
aadA1
sul1
blaSHV–11
46k
fosA
oqxA
blaSHV–1
48k
fosA
oqxA
catA1
strA, strB
sul2
tet(A)
blaSHV–11, blaTEM–1B
49k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
blaOXA–2, blaSHV–129, blaTEM–1B
50k
fosA
oqxA
catA1
aac(6′), aadA1, aadA5
sul1
blaOXA–2, blaTEM–1B
Positives
38
37
19
30
26
15
17
38
Prevalence
1.00
0.97
0.50
0.79
0.68
0.39
0.45
1.00
MIC (μg/mL) of the 38 clinical isolates of K. pneumoniae according to CLSI.
Strain
FM
CF
MF
CM
AK
GM
TM
T/S
TC
2K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
>4.0R
>16R
4K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
5K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
6K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
7K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
<2.0
<4.0
8K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
<2.0
<4.0
9K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
10K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
11K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
<2.0
<4.0
12K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
13K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
15K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
16K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
17K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
18K
>256R
<0.25
<0.25
<8.0
<16
<4.0
<4.0
<2.0
<4.0
19K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
20K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
21K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
>4.0R
>16R
22K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
>4.0R
>16R
23K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
>16R
24K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
>16R
<2.0
>16R
25K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
<4.0
<2.0
<4.0
27K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
<4.0
28K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
<4.0
<2.0
<4.0
30K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
31K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
<4.0
<2.0
<4.0
34K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
35K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
36K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
37K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
40K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
42K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
43K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
>16R
>4.0R
>16R
45K
>256R
>1.0R
>1.0R
<8.0
>64R
>16R
<4.0
<2.0
<4.0
46K
>256R
>1.0R
>1.0R
<8.0
<16
<4.0
<4.0
<2.0
<4.0
48K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
>16R
49K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
50K
>256R
>1.0R
>1.0R
>32R
>64R
>16R
<4.0
<2.0
<4.0
MIC, minimum inhibitory concentration; CLSI, Clinical and Laboratory Standards Institute; FM, fosfomycin. Quinolones: CF, ciprofloxacin; MF, moxifloxacin. Phenicols: CM, chloramphenicol. Aminoglycosides: AK, amikacin; GM, gentamicin; TM, tobramycin. Trimethoprim/Sulphonamides: T/S, trimethoprim/sulfamethoxazole. Tetracyclines: TC, tetracycline. R, resistant.The K. pneumoniae Genomes Revealed a Myriad of Virulence Factors
The prevalence of virulence genes in the genomes of K. pneumoniae strains is shown in Figure 2. For fimbrial genes, the prevalence was very high showing that, with exception of strain 27k that lacks the fim cluster, all the strains possess the fim (Type 1 fimbriae), mrk (Type 3 fimbriae) and ecp (E. coli common pilus) genes. In terms of genes that code for outer membrane proteins, ompK35 and ompA were present in all K. pneumoniae genomes, while ompK36 gene presented in 63% of the strains. Interestingly, seven (18%) K. pneumoniae strains were positive for vir genes, which code for a putative T4SS. Regarding siderophores, fes (enterobactin) and fyuA (yersiniabactin) genes were present in 100 and 18% of the K. pneumoniae strains, respectively. Of note, only strain 16k carried genes that code for the four siderophores: enterobactin, yersiniabactin, aerobactin, and salmochelin. In addition, this strain was the only one to possess the hypermucoviscous associated gene rmpA. None of the strains were positive for magA gene. Furthermore, a correlation is depicted between the presence of virulence factors and the clustering by the VGF in Figure 2. Strains 16k and the NTUH K-2044, which possess the higher number of the virulence factors, were grouped into a common clade. Strains 6k, 22k, 23k, and 9k that contain the presence of fyuA and vir genes, formed an independent clade. Interestingly, strains 13k and 20k that belong to the O4:K58 serotype were allocated into a different clade although they shared the fyuA and vir genes.
Virulome of K. pneumoniae strains. Color filled box mean presence of virulence factor. NTUH-K2044 was added as reference strain.
The K. pneumoniae K2 Serotype Is More Resistant to Phagocytosis and Polymyxin B
The analysis of the O groups among the strains showed that the most prevalent serogroup was O1 (39%), followed by O3 (26%) (Table 4). Seventeen different K-serotypes were identified in the all strains, being serotype K10 the most prevalent (16%), followed by K136 (13%) (Table 4). Three strains were positive for serotype K2 (6k, 7k, and 16k), which is associated to hypervirulence phenotype due to capsule hyperproduction (Yeh et al., 2007; Yu et al., 2008). We analyzed the capsule production and phagocytosis resistance in K. pneumoniae K2 serotype strains and compared them to the K10 and K136 serotypes, which were the most prevalent (Figure 3). In addition, polymyxin B susceptibility assays were performed to evaluate the protective capacity of the capsule against this antimicrobial cationic peptide (Campos et al., 2004). As controls, we used the wild-type K. pneumoniae 123/01 strain (K39 serotype) and an Δhns-derivative mutant, which show low and high level of capsule production, respectively (Ares et al., 2016). Indeed, the K2 serotype strains produced more capsule than K10/K136 serotypes to similar levels as the Δhns mutant, which is a hypercapsulated K. pneumoniae. Therefore, K. pneumoniae K2 serotype strains were less phagocytized by THP-1 macrophages than K10 and K136 serotypes strains. In addition, the action of polymyxin B was lower in the three K2 serotype K. pneumoniae.
K2 serotype K. pneumoniae is highly resistant to phagocytosis and polymyxin B. (A) Capsule quantification by determination of the glucuronic acid concentration from capsular polysaccharides extracted from 0.5 ml bacterial cultures. (B) Comparison of phagocytic uptake of K. pneumoniae grown under the indicated different conditions by THP-1 macrophages. (C) Polymyxin B activity (4 mg/ml) against K. pneumoniae strains. An Δhns mutant was evaluated as control of phagocytosis and polymyxin B. ns, not significant; statistically significant with respect to bacteria grown in LB medium: ***p < 0.001.
Discussion
Gram-negative bacilli infections are the most serious and potentially life-threatening causes of infectious disease in hospitalized patients (Joram et al., 2012; Ares et al., 2013). K. pneumoniae is an opportunistic pathogen associated with both community-acquired and nosocomial infections, causing pneumonia, abscess, bacteremia, and urinary tract infections (Gassama-Sow et al., 2010). Due to its ability to rapidly acquire antimicrobial resistance, K. pneumoniae has become a global concern calling for actions to prevent the spread of multidrug resistant bacteria (Jansen et al., 2018). Over the last few years, whole-genome sequencing (WGS) has become more accessible and affordable, resulting in its increased use in many fields, including clinical microbiology (Mathers et al., 2015). One of the main benefits of using WGS is to characterize the genomic content of clinically relevant bacteria and relate it to virulence-associated phenotypes to allow the understanding of their transmission within the hospital and to provide accurate diagnostics and to apply prompt therapies. The present study aimed to characterize K. pneumoniae isolates from pediatric patients hospitalized between 2010 and 2013 in Mexico City. This is the first study in Mexico that uses WGS to deeply characterize this opportunistic pathogen. In this work, 24 STs were identified in the bacterial genomes analyzed, showing in general, a heterogeneous distribution among which ST-76 (16%) and ST-70 (11%) were the most prevalent STs. Of note, the ST-29 represented by strains 36k and 37k, has previously been associated with hypervirulence (Shon et al., 2013; Lev et al., 2018). This same ST-29 was previously reported in isolates from community-acquired infections in the State of Guerrero, México, presenting a hypermucoviscosity phenotype (Garza-Ramos et al., 2018).
Regarding the prevalence of antibiotic resistance, K. pneumoniae genomes showed high prevalence (97%) for oqxA or qnrB genes, which conferred resistance to quinolones. Only one strain (18k) was sensitive to ciprofloxacin and moxifloxacin, which corresponded to the lack of oqxA/qnrB genes in the bacterial genome. This high resistance to this class of antibiotics has been reported by other researchers in both Mexico and worldwide (Rodríguez-Zulueta et al., 2013; Osei Sekyere and Amoako, 2017; Salah et al., 2019; Geetha et al., 2020). Sulfonamide resistance genes such as sul1 or sul2 were highly prevalent (68%) in the K. pneumoniae strains analyzed in this study, which is in agree with other reports (Khamesipour and Tajbakhsh, 2016; Taitt et al., 2017). Genes that confer tetracycline [tet(A), tet(B), tet(C), tet(D), or tet(E)] were present in 45% of the strains, which is low in comparison to other studies (Bokaeian et al., 2015; Lev et al., 2018). All strains presented at least one type of gene coding for β-lactamases and the SHV-type was the most prevalent in the sequenced genomes (53%). However, for the blaTEM-type, only blaTEM–1B gene was identified in 20 strains (53%), which has been reported as the most prevalent β-lactamases genes in K. pneumoniae (Chang et al., 2001; Feizabadi et al., 2010; Shaikh et al., 2015). The β-lactamases genes blaOXA–1 and blaOXA–2 were found in 20 strains (53%), while for the blaCTX–M–15 gene was found in 10 bacterial genomes analyzed (26%), which has been reported by different authors as common in nosocomial outbreaks (Zhou et al., 2015; Bevan et al., 2017; Gautam et al., 2019). β-lactamases genes blaOKP and blaLEN were not found in the strains.
Seventeen different types of CPS were found among the clinical K. pneumoniae collection. Three strains (6k, 7k, and 16k) were identified as belonging to the K2 capsular serotype, which is the main antigen associated with hypervirulence in K. pneumoniae (Brisse et al., 2013; Lin et al., 2014; Lev et al., 2018). These three strains produced more capsule and consequently they were more resistant to both macrophage-mediated phagocytosis and polymyxin, than the K10 and K136 capsular types, which were the most prevalent serotypes in our study. Our results are consistent with previous studies confirming that the K2 serotype confer K. pneumoniae higher resistance to both macrophage-mediated phagocytosis and polymyxin B due to increased capsule production (Yu et al., 2015). In the case of the LPS O-antigen, five different serotypes were found. In agreement with other studies (Shankar-Sinha et al., 2004; Hsieh et al., 2014; Follador et al., 2016), the O1 antigen, which confers greater resistance to the serum’s bactericidal activity, was the most prevalent O1 serogroup (39%) in the K. pneumoniae strains isolated from blood. Whilst O2 serogroup is more sensitive to human serum killing than O1 (Pennini et al., 2017), there are not reports about the O3 antigen immunogenicity, which was the second serogroup most prevalent (24%) in the K. pneumoniae strains collection.
In terms of adhesive structures, type 3 fimbria is the most important organelle of K. pneumoniae, playing an important role in adherence to both abiotic and biotic surfaces (Langstraat et al., 2001; Di Martino et al., 2003; Struve et al., 2009; Schroll et al., 2010; Ares et al., 2016). Virulence gene analysis revealed that the mrk cluster was present in all the strains, supporting the notion that type 3 fimbriae are a hallmark of K. pneumoniae virulence. Porins have been shown to be important for bacterial survival due to their role in the exchange of substances, including nutrients and toxic metabolites (Tsai et al., 2011). All the strains presented both the ompA and ompK35 genes. In contrast, only 25 (66%) K. pneumoniae strains were positive for ompK36 gene. The high prevalence of ompK35 versus ompK36, supports the structural role of OmpK35 as a porin related to the outer membrane integrity, whilst OmpK36 provides increased antibiotic resistance with no significant loss of fitness in the host (Fajardo-Lubián et al., 2019). Interestingly, K. pneumoniae strains lacking ompK36 presented over 70% prevalence of blaOXA carbapenemase gene, which support the notion that the resistance to carbapenem antibiotics is related to the absence of ompK36 porin gene (Landman et al., 2009; Wong et al., 2019). The combination of siderophores secreted by K. pneumoniae during infection affect tissue localization, systemic dissemination, and host survival (Holden et al., 2016). Enterobactin genes were present in 100% of the K. pneumoniae genomes, which is understandable since the catecholate is typically encoded in the core genome of K. pneumoniae (Martin and Bachman, 2018). Genes that encodes yersiniabactin siderophore were present in seven (18%) K. pneumoniae, and two of these genomes presented the K2 capsular serotype, reflecting the virulence properties of these strains, since yersiniabactin is significantly associated with invasive infections in humans (Holt et al., 2015).
With exception of the strain 46k, at least one plasmid sequence was detected in 37 of the 38 genomes analyzed. IncFIB plasmid was present in 36 of 37 of the bacterial isolates. This conjugative plasmid has been associated with the dissemination of the blaNDM–1, blaSHV–12, blaCTXM–15, and blaOXA–1 genes in K. pneumoniae (Oliveira et al., 2020). Methods of circularization and annotation of the plasmid DNA using sequencing techniques that generate long-read sequences will offer insights into the content of those plasmids and its evolution and spread (Ramsamy et al., 2020).
Traditional phylogenetic approaches require the calculation of multiple sequence alignment of homologous sequences for finding informative characters or for the calculation of sequence differences that are used as overall measures to estimate evolutionary distances. Thus, the time required for a phylogenetic study is distributed into that which is needed for calculating the alignment and that which is needed for the phylogenetic method itself (Yang and Rannala, 2012). The VAMPhyRE method is an alignment free tool based in the calculation of Virtual Genomic Fingerprints by a virtual hybridization with predesigned probe sets. Additionally, SNPs methods usually require a collection of SNPs relative to a reference sequence which is not required by the VAMPhyRE approach. The phylogenomic reconstructions obtained by VAMPhyRE correlated to the MLST determination. These results demonstrated a strong relation between MLST classification scheme and the phylogenomic analysis in the bacterial strains. This same behavior has been reported by Muñoz-Ramírez et al. (2017) in Helicobacter pylori strains. Furthermore, the VAMPhyRE methodology was able to cluster together strains with high number of virulence factors and hypervirulent serotypes from those lacking them. Interestingly, strain 16k presented a strong phylogenetic relation with NTUH-K2044 reference strain, carrying all virulence factors analyzed in our study, including the four siderophores and the T4SS genes. Ongoing work in our laboratory aims to characterize the pathogenic features of this strain. In summary, the search of virulence factors by using whole genome sequencing and virtual genome fingerprinting is an alternative method to classify clinical K. pneumoniae strains.
Data Availability Statement
The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
MD and AM-T conceived and designed the experiments. MA and RR-R performed the experiments. MF-V, MA, JT, JG, BW, AM-T, and MD analyzed the data. MF-V, JG, and MD wrote the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Funding. This work was supported by grants FIS/IMSS/PROT/G16/1617 (to MD) from the Fondo de Investigación en Salud (FIS)-Instituto Mexicano del Seguro Social (IMSS), Mexico. MF-V was supported by a pre-doctoral fellowship from CONACYT (586925).
The genomes were sequenced using the 100K Pathogen Genome Project resources. We thank the team of curators of the Institut Pasteur MLST and whole genome MLST databases for curating the data and making them publicly available at http://bigsdb.pasteur.fr.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.711577/full#supplementary-material
ReferencesAlcántar-CurielM. D.BlackburnD.SaldañaZ.Gayosso-VázquezC.IovineN. M.De la CruzM. A. (2013). Multi-functional analysis of Klebsiella pneumoniae fimbrial types in adherence and biofilm formation.4129–138. 10.4161/viru.2297423302788AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool.215403–410. 10.1016/S0022-2836(05)80360-80362AresM. A.Fernández-VázquezJ. L.Rosales-ReyesR.Jarillo-QuijadaM. D.von BargenK.TorresJ. (2016). H-NS nucleoid protein controls virulence features of Klebsiella pneumoniae by regulating the expression of Type 3 pili and the capsule polysaccharide.6:13. 10.3389/fcimb.2016.0001326904512AresM. A.SansabasA.Rodríguez-ValverdeD.Siqueiros-CendónT.Rascón-CruzQ.Rosales-ReyesR. (2019). The interaction of Klebsiella pneumoniae with lipid rafts-associated cholesterol increases macrophage-mediated phagocytosis due to down regulation of the capsule polysaccharide.9:255. 10.3389/fcimb.2019.0025531380298AresM. ÁAlcántar-CurielM. D.Jiménez-GaliciaC.Rios-SarabiaN.PachecoS.De la CruzM. Á (2013). Antibiotic resistance of gram-negative bacilli isolated from pediatric patients with nosocomial bloodstream infections in a Mexican tertiary care hospital.59361–368. 10.1159/00036208524821320BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.19455–477. 10.1089/cmb.2012.002122506599BevanE. R.JonesA. M.HawkeyP. M. (2017). Global epidemiology of CTX-M β-lactamases: temporal and geographical shifts in genotype.722145–2155. 10.1093/jac/dkx14628541467BokaeianM.Shahraki ZahedaniS.Soltanian BajgiranM.Ansari MoghaddamA. (2015). Frequency of PER, VEB, SHV, TEM and CTX-M genes in resistant strains of Pseudomonas aeruginosa producing extended spectrum β-Lactamases.8:e13783. 10.5812/jjm.1378325789123BrisseS.PassetV.HaugaardA. B.BabosanA.Kassis-ChikhaniN.StruveC. (2013). wzi gene sequencing, a rapid method for determination of capsular type for Klebsiella strains.514073–4078. 10.1128/JCM.01924-1913CamachoC.CoulourisG.AvagyanV.MaN.PapadopoulosJ.BealerK. (2009). BLAST+: architecture and applications.10:421. 10.1186/1471-2105-10-42120003500CamposM. A.VargasM. A.RegueiroV.LlompartC. M.AlbertíS.BengoecheaJ. A. (2004). Capsule polysaccharide mediates bacterial resistance to antimicrobial peptides.727107–7114. 10.1128/iai.72.12.7107-7114.200415557634CarattoliA.ZankariE.García-FernándezA.Voldby LarsenM.LundO.VillaL. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing.583895–3903. 10.1128/AAC.02412-1424777092ChangF.-Y.SiuL. K.FungC.-P.HuangM.-H.HoM. (2001). Diversity of SHV and TEM β-Lactamases in Klebsiella pneumoniae: gene evolution in northern taiwan and two novel β-Lactamases, SHV-25 and SHV-26.452407–2413. 10.1128/AAC.45.9.2407-2413.200111502506ChongY.ShimodaS.ShimonoN. (2018). Current epidemiology, genetic evolution and clinical impact of extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae.61185–188. 10.1016/j.meegid.2018.04.00529626676ChungD. R.ParkM. H.KimS. H.KoK. S.KangC.-I.PeckK. R. (2012). Prevalence and molecular characterization of serotype K1 Klebsiella pneumoniae strains from various clinical specimen sources in 11 Asian countries.64622–625. 10.1016/j.jinf.2012.02.00722343125ClementsA.GaboriaudF.DuvalJ. F. L.FarnJ. L.JenneyA. W.LithgowT. (2008). The major surface-associated saccharides of Klebsiella pneumoniae contribute to host cell association.3:e3817. 10.1371/journal.pone.000381719043570Di MartinoP.CafferiniN.JolyB.Darfeuille-MichaudA. (2003). Klebsiella pneumoniae type 3 pili facilitate adherence and biofilm formation on abiotic surfaces.1549–16. 10.1016/S0923-2508(02)00004-9DiancourtL.PassetV.VerhoefJ.GrimontP. A. D.BrisseS. (2005). Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates.434178–4182. 10.1128/JCM.43.8.4178-4182.200516081970EvrardB.BalestrinoD.DosgilbertA.Bouya-GachancardJ.-L. J.CharbonnelN.ForestierC. (2010). Roles of capsule and lipopolysaccharide O antigen in interactions of human monocyte-derived dendritic cells and Klebsiella pneumoniae.78210–219. 10.1128/IAI.00864-869Fajardo-LubiánA.Ben ZakourN. L.AgyekumA.QiQ.IredellJ. R. (2019). Host adaptation and convergent evolution increases antibiotic resistance without loss of virulence in a major human pathogen.15:e1007218. 10.1371/journal.ppat.100721830875398FangC.-T.ShihY.-J.CheongC.-M.YiW.-C. (2016). Rapid and accurate determination of lipopolysaccharide O-Antigen types in Klebsiella pneumoniae with a novel PCR-Based O-Genotyping method.54666–675. 10.1128/JCM.02494-2415FeizabadiM. M.Mahamadi-YeganehS.MirsalehianA.MirafsharS.-M.MahboobiM.NiliF. (2010). Genetic characterization of ESBL producing strains of Klebsiella pneumoniae from Tehran hospitals.4609–615. 10.3855/jidc.105921045352FolladorR.HeinzE.WyresK. L.EllingtonM. J.KowarikM.HoltK. E. (2016). The diversity of Klebsiella pneumoniae surface polysaccharides.2:e000073. 10.1099/mgen.0.00007328348868Garza-RamosU.Barrios-CamachoH.Moreno-DomínguezS.Toribio-JiménezJ. (2018). Phenotypic and molecular characterization of Klebsiella spp. isolates causing community-acquired infections.2317–27. 10.1016/j.nmni.2018.02.00229692906Gassama-SowA.DialloM. H.WaneA. A.SeckA.Samb-BaB.SowP. S. (2010). Genetic determinants of antibiotic resistance in diarrheagenic Klebsiella pneumoniae subspecies ozaenae: an emerging enteropathogen in Senegal.50453–455. 10.1086/64989220064035GautamV.ThakurA.SharmaM.SinghA.BansalS.SharmaA. (2019). Molecular characterization of extended-spectrum β-lactamases among clinical isolates of Escherichia coli & Klebsiella pneumoniae: a multi-centric study from tertiary care hospitals in India.149208–215. 10.4103/ijmr.IJMR_172_18GeethaP. V.AishwaryaK. V. L.MariappanS.SekarU. (2020). Fluoroquinolone resistance in clinical isolates of e.12121–125.HoldenV. I.BreenP.HouleS.DozoisC. M.BachmanM. A. (2016). Klebsiella pneumoniae siderophores induce inflammation, bacterial dissemination, and HIF-1α stabilization during pneumonia.7:e01397-16. 10.1128/mBio.01397-1316HoltK. E.WertheimH.ZadoksR. N.BakerS.WhitehouseC. A.DanceD. (2015). Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health.112E3574–E3581.HsiehP.-F.WuM.-C.YangF.-L.ChenC.-T.LouT.-C.ChenY.-Y. (2014). D-galactan II is an immunodominant antigen in O1 lipopolysaccharide and affects virulence in Klebsiella pneumoniae: implication in vaccine design.5:608. 10.3389/fmicb.2014.0060825477867JansenK. U.KnirschC.AndersonA. S. (2018). The role of vaccines in preventing bacterial antimicrobial resistance.2410–19. 10.1038/nm.446529315295JolleyK. A.MaidenM. C. J. (2010). BIGSdb: scalable analysis of bacterial genome variation at the population level.11:595. 10.1186/1471-2105-11-59521143983JoramN.de Saint BlanquatL.StammD.LaunayE.Gras-Le GuenC. (2012). Healthcare-associated infection prevention in pediatric intensive care units: a review.312481–2490. 10.1007/s10096-012-1611-1610KhamesipourF.TajbakhshE. (2016). Analyzed the genotypic and phenotypic antibiotic resistance patterns of Klebsiella pneumoniae isolated from clinical samples in Iran.27, 0970–938X.KoK. S. (2017). The contribution of capsule polysaccharide genes to virulence of Klebsiella pneumoniae.8485–486. 10.1080/21505594.2016.124086227715471KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms.351547–1549. 10.1093/molbev/msy09629722887LandmanD.BratuS.QualeJ. (2009). Contribution of OmpK36 to carbapenem susceptibility in KPC-producing Klebsiella pneumoniae.58(Pt 10), 1303–1308. 10.1099/jmm.0.012575-019556371LangstraatJ.BohseM.CleggS. (2001). Type 3 fimbrial shaft (MrkA) of Klebsiella pneumoniae, but not the fimbrial adhesin (MrkD), facilitates biofilm formation.695805–5812. 10.1128/iai.69.9.5805-5812.200111500458LevA. I.AstashkinE. I.KislichkinaA. A.SolovievaE. V.KombarovaT. I.KorobovaO. V. (2018). Comparative analysis of Klebsiella pneumoniae strains isolated in 2012-2016 that differ by antibiotic resistance genes and virulence genes profiles.112142–151. 10.1080/20477724.2018.146094929708041LiB.ZhaoY.LiuC.ChenZ.ZhouD. (2014). Molecular pathogenesis of Klebsiella pneumoniae.91071–1081. 10.2217/fmb.14.4825340836LinJ.-C.KohT. H.LeeN.FungC.-P.ChangF.-Y.TsaiY.-K. (2014). Genotypes and virulence in serotype K2 Klebsiella pneumoniae from liver abscess and non-infectious carriers in Hong Kong, Singapore and Taiwan.6:21. 10.1186/1757-4749-6-2124987462MarchC.CanoV.MorantaD.LlobetE.Pérez-GutiérrezC.TomásJ. M. (2013). Role of bacterial surface structures on the interaction of Klebsiella pneumoniae with phagocytes.8:e56847. 10.1371/journal.pone.005684723457627MartinR. M.BachmanM. A. (2018). Colonization, infection, and the accessory genome of Klebsiella pneumoniae.8:4. 10.3389/fcimb.2018.0000429404282MathersA. J.StoesserN.SheppardA. E.PankhurstL.GiessA.YehA. J. (2015). Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae at a single institution: insights into endemicity from whole-genome sequencing.591656–1663. 10.1128/AAC.04292-4214MeatherallB. L.GregsonD.RossT.PitoutJ. D. D.LauplandK. B. (2009). Incidence, risk factors, and outcomes of Klebsiella pneumoniae bacteremia.122866–873. 10.1016/j.amjmed.2009.03.03419699383Muñoz-RamírezZ. Y.Mendez-TenorioA.KatoI.BravoM. M.RizzatoC.ThorellK. (2017). Whole genome sequence and phylogenetic analysis show strains from latin america have followed a unique evolution pathway.7:50. 10.3389/fcimb.2017.0005028293542NangS. C.HanM.-L.YuH. H.WangJ.TorresV. V. L.DaiC. (2019). Polymyxin resistance in Klebsiella pneumoniae: multifaceted mechanisms utilized in the presence and absence of the plasmid-encoded phosphoethanolamine transferase gene mcr-1.743190–3198. 10.1093/jac/dkz31431365098NeiM.LiW. H. (1979). Mathematical model for studying genetic variation in terms of restriction endonucleases.765269–5273. 10.1073/pnas.76.10.5269291943OliveiraÉ. M.BeltrãoE.ScavuzziA.BarrosJ. F.LopesA. (2020). High plasmid variability, and the presence of IncFIB, IncQ, IncA/C, IncHI1B, and IncL/M in clinical isolates of Klebsiella pneumoniae with bla KPC and bla NDM from patients at a public hospital in Brazil.53: e20200397. 10.1590/0037-8682-0397-202033111914Osei SekyereJ.AmoakoD. G. (2017). Genomic and phenotypic characterisation of fluoroquinolone resistance mechanisms in Enterobacteriaceae in Durban. South Africa.12:e0178888. 10.1371/journal.pone.017888828636609PaczosaM. K.MecsasJ. (2016). Klebsiella pneumoniae: going on the offense with a strong defense.80629–661. 10.1128/MMBR.00078-1527307579PageA. J.CumminsC. A.HuntM.WongV. K.ReuterS.HoldenM. T. G. (2015). Roary: rapid large-scale prokaryote pan genome analysis.313691–3693. 10.1093/bioinformatics/btv42126198102PenniniM. E.De MarcoA.PelletierM.BonnellJ.CvitkovicR.BeltramelloM. (2017). Immune stealth-driven O2 serotype prevalence and potential for therapeutic antibodies against multidrug resistant Klebsiella pneumoniae.8:1991. 10.1038/s41467-017-02223-2227PiperakiE.-T.SyrogiannopoulosG. A.TzouvelekisL. S.DaikosG. L. (2017). Klebsiella pneumoniae: virulence, biofilm and antimicrobial resistance.361002–1005. 10.1097/INF.000000000000167528914748PodschunR.UllmannU. (1998). Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors.11589–603. 10.1128/CMR.11.4.5899767057RamsamyY.MlisanaK. P.AllamM.AmoakoD. G.AbiaA. L. K.IsmailA. (2020). Genomic analysis of carbapenemase-producing extensively drug-resistant Klebsiella pneumoniae isolates reveals the horizontal spread of p18-43_01 plasmid encoding blaNDM-1 in South Africa.8:137. 10.3390/microorganisms801013731963608Rodríguez-ZuluetaP.Silva-SánchezJ.BarriosH.Reyes-MarJ.Vélez-PérezF.Arroyo-EscalanteS. (2013). First outbreak of KPC-3-producing Klebsiella pneumoniae (ST258) clinical isolates in a Mexican Medical Center.574086–4088. 10.1128/AAC.02530-2512SalahF. D.SoubeigaS. T.OuattaraA. K.SadjiA. Y.Metuor-DabireA.Obiri-YeboahD. (2019). Distribution of quinolone resistance gene (qnr) in producing and in Lomé.8:104.SchrollC.BarkenK. B.KrogfeltK. A.StruveC. (2010). Role of type 1 and type 3 fimbriae in Klebsiella pneumoniae biofilm formation.10:179. 10.1186/1471-2180-10-17920573190SeemannT. (2014). Prokka: rapid prokaryotic genome annotation.302068–2069. 10.1093/bioinformatics/btu15324642063ShaikhS.FatimaJ.ShakilS.RizviS. M.KamalM. A. (2015). Risk factors for acquisition of extended spectrum beta lactamase producing Escherichia coli and Klebsiella pneumoniae in North-Indian hospitals.2237–41. 10.1016/j.sjbs.2014.05.00625561881Shankar-SinhaS.ValenciaG. A.JanesB. K.RosenbergJ. K.WhitfieldC.BenderR. A. (2004). The Klebsiella pneumoniae O antigen contributes to bacteremia and lethality during murine pneumonia.721423–1430. 10.1128/iai.72.3.1423-1430.200414977947ShonA. S.BajwaR. P. S.RussoT. A. (2013). Hypervirulent (hypermucoviscous) Klebsiella pneumoniae: a new and dangerous breed.4107–118. 10.4161/viru.2271823302790StojkovicK.SzijártóV.KaszowskaM.NiedzielaT.HartlK.NagyG. (2017). Identification of d-Galactan-III as part of the lipopolysaccharide of serotype O1.8:684. 10.3389/fmicb.2017.0068428487676StruveC.BojerM.KrogfeltK. A. (2009). Identification of a conserved chromosomal region encoding Klebsiella pneumoniae type 1 and type 3 fimbriae and assessment of the role of fimbriae in pathogenicity.775016–5024. 10.1128/IAI.00585-589TaittC. R.LeskiT. A.ErwinD. P.OdundoE. A.KipkemoiN. C.NdonyeJ. N. (2017). Antimicrobial resistance of Klebsiella pneumoniae stool isolates circulating in Kenya.12:e0178880. 10.1371/journal.pone.017888028575064TsaiY.-K.FungC.-P.LinJ.-C.ChenJ.-H.ChangF.-Y.ChenT.-L. (2011). Klebsiella pneumoniae outer membrane porins OmpK35 and OmpK36 play roles in both antimicrobial resistance and virulence.551485–1493. 10.1128/AAC.01275-1210WickR. R.HeinzE.HoltK. E.WyresK. L. (2018). Kaptive web: user-friendly capsule and lipopolysaccharide serotype prediction for klebsiella genomes.56:e00197-18. 10.1128/JCM.00197-118WongJ. L. C.RomanoM.KerryL. E.KwongH.-S.LowW.-W.BrettS. J. (2019). OmpK36-mediated carbapenem resistance attenuates ST258 Klebsiella pneumoniae in vivo.10:3957.WuK.-M.LiL.-H.YanJ.-J.TsaoN.LiaoT.-L.TsaiH.-C. (2009). Genome sequencing and comparative analysis of Klebsiella pneumoniae NTUH-K2044, a strain causing liver abscess and meningitis.1914492–4501. 10.1128/jb.00315-0919447910WyresK. L.WickR. R.GorrieC.JenneyA.FolladorR.ThomsonN. R. (2016). Identification of capsule synthesis loci from whole genome data.2:e000102. 10.1099/mgen.0.00010228348840YangZ.RannalaB. (2012). Molecular phylogenetics: principles and practice.13303–314. 10.1038/nrg318622456349YehK.-M.KurupA.SiuL. K.KohY. L.FungC.-P.LinJ.-C. (2007). Capsular serotype K1 or K2, rather than magA and rmpA, is a major virulence determinant for Klebsiella pneumoniae liver abscess in Singapore and Taiwan.45466–471. 10.1128/JCM.01150-1156YuW.-L.KoW.-C.ChengK.-C.LeeC.-C.LaiC.-C.ChuangY.-C. (2008). Comparison of prevalence of virulence factors for Klebsiella pneumoniae liver abscesses between isolates with capsular K1/K2 and non-K1/K2 serotypes.621–6. 10.1016/j.diagmicrobio.2008.04.00718486404YuZ.QinW.LinJ.FangS.QiuJ. (2015). Antibacterial mechanisms of polymyxin and bacterial resistance.2015:679109. 10.1155/2015/67910925664322ZankariE.HasmanH.CosentinoS.VestergaardM.RasmussenS.LundO. (2012). Identification of acquired antimicrobial resistance genes.672640–2644. 10.1093/jac/dks26122782487ZhouK.LokateM.DeurenbergR. H.ArendsJ.Lo-Ten FoeJ.GrundmannH. (2015). Characterization of a CTX-M-15 producing Klebsiella pneumoniae outbreak strain assigned to a novel sequence Type (1427).6:1250. 10.3389/fmicb.2015.0125026617589