Comprehensive Genomic Investigation of Adaptive Mutations Driving the Low-Level Oxacillin Resistance Phenotype in Staphylococcus aureus

Antistaphylococcal penicillins such as oxacillin are the key antibiotics in the treatment of invasive methicillin-susceptible Staphylococcus aureus (MSSA) infections; however, mec gene-independent resistance adaptation can cause treatment failure. Despite its clinical relevance, the basis of this phenomenon remains poorly understood. Here, we investigated the genomic adaptation to oxacillin at an unprecedented scale using a large collection of 503 clinical mec-negative isolates and 30 in vitro-adapted isolates from independent oxacillin exposures. By combining comparative genomics, evolutionary convergence, and genome-wide association analysis, we found 21 genetic loci associated with low-level oxacillin resistance, underscoring the polygenic nature of this phenotype. Evidence of adaptation was particularly strong for the c-di-AMP signal transduction pathways (gdpP and dacA) and in the clpXP chaperone-protease complex. The role of mutations in gdpP in conferring low-level oxacillin resistance was confirmed by allele-swapping experiments. We found that resistance to oxacillin emerges at high frequency in vitro (median, 2.9 × 10(−6); interquartile range [IQR], 1.9 × 10(−6) to 3.9 × 10(−6)), which is consistent with a recurrent minimum inhibitory concentration (MIC) increase across the global phylogeny of clinical isolates. Nevertheless, adaptation in clinical isolates appears sporadically, with no stably adapted lineages, suggesting a high fitness cost of resistance, confirmed by growth assessment of mutants in rich media. Our data provide a broader understanding of the emergence and dynamics of oxacillin resistance adaptation in S. aureus and a framework for future surveillance of this clinically important phenomenon.