Introduction
Sequence type 398 (ST398) Staphylococcus aureus emerged as a novel type globally in the early 2000s. Reference Bhat, Dumortier and Taylor1,Reference Bouiller, Bertrand, Hocquet and Chirouze2 ST398 with methicillin resistance (MRSA) was commonly associated with livestock contact, while ST398 with methicillin susceptibility (MSSA) was not. Reference Uhlemann, Porcella and Trivedi3 Genomic analysis indicates further contrasts between livestock-associated ST398 MRSA and human-associated ST398 MSSA, such as the varying presence of genes encoding resistance to tetracycline and erythromycin. Reference Bouiller, Bertrand, Hocquet and Chirouze2,Reference Uhlemann, McAdam and Sullivan4 Although ST398 MSSA was widely detected in European and Asian countries, in the United States it was first detected in New York City (NYC) and the subsequent region of greatest prevalence was New York state and the Northeast. Reference Uhlemann and Lowy5–Reference Uhlemann, Hafer and Miko8 Utilizing the phenotypic characteristics commonly associated with ST398 MSSA, prior analysis of MSSA collected nationwide in US hospitals documented an increased prevalence of potential ST398 from 2003 to 2014 with a significant negative association with distance from NYC. Reference Uhlemann, Hafer and Miko8,Reference Carrel, Goto, Schweizer, David, Livorsi and Perencevich9 Molecular characterization of MSSA isolates across the United States from 2004 to 2010 similarly found a low prevalence of this phenotype outside of the Northeast. Reference Miko, Hafer and Lee10
This phenotypic profile has been shown to be associated, but perhaps not exclusively, with ST398 MSSA. However, in the absence of widespread genotyping of S. aureus isolates, spatiotemporal surveillance of phenotypic profiles can indicate the prevalence and diffusion patterns of novel antimicrobial resistance varieties. Utilizing a retrospective cohort of all S. aureus positive cultures in outpatient Veterans Health Administration (VHA) settings from 2010 to 2019, we explored whether previously observed rates of potential-ST398 MSSA held steady or changed over time and whether the Northeast continued to serve as the epicenter for this phenotype.
Data and methods
All S. aureus isolates from clinical specimens collected in outpatient VHA settings from January 1, 2010 to December 31, 2019 were analyzed. Microbiology reports were extracted from the Corporate Data Warehouse, the integrated electronic medical record system data repository of the VHA health system. We included all isolates from patients with a residential address in the contiguous United States who were ≥18 yr old at the time of specimen collection. Isolates were excluded from analysis if they were obtained at institutionalized care settings (e.g., acute care, nursing homes, rehabilitation facilities) or if they were taken after either 48 h of an inpatient admission or 72 h of an inpatient discharge. MSSA was defined as susceptibility to first-generation cephalosporins, cefoxitin, nafcillin, oxacillin, or carbapenems. Only isolates that were tested for resistance against all three relevant antimicrobial classes, clindamycin, erythromycin, and tetracycline, were included in the analysis. The first record per unique patient per month was retained. Potential ST398 was defined as resistance, including intermediate resistance, to clindamycin and erythromycin with susceptibility to tetracycline. Reference Gandra, Braykov and Laxminarayan7,Reference Uhlemann, Hafer and Miko8 When reported, inducible clindamycin resistance (commonly known as D-test) was also considered to be resistance.
Based on patient county of residence, MSSA isolates were assigned to Census regions: Northeast, Midwest, South, or West. Prevalence of potential ST398 MSSA was determined at county and regional scales for each month and year of the study.
With the county-level monthly counts of potential ST398 as the outcome variable, a Poisson regression model with robust standard errors was fitted to account for the variation in the trends of potential ST398 between 2010 and 2019 as repeated measurements within each region. Region was included in the model as a fixed effect and an interaction term between time and region was included in the model to detect regional variability in temporal trends. The total number of MSSA infections was included as an offset variable. Baseline incidence rates and temporal trends were estimated as incidence rate ratios (IRRs), using Northeast and the year 2010 as references.
The study was approved by the IRB of the University of Iowa.
Results
In 2010–2019, 235,033 outpatient specimens were positive for MSSA and 161,816 (68.85%) were tested against all three relevant antimicrobial classes. Among the first 145,304 isolates per unique patient per month, 22,919 (15.77%) were classified as potential ST398. Rates of potential ST398 rose slightly from ∼15% in the early years of the study to ∼17% in the final years (Figure 1A). Prevalence differed by Census region, with the Northeast having the highest prevalence of potential ST398 and the West having the greatest increases in rates over time (Figure 1B).
Models accounting for region and time showed similar patterns of potential ST398 MSSA (Table 1). No significant change in prevalence over time was observed (IRR: 1.01 [0.99–1.02]). The Midwest, South, and West all had significantly lower potential ST398 MSSA than did the Northeast. When the interaction between time and region was examined, only the West had significantly increased potential ST398 from 2010 to 2019 (IRR: 1.02 [1.01–1.05]) compared to the Northeast.
Note. Bold is statistically significant (P < 0.05).
Discussion
Prior analysis of sterile site MSSA cultures within VHA inpatient settings indicated an increase in potential ST398 MSSA from 2003 to 2010, with prevalence plateauing from 2010 to 2014 around 15%. Reference Carrel, Goto, Schweizer, David, Livorsi and Perencevich9 In this study which includes data through 2019, similar prevalence is observed in cultures from all sites obtained in VHA outpatient settings. This suggests that potential ST398 MSSA has established stable levels of circulation in the US population.
While prevalence over time is relatively static, regional differences persist. Potential ST398 MSSA prevalence remained highest in outpatients residing in the Northeast, consistent with earlier studies that documented early detection of ST398 MSSA in NYC. Significant increases in potential ST398 MSSA are also observed in outpatients residing in the Western region of the United States. Long-distance spread and intense short-distance circulation of ST398 MSSA have been documented using phylogeographic analysis, pointing to migration routes and human circulation as responsible for the diffusion and establishment of this novel type. Reference Uhlemann, McAdam and Sullivan4
Understanding these geographic differences in potential ST398 prevalence, and its associated phenotypic profile, has important implications for empiric therapy, particularly in outpatient settings. In the Northeast, over 20% of the MSSA infections observed in outpatient settings were resistant to clindamycin and erythromycin, and the rising rates of resistance may influence clinicians in that region to change their empiric prescribing patterns based on their risk tolerance. At the same time, clinicians in the West need to be aware of the novel resistance profile that is slowly increasing in prevalence among their patients. While beta-lactam agents (e.g., first-generation cephalosporins) are generally preferred options for MSSA, most empiric therapy decisions, particularly in outpatient settings, are made without susceptibility profile information. As such, epidemiologic trends of phenotypic change can be useful information to guide clinicians’ prescribing patterns.
Sample collection was limited to the end of 2019 to avoid COVID pandemic confounding factors such as decreased outpatient visits and increased antimicrobial usage. Future work will explore potential ST398 MSSA prevalence to understand if spatiotemporal patterns have changed in VHA populations following the pandemic.
Though its patient population is older with higher male proportions than the rest of the US population, the VHA represents one of the only sources for nationwide retrospective phenotypic analysis of bacterial isolates in the United States. While the phenotypic profile used in this study may incorrectly identify ST398 MSSA, it is still useful to observe spatiotemporal trends in resistance profiles. Moreover, many studies focus on MRSA or invasive infections rather than more frequently experienced outpatient illnesses. Understanding the resistance patterns and associated genotypes occurring in MSSA and the community rather than in healthcare settings is an important indicator of the successful spread and persistence of novel S. aureus types in the non-hospitalized population.
Financial support
M.C., Q.S., and G.S.C. received support from the Centers for Disease Control and Prevention Epicenter program award CK000613. M.G. received support from the Agency for Healthcare Research and Quality under award number K08HS027472. E.N.P. received support from the Veterans Affairs Health Services Research and Development (HSR&D) Service Center of Innovation under award number 13-412. M.S., S.H., and M.Z.D. report no relevant funding sources. The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs or the United States Government.
Competing interests
M.Z.D. has received funding for a research contract from Johnson and Johnson and served as local PI for clinical trials supported by GSK and Covance. He also served on an advisory board for GSK. All other authors report no conflicts of interest relevant to this article.