Abstract
Pseudomonas aeruginosa (P. aeruginosa) is one of the most common nosocomial pathogens worldwide. Although the emergence of multidrug-resistant (MDR) P. aeruginosa is a critical problem in medical practice, the key features involved in the emergence and spread of MDR P. aeruginosa remain unknown. This study utilized whole genome sequence (WGS) analyses to define the population structure of 185 P. aeruginosa clinical isolates from several countries. Of these 185 isolates, 136 were categorized into sequence type (ST) 235, one of the most common types worldwide. Phylogenetic analysis showed that these isolates fell within seven subclades. Each subclade harbors characteristic drug resistance genes and a characteristic genetic background confined to a geographic location, suggesting that clonal expansion following antibiotic exposure is the driving force in generating the population structure of MDR P. aeruginosa. WGS analyses also showed that the substitution rate was markedly higher in ST235 MDR P. aeruginosa than in other strains. Notably, almost all ST235 isolates harbor the specific type IV secretion system and very few or none harbor the CRISPR/CAS system. These findings may help explain the mechanism underlying the emergence and spread of ST235 P. aeruginosa as the predominant MDR lineage.