Skip to main content

Genome analyses of colistin-resistant high-risk blaNDM-5 producing Klebsiella pneumoniae ST147 and Pseudomonas aeruginosa ST235 and ST357 in clinical settings

Abstract

Background

Colistin is a last-resort antibiotic used in extreme cases of multi-drug resistant (MDR) Gram-negative bacterial infections. Colistin resistance has increased in recent years and often goes undetected due to the inefficiency of predominantly used standard antibiotic susceptibility tests (AST). To address this challenge, we aimed to detect the prevalence of colistin resistance strains through both Vitek®2 and broth micro-dilution. We investigated 1748 blood, tracheal aspirate, and pleural fluid samples from the Intensive Care Unit (ICU), Neonatal Intensive Care Unit (NICU), and Tuberculosis and Respiratory Disease centre (TBRD) in an India hospital. Whole-genome sequencing (WGS) of extremely drug-resitant (XDR) and pan-drug resistant (PDR) strains revealed the resistance mechanisms through the Resistance Gene Identifier (RGI.v6.0.0) and Snippy.v4.6.0. Abricate.v1.0.1, PlasmidFinder.v2.1, MobileElementFinder.v1.0.3 etc. detected virulence factors, and mobile genetic elements associated to uncover the pathogenecity and the role of horizontal gene transfer (HGT).

Results

This study reveals compelling insights into colistin resistance among global high-risk clinical isolates: Klebsiella pneumoniae ST147 (16/20), Pseudomonas aeruginosa ST235 (3/20), and ST357 (1/20). Vitek®2 found 6 colistin-resistant strains (minimum inhibitory concentrations, MIC = 4 μg/mL), while broth microdilution identified 48 (MIC = 32–128 μg/mL), adhering to CLSI guidelines. Despite the absence of mobile colistin resistance (mcr) genes, mechanisms underlying colistin resistance included mgrB deletion, phosphoethanolamine transferases arnT, eptB, ompA, and mutations in pmrB (T246A, R256G) and eptA (V50L, A135P, I138V, C27F) in K. pneumoniae. P. aeruginosa harbored phosphoethanolamine transferases basS/pmrb, basR, arnA, cprR, cprS, alongside pmrB (G362S), and parS (H398R) mutations. Both strains carried diverse clinically relevant antimicrobial resistance genes (ARGs), including plasmid-mediated blaNDM-5 (K. pneumoniae ST147) and chromosomally mediated blaNDM-1 (P. aeruginosa ST357).

Conclusion

The global surge in MDR, XDR and PDR bacteria necessitates last-resort antibiotics such as colistin. However, escalating resistance, particularly to colistin, presents a critical challenge. Inefficient colistin resistance detection methods, including Vitek2, alongside limited surveillance resources, accentuate the need for improved strategies. Whole-genome sequencing revealed alarming colistin resistance among K. pneumoniae and P. aeruginosa in an Indian hospital. The identification of XDR and PDR strains underscores urgency for enhanced surveillance and infection control. SNP analysis elucidated resistance mechanisms, highlighting the complexity of combatting resistance.

Peer Review reports

Background

The antibiotic development pipeline has run dry, failing to compete with the constantly evolving antimicrobial resistance (AMR) crisis [1]. This ‘silent pandemic’ led to approximately 5 million mortalities in 2019 and is predicted to escalate drastically in 2050 with 10 million annual deaths [2]. Multidrug resistance (MDR) is a growing global public health concern, with multidrug-resistant bacterial pathogens emerging from various sources. This underscores the crucial need for judicious antibiotic use. Routine antimicrobial susceptibility testing is vital for identifying the appropriate antibiotics and screening for emerging MDR strains [3,4,5,6,7,8]. Global dissemination of MDR, hypervirulent, and high-risk clones poses a significant healthcare challenge, particularly in low and middle-income countries (LMICs) [2, 9]. Colistin, a crucial last-resort antibiotic for MDR gram-negative infections, has witnessed a surge in resistance over the past decade. The mechanisms of colistin resistance, though not fully elucidated, involve plasmid-mediated mcr genes, intrinsic factors, and genetic mutations [10, 11]. Since colistin’s bactericidal activity is associated with penetrating and disrupting the bacterial cell membrane, reduced colistin susceptibility is mediated by adding cationic molecules 4-deoxyaminoarabinose (LAra4N) and phosphoethanolamine (pEtN) to the lipopolysaccharide surface (LPS), particularly lipid A, extrusion of colistin through efflux pumps, and overexpression of outer membrane proteins. In Klebsiella pneumoniae, a transmembrane protein mgrB (47 amino acids) is a crucial factor intricately involved in the colistin resistance pathway. Insertional inactivation by interrupting IS element, truncation, or complete deletion of mgrB alters colistin susceptibility [12,13,14,15,16]. mgrB is the negative regulator of the two-component system PhoPQ. PhoQ is a sensor kinase activated by the low extracellular magnesium (Mg2+), acidic pH (5.5), or the presence of cationic antimicrobial peptides (AMPs). PhoQ activates PhoP through phosphorylation, which triggers the expression of the pmrHFIJKLM operon, also known as pmrF, pbg, or arn operon. The pmrHFIJKLM operon adds LAra4N to the lipid A of LPS. Several other factors contribute towards colistin resistance, notably, pEtN transferases such as eptB and arnT, alterations in two-component system pmrAB, phoPQ, and crrAB, parRS, colRS, upregulation of acrA, tolC, and soxRS, efflux pumps including KpnEF, KpnG, acrAB and sap protein systems, and porin ompA [12, 13, 15,16,17,18,19,20]. In Pseudomonas aeruginosa, LPS modifications reducing the surface anionic charge, the arnBCADTEF-pmrE (pmrHFIJKLME) encoding enzymes, which add LAra4N to lipid A, and alterations in PmrAB, PhoPQ, ParRS, and CprRS two-component regulatory systems are important mechanisms inducing colistin resistance [21].

The World Health Organisation (WHO) recognizes carbapenem-resistant K. pneumoniae and P. aeruginosa as ‘critical priority’ pathogens [22]. K. pneumoniae is a typical perpetrator of nosocomial, community-acquired, and opportunistic infections such as pneumonia, urinary tract infections (UTIs), bloodstream infections, and sepsis [23,24,25,26]. Klebsiella pneumoniae’s antimicrobial resistance primarily stems from: i. Enzyme production, notably β-lactamases like blaNDM, blaKPC, blaOXA. ii. Reduced cell permeability due to loss of outer membrane proteins (OmpA, OmpK35, OmpK36, OmpK37). iii. Overexpression of efflux pumps, expelling antibiotics from the cell. iv. Target modification, altering the antimicrobial agents’ targets [27]. β-lactamase emergence has significantly reduced the effectiveness of multiple antibiotics like carbapenems, cephems, clavams, monobactams, and penems, posing major public health challenges [28, 29]. The clonal type K. pneumoniae ST147 is a high-risk clone associated with MDR and pan-drug resistant chronic infections [9]. It was first reported in Hungary and Spain, genotypically characterized with fluoroquinolone resistance induced by gyrAS83I and parCS80I QRDR mutations and extended-spectrum β-lactamase gene blaCTX-M-15 [30, 31]. Later, it expanded globally, acquiring resistance against major antimicrobial classes, including ESBLs and carbapenems [9, 32, 33]. K. pneumoniae ST147 is endemic to India, Italy, Greece, and certain North African countries [9]. The high genomic plasticity, virulence, and extensive resistance acquisition make K. pneumoniae ST147 a potential epidemic trigger. K. pneumoniae ST147 carrying blaNDM-1 caused prolonged large outbreaks in Poland and the USA (2015–19) [34, 35]. Hypervirulent K. pneumoniae ST147 armed with a chimeric plasmid carrying resistance and virulence determinants caused an outbreak in Tuscany, Italy [36]. Beta-lactamase blaNDM-5 is endemic to India, and reportedly, K. pneumoniae ST147 acquired blaNDM-5 and blaOXA-181 in India, later disseminating it globally [37,38,39]. In India, colistin-resistant K. pneumoniae ST147 isolates have been reported in sepsis, acute kidney injury, and pulmonary tuberculosis [40,41,42].

P. aeruginosa is an urgent threat associated with severe nosocomial infections, bacteremia, sepsis, hospital-acquired pneumonia, and respiratory failure, especially in immunocompromised patients, diabetic ulcers, burn wounds, corneal ulcers, and surgical wounds [43]. It colonizes the lungs of cystic fibrosis (CF) patients, leading to pulmonary failure [44]. In Pseudomonas aeruginosa, antibiotic resistance is driven by three main mechanisms: intrinsic, acquired, and adaptive. i. Intrinsic resistance is due to reduced outer membrane permeability and presence of antibiotic efflux pumps. ii. Acquired resistance arises from horizontal gene transfer of antibiotic resistance genes (ARGs), mutations leading to increased expression of efflux pumps and β-lactamases, and alterations in antibiotic target sites. Iii. Adaptive resistance, which is transient and induced by external factors like stress or specific antibiotics, involves regulatory pathways, changes in gene expression, protein production, or modifications to antibiotic target [45]. P. aeruginosa ST235 and ST357 are widespread high-risk clones, with ST235 being the most prevalent. Their capability to acquire resistance against major antibiotics, including beta-lactamases, carbapenemases, and aminoglycosides, is a significant factor associated with high morbidity and mortality. The production of ExoU cytotoxin is a major virulence factor contributing to high lethal potential [46, 47].

Virulence factors (V.F.s) are crucial for enhancing bacterial infectivity by promoting strong attachment and colonization within host cells, disrupting tissue integration, evading the host immune response, and exploiting nutrient availability [48]. Classical and Hypervirulent Klebsiella pneumoniae (HVKP) strains differ in virulence, identified through a string test and gene-carrying plasmids. Key V.F.s include fimbriae adhesion, capsule, lipopolysaccharide (LPS), and siderophores, aiding the attachment, evasion from immunity, enhanced survival and colonization [27].

In P. aeruginosa, virulence factors (V.F.s) are categorized into three main categories: i. Bacterial surface structures including pili, flagella, secretion systems (Types 1, 2, 3, 5, and 6), and outer membrane proteins (e.g., lipopolysaccharide). These aid in attachment, motility, immune evasion, toxin delivery, and biofilm formation. ii. Secreted factors such as exopolysaccharides (Alginate, Pel, Psl), siderophores (Pyoverdine, Pyochelin), proteases (AprA, LasA, LasB, Protease IV), and toxins (ExoS, ExoT, ExoU, ExoY). These contribute to adhesion, immune evasion, growth promotion, and virulence induction. Iii. Virulence factor-mediated interactions between bacterial cells regulate quorum sensing and biofilm production, aiding in evasion of host immune responses [49].

Biofilms, intricate three-dimensional structures formed by P. aeruginosa and K. pneumoniae, enhance their pathogenicity and resistance to antimicrobial agents. The extracellular matrix, comprising polysaccharides, proteins, nucleic acids, lipids, and extracellular DNA (eDNA), supports bacterial growth within biofilms. In P. aeruginosa, virulence factors like alginate, Psl, and Pel aid in biofilm formation [49, 50]. Similarly, in K. pneumoniae, biofilm formation is influenced by factors such as the capsule, type 1 and type 3 fimbriae, and the Escherichia coli common pilus (ECP) fimbriae gene cluster [51]. Quorum sensing mechanisms regulate biofilm formation by controlling the synthesis of fimbriae, exopolysaccharides, and adhesins through signaling molecules, profoundly impacting both inter- and intra-species communication [52, 53].

Indian hospitals are one of the largest consumers of polymyxins, driven by the increasing cases of carbapenem-resistant K. pneumoniae, P. aeruginosa, and Acinetobacter baumannii infections [54]. Colistin resistance often goes unnoticed and is underreported in hospitals in India due to the incompetence of predominantly used standard antibiotic susceptibility tests such as Vitek and disk diffusion [55]. In this comprehensive whole-genome study, we investigated clinical isolates sourced from blood cultures of sepsis patients and pleural fluid, as well as tracheal aspirate samples from tuberculosis patients, collected during 2020–21. Our primary objective was to elucidate the factors contributing to colistin resistance within an Indian hospital setting, focusing on predominant bacterial species and the prevalence of other antibiotic resistance determinants. Through thorough analysis, we aimed to identify mechanisms of colistin resistance, including horizontal gene transfer facilitated by plasmid-mediated antibiotic resistance genes (ARGs), insertion sequence (IS) elements, virulence determinants, genomic islands, and bacteriophages. Additionally, we investigated the presence of CRISPR elements to assess bacterial immune responses against foreign genetic materials. Furthermore, we employed pangenome analysis to characterize the phylogenetic diversity of the bacterial isolates under study. Moreover, we sought to observe the differences between colistin resistance detection using Vitek2, which has previously yielded false negatives, and antimicrobial susceptibility testing (AST) through broth micro-dilution, aiming to provide a more accurate assessment of colistin resistance in our sample population.

Methods

Wet lab section

Isolation and antimicrobial susceptibility tests (AST) through VITEK2 and broth micro-dilution

The blood, pleural fluid, and tracheal aspirate (1748 samples) were inoculated on Luria Bertani agar (Hi Media, India) and incubated for 18 h at 37 °C. The pure colonies were picked by a sterile loop, inoculated in 5 mL LB broth and incubated on a shaker at 37 °C for 18 h. 0.3 mL of overnight culture was mixed with 0.7 mL of glycerol in sterile vials and stored at − 80 °C for future use. Gram-staining was performed to detect gram-positive and gram-negative isolates [49, 56]. The preliminary identification and antimicrobial susceptibility test against colistin and other antimicrobials (amikacin, amoxicillin, ampicillin, cefepime, cefoperazone/sulbactam, ceftazidime, ceftriaxone, ciprofloxacin, colistin, gentamicin, imipenem, levofloxacin, meropenem, piperacillin/tazobactam, ticarcillin/clavulanic acid+ cefotaxime, tigecycline, and trimethoprim/sulfamethoxazole) on isolates was performed using a semi-automated commercial system VITEK2 Advanced Expert system (BioMérieux, Marcy l’Étoile, France), using the European Committee on Antimicrobial Susceptibility Testing (EUCAST) breakpoints (http://www.eucast.org/clinical_breakpoints) [57]. The minimum inhibitory concentration (MIC) of colistin (Hi Media, India) and other antimicrobials including amikacin (Hi Media, India), aztreonam (Sigma-Aldrich, USA), cefazolin (Hi Media, India), cefepime (Hi Media, India), cefoxitin (Sigma-Aldrich, USA), ceftazidime (Sigma-Aldrich, USA), ceftriaxone (Sigma-Aldrich, USA), ciprofloxacin (Sigma-Aldrich, USA), colistin (Hi Media, India), Fosfomycin (Hi Media, India), gentamicin (Hi Media, India), imipenem (Hi Media, India), levofloxacin (Hi Media, India), meropenem (Hi Media, India), piperacillin (Hi Media, India), tetracycline (Hi Media, India), and ticarcillin-clavulanate (Hi Media, India) against the isolates was determined by broth microdilution as well, using cation-adjusted Mueller-Hinton broth (MHCB) [55, 58] (Additional Table 1). The results were interpreted according to the Clinical and Laboratory Standards Institute (CLSI) guidelines [59]. The classification of isolates into multi-drug resistant (MDR), extensively drug-resistant (XDR) and pandrug-resistant (PDR) was done according to Magiorakos et al. [60]. The MAR (Multiple Antibiotic Resistance) index was computed and analyzed following the methodology outlined by Krumperman, employing the formula: a/b. Here, ‘a’ denotes the count of antibiotics for which an isolate exhibited resistance, while ‘b’ stands for the total number of antibiotics subjected to testing [61]. The whole-cell DNA from the strains was obtained by collecting pure colonies from a culture plate and suspending them in 100 μL of nuclease-free water. Subsequently, the mixture was heated at 95 °C for 15 minutes in a thermal block and then centrifuged at 10,000 rpm at 4 °C for 10 minutes [62]. Gene amplification was performed using polymerase chain reaction (PCR) with gene-specific primers targeting mcr-1. The primers used were mcr-F(forward) 5′-CGG TCA GTC CGT TTG TTC-3′ and mcr-R (reverse) 5′-CTT GGT CGG TCT GTA GGG-3′ [63]. The PCR products were mixed with gel loading dye and run on 1% agarose gels in 1x TAE buffer. Ethidium bromide was added for visualization, and gel images were captured using a gel documentation system (Gel Doc) .

Dry lab section

DNA extraction, library preparation, sequencing, and quality check

The DNA of colistin-resistant isolates was extracted using a QIAGEN DNeasy PowerSoil kit (Qiagen, Germany) according to the manufacturer’s protocol and subjected to a Nano Drop™ 2000 Spectrophotometer (ThermoFisher Scientific, USA) for checking RNA and protein contaminants [64]. The standard set for quality DNA is approximately 1.8 to 2 at 260/280 nm. QUBIT 3.0 Fluorometer (ThermoFisher Scientific, USA) quantified the extracted DNA using dS DNA HS (High Sensitivity) Dye. One hundred nanograms of intact DNA was enzymatically fragmented using Covaris, targeting 200–300 bp fragment size. The DNA fragments undergo end repair, where the mix converts the overhangs resulting from fragmentation into blunt ends. The end-repair mix’s 3′ to 5′ exonuclease activity removes the 3′ overhangs, and polymerase activity fills in the 5′ overhangs. The blunt-ended fragments are adenylated, adding a single ‘A’ nucleotide to the 3′ ends. The adenylated fragments loop adapters were ligated and cleaved with a uracil-specific excision reagent (USER) enzyme. The AMPure beads (Beckman Coulter Life Sciences, USA) further purified the DNA. The polymerase chain reaction (PCR) with 6 cycles using NEBNext Ultra II Q5 master mix (New England BioLabs®Inc., USA), Illumina universal primer (Illumina, USA), and sample-specific octamer primers amplified the DNA. The AMPure beads cleaned the amplified DNA, and the final DNA library was eluted in 15 μL of 0.1X Tris EDTA buffer. QUBIT 3.0 Fluorometer quantified 1ul of the library using dS DNA HS reagent [65]. 1 μL of the library was loaded on an Agilent DNA 7500 chip, and fragments were analyzed by Agilent 2100 Bioanalyzer. Illumina HiSeq 4000 (Illumina, USA) sequenced the DNA with a 2 × 150 bp paired-end run.

FASTQC v0.11.9 and MultiQC v1.10.1 checked the quality of raw reads [66, 67]. Trimgalore v0.6.6 removed the adapter and other contaminations (https://github.com/FelixKrueger/TrimGalore). The integration of both dry lab (in silico) and wet lab (in vitro) studies, utilizing a variety of databases, software, and tools (as illustrated in the following flowchart) enhance accuracy and minimize biases [68].

figure a

Assembly and annotation of raw reads

The raw reads were assembled de novo (primary assembly) using Unicycler v0.4.9 [69]. The identity of isolate was further confirmed using 16S rRNA, and the average nucleotide identity (ANI) which was calculated from pair-wise comparisons of all sequences shared between the query and reference strain [70, 71]. The PubMLST web server (https://pubmlst.org/bigsdb?db=pubmlst_rmlst_seqdef_kiosk) identified the bacterial isolates [72]. Contiguator assembled the raw reads by mapping the contigs over a reference genome (secondary assembly) retrieved from the NCBI database, NZ_CP042858.1 for K. pneumoniae, NZ_CP053706.1, and NZ_LN831024.1 for P. aeruginosa [73]. PROKKA (Prokaryotic Genome Annotation) v1.12 annotated all the genomic features of the de novo assembled genome [74]. Mob-suite v3.0.3 distinguished and partitioned plasmid and nuclear sequence [75]. Plasmid types were further identified using PlasmidFinder2.1 (https://cge.food.dtu.dk/services/PlasmidFinder/) [76, 77].

MLST, capsular, and antigen typing

The multi-locus sequence typing of K. pneumoniae for seven housekeeping genes rpoB, gapA, mdh, pgi, phoE, infB, and tonB, and P. aeruginosa for acsA, aroE, guaA, mutL, nuoD, ppsA and trpE was performed using the command-line version of mlst (https://github.com/tseemann/mlst) against PubMLST typing schemes based on Bacterial Isolate Genome Sequence database (BIGSdb) [72]. Kleborate screened the K. pneumoniae strains and predicted the K (capsule) and O antigen (LPS) serotypes [78]. The Pseudomonas aeruginosa serotyper (PAst) v1.0 (https://cge.food.dtu.dk/services/PAst/) identified the serogroup of P. aeruginosa by running a BLAST search on the O-specific antigen (OSA) gene cluster [77, 79].

SNP and Indel detection

Snippy v4.6.0 (https://github.com/tseemann/snippy) identified the SNPs and indels conferring colistin resistance by mapping against colistin-susceptible reference whole genomes extracted from the NCBI database (https://www.ncbi.nlm.nih.gov/), NC_009648.1 K. pneumoniae subsp. pneumoniae MGH 78578 (ATCC® 700721™) and NZ_CP017149.1 P. aeruginosa strain (ATCC® 15692™). The SNPs rendering deleterious amino acid changes were assigned by SIFT (Sorting Intolerant From Tolerant) [80].

ARGs identification

Resistance Gene Identifier (RGI) v6.0.0 identified the antibiotic resistance genes with 80% similarity and 90% identity in the isolates against using the Comprehensive Antibiotic Resistance Database (CARD) [81].

Plasmids and IS elements

PlasmidFinder v2.1(https://cge.food.dtu.dk/services/PlasmidFinder/) identified the plasmids and mobile genetic elements were identified by Mobile Element Finder v1.0.3 (https://cge.food.dtu.dk/services/MobileElementFinder/).

Virulence factors and pathogenicity toward humans

Virulence factors were identified by ABRicate v1.0.1 (https://github.com/tseemann/abricate) using the Virulence Factor Database (VFDB) with a cut-off coverage of 80 and 70% identification. The pathogenicity towards humans was predicted by PathogenFinder v1.1 (https://cge.food.dtu.dk/services/PathogenFinder/).

CRISPRs, bacteriophages, and genomic islands interspersed in the bacterial genome

CRISPR (Clustered regularly interspaced short palindromic repeat) arrays and Cas proteins were identified by CRISPRCasFinder (https://crisprcas.i2bc.paris-saclay.fr/CrisprCasFinder/Index) [82]. Phigaro v2.2.6 detected and annotated prophage sequences within the bacterial genome. Genomic islands were individually identified using IslandViewer 4 (https://www.pathogenomics.sfu.ca/islandviewer/) and compared using IslandCompare (https://islandcompare.ca/).

Phylogenetic analyses

Phylogenetic relationship was obtained using Roary v3.13.0 (https://github.com/sanger-pathogens/Roary) with a blastP identity cut-off criterion of 90%. The genes present in more than 99% of the isolates were assigned as core genes. Gubbins v 3.3 removed the recombination events (https://github.com/nickjcroucher/gubbins).

Comparative genome analyses

The genome of K. pneumoniae strains and P. aeruginosa strains were compared against NC_009648.1 K. pneumoniae subsp. pneumoniae MGH 78578 (ATCC® 700,721™) and NZ_CP017149.1 P. aeruginosa strain (ATCC® 15,692™) respectively using Proksee (https://proksee.ca/).

The figures were generated in R using cowplot, ggplot2, pheatmap and Easyfig v.2.1 (https://mjsull.github.io/Easyfig/).

Statistical analyses

The chi-square test was conducted to analyze the results using R v.4.1.2. The significance level was p < 0.05.

Results

The primary objective of our study was to explore colistin resistance in XDR and PDR clinical isolates obtained from sepsis and tuberculosis patients in an Indian hospital. Notable observations include the discordance in colistin resistance detection between Vitek®2 and broth microdilution methods. Whole-genome sequencing unveiled a substantial prevalence of colistin-resistant strains, with 17 demonstrating extensive drug resistance and 3 displaying pan-drug resistance. Detailed analyses encompassing genomic features, plasmid sequences, molecular resistance determinants, virulence factors distribution, and CRISPR/Cas systems shed light on the underlying mechanisms of colistin resistance in K. pneumoniae and P. aeruginosa strains.

Phenotypic characteristics of the recovered isolates

Clinical data and phenotypic characterisation of colistin-resistant bacterial isolates

We investigated colistin resistance in the isolates retrieved from blood culture of sepsis patients, pleural fluid, and tracheal aspirate culture of patients suffering from tuberculosis in Jawaharlal Nehru Medical College and Hospital (J.N.M.CH.), Aligarh, Uttar Pradesh, India, during the year 2020–2021. During this period, 1784 blood, tracheal aspirate, and pleural fluid samples were collected from Tuberculosis and Respiratory Disease centre (TBRD), Intensive Care Unit (ICU), and Neonatal Intensive Care Unit (NICU) for antibiotic susceptibility tests (AST), out of which 214 samples showed positive growth. AST through automated Vitek®2, detected only six colistin-resistant strains (MIC = 4 μg/mL). AST, through the broth microdilution method according to the latest CLSI guidelines, demonstrated 48 strains as colistin-resistant, out of which 17 strains were identified as extremely drug-resistant (XDR) and three as pan-drug resistant (PDR) [60]. Colistin resistance was high for these 20 strains, with a MIC range between 32 and 128 μg/mL (Additional Table 1). According to MAR index, K. pneumoniae ST147 isolates (MAR index = 0.846) were resistant to approximately 84.6% of the antibiotics included in the study. P. aeruginosa ST235 (MAR index = 1.0) exhibited resistance against all the tested antibiotics whereas P. aeruginosa ST357 (MAR index = 0.923) displayed resistance against approximately 92.3% of the tested antibiotics [61]. PCR amplification for detecting mcr-1 was negative for all the isolates.

Genomic features identical in colistin-resistant strains

The whole genome of 20 XDR strains was sequenced for analysis. The GC content was comparable across all the strains in each ST lineage (Additional Table 2). All the sequenced strains were identified as AMR-high-risk clones, 16 as K. pneumoniae ST 147 (Additional Table 3), three as P. aeruginosa ST235 (AK-624, 625, 628), and one as P. aeruginosa ST357 (AK-631). All the strains were predicted as human pathogens. The K. pneumoniae ST147 belongs to clonal group CG147, which has disseminated globally, causing several outbreaks in recent years. The identified capsular types were KL10, wzi420, and O3/O3a O locus in all the K. pneumoniae strains. None of the strains possessed hypervirulence genes rmpA and rmpA2. P. aeruginosa ST235 is the founder of the CC235 clonal complex, and P. aeruginosa ST357 is the founder of the CC357 clonal complex. All the P. aeruginosa strains displayed O11 serogroup (Additional Table 3).

Identification of nuclear and plasmid sequences

Eight plasmids were reconstructed in K. pneumoniae strains, further identified as IncFIB, IncFIIA, IncR, IncFIA, and ColRNAI. In addition to the plasmids mentioned above, PlasmidFinder-2.1 identified Col(phAD28) in AK-618, 620, 626, and 630 (Additional Table 4). In P. aeruginosa strains (AK-624, 625, 628 and 631), the reconstructed plasmid was identified as IncP by Mob-suite v3.0.3. PlasmidFinder-2.1 didn’t detect any plasmid in P. aeruginosa strains (Additional Table 4).

Molecular determinants of colistin resistance indicate mgrB deletion and mutations in chromosomally mediated factors

Although no mcr genes were detected, several factors inducing colistin resistance were identified in K. pneumoniae and P. aeruginosa strains. The pmr phosphoethanolamine transferase encoding genes arnT, eptB, ompA bacterial porin with reduced permeability to peptide antibiotics were detected in K. pneumoniae (Fig. 1A). None of the K. pneumoniae strains possessed mgrB, indicating deletion. The SNP/indel detection against susceptible strain confirmed the complete deletion of mgrB due to bidirectional gene fusion.

Fig. 1
figure 1

A Antibiotic resistance genes in K. pneumoniae strains. The phylogenetic relatedness among the strains was determined by recombination free alignment of single-nucleotide variants maximum-likelihood phylogenetic tree rooted at K. pneumoniae MGH 78578 (ATCC® 700,721™). All strains possessed KL10 (wzi420) capsular type and O3/O3a antigen. B Genetic characterization of blaNDM-5. Bleomycin resistance gene BRP(MBL), Phosphoribosylanthranilate (PRA) isomerase (TrpF) and disulfide isomerase-like protein (DsbD) were flanked downstream of blaNDM-5

There were no mutations detected in the PhoPQ two-component system. Mutations at two amino acid positions in pmrB (T246A and R256G) and four amino acids in eptA (V50L, A135P, I138V, and C27F) were found in all the K. pneumoniae isolates. Missense mutations in arnA, arnB, arnT, and kpnF were also detected (Additional Table 5A). Frameshift mutations in ompK35 increase the carbapenem resistance [83].

Five pmr phosphoethanolamine transferase genes basS (also known as pmrB), basR, arnA, cprR, and cprS were identified in P. aeruginosa strains (Fig. 2A). These genes modulate the addition of phosphoethanolamine to lipid A, creating a cationic charge on the cell membrane inhibiting the binding of positively charged colistin. Mutations in arnA, arnB, arnC, arnD, arnF, arnT, pmrA, pmrB, and parS were identified in P. aeruginosa (Additional Table 5B and C).

Fig. 2
figure 2

A Antibiotic resistance genes in P. aeruginosa strains. The phylogenetic relatedness among the strains was determined by recombination free alignment of single-nucleotide variants maximum-likelihood phylogenetic tree rooted at P. aeruginosa strain (ATCC® 15,692™). Three strains were ST235 and one was ST357. All strains possessed O11 serogroup. B Genetic characterization of blaNDM-1. No adjacent genes were identified in the contig carrying blaNDM-1, owing to limited sequencing depth

Distribution of virulence factors and antimicrobial resistance genes (ARGs) among the recovered isolates

Prevalence of ARGs against all the major antibiotic classes

All the strains were detected as extended-spectrum-lactamase (ESBL)- and carbapenem-resistant strains. The resistance profile and total number of ARGs in all the K. pneumoniae strains were comparatively similar, conferring resistance against approximately all the major classes of antibiotics, including aminocoumarin, aminoglycoside, carbapenem, cephalosporin, macrolide, monobactam, peptide and sulfonamide (Fig. 1A and Additional Table 7). The predominant resistance mechanisms were antibiotic efflux and antibiotic inactivation. All the K. pneumoniae strains carried metallo-β-lactamase blaNDM-5, extended-spectrum β-lactamases blaSHV-11, blaCTX-M-15 and carbapenemases blaOXA-1 and blaOXA-232. The blaNDM-5 varies from blaNDM-1 at two positions, i.e., V88L and M154L, and is endemic to India [37, 39]. blaNDM-5 was mediated by the IncFII plasmid (Fig. 1B). K. pneumoniae ST147 is often characterized by gyrA S831 and parC S80I QRDR mutations conferring quinolone resistance and blaCTX-M-15 mediated by ISEcp1. K. pneumoniae strains in this study possessed gyrA D87A and parC with S80I and N304S mutation. blaCTX-M-15 was mediated by the IncR plasmid.

The resistome of P. aeruginosa strains was distinct from that of K. pneumoniae. The total number of ARGs was comparatively higher (Fig. 2A). In P. aeruginosa ST235, we detected metallo-β-lactamase blaVIM-2 and carbapenemase blaOXA-488. In P. aeruginosa ST357, blaNDM-1, blaOXA-10 and blaOXA-846 were identified. blaNDM-1 was chromosomally mediated, inferred through its presence in the nuclear partition of the genome, and further confirmed with PlasmidFinder2.1 (Fig. 2B).

Virulence factors (V.F.S)

All K. pneumoniae strains possessed 19 virulence factors, notably, entB and entA encoding siderophores system enterobactin, ompA, pilus chaperones encoding genes yagZ/ecpA, yagY/ecpB, yagV/ecpE and yagX/ecpC, polymerized tip adhesin yagW/ecpD, iron-enterobactin ABC transporter permease fepG, and ferrienterobactin ABC transporter ATPase fepC (Additional Table 6).

P. aeruginosa AK-624, AK-625, and AK-628 had an extensive array of 224 virulence factors contributing to pathogenesis (Additional Table 6). AK-631 carried 216 virulence factors. Alginate biosynthesis proteins encoding genes (alg44, alg8, algA, algB, algC, algD, algE, algF, algG, algI, algJ, algK, algL, algP/algR3, algQ, algR, algU, algW, algX, algZ), type III, IV and VI secretion system associated virulence factors were prevalent. Exoenzymes Exo T, ExoU, and ExoY of the type III secretion system were present. ExoU is a significant virulence factor causing acute epithelial injury, rapid cell lysis, and destruction of the host cell plasma membrane [46].

IS elements associated with ARGs

IS5075, ISEc29, ISKpn1, and ISKpn14 were common in all the ColRKp (colistin-resistant K. pneumoniae) strains (Additional Table 4). The beta-lactamase blaTEM-1, aminoglycoside resistance genes aph(6)-Id and aph(3″)-Ib, and sulfonamide resistance gene sul2 were present adjacent to each other, flanked by IS5075. Aminoglycoside resistance gene armA and macrolide resistance genes msrE and mphE were flanked by ISEc29. In P. aeruginosa, ISPa6, ISPa86, IS6100 and ISpa7 were identified. In AK-631, the phosphoethanolamine transferases cprR and cprS were flanked by ISPre2, aiding in the HGT of these genes conferring colistin resistance. Only the PmpM gene encoding for fluoroquinolone was flanked by ISPa86. Tn6196 was associated with aph(3′)-VI, ISPa6 with pmpM, phosphonic acid fosA, and phenicol resistance gene mexM and mexN, and IS26 with fluoroquinolone resistance gene qnrVC1.

CRISPR/Cas systems

All the K. pneumoniae strains were accompanied by CRISPR/Cas Class I with 43 spacers. Class I CRISPR with 9, 17, and 7 spacers were detected in P. aeruginosa ST357 (Additional Table 3).

Genomic Islands (GIs)

K. pneumoniae genomes screened against K. pneumoniae subsp. pneumoniae MGH 78578 predicted common genomic islands mostly consisting of hypothetical proteins with unidentified functions (694), and protein-encoding genes were 254. Tn3 family transposase genes and prophage integrase IntA and IntS were represented in the GIs. GIs were comparable across K. pneumoniae strains. AK-624 had GIs comprised of 339 hypothetical proteins and 113 protein-encoding genes. IS1182 family transposase ISCfr1, IS1182 family transposase ISPa7, IS3 family transposase ISPa86, IS3 family transposase ISPosp5, ISL3 family transposase ISStma11, and ISL3 family transposase ISStma11 were present. The distribution of GIs was comparable in AK-625 and AK-631 but significantly varied in AK-628 with 5786 GIs, out of which 3446 were encoded for proteins. A comparative analysis revealed common GIs in K. pneumoniae and P. aeruginosa strains, respectively (Fig. 3A and B).

Fig. 3
figure 3

A Genomic islands in K. pneumoniae strains. GIs screened against K. pneumoniae subsp. pneumoniae MGH 78578 predicted common genomic islands mostly consisting of hypothetical proteins with unidentified functions (694), and protein-encoding genes were 254. Tn3 family transposase genes and prophage integrase IntA and IntS were represented in the GIs. B Genomic islands in P. aeruginosa strains. P. aeruginosa ST235 carried GIs comprised of 339 hypothetical proteins and 113 protein-encoding genes comprising of IS1182 family transposase ISCfr1, IS1182 family transposase ISPa7, IS3 family transposase ISPa86, IS3 family transposase ISPosp5, ISL3 family transposase ISStma11, and ISL3 family transposase ISStma11 were present. P. aeruginosa ST357 was distinct with 5786 GIs, out of which 3446 were encoded for proteins

Prophage sequences within the bacterial genome

Sequences of myoviridae and siphoviridae were incorporated into the genome of all the strains.

Phylogenetic analyses

The phylogenetic relationship was explored based on concatenating 4528 core genes for K. pneumoniae strains and 5376 core genes for P. aeruginosa strains. The final recombination-free alignment of 21,687 and 95,729 single-nucleotide variants (SNVs) for K. pneumoniae and P. aeruginosa, respectively, were used to construct a maximum-likelihood phylogenetic tree using Interactive Tree of Life iTOL v 6.7.3 (https://itol.embl.de/). The tree was rooted with NC_009648.1 K. pneumoniae subsp. pneumoniae MGH 78578 (ATCC® 700,721™) and NZ_CP017149.1 P. aeruginosa strain (ATCC® 15,692™) (Figs. 1 and 2).

Discussion

The global spread of multi-drug resistant pathogenic bacterial species has become detrimental to the available antimicrobial treatment regimen. The emergence of beta-lactamases such as blaNDM, extended spectrum beta lactamases (ESBLs) such as blaCTX-M, blaSHV, and blaTEM etc. has led to the use of last resort antibiotics such as tigecycline, polymyxin E, daptomycin, vancomycin and linezolid [84,85,86]. However, in the last decade several ARGs against the last-line antibiotics have been reported [87,88,89,90,91,92,93]. Colistin is one of the last-resort antibiotics which is administered in cases of MDR Gram-negative bacterial infections. Colistin resistance is a huge challenge that may not be as prevalent as beta-lactamases such as blaNDM at present but its emergence and dissemination, especially in hospital settings, is a premonition of impending outbreaks [94]. The inefficiency of standard disk diffusion methods and automated ASTs such as Vitek2 for accurate detection of colistin resistance unified with limited resources and expertise for broth micro-dilution at most hospitals in LMICs such as India demonstrates a low presence of colistin resistance in the country [55]. Next-generation sequencing is an efficient, non-PCR-biased technology that allows for the accurate surveillance of antimicrobial resistance genes and the exploration of their underlying mechanisms. Hospitals are often a hotspot for nosocomial infections and outbreaks. In this study, we analyzed clinical isolates obtained from blood, tracheal aspirate, and pleural fluid of patients admitted to an Indian hospital to assess the prevalence of XDR or PDR bacterial strains exhibiting colistin resistance. While Vitek2 identified only six colistin-resistant isolates, antimicrobial susceptibility testing (AST) using broth micro-dilution revealed 48 isolates as colistin-resistant. Among them, 17/48 (16 K. pneumoniae and 1 P. aeruginosa) isolates were identified as XDR, and 3/48 (all P. aeruginosa) were PDR. Notably, all isolates tested negative for the mcr-1 gene. Subsequently, whole-genome sequencing and pan-genomic analyses were conducted to mitigate PCR biases and facilitate a comprehensive exploration of the underlying mechanisms driving colistin resistance. Interestingly, although the strains were isolated from different patients and wards in a hospital (Additional Table 3), all the K. pneumoniae strains and P. aeruginosa strains displayed remarkable genomic similarities in each ST lineage (Fig. 4A and B). The identical clonal types, capsular type, resistance profile, and plasmids in colistin-resistant K. pneumoniae signal a potential XDR K. pneumoniae ST147 outbreak in hospital settings (Additional Table 3). The same goes for circulating colistin-resistant P. aeruginosa ST235 strains. K. pneumoniae ST147 and P. aeruginosa ST235 and ST357 are global high-risk clones carrying multiple ARGs [9, 46, 95].

Fig. 4
figure 4

Comparison of whole genome against colistin susceptible strains. A Genomes of K. pneumoniae strains were compared against colistin susceptible ATCC strain, NC_009648.1 K. pneumoniae subsp. pneumoniae MGH 78578 (ATCC® 700,721™). B Genomes of P. aeruginosa strains were compared against colistin susceptible ATCC strain NZ_CP017149.1 P. aeruginosa (ATCC® 15,692™)

Earlier, K. pneumoniae ST147 ColRKp displaying colistin resistance devoid of mcr genes has been reported in India. Mutation, truncation, or complete deletion of mgrB, a mutation in pmrB, phoQ, and insertional inactivation have been identified as significant factors inducing colistin resistance [20, 41, 42, 96,97,98]. Colistin displays bactericidal activity in gram-negative bacteria by binding to the lipid A of LPS. Colistin competitively destabilizes the outer membrane’s surface charge by displacing calcium (Ca2+) and magnesium (Mg2+) divalent cations. Consequently, the three-dimensional structure of LPS is impaired, and permeability is increased, allowing colistin to insert its hydrophobic terminal acyl fat chain. The chromosomally mediated mgrB gene is a crucial factor that modulates the surface charge stability to facilitate colistin susceptibility. mgrB negatively regulates the functioning of PhoPQ, a two-component system involved in LPS modification. The deletion of mgrB in all ColRKp might be the major factor driving colistin resistance (Additional Table 5) [16]. Without mgrB, PhoPQ actively adds the positively charged LAra4N to lipid A, increasing the cationic charge on the surface and prohibiting colistin binding. The phosphoethanolamine transferases add pEtN to lipid A, inhibiting colistin. We detected multiple phosphoethanolamine transferases, namely, arnT and eptB in K. pneumoniae strains and basS, basR, arnA, cprR, and cprS in P. aeruginosa (Figs. 1 and 2). All the strains exhibited amino acid substitutions T246A and R256G in pmrB (Additional Table 5). Although earlier studies have also reported these mutations in colistin-resistant strains, they are identified as neutral mutations and believed to have no direct correlation with colistin resistance [99]. As the colistin resistance mechanism is complex, these predictions are not absolutely conclusive, and more investigations, especially at the transcriptional and translational level, are required to understand their relationship with colistin resistance. lpxM adds a secondary acyl chain to lipid A, and mutation in lpxM can reduce colistin susceptibility in K. pneumoniae [18]. The identified mutations in eptA, arnA, arnB, arnT, and kpnF in K. pneumoniae have been reported to increase colistin resistance [100]. The mutations in arnA, arnB, arnC, arnD, arnF, and arnT modify the LPS. The pmrB and parS mutations in P. aeruginosa also increase colistin resistance. His398Arg mutation in parS is deleterious, leading to alteration in the functional domain of histidine kinase-like ATPase, which often develops under external colistin stress and ultimately causes colistin resistance [21]. The complete and accurate mechanism of colistin resistance remains unknown. Although no mcr gene was detected, the cumulative effect of all these factors reduces colistin susceptibility.

The prevalence of carbapenemases (blaNDM-5, blaOXA-1, blaOXA-232, blaSHV-11, blaCTX-M-15 and LptD), aminoglycoside resistance genes (aadA2, acrD, APH(3″)-Ib, APH(6)-Id, armA, and baeR), and sulfonamide resistance genes (sul1, sul2) along with ARGs against other major antibiotic classes in K. pneumoniae strains poses a severe threat (Fig. 1A). The porins ompK35, ompK36, ompC, and ompA mutations reduce susceptibility towards carbapenem [83]. ompA mutation is also involved in colistin resistance. The association of beta-lactamase blaTEM-1, aminoglycosides APH(6)-Id and APH(3″)-Ib and armA, macrolide resistance genes msrE and mphE, and sulfonamide resistance gene sul2 with IS elements in K. pneumoniae strains may expedite the horizontal gene transmission. The phosphoethanolamine transferases cprR and cprS flanked by ISPre2 in P. aeruginosa can disseminate colistin resistance-carrying genes horizontally. The resistance of all these clinical strains towards colistin and approximately every other class of antibiotics in clinical settings poses significant challenges due to colistin’s critical status as a last-resort antibiotic. These bacteria’s adaptability and rapid spread of resistance determinants worsen the global healthcare burden [92, 101].

The chi-square test (p = 0.007) revealed a significant association between K. pneumoniae and P. aeruginosa isolates and their prevalence rates. However, an undefined odds ratio was observed due to a lack of representation in the comparison group. Fisher’s Exact Test yielded a p-value of 1, indicating no significant association between resistance genes and phenotypic resistance (Additional Table 7B and C). Despite the absence of strains without corresponding resistance genes, further research is needed to fully understand this relationship and develop effective interventions against antimicrobial resistance.

Virulence factors in K. pneumoniae and P. aeruginosa enhance pathogenicity, attachement to the host cell and colonization, protection against host immune response and biofilm production [48, 53, 102,103,104]. Globally, K. pneumoniae ST147 comprises two major lineages, one with KL10 (wzi420) capsular type and O3 or O3a locus, and another with KL64 (wzi64) capsular type and O2v1 locus [36] Unlike previously reported hospital outbreaks of K. pneumoniae ST147 with KL64 capsular typing in Tunisia and Poland, in our study, KL10 (wzi420) was associated [17, 34] (Additional Table 3). KL10 capsulate serotype alongwith V.F.s encoding for Type 1 fimbriae, pili, siderophores and ompA increases the K. pneumoniae adherence to the host cell and biofilm formation. Biofilm formation protects the underlying bacterial colonies against exogenous stressors and antimicrobial agents [49, 50, 105]. Even though the hypervirulence genes rmpA/rmpA2 were absent, the virulence factors increasing adherence and penetration were common in all K. pneumoniae strains (Additional Table 6).

In P. aeruginosa, exoU is a critical cytotoxic virulence factor that damages the host cells and is involved in mortality due to P. aeruginosa infections [106]. The cascade of alginate proteins in P. aeruginosa facilitates biofilm formation, protecting bacteria from antibiotics and immune responses, and hence aiding chronic infections [49]. Further, the presence of Type 4 pilus proteins, type 4 fimbrieae, and twitching motility proteins contributes to its infectivity. These structures play crucial roles in the pathogen’s ability to adhere to and invade host cells, thereby facilitating the initiation and progression of infection. Type 4 pilus proteins enable P. aeruginosa to firmly attach to host tissues, promoting colonization and the formation of biofilms. The twitching motility proteins allow P. aeruginosa to actively move across surfaces and penetrate deeper into host tissues, facilitating the spread of infection and exacerbating disease severity. This expanded repertoire of virulence factors underscores the formidable pathogenicity of P. aeruginosa and its ability to cause various infections, particularly in immunocompromised individuals or those with underlying health conditions. The Type VI Secretion System (T6SS) is a molecular weapon used by P. aeruginosa and K. pneumoniae, to deliver toxic effector proteins directly into neighboring cells. In P. aeruginosa, the T6SS is essential for virulence, interbacterial interactions, and outcompeting other bacteria by delivering toxins into neighboring cells, thereby killing or inhibiting their growth. This competitive advantage is particularly important in environments such as the human body, where P. aeruginosa often encounters other microbial species [48, 49, 52, 53]. (Additional Table 6). Biofilm-associated infections in healthcare present challenges due to their resilience, leading to chronicity and increased costs. Biofilms on implanted devices often require removal, while chronic wounds, like diabetic ulcers, face healing delays and complications [107]. In conditions like cystic fibrosis, biofilm-forming P. aeruginosa worsen lung function, emphasizing the severity of associated respiratory infections [44, 50]. CRISPR regions in K. pneumoniae and P. aeruginosa unveiled genomic diversity and assessed immunity against foreign invaders. Presence of CRISPR I enhances survival.

As we explored a limited number of samples from one hospital, it doesn’t give a definite prediction regarding the prevalence of colistin resistance throughout the country. Still, it is a problematic snapshot, with approximately 22% (48/214) colistin-resistant strains. It is in congruence with colistin prevalence in Asia [108]. Our study demonstrates a remarkable similarity among the genomes of colistin-resistant K. pneumoniae and P. aeruginosa strains, respectively, indicating a potential outbreak from the same source of origin (Figs. 14). Due to a lack of real-time surveillance strategies, these outbreaks often go unnoticed in low-middle-income countries. Plasmid-associated ARGs, ARGs flanked by IS elements, multiple virulence factors, genomic islands, and prophage sequences incorporated in bacterial genome ensure a high rate of horizontal gene transfer. Co-resistance of colistin and carbapenem resistance poses a colossal challenge, leaving no alternative antibiotic treatment regimen to the rescue.

Limitations of the study

The primary limitation of this study lies in its reliance on samples collected solely from a single hospital setting. Consequently, the findings may not fully capture the overall prevalence of colistin resistance and the dissemination of XDR and PDR strains across wider geographic regions or diverse patient populations. While the study coincidentally observed XDR and PDR colistin-resistant strains predominantly among K. pneumoniae and P. aeruginosa isolates, the exclusive focus on these pathogens inadvertently overlooks other priority pathogens, such as Acinetobacter baumannii, which also possess significant resistance capabilities. Consequently, the generalizability of the study’s findings to other healthcare settings or regions may be constrained, necessitating further research encompassing a broader spectrum of pathogens and healthcare contexts for a comprehensive understanding of colistin resistance dissemination.

Strength of the study

Despite its limitations, this study provides valuable insights into the prevalence and genomic characteristics of colistin-resistant strains, especially among K. pneumoniae and P. aeruginosa in a hospital setting. Comprehensive genomic analyses, including whole-genome sequencing and pan-genomic analyses, enhance our understanding of colistin resistance dissemination and identify MDR, XDR and PDR strains. Utilizing broth microdilution for antimicrobial susceptibility testing yields more accurate results than Vitek2 systems, improving surveillance and management of antimicrobial resistance. Additionally, despite the absence of the mobile colistin resistance (mcr) gene, SNP analysis through whole-genome sequencing reveals novel genetic mechanisms underlying colistin resistance, emphasizing the importance of advanced genomic techniques in combating antimicrobial resistance.

Conclusions

The global dissemination of MDR pathogenic bacterial species presents a critical challenge to current antimicrobial treatment regimens. The emergence of beta-lactamases such as blaNDM and extended-spectrum beta-lactamases (ESBLs) like blaCTX-M, blaSHV, and blaTEM has necessitated the utilization of last-resort antibiotics such as colistin. However, the rise of resistance against these last-line antibiotics, exemplified by colistin, poses a significant threat, particularly in healthcare settings. The inefficiency of standard detection methods for colistin resistance, coupled with limited resources for comprehensive surveillance, underscores the urgency for improved detection strategies. Through whole-genome sequencing and pan-genomic analyses, we identified a concerning prevalence of colistin-resistant strains, particularly among K. pneumoniae and P. aeruginosa isolates from a single hospital setting in India. Notably, despite originating from different patients and wards, the genomic similarity among colistin-resistant strains suggests a potential outbreak scenario, emphasizing the importance of real-time surveillance and infection control measures. The detection of extensively drug-resistant (XDR) and pandrug-resistant (PDR) strains among colistin-resistant isolates is concerning. XDR and PDR bacteria resist multiple antibiotic classes, severely limiting treatment options and increasing treatment failure risks. Failure to address XDR and PDR spread could profoundly impact patient care, healthcare resources, and public health. The identification of diverse resistance mechanisms, including mutations in key regulatory genes and plasmid-mediated transfer of resistance determinants, underscores the complexity of combating colistin resistance. Addressing the challenges posed by colistin-resistant MDR, XDR and PDR strains requires a multifaceted approach, encompassing enhanced surveillance, infection control strategies, and the development of novel therapeutic interventions.

Availability of data and materials

Data is provided within the manuscript or supplementary information files.

The whole genome sequencing data are submitted to the National Center for Biotechnology Information (NCBI) database with BioProject Accession number PRJNA905193.

Abbreviations

ARG:

Antibiotic resistance genes

AMR:

Antimicrobial resistance

ST:

Sequence type

MLST:

Multi locus sequence typing

ColRKp:

Colistin resistant Klebsiella pneumoniae

References

  1. Walesch S, Birkelbach J, Jézéquel G, Haeckl FPJ, Hegemann JD, Hesterkamp T, et al. Fighting antibiotic resistance—strategies and (pre) clinical developments to find new antibacterials. EMBO Rep. 2023;24.

  2. Murray CJ, Ikuta KS, Sharara F, Swetschinski L, Robles Aguilar G, Gray A, et al. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet. 2022;399:629–55.

    Article  CAS  Google Scholar 

  3. Algammal AM, Wahdan A, Elhaig MM. Potential efficiency of conventional and advanced approaches used to detect Mycobacterium bovis in cattle. Microb Pathog. 2019;134.

  4. Elbehiry A, Marzouk E, Aldubaib M, Moussa I, Abalkhail A, Ibrahem M, et al. Pseudomonas species prevalence, protein analysis, and antibiotic resistance: an evolving public health challenge. AMB Express. 2022;12.

  5. Algammal AM, Hashem MEA, Alfifi KJ, Al-Otaibi AS, Alatawy M, Eltarabili RM, et al. Sequence analysis, Antibiogram profile, virulence and antibiotic resistance genes of XDR and MDR Gallibacterium anatis isolated from layer chickens in Egypt. Infect Drug Resist. 2022;15:4321–34.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Algammal AM, Ibrahim RA, Alfifi KJ, Ghabban H, Alghamdi S, Kabrah A, et al. A first report of molecular typing, virulence traits, and phenotypic and genotypic resistance patterns of newly emerging XDR and MDR Aeromonas veronii in Mugil seheli. Pathogens. 2022;11.

  7. Kareem SM, Al-Kadmy IMS, Kazaal SS, Ali ANM, Aziz SN, Makharita RR, et al. Detection of gyrA and parC mutations and prevalence of plasmid-mediated quinolone resistance genes in Klebsiella pneumoniae. Infect Drug Resist. 2021;14:555–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Shafiq M, Zeng M, Permana B, Bilal H, Huang J, Yao F, et al. Coexistence of blaNDM–5 and tet(X4) in international high-risk Escherichia coli clone ST648 of human origin in China. Front Microbiol. 2022;13.

  9. Peirano G, Chen L, Kreiswirth BN, Pitout JDD. Emerging Antimicrobial-Resistant High-Risk Klebsiella pneumoniae Clones ST307 and ST147. Antimicrob Agents Chemother. 2020;64(10):e01148–20. https://doi.org/10.1128/AAC.01148-20.

  10. Liu YY, Wang Y, Walsh TR, Yi LX, Zhang R, Spencer J, et al. Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect Dis. 2016;16:161–8.

    Article  PubMed  Google Scholar 

  11. Gogry FA, Siddiqui MT, Sultan I, Haq QMR. Current update on intrinsic and acquired Colistin resistance mechanisms in Bacteria. Front Med. 2021;8.

  12. Haeili M, Javani A, Moradi J, Jafari Z, Feizabadi MM, Babaei E. MgrB alterations mediate colistin resistance in Klebsiella pneumoniae isolates from Iran. Front Microbiol. 2017;8.

  13. Yap PSX, Cheng WH, Chang SK, Lim SHE, Lai KS. MgrB mutations and altered cell permeability in Colistin resistance in Klebsiella pneumoniae. Cells. 2022;11.

  14. Giordano C, Barnini S, Tsioutis C, Chlebowicz MA, Scoulica EV, Gikas A, et al. Expansion of KPC-producing Klebsiella pneumoniae with various mgrB mutations giving rise to colistin resistance: the role of ISL3 on plasmids. Int J Antimicrob Agents. 2018;51:260–5.

    Article  CAS  PubMed  Google Scholar 

  15. Kocsis B, Kádár B, Tóth Á, Fullár A, Szabó D. MgrB variants in colistin-susceptible and colistin-resistant Klebsiella pneumoniae ST258. J Microbiol Immunol Infect. 2017;50:735–6.

    Article  CAS  PubMed  Google Scholar 

  16. Song Y, Yang J, Ju Y, Liu D, Wang Y, Sun X, et al. Genomic insights into a complete deletion of the mgrB locus in colistin-resistant Klebsiella pneumoniae ST2268 isolated from a human infection. J Glob Antimicrob Resist. 2021;27:75–8.

    Article  PubMed  Google Scholar 

  17. Jaidane N, Bonnin RA, Mansour W, Girlich D, Creton E, Cotellon G, et al. Genomic insights into colistin-resistant Klebsiella pneumoniae from a tunisian teaching hospital. Antimicrob Agents Chemother. 2020;64(10):e01148-20. https://doi.org/10.1128/AAC.01148-20.

  18. Halaby T, Kucukkose E, Janssen AB, Rogers MRC, Doorduijn DJ, van der Zanden AGM, et al. Genomic characterization of colistin heteroresistance in Klebsiella pneumoniae during a nosocomial outbreak. Antimicrob Agents Chemother. 2016;60:6837–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ghafur A, Shankar C, GnanaSoundari P, Venkatesan M, Mani D, Thirunarayanan MA, et al. Detection of chromosomal and plasmid-mediated mechanisms of colistin resistance in Escherichia coli and Klebsiella pneumoniae from Indian food samples. J Glob Antimicrob Resist. 2019;16:48–52.

    Article  PubMed  Google Scholar 

  20. Nirwan PK, Chatterjee N, Panwar R, Dudeja M, Jaggi N. Mutations in two component system (PhoPQ and PmrAB) in colistin resistant Klebsiella pneumoniae from north Indian tertiary care hospital. J Antibiot. 2021;74:450–7.

    Article  CAS  Google Scholar 

  21. Lee JY, Chung ES, Na IY, Kim H, Shin D, Ko KS. Development of colistin resistance in pmrA-, phoP-, parR- and cprR-inactivated mutants of Pseudomonas aeruginosa. J Antimicrob Chemother. 2014;69:2966–71.

    Article  CAS  PubMed  Google Scholar 

  22. Tacconelli E, Carrara E, Savoldi A, Harbarth S, Mendelson M, Monnet DL, et al. Discovery, research, and development of new antibiotics: the WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect Dis. 2018;18:318–27.

    Article  PubMed  Google Scholar 

  23. Chen Q, Wang M, Han M, Xu L, Zhang H. Molecular basis of Klebsiella pneumoniae colonization in host. Microb Pathog. 2023;177.

  24. Karami-Zarandi M, Rahdar HA, Esmaeili H, Ranjbar R. Klebsiella pneumoniae : an update on antibiotic resistance mechanisms. Future Microbiol. 2023; https://doi.org/10.2217/FMB-2022-0097.

  25. le Guern R, Grandjean T, Stabler S, Bauduin M, Gosset P, Kipnis É, et al. Gut colonisation with multidrug-resistant Klebsiella pneumoniae worsens Pseudomonas aeruginosa lung infection. Nat Commun. 2023;14:78.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Xu Q, Zheng B, Shen P, Xiao Y. Protective efficacy of statins in patients with Klebsiella pneumoniae bloodstream infection. Front Cell Infect Microbiol. 2023;12.

  27. Karampatakis T, Tsergouli K, Behzadi P. Carbapenem-resistant Klebsiella pneumoniae: virulence factors, molecular epidemiology and latest updates in treatment options. Antibiotics (Basel). 2023;12.

  28. Behzadi P, García-Perdomo HA, Karpiński TM, Issakhanian L. Metallo-ß-lactamases: a review. Mol Biol Rep. 2020;47:6281–94.

    Article  CAS  PubMed  Google Scholar 

  29. Algammal A, Hetta HF, Mabrok M, Behzadi P. Editorial: emerging multidrug-resistant bacterial pathogens “superbugs”: a rising public health threat. Front Microbiol. 2023;14.

  30. Damjanova I, Tóth A, Pászti J, et al. Expansion and countrywide dissemination of ST11, ST15 and ST147 ciprofloxacin-resistant CTX-M-15-type beta-lactamase-producing Klebsiella pneumoniae epidemic clones in Hungary in 2005--the new 'MRSAs'?. J Antimicrob Chemother. 2008;62(5):978–85. https://doi.org/10.1093/jac/dkn287.

  31. Coelho A, Mirelis B, Alonso-Tarrés C, Larrosa MN, Miró E, Abad RC, et al. Detection of three stable genetic clones of CTX-M-15-producing Klebsiella pneumoniae in the Barcelona metropolitan area, Spain. J Antimicrob Chemother. 2009;64:862–4.

    Article  CAS  PubMed  Google Scholar 

  32. Rodrigues C, Desai S, Passet V, Gajjar D, Brisse S. Genomic evolution of the globally disseminated multidrug-resistant Klebsiella pneumoniae clonal group 147. Microb Genom. 2022;8.

  33. Avgoulea K, Pilato V di, Zarkotou O, Sennati S, Politi L, Cannatelli A, et al. Characterization of extensively drug-resistant or Pandrug-resistant sequence type 147 and 101 OXA-48-producing Klebsiella pneumoniae causing bloodstream infections in patients in an intensive care unit. Antimicrob Agents Chemother 2018;62.

  34. Biedrzycka M, Urbanowicz P, Guzek A, Brisse S, Gniadkowski M, Izdebski R. Dissemination of Klebsiella pneumoniae ST147 NDM-1 in Poland, 2015-19. J Antimicrob Chemother. 2021;76:2538–45.

    Article  CAS  PubMed  Google Scholar 

  35. Lapp Z, Crawford R, Miles-Jay A, Pirani A, Trick WE, Weinstein RA, et al. Regional spread of blaNDM-1-containing Klebsiella pneumoniae ST147 in post-acute care facilities. Clin Infect Dis. 2021;73:1431–9.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Di Pilato V, Henrici De Angelis L, Aiezza N, Baccani I, Niccolai C, Parisio EM, et al. Resistome and virulome accretion in an NDM-1-producing ST147 sublineage of Klebsiella pneumoniae associated with an outbreak in Tuscany, Italy: a genotypic and phenotypic characterisation. Lancet Microbe. 2022;3:e224–34.

    Article  PubMed  Google Scholar 

  37. Rojas LJ, Hujer AM, Rudin SD, Wright MS, Domitrovic TN, Marshall SH, et al. NDM-5 and OXA-181 beta-lactamases, a significant threat continues to spread in the Americas. Antimicrob Agents Chemother. 2017;61.

  38. Peirano G, Pillai DR, Pitondo-Silva A, Richardson D, Pitout JDD. The characteristics of NDM-producing Klebsiella pneumoniae from Canada. Diagn Microbiol Infect Dis. 2011;71:106–9.

    Article  CAS  PubMed  Google Scholar 

  39. Hornsey M, Phee L, Wareham DW. A novel variant, NDM-5, of the New Delhi Metallo-β-lactamase in a multidrug-resistant Escherichia coli ST648 isolate recovered from a patient in the United Kingdom. Antimicrob Agents Chemother. 2011;55:5952–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Bean DC, Agarwal A, Cherian BP, Wareham DW. Hypermucoviscous polymyxin-resistant Klebsiella pneumoniae from Kolkata, India: genomic and phenotypic analysis. J Glob Antimicrob Resist. 2019;17:1–2.

    Article  PubMed  Google Scholar 

  41. Dey S, Gaur M, Sahoo RK, Das A, Jain B, Pati S, et al. Genomic characterization of XDR Klebsiella pneumoniae ST147 co-resistant to carbapenem and colistin – the first report in India. J Glob Antimicrob Resist. 2020;22:54–6.

    Article  PubMed  Google Scholar 

  42. Pragasam AK, Shankar C, Veeraraghavan B, Biswas I, Nabarro LEB, Inbanathan FY, et al. Molecular mechanisms of Colistin resistance in Klebsiella pneumoniae causing bacteremia from India—a first report. Front Microbiol. 2016;7:2135.

    PubMed  Google Scholar 

  43. Wood SJ, Kuzel TM, Shafikhani SH. Pseudomonas aeruginosa: infections, animal modeling, and therapeutics. Cells. 2023;12.

  44. Beg AZ, Rashid F, Talat A, Haseen MA, Raza N, Akhtar K, et al. Functional amyloids in Pseudomonas aeruginosa are essential for the proteome modulation that leads to Pathoadaptation in pulmonary niches. Microbiol Spectr. 2023;11.

  45. Langendonk RF, Neill DR, Fothergill JL. The building blocks of antimicrobial resistance in Pseudomonas aeruginosa: implications for current resistance-breaking therapies. Front Cell Infect Microbiol. 2021;11.

  46. Treepong P, Kos VN, Guyeux C, Blanc DS, Bertrand X, Valot B, et al. Global emergence of the widespread Pseudomonas aeruginosa ST235 clone. Clin Microbiol Infect. 2018;24:258–66.

    Article  CAS  PubMed  Google Scholar 

  47. del Barrio-Tofiño E, López-Causapé C, Oliver A. Pseudomonas aeruginosa epidemic high-risk clones and their association with horizontally-acquired β-lactamases: 2020 update. Int J Antimicrob Agents. 2020;56.

  48. Karampatakis T, Tsergouli K, Behzadi P. Pan-genome plasticity and virulence factors: a natural treasure trove for Acinetobacter baumannii. Antibiotics (Basel). 2024;13:257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Behzadi P, Gajdács M, Pallós P, Ónodi B, Stájer A, Matusovits D, et al. Relationship between biofilm-formation, phenotypic virulence factors and antibiotic resistance in environmental Pseudomonas aeruginosa. Pathogens. 2022;11.

  50. Tuon FF, Dantas LR, Suss PH, Tasca Ribeiro VS. Pathogenesis of the Pseudomonas aeruginosa biofilm: a review. Pathogens. 2022;11.

  51. Li Y, Ni M. Regulation of biofilm formation in Klebsiella pneumoniae. Front Microbiol. 2023;14.

  52. Qin S, Xiao W, Zhou C, Pu Q, Deng X, Lan L, et al. Pseudomonas aeruginosa: pathogenesis, virulence factors, antibiotic resistance, interaction with host, technology advances and emerging therapeutics. Signal Transduct Target Ther. 2022;7:1–27.

    PubMed  PubMed Central  Google Scholar 

  53. Santajit S, Sookrung N, Indrawattana N. Quorum sensing in ESKAPE bugs: a target for combating antimicrobial resistance and bacterial virulence. Biology. 2022;11:1466.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Laxminarayan R, Chaudhury RR. Antibiotic resistance in India: drivers and opportunities for action. PLoS Med. 2016;13.

  55. Khurana S, Malhotra R, Mathur P. Evaluation of Vitek®2 performance for colistin susceptibility testing for gram-negative isolates. JAC Antimicrob Resist. 2020;2.

  56. Tripathi N, Sapra A. Gram Staining. StatPearls; 2023.

    Google Scholar 

  57. Pages Monteiro L, Von Allmen N, Friesen I, Huth K, Zambardi G. Performance of the VITEK®2 advanced expert system™ for the validation of antimicrobial susceptibility testing results. Eur J Clin Microbiol Infect Dis. 2021;40:1333–5.

    Article  CAS  PubMed  Google Scholar 

  58. Kowalska-Krochmal B, Dudek-Wicher R. The minimum inhibitory concentration of antibiotics: methods, interpretation, clinical relevance. Pathogens. 2021;10:1–21.

    Article  Google Scholar 

  59. Weinstein MP, Lewis JS. The clinical and laboratory standards institute subcommittee on antimicrobial susceptibility testing: background, organization, functions, and processes. J Clin Microbiol. 2020;58.

  60. Magiorakos AP, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect. 2012;18:268–81.

    Article  CAS  PubMed  Google Scholar 

  61. Krumperman PH. Multiple antibiotic resistance indexing of Escherichia coli to identify high-risk sources of fecal contamination of foods. Appl Environ Microbiol. 1983;46:165–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Khalid S, Ahmad N, Ali SM, Khan AU. Outbreak of efficiently transferred Carbapenem-resistant blaNDM-producing gram-negative Bacilli isolated from neonatal intensive care unit of an Indian hospital. Microb Drug Resist. 2020;26:284–9.

    Article  CAS  PubMed  Google Scholar 

  63. Moawad AA, Hotzel H, Awad O, Tomaso H, Neubauer H, Hafez HM, et al. Occurrence of salmonella enterica and Escherichia coli in raw chicken and beef meat in northern Egypt and dissemination of their antibiotic resistance markers. Gut Pathog. 2017;9.

  64. Feng S, DeKlotz M, Taş N. Comparison of three DNA extraction methods for recovery of microbial DNA from Arctic permafrost. MicroPubl Biol. 2023;2023.

  65. Fisher Scientific T. Qubit ® 3.0 Fluorometer Catalog Number Q33216.

  66. Babraham Bioinformatics - FastQC A Quality Control tool for High Throughput Sequence Data.

  67. Ewels P, Ns Magnusson M, Lundin S, Aller MK. Data and text mining MultiQC: summarize analysis results for multiple tools and samples in a single report. https://doi.org/10.1093/bioinformatics/btw354.

  68. Behzadi P, Ranjbar R. DNA microarray technology and bioinformatic web services. Acta Microbiol Immunol Hung. 2019;66:19–30.

    Article  CAS  PubMed  Google Scholar 

  69. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13:e1005595.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Han N, Qiang Y, Zhang W. ANItools web: a web tool for fast genome comparison within multiple bacterial strains. Database (Oxford). 2016;2016.

  71. Clarridge JE. Impact of 16S rRNA gene sequence analysis for identification of Bacteria on clinical microbiology and infectious diseases. Clin Microbiol Rev. 2004;17:840.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;3.

  73. Galardini M, Biondi EG, Bazzicalupo M, Mengoni A. CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes. Source Code Biol Med. 2011;6:1–5.

    Article  Google Scholar 

  74. Seemann T. Genome analysis Prokka: rapid prokaryotic genome annotation. 2014;30:2068–9.

  75. Robertson J, Nash JHE. MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb Genom. 2018;4(8):e000206. https://doi.org/10.1099/mgen.0.000206.

  76. Carattoli A, Zankari E, Garciá-Fernández A, Larsen MV, Lund O, Villa L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58:3895–903.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10.

  78. Lam MMC, Wick RR, Watts SC, Cerdeira LT, Wyres KL, Holt KE. A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex. Nat Commun. 2021;12.

  79. Thrane SW, Taylor VL, Lund O, Lam JS, Jelsbak L. Application of whole-genome sequencing data for O-specific antigen analysis and in Silico serotyping of Pseudomonas aeruginosa isolates. J Clin Microbiol. 2016;54:1782–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31:3812.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Alcock BP, Raphenya AR, Lau TT, Tsang KK, Egane B, Edalatmand A, et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2019;48.

  82. Couvin D, Bernheim A, Toffano-Nioche C, Touchon M, Michalik J, Néron B, et al. CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res. 2018;46:W246–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Shankar C, Kumar S, Venkatesan M, Veeraraghavan B. Emergence of ST147 Klebsiella pneumoniae carrying blaNDM-7 on IncA/C2 with ompK35 and ompK36 mutations in India. J Infect Public Health. 2019;12:741–3.

    Article  PubMed  Google Scholar 

  84. Farhat N, Khan AU. Therapeutic approaches to combat the global antibiotic resistance challenge. Future Microbiol. 2022;17:1515–29.

    Article  CAS  PubMed  Google Scholar 

  85. Farhat N, Khan AU. Evolving trends of New Delhi Metallo-betalactamse (NDM) variants: a threat to antimicrobial resistance. Infect Genet Evol. 2020;86.

  86. Talat A, Blake KS, Dantas G, Khan AU. Metagenomic insight into microbiome and antibiotic resistance genes of high clinical concern in urban and rural hospital wastewater of northern India origin: a major reservoir of antimicrobial resistance. Microbiol Spectr. 2023;11.

  87. Gharaibeh MH, Shatnawi SQ. An overview of colistin resistance, mobilized colistin resistance genes dissemination, global responses, and the alternatives to colistin: a review. Veterinary World. 2019;12:1735–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Bardhan T, Chakraborty M, Bhattacharjee B. Prevalence of colistin-resistant, carbapenem-hydrolyzing proteobacteria in hospital water bodies and out-falls of West Bengal, India. Int J Environ Res Public Health. 2020;17.

  89. Guducuoglu H, Gursoy NC, Yakupogullari Y, Parlak M, Karasin G, Sunnetcioglu M, et al. Hospital outbreak of a Colistin-resistant, NDM-1- and OXA-48-producing Klebsiella pneumoniae: high mortality from Pandrug resistance. Microb Drug Resist. 2018;24:966–72.

    Article  CAS  PubMed  Google Scholar 

  90. Carroll LM, Gaballa A, Guldimann C, Sullivan G, Henderson LO, Wiedmann M. Identification of novel mobilized colistin resistance gene mcr-9 in a multidrug-resistant, colistin-susceptible salmonella enterica serotype typhimurium isolate. mBio. 2019;10.

  91. Fukuda A, Nakano H, Suzuki Y, Nakajima C, Usui M. Conjugative IncHI2/HI2A plasmids harbouring mcr-9 in colistin-susceptible Escherichia coli isolated from diseased pigs in Japan. Access Microbiol. 2022;4.

  92. Talat A, Usmani A, Khan AU. Detection of E. Coli IncX1 plasmid-mediated mcr-5.1 gene in an Indian hospital sewage water using shotgun metagenomic sequencing: a first report. Microb Drug Resist. 2022;28.

  93. Talat A, Miranda C, Poeta P, Khan AU. Farm to table: colistin resistance hitchhiking through food. Arch Microbiol. 2023;205.

  94. Laxminarayan R. The overlooked pandemic of antimicrobial resistance. Lancet. 2022;399:606–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Mihara T, Kimura T, Momiyama K, Kainuma A, Akiyama K, Ohara J, et al. Secondary in-hospital epidemiological investigation after an outbreak of Pseudomonas aeruginosa ST357. J Infect Chemother. 2020;26:257–65.

    Article  CAS  PubMed  Google Scholar 

  96. Naha S, Sands K, Mukherjee S, Roy C, Rameez MJ, Saha B, et al. KPC-2-producing Klebsiella pneumoniae ST147 in a neonatal unit: clonal isolates with differences in colistin susceptibility attributed to AcrAB-TolC pump. Int J Antimicrob Agents. 2020;55:105903.

    Article  CAS  PubMed  Google Scholar 

  97. Yusof NY, Norazzman NII, Hakim SNWA, Azlan MM, Anthony AA, Mustafa FH, et al. Prevalence of mutated Colistin-resistant Klebsiella pneumoniae: a systematic review and Meta-analysis. Trop Med Infect Dis. 2022;7.

  98. Torres DA, Seth-Smith HMB, Joosse N, Lang C, Dubuis O, Nüesch-Inderbinen M, et al. Colistin resistance in gram-negative bacteria analysed by five phenotypic assays and inference of the underlying genomic mechanisms. BMC Microbiol. 2021;21.

  99. Liu S, Ding Y, Xu Y, Li Z, Zeng Z, Liu J. An outbreak of extensively drug-resistant and hypervirulent Klebsiella pneumoniae in an intensive care unit of a teaching hospital in Southwest China. Front Cell Infect Microbiol. 2022;12.

  100. Urooj M, Ullah R, Ali S, Mohyuddin A, Mirza HM, Faryal R. Elucidation of molecular mechanism for colistin resistance among gram-negative isolates from tertiary care hospitals. J Infect Chemother. 2022;28:602–9.

    Article  CAS  PubMed  Google Scholar 

  101. Colistin M, Genes R. Global Burden of Colistin-Resistant Bacteria: Mobilized Colistin Resistance Genes Study. 2019;:1–18.

  102. Lu J, Chen J, Liu C, Zeng Y, Sun Q, Li J, et al. Identification of antibiotic resistance and virulence-encoding factors in Klebsiella pneumoniae by Raman spectroscopy and deep learning. Microb Biotechnol. 2022;15:1270–80.

    Article  CAS  PubMed  Google Scholar 

  103. Ghaffoori Kanaan MH, Al-Shadeedi SMJ, Al-Massody AJ, Ghasemian A. Drug resistance and virulence traits of Acinetobacter baumannii from Turkey and chicken raw meat. Comp Immunol Microbiol Infect Dis. 2020;70.

  104. Kubelová M, Koláčková I, Gelbíčová T, Florianová M, Kalová A, Karpíšková R. Virulence properties of mcr-1-positive Escherichia coli isolated from retail poultry meat. Microorganisms. 2021;9(2):308. https://doi.org/10.3390/microorganisms9020308.

  105. Zaidi S, Ali K, Khan AU. It’s all relative: analyzing microbiome compositions, its significance, pathogenesis and microbiota derived biofilms: challenges and opportunities for disease intervention. Arch Microbiol. 2023;205.

  106. Kainuma A, Momiyama K, Kimura T, Akiyama K, Inoue K, Naito Y, et al. An outbreak of fluoroquinolone-resistant Pseudomonas aeruginosa ST357 harboring the exoU gene. J Infect Chemother. 2018;24:615–22.

    Article  CAS  PubMed  Google Scholar 

  107. Akhtar F, Khan AU, Qazi B, Kulanthaivel S, Mishra P, Akhtar K, et al. A nano phototheranostic approach of toluidine blue conjugated gold silver core shells mediated photodynamic therapy to treat diabetic foot ulcer. Sci Rep. 2021;11.

  108. Uzairue LI, Rabaan AA, Adewumi FA, Okolie OJ, Folorunso JB, Bakhrebah MA, et al. Global prevalence of Colistin resistance in Klebsiella pneumoniae from bloodstream infection: a systematic review and Meta-analysis. Pathogens. 2022;11.

Download references

Acknowledgements

We acknowledge the facility and support provided by the Department of Biotechnology, Ministry of Science and Technology, Government of India. We acknowledge Nabeela Farhat and Shamsi Khalid for their assistance in the antibiotic sensitivity test and for procuring pure cultures of strains.

Funding

We acknowledge funding from the DBT, Government of India, grant no. BT/PR40148/BTIS/137/20/2021 and Tata Innovation Fellowship, HRD-16012/6/2020-AFS-DBT to AUK. A.T. is supported by the Council of Scientific and Industrial Research in the form of CSIR-NET/SRF (09/112(0648)/2019-EMR-I). The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies.

Author information

Authors and Affiliations

Authors

Contributions

AT performed sequence analysis, MIC analysis and wrote draft manuscript, FK collected clinical samples and did primary screening, and AUK conceived the idea, guided the study and checked draft manuscript.

Corresponding author

Correspondence to Asad U. Khan.

Ethics declarations

Ethics approval and consent to participate

The Institutional Ethics Committee approved the study (D. No. 269/FM/IEC) under National Ethics Committee Registry for Biomedical and Health Research-NECRBHR DHR-ICMR, Jawaharlal Nehru Medical College and Hospital, Faculty of Medicine, Aligarh Muslim University. The samples were collected after obtaining informed consent from patients.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Talat, A., Khan, F. & Khan, A.U. Genome analyses of colistin-resistant high-risk blaNDM-5 producing Klebsiella pneumoniae ST147 and Pseudomonas aeruginosa ST235 and ST357 in clinical settings. BMC Microbiol 24, 174 (2024). https://doi.org/10.1186/s12866-024-03306-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03306-4

Keywords