Despite its importance as a human pathogen, information on population structure and global epidemiology of Staphylococcus epidermidis is scarce and the relative importance of the mechanisms contributing to clonal diversification is unknown. In this study, we addressed these issues by analyzing a representative collection of S. epidermidis isolates from diverse geographic and clinical origins using multilocus sequence typing (MLST). Additionally, we characterized the mobile element (SCC mec ) carrying the genetic determinant of methicillin resistance. The 217 S. epidermidis isolates from our collection were split by MLST into 74 types, suggesting a high level of genetic diversity. Analysis of MLST data using the eBURST algorithm revealed the existence of nine epidemic clonal lineages that were disseminated worldwide. One single clonal lineage (clonal complex 2) comprised 74% of the isolates, whereas the remaining isolates were clustered into 8 minor clonal lineages and 13 singletons. According to our evolutionary model, SCC mec was acquired at least 56 times by S. epidermidis . Although geographic dissemination of S. epidermidis strains and the value of the index of association between the alleles, 0.2898 ( P < 0.05), support the clonality of S. epidermidis species, examination of the sequence changes at MLST loci during clonal diversification showed that recombination gives rise to new alleles approximately twice as frequently as point mutations. We suggest that S. epidermidis has a population with an epidemic structure, in which nine clones have emerged upon a recombining background and evolved quickly through frequent transfer of genetic mobile elements, including SCC mec .

We evaluated three multilocus sequence typing (MLST) schemes for Staphylococcus epidermidis and selected the seven most discriminatory loci for the formation of a new, more powerful MLST scheme. This improved scheme gave 31 sequence types (STs) and 5 clonal complexes (CCs), whereas the other schemes delineate 16 to 24 STs and 1 to 3 CCs.

Some of the most common infectious diseases are caused by bacteria that naturally colonise humans asymptomatically. Combating these opportunistic pathogens requires an understanding of the traits that differentiate infecting strains from harmless relatives. Staphylococcus epidermidis is carried asymptomatically on the skin and mucous membranes of virtually all humans but is a major cause of nosocomial infection associated with invasive procedures. Here we address the underlying evolutionary mechanisms of opportunistic pathogenicity by combining pangenome-wide association studies and laboratory microbiology to compare S. epidermidis from bloodstream and wound infections and asymptomatic carriage. We identify 61 genes containing infection-associated genetic elements (k-mers) that correlate with in vitro variation in known pathogenicity traits (biofilm formation, cell toxicity, interleukin-8 production, methicillin resistance). Horizontal gene transfer spreads these elements, allowing divergent clones to cause infection. Finally, Random Forest model prediction of disease status (carriage vs. infection) identifies pathogenicity elements in 415 S. epidermidis isolates with 80% accuracy, demonstrating the potential for identifying risk genotypes pre-operatively.

The understanding of the mechanisms of antibiotic resistance development are fundamental to alert and preview beforehand, the large scale dissemination of resistance to antibiotics, enabling the design of strategies to prevent its spread. The mecA-mediated methicillin resistance conferring resistance to broad-spectrum β-lactams is globally spread in staphylococci including hospitals, farms and community environments, turning ineffective the most widely used and efficient class of antibiotics to treat staphylococcal infections. The use of whole genome sequencing (WGS) technologies at a bacterial population level has provided a considerable progress in the identification of key steps that led to mecA-mediated β-lactam resistance development and dissemination. Data obtained from multiple studies indicated that mecA developed from a harmless core gene (mecA1) encoding the penicillin-binding protein D (PbpD) from staphylococcal species of animal origin (S. sciuri group) due to extensive β-lactams use in human created environments. Emergence of the resistance determinant involved distortion of PbpD active site, increase in mecA1 expression, addition of regulators (mecR1, mecI) and integration into a mobile genetic element (SCCmec). SCCmec was then transferred into species of coagulase-negative staphylococci (CoNS) that are able to colonize both animals and humans and subsequently transferred to S. aureus of human origin. Adaptation of S. aureus to the exogenously acquired SCCmec involved, deletion and mutation of genes implicated in general metabolism (auxiliary genes) and general stress response and the adjustment of metabolic networks, what was accompanied by an increase in β-lactams minimal inhibitory concentration and the transition from a heterogeneous to homogeneous resistance profile. Nowadays, methicillin-resistant S. aureus (MRSA) carrying SCCmec constitutes one of the most important worldwide pandemics. The stages of development of mecA-mediated β-lactam resistance described here may serve as a model for previewing and preventing the emergence of resistance to other classes of antibiotics.

The opportunistic pathogens Staphylococcus aureus and Staphylococcus epidermidis represent major causes of severe nosocomial infection, and are associated with high levels of mortality and morbidity worldwide. These species are both common commensals on the human skin and in the nasal pharynx, but are genetically distinct, differing at 24% average nucleotide divergence in 1,478 core genes. To better understand the genome dynamics of these ecologically similar staphylococcal species, we carried out a comparative analysis of 324 S. aureus and S. epidermidis genomes, including 83 novel S. epidermidis sequences. A reference pan-genome approach and whole genome multilocus-sequence typing revealed that around half of the genome was shared between the species. Based on a BratNextGen analysis, homologous recombination was found to have impacted on 40% of the core genes in S. epidermidis, but on only 24% of the core genes in S. aureus. Homologous recombination between the species is rare, with a maximum of nine gene alleles shared between any two S. epidermidis and S. aureus isolates. In contrast, there was considerable interspecies admixture of mobile elements, in particular genes associated with the SaPIn1 pathogenicity island, metal detoxification, and the methicillin-resistance island SCCmec. Our data and analysis provide a context for considering the nature of recombinational boundaries between S. aureus and S. epidermidis and, the selective forces that influence realized recombination between these species.

Several studies have addressed the epidemiology of community-associated Staphylococcus aureus (CA-SA) in Europe; nonetheless, a comprehensive perspective remains unclear. In this study, we aimed to describe the population structure of CA-SA and to shed light on the origin of methicillin-resistant S. aureus (MRSA) in this continent. A total of 568 colonization and infection isolates, comprising both MRSA and methicillin-susceptible S. aureus (MSSA), were recovered in 16 European countries, from community and community-onset infections. The genetic background of isolates was characterized by molecular typing techniques (spa typing, pulsed-field gel electrophoresis and multilocus sequence typing) and the presence of PVL and ACME was tested by PCR. MRSA were further characterized by SCCmec typing. We found that 59% of all isolates were associated with community-associated clones. Most MRSA were related with USA300 (ST8-IVa and variants) (40%), followed by the European clone (ST80-IVc and derivatives) (28%) and the Taiwan clone (ST59-IVa and related clonal types) (15%). A total of 83% of MRSA carried Panton-Valentine leukocidin (PVL) and 14% carried the arginine catabolic mobile element (ACME). Surprisingly, we found a high genetic diversity among MRSA clonal types (ST-SCCmec), Simpson’s index of diversity = 0.852 (0.788–0.916). Specifically, about half of the isolates carried novel associations between genetic background and SCCmec. Analysis by BURP showed that some CA-MSSA and CA-MRSA isolates were highly related, suggesting a probable local acquisition/loss of SCCmec. Our results imply that CA-MRSA origin, epidemiology and population structure in Europe is very dissimilar from that of USA.