Abstract
Ciprofloxacin is an important antibacterial drug targeting Type II topoisomerases, highly active against Gram-negatives including Escherichia coli. The evolution of resistance to ciprofloxacin in E. coli always requires multiple genetic changes, usually including mutations affecting two different drug target genes, gyrA and parC. Resistant mutants selected in vitro or in vivo can have many different mutations in target genes and efflux regulator genes that contribute to resistance. Among resistant clinical isolates the genotype, gyrA S83L D87N, parC S80I is significantly overrepresented suggesting that it has a selective advantage. However, the evolutionary or functional significance of this high frequency resistance genotype is not fully understood. By combining experimental data and mathematical modeling, we addressed the reasons for the predominance of this specific genotype. The experimental data were used to model trajectories of mutational resistance evolution under different conditions of drug exposure and population bottlenecks. We identified the order in which specific mutations are selected in the clinical genotype, showed that the high frequency genotype could be selected over a range of drug selective pressures, and was strongly influenced by the relative fitness of alternative mutations and factors affecting mutation supply. Our data map for the first time the fitness landscape that constrains the evolutionary trajectories taken during the development of clinical resistance to ciprofloxacin and explain the predominance of the most frequently selected genotype. This study provides strong support for the use of in vitro competition assays as a tool to trace evolutionary trajectories, not only in the antibiotic resistance field.