Genomic Characterization of Colistin-Resistant Isolates from the King Fahad Medical City, Kingdom of Saudi Arabia

Background: Whole-genome sequencing is one of the best ways to investigate resistance mechanisms of clinical isolates as well as to detect and identify circulating multi-drug-resistant (MDR) clones or sub-clones in a given hospital setting. Methods: Here, we sequenced 37 isolates of Acinetobacter baumannii, 10 Klebsiella pneumoniae, and 5 Pseudomonas aeruginosa collected from the biobank of the hospital setting of the King Fahad Medical City. Complete phenotypic analyses were performed, including MALDI-TOF identification and antibiotic susceptibility testing. After the genome assembly of raw data, exhaustive genomic analysis was conducted including full resistome determination, genomic SNP (gSNP) analysis, and comparative genomics. Results: Almost all isolates were highly resistant to all tested antibiotics, including carbapenems and colistin. Resistome analysis revealed many antibiotic resistance genes, including those with resistance to β-lactams, aminoglycosides, macrolides, tetracyclines, sulfamids, quinolones, and phenicols. In A. baumannii isolates, the endemic carbapenemase bla(OXA-23) gene was detected in 36 of the 37 isolates. Non-synonymous mutations in pmrB were detected in almost all of the isolates and likely mediated colistin resistance. Interestingly, while classical analyses, such as MLST, revealed the predominance of an ST2 clone in A. baumannii isolates, the genomic analysis revealed the presence of five circulating sub-clones and identified several isolate transmissions between patients. In the 10 K. pneumoniae isolates, several resistance genes were identified, and the observed carbapenem resistance was likely mediated by overexpression of the detected extended-spectrum-β-lactamase (ESBL) genes associated with low membrane permeability as few carbapenemase genes were detected with just bla(OXA-48) in three isolates. Colistin resistance was mediated either by non-synonymous mutations in the MgrB regulator, PmrA, PmrB, and PhoQ proteins or the presence of the MCR-1 protein. Here, gSNP analysis also revealed the existence of bacterial clones and cases of isolate transmissions between patients. The five analyzed P. aeruginosa isolates were highly resistant to all tested antibiotics, including carbapenems mediated by loss or truncated OprD porin, and colistin resistance was associated with mutations in the genes encoding the PmrA, PmrB, or PhoQ proteins. Conclusion: We demonstrate here the usefulness of whole-genome sequencing to exhaustively investigate the dissemination of MDR isolates at the sub-clone level. Thus, we suggest implementing such an approach to monitor the emergence and spread of new clones or sub-clones, which classical molecular analyses cannot detect. Moreover, we recommend increasing the surveillance of the endemic and problematic colistin resistance mcr-1 gene to avoid extensive dissemination.