Back to Journals » Infection and Drug Resistance » Volume 17
Global Emergence and Genomic Epidemiology of blaNDM-Carrying Klebsiella variicola
Authors Li L, Zhang Y, Wang W, Chen Y, He F , Yu Y
Received 26 February 2024
Accepted for publication 1 May 2024
Published 14 May 2024 Volume 2024:17 Pages 1893—1901
DOI https://doi.org/10.2147/IDR.S460569
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 2
Editor who approved publication: Professor Suresh Antony
Lirong Li,1,* Yawen Zhang,1,* Weizhong Wang,1 Yanmin Chen,1 Fang He,1 Yan Yu2
1Laboratory Medicine Center, Department of Clinical Laboratory, Zhejiang Provincial People’s Hospital (Affiliated People’s Hospital), Hangzhou Medical College, Hangzhou, Zhejiang, People’s Republic of China; 2Center for Rehabilitation Medicine, Department of Ophthalmology, Zhejiang Provincial People’s Hospital (Affiliated People’s Hospital), Hangzhou Medical College, Hangzhou, Zhejiang, People’s Republic of China
*These authors contributed equally to this work
Correspondence: Yan Yu; Fang He, Email [email protected]; [email protected]
Purpose: Klebsiella variicola has emerged as a human pathogen in the past decade. Here, we present findings related to a K. variicola strain carrying the blaNDM-1 gene, which was isolated from a urinary tract infection in China. Global transmission dynamics and genomic epidemiology of blaNDM-carrying K. variicola were further investigated.
Material and Methods: The complete genome sequence of the strain was determined using the Illumina NovaSeq 6000 and Nanopore MinION sequencer. Genomic features and resistance mechanisms were analyzed through diverse bioinformatics approaches. Additionally, genome sequences of K. variicola strains carrying blaNDM were retrieved from the NCBI database, and a comprehensive analysis of the global dissemination trends of these strains was conducted.
Results: K. variicola strain 353 demonstrated resistance to multiple antimicrobials, including carbapenems. Within its genome, we identified fourteen antimicrobial resistance genes associated with β-lactam, aminoglycoside, fosfomycin, quinolone, trimethoprim, rifamycin, and sulfonamide resistance. The carbapenem-resistant gene blaNDM-1 was located on an IncU-type plasmid spanning 294,608 bp and flanked by ISCR1 and IS 26. Downstream of blaNDM-1, we identified an Intl1 element housing numerous antibiotic resistance genes. A comprehensive search of the NCBI database revealed 72 K. variicola strains carrying blaNDM from twelve different countries, predominantly from clinical sources, with the highest prevalence observed in the USA and China. A total of 28 distinct sequence types (STs) were identified, with ST115 being the most prevalent, followed by ST60.
Conclusion: In summary, this study presents the genomic characterization of a K. variicola strain carrying blaNDM-1 on an IncU-type plasmid. The research highlights the global dissemination of blaNDM-carrying K. variicola, observed in both healthcare settings and natural environments. Our data have revealed a diverse array of antimicrobial resistance determinants in K. variicola, providing valuable insights that could aid in the development of strategies for the prevention, diagnosis, and treatment of K. variicola infections.
Keywords: Klebsiella variicola, blaNDM-1, whole-genome sequencing, IncU type plasmid, urinary tract infection
Introduction
Klebsiella variicola is a gram-negative, facultative anaerobe, nonsporogenic, and nonmotile rod-shaped bacterium that forms round, convex, and smooth colonies.1 Initially discovered in bananas in 2004, K. variicola is commonly present in agriculturally sourced soils, plants, freshwater, and sewage.2 It frequently contributes to nitrogen fixation and promotes plant growth. Moreover, it is recognized as a commensal of insects and can act as a pathogen for both plants and animals.3–5 Belonging to the Klebsiella genus, K. variicola has historically been prone to misidentification as Klebsiella pneumoniae through traditional biochemical methods. This challenge persisted until the adoption of advanced techniques such as matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and whole-genome sequencing for accurate species identification.6,7 The past misidentification, particularly within the K. pneumoniae complex, has hindered comprehensive research on K. variicola within the healthcare system.
K. variicola has emerged as a noteworthy human pathogen, witnessing a consistent rise in clinical infections in recent years.8 Reports have extensively documented infections involving isolates from various sources, including fecal, blood, sputum, vaginal, and urine samples.9,10 Notably, hypermucoviscous strains of K. variicola were initially identified in 2015, intensifying the concerns surrounding its clinical impact.11,12 Adding to these apprehensions, a study highlighted that among hospitalized adults with bloodstream infections, the mortality rate associated with K. variicola exceeded that linked to K. pneumoniae.13 This underscores the increasing significance of K. variicola as a clinically relevant pathogen, prompting a closer examination of its implications in healthcare settings.
In this investigation, the K. variicola strain 353, which carries the blaNDM-1 gene, was isolated from a urinary tract infection in a male patient hospitalized in the department of rehabilitation of a teaching hospital in China. The isolate was preliminarily identified using the VITEK MS system (bioMérieux, France) and was further confirmed by whole-genome sequencing. Bioinformatics analysis was undertaken to delve into the genetic characteristics of both the strain and the plasmid carrying the blaNDM-1 gene. The global transmission dynamics and genomic epidemiology of blaNDM-carrying K. variicola were further investigated.
Materials and Methods
Antimicrobial Susceptibility Test
Antimicrobial susceptibility testing was performed using the VITEK 2 system (bioMérieux, France) with Gram-negative antimicrobial susceptibility testing cards (AST-GN13) and the Etest method. The testing procedures adhered to the guidelines established by the Clinical and Laboratory Standards Institute (CLSI) M100, 33rd edition. Breakpoints were interpreted in accordance with the recommendations outlined in the CLSI guidelines. In cases where CLSI breakpoints were unavailable for colistin and tigecycline, interpretations for colistin minimum inhibitory concentration (MIC) followed the guidelines of the European Committee on Antimicrobial Susceptibility Testing (EUCAST), while standards set by the US Food and Drug Administration (FDA) were utilized for tigecycline.
Whole-Genome Sequencing
Whole-genome sequence of the strain was determined utilizing the Illumina NovaSeq 6000 platform (Illumina Inc., San Diego, CA, USA) in the 150-bp paired-end sequencing mode, with an average sequencing depth of ≥100×. Additionally, long-read sequencing was performed using a Nanopore MinION sequencer (Nanopore, Oxford, UK). Both short Illumina reads and long MinION reads underwent hybrid assembly using Unicycler (v0.4.7) in the conservative mode. This process resulted in complete circular contigs, which underwent further refinement and correction using Pilon with Illumina reads through multiple rounds of iteration until no further changes were detected. The resultant complete genome sequence was subsequently automatically annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP) server.
Genomic Features and Plasmid Characterization
The investigation into the antimicrobial resistance genes and plasmid replicons of the strain was carried out using the BacWGSTdb server.14,15 MLST (Multi Locus Sequence Typing) analysis was conducted utilizing the database available at http://mlstkv.insp.mx/.16 Circular comparisons, represented by concentric rings, were performed to examine the blaNDM-1-carrying plasmid and its similarity to analogous plasmids. This comparative analysis was visualized using the BLAST Ring Image Generator (BRIG).17
Phylogenetic Analysis
The phylogenetic relationship between K. variicola 353 and other K. variicola strains obtained from the NCBI GenBank database was assessed using the BacWGSTdb server.14,15 This server employs single nucleotide polymorphism (SNP) approaches to analyze the phylogenetic relationship of the uploaded genome sequence with sequences available in the database. The resulting phylogenetic tree was enhanced for visual clarity using iTOL.18
Nucleotide Sequence Accession Numbers
The complete genome sequence of K.variicola 353 has been submitted to the NCBI GenBank database and is assigned the accession number CP141632-CP141634.
Ethical Approval
This study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Zhejiang Provincial People’s Hospital (Ethics approval number 2019KY244).
Results and Discussion
K. variicola strain 353 was isolated from a urine sample obtained from a male patient who was admitted to a tertiary hospital in China following a cerebral hemorrhage. The genome of the K. variicola strain 353 is composed of three contigs, totaling 6,213,895 bp. Notably, one of these contigs, designated as contig 1 and spanning 5,727,068 bp, is associated with the chromosome, while the remaining two contigs correspond to plasmids (contig pCRKP353-NDM1: 294,608 bp and contig 3: 192,219 bp). The genome harbors two distinct plasmid replicons: one on plasmid pCRKP353-NDM1 (IncU) and another on contig 3 (IncFIB(K)). Analysis conducted through the PGAP server yielded predictions for 6001 protein-coding sequences, 88 tRNA genes, and 25 rRNA operons.
The antibiotic susceptibility profiles, as depicted in Table S1, reveal that K. variicola 353 exhibited resistance to a broad spectrum of antibiotics, including ceftazidime, ceftriaxone, cefotetan, cefazolin, cefepime, ciprofloxacin, levofloxacin, ertapenem, imipenem, meropenem, ampicillin/sulbactam, and sulfamethoxazole/trimethoprim. However, it demonstrated susceptibility to amikacin, gentamicin, aztreonam, tigecycline, and colistin. Table 1 outlines the resistance genes identified in the genome of the isolate. Notably, β-lactam resistance genes blaLEN-16 and blaNDM-1 were identified, along with aminoglycoside resistance genes aph(3’)-Ia, aadA16, aadA2 and aac(6’)-Ib-cr, fosfomycin resistance gene fosA, quinolone resistance genes oqxB, oqxA, and qnrS1, trimethoprim resistance gene dfrA27, rifamycin resistance gene arr-3, and sulfonamide resistance gene sul1. Except for oqxA, oqxB, blaLEN-16, and fosA, which are located on the chromosome, all other resistance genes, including blaNDM-1, are situated on plasmid pCRKP353-NDM1.
![]() |
Table 1 Antimicrobial Resistance Genes in Klebsiella Variicola 353 |
The plasmid pCRKP353-NDM1, carrying the blaNDM gene, was identified as an IncU-type plasmid. Two copies of IS26 are situated upstream of blaNDM-1, while two ISCR1 (insertion sequence common region 1) elements flank blaNDM-1. Downstream of blaNDM-1, an Intl1 (Class 1 Integron) is present, housing numerous antibiotic resistance genes such as sul1, aadA16, dfrA27, arr-3 and aac(6’)-Ib-cr. IS26 typically plays a pivotal role in disseminating antibiotic resistance genes within Gram-negative bacteria, contributing to their widespread distribution.19–21 Additionally, ISCR1 is closely associated with various antibiotic resistance determinants, underscoring its significance in this context.22 This implies that both IS26 and ISCR1 are crucial players in the propagation of blaNDM. The presence of Intl1 (Class 1 Integron) downstream of blaNDM-1 enables the plasmid to harbor multiple drug-resistant genes, transforming it into a multi-drug-resistant plasmid. This emphasizes the role of Intl1 in facilitating the accumulation of diverse drug resistance determinants within the plasmid.
Comparative analysis of pCRKP353-NDM1 with similar plasmids obtained from the NCBI database was performed using the Basic Local Alignment Search Tool (Figure 1). These plasmids, all falling under the IncU type, include three that harbor the blaNDM gene (Table S2). IncU plasmids are distinguished for their extensive host ranges, enabling the efficient transmission of antibiotic resistance genes owing to their robust binding and mobility features.23 Notably, characterized by their high binding and mobility attributes, IncU plasmids demonstrate resilience and adept replication capabilities across diverse bacterial species. Consequently, they serve as crucial vectors in the dissemination of antibiotic resistance genes. It is noteworthy that IncU plasmids carrying blaNDM-1 genes have begun to proliferate across diverse genera, encompassing Raoultella ornithinolytica, Enterobacter cloacae, and K. variicola.
In order to gain a more profound understanding of global transmission dynamics and genomic epidemiology of blaNDM-carrying K. variicola, an extensive search was conducted within the NCBI pathogen database for strains harboring the blaNDM gene until December 1, 2023. A total of 72 blaNDM-carrying K. variicola strains were identified from the NCBI database (Table 2). These strains were discovered in 12 countries worldwide, with the highest prevalence observed in the USA (32 strains), followed by China (18 strains), Bangladesh (4 strains), and three strains each from the United Kingdom, Switzerland, and South Korea. Additionally, two strains each were found in Australia and Brazil, as well as one strain each from Lebanon, Germany and Canada. This global distribution highlights the widespread dissemination of blaNDM-carrying K. variicola. Among these strains, five sub-types of blaNDM were identified: 38 strains carried blaNDM-1, 13 strains carried blaNDM-5, 11 strains carried blaNDM-4, 3 strains carried blaNDM-9, and 1 strain carried blaNDM-18. According to the MLST (Multi Locus Sequence Typing) results available at http://mlstkv.insp.mx/,16 KV353 is assigned to a newly identified, unnamed sequence type (ST) distinguished by leuS60, pgi39, pgk1, phoE29, pyrG1, rpoB1 and fusA4. Among the 72 blaNDM-carrying K. variicola strains, excluding 20 strains with unclassified sequence types, a total of 28 distinct ST types were identified (Figure 2). ST115 was the most prevalent, comprising 12 strains, followed by ST60 (4 strains), ST10 (3 strains), ST20 (3 strains), ST64 (3 strains), and ST277 (3 strains). Phylogenetic analysis revealed that K. variicola 353 belongs to a distinct clone. Furthermore, eleven strains originating from the United States were identified to be part of a single clone, exhibiting single nucleotide polymorphism (SNP) differences of ≤20. This observation suggests that clonal transmission, possibly associated with nosocomial infection, could be contributing to the epidemic of blaNDM-carrying K. variicola.
![]() |
Table 2 Clinical Metadata for 72 Klebsiella Variicola Strains Carrying blaNDM Retrieved from the NCBI Database |
Relevant analyses were conducted on the blaNDM-carrying plasmid. Subsequent investigation into these blaNDM-carrying strains, among which 34 had complete genome sequences, revealed that all blaNDM-carrying plasmids were located on plasmids. Specifically, 13 were classified as IncX3 type, 10 as IncA/C2 type, 5 as IncFII type, 3 strains as IncHI1B type, with one each belonging to IncFIA, IncFIB, and IncN types. The occurrence of IncU type blaNDM-carrying plasmid in K. variicola is uncommon.
Remarkably, seven strains from distinct regions were identified to carry multiple carbapenemases, with three strains co-harboring blaNDM and blaIMP, and four strains co-harboring blaNDM and blaKPC. Over the past decade, there has been a substantial increase in the proportion of K. pneumoniae strains carrying multiple-carbapenemase genes, presenting an elevated threat to public health.24 The current study suggests that K. variicola strains carrying multiple-carbapenemase genes are also undergoing global dissemination. The existence and prevalence of K. variicola with multiple carbapenemases pose formidable challenges for clinical treatment.
Among the 72 strains of K. variicola obtained from NCBI, there was a notably higher prevalence in urine (18/72), blood (10/72), and sputum (6/72) samples. Additionally, K. variicola strains were detected in environmental sewage and fields. It is noteworthy, that one strain carrying the blaNDM-1 gene was detected in hospital wastewater in China, and three strains carrying the blaNDM-9 gene from South Korea were found in river water, while the remaining strains were of clinical origin. This implies that bacteria carrying the blaNDM gene have transferred between clinical settings and the natural environment.
In summary, our findings reveal the global dissemination of K. variicola carrying the blaNDM gene in both healthcare settings and natural environments. These data have unveiled a diverse array of antimicrobial resistance determinants in K. variicola, offering valuable insights that may contribute to the development of strategies for the prevention, diagnosis, and treatment of K. variicola infections.
Acknowledgments
This work was supported by the Traditional Chinese Medicine Science and Technology Project of Zhejiang Province, China (Grant number 2024ZL253), and the Medicine and Health Science and Technology Project of Zhejiang Province, China (Grant number 2021KY475).
Disclosure
The authors report no conflicts of interest in this work.
References
1. Lin L, Wei C, Chen M, et al. Complete genome sequence of endophytic nitrogen-fixing Klebsiella variicola strain DX120E. Stand Genomic Sci. 2015;10:22. doi:10.1186/s40793-015-0004-2
2. Rosenblueth M, Martinez L, Silva J, Martinez-Romero E. Klebsiella variicola, a novel species with clinical and plant-associated isolates. Syst Appl Microbiol. 2004;27(1):27–35. doi:10.1078/0723-2020-00261
3. Martinez-Romero E, Silva-Sanchez J, Barrios H, et al. Draft Genome Sequences of Klebsiella variicola Plant Isolates. Genome Announc. 2015;3(5):15. doi:10.1128/genomeA.01015-15
4. Martinez-Romero E, Rodriguez-Medina N, Beltran-Rojel M, et al. Genome misclassification of Klebsiella variicola and Klebsiella quasipneumoniae isolated from plants, animals and humans. Salud Publica Mex. 2018;60(1):56–62. doi:10.21149/8149
5. Pinto-Tomas AA, Anderson MA, Suen G, et al. Symbiotic nitrogen fixation in the fungus gardens of leaf-cutter ants. Science. 2009;326(5956):1120–1123. doi:10.1126/science.1173036
6. Ohama Y, Nomura Y, Mizoguchi M, Higurashi Y, Okamoto K, Harada S. Accurate Identification of Klebsiella variicola by MALDI-TOF Mass Spectrometry in Clinical Microbiology Laboratories. Microbiol Spectr. 2022;10(5):e0284422. doi:10.1128/spectrum.02844-22
7. Potter RF, Lainhart W, Twentyman J, et al. Population Structure, Antibiotic Resistance, and Uropathogenicity of Klebsiella variicola. mBio. 2018;9(6):2481. doi:10.1128/mBio.02481-18
8. Rodriguez-Medina N, Barrios-Camacho H, Duran-Bedolla J, Garza-Ramos U. Klebsiella variicola: an emerging pathogen in humans. Emerg Microbes Infect. 2019;8(1):973–988. doi:10.1080/22221751.2019.1634981
9. Holt KE, Wertheim H, Zadoks RN, et al. Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc Natl Acad Sci U S A. 2015;112(27):E3574–81. doi:10.1073/pnas.1501049112
10. Alves MS, Dias RC, de Castro AC, Riley LW, Moreira BM. Identification of clinical isolates of indole-positive and indole-negative Klebsiella spp. J Clin Microbiol. 2006;44(10):3640–3646. doi:10.1128/JCM.00940-06
11. Garza-Ramos U, Silva-Sanchez J, Barrios H, Rodriguez-Medina N, Martinez-Barnetche J, Andrade V. Draft Genome Sequence of the First Hypermucoviscous Klebsiella variicola Clinical Isolate. Genome Announc. 2015;3(2):14. doi:10.1128/genomeA.01352-14
12. Rodriguez-Medina N, Martinez-Romero E, De la Cruz MA, et al. A Klebsiella variicola Plasmid Confers Hypermucoviscosity-Like Phenotype and Alters Capsule Production and Virulence. Front Microbiol. 2020;11:579612. doi:10.3389/fmicb.2020.579612
13. Maatallah M, Vading M, Kabir MH, et al. Klebsiella variicola is a frequent cause of bloodstream infection in the Stockholm area, and associated with higher mortality compared to K. pneumoniae. PLoS One. 2014;9(11):e113539. doi:10.1371/journal.pone.0113539
14. Ruan Z, Feng Y. BacWGSTdb, a database for genotyping and source tracking bacterial pathogens. Nucleic Acids Res. 2016;44(D1):D682–7. doi:10.1093/nar/gkv1004
15. Ruan Z, Yu Y, Feng Y. The global dissemination of bacterial infections necessitates the study of reverse genomic epidemiology. Briefings Bioinf. 2020;21(2):741–750. doi:10.1093/bib/bbz010
16. Barrios-Camacho H, Aguilar-Vera A, Beltran-Rojel M, et al. Molecular epidemiology of Klebsiella variicola obtained from different sources. Sci Rep. 2019;9(1):10610. doi:10.1038/s41598-019-46998-9
17. Alikhan NF, Petty NK, Ben Zakour NL, Beatson SA. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:402. doi:10.1186/1471-2164-12-402
18. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–W296. doi:10.1093/nar/gkab301
19. Partridge SR, Kwong SM, Firth N, Jensen SO. Mobile Genetic Elements Associated with Antimicrobial Resistance. Clin Microbiol Rev. 2018;31(4):17. doi:10.1128/CMR.00088-17
20. Harmer CJ, Moran RA, Hall RM. Movement of IS26-associated antibiotic resistance genes occurs via a translocatable unit that includes a single IS26 and preferentially inserts adjacent to another IS26. mBio. 2014;5(5):e01801–14. doi:10.1128/mBio.01801-14
21. He S, Hickman AB, Varani AM, et al. Insertion Sequence IS26 Reorganizes Plasmids in Clinically Isolated Multidrug-Resistant Bacteria by Replicative Transposition. mBio. 2015;6(3):e00762. doi:10.1128/mBio.00762-15
22. Toleman MA, Bennett PM, Walsh TR. Common regions e.g. orf513 and antibiotic resistance: IS91-like elements evolving complex class 1 integrons. J Antimicrob Chemother. 2006;58(1):1–6. doi:10.1093/jac/dkl204
23. Dang B, Xu Y, Mao D, Luo Y. Complete nucleotide sequence of plasmid pNA6 reveals the high plasticity of IncU family plasmids. Gene. 2016;591(1):74–79. doi:10.1016/j.gene.2016.06.060
24. Guo H, Wu Y, Li L, Wang J, Xu J, He F. Global emergence of carbapenem-resistant Klebsiella pneumoniae co-carrying multiple carbapenemases. Comput Struct Biotechnol J. 2023;21:3557–3563. doi:10.1016/j.csbj.2023.07.013
© 2024 The Author(s). This work is published and licensed by Dove Medical Press Limited. The
full terms of this license are available at https://www.dovepress.com/terms.php
and incorporate the Creative Commons Attribution
- Non Commercial (unported, v3.0) License.
By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted
without any further permission from Dove Medical Press Limited, provided the work is properly
attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.