- Open Access
Whole-genome epidemiology, characterisation, and phylogenetic reconstruction of Staphylococcus aureus strains in a paediatric hospital
Genome Medicine volume 10, Article number: 82 (2018)
Staphylococcus aureus is an opportunistic pathogen and a leading cause of nosocomial infections. It can acquire resistance to all the antibiotics that entered the clinics to date, and the World Health Organization defined it as a high-priority pathogen for research and development of new antibiotics. A deeper understanding of the genetic variability of S. aureus in clinical settings would lead to a better comprehension of its pathogenic potential and improved strategies to contrast its virulence and resistance. However, the number of comprehensive studies addressing clinical cohorts of S. aureus infections by simultaneously looking at the epidemiology, phylogenetic reconstruction, genomic characterisation, and transmission pathways of infective clones is currently low, thus limiting global surveillance and epidemiological monitoring.
We applied whole-genome shotgun sequencing (WGS) to 184 S. aureus isolates from 135 patients treated in different operative units of an Italian paediatric hospital over a timespan of 3 years, including both methicillin-resistant S. aureus (MRSA) and methicillin-sensitive S. aureus (MSSA) from different infection types. We typed known and unknown clones from their genomes by multilocus sequence typing (MLST), Staphylococcal Cassette Chromosome mec (SCCmec), Staphylococcal protein A gene (spa), and Panton-Valentine Leukocidin (PVL), and we inferred their whole-genome phylogeny. We explored the prevalence of virulence and antibiotic resistance genes in our cohort, and the conservation of genes encoding vaccine candidates. We also performed a timed phylogenetic investigation for a potential outbreak of a newly emerging nosocomial clone.
The phylogeny of the 135 single-patient S. aureus isolates showed a high level of diversity, including 80 different lineages, and co-presence of local, global, livestock-associated, and hypervirulent clones. Five of these clones do not have representative genomes in public databases. Variability in the epidemiology is mirrored by variability in the SCCmec cassettes, with some novel variants of the type IV cassette carrying extra antibiotic resistances. Virulence and resistance genes were unevenly distributed across different clones and infection types, with highly resistant and lowly virulent clones showing strong association with chronic diseases, and highly virulent strains only reported in acute infections. Antigens included in vaccine formulations undergoing clinical trials were conserved at different levels in our cohort, with only a few highly prevalent genes fully conserved, potentially explaining the difficulty of developing a vaccine against S. aureus. We also found a recently diverged ST1-SCCmecIV-t127 PVL− clone suspected to be hospital-specific, but time-resolved integrative phylogenetic analysis refuted this hypothesis and suggested that this quickly emerging lineage was acquired independently by patients.
Whole genome sequencing allowed us to study the epidemiology and genomic repertoire of S. aureus in a clinical setting and provided evidence of its often underestimated complexity. Some virulence factors and clones are specific of disease types, but the variability and dispensability of many antigens considered for vaccine development together with the quickly changing epidemiology of S. aureus makes it very challenging to develop full-coverage therapies and vaccines. Expanding WGS-based surveillance of S. aureus to many more hospitals would allow the identification of specific strains representing the main burden of infection and therefore reassessing the efforts for the discovery of new treatments and clinical practices.
Staphylococcus aureus is a bacterium commonly found on the skin (15%), in the nostrils (27%), and in the pharynx (10–20%) of healthy adults [1,2,3], but it is also the cause of a number of diseases, whose severity ranges from common community-associated skin infections to fatal bacteraemia [3,4,5]. S. aureus is a leading cause of surgical, device-related, and pleuropulmonary infections, which can result into life-threatening infective endocarditis or even sepsis . The mortality of S. aureus invasive infections was extremely high (> 80%) in the pre-antibiotic era [6, 7], and only the introduction of penicillin at the beginning of the 1940s was able to contain it. However, resistant strains carrying a penicillinase/beta-lactamase quickly emerged [8,9,10], and more than 90% of current human-associated isolates are resistant to penicillin . Similarly, the introduction of the penicillinase-resistant antibiotic methicillin was quickly followed by the emergence of methicillin-resistant S. aureus (MRSA) clones [11,12,13]. S. aureus is capable of acquiring resistance to virtually every antibiotic that has entered clinical use [14, 15], including recently developed agents like daptomycin and linezolid [16, 17] and the last resort antibiotic vancomycin [18, 19]. In 2017, the World Health Organization has listed vancomycin-intermediate and vancomycin-resistant MRSA among the high priority pathogens for research and development of new antibiotics .
S. aureus’s ability to spread worldwide and to cause outbreaks in both hospitals and the community [21, 22] has fostered the study of its global epidemiology [3, 15, 23,24,25]. Some lineages are very prevalent worldwide (e.g. CC5 and CC8) , whereas others have a more localised spreading range, like the CC5 ST612 clone, which has been found only in South Africa and Australia [24, 26]. MRSA prevalence is also highly geographically variable, ranging from < 1% in some Northern European countries to > 50% in some American and Asian countries, with livestock-associated MRSA disseminating in the last two decades . Newly emerging highly pathogenic and pandemic clones have also been globally characterised [27, 28] and are often the results of recombination events as in the case of the ST239-SCCmecIII clone [25, 27, 29]. S. aureus investigations have however often underestimated the importance of non-MRSA clones, usually considering only hypervirulent or specifically relevant methicillin-sensitive S. aureus (MSSA) lineages , even though MSSA is the most common cause of surgical site infection [30, 31] and one of the major nosocomial pathogens .
Untargeted profiling of the entire S. aureus population in a given site or area is as important as its global epidemiology, and it is crucial for surveillance and prevention of local outbreaks. Some studies have for instance unbiasedly assessed the local epidemiology of nosocomial S. aureus, suggesting that this pathogen is only rarely transmitted from nurses to hospitalised patients in presence of adequate infection prevention measures  and that the community acts as major source of nosocomial MRSA . Studies surveying the whole S. aureus population in hospitals have however focused on single aspects, like the diversity of the population, its virulence and resistance traits, and its transmission in presence of an outbreak [34,35,36,37,38,39] or in non-emergency conditions [40,41,42]. Despite the large body of researches on S. aureus, studies addressing a whole S. aureus infective population at a given site through whole genome sequencing to simultaneously look at the epidemiology, phylogenetic reconstruction, genomic characterisation, and transmission pathways of infective clones are currently limited . Expanding these types of studies will be crucial for an in-depth global monitoring of S. aureus.
Here we report an in-depth epidemiological and genomic investigation of S. aureus infections in a paediatric hospital in Italy. With a whole-genome sequencing approach, we reconstructed the phylogenies of the clones in the cohort, characterised known clones and variants, screened for resistance and virulence genes, and tested for the presence of an outbreak. This allowed us to appreciate the high diversity of the S. aureus community, with 80 different lineages, variability of the resistance cassettes, and uneven conservation of various antigens previously clinically tested for vaccine development. We further report an increased prevalence of highly resistant and lowly virulent clones in chronic infections, and the rise of a newly emerging clone already reported in other hospitals. Overall, our results highlight the complexity of S. aureus epidemiology and advocate the need for wider genome-based analysis.
Materials and methods
Sample collection and S. aureus isolation
Samples were collected at Anna Meyer Children’s University Hospital (Florence, Italy) from 160 patients from January 2013 to December 2015. Metadata were also collected (Additional file 1: Table S1). We analysed samples obtained from the most common sites of infection for S. aureus, namely airways (bronchial aspirates, sputum or oropharyngeal and nasal swabs) or from soft-tissue and skin lesions. All samples were processed for the detection of bacteria using selective (Mannitol Salt Agar 2, bioMérieux) and chromogenic culture media for MRSA (BBL™ CHROMagar™ MRSA II, Becton Dickinson). In order to confirm species-level identification, mass spectrometry analysis was performed using matrix-assisted laser desorption/ionisation time of flight (MALDI-TOF) (VITEK® MS, bioMérieux). Antibiotic susceptibility was evaluated using the automated system VITEK®2 (bioMerieux) with the card AST-P632 (see Additional file 1: Table S1 for antibiograms). All identified strains were stored at − 80 °C for the following molecular analysis.
Molecular characterisation of S. aureus and MRSA isolates
DNA extraction was performed from pure S. aureus cultures after 24 h of incubation at 37 °C on Columbia agar + 5% sheep blood (bioMérieux) using QIAamp DNA Mini Kit (cat. num. 51306, QIAGEN, Netherlands) according to the manufacturer’s specifications. DNA was purified using Agencourt AMPure XP (Beckman Coulter, California, USA) according to the manufacturer’s specifications. Extracted DNA was stored at − 20 °C for further analyses.
In order to determine the potential virulence of SA/MRSA strains, a specific PCR assay for the presence of the gene (lukS-lukF) encoding for the Panton-Valentine Leukocidin (PVL) was set up following a previously published protocol . The mecA gene and other loci of the SCCmec cassette were analysed using different multiplex PCR. The protocol suggested by Milheirico et al.  has been used as a screening test for most frequent SCCmec cassettes types (types I, II, III, IV, V, and VI) and then confirmed with other methods in equivocal cases [45,46,47,48].
PCR-based multilocus sequence typing (MLST) was carried out with 25 μl reaction volumes containing 2 μl of chromosomal DNA, 20 μM of each primer, 1 U of Taq DNA polymerase (Super AB Taq, AB analitica), 2.5 μl of 10× PCR buffer (supplied with the Taq polymerase), 1.5 μM MgCl2, and 250 μM each deoxynucleoside triphosphates. The PCR was performed with an initial 5-min denaturation at 95 °C, followed by 30 cycles of annealing at 55 °C for 1 min, extension at 72 °C for 1 min, and denaturation at 95 °C for 1 min, followed by a final extension step of 72 °C for 5 min. The amplified products were purified and then amplified with the BigDye® Terminator v3.1 Cycle Sequencing Kits (Applied Biosystem) with the primers used in the initial PCR amplification. The sequences of both strands were determined with an ABI Prism 310 DNA sequencer. Isolates with the same ST have identical sequences at all seven MLST loci.
Isolates sequencing and data pre-processing
DNA libraries were prepared with Nextera XT DNA Library Preparation Kit (Illumina, California, USA). Quality control was performed with Caliper LabChip GX (Perkin Elmer) prior to shotgun sequencing with MiSeq (Illumina, California, USA), with an expected sequencing depth of 260 Mb/library (expected coverage > 80×). One hundred twenty-nine million reads were generated (704 thousand reads/sample s.d. 349 thousand).
Sequences were pre-processed by removing low-quality (mean quality lower than 25) or low-complexity reads, reads mapping to human genome or to large and small ribosomal units of bacteria, fungi and human, and known contaminants (e.g. phiX174, Illumina spike-in). All genomes are available at the NCBI Sequence Read Archive (BioProject accession number PRJNA400143).
Genome assembly and annotation
Pre-processed reads were de novo assembled using SPAdes version 3.6.1  and discarding contigs shorter than 1000 nt. We selected for our analysis only reconstructed genomes with an N50 > 50,000. We obtained high-quality genomes (N50 > 50,000 and less than 150 contigs) for 135 of the 160 patients enrolled. Genomes belonging to the remaining 25 patients were excluded from further analyses. Genomes were annotated with Prokka version 1.11  using default parameters and adding --addgenes and --usegenus options.
Genome alignment/phylogenetic analysis
The sets of 1464 concatenated genes used as input for constructing whole cohort (Fig. 1) and strain (Fig. 2) phylogenetic trees were generated using Roary version 3.4.2 . Maximum likelihood trees were inferred with RAxML version 8.0.26  using a GTR replacement model with four discrete categories of Gamma. Support at nodes was estimated using 100 bootstrap pseudo-replicates (option “-f a”). The phylogenetic tree in Additional file 2: Figure S1 was inferred using the presence-absence binary matrix of the core and accessory genes computed with Roary version 3.4.2  in RAxML version 8.0.26  with option “-m BINGAMMA”. Phylogenetic analyses were conducted using only one single isolate per patient; when multiple isolates from different timepoints of the same patient were available, the reconstructed genome with the highest N50 and the lowest number of contigs was selected. In most cases (n = 30), patients maintained the same ST over time; in discrepant cases (n = 2), we selected the most prevalent clone.
In silico sequence type (ST), SCCmec, and spa-type identification
In order to assign SCCmec type also to equivocal cases and to confirm PCR-based SCCmec typing, the same set of primers  and other primer sets [53, 54] were mapped to reconstructed genomes by BLAST . In most cases, the two methods were consistent. In discordant cases, PCR was repeated. Sequence typing and spa-typing were conducted using MetaMLST  and the DNAGear software  respectively.
Many isolates were not assigned a spa-type because of the limitations of short-read shotgun sequencing in repeated regions, which cause problems in genome assembly.
Virulence factors and resistance gene analysis
Selected virulence factors and resistance genes (as in ) were searched for by mapping reference genes (Additional file 3: Table S2) to all reconstructed genomes with BLAST  with the following parameters [−evalue 1e−10 -perc_identity 90 -gapopen 5 -gapextend 5] with a match > 75%. Virulence genes to be searched for were selected on the basis of a careful literature review for their clinical relevance [59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84].
Analysis of genes with available vaccine targets
Genes of interest were identified as those S. aureus vaccine candidates that had already entered clinical trials (according to http://clinicaltrials.gov as of January 2018), and those candidates that showed promising results in preclinical trials. For each genome, we extracted the reference sequences using BLAST  with default parameters. Extracted genes were pairwise globally aligned with the reference and evaluated for synonymous and non-synonymous single-nucleotide variation (SNVs), insertions, and/or deletions.
Bayesian divergence estimates
We estimated divergence times of ST1 SCCmecIV t127 PVL− clones using BEAST2  and the core genome (core genes = 1464). We defined the best fitting model priors by testing the combination of three clock models (uncorrelated relaxed exponential, uncorrelated relaxed lognormal, and strict), three demographic models (birth-death, coalescent Bayesian skyline, and constant), and two substitution models (HKY - Hasegawa, Kishino, Yano and generalised time reversible). Bayesian Markov chain Monte Carlo were run for 500 Mio. generations and sampled every 1000 generations. We chose the combination of models that resulted in the highest Bayes factor after parameter correction using AICM in Tracer (see Additional file 4: Table S3).
Associations between STs/virulence genes/antibiotic resistance markers and sample/operative unit types were found by performing Fisher’s exact test between the class of interest and the remaining set of samples.
Results and discussion
We investigated the epidemiology and the whole-genome genetics of S. aureus isolated from multiple operative units of the same paediatric hospital in Italy (Meyer’s Children Hospital, Florence). Two hundred thirty-four S. aureus isolates from 160 patients were retrieved from diverse clinical specimens, tested for antibiotic susceptibility, and subjected to whole-genome sequencing (see Materials and methods). The study produced 184 high-quality reconstructed S. aureus genomes with a N50 larger than 50,000 and less than 250 contigs (Additional file 1: Table S1). Downstream analyses are focused on the 135 high-quality strains recovered from distinct patients.
Genome sequencing highlights the presence of common clonal complexes and five newly sequenced clones
We first performed a whole-genome phylogenetic analysis to investigate the population structure of S. aureus in our cohort. The phylogeny was built using one isolate for each patient (n = 135) and using the 1464 core genes representing a core genome of > 1.19 M bases (see Materials and methods and Fig. 1). The genomic diversity of S. aureus is highlighted by the relatively large number of accessory genes even in a limited cohort of clinical isolates (n = 6909 from a pangenome of 8373 (Additional file 2: Figure S2), in concordance with a recent study based on the pangenome of 64 strains from different ecological niches . The gene presence/absence phylogenetic model considering both core and genes confirmed the structure of the one built on the core genome alone, with however a slightly higher strain-diversity for isolates belonging to the same ST (Additional file 2: Figure S1). Despite this diversity, we found the presence of a reduced set of closely related strains in the cohort (Fig. 1) mostly associated with distinct multilocus sequence typing clones (STs)  (see Materials and methods). We identified a total of 29 different STs, with five of them—ST228, ST22, ST5, ST121, and ST1—found in at least 12 patients (Table 1 and Additional file 1: Table S1) with evidence of ST replacement in only one patient (Patient 091 switching from ST228 to ST22) of the 32 patients sampled at multiple timepoints. This longitudinal strain consistency was confirmed by whole-genome analysis (mean intra-patient variability = 56.42 SNVs), for which the replacing event in Patient 091 accounted for 6238 SNVs between the 2013 and 2016 isolates, 0.22% of the genome. The 29 identified STs belong to 14 clonal complexes (CCs), with the five most prevalent CCs (CC5, CC22, CC8, CC1, and CC121) comprising more than 60% of the isolates. Spa-typing  further refined the typing resolution: we found 44 distinct spa-types (Additional file 1: Table S1), with t001, t002, t008, and t127 being the most prevalent (i.e. present in > 4 isolates, Table 1). We also investigated the presence of the Panton-Valentine Leukocidin (PVL), a two-component prophage virulence factor allowing S. aureus to escape from the host immune system, that was found in 27.4% of the samples (Additional file 1: Table S1).
According to both antibiotic susceptibility testing (oxacillin and cefoxitin susceptibility, Additional file 1: Table S1) and genome analysis (presence of the SCCmec cassette, see Materials and methods), 63.7% of the isolates were classified as methicillin-resistant S. aureus (MRSA). Most strains (n = 54) belonged to SCCmecIV; type I cassettes were also abundant (n = 19), whereas cassettes type V (n = 8) and II (n = 1) were less represented. Methicillin resistance was unevenly distributed across the phylogenetic tree (Fig. 1) and partially independent from the STs. All CC1 isolates (n = 14, ST1 and ST772) were MRSA, and so were the isolates belonging to CC5 (n = 30, ST5 and ST228) and CC22 (n = 16, ST22 and ST1327). All CC121 (n = 12, ST121) and CC10 (n = 3, ST10 and ST1162) isolates were instead methicillin-sensitive (MSSA), and other clonal complexes (CC8, CC30, CC45) showed balanced proportions of sensitive and resistant strains. SCCmecI (n = 19) was the most CC-specific cassette, as it was found almost exclusively in CC5 isolates (ST5 and ST228), with the exception of one ST15 and one ST8 isolates, while neither SCCmecIV nor SCCmecV were associated with specific STs.
For five of the recovered STs, namely ST241, ST942, ST1162, ST1327, and ST1866, no sequenced genome is publicly available (as genomes of S. aureus in RefSeq  version 2017 ). Although a large number of S. aureus genome sequences are available in NCBI, these are biased toward a limited set of clinically relevant STs [43, 90], with many others being neglected. This underrepresentation of less-pathogenic or less-known strains may lead to a poor understanding of the host–pathogen interactions at the genomic level, and to an underestimation of emerging or re-emerging pathogenic strains [25, 43].
Co-presence of local, global, animal-associated, and hypervirulent clones
We combined the four characterisation methods (MLST, SCCmec-, spa-, and PVL typing) to identify specific known clones in the cohort, yielding 80 different lineages. The most prevalent were the South German/Italian ST228-SCCmecI clone (n = 16, 11.85%) and the E-MRSA-15 ST22-SCCmecIV clone (n = 13, 9.63%), followed by the USA400 ST1-SCCmecIV t127 (n = 11, 8.15%) clone, the USA800 paediatric clone ST5-SCCmecIV t002 (n = 10, 7.41%), and the USA500 E-MRSA-2/6 clone ST8-SCCmecIV t008 PVL− (n = 4, 2.96%) (Additional file 5: Table S4). Several other clones, including the highly virulent USA300 ST8-SCCmecIV PVL+ clone (n = 2, 1.48%), were also found, confirming a heterogeneous clone composition in Italian hospitals [91, 92]. Surprisingly, we did not isolate any ST80, the most prevalent community-associated MRSA clone in Europe .
We moreover identified two isolates (1.48%) belonging to the livestock-associated MRSA (LA-MRSA) ST398 clone [24, 94] (Table 1). This clone has already been reported in patients that had regular exposure to livestock in several countries [24, 95, 96] including Italy [97,98,99], but our results and other reports [96, 100,101,102] of infections in non-exposed subjects suggest that the within-subject transmission for these clones is not rare. Similar conclusions can be drawn for another LA-MRSA, namely ST97 (n = 2, 1.48%, Table 1), which is the leading cause of bovine mastitis, but is only rarely reported in humans [103,104,105,106]. This growing incidence of LA-MRSA strains (n = 4, 2.96% in our cohort) causing zoonotic infections highlights the existence of underestimated reservoirs of S. aureus strains that could become epidemic [28, 107, 108].
One isolate was assigned to ST395, which is an unusual strain unable to exchange DNA via bacteriophages with other S. aureus strains because of a modification in the wall teichoic acid (WTA) [109, 110]. The same modification, however, enables ST395 to exchange DNA with coagulase-negative staphylococci (CoNS) , making it particularly prone to exchange SCCmec elements and others with other commonly found staphylococci, e.g. S. epidermidis.
Genomic signatures of chronic versus acute S. aureus infections
In order to investigate the potential association of clones and antibiotic resistance with specific hospital operative units, we cross-checked the prevalence of SCCmec types, STs, and PVL+ clones with both OUs and sample types (see Materials and methods). Strains from the cystic fibrosis (CF, n = 76) unit were positively associated with the presence of SCCmecI (n = 19, ten from CF unit; p value = 0.03), a cassette known to be hospital-associated [111, 112]. Strains from the same unit were also associated with ST1 (n = 12, seven from CF unit; p value = 0.04), whereas we noted a reduced prevalence of the PVL genes (n = 37, only two from CF unit; p value = 0.0002) and of ST121 (n = 12, none from CF unit; p value = 0.02). This reflects the relatively attenuated virulence which is a well-known phenomenon in long-term S. aureus infections [113,114,115,116]. Similarly, sputum samples (n = 33; 88.7% from CF unit) were associated with ST228 (n = 16, nine from sputum; p value = 0.004) and SCCmecI (n = 19, 11 from sputum; p value = 0.0008), and negatively correlated with PVL (n = 37, only two from sputum; p value = 0.001). The high correlation of ST228 with lung isolates and specifically with CF has already been observed in Spain . A similar pattern of increased resistance and lowered virulence has been observed for another sample type linked with long-term lung infections, namely broncho-aspiration material (n = 23; 78.2% from intensive care unit). Strains from this sample type were associated to SCCmecIV (n = 54, 14 from broncho-aspiration material; p value = 0.008), and with PVL− (n = 98, 23 from broncho-aspiration material; p value = 0.0005) and MRSA clones (n = 83, 21 from broncho-aspiration material; p value = 0.002), highlighting once again the loss of virulence and the acquisition of resistance in long-term lung infections [113,114,115,116].
On the contrary, patients from both emergency room (n = 5) and the infectious disease unit (n = 15) show an overrepresentation of PVL+ clones (n = 37, four from emergency room and nine from infectious diseases; p values = 0.02 and 0.005, respectively), indicative of acute rather than chronic infections. Lesion swabs (n = 31) are strongly associated with MSSA (n = 49, 31 from lesion swabs; p value = 3e−08). This sample type was also associated to the hypervirulent ST121 clone [117, 118] (n = 12, 11 from lesion swabs; p value = 2e−05) and to the presence of the PVL (n = 37, 14 from lesion swabs; p value = 3e−07), suggesting that in our cohort skin and soft tissue infections (SSTIs) are predominantly caused by hypervirulent MSSA strains. Lesion swabs from children in care at the infectious diseases unit (n = 12, 80% of the samples from this operative unit) are also characterised by high prevalence of the virulent ST45 clone [119, 120] (n = 8, three from lesion swabs; p value = 0.04) that is known to be associated with SSTIs [121,122,123,124]. The expected  association between PVL (n = 37) and ST121 (n = 12, nine PVL+; p value = 0.001) and ST30 (n = 7, five PVL+; p value = 0.003) supports once again the observed increased virulence of these STs [117, 118, 126, 127], which is partially in conflict with the hypothesis of lesion colonisation by commensal strains present in the skin microbiome [128, 129].
Discovery of novel variants of SCCmecIV with kanamycin, trimethoprim, and bleomycin resistance
We next investigated the specific genetic variants of the four types of SCCmec cassettes identified and discussed above. This is relevant because the epidemiology of this genetic element is disentangled by that of the rest of the genome by virtue of its high horizontal mobility [130, 131]. Moreover, the SCCmec can host genes encoding not only for resistance to beta-lactams [132, 133], but also for other antibiotic resistances or virulence factors .
More than a half of the MRSA isolates in our collection (n = 86) carried SCCmecIV (62.8%). This cassette type has spread widely in the last decades, often substituting the previously more prevalent nosocomial SCCmec types I and II [24, 134], and it is now common especially in European clinical isolates [24, 92]. Another cassette that has spread in recent years following a similar path is SCCmecV [134, 135], the third most prevalent cassette type in our cohort (10.5% of the MRSA isolates) after the more traditionally hospital-associated SCCmecI [24, 112] (22.1% of the MRSA isolates). We moreover isolated one MRSA carrying SCCmecII, which is widely diffused in the USA but only rarely found in Italy/Europe [25, 136]. Consistently, the SCCmecII isolate was recovered from Patient 115, which is consistent with the personal history of the patient. For two isolates, it was not possible to classify the cassette neither with PCR nor with in silico PCR using standard primers .
By aligning reconstructed SCCmec with reference cassettes (see Materials and methods), we observed a certain degree of variability inside the same cassette type, specifically in type IV (Fig. 3). Subtypes IVa, IVb, and IVc were identified, with some SCCmec elements showing insertions. Two cassettes in particular were not consistent with the already described subtypes: the SCCmec type IVc carried by MF062, which was enriched with genes for kanamycin  and bleomycin [138, 139] resistance, and the type IVa carried by MR090 that showed insertion of genes involved in resistance to trimethoprim [140, 141] (Fig. 3).
Non-SCCmec resistance profiles show different patterns in chronic and acute infections
S. aureus can easily acquire a number of resistances, including those to the last resort antibiotics vancomycin [142, 143] and daptomycin . According to results presented in previous paragraphs and elsewhere , resistances can occur by gene acquisition in the SCCmec cassette. Most resistances are however encoded by genes that are found in other parts of the genome or that have been horizontally transferred through different genetic elements . Given the high importance of multi-drug resistance in S. aureus , we therefore tested the presence or absence of specific resistance genes in our cohort  (Fig. 4 and Additional file 3: Table S2). Consistently with previous literature , most of the isolates tested positive for blaZ (81.5%), responsible for penicillin resistance (96.3% concordance with antibiotic susceptibility test, as per presence of the pbp and/or mecA genes). No isolates were found positive for genes encoding resistance to vancomycin (van, 100% concordance with antibiotic susceptibility test) and to fusidic acid (fusB and far, 94.1% concordance with antibiotic susceptibility test). Antibiotic resistances were sometimes associated with specific CCs, as for the increased representation of aacA.aphD (gentamicin resistance, 92.6% concordance with antibiotic susceptibility test) and ermA (erythromycin resistance, phenotypic resistance not tested) in CC5 isolates, whose genomes tended to lack instead the blaZ gene (penicillin resistance) (Fig. 4). Overall, two isolates from acute skin infections were negative for all the resistance genes tested, while six CF and intensive care unit isolates were positive for six (33.3%) of them. This pattern of increased resistance in long-term infections, together with their observed reduced virulence, completes the scenario of reduced virulence and increased resistance that has been observed in this and previous studies [113,114,115,116].
Emergence and disease-associated diversity of clinically relevant virulence factors
S. aureus has a large repertoire of virulence genes, and it is able to evade the host immune system through a variety of strategies. Some of the genes usually involved in immune evasion were present in almost all our isolates (Fig. 4 and Additional file 3: Table S2). These include genes encoding the phenol-soluble modulin alpha and beta and the delta-haemolysin Hld, responsible for leukocytes and erythrocyte lysis respectively ; the immunoglobulin-binding protein Sbi that inhibits IgG and IgA [61, 62]; and some genes part of the GIɑ genomic island (ssl6 and ssl9).
Other genes belonging to the immune evasion island IEC2 were present in many but not all isolates, for example, the one encoding for the antiplatelet extracellular fibrinogen binding protein Efb [63, 64] and those encoding various haemolysins (hla, hlg) [59, 60] (Fig. 4 and Additional file 3: Table S2). In addition to the 27.4% prevalence of the lukF and lukS PVL genes discussed above, one sample (MR029, from emergency room) was positive for the epidermal cell differentiation inhibitor Edin, which has been found to promote the translocation of S. aureus into the bloodstream . One of the two USA300 isolates (MR047, from nasal swab) tested positive for the arginine catabolic mobile element (ACME), another important virulence factor (gene arcA) that has been shown to be responsible for the increased pathogenicity of S. aureus and specifically of USA300 clones [66, 67].
Many virulence genes were associated to specific STs (Fig. 4 and Additional file 3: Table S2). ST22 (n = 15), for instance, was associated with the toxic shock syndrome toxin TSST-1 (n = 8, three from ST22; p value = 0.04; present in 20% of the ST22 clones) [60, 68, 69], other pyrogenic toxin superantigens known as staphylococcal enterotoxins (SEs, mean n = 34.9 ± 28.1 s.d.; p value < 0.02 for seg, sei, sem, sen, seo, present on average in 86.7% of ST22 and 49.5% of non-ST22), and various ssl immune evasion genes (mean n = 57.4 ± 40.2 s.d.; p value < 0.01 for ssl1, ssl3, ssl4, ssl7, ssl11, ssl12, present on average in 93.3% of ST22 and 21.1% of non-ST22) [70, 71]. ST22-IV EMRSA-15 clones positive for tst1 are usually described as “Middle Eastern variant” [72,73,74], but a high prevalence in an Italian neonatal intensive care unit  and pre-school children living in Palermo, Italy , has been observed. Authors suggested that the Middle Eastern clone might be more widely spread than estimated and might have diffused in the Mediterranean populations as a community-acquired MRSA [75, 76], as suggested by our analysis. TSST-1 is responsible for an increased pyrogenic, emetic, and superantigen activity, together with SEs (10627489; 11544350). SEs (mean n = 34.9 ± 28.1 s.d.) were associated with all “virulent” STs, such as ST5 (n = 15, p value < 0.04 for sed, seg, sei, sej, sem, sen, seo, sep, present on average in 79.2% of ST5 and 32.3% of non-ST5), ST45 (n = 8; p value < 0.02 for sec, seg, sei, sel, sem, seo, present on average in 98.2% of ST45 and 38.9% of non-ST45), ST121 (n = 12; p value < 0.02 for seb, seg, sei, sem, sen, seo, present on average in 86.1% of ST121 and 41.7% of non-ST121), and—to a lower extent—ST30 (n = 7; p value < 0.02 for sei, sem, sen, present on average in 100% of ST30 and 50% of non-ST30).
The hypervirulent ST121 MSSA isolates obtained from lesion swabs (n = 12) were instead associated with the genes encoding for the exfoliative toxins Eta and Etb (n = 3 from ST121 swabs, and n = 0 for non-ST121, p value = 0.0006 for both genes), responsible for the skin manifestations of bullous impetigo and Staphylococcal scalded skin syndrome [77,78,79], the gene bbp (n = 12 from ST121, n = 7 from non-ST121; p value = 1.35e−08) that interacts with the extracellular matrix bone sialoprotein and contributes to staphylococcal arthritis and osteomyelitis , and the immune evasion gene ecb (n = 12 from ST121, n = 36 from non-ST121; p value = 1.51e06), which is required for the persistence of S. aureus in host tissues and the formation of abscesses . The latter was also present in all and only the isolates belonging to ST1, ST7, ST10, ST15, ST30, ST34, and ST398, suggesting a strong dependence on ST (Fig. 4 and Additional file 3: Table S2).
Isolates retrieved from sputum samples of CF patients (n = 38) showed a positive association with the adhesin-encoding genes sdrD (n = 34 from CF, n = 69 from non-CF; p value = 0.03) and sdrE (n = 27 from CF, n = 48 from non-CF; p value = 0.03), and a negative association with bbp (n = 1 from CF, n = 18 from non-CF; p value = 0.01), contrary to samples from the infectious disease unit (n = 15, four positive for bbp gene). This finding is consistent with the increased need for adhesins in chronic lung infections [82, 83, 116], including in CF .
Conservation of genes encoding vaccine candidates
Unlike other bacterial infections, prior exposure to S. aureus does not seem to provide protective immunity ; therefore, vaccines are an attractive yet challenging option to prevent disease. Researchers have long attempted to produce an effective vaccine against S. aureus, but even though few have proved promising in animal models, the two vaccines so far tested in efficacy clinical trials have failed [147,148,149,150,151,152,153]. Since the main issue is the polymorphic expression of S. aureus surface antigens and the redundancy of its virulence proteins [147, 154, 155], we tested the prevalence and conservation of a number of genes encoding vaccine candidates described in the literature (Table 2).
Among antigens that have been proposed as targets for vaccine development, the alpha haemolysin toxin gene hla [147, 156, 157] and the genes coding for capsular biosynthesis cap5 and cap8 [150, 151] are highly prevalent in our cohort (91.9% and 97.8% of the isolates respectively). Nevertheless, these genes showed a larger degree of variability compared to the others we considered, which may explain the poor results obtained in clinical trials [147, 150, 151, 156, 157]. Other genes that code for proteins used alone or in combination in vaccine formulations, such as the virulence determinant SpA  and the fibronectin-binding protein ClfA [159,160,161], are present in most of our strain collection. In some of these genes, indels are prevalent (> 90%, Table 2), but they are frequently found in repeated regions that may not critically impact the protein structure, as in the case of the spa gene.
Vaccines have also been proposed for S. aureus strains with specific characteristics. For instance, targeting the toxicity determinant TSST-1 (5.9% prevalence of tst1) [162, 163] or the PVL proteins LukF-LukS (27.4% prevalence of lukF-lukS) [164, 165] aims at selectively preventing the most virulent or lethal infections. In our cohort, despite their low prevalence, both tst1 and the PVL genes were conserved at 99%, except for a few isolates that had indels in the latter (Table 2). The gamma-haemolysins HlgAB and HlgCB genes [164, 165] were instead highly prevalent (97.8–100%) and quite conserved (69.6–94.8%). The opposite approach is targeting genes with a lower virulence profile, which may be more prevalent and conserved than those coding for highly toxic factors. Among them, the genes encoding for the manganese uptake receptor (mntC) [159,160,161] and for the iron acquisition factor (isdB) [152, 166], which are indeed present in all or all but one the isolates of our cohort. Non-synonymous mutations are rare in mntC (20.7% of the isolates, with only one non-synonymous SNV), and, whenever not affected by indels that may or may not affect the protein structure, also the isdB gene is highly conserved (> 99% identity, Table 2).
Finally, we also analysed the conservation of csa1A, csa1B, fhuD2, and esxA, genes recently described as being promising vaccine candidates in preclinical studies [167, 168]. The two genes encoding for the conserved antigen Csa (csa1A and csa1B) are present in 51.9% and 26.7% of the isolates respectively and are conserved in only a fraction of the cases (Table 2). By contrast, the iron uptake gene fhuD2 is present in all isolates, with a maximum of 1% non-synonymous variation in sequence (Table 2). Also the genes encoding for the ESAT-6-like secretion system (esxA, esxB, esxC, esxD) are well represented in the cohort, but only esxA is present in all but one isolate and has no non-synonymous mutations in 85.1% of the isolates (Table 2). Therefore, on the basis of their conservation, both FhuD2 and EsxA appear to be promising targets for vaccine formulations.
Phylogenetics of specific STs highlights the aggressive spread of a novel independently acquired ST1 clone
We investigated the hypothesis that some of the prevalent STs could be hospital-associated clones. We estimated the ST phylogenies using a whole-genome maximum likelihood approach (see Materials and methods). In most cases, we observed that isolates in our cohort, despite sharing the same ST, SCCmec, and spa types, were not monophyletic subtrees when considering external reference genomes for the same STs. This is the case, for example, of the ST228 and ST5 clones (Fig. 2). This suggests independent acquisition of the clones and no evidence of transmission among the selected hospitalised patients, while person-to-person transmission from healthy carriers or non-selected patients cannot be ruled out [21, 22]. Only two ST121 MSSA isolates were found to be almost identical and both were retrieved in the same time window from patients 096 and 098 (8 SNVs). For ST1, instead, all but two isolates belonged to the same sub-lineage, typed as SCCmecIV t127 PVL−.
We further estimated divergence times for all the 16 isolates belonging to the ST1 SCCmecIV t127 PVL− clone, including those obtained from earlier or later time points of the same patients. We used a Bayesian approach  (see Materials and methods) integrating all the reference genomes publicly available for ST1 and the two ST1 SCCmecV isolates from our cohort (Additional file 4: Table S3). These analyses were performed to test the hypothesis that all ST1 SCCmecIV t127 belong to a clone specific of Meyer’s hospital. The relaxed exponential clock model with constant coalescent prior and GTR substitution model resulted to the most appropriate model (Additional file 6: Table S5). This model estimated that Meyer’s clone has emerged approximately 6 to 28 years ago as a specific branch of the ST1 tree, which has been estimated to be 26–160 years old (Fig. 5). However, age of Meyer’s clone does not match with the time of emergence of the clone in the hospital. Moreover, an isolate obtained in a recent study investigating the spread of a ST1 SCCmecIV t127 clone in Irish hospitals  and carrying a virulence and resistance profile very close to the one of our cohort (differences in gene presence: 2/79 and 0/18 respectively) is phylogenetically rooted inside Meyer’s cluster (161 SNVs intra-cluster; 412 SNVs inter-cluster). These two findings suggest that ST1 SCCmecIV t127 is not specific of the Meyer Children’s hospital but might represent a newly arising community clone that is now spreading in the nosocomial environment of different countries [169, 170].
In this study, we investigated the epidemiology of S. aureus in different operative units of Anne Meyer’s Children’s University Hospital (Florence, Italy) over a timespan of 3 years by whole genome isolate sequencing. Our analyses highlighted a high diversity of STs, SCCmec, and spa-types, resulting into a wide number of clones. Some of these clones had been previously described in the literature as livestock-associated, and we described them in non-exposed children thus supporting the spreading of such clones in the non-at-risk community. We moreover described the presence of hypervirulent and geographically unusual clones, and of five STs for which no sequenced genome was available in public databases. Our refined analysis of the SCCmec cassettes highlighted the presence of further resistances and diversity within the same cassette type. On the contrary, when considering single infection types or specific STs or clones as it is usual in S. aureus epidemiological studies, the genomic diversity was limited, with an increased pattern of resistance genes in chronic patients and a larger number of virulence factors in acute infections. Altogether, these observations shed more light on the complexity of S. aureus epidemiology and on the need for a more unbiased survey of the commensal and pathogenic S. aureus community, to avoid the misrepresentation of specific genomic traits.
Whole-genome-based routine surveillance of S. aureus and other hospital-related pathogens would further allow to get a more unbiased idea of the rising clones and better informing clinical practices, which usually focused on the most dangerous or well-known strains. Performing such epidemiological studies as soon as a new putative nosocomial clone arises could allow us to conclude whether the new clone has arisen in that very hospital or it is a recent sub-clone spreading also in the non-hospitalised population and therefore more frequently isolated also in the clinics. These wider-focus studies would not only allow the assessment of the epidemiology of specific pathogens and clones in the hospital setting, but also the survey of the prevalence and conservation of their virulence and resistance traits. This could lead to the identification of antigens of interest for vaccine development and of specific sub-clones representing the main burden of infection, and therefore reassessing the efforts for the discovery of new treatments.
Whole genome sequencing studies are crucial to survey the global epidemiology of infectious agents, including S. aureus, as genome-based data are reproducible and can be easily meta-analysed without the confounding of batch effects. The meta-analysis of pathogenic, commensal, and environmental S. aureus isolates could lead to a deeper knowledge of the epidemiology of this bacterium and may help in understanding how to prevent and treat infections without boosting antibiotic resistance.
Insertions and deletions
Multilocus sequence typing
Methicillin-resistant S. aureus
Methicillin-sensitive S. aureus
Staphylococcal Cassette Chromosome mec
Staphylococcal protein A
Whole-genome shotgun sequencing
Wertheim HF, Melles DC, Vos MC, van Leeuwen W, van Belkum A, Verbrugh HA, Nouwen JL. The role of nasal carriage in Staphylococcus aureus infections. Lancet Infect Dis. 2005;5(12):751–62.
Mainous AG 3rd, Hueston WJ, Everett CJ, Diaz VA. Nasal carriage of Staphylococcus aureus and methicillin-resistant S aureus in the United States, 2001-2002. Ann Fam Med. 2006;4(2):132–7.
Tong SY, Davis JS, Eichenberger E, Holland TL, Fowler VG Jr. Staphylococcus aureus infections: epidemiology, pathophysiology, clinical manifestations, and management. Clin Microbiol Rev. 2015;28(3):603–61.
Esposito S, Noviello S, Leone S. Epidemiology and microbiology of skin and soft tissue infections. Curr Opin Infect Dis. 2016;29(2):109–15.
Rhee Y, Aroutcheva A, Hota B, Weinstein RA, Popovich KJ. Evolving epidemiology of Staphylococcus aureus bacteremia. Infect Control Hosp Epidemiol. 2015;36(12):1417–22.
Peacock SJ, Paterson GK. Mechanisms of methicillin resistance in Staphylococcus aureus. Annu Rev Biochem. 2015;84(1):577–601.
Skinner D. Significance of bacteremia caused by Staphylococcus aureus. Arch Intern Med. 1941;68(5):851.
Rammelkamp CH, Maxon T. Resistance of Staphylococcus aureus to the action of penicillin. Exp Biol Med. 1942;51(3):386–9.
Kirby WM. Extraction of a highly potent penicillin Inactivator from penicillin resistant staphylococci. Science. 1944;99(2579):452–3.
Bondi A Jr, Dietz CC. Penicillin resistant staphylococci. Proc Soc Exp Biol Med. 1945;60:55–8.
Barber M. Methicillin-resistant staphylococci. J Clin Pathol. 1961;14:385–93.
Jevons MP, Coe AW, Parker MT. Methicillin resistance in staphylococci. Lancet. 1963;1(7287):904–7.
Parker MT, Jevons MP. A survey of methicillin resistance in Staphylococcus aureus. Postgrad Med J. 1964;40(SUPPL):170–8.
Pantosti A, Sanchini A, Monaco M. Mechanisms of antibiotic resistance in Staphylococcus aureus. Future Microbiol. 2007;2(3):323–34.
Monaco M, Pimentel de Araujo F, Cruciani M, Coccia EM, Pantosti A. Worldwide epidemiology and antibiotic resistance of Staphylococcus aureus. Curr Top Microbiol Immunol. 2017;409:21–56.
Liu C, Bayer A, Cosgrove SE, Daum RS, Fridkin SK, Gorwitz RJ, Kaplan SL, Karchmer AW, Levine DP, Murray BE, et al. Clinical practice guidelines by the infectious diseases society of america for the treatment of methicillin-resistant Staphylococcus aureus infections in adults and children. Clin Infect Dis. 2011;52(3):e18–55.
Raad I, Hanna H, Jiang Y, Dvorak T, Reitzel R, Chaiban G, Sherertz R, Hachem R. Comparative activities of daptomycin, linezolid, and tigecycline against catheter-related methicillin-resistant Staphylococcus bacteremic isolates embedded in biofilm. Antimicrob Agents Chemother. 2007;51(5):1656–60.
Kos VN, Desjardins CA, Griggs A, Cerqueira G, Van Tonder A, Holden MT, Godfrey P, Palmer KL, Bodi K, Mongodin EF, et al. Comparative genomics of vancomycin-resistant Staphylococcus aureus strains and their positions within the clade most commonly associated with Methicillin-resistant S. aureus hospital-acquired infection in the United States. MBio. 2012;3(3):e00112-12.
Hiramatsu K. Vancomycin-resistant Staphylococcus aureus: a new model of antibiotic resistance. Lancet Infect Dis. 2001;1(3):147–55.
Tacconelli E, Carrara E, Savoldi A, Harbarth S, Mendelson M, Monnet DL, Pulcini C, Kahlmeter G, Kluytmans J, Carmeli Y, et al. Discovery, research, and development of new antibiotics: the WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect Dis. 2018;18(3):318–27.
Tosas Auguet O, Stabler RA, Betley J, Preston MD, Dhaliwal M, Gaunt M, Ioannou A, Desai N, Karadag T, Batra R, et al. Frequent undetected ward-based methicillin-resistant Staphylococcus aureus transmission linked to patient sharing between hospitals. Clin Infect Dis. 2018;66(6):840–8.
Coll F, Harrison EM, Toleman MS, Reuter S, Raven KE, Blane B, Palmer B, Kappeler ARM, Brown NM, Torok ME, et al. Longitudinal genomic surveillance of MRSA in the UK reveals transmission patterns in hospitals and the community. Sci Transl Med. 2017;9(413):eaak9745.
Voss A, Doebbeling BN. The worldwide prevalence of methicillin-resistant Staphylococcus aureus. Int J Antimicrob Agents. 1995;5(2):101–6.
Stefani S, Chung DR, Lindsay JA, Friedrich AW, Kearns AM, Westh H, Mackenzie FM. Methicillin-resistant Staphylococcus aureus (MRSA): global epidemiology and harmonisation of typing methods. Int J Antimicrob Agents. 2012;39(4):273–82.
Chambers HF, Deleo FR. Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat Rev Microbiol. 2009;7(9):629–41.
Jansen van Rensburg MJ, Eliya Madikane V, Whitelaw A, Chachage M, Haffejee S, Gay Elisha B. The dominant methicillin-resistant Staphylococcus aureus clone from hospitals in Cape Town has an unusual genotype: ST612. Clin Microbiol Infection. 2011;17(5):785–92.
Harris SR, Feil EJ, Holden MT, Quail MA, Nickerson EK, Chantratita N, Gardete S, Tavares A, Day N, Lindsay JA, et al. Evolution of MRSA during hospital transmission and intercontinental spread. Science. 2010;327(5964):469–74.
Mediavilla JR, Chen L, Mathema B, Kreiswirth BN. Global epidemiology of community-associated methicillin resistant Staphylococcus aureus (CA-MRSA). Curr Opin Microbiol. 2012;15(5):588–95.
Deurenberg RH, Stobberingh EE. The evolution of Staphylococcus aureus. Infect Genet Evol. 2008;8(6):747–63.
Point prevalence survey of healthcare-associated infections and antimicrobial use in European long-term care facilities, May–September 2010 [https://publications.europa.eu/en/publication-detail/-/publication/d6355bbd-f651-4770-bcad-0456b1dc2515 Accessed 11 Apr 2018].
Sievert DM, Ricks P, Edwards JR, Schneider A, Patel J, Srinivasan A, Kallen A, Limbago B, Fridkin S, National Healthcare Safety Network T, et al. Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2009-2010. Infect Control Hosp Epidemiol. 2013;34(1):1–14.
Price JR, Cole K, Bexley A, Kostiou V, Eyre DW, Golubchik T, Wilson DJ, Crook DW, Walker AS, Peto TEA, et al. Transmission of Staphylococcus aureus between health-care workers, the environment, and patients in an intensive care unit: a longitudinal cohort study based on whole-genome sequencing. Lancet Infect Dis. 2017;17(2):207–14.
Prosperi M, Veras N, Azarian T, Rathore M, Nolan D, Rand K, Cook RL, Johnson J, Morris JG Jr, Salemi M. Molecular epidemiology of community-associated methicillin-resistant Staphylococcus aureus in the genomic era: a cross-sectional study. Sci Rep. 2013;3:1902.
Harris SR, Cartwright EJ, Torok ME, Holden MT, Brown NM, Ogilvy-Stuart AL, Ellington MJ, Quail MA, Bentley SD, Parkhill J, et al. Whole-genome sequencing for analysis of an outbreak of meticillin-resistant Staphylococcus aureus: a descriptive study. Lancet Infect Dis. 2013;13(2):130–6.
Stein M, Navon-Venezia S, Chmelnitsky I, Kohelet D, Schwartz O, Agmon O, Somekh E. An outbreak of new, nonmultidrug-resistant, methicillin-resistant Staphylococcus aureus strain (sccmec type iiia variant-1) in the neonatal intensive care unit transmitted by a staff member. Pediatr Infect Dis J. 2006;25(6):557–9.
Bertin ML, Vinski J, Schmitt S, Sabella C, Danziger-Isakov L, McHugh M, Procop GW, Hall G, Gordon SM, Goldfarb J. Outbreak of methicillin-resistant Staphylococcus aureus colonization and infection in a neonatal intensive care unit epidemiologically linked to a healthcare worker with chronic otitis. Infect Control Hosp Epidemiol. 2006;27(6):581–5.
Coombs GW, Van Gessel H, Pearson JC, Godsell MR, O'Brien FG, Christiansen KJ. Controlling a multicenter outbreak involving the New York/Japan methicillin-resistant Staphylococcus aureus clone. Infect Control Hosp Epidemiol. 2007;28(7):845–52.
Saiman L, Cronquist A, Wu F, Zhou J, Rubenstein D, Eisner W, Kreiswirth BN, Della-Latta P. An outbreak of methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit. Infect Control Hosp Epidemiol. 2003;24(5):317–21.
Wang JT, Chang SC, Ko WJ, Chang YY, Chen ML, Pan HJ, Luh KT. A hospital-acquired outbreak of methicillin-resistant Staphylococcus aureus infection initiated by a surgeon carrier. J Hosp Infect. 2001;47(2):104–9.
Givney R, Vickery A, Holliday A, Pegler M, Benn R. Methicillin-resistant Staphylococcus aureus in a cystic fibrosis unit. J Hosp Infect. 1997;35(1):27–36.
Blok HE, Troelstra A, Kamp-Hopmans TE, Gigengack-Baars AC, Vandenbroucke-Grauls CM, Weersink AJ, Verhoef J, Mascini EM. Role of healthcare workers in outbreaks of methicillin-resistant Staphylococcus aureus: a 10-year evaluation from a Dutch university hospital. Infect Control Hosp Epidemiol. 2003;24(9):679–85.
Tong SY, Holden MT, Nickerson EK, Cooper BS, Koser CU, Cori A, Jombart T, Cauchemez S, Fraser C, Wuthiekanun V, et al. Genome sequencing defines phylogeny and spread of methicillin-resistant Staphylococcus aureus in a high transmission setting. Genome Res. 2015;25(1):111–8.
Copin R, Shopsin B, Torres VJ. After the deluge: mining Staphylococcus aureus genomic data for clinical associations and host-pathogen interactions. Curr Opin Microbiol. 2017;41:43–50.
Lina G, Piemont Y, Godail-Gamot F, Bes M, Peter MO, Gauduchon V, Vandenesch F, Etienne J. Involvement of Panton-Valentine leukocidin-producing Staphylococcus aureus in primary skin infections and pneumonia. Clin Infect Dis. 1999;29(5):1128–32.
Milheirico C, Oliveira DC, de Lencastre H. Update to the multiplex PCR strategy for assignment of mec element types in Staphylococcus aureus. Antimicrob Agents Chemother. 2007;51(9):3374–7.
Oliveira DC, de Lencastre H. Multiplex PCR strategy for rapid identification of structural types and variants of the mec element in methicillin-resistant Staphylococcus aureus. Antimicrob Agents Chemother. 2002;46(7):2155–61.
Oliveira DC, Milheirico C, de Lencastre H. Redefining a structural variant of staphylococcal cassette chromosome mec, SCCmec type VI. Antimicrob Agents Chemother. 2006;50(10):3457–9.
Milheirico C, Oliveira DC, de Lencastre H. Multiplex PCR strategy for subtyping the staphylococcal cassette chromosome mec type IV in methicillin-resistant Staphylococcus aureus: 'SCCmec IV multiplex. J Antimicrob Chemother. 2007;60(1):42–8.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, Fookes M, Falush D, Keane JA, Parkhill J. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(22):3691–3.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Zhang K, McClure JA, Elsayed S, Louie T, Conly JM. Novel multiplex PCR assay for characterization and concomitant subtyping of staphylococcal cassette chromosome mec types I to V in methicillin-resistant Staphylococcus aureus. J Clin Microbiol. 2005;43(10):5026–33.
Boye K, Bartels MD, Andersen IS, Moller JA, Westh H. A new multiplex PCR for easy screening of methicillin-resistant Staphylococcus aureus SCCmec types I-V. Clin Microbiol Infect. 2007;13(7):725–7.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Zolfo M, Tett A, Jousson O, Donati C, Segata N. MetaMLST: multi-locus strain-level bacterial typing from metagenomic samples. Nucleic Acids Res. 2017;45(2):e7.
AL-T F, Brunel AS, Bouzinbi N, Corne P, Banuls AL, Shahbazkia HR. DNAGear—a free software for spa type identification in Staphylococcus aureus. BMC Res Notes. 2012;5:642.
Gordon NC, Price JR, Cole K, Everitt R, Morgan M, Finney J, Kearns AM, Pichon B, Young B, Wilson DJ, et al. Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing. J Clin Microbiol. 2014;52(4):1182–91.
Inoshima I, Inoshima N, Wilke GA, Powers ME, Frank KM, Wang Y, Bubeck Wardenburg J. A Staphylococcus aureus pore-forming toxin subverts the activity of ADAM10 to cause lethal infection in mice. Nat Med. 2011;17(10):1310–4.
Dinges MM, Orwin PM, Schlievert PM. Exotoxins of Staphylococcus aureus. Clin Microbiol Rev. 2000;13(1):16–34 table of contents.
Smith EJ, Visai L, Kerrigan SW, Speziale P, Foster TJ. The Sbi protein is a multifunctional immune evasion factor of Staphylococcus aureus. Infect Immun. 2011;79(9):3801–9.
Haupt K, Reuter M, van den Elsen J, Burman J, Halbich S, Richter J, Skerka C, Zipfel PF. The Staphylococcus aureus protein Sbi acts as a complement inhibitor and forms a tripartite complex with host complement factor H and C3b. PLoS Pathog. 2008;4(12):e1000250.
Shannon O, Flock JI. Extracellular fibrinogen binding protein, Efb, from Staphylococcus aureus binds to platelets and inhibits platelet aggregation. Thromb Haemost. 2004;91(4):779–89.
Shannon O, Uekotter A, Flock JI. Extracellular fibrinogen binding protein, Efb, from Staphylococcus aureus as an antiplatelet agent in vivo. Thromb Haemost. 2005;93(5):927–31.
Courjon J, Munro P, Benito Y, Visvikis O, Bouchiat C, Boyer L, Doye A, Lepidi H, Ghigo E, Lavigne JP, et al. EDIN-B promotes the translocation of Staphylococcus aureus to the bloodstream in the course of pneumonia. Toxins (Basel). 2015;7(10):4131–42.
Diep BA, Stone GG, Basuino L, Graber CJ, Miller A, des Etages SA, Jones A, Palazzolo-Ballance AM, Perdreau-Remington F, Sensabaugh GF, et al. The arginine catabolic mobile element and staphylococcal chromosomal cassette mec linkage: convergence of virulence and resistance in the USA300 clone of methicillin-resistant Staphylococcus aureus. J Infect Dis. 2008;197(11):1523–30.
Ellington MJ, Yearwood L, Ganner M, East C, Kearns AM. Distribution of the ACME-arcA gene among methicillin-resistant Staphylococcus aureus from England and Wales. J Antimicrob Chemother. 2008;61(1):73–7.
Kim J, Urban RG, Strominger JL, Wiley DC. Toxic shock syndrome toxin-1 complexed with a class II major histocompatibility molecule HLA-DR1. Science. 1994;266(5192):1870–4.
Miethke T, Duschek K, Wahl C, Heeg K, Wagner H. Pathogenesis of the toxic shock syndrome: T cell mediated lethal shock caused by the superantigen TSST-1. Eur J Immunol. 1993;23(7):1494–500.
Rooijakkers SH, van Strijp JA. Bacterial complement evasion. Mol Immunol. 2007;44(1-3):23–32.
Thammavongsa V, Kim HK, Missiakas D, Schneewind O. Staphylococcal manipulation of host immune responses. Nat Rev Microbiol. 2015;13(9):529–43.
Udo EE, Boswihi SS, Al-Sweih N. High prevalence of toxic shock syndrome toxin-producing epidemic methicillin-resistant Staphylococcus aureus 15 (EMRSA-15) strains in Kuwait hospitals. New Microbes New Infect. 2016;12:24–30.
Al Laham N, Mediavilla JR, Chen L, Abdelateef N, Elamreen FA, Ginocchio CC, Pierard D, Becker K, Kreiswirth BN. MRSA clonal complex 22 strains harboring toxic shock syndrome toxin (TSST-1) are endemic in the primary hospital in Gaza, Palestine. PLoS One. 2015;10(3):e0120008.
Biber A, Abuelaish I, Rahav G, Raz M, Cohen L, Valinsky L, Taran D, Goral A, Elhamdany A, Regev-Yochay G, et al. A typical hospital-acquired methicillin-resistant Staphylococcus aureus clone is widespread in the community in the Gaza strip. PLoS One. 2012;7(8):e42864.
Geraci DM, Bonura C, Giuffre M, Aleo A, Saporito L, Graziano G, Valenti RM, Mammina C. tst1-positive ST22-MRSA-IVa in healthy Italian preschool children. Infection. 2014;42(3):535–8.
Geraci DM, Giuffre M, Bonura C, Matranga D, Aleo A, Saporito L, Corsello G, Larsen AR, Mammina C. Methicillin-resistant Staphylococcus aureus colonization: a three-year prospective study in a neonatal intensive care unit in Italy. PLoS One. 2014;9(2):e87760.
Amagai M, Matsuyoshi N, Wang ZH, Andl C, Stanley JR. Toxin in bullous impetigo and staphylococcal scalded-skin syndrome targets desmoglein 1. Nat Med. 2000;6(11):1275–7.
Hanakawa Y, Schechter NM, Lin C, Garza L, Li H, Yamaguchi T, Fudaba Y, Nishifuji K, Sugai M, Amagai M, et al. Molecular mechanisms of blister formation in bullous impetigo and staphylococcal scalded skin syndrome. J Clin Invest. 2002;110(1):53–60.
Ladhani S. Understanding the mechanism of action of the exfoliative toxins of Staphylococcus aureus. FEMS Immunol Med Microbiol. 2003;39(2):181–9.
Tung H, Guss B, Hellman U, Persson L, Rubin K, Ryden C. A bone sialoprotein-binding protein from Staphylococcus aureus: a member of the staphylococcal Sdr family. Biochem J. 2000;345(Pt 3):611–9.
Jongerius I, von Kockritz-Blickwede M, Horsburgh MJ, Ruyken M, Nizet V, Rooijakkers SH. Staphylococcus aureus virulence is enhanced by secreted factors that block innate immune defenses. J Innate Immun. 2012;4(3):301–11.
Sanchez CJ Jr, Mende K, Beckius ML, Akers KS, Romano DR, Wenke JC, Murray CK. Biofilm formation by clinical isolates and the implications in chronic infections. BMC Infect Dis. 2013;13:47.
Paharik AE, Horswill AR. The staphylococcal biofilm: adhesins, regulation, and host response. Microbiol Spectr. 2016;4(2):VMBF-0022-2015.
Schwab UE, Wold AE, Carson JL, Leigh MW, Cheng PW, Gilligan PH, Boat TF. Increased adherence of Staphylococcus aureus from cystic fibrosis lungs to airway epithelial cells. Am Rev Respir Dis. 1993;148(2):365–9.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
Bosi E, Monk JM, Aziz RK, Fondi M, Nizet V, Palsson BO. Comparative genome-scale modelling of Staphylococcus aureus strains identifies strain-specific metabolic capabilities linked to pathogenicity. Proc Natl Acad Sci U S A. 2016;113(26):E3801–9.
Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A. 1998;95(6):3140–5.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Database issue):D61–5.
Haft DH, DiCuccio M, Badretdin A, Brover V, Chetvernin V, O'Neill K, Li W, Chitsaz F, Derbyshire MK, Gonzales NR, et al. RefSeq: an update on prokaryotic genome annotation and curation. Nucleic Acids Res. 2018;46(D1):D851–60.
Planet PJ, Narechania A, Chen L, Mathema B, Boundy S, Archer G, Kreiswirth B. Architecture of a species: Phylogenomics of Staphylococcus aureus. Trends Microbiol. 2017;25(2):153–66.
Mato R, Campanile F, Stefani S, Crisostomo MI, Santagati M, Sanches SI, de Lencastre H. Clonal types and multidrug resistance patterns of methicillin-resistant Staphylococcus aureus (MRSA) recovered in Italy during the 1990s. Microb Drug Resist. 2004;10(2):106–13.
Campanile F, Bongiorno D, Perez M, Mongelli G, Sessa L, Benvenuto S, Gona F, Participants ASS, Varaldo PE, Stefani S. Epidemiology of Staphylococcus aureus in Italy: first nationwide survey, 2012. J Glob Antimicrob Resist. 2015;3(4):247–54.
Vandenesch F, Naimi T, Enright MC, Lina G, Nimmo GR, Heffernan H, Liassine N, Bes M, Greenland T, Reverdy ME, et al. Community-acquired methicillin-resistant Staphylococcus aureus carrying Panton-Valentine leukocidin genes: worldwide emergence. Emerg Infect Dis. 2003;9(8):978–84.
Witte W, Strommenger B, Stanek C, Cuny C. Methicillin-resistant Staphylococcus aureus ST398 in humans and animals, Central Europe. Emerg Infect Dis. 2007;13(2):255–8.
van Cleef BA, Graveland H, Haenen AP, van de Giessen AW, Heederik D, Wagenaar JA, Kluytmans JA. Persistence of livestock-associated methicillin-resistant Staphylococcus aureus in field workers after short-term occupational exposure to pigs and veal calves. J Clin Microbiol. 2011;49(3):1030–3.
Cuny C, Kock R, Witte W. Livestock associated MRSA (LA-MRSA) and its relevance for humans in Germany. Int J Med Microbiol. 2013;303(6-7):331–7.
Monaco M, Pedroni P, Sanchini A, Bonomini A, Indelicato A, Pantosti A. Livestock-associated methicillin-resistant Staphylococcus aureus responsible for human colonization and infection in an area of Italy with high density of pig farming. BMC Infect Dis. 2013;13:258.
Pan A, Battisti A, Zoncada A, Bernieri F, Boldini M, Franco A, Giorgi M, Iurescia M, Lorenzotti S, Martinotti M, et al. Community-acquired methicillin-resistant Staphylococcus aureus ST398 infection, Italy. Emerg Infect Dis. 2009;15(5):845–7.
Soavi L, Stellini R, Signorini L, Antonini B, Pedroni P, Zanetti L, Milanesi B, Pantosti A, Matteelli A, Pan A, et al. Methicillin-resistant Staphylococcus aureus ST398, Italy. Emerg Infect Dis. 2010;16(2):346–8.
Mammina C, Cala C, Plano MR, Bonura C, Vella A, Monastero R, Palma DM. Ventilator-associated pneumonia and MRSA ST398, Italy. Emerg Infect Dis. 2010;16(4):730–1.
van Cleef BA, Monnet DL, Voss A, Krziwanek K, Allerberger F, Struelens M, Zemlickova H, Skov RL, Vuopio-Varkila J, Cuny C, et al. Livestock-associated methicillin-resistant Staphylococcus aureus in humans, Europe. Emerg Infect Dis. 2011;17(3):502–5.
Verkade E, Bergmans AM, Budding AE, van Belkum A, Savelkoul P, Buiting AG, Kluytmans J. Recent emergence of Staphylococcus aureus clonal complex 398 in human blood cultures. PLoS One. 2012;7(10):e41855.
Sung JM, Lloyd DH, Lindsay JA. Staphylococcus aureus host specificity: comparative genomics of human versus animal isolates by multi-strain microarray. Microbiology. 2008;154(Pt 7):1949–59.
Spoor LE, McAdam PR, Weinert LA, Rambaut A, Hasman H, Aarestrup FM, Kearns AM, Larsen AR, Skov RL, Fitzgerald JR. Livestock origin for a human pandemic clone of community-associated methicillin-resistant Staphylococcus aureus. MBio. 2013;4(4):e00356–13.
Cuny C, Wieler LH, Witte W. Livestock-associated MRSA: the impact on humans. Antibiotics (Basel). 2015;4(4):521–43.
Feltrin F, Alba P, Kraushaar B, Ianzano A, Argudin MA, Di Matteo P, Porrero MC, Aarestrup FM, Butaye P, Franco A, et al. A livestock-associated, multidrug-resistant, methicillin-resistant Staphylococcus aureus clonal complex 97 lineage spreading in dairy cattle and pigs in Italy. Appl Environ Microbiol. 2015;82(3):816–21.
Fitzgerald JR. Livestock-associated Staphylococcus aureus: origin, evolution and public health threat. Trends Microbiol. 2012;20(4):192–8.
Harrison EM, Coll F, Toleman MS, Blane B, Brown NM, Torok ME, Parkhill J, Peacock SJ. Genomic surveillance reveals low prevalence of livestock-associated methicillin-resistant Staphylococcus aureus in the East of England. Sci Rep. 2017;7(1):7406.
Winstel V, Liang C, Sanchez-Carballo P, Steglich M, Munar M, Broker BM, Penades JR, Nubel U, Holst O, Dandekar T, et al. Wall teichoic acid structure governs horizontal gene transfer between major bacterial pathogens. Nat Commun. 2013;4:2345.
Larsen J, Andersen PS, Winstel V, Peschel A. Staphylococcus aureus CC395 harbours a novel composite staphylococcal cassette chromosome mec element. J Antimicrob Chemother. 2017;72(4):1002–5.
Molina A, Del Campo R, Maiz L, Morosini MI, Lamas A, Baquero F, Canton R. High prevalence in cystic fibrosis patients of multiresistant hospital-acquired methicillin-resistant Staphylococcus aureus ST228-SCCmecI capable of biofilm formation. J Antimicrob Chemother. 2008;62(5):961–7.
Asghar AH. Molecular characterization of methicillin-resistant Staphylococcus aureus isolated from tertiary care hospitals. Pak J Med Sci. 2014;30(4):698–702.
McAdam PR, Holmes A, Templeton KE, Fitzgerald JR. Adaptive evolution of Staphylococcus aureus during chronic endobronchial infection of a cystic fibrosis patient. PLoS One. 2011;6(9):e24301.
Goerke C, Wolz C. Adaptation of Staphylococcus aureus to the cystic fibrosis lung. Int J Med Microbiol. 2010;300(8):520–5.
Kahl BC. Impact of Staphylococcus aureus on the pathogenesis of chronic cystic fibrosis lung disease. Int J Med Microbiol. 2010;300(8):514–9.
Cullen L, McClean S. Bacterial adaptation during chronic respiratory infections. Pathogens. 2015;4(1):66–89.
Rao Q, Shang W, Hu X, Rao X. Staphylococcus aureus ST121: a globally disseminated hypervirulent clone. J Med Microbiol. 2015;64(12):1462–73.
Goering RV, Shawar RM, Scangarella NE, O'Hara FP, Amrine-Madsen H, West JM, Dalessandro M, Becker JA, Walsh SL, Miller LA, et al. Molecular epidemiology of methicillin-resistant and methicillin-susceptible Staphylococcus aureus isolates from global clinical trials. J Clin Microbiol. 2008;46(9):2842–7.
Moore CL, Osaki-Kiyan P, Perri M, Donabedian S, Haque NZ, Chen A, Zervos MJ. USA600 (ST45) methicillin-resistant Staphylococcus aureus bloodstream infections in urban Detroit. J Clin Microbiol. 2010;48(6):2307–10.
Roberts JC. Classification of epidemic community-acquired methicillin-resistant Staphylococcus aureus by anatomical site of isolation. Biomed Res Int. 2014;2014:904283.
David MZ, Boyle-Vavra S, Zychowski DL, Daum RS. Methicillin-susceptible Staphylococcus aureus as a predominantly healthcare-associated pathogen: a possible reversal of roles? PLoS One. 2011;6(4):e18217.
Tinelli M, Monaco M, Vimercati M, Ceraminiello A, Pantosti A. Methicillin-susceptible Staphylococcus aureus in skin and soft tissue infections, Northern Italy. Emerg Infect Dis. 2009;15(2):250–7.
Baranovich T, Zaraket H, Shabana II, Nevzorova V, Turcutyuicov V, Suzuki H. Molecular characterization and susceptibility of methicillin-resistant and methicillin-susceptible Staphylococcus aureus isolates from hospitals and the community in Vladivostok, Russia. Clin Microbiol Infect. 2010;16(6):575–82.
Wu D, Wang Q, Yang Y, Geng W, Wang Q, Yu S, Yao K, Yuan L, Shen X. Epidemiology and molecular characteristics of community-associated methicillin-resistant and methicillin-susceptible Staphylococcus aureus from skin/soft tissue infections in a children's hospital in Beijing, China. Diagn Microbiol Infect Dis. 2010;67(1):1–8.
Rasigade JP, Laurent F, Lina G, Meugnier H, Bes M, Vandenesch F, Etienne J, Tristan A. Global distribution and evolution of Panton-Valentine leukocidin-positive methicillin-susceptible Staphylococcus aureus, 1981-2007. J Infect Dis. 2010;201(10):1589–97.
Isobe H, Takano T, Nishiyama A, Hung WC, Kuniyuki S, Shibuya Y, Reva I, Yabe S, Iwao Y, Higuchi W, et al. Evolution and virulence of Panton-Valentine leukocidin-positive ST30 methicillin-resistant Staphylococcus aureus in the past 30 years in Japan. Biomed Res. 2012;33(2):97–109.
Fernandez S, Ledo C, Lattar S, Noto Llana M, Bertelli AM, Di Gregorio S, Sordelli DO, Gomez MI, Mollerach ME. High virulence of methicillin resistant Staphylococcus aureus ST30-SCCmecIVc-spat019, the dominant community-associated clone in Argentina. Int J Med Microbiol. 2017;307(4-5):191–9.
Byrd AL, Belkaid Y, Segre JA. The human skin microbiome. Nat Rev Microbiol. 2018;16(3):143–55.
Tett A, Pasolli E, Farina S, Truong DT, Asnicar F, Zolfo M, Beghini F, Armanini F, Jousson O, De Sanctis V, et al. Unexplored diversity and strain-level structure of the skin microbiome associated with psoriasis. NPJ Biofilms Microbiomes. 2017;3:14.
Hiramatsu K, Cui L, Kuroda M, Ito T. The emergence and evolution of methicillin-resistant Staphylococcus aureus. Trends Microbiol. 2001;9(10):486–93.
Hanssen AM, Ericson Sollid JU. SCCmec in staphylococci: genes on the move. FEMS Immunol Med Microbiol. 2006;46(1):8–20.
Hartman B, Tomasz A. Altered penicillin-binding proteins in methicillin-resistant strains of Staphylococcus aureus. Antimicrob Agents Chemother. 1981;19(5):726–35.
Archer GL, Niemeyer DM, Thanassi JA, Pucci MJ. Dissemination among staphylococci of DNA sequences associated with methicillin resistance. Antimicrob Agents Chemother. 1994;38(3):447–54.
Campanile F, Bongiorno D, Falcone M, Vailati F, Pasticci MB, Perez M, Raglio A, Rumpianesi F, Scuderi C, Suter F, et al. Changing Italian nosocomial-community trends and heteroresistance in Staphylococcus aureus from bacteremia and endocarditis. Eur J Clin Microbiol Infect Dis. 2012;31(5):739–45.
Valsesia G, Rossi M, Bertschy S, Pfyffer GE. Emergence of SCCmec type IV and SCCmec type V methicillin-resistant Staphylococcus aureus containing the Panton-Valentine leukocidin genes in a large academic teaching hospital in Central Switzerland: external invaders or persisting circulators? J Clin Microbiol. 2010;48(3):720–7.
Tenover FC, Goering RV. Methicillin-resistant Staphylococcus aureus strain USA300: origin and epidemiology. J Antimicrob Chemother. 2009;64(3):441–6.
Pedersen LC, Benning MM, Holden HM. Structural investigation of the antibiotic and ATP-binding sites in kanamycin nucleotidyltransferase. Biochemistry. 1995;34(41):13305–11.
Gennimata D, Davies J, Tsiftsoglou AS. Bleomycin resistance in Staphylococcus aureus clinical isolates. J Antimicrob Chemother. 1996;37(1):65–75.
Dortet L, Girlich D, Virlouvet AL, Poirel L, Nordmann P, Iorga BI, Naas T. Characterization of BRPMBL, the bleomycin resistance protein associated with the Carbapenemase NDM. Antimicrob Agents Chemother. 2017;61(3):AAC–02413.
Burdeska A, Ott M, Bannwarth W, Then RL. Identical genes for trimethoprim-resistant dihydrofolate reductase from Staphylococcus aureus in Australia and Central Europe. FEBS Lett. 1990;266(1-2):159–62.
Rouch DA, Messerotti LJ, Loo LS, Jackson CA, Skurray RA. Trimethoprim resistance transposon Tn4003 from Staphylococcus aureus encodes genes for a dihydrofolate reductase and thymidylate synthetase flanked by three copies of IS257. Mol Microbiol. 1989;3(2):161–75.
Whitener CJ, Park SY, Browne FA, Parent LJ, Julian K, Bozdogan B, Appelbaum PC, Chaitram J, Weigel LM, Jernigan J, et al. Vancomycin-resistant Staphylococcus aureus in the absence of vancomycin exposure. Clin Infect Dis. 2004;38(8):1049–55.
Tenover FC, Weigel LM, Appelbaum PC, McDougal LK, Chaitram J, McAllister S, Clark N, Killgore G, O'Hara CM, Jevitt L, et al. Vancomycin-resistant Staphylococcus aureus isolate from a patient in Pennsylvania. Antimicrob Agents Chemother. 2004;48(1):275–80.
Stefani S, Campanile F, Santagati M, Mezzatesta ML, Cafiso V, Pacini G. Insights and clinical perspectives of daptomycin resistance in Staphylococcus aureus: a review of the available evidence. Int J Antimicrob Agents. 2015;46(3):278–89.
Noto MJ, Kreiswirth BN, Monk AB, Archer GL. Gene acquisition at the insertion site for SCCmec, the genomic island conferring methicillin resistance in Staphylococcus aureus. J Bacteriol. 2008;190(4):1276–83.
Wright GD. The antibiotic resistome: the nexus of chemical and genetic diversity. Nat Rev Microbiol. 2007;5(3):175–86.
Giersing BK, Dastgheyb SS, Modjarrad K, Moorthy V. Status of vaccine research and development of vaccines for Staphylococcus aureus. Vaccine. 2016;34(26):2962–6.
Verkaik NJ, van Wamel WJ, van Belkum A. Immunotherapeutic approaches against Staphylococcus aureus. Immunotherapy. 2011;3(9):1063–73.
Salgado-Pabon W, Schlievert PM. Models matter: the search for an effective Staphylococcus aureus vaccine. Nat Rev Microbiol. 2014;12(8):585–91.
Fattom AI, Horwith G, Fuller S, Propst M, Naso R. Development of StaphVAX, a polysaccharide conjugate vaccine against S. aureus infection: from the lab bench to phase III clinical trials. Vaccine. 2004;22(7):880–7.
Fattom A, Matalon A, Buerkert J, Taylor K, Damaso S, Boutriau D. Efficacy profile of a bivalent Staphylococcus aureus glycoconjugated vaccine in adults on hemodialysis: phase III randomized study. Hum Vaccin Immunother. 2015;11(3):632–41.
Fowler VG, Allen KB, Moreira ED, Moustafa M, Isgro F, Boucher HW, Corey GR, Carmeli Y, Betts R, Hartzel JS, et al. Effect of an investigational vaccine for preventing Staphylococcus aureus infections after cardiothoracic surgery: a randomized trial. JAMA. 2013;309(13):1368–78.
Fowler VG Jr, Proctor RA. Where does a Staphylococcus aureus vaccine stand? Clin Microbiol Infect. 2014;20(Suppl 5):66–75.
Golubchik T, Batty EM, Miller RR, Farr H, Young BC, Larner-Svensson H, Fung R, Godwin H, Knox K, Votintseva A, et al. Within-host evolution of Staphylococcus aureus during asymptomatic carriage. PLoS One. 2013;8(5):e61319.
Dreisbach A, van der Kooi-Pol MM, Otto A, Gronau K, Bonarius HP, Westra H, Groen H, Becher D, Hecker M, van Dijl JM. Surface shaving as a versatile tool to profile global interactions between human serum proteins and the Staphylococcus aureus cell surface. Proteomics. 2011;11(14):2921–30.
Hua L, Cohen TS, Shi Y, Datta V, Hilliard JJ, Tkaczyk C, Suzich J, Stover CK, Sellman BR. MEDI4893* promotes survival and extends the antibiotic treatment window in a Staphylococcus aureus immunocompromised pneumonia model. Antimicrob Agents Chemother. 2015;59(8):4526–32.
Bagnoli F. Staphylococcus aureus toxin antibodies: good companions of antibiotics and vaccines. Virulence. 2017;8(7):1037–42.
Yang L, Zhou H, Cheng P, Yang Y, Tong Y, Zuo Q, Feng Q, Zou Q, Zeng H. A novel bivalent fusion vaccine induces broad immunoprotection against Staphylococcus aureus infection in different murine models. Clin Immunol. 2018;188:85–93.
Frenck RW Jr, Creech CB, Sheldon EA, Seiden DJ, Kankam MK, Baber J, Zito E, Hubler R, Eiden J, Severs JM, et al. Safety, tolerability, and immunogenicity of a 4-antigen Staphylococcus aureus vaccine (SA4Ag): results from a first-in-human randomised, placebo-controlled phase 1/2 study. Vaccine. 2017;35(2):375–84.
Begier E, Seiden DJ, Patton M, Zito E, Severs J, Cooper D, Eiden J, Gruber WC, Jansen KU, Anderson AS, et al. SA4Ag, a 4-antigen Staphylococcus aureus vaccine, rapidly induces high levels of bacteria-killing antibodies. Vaccine. 2017;35(8):1132–9.
Anderson AS, Miller AA, Donald RG, Scully IL, Nanra JS, Cooper D, Jansen KU. Development of a multicomponent Staphylococcus aureus vaccine designed to counter multiple bacterial virulence factors. Hum Vaccin Immunother. 2012;8(11):1585–94.
Roetzer A, Jilma B, Eibl MM. Vaccine against toxic shock syndrome in a first-in-man clinical trial. Expert Rev Vaccines. 2017;16(2):81–3.
Narita K, Hu DL, Asano K, Nakane A. Vaccination with non-toxic mutant toxic shock syndrome toxin-1 induces IL-17-dependent protection against Staphylococcus aureus infection. Pathog Dis. 2015;73(4):ftv023.
Rouha H, Badarau A, Visram ZC, Battles MB, Prinz B, Magyarics Z, Nagy G, Mirkina I, Stulik L, Zerbs M, et al. Five birds, one stone: neutralization of alpha-hemolysin and 4 bi-component leukocidins of Staphylococcus aureus with a single human monoclonal antibody. MAbs. 2015;7(1):243–54.
Badarau A, Rouha H, Malafa S, Battles MB, Walker L, Nielson N, Dolezilkova I, Teubenbacher A, Banerjee S, Maierhofer B, et al. Context matters: the importance of dimerization-induced conformation of the LukGH leukocidin of Staphylococcus aureus for the generation of neutralizing antibodies. MAbs. 2016;8(7):1347–60.
Moustafa M, Aronoff GR, Chandran C, Hartzel JS, Smugar SS, Galphin CM, Mailloux LU, Brown E, Dinubile MJ, Kartsonis NA, et al. Phase IIa study of the immunogenicity and safety of the novel Staphylococcus aureus vaccine V710 in adults with end-stage renal disease receiving hemodialysis. Clin Vaccine Immunol. 2012;19(9):1509–16.
Schluepen C, Malito E, Marongiu A, Schirle M, McWhinnie E, Lo Surdo P, Biancucci M, Falugi F, Nardi-Dei V, Marchi S, et al. Mining the bacterial unknown proteome: identification and characterization of a novel family of highly conserved protective antigens in Staphylococcus aureus. Biochem J. 2013;455(3):273–84.
Bagnoli F, Fontana MR, Soldaini E, Mishra RP, Fiaschi L, Cartocci E, Nardi-Dei V, Ruggiero P, Nosari S, De Falco MG, et al. Vaccine composition formulated with a novel TLR7-dependent adjuvant induces high and broad protection against Staphylococcus aureus. Proc Natl Acad Sci U S A. 2015;112(12):3680–5.
Earls MR, Kinnevey PM, Brennan GI, Lazaris A, Skally M, O'Connell B, Humphreys H, Shore AC, Coleman DC. The recent emergence in hospitals of multidrug-resistant community-associated sequence type 1 and spa type t127 methicillin-resistant Staphylococcus aureus investigated by whole-genome sequencing: implications for screening. PLoS One. 2017;12(4):e0175542.
David MD, Kearns AM, Gossain S, Ganner M, Holmes A. Community-associated meticillin-resistant Staphylococcus aureus: nosocomial transmission in a neonatal unit. J Hosp Infect. 2006;64(3):244–50.
Creech CB, Frenck RW Jr, Sheldon EA, Seiden DJ, Kankam MK, Zito ET, Girgenti D, Severs JM, Immermann FW, McNeil LK, et al. Safety, tolerability, and immunogenicity of a single dose 4-antigen or 3-antigen Staphylococcus aureus vaccine in healthy older adults: results of a randomised trial. Vaccine. 2017;35(2):385–94.
Torre A, Bacconi M, Sammicheli C, Galletti B, Laera D, Fontana MR, Grandi G, De Gregorio E, Bagnoli F, Nuti S, et al. Four-component Staphylococcus aureus vaccine 4C-staph enhances Fcgamma receptor expression in neutrophils and monocytes and mitigates S. aureus infection in neutropenic mice. Infect Immun. 2015;83(8):3157–63.
Yu XQ, Robbie GJ, Wu Y, Esser MT, Jensen K, Schwartz HI, Bellamy T, Hernandez-Illas M, Jafri HS. Safety, tolerability, and pharmacokinetics of MEDI4893, an investigational, extended-half-life, anti-Staphylococcus aureus alpha-toxin human monoclonal antibody, in healthy adults. Antimicrob Agents Chemother. 2017;61(1):e01020–16.
Delfani S, Mohabati Mobarez A, Imani Fooladi AA, Amani J, Emaneini M. Protection of mice against Staphylococcus aureus infection by a recombinant protein ClfA-IsdB-Hlg as a vaccine candidate. Med Microbiol Immunol. 2016;205(1):47–55.
Landrum ML, Lalani T, Niknian M, Maguire JD, Hospenthal DR, Fattom A, Taylor K, Fraser J, Wilkins K, Ellis MW, et al. Safety and immunogenicity of a recombinant Staphylococcus aureus alpha-toxoid and a recombinant Panton-Valentine leukocidin subunit, in healthy adults. Hum Vaccin Immunother. 2017;13(4):791–801.
Schwameis M, Roppenser B, Firbas C, Gruener CS, Model N, Stich N, Roetzer A, Buchtele N, Jilma B, Eibl MM. Safety, tolerability, and immunogenicity of a recombinant toxic shock syndrome toxin (rTSST)-1 variant vaccine: a randomised, double-blind, adjuvant-controlled, dose escalation first-in-man trial. Lancet Infect Dis. 2016;16(9):1036–44.
The authors would like to acknowledge the laboratory of Computational Metagenomics for the valuable input. We would also like to thank Sander Wuyts and Stijn Wittouck for assistance in the plotting of the SCCmecIV graphs, and Dr. Megan Earls and Prof. David Coleman for sharing the ST1 SCCmecIV t127 genome.
This work was supported in part by the Italian Ministry of Health “Ricerca finalizzata” RF-2010-2316179 to GT. This work was also partially supported by the European Union H2020 Marie Curie grant (707345) to E.P. and by a European Union FP7 Marie Curie grant (PCIG13-618833), H2020-ERC-STG grant (MetaPG-716575), and by MIUR project FIR-RBFR13EWWI to NS.
Availability of data and materials
The datasets generated and analysed during the current study are available in the NCBI Sequence Read Archive (BioProject accession number PRJNA400143, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA400143/).
Ethics approval and consent to participate
The study received ethical approval from the Paediatric Ethics Committee, Autonomous Section of the Regional Ethics Committee for Clinical Trials at the Children’s Meyer Hospital, Florence on July 1, 2014. Data collection has been performed in accordance with the Declaration of Helsinki. All patients gave written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics of the single isolates, including collection details, genome assembly statistics, genomic features, and results of antibiotic susceptibility testing. (XLSX 69 kb)
Figure S1. Pangenome analysis statistics. Figure S2. Phylogenetic model based on gene presence/absence. (PDF 580 kb)
Table S2. Presence/absence profile of virulence and antibiotic resistance genes in the cohort. (XLSX 6063 kb)
Table S3. Characteristics of the genomes included in the ST1 analyses. (XLSX 11 kb)
Table S4. Most relevant clones represented in the cohort and their abundance. (XLSX 9 kb)
Table S5. Bayesian clock models tested and their results. (XLSX 10 kb)