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.
IMPORTANCE The majority of Staphylococcus aureus strains causing human disease are methicillin-susceptible (MSSA) and can be treated with antistaphylococcal penicillins (such as oxacillin). While acquisition of the mec gene represents the main resistance mechanism to oxacillin, S. aureus can acquire low-level resistance through adaptive mutations in other genes. In this study, we used genomic approaches to understand the basis of S. aureus adaption to oxacillin and its dynamic at the population level. By combining a genome analysis of clinical isolates from persistent MSSA infections, in vitro selection of oxacillin resistance, and genome-wide association analysis on a large collection of isolates, we identified 21 genes linked to secondary oxacillin resistance. Adaptive mutations in these genes were easy to select when S. aureus was exposed to oxacillin, but they also came at a substantial cost in terms of bacterial fitness, suggesting that this phenotype emerges preferentially in the setting of sustained antibiotic exposure.