Methicillin resistant Staphylococcus aureus (MRSA) is currently a major cause of skin and soft tissue infections (SSTI) in the United States. Seasonal variation of MRSA infections in hospital settings has been widely observed. However, systematic time-series analysis of incidence data is desirable to understand the seasonality of community acquired (CA)-MRSA infections at the population level. In this paper, using data on monthly SSTI incidence in children aged 0–19 years and enrolled in Medicaid in Maricopa County, Arizona, from January 2005 to December 2008, we carried out time-series and nonlinear regression analysis to determine the periodicity, trend, and peak timing in SSTI incidence in children at different age: 0-4 years, 5-9 years, 10-14 years, and 15-19 years. We also assessed the temporal correlation between SSTI incidence and meteorological variables including average temperature and humidity. Our analysis revealed a strong annual seasonal pattern of SSTI incidence with peak occurring in early September. This pattern was consistent across age groups. Moreover, SSTIs followed a significantly increasing trend over the 4-year study period with annual incidence increasing from 3.36% to 5.55% in our pediatric population of approximately 290,000. We also found a significant correlation between the temporal variation in SSTI incidence and mean temperature and specific humidity. Our findings could have potential implications on prevention and control efforts against CA-MRSA.

Community associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has become a major cause of skin and soft tissue infections (SSTIs) in the US. We developed an age-structured compartmental model to study the spread of CA-MRSA at the population level and assess the effect of control intervention strategies. We used Monte-Carlo Markov Chain (MCMC) techniques to parameterize our model using monthly time series data on SSTIs incidence in children (≤19 years) during January 2004 -December 2006 in Maricopa County, Arizona. Our model-based forecast for the period January 2007–December 2008 also provided a good fit to data. We also carried out an uncertainty and sensitivity analysis on the control reproduction number, R_c which we estimated at 1.3 (95% CI [1.2,1.4]) based on the model fit to data. Using our calibrated model, we evaluated the effect of typical intervention strategies namely reducing the contact rate of infected individuals owing to awareness of infection and decolonization strategies targeting symptomatic infected individuals on both and the long-term disease dynamics. We also evaluated the impact of hypothetical decolonization strategies targeting asymptomatic colonized individuals. We found that strategies focused on infected individuals were not capable of achieving disease control when implemented alone or in combination. In contrast, our results suggest that decolonization strategies targeting the pediatric population colonized with CA-MRSA have the potential of achieving disease elimination.