Abstract
Pseudomonas aeruginosa is characterized by a high adaptive potential, developing resistance in response to antimicrobial pressure. We employed a spatiotemporal evolution model to disclose the pathways of adaptation to colistin, a last-resort polymyxin antimicrobial, among three unrelated P. aeruginosa lineages. The P. aeruginosa ATCC-27833 reference strain (Pa_ATCC), an environmental P. aeruginosa isolate (Pa_Environment), and a clinical isolate with multiple drug resistance (Pa_MDR) were grown over an increasing 5-step colistin concentration gradient from 0 to 400 mg/L. Pa_Environment demonstrated the highest growth pace, achieving the 400 mg/L band in 15 days, whereas it took 37 and 60 days for Pa_MDR and Pa_ATCC, respectively. To identify the genome changes that occurred during adaptation to colistin, the isolates selected during the growth of the bacteria (n = 185) were subjected to whole genome sequencing. In total, 17 mutation variants in eight lipopolysaccharide-synthesis-associated genes were detected. phoQ and lpxL/PA0011 were affected in all three lineages, whereas changes in pmrB were found in Pa_Environment and Pa_MDR but not in Pa_ATCC. In addition, mutations were detected in 34 general metabolism genes, and each lineage developed mutations in a unique set of such genes. Thus, the three examined distinct P. aeruginosa strains demonstrated different capabilities and genetic pathways of colistin adaptation.