Genomic Surveillance and Resistome in ICUs: Unraveling Klebsiella pneumoniae and Acinetobacter baumannii with an Interactive Dashboard

Bacterial Genomics Clinical Bioinformatics R Shiny Antimicrobial Resistance (AMR)

Data repository: NCBI BioProject PRJNA1185159 | SRA Study SRP545647

1. Introduction: The Challenge of Antimicrobial Resistance in Brazilian ICUs

Antimicrobial resistance (AMR) in opportunistic bacterial pathogens is one of the most severe threats to global public health. In the Intensive Care Unit (ICU), where immunocompromised patients with invasive devices face serious hospital-acquired infections, multidrug-resistant Gram-negative bacteria such as carbapenemase-producing (KPC) Klebsiella pneumoniae and carbapenem-resistant Acinetobacter baumannii (CRAB) lead to challenging clinical outcomes and a shortage of therapeutic options.

However, monitoring the spread of these pathogens with traditional phenotypic tests alone is no longer enough. We need to investigate the genomic architecture of circulating strains: which clones are expanding? Which resistance determinants sit on plasmids or integrons? How do chromosomal mutations affect last-line antimicrobials such as colistin?

In this article, we share the behind-the-scenes work and genomic findings of a large prospective multicenter study conducted in six tertiary hospitals in São Paulo, coordinated by Prof. Dr. Eduardo Alexandrino Servolo de Medeiros (Infectious Diseases Division / LEMC / UNIFESP), in which I worked directly on genomic curation, bioinformatic auditing and the development of analytical tools. To make this large dataset accessible, we built a complete analytical dashboard, available interactively right here on the blog!

Felipe Alberto Lei

Founder of LogikBioinfo. Bioinformatician specialized in bacterial genomics, antimicrobial resistance and the development of analytical tools.

2. The Multicenter Cohort and the Bioinformatics Pipeline

The study followed ICUs at six major hospitals in São Paulo State: Hospital São Paulo (HSP), Hospital Estadual de Diadema (HED), Hospital Municipal de Barueri (HMB), Hospital Geral de Pirajussara (HGPI), Hospital Geral de Guarulhos (HGG) and Hospital Geral de Pedreira (HGPE).

Starting from 290 initial raw samples from active surveillance (rectal swabs) and diagnostic cultures, we applied strict quality control (QC) and epidemiological inclusion criteria, consolidating a definitive cohort of 236 isolates with Whole Genome Sequencing (WGS) on the Illumina NextSeq 1000 platform:

236

WGS isolates

144

K. pneumoniae

92

A. baumannii

The Decisive Role of Bioinformatic Curation: Practical Lessons

In routine clinical bioinformatics, automated "out-of-the-box" pipelines often produce pitfalls that only expert curation can uncover:

  1. The mgrB trap in colistin resistance: in K. pneumoniae, mgrB inactivation is the main cause of polymyxin resistance. However, standard point-SNP alignment tools reported that most resistant isolates had no mgrB mutation. By inspecting base coverage, we found that 28 resistant isolates had 0% coverage of the mgrB gene — that is, they had undergone complete structural deletion or disruption by insertion sequences (IS). Once the artifact was corrected, 77.8% (35/45) of colistin-resistant isolates turned out to carry functional mgrB inactivation!
  2. Detection and exclusion of chimeras/contamination: isolate HED-117-ACB stood out with 304 virulence genes (versus an average of 111 for Acinetobacter). An in-depth audit revealed complete biosynthesis pathways for pyocyanin, pyoverdine and type III secretion systems exclusive to Pseudomonas aeruginosa. It was pre-sequencing microbiological contamination that, had it not been audited, would have seriously distorted the group's publications.
  3. International breakpoint harmonization (BrCAST/EUCAST vs. CLSI): we aligned reports between the Brazilian ANVISA standard (BrCAST v14.0) and the international standard (CLSI M100). For colistin in A. baumannii, for example, the "Susceptible" category does not exist in CLSI (MIC ≤ 2 µg/mL is classified as Intermediate), requiring dual reporting for compatibility with high-impact international journals.

3. Key Biological Findings Revealed by WGS

Klebsiella pneumoniae: the ST11 vs. ST258 clonal dichotomy

  • ST11 (41.0%): the most prevalent clone in the study, widely spread across HGPI, HED and HMB. It is characterized by a dangerous combination of multidrug resistance and hypervirulence: 61% of ST11 isolates are colistin-resistant and carry the mobile element ICEKp10, encoding yersiniabactin (ybt) and colibactin (clb).
  • ST258 (18.1%): a classic lineage associated with the blaKPC-2* variant (with the conserved G822T mutation). Notably, 100% of ST258 isolates remained fully susceptible to colistin (MIC = 1 µg/mL).
  • Carbapenemases: near-absolute predominance of blaKPC-2 / blaKPC-2* (95.8%), with blaNDM-1 detected in 2.1% of isolates.

Acinetobacter baumannii: the NDM-1 + armA cluster alert in Diadema

  • Dominant oxacillinases: blaOXA-23 was identified in 78.3% of isolates and blaOXA-72 in 20.7%.
  • Emerging cluster at HED: we detected a high-risk clonal cluster at Hospital Estadual de Diadema made up of 6 isolates of the international clone ST2 (Global Clone 2 - GC2) co-producing blaNDM-1 + blaOXA-23/66.
  • 100% co-occurrence with armA: all NDM-1 isolates carry the 16S rRNA methyltransferase armA, conferring pan-resistance to aminoglycosides (gentamicin, amikacin and tobramycin with MIC > 64 µg/mL).

4. The Interactive Dashboard

To turn these analyses into a visual and dynamic experience, we developed an analytical dashboard. The full version was built in R Shiny with the modern bslib framework (Flatly theme), interactive plotly charts and reactive DT tables. Below is a lightweight web version of the same panel, with data for all 236 isolates, running directly in your browser.

Dashboard features

  1. Reactive dynamic filters: instantly filter by species (K. pneumoniae, A. baumannii) or by specific hospitals.
  2. Epidemiological indicators: total filtered samples, KPC, oxacillinase and NDM-1 cases, and colistin-resistant isolates.
  3. Species modules: ST distribution, colistin by clone and detailed mgrB status.
  4. Searchable data table: individual access to each of the 236 isolates, with global search and sorting by any column.

Try the interactive dashboard

Data: Consolidated General Manifest (N = 236), UNIFESP/LEMC Multicenter Project. No identifiable patient data is displayed. Colistin phenotype in A. baumannii: S (BrCAST) / I (CLSI) for MIC ≤ 2 µg/mL.

5. How to Run and Publish the Shiny Application

If you want to run the R Shiny version of the dashboard on your own machine or deploy it to shinyapps.io, the code is fully prepared and modularized.

Running locally in R / RStudio

# Install the required packages
install.packages(c("shiny", "bslib", "DT", "plotly", "rsconnect", "bsicons"))

# Run the application
shiny::runApp("app.R")

Publishing to shinyapps.io

Just insert your shinyapps.io credentials (Account → Tokens → Show):

library(rsconnect)

rsconnect::setAccountInfo(
  name   = "logikbioinfo",
  token  = "YOUR_TOKEN",
  secret = "YOUR_SECRET"
)

rsconnect::deployApp(
  appDir        = getwd(),
  appPrimaryDoc = "app.R",
  appName       = "amr-genomics-multicenter-unifesp"
)

6. Conclusion and LogikBioinfo Services

Integrating classical microbiology, next-generation sequencing and curated bioinformatics is today's gold standard for high-impact research and health surveillance. Cases such as elucidating the mgrB mechanism and identifying the NDM-1 cluster in Diadema show that detailed bioinformatic analysis goes far beyond generating raw tables: it turns nucleotide sequences into actionable epidemiological knowledge.

Need bioinformatics support for your research project?

At LogikBioinfo, we offer complete scientific consulting for researchers, research groups and clinical laboratories:

  • Genome assembly and resistome/virulome analysis (WGS)
  • Phylogenomic and pangenomic analyses
  • Custom dashboards in R Shiny and Python (Streamlit/Dash)
  • Data submission and curation in public databases (NCBI BioProject, BioSample, SRA, GenBank)
  • Hands-on bioinformatics training with real data from your project
Previous post Back to Blog