Complex pathways to ceftolozane-tazobactam resistance in clinical Pseudomonas aeruginosa isolates: a genomic...

We aimed to conduct a comprehensive genomic analysis of ceftolozane/tazobactam (C/T) resistance mechanisms in Pseudomonas aeruginosa by combining novel institutional data with publicly-available sequencing data.

We analysed 1682 P. aeruginosa isolates, comprising 339 isolates from Alfred Hospital (Melbourne, Australia) and 1343 isolates from six public datasets. All isolates underwent whole-genome sequencing and C/T broth microdilution susceptibility testing. We assessed previously reported intrinsic and acquired resistance mechanisms. We then conducted a genome-wide association study and machine learning analysis to identify novel genes associated with resistance. We then evaluated the impact of mutations in these genes on MIC values and ceftolozane binding affinity.

Among 1682 P. aeruginosa isolates representing 527 distinct sequence types, 343 of 1682 (20.4%) were C/T-resistant. Carbapenemase genes were detected in 206 of 1682 (12.2%) isolates. Mutations in previously reported resistance-associated genes (ftsI, mpl, ampD, ampC, ampR, and oprD) were more frequent in resistant isolates but were also found in almost all susceptible isolates. Successive mutations conferred additive increases in MIC. A combined genome-wide association study and machine learning analyses a priori identified five key genes significantly associated with resistance: ftsI, ampR, ampC, PA3329, and PA4311. Molecular docking simulation revealed that the R504C mutation in penicillin-binding protein 3, which is encoded by ftsI, reduced binding contacts and hydrogen bonds with ceftolozane, significantly decreasing binding affinity (p 0.016).

Our analysis of 1682 P. aeruginosa genomes demonstrated complex pathways to C/T resistance and showed that ftsI may play an underappreciated role. We discovered two previously unidentified genes associated with C/T resistance, whose function remains to be determined.