Genomic Epidemiology and Characterization of Methicillin-Resistant Staphylococcus aureus from Bloodstream Infections in China

Since 2010, methicillin-resistant Staphylococcus aureus (MRSA) ST59 began to increase in prevalence in China, gradually replacing ST239 and has become the dominant clone in most hospitals in China. Here, we investigated the changing epidemiology, phylogenetic reconstruction, and genomic characterization of MRSA clones in China to identify the genomic driving factors in the prevalence of ST59. Most MRSA isolates were identified as ST59 (36.98%; 277/749), which increased from 25.09% in 2014 to 35.53% in 2019. The phylogenetic analysis of the 749 MRSA isolates showed a high level of diversity and the copresence of hospital-associated, community-associated, livestock-associated, and hypervirulent clones. Furthermore, minimum spanning trees revealed that ST59 MRSA clones from different hospitals and regions were integrated, suggesting that frequent exchanges had occurred between regions and hospitals. ST59 clones displayed higher susceptibility to antimicrobials than did ST239 and ST5 MRSA clones, indicating that resistance to non-β-lactam and fluoroquinolone antibiotics may be not critical for the epidemic success of ST59 clones. Virulence factors detection showed that sak and chp genes enriched in MRSA ST59 may be associated with the enhanced spreading success of ST59, whereas qacA may have contributed to the predominance of ST5 in East China. Our refined analysis of different clones among ST239, ST5, ST59, and ST398 demonstrated the existence of potential driving factors for the evolution of nosocomial MRSA populations and diversity of the evolutionary events surrounding clonal replacement. IMPORTANCE As a developing country, China has an unbalanced health care system due to regional differences in economic development. However, China is also a country worthy of study with regard to the population dynamics of MRSA within the more resource-rich health care systems. In this study, we carried out genomic analysis to investigate the genomic epidemiology and characterization of MRSA isolated from bloodstream infections over a timespan of 6 years. Our refined analysis of different MRSA clones among ST59, ST5, ST239, and ST398 demonstrated the existence of driving factors for the evolution of nosocomial MRSA populations and diversity of the evolutionary events surrounding clonal replacement.