Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article

Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe

David M. Aanensen, Edward J. Feil, Matthew T. G. Holden, Janina Dordel, Corin A. Yeats, Artemij Fedosejev, Richard Goater, Santiago Castillo-Ramírez, Jukka Corander, Caroline Colijn, Monika A. Chlebowicz, Leo Schouls, Max Heck, Gerlinde Pluister, Raymond Ruimy, Gunnar Kahlmeter, Jenny Åhman, Erika Matuschek, Alexander W. Friedrich, Julian Parkhill, Stephen D. Bentley, Brian G. Spratt, Hajo Grundmann
Keith P. Klugman, Editor
David M. Aanensen
aDepartment of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
bThe Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edward J. Feil
cThe Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Bath, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthew T. G. Holden
dSchool of Medicine, University of St. Andrews, St. Andrews, United Kingdom
ePathogen Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Janina Dordel
ePathogen Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
fDepartment of Biology, Drexel University, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Corin A. Yeats
aDepartment of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
bThe Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Artemij Fedosejev
aDepartment of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
bThe Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard Goater
bThe Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Santiago Castillo-Ramírez
gPrograma de Genómica Evolutiva, Centro de Ciencias Genómicas, Universidad Nacional Autónoma de Mexico, Cuernavaca, Morelos, Mexico
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jukka Corander
hHelsinki Institute for Information Technology HIIT, Aalto, Finland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Caroline Colijn
iDepartment of Mathematics, Imperial College London, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Monika A. Chlebowicz
jDepartment of Medical Microbiology, University Medical Center Groningen, Rijksuniversteit Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leo Schouls
kNational Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Max Heck
kNational Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gerlinde Pluister
kNational Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Raymond Ruimy
lCentre Hospitalier Universitaire de Nice, Nice, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gunnar Kahlmeter
mEUCAST Development Laboratory, Växjö, Sweden
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jenny Åhman
mEUCAST Development Laboratory, Växjö, Sweden
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Erika Matuschek
mEUCAST Development Laboratory, Växjö, Sweden
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander W. Friedrich
jDepartment of Medical Microbiology, University Medical Center Groningen, Rijksuniversteit Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julian Parkhill
ePathogen Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Julian Parkhill
Stephen D. Bentley
ePathogen Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brian G. Spratt
aDepartment of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hajo Grundmann
jDepartment of Medical Microbiology, University Medical Center Groningen, Rijksuniversteit Groningen, Groningen, The Netherlands
nDepartment of Infection Prevention and Hospital Hygiene, Faculty of Medicine, University of Freiburg, Freiburg, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keith P. Klugman
Department of Global Health, Emory University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jennifer Gardy
B.C. Centre for Disease Control
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geoffrey Coombs
Murdoch University, Australia
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.00444-16
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The implementation of routine whole-genome sequencing (WGS) promises to transform our ability to monitor the emergence and spread of bacterial pathogens. Here we combined WGS data from 308 invasive Staphylococcus aureus isolates corresponding to a pan-European population snapshot, with epidemiological and resistance data. Geospatial visualization of the data is made possible by a generic software tool designed for public health purposes that is available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc). Our analysis demonstrates that high-risk clones can be identified on the basis of population level properties such as clonal relatedness, abundance, and spatial structuring and by inferring virulence and resistance properties on the basis of gene content. We also show that in silico predictions of antibiotic resistance profiles are at least as reliable as phenotypic testing. We argue that this work provides a comprehensive road map illustrating the three vital components for future molecular epidemiological surveillance: (i) large-scale structured surveys, (ii) WGS, and (iii) community-oriented database infrastructure and analysis tools.

IMPORTANCE The spread of antibiotic-resistant bacteria is a public health emergency of global concern, threatening medical intervention at every level of health care delivery. Several recent studies have demonstrated the promise of routine whole-genome sequencing (WGS) of bacterial pathogens for epidemiological surveillance, outbreak detection, and infection control. However, as this technology becomes more widely adopted, the key challenges of generating representative national and international data sets and the development of bioinformatic tools to manage and interpret the data become increasingly pertinent. This study provides a road map for the integration of WGS data into routine pathogen surveillance. We emphasize the importance of large-scale routine surveys to provide the population context for more targeted or localized investigation and the development of open-access bioinformatic tools to provide the means to combine and compare independently generated data with publicly available data sets.

INTRODUCTION

Most bacterial infections of humans are caused by organisms that have a large population size and short generation times, and lineages with novel properties emerge and expand within observable time scales (1). Lineages that present a serious threat to public health are designated high-risk clones (HRCs) (2); these often combine enhanced virulence or transmission potential with multiple-antibiotic resistance. As HRCs are often difficult to treat and are associated with significant morbidity, mortality, and economic cost, they require targeted surveillance and containment at the population level. The considerable management challenge for public health microbiologists in tackling HRCs can be broken down into three tasks: (i) identification of public health risks posed by emerging and/or expanding HRCs, (ii) assessment of this risks by predicting important clinical and epidemiological consequences, and (iii) risk management through the implementation of prevention and control strategies. With the cost of sequencing entire bacterial genomes in steady decline and the development of powerful bioinformatic tools gathering pace, whole-genome sequencing (WGS) will inevitably be widely implemented for routine epidemiological surveillance. The dissemination of genomic data through established national and international networks of collaborating specialist laboratories will lead to increased awareness and shorter response times for HRCs. Ideally, these networks will be based on a shared bioinformatic infrastructure that links molecular data (genome sequences) and metadata (time, place, clinical details, and additional variables) with tools that help appraise the clinical and public health relevance of any given entry. This poses significant technical, ethical, and political challenges; however, the benefits of this fundamental shift in molecular epidemiology cannot be overstated. The efficient management and interrogation of genome data by public health and medical audiences may ultimately lead to the erosion of conventional reference diagnostic tasks such as identification to the species level, characterization of clinically important virulence or resistance phenotypes, and identification of outbreaks and inference of national or international transmission.

Staphylococcus aureus represents an epidemiological paradigm because of the undisputed public health relevance of this species, which is characterized by multiple-antibiotic resistance and a potential for swift dissemination through health care, social, and farm animal production networks. Here we demonstrate the utility of WGS when applied to a continental-scale representative “population snapshot” by using a novel data visualization platform. We sequenced 308 S. aureus isolates responsible for invasive infections that were recovered from 186 hospitals in 21 countries across Europe in a 6-month period. We consider three analytical strands: (i) a representative phylogeographic analysis that defines HRCs on a population level, (ii) an analysis of the dynamics of virulence and resistance carried by mobile genetic elements (MGEs), and (iii) an in silico ascertainment of antibiotic resistance encompassing 19 antibiotic compounds of clinical relevance.

To underline the added value for public health decision making when WGS data are supplemented with epidemiological metadata collected through structured surveys, we developed a web application (Microreact) that allows easy access and visualization of our data by medical and public health audiences and is available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc). This tool can be used for any appropriate data set where a phylogenetic tree and associated metadata are available (see http://www.microreact.org).

RESULTS

We chose a random sample (n = 308; 10.6%) of isolates collected as part of a European structured survey of S. aureus from invasive diseases (3). Sixty percent (n = 185) of these isolates were methicillin-sensitive S. aureus (MSSA), and 40% (n = 123) were methicillin-resistant S. aureus (MRSA). A total of 235,226 SNP sites within the core genome were identified by mapping against a single reference genome, HO 5096 0412 (sequence type 22 [ST22]). We have divided the analysis and interpretation of the results into three parts: (i) a broad phylogenetic analysis and more fine-scaled consideration of exemplar lineages of high public health relevance, (ii) the distribution and dynamics of the accessory genes conferring virulence and resistance traits in the context of individual HRCs, and (iii) a comparison of antibiotic susceptibility profiles ascertained by in silico prediction from genome data with conventional susceptibility testing carried out in Staphylococcal Reference Laboratories (SRLs) and the European Committee on Antimicrobial Susceptibility Testing (EUCAST) Development Laboratory (EDL).

Broad and fine-scale phylogenetic analyses for the identification of high-risk clones.Analysis of core SNPs resolved the population into 6 major (CC5, CC22, CC8, CC30, CC45, CC15) and 10 minor clonal complexes (CCs) and STs (CC1, ST20, ST25, ST7, CC121, ST88, CC12, CC398, ST101, ST72), as previously defined by multilocus sequence typing (MLST)/eBURST (4) (Table 1; Fig. 1 and 2). The six major CCs are each represented by at least 24 isolates, whereas the minor lineages are represented by a minimum of four isolates and a maximum of 14. Subdivisions within a given CC and within individual STs are resolved by the data (Fig. 2E), which identify very closely related isolates resulting from recent expansion. We refer to these groups within CCs simply as clusters. These include many well-known MRSA “clones” previously defined by MLST, pulsed-field gel electrophoresis, and spa typing.

View this table:
  • View inline
  • View popup
TABLE 1 

Abundance, diversity, and proportion of MRSA isolates in each major or minor CC detected in the sample

FIG 1 
  • Open in new tab
  • Download powerpoint
FIG 1 

Phylogenetic relationship of the invasive S. aureus population circulating in Europe in 2006. A rooted neighbor-joining tree based on 235,226 genomewide core SNPs is shown. Lineages are highlighted and named according to the corresponding CC or ST.

FIG 2 
  • Open in new tab
  • Download powerpoint
FIG 2 

Phylogeny decoration. Colors of branches indicate MSSA (green) and MRSA (red) states. Each isolate is annotated by affiliation with a CC or ST (A), SSCmec type and MSSA (green) or MRSA (red) state (B), or antibiotic resistance profile (C). Red boxes indicate the presence of genetic resistance markers, black dots indicate phenotypic resistance, and gray boxes highlight resistance in reference genomes. (D) Size and composition of the accessory genome based on the number of noncore homologous groups with further categorization according to MGE type. (E) Close-up of phylogenetic trees of the six major lineages.

The distribution of MRSA and MSSA isolates is not random with respect to CCs (Table 1). Ninety-three percent of the MRSA isolates belonged to only three CCs, CC22, CC5, and CC8, although these CCs represent only approximately half of the isolates (Table 1). Over 70% of the isolates in CC5, CC8, and CC22 are MRSA, whereas all of the isolates in CC15 are MSSA. CC45 shows moderate levels of MRSA (23%), and CC30 shows low levels (6%). The overall topology of the tree is highly consistent with previous studies (5) (Fig. 1 and 2). The impact of recombination on tree topology was investigated by computing the frequency of homoplasy (phylogenetic conflicts within the data). This confirmed previous work that revealed a high rate of homoplasy between the CCs (37.8% of the core SNPs) but very low levels when considering the variation within each of the CCs (0.62% of the core SNPs) (5). We are therefore confident that the fine-scaled phylogeographic inferences within individual CCs have not been seriously compromised by recombination. Intra-CC diversity was calculated simply as the mean number of SNP differences in all pairwise comparisons, and the average time since each pair of genomes shared a common ancestor was calculated on the basis of the published S. aureus mutation rate (Table 1) (6–9). This revealed similar levels of diversity between CCs and that, on average, the common ancestor of any pair of genomes belonging to the same major CC existed from the mid-1930s (CC45) to the mid-1970s (CC15).

A key observation from our previous work on spa typing was that MRSA variants tend to show geographic clustering, whereas MSSA variants do not (3). The WGS data permit more detailed comparisons of spatial signals within and between MRSA and MSSA clusters belonging to individual CCs. We examined three major CCs that illustrate the usefulness of WGS for public health and infection control; other CCs (CC45, CC8, ST239, and CC15) are discussed in the supplemental material (see Fig. S1 to S4, respectively). All of the phylogenetic analyses, presence/absence of accessory genes, and associated metadata are available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc) by use of the Microreact tool; a detailed guide is provided in Text S1 the supplemental material, with CC45 as an example.

Text S1 

Appendix: Walkthrough of web application for genomic pathogen surveillance. Download Text S1, DOCX file, 1.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure S1 

Phylogenetic reconstruction of CC45. Branch color indicates MSSA (green) or MRSA (red). The symbols at branch tips indicate the geographic origins of these isolates. SCCmec IV subtypes are shown. Download Figure S1, PDF file, 0.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure S2 

Phylogenetic reconstruction of CC8 (excluding ST239). Branch color indicates MSSA (green) or MRSA (red). The symbols at branch tips indicate the geographic origins of these isolates. Close-ups of the two SCCmec IV clusters show their starlike radiation. SCCmec IV subtypes are shown. The third cluster represents USA300 isolates consisting of two reference isolates and one Belgian isolate. Download Figure S2, PDF file, 0.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure S3 

(A) Phylogenetic reconstruction of ST239 isolates in the sample. Red branches depict all of the MRSA isolates. The symbols at the tips of the branches indicate the geographic origins of the isolates. Isolates carrying the sasX gene are indicated by the letter X. Blue shading highlights isolates belonging to the Asian clade according to Harris et al. (S. R. Harris, et al., Science 327:469-474, 2010). (B) Phylogenetic tree of combined ST239 data obtained from isolates from the reports of Harris et al. and Castillo-Ramírez et al. (S. Castillo-Ramírez, et al., Genome Biol 13:R126, 2012) and isolates in panel A. Clusters are named in accordance with the reports of Harris et al. and Castillo-Ramírez et al. Download Figure S3, PDF file, 0.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure S4 

Phylogenetic reconstruction of CC15. Starlike phylogeny indicates the recent rapid expansion of this MSSA-only lineage. Symbols at the tips indicate the geographic origins of these isolates. The branch leading to an ST582 isolate—a single-locus variant of CC15—is shortened, and the number of SNPs defining that branch is shown. Download Figure S4, PDF file, 0.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

CC5.CC5 is the most abundant CC in this study, being represented by 70 isolates, 80% of which are MRSA (red branches) (Fig. 3). Three MRSA clusters (ST225, ST228, and ST125) exhibit strong geographic structuring indicative of regional expansion. One MSSA (green branches) cluster is resolved on the tree but is not geographically restricted. This cluster largely corresponds to spa type t002, and this observation is consistent with our previous work showing that MSSA genotypes are widely distributed across Europe (3). The 13 isolates within the ST225 cluster all contain SCCmec type II elements. Seven of these originate from Germany, six are from the Czech Republic, and the partitioning within the phylogenetic tree coincides perfectly with the country of origin, although some of these isolates originate from areas close to the shared border (e.g., Ústí nad Labem in the Czech Republic and Grossenhain in Germany are just 122 km apart by road). The data thus point to the border between Germany and the Czech Republic as being a barrier to health care referral practices and thus MRSA spread.

FIG 3 
  • Open in new tab
  • Download powerpoint
FIG 3 

Phylogenetic reconstruction of CC5. Branch color indicates MSSA (green) or MRSA (red). Clusters described in Results are shaded gray. Symbols at the tips indicate the geographic origins of these isolates. SCCmec IV subtypes are shown for ST125.

The MRSA cluster corresponding to ST228 and related STs consists of 15 isolates all harboring a type I SCCmec element. Nine of these isolates define a very tight subcluster corresponding to ST111 (n = 8) and ST1481 (n = 1). Eight of these nine isolates originate from Croatia, which is indicative of rapid epidemic radiation in that country, and three isolates from Split cluster together on the tree, reflecting geographic structuring on a national level. One ST111 isolate originates from Sweden, consistent with the probable importation of ST111 into Sweden from Croatia by travel. The MRSA ST125 cluster consists of eight MRSA isolates restricted to central and northern Spain. For further investigation, see the Microreact tool, where supplementary analysis and geographic detail can be explored.

CC22.CC22 contains epidemic MRSA-15 (EMRSA-15) (Fig. 4A), which is currently the most abundant and fastest growing health care-associated MRSA (HA-MRSA) clone in Europe (10). This clone accounts for the majority of the CC22 isolates (31/41; 76%) in our sample, which are characterized by an SCCmec type IVh element and form a tightly clustered starlike phylogeny (Fig. 4A). Fourteen of these isolates originated from the United Kingdom, 11 were from Germany, and 7 were from Portugal. The topology of the tree supports a United Kingdom origin of EMRSA-15, followed by separate introductions from the United Kingdom into Germany and Portugal, consistent with the work of Holden et al. (8). In the German cluster, isolates from Berlin show the fewest SNP differences from the United Kingdom isolates, implicating that city as the initial point of introduction from the United Kingdom (Fig. 4B). This clone subsequently spread to the surrounding cities of Kiel, Frankfurt (Oder), Hanover, and Magdeburg, supporting epidemiological observations made by the Robert Koch Institute during investigations in the 1990s (11). Similarly, the Portuguese CC22 cluster points toward Lisbon as the point of introduction, with subsequent spread to Braga (in the far north) and Coimbra (central). All of the CC22 isolates located basal to the EMRSA-15 subtree are MSSA. They form a more diverse population, genetically and geographically, and represent the MSSA reservoir from which EMRSA-15 emerged (8).

FIG 4 
  • Open in new tab
  • Download powerpoint
FIG 4 

Phylogenetic reconstruction of CC22. Branch color indicates MSSA (green) or MRSA (red). The EMRSA-15 cluster is shaded gray. Symbols at the tips of the branches indicate the geographic origins of these isolates. A cluster consisting of isolates from Berlin indicating the possible point of EMRSA-15 introduction into Germany from the United Kingdom is shaded a darker gray. The position of an isolate from Lisbon is shown indicating the possible location of its entry into Portugal.

CC30.CC30 contains successful HA-MRSA and community-associated MRSA (CA-MRSA) epidemic clones, including EMRSA-16 and the Panton-Valentine leukocidin (PVL) toxin-positive Southwest Pacific (SWP) clone. Figure 5A shows the tree for the 32 CC30 isolates present in the current sample and illustrates that only two (6.5%) were MRSA. One of the MRSA isolates is closely related to the MRSA252 reference, which is a HA-EMRSA-16 isolate (ST36), and the other (ST1829) is related to the TCH60 reference, a representative isolate of the SWP clone (14). The majority of 28 MSSA isolates form a striking starlike phylogeny with little geographic structure (for detailed investigation of geographic origin, see the Microreact tool). To place this lineage in a broader context, we combined our data with those of McAdam et al. (7), who recently sequenced representatives of the three major CC30 epidemic lineages; EMRSA-16, SWP, and the historical MSSA phage type 80/81 clone (Fig. 5B). This revealed that the large MSSA radiation corresponds to a successful progenitor population from which HA-EMRSA-16 emerged. This MSSA cluster, referred to by MacAdam et al. as “other epidemic,” is responsible for a significant disease burden in the community and in hospitals and encompasses a variety of spa types. Given the public health significance of this cluster, we tentatively suggest the designation EMSSA-ST30. All of the isolates in this cluster contain the SNPs in hla and agrC previously reported to attenuate virulence (18) (Fig. 5B). The presence of these SNPs among community-associated invasive disease isolates further challenges the suggestion that they play a role in nosocomial adaptation (18).

FIG 5 
  • Open in new tab
  • Download powerpoint
FIG 5 

(A) Phylogenetic reconstruction of CC30 isolates in the sample. Branch color indicates MSSA (green) or MRSA (red). Reference genomes are named and clusters are highlighted and named in accordance with the report of McAdam et al. (7). The dashed line indicates the long branch leading to three ST34 isolates that evolved through the acquisition of a 200-kb homologous replacement within the chromosome from an ST10/ST145-like parent (66). (B) Phylogenetic tree of combined CC30 data obtained from isolates from the study of McAdam et al. (7) and isolates from panel A. Colors and cluster names are as in panel A. Light gray shading indicates isolates carrying SNPs in the hla and agrC genes thought to restrict these lineages to health care settings. (C) Representation of successive clonal radiations within the recent evolutionary history of CC30. These radiations correspond to recognized HRCs, both contemporary and historic. The SWP clone is a historic diversified starlike expansion with relatively long branches. Phage type 80/81 probably emerged from within this starlike expansion, as did the current MSSA HRC, which we have termed EMSSA-ST30. Finally, EMRSA-16 emerged from within the successful EMSSA-ST30 lineage, resulting in a more recent, and tightly clustered, starlike expansion.

Distribution and dynamics of accessory genes.The accessory genome consists of genes that are variably present in a defined population and can be major drivers of adaptation (15). An understanding of the dynamics of key MGEs such as phages and the SCCmec elements, which carry determinants of the virulence and resistance phenotypes, is critical for monitoring the emergence and diversification of successful HRCs. We explored the diversity and distribution of accessory genes in our de novo assembled genomes by assigning genes to noncore homology groups (ncHGs). We categorized these groups according to the types of mobile elements with which they are associated, specifically, phages, SCCmec elements, plasmids, S. aureus pathogenicity islands, and transposons (Fig. 2). Fifty-seven percent of the ncHGs in the sample were phage associated, representing the most dominant category by far. The total number of ncHGs varies between different CCs; CC8 harbors the highest average number of ncHGs per isolate at 197, and CC15 has the lowest at 71 (see Fig. S5 in the supplemental material).

Figure S5 

Number of ncHGs per major lineage. The bar chart shows the average number of ncHGs per isolate split into MGE type for each of the six major lineages. Each pie chart indicates the proportion of MGE types for each major lineage. Download Figure S5, PDF file, 0.2 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Figure 2D illustrates the variation in the accessory genome within single CCs, and Fig. 6 shows the total number of ncHGs shared by each pair of genomes. Many more ncHGs are shared by isolates belonging to the same CC than by isolates belonging to different CCs. This supports the concept of a “core variable” genome (16), representing genes that are universally present at the scale of an individual CC but are variably present or absent with respect to the species as a whole. We also noted cases of MRSA clusters within CCs that have acquired a distinctive repertoire of ncHGs, evident as dark squares in Fig. 6; for example, the ST239 cluster (CC8), the ST228 cluster (CC5), and the ST225 cluster (also CC5). The high level of consistency between phylogenetic relatedness and the ncHG repertoire is expected over these very fine phylogenetic scales because of common inheritance. However, this consistency rapidly decays with increasing core genome divergence, particularly for highly dynamic elements such as phages. It is evident from Fig. 6 that CC5 and CC8 have a high number of ncHGs in common, resulting in a single large square that encompasses both CCs. The large proportion of MRSA strains within CC5 and CC8 does not, by itself explain, this observation, as it is still apparent when SCCmec is excluded (data not shown).

FIG 6 
  • Open in new tab
  • Download powerpoint
FIG 6 

Conservation of ncHGs across the invasive S. aureus population in Europe. Isolates are arranged in the tree order of Fig. 2 along the left and top. The ncHGs in the pairwise comparison are displayed as a heat map matrix, and colors represent the total numbers of ncHGs that genome pairs have in common; darker shading indicates more ncHGs in common (see the scale at the bottom). The dark squares corresponding to (from the top) CC45, CC30, CC22, CC8, CC5, and (at the bottom) CC15.

Comparing the total number of ncHGs within MRSA and MSSA isolates in the same CC revealed that MRSA clusters contain more ncHGs than closely related MSSA genomes, even when the SCCmec elements were excluded (see Fig. S6 in the supplemental material). For example, considering CC5 and CC22, MRSA isolates contain approximately 50% more ncHGs than MSSA isolates of the same complex. Although it is unclear whether the acquisition of SCCmec is the cause or the consequence of an increased propensity to integrate horizontally acquired DNA, this observation may be relevant to the finding that MRSA strains tend to contain many more resistance determinants than MSSA strains, conferring combined resistance to multiple classes of antibiotics.

Figure S6 

MRSA or MSSA genome size. The bar chart shows the average number of ncHGs per MSSA or MRSA isolate, respectively, excluding the HGs associated with SCCmec, for each major lineage. Asterisks indicate lineages with a significant size difference of the accessory genome in MSSA and MRSA isolates. CC8, P = 0.039; CC5, P = 0.0325 (statistically significant); CC22, P = 0.0539 (not quite statistically significant). Download Figure S6, PDF file, 0.03 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Phages are important from a public health perspective, as they are carriers of virulence genes and resistance determinants. We assayed the distribution of seven phage types on the basis of the presence/absence of their integrase genes as described previously (17). The overall number of prophages per genome ranged from zero (9%) to five (0.3%), with a median of two (46%). Sa3int was the most commonly observed phage type, being present in 82% of our isolates, followed by Sa2int (32%), and Sa1int (27%). Sa4int was very rare, being observed in only four isolates (1.3%). The distribution of these phage types and associated cargo genes across the tree reveals important differences between the CCs and between different clusters belonging to the same CC (see Fig. S7 in the supplemental material). For example, the Sa3int prophage was found in all of the CCs, with the exception of CC15. This prophage can harbor four genes belonging to the immune evasion cluster, chp, sak, scn, and sea/sep, and in our data mostly harbored sak and scn, with the addition of either chp (e.g., CC22, CC45) or sea/sep (e.g., CC8, except USA300). Notably, chp and scn are common in CC15 even though the prophage itself is absent. The most likely explanation for this pattern is that the prophage was once integrated and then lost from the genome, leaving chp and scn behind.

Figure S7 

Prophage distribution. Rooted neighbor-joining tree like that in Fig. 1. Colors of branches indicate MSSA (green) and MRSA (red). Each isolate is annotated for affiliation to CCs or STs (A), SSCmec type and MSSA (green) or MRSA (red) (B), and seven prophage types classified on the basis of the presence or absence of their integrase genes (C). Colored boxes indicate presence. Black boxes indicate the presence of virulence genes associated with the prophage. (D) Colored boxes indicate the presence of two virulence genes (sasX and tstH) and acr—the gene commonly used to detect ACMEs. (E) Size and composition of the accessory genome based on the number of ncHGs with further categorization according to MGE type. (F) Close-up of phylogenetic trees of the six major lineages. Download Figure S7, PDF file, 0.8 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

SasX is a cell wall-anchored surface protein that is linked to the epidemiological success of ST239 in China and Southeast Asia (19). The ϕSPβ-like prophage that encodes SasX was identified among three ST239 isolates from Poland (see Fig. S7D). Two of these isolates cluster with the TW20 reference, which is representative of the “Asian clade” of ST239 (see Fig. S3), pointing to the introduction of these isolates into Poland from Asia. The third sasX-positive ST239 isolate in our data is more likely to be a native Polish variant and most likely acquired sasX from the imported Asian strains cocirculating in Poland. The data thus point to the horizontal dissemination, presumably via lysogenic conversion, of the phage-borne sasX virulence factor within Europe. Other examples of the horizontal dissemination of virulence genes are evident in our data. We note only a single representative of USA300, an MRSA strain from Belgium that closely clusters with the reference USA300 genome and contains a type I arginine catabolic MGE (ACME) (see Fig. S2). However, we also note a single CC5 SCCmec type V MRSA isolate from Portugal that has also acquired a type I ACME region (Fig. 3; see Fig. S7D). The distribution of other toxin and virulence genes, such as the presence of tstH in the cluster we designated EMSSA-ST30, is illustrated in Fig. S7D and incorporated into the visualization platform at the project URL.

Multiple SCCmec types (I, II, III, IV, V, and IV) were present among the 120 MRSA isolates. Type IV elements are the most widely distributed among different lineages, being present in all five of the major CCs that contain MRSA strains. From the phylogenetic analyses, we estimate that a minimum of 20 independent acquisitions of this element have occurred in our sample. However, the rate of SCCmec acquisition clearly varies between CCs, whereas only a single acquisition of the SCCmec type IVh element is observed in CC22, in CC5, there have been a minimum of seven different acquisition events encompassing five different SCCmec types. For example, the Spanish ST125 cluster harbors both SCCmec type IVc and IVs elements, indicating multiple acquisitions. In contrast, no CC15 isolates have acquired SCCmec. The reversion of MRSA to MSSA because of the deletion of SCCmec appears to be a relatively rare event and was identified in only two CC8 isolates in our sample (see Fig. S2).

WGS as a tool for predicting antibiotic resistance.The EDL tested all of the isolates against 16 antibiotics. In addition, all MRSA isolates were tested against a further three antibiotics that are prescribed mainly for the treatment of infections caused by MRSA. A total of 5,288 in silico predictions of resistance/susceptibility were made and compared against the EUCAST reference results in a blinded fashion. Of these, 1,075 were predicted to be resistant, compared with 1,050 identified phenotypically as resistant by the EDL, with concordance noted for 5,213 (98.6%) of the individual tests (Fig. 2; Table 2). However, there were some noteworthy exceptions. For the aminoglycoside amikacin, 21 isolates were falsely predicted to be resistant on the basis of the presence of the aphA-3 kanamycin resistance gene and aadD. It has previously been noted that the presence of these two aminoglycoside-modifying enzyme genes is a poor indicator of amikacin resistance (20, 21). Resistance to mupirocin was not predicted in five isolates that had inhibition zones (29 mm) close to the susceptibility breakpoint (≥30 mm). The reason for this discrepancy is likely to be normal variation or the choice of medium. There were also eight incorrect predictions of erythromycin susceptibility. In these cases, the result could readily be explained by loss of the ermC-carrying plasmid after sequencing but prior to phenotypic testing. Similarly, the loss of a BlaZ β-lactamase-carrying plasmid could explain the three false-negative predictions of penicillin resistance.

View this table:
  • View inline
  • View popup
TABLE 2 

Comparison of antibiotic resistances predicted by in silico and SRL test results against the EDL reference

We also determined the concordance of the antimicrobial susceptibility test (AST) results provided by the European Staphylococcus Reference Laboratories (SRLs), which contributed the isolates to our sample, with the EDL reference data. The SRLs submitted a total of 2,252 AST results. Of these, 2,203 (97.8%) were concordant with the EUCAST test results (Fig. 2; Table 2). Thus, it is clear that a small amount of discordance between phenotypic tests can be expected, even when they are carried out by experienced reference laboratories. We note that the degree of discordance between the in silico predictions and the EUCAST data is comparable to the degree of discordance between the two sets of phenotypic data. This demonstrates that in silico predictions were at least as reliable, in terms of matching the gold standard, as the AST results generated by independent reference laboratories.

DISCUSSION

A consistent and comprehensive sampling frame is a crucial component of effective pathogen surveillance in public health, though in reality, this may involve significant logistical and political challenges. By taking a random sample of a larger collection assembled as part of a pan-European structured survey, we have minimized the problems associated with phylogenetic discovery bias (22) and are confident that our data capture a meaningful snapshot of S. aureus invasive disease isolates in Europe at the time of sampling. However, as with all sampling frameworks, important caveats remain. First, all of the isolates were recovered from invasive infections, primarily blood, meaning that we have underrepresented isolates causing skin and soft tissue infections. Second, laboratories submitted the first five MSSA isolates and the first five MRSA isolates from individual patients. A small number of laboratories did not receive five MRSA isolates during the 6-month survey because of a very low frequency of MRSA disease in their regions, for example in Scandinavia, and these laboratories completed their quota with MSSA isolates.

Despite these caveats, the combined analysis of WGS with epidemiological metadata addresses the following three key elements of managing infectious disease threats in public health.

Genetic population structure and identification of HRCs.The WGS data demonstrate that the disease-causing population of S. aureus is readily partitioned into highly discrete subpopulations or CCs. These CCs vary in their potential to spawn HRCs such as EMRSA-15 (8) or USA300 (9, 22), and our analysis suggests that these differences may be related to the varying propensity of each CC to acquire exogenous DNA (see Fig. S5 in the supplemental material). Analysis of fine-scale genetic variation and spatial structuring of the well-characterized HRCs circulating in Europe during the sampling period helped identify key population level properties, or HRC signatures, that can be used to recognize candidate HRCs, even in the absence of detailed phenotypic data (Fig. 3 to 5; see Fig. S1 to S4). Central to this is a consideration of phylogenetic tree topology, relative abundance, and geographic structuring.

Perhaps the most striking example is EMRSA-15, which emerged from the CC22 population approximately 30 years ago (8) and has subsequently become pandemic. Phylogenetic analysis of ST22 reveals a comet-shaped phylogeny, with two distinct parts; the comet head consists of the starlike radiation of EMRSA-15, reflecting recent rapid expansion, while the tail represents the more diverse MSSA progenitor population (Fig. 4A).

We find multiple similar signatures of rapid clonal expansions that coincide with recognized successful clones such as ST225, ST228, and ST125 within CC5 (Fig. 3). We also note previously unrecognized clusters in Europe that constitute candidate HRCs in the MSSA population, the most notable example being the CC30 MSSA designated EMSSA-ST30 here. By combining our data with the WGS data from a previous study of well-known CC30 HRCs (Fig. 5B), we placed EMSSA-ST30 within a broader evolutionary context. This revealed that the EMRSA-16 HRC emerged from EMSSA-ST30, as hypothesized originally from MLST data (24). This highlights the importance of recognizing successful MSSA lineages not only as HRCs in their own right but also in their role as likely progenitors of EMRSA (25).

Assessing the risks posed by virulence and resistance determinants for public health.CCs differ in their propensities to acquire and maintain accessory genome elements (see Fig. S1 in the supplemental material), a factor that may influence the emergence of HRCs. CC30 represents a good example of a lineage that has acquired elements carrying important virulence and resistance determinants. In addition to the phage-borne PVL toxin gene harbored in phage type 80/81 and the SWP clone, the cluster of isolates that we designated EMSSA-ST30 maintains the tstH gene encoding a superantigen that can cause toxic shock syndrome at high frequency (26) (see Fig. S7). This cluster is therefore reminiscent of, and likely descended from, the ST30 variants responsible for the tampon-associated toxic shock syndrome that emerged in the 1980s (27). This provides another line of evidence for the public health risk posed by this widespread clone.

The success of ST239 in China and much of Southeast Asia has been largely attributed to SasX, a surface-anchored protein that modulates host interactions and transmissibility and is carried on the ϕSPβ-like prophage (19). In Poland, we found evidence that ST239 isolates from the successful Asian lineage spread repeatedly to Europe (12) but also that this sasX virulence determinant has been transduced from the European ST239 genetic background (see Fig. S3). This observation justifies heightened epidemiological vigilance with respect to this HRC, not only for particular epidemiologically successful clones but also for virulence determinants they can spread.

A key characteristic of HRCs has been the emergence of antibiotic resistance. WGS provides the opportunity to scan genomes for all known genetic determinants of antibiotic resistance. We examined the distribution of known resistance determinants and the extent to which it is possible to predict resistance profiles from sequence data in a blinded fashion. Whereas previous studies have shown a high concordance between phenotypic data and genotypic predictions (8, 13, 28), we sought to extend these studies by including a greater number of antibiotics (19 in total) and by including a representative sample of isolates that cause invasive disease among patients in European hospitals and that we consider the clinically most relevant representatives of the S. aureus population. Further, as discrepancies between the two methods could arise from uncertainties in the phenotypic data, as well as inaccurate predictions from genotyping, we compared two independent sets of phenotypic data, one contributed by various SRLs and the second generated at the EDL. Taking the EDL data to represent the “gold standard,” our genotypic predictions show higher rates of concordance with these data (EDL versus genotype, 98.6%) than do the data generated by the SRLs (EDL versus SRL, 97.8%, Table 2). We conclude that genotypic prediction is at least as reliable as routine phenotypic testing and that any discrepancies between the two approaches are just as likely to represent inaccuracies in the phenotypic testing as inaccurate genotypic predictions.

We do not, however, propose that genotypic prediction should replace phenotypic testing, as such a strategy would be vulnerable to the emergence of new and uncharacterized antibiotic resistance mechanisms. Moreover, in silico prediction for other organisms has proved more challenging, especially for Gram-negative bacteria (29), where the present understanding of the genetic basis of resistance is less comprehensive. Despite these limitations, WGS data clearly provide additional objective evidence that is not prone to heterogeneities inherent to conventional phenotypic test methods deployed in different reference laboratories. Moreover, in the longer term, WGS will deliver a digitized and cumulative record that will address the need for internationally agreed standards for collection of data and reporting on antibacterial resistance in human health and for harmonizing standards across medical, veterinary, and agricultural sectors as required by the WHO, the World Organization for Animal Health, and the Food and Agriculture Organization of the United Nations (30).

Deployment of informed and targeted prevention and control strategies.We have argued for a two-pronged approach to the identification of emerging HRCs and an assessment of the threats they pose to public health, (i) a consideration of clonal relatedness, abundance, and phylogeographic structure at the population level and (ii) mining of the accessory genome repertoire to ascertain likely virulence and resistance properties. These two perspectives should also be combined when considering appropriate containment measures. The recognition of an HRC through large-scale surveillance would enable the development of tailored rapid diagnostic tests based on distinguishing SNPs or accessory gene signatures. This will allow for selective screening and targeted containment by decontamination and isolation strategies that are easier to implement, more economical, and more likely to be effective than more generic procedures. The expedition of specifically tailored tests during outbreaks would allow health authorities or infection control practitioners to screen potential hosts for colonization or infection, thereby reducing the chance of onward transmission. WGS allows for rapid appraisal of outbreak sources and transmission pathways that will also help in weighing up infection control priorities. Moreover, information regarding the virulence gene repertoire associated with an emerging or outbreak HRC will inform clinicians and medical microbiologists of likely clinical manifestations.

Ultimately, the key to managing HRCs lies in making the data available in an open and intuitive format for infection control and public health audiences. This democratization of the data increases the collective power of interpretation (i.e., the identification of HRCs) while decreasing the necessity for local expertise in bioinformatics. The coupling of large-scale population sampling by WGS within open-access and freely available web resources empowers the community to identify clone-specific signatures (canonical SNPs and/or accessory genes), promoting the design of HRC-specific rapid and cost-effective molecular diagnostic tests (31).

With this communication, we aimed to demonstrate how the integration of WGS into epidemiological surveillance programs provides the means for both the early warning of emerging HRCs and a robust assessment of associated public health threats. However, the advent of the underlying sequencing technology addresses only a small part of the challenge of managing the threat from infectious disease. In order to exploit this technical advance to its maximum potential, two things must happen. First, national reference laboratories need to agree to, and abide by, common standards for the contribution of isolates to structured surveys, ideally at 5-year or shorter intervals, on pancontinental and, ideally, global scales. Second, platforms for bioinformatics need to be developed to allow intuitive management and exploration of the data. These data sets should conform to internationally curated standards for sets of genes and mutations that are recognized as key virulence or resistance determinants, and we propose the initiation of internationally curated data sets to act as a gold standard resource for genomic antimicrobial resistance determinants.

The establishment of cumulative databases will engender a far richer understanding of the detailed dynamics underpinning clonal emergence and replacement, national and international transmission, and the horizontal transfer of core genes and MGEs. The increasing public health threat from bacterial pathogens is, in large part, down to the ability of these organisms to rapidly adapt through the dissemination of genes and mobile elements. Our best chance of managing these threats in the future is to emulate, as far as possible, this resource and data sharing through the development of international surveillance networks and a common data exchange infrastructure.

MATERIALS AND METHODS

Sampling framework and bacterial isolates.From September 2006 to February 2007, 357 laboratories serving 450 hospitals in 26 countries collected nearly 3,000 MSSA and MRSA isolates from patients with invasive S. aureus infections, as described previously (3). Approximately 10% (n = 308) of these were randomly selected for sequencing. Isolate details are available for download at the project URL.

Genomic library creation and sequencing.For each sample, index-tagged libraries were prepared and sequenced in Illumina Genome Analyzer GAII cells with 54-base paired-end reads or Illumina HiSeq with 75-bp paired-end reads. Downstream processing utilized the index tag sequence information to assign reads to individual samples.

Reference genomes.To provide a wider context to the data, we utilized 26 fully annotated complete S. aureus reference genomes. The strain names and accession numbers are as follows: 04-02981, CP001844 (6); COL, CP000046 (32); ECT-R2, FR714927 (33); ED133, CP001996 (34); ED98, CP001781 (35); HO 5096 0412, HE681097 (8); JH1, CP000736 (36); JH9, CP000703 (36); JDK6008, CP002120 (37); JKD6159, CP002114 (38); LGA251, FR821779 (39); MRSA252, BX571856; MSSA476, BX571857 (40); Mu3, AP009324 (41); Mu50, BA000017 (42); MW2, BA000033 (43); N315, BA000018 (42); NCTC8325, CP000253 (44); Newman, AP009351 (45); RF122, AJ938182 (46); BB155, LN854556; ST398, AM990992 (47); TCH60, CP002110; TW20, FN433596 (48); USA300_FPR3757, CP000255 (23); USA300_TCH1516, CP000730 (49).

Detection of SNPs in the core genome.The paired-end reads for the survey isolates and the 26 publicly available genomes were mapped with SMALT (http://www.sanger.ac.uk/resources/software/smalt/) against the chromosome of S. aureus strain HO 5096 0412 (EMRSA-15, ST22; accession number HE681097) (8). The reference genomes were mapped by using simulated paired fastq data (54-bp paired ends with a 300-bp insert). SNPs were identified as previously described (50). Indels were identified with the Genome Analysis Toolkit (https://www.broadinstitute.org/gatk/) (51). Unmapped reads and sequences that were not present in all of the genomes were not considered part of the core genome, and SNPs from these regions were not included in the phylogenetic analysis. SNPs falling within MGE regions were also excluded from the phylogenetic analysis, as described previously (8).

Phylogenetic analysis.Phylogenetic analysis of the isolates combined was carried out by the neighbor-joining method as implemented in FastTree (52). In cases where data were combined with already published data, RAxML v8.0 (53) was used for tree reconstruction based on maximum likelihood. Only SNPs corresponding to the core genomes were used for phylogenetic analysis. In silico reconstruction of STs conventionally generated by MLST for each isolate was carried out with the sequence data as described by Croucher et al. (50). To root the tree, we used the sequence of BB155, an S. aureus reference genome that belongs to ST152. This was considered an appropriate outgroup as, according to MLST data, ST152 is a divergent lineage approximately equidistant from each of the major S. aureus STs and CCs previously recorded in Europe (54).

Pangenome analysis.Assemblies were produced with a pipeline comprising three steps: (i) correction of sequencing errors in reads with Quake v0.3 (55), (ii) de novo assembly of corrected reads into scaffolds with SOAPdenovo2 (56) by using a kmer of 31, and (iii) remapping of the corrected reads to fill scaffold gaps with GapCloser v1.12 (56). Prediction of protein coding sequences (CDSs) in assemblies was done with Prodigal (57), and FASTA sequences and coordinate positions were extracted. An all-versus-all BlastP (58) search of translated sequences was used as the input for the definition of homology groups (HGs) with TribeMCL (59) (cutoff, 1e-50) with group membership stored for querying. From a total of 749,089 putative CDSs, 4,281 HGs were defined with various numbers of members (1 to 1,827). Core HGs were defined as being present in at least a single copy in each assembly, and all others were considered ncHGs. SCCmec elements were typed by mapping reads to a pseudomolecule of the known SCCmec variants and also by comparative analysis using the assemblies.

Detection of prophages and virulence genes.Prophages were classified by using conserved areas of the integrase (int) genes as described previously (17) by an in silico PCR approach. Virulence and toxin genes of public health interest were similarly detected through in silico searches for previously published primer sequences (eta [60]; luk-PV [61]; chp, sak, scn, sea, sep, hlb [62]; lukM [63]; sasX [18]; tstH [64]; acr [65]).

Genomic prediction of antibiotic resistance.Resistance profiles for 19 antibiotics were predicted in silico from the sequence data as previously described (8). In brief, the literature was comprehensively mined for the known genetic mechanisms of antibiotic resistance in S. aureus (see Table S1 in the supplemental material). Antibiotic resistance conferred by SNPs in components of the core chromosome were identified from the mapping data. Antibiotic resistance conferred by accessory genes was identified by comparing de novo assemblies against a database of S. aureus resistance genes with BlastN (58) and by mapping sequence reads to a pseudomolecule consisting of concatenated antibiotic resistance genes as detailed by Holden et al. (8).

Table S1 

Known genetic mechanisms of antibiotic resistance in S. aureus used to predict genotypic resistance. Table S1, XLSX file, 0.04 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Phenotypic testing of antibiotic resistance.The antibiotic susceptibility of all S. aureus isolates was tested by the standardized EUCAST disk diffusion method in the EDL, Växjö, Sweden, without prior knowledge of the sequence data or metadata. Full methodological details are available in the EUCAST Disk Diffusion Test Manual, v 3.0, 2013 (http://www.eucast.org). Phenotypic resistance was defined by applying inhibition zone diameter epidemiological cutoff values (ECOFFs) and EUCAST clinical breakpoints (EUCAST breakpoint tables for interpretation of MICs and zone diameters, version 3.1, 2013 [http://www.eucast.org]). We pooled the intermediate and resistant categories into a single nonsusceptible category. For our sample, this partition was fully consistent with the classification into wild-type and non-wild-type isolates defined by the EUCAST ECOFFs. The antibiotics tested were penicillin, cefoxitin, ciprofloxacin, moxifloxacin, amikacin, gentamicin, tobramycin, erythromycin, clindamycin, tetracycline, tigecycline, fusidic acid, linezolid, mupirocin, rifampin, and trimethoprim. MRSA resistance to clinically relevant reserve antibiotics (teicoplanin, vancomycin, and daptomycin) not suitable for disk diffusion testing was determined by using the respective MIC ECOFFs (2.0, 2.0, and 1.0 mg/liter, respectively) with E tests (bioMérieux Clinical Diagnostics, Marcy l’Etoile, France). Phenotypic test results were dichotomized and grouped into wild-type and non-wild-type categories. Results of genomic predictions were compared with the phenotypic classifications as resistant or susceptible as defined by the ECOFFS. Moreover, results that were made available by the SRLs were equivalently dichotomized and also compared to those generated by the EDL in a consistent manner.

Web application.Microreact is a Node.js application written in JavaScript. Locational data are displayed by using the Google Maps Application Programming Interface, and phylogenetic data (in Newick tree format) are displayed by using a custom tree viewer developed for the HTML5 Canvas element (PhyloCanvas). Image files (.png) of annotated trees and subtrees can be saved. Instructions and examples are shown at the microreact.org website. The project URL for exploring the data described in this article is http://www.microreact.org/project/EkUvg9uY?tt=rc.

Nucleotide accession numbers.WGS data for all of the isolates tested in this study have been deposited in the Sequence Read Archive, and accession numbers are included in the metadata available at the project URL.

Table S2 

Affiliations of and contact information for members of the European SRL Working Group. Table S2, XLSX file, 0.02 MB.
Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

D.M.A., C.Y., and A.F. were supported by Wellcome Trust grant 099202 and by The Centre for Genomic Pathogen Surveillance (R.G.). J.D., S.D.B., J.P., and M.T.G.H. were supported by Wellcome Trust grant 098051. B.G.S. was supported by Wellcome Trust grant 089472. S.C.-R. was funded by UKCRC grant MRC G1000803.

We thank the diagnostic laboratories across Europe that contributed their isolates during the structured surveys and whose efforts are crucial to the identification of HRCs and understanding of their spread through European populations. We also thank the sequencing and informatics teams at the Sanger Institute for their assistance and the Wellcome Trust for their support of the Sanger Institute Pathogen Genomics and Biology groups. We thank Carola Schinkel of Tomorrow’s Events for organizing the annual training and planning of capacity building workshops and Corinna Glasner for lab assistance.

D.M.A., E.J.F., M.T.G.H., and H.G. designed the experiments/the study. Members of the SRL Working Group collected isolates and epidemiological data. Members of the SRL Working Group, L.S., M.H., G.P., and H.G. managed isolate and data collection. G.K., J.Å., and E.M. carried out reference susceptibility testing. D.M.A., E.J.F., M.T.G.H., J.D., M.A.C., S.C.-R., J.C., C.C., R.R., and H.G. analyzed the data. D.M.A., E.J.F., M.T.G.H., J.D., and H.G. wrote the paper. D.M.A., C.A.Y., A.F., and R.G. developed the public-domain web-based interactive mapping tool. J.P., S.D.B., A.F., and B.G.S. provided scientific advice and coedited the manuscript.

The members of the European SRL Working Group are Helmut Mittermayer, Karina Krziwanek, Sabine Stumvoll, Walter Koller, Olivier Denis, Marc Struelens, Dimitr Nashev, Ana Budimir, Smilja Kalenic, Despo Pieridou-Bagatzouni, Vladislav Jakubu, Helena Zemlickova, Henrik Westh, Anders Rhod Larsen, Robert Skov, Frederic Laurent, Jerome Ettienne, Birgit Strommenger, Wolfgang Witte, Sofia Vourli, Alkis Vatopoulos, Anni Vainio, Jaana Vuopio-Varkila, Miklos Fuzi, Erika Ungvári, Stephan Murchan, Angela Rossney, Edvins Miklasevics, Arta Balode, Gunnsteinn Haraldsson, Karl G. Kristinsson, Monica Monaco, Annalisa Pantosti, Michael Borg, Marga van Santen-Verheuvel, Xander Huijsdens, Lillian Marstein, Trond Jacobsen, Gunnar Skov Simonsen, Marta Aires-de-Sousa, Herminia de Lencastre, Agnieszka Luczak-Kadlubowska, Waleria Hryniewicz, Monica Straut, Irina Codita, Maria Perez-Vazquez, Jesus Oteo Iglesias, Vesna Cvitkovic Spik, Manica Mueller-Premru, Sara Haeggman, Barbro Olsson-Liljequist, Matthew Ellington, and Angela Kearns. For their affiliations, contact information, and contributions, see Table S2 in the supplemental material.

FOOTNOTES

    • Received 5 April 2016
    • Accepted 12 April 2016
    • Published 5 May 2016
  • Copyright © 2016 Aanensen et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Smith JM,
    2. Feil EJ,
    3. Smith NH
    . 2000. Population structure and evolutionary dynamics of pathogenic bacteria. Bioessays 22:1115–1122. doi:10.1002/1521-1878(200012)22:12<1115::AID-BIES9>3.3.CO;2-I.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Cantón R,
    2. Ruiz-Garbajosa P
    . 2011. Co-resistance: an opportunity for the bacteria and resistance genes. Curr Opin Pharmacol 11:477–485. doi:10.1016/j.coph.2011.07.007.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Grundmann H,
    2. Aanensen DM,
    3. van den Wijngaard CC,
    4. Spratt BG,
    5. Harmsen D,
    6. Friedrich AW, European Staphylococcal Reference Laboratory Working Group
    . 2010. Geographic distribution of Staphylococcus aureus causing invasive infections in Europe: a molecular-epidemiological analysis. PLoS Med 7:e1000215. doi:10.1371/journal.pmed.1000215.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Feil EJ,
    2. Li BC,
    3. Aanensen DM,
    4. Hanage WP,
    5. Spratt BG
    . 2004. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol 186:1518–1530. doi:10.1128/JB.186.5.1518-1530.2004.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Everitt RG,
    2. Didelot X,
    3. Batty EM,
    4. Miller RR,
    5. Knox K,
    6. Young BC,
    7. Bowden R,
    8. Auton A,
    9. Votintseva A,
    10. Larner-Svensson H,
    11. Charlesworth J,
    12. Golubchik T,
    13. Ip CL,
    14. Godwin H,
    15. Fung R,
    16. Peto TE,
    17. Walker AS,
    18. Crook DW,
    19. Wilson DJ
    . 2014. Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus. Nat Commun 5:3956. doi:10.1038/ncomms4956.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Nübel U,
    2. Dordel J,
    3. Kurt K,
    4. Strommenger B,
    5. Westh H,
    6. Shukla SK,
    7. Žemličková H,
    8. Leblois R,
    9. Wirth T,
    10. Jombart T,
    11. Balloux F,
    12. Witte W
    . 2010. A timescale for evolution, population expansion, and spatial spread of an emerging clone of methicillin-resistant Staphylococcus aureus. PLoS Pathog 6:e1000855. doi:10.1371/journal.ppat.1000855.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. McAdam PR,
    2. Templeton KE,
    3. Edwards GF,
    4. Holden MTG,
    5. Feil EJ,
    6. Aanensen DM,
    7. Bargawi HJA,
    8. Spratt BG,
    9. Bentley SD,
    10. Parkhill J,
    11. Enright MC,
    12. Holmes A,
    13. Girvan EK,
    14. Godfrey PA,
    15. Feldgarden M,
    16. Kearns AM,
    17. Rambaut A,
    18. Robinson DA,
    19. Fitzgerald JR
    . 2012. Molecular tracing of the emergence, adaptation, and transmission of hospital-associated methicillin-resistant Staphylococcus aureus. Proc Natl Acad Sci USA 109:9107–9112. doi:10.1073/pnas.1202869109.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Holden MTG,
    2. Hsu LY,
    3. Kurt K,
    4. Weinert LA,
    5. Mather AE,
    6. Harris SR,
    7. Strommenger B,
    8. Layer F,
    9. Witte W,
    10. de Lencastre H,
    11. Skov R,
    12. Westh H,
    13. Zemlickova H,
    14. Coombs G,
    15. Kearns AM,
    16. Hill RLR,
    17. Edgeworth J,
    18. Gould I,
    19. Gant V,
    20. Cooke J,
    21. Edwards GF,
    22. McAdam PR,
    23. Templeton KE,
    24. McCann A,
    25. Zhou Z,
    26. Castillo-Ramírez S,
    27. Feil EJ,
    28. Hudson LO,
    29. Enright MC,
    30. Balloux F,
    31. Aanensen DM,
    32. Spratt BG,
    33. Fitzgerald JR,
    34. Parkhill J,
    35. Achtman M,
    36. Bentley SD,
    37. Nubel U
    . 2013. A genomic portrait of the emergence, evolution, and global spread of a methicillin-resistant Staphylococcus aureus pandemic. Genome Res 23:653–664. doi:10.1101/gr.147710.112.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Uhlemann AC,
    2. Dordel J,
    3. Knox JR,
    4. Raven KE,
    5. Parkhill J,
    6. Holden MTG,
    7. Peacock SJ,
    8. Lowy FD
    . 2014. Molecular tracing of the emergence, diversification, and transmission of S. aureus sequence type 8 in a New York community. Proc Natl Acad Sci USA 111:6738–6743. doi:10.1073/pnas.1401006111.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Grundmann H,
    2. Schouls L,
    3. Aanensen D,
    4. Pluister G,
    5. Tami A,
    6. Chlebowicz M,
    7. Glasner C,
    8. Sabat A,
    9. Weist K,
    10. Heuer O,
    11. Friedrich A, ESCMID Study Group on Molecular Epidemiological Markers, European Staphylococcal Reference Laboratory Working Group
    . 2014. The dynamic changes of dominant clones of Staphylococcus aureus causing bloodstream infections in the European region: results of a second structured survey. Euro Surveill 19:20987. doi:10.2807/1560-7917.ES2014.19.49.20987.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Witte W,
    2. Enright M,
    3. Schmitz FJ,
    4. Cuny C,
    5. Braulke C,
    6. Heuck D
    . 2001. Characteristics of a new epidemic MRSA in Germany ancestral to United Kingdom EMRSA 15. Int J Med Microbiol 290:677–682. doi:10.1016/S1438-4221(01)80006-0.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Harris SR,
    2. Feil EJ,
    3. Holden MTG,
    4. Quail MA,
    5. Nickerson EK,
    6. Chantratita N,
    7. Gardete S,
    8. Tavares A,
    9. Day N,
    10. Lindsay JA,
    11. Edgeworth JD,
    12. de Lencastre H,
    13. Parkhill J,
    14. Peacock SJ,
    15. Bentley SD
    . 2010. Evolution of MRSA during hospital transmission and intercontinental spread. Science 327:469–474. doi:10.1126/science.1182395.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Köser CU,
    2. Holden MTG,
    3. Ellington MJ,
    4. Cartwright EJP,
    5. Brown NM,
    6. Ogilvy-Stuart AL,
    7. Hsu LY,
    8. Chewapreecha C,
    9. Croucher NJ,
    10. Harris SR,
    11. Sanders M,
    12. Enright MC,
    13. Dougan G,
    14. Bentley SD,
    15. Parkhill J,
    16. Fraser LJ,
    17. Betley JR,
    18. Schulz-Trieglaff OB,
    19. Smith GP,
    20. Peacock SJ
    . 2012. Rapid whole-genome sequencing for investigation of a neonatal MRSA outbreak. N Engl J Med 366:2267–2275. doi:10.1056/NEJMoa1109910.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    1. Collignon P,
    2. Gosbell I,
    3. Vickery A,
    4. Nimmo G,
    5. Stylianopoulos T,
    6. Gottlieb T
    . 1998. Community-acquired meticillin-resistant Staphylococcus aureus in Australia. Australian Group on Antimicrobial Resistance. Lancet 352:145–146. doi:10.1016/S0140-6736(98)85051-4.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Lindsay JA,
    2. Holden MTG
    . 2004. Staphylococcus aureus: superbug, super genome? Trends Microbiol 12:378–385. doi:10.1016/j.tim.2004.06.004.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Lindsay JA,
    2. Moore CE,
    3. Day NP,
    4. Peacock SJ,
    5. Witney AA,
    6. Stabler RA,
    7. Husain SE,
    8. Butcher PD,
    9. Hinds J
    . 2006. Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J Bacteriol 188:669–676. doi:10.1128/JB.188.2.669-676.2006.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Goerke C,
    2. Pantucek R,
    3. Holtfreter S,
    4. Schulte B,
    5. Zink M,
    6. Grumann D,
    7. Broker BM,
    8. Doskar J,
    9. Wolz C
    . 2009. Diversity of prophages in dominant Staphylococcus aureus clonal lineages. J Bacteriol 191:3462–3468. doi:10.1128/JB.01804-08.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. DeLeo FR,
    2. Kennedy AD,
    3. Chen L,
    4. Wardenburg JB,
    5. Kobayashi SD,
    6. Mathema B,
    7. Braughton KR,
    8. Whitney AR,
    9. Villaruz AE,
    10. Martens CA,
    11. Porcella SF,
    12. McGavin MJ,
    13. Otto M,
    14. Musser JM,
    15. Kreiswirth BN
    . 2011. Molecular differentiation of historic phage-type 80/81 and contemporary epidemic Staphylococcus aureus. Proc Natl Acad Sci U S A 108:18091–18096. doi:10.1073/pnas.1111084108.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Li M,
    2. Du X,
    3. Villaruz AE,
    4. Diep BA,
    5. Wang D,
    6. Song Y,
    7. Tian Y,
    8. Hu J,
    9. Yu F,
    10. Lu Y,
    11. Otto M
    . 2012. MRSA epidemic linked to a quickly spreading colonization and virulence determinant. Nat Med 18:816–819. doi:10.1038/nm.2692.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Davies J,
    2. Courvalin P
    . 1977. Mechanisms of resistance to aminoglycosides. Am J Med 62:868–872. doi:10.1016/0002-9343(77)90654-4.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Courvalin P,
    2. Davies J
    . 1977. Plasmid-medicated [sic] aminoglycoside phosphotransferase of broad substrate range that phosphorylates amikacin. Antimicrob Agents Chemother 11:619–624. doi:10.1128/AAC.11.4.619.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Pearson T,
    2. Busch JD,
    3. Ravel J,
    4. Read TD,
    5. Rhoton SD,
    6. U’Ren JM,
    7. Simonson TS,
    8. Kachur SM,
    9. Leadem RR,
    10. Cardon ML,
    11. Van Ert MN,
    12. Huynh LY,
    13. Fraser CM,
    14. Keim P
    . 2004. Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing. Proc Natl Acad Sci U S A 101:13536–13541. doi:10.1073/pnas.0403844101.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Diep BA,
    2. Gill SR,
    3. Chang RF,
    4. Phan TH,
    5. Chen JH,
    6. Davidson MG,
    7. Lin F,
    8. Lin J,
    9. Carleton HA,
    10. Mongodin EF,
    11. Sensabaugh GF,
    12. Perdreau-Remington F
    . 2006. Complete genome sequence of USA300, an epidemic clone of community-acquired meticillin-resistant Staphylococcus aureus. Lancet 367:731–739. doi:10.1016/S0140-6736(06)68231-7.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Enright MC,
    2. Day NP,
    3. Davies CE,
    4. Peacock SJ,
    5. Spratt BG
    . 2000. Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J Clin Microbiol 38:1008–1015.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Enright MC,
    2. Robinson DA,
    3. Randle G,
    4. Feil EJ,
    5. Grundmann H,
    6. Spratt BG
    . 2002. The evolutionary history of methicillin-resistant Staphylococcus aureus (MRSA). Proc Natl Acad Sci U S A 99:7687–7692. doi:10.1073/pnas.122108599.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Stach CS,
    2. Herrera A,
    3. Schlievert PM
    . 2014. Staphylococcal superantigens interact with multiple host receptors to cause serious diseases. Immunol Res 59:177–181. doi:10.1007/s12026-014-8539-7.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Schlievert PM,
    2. Shands KN,
    3. Dan BB,
    4. Schmid GP,
    5. Nishimura RD
    . 1981. Identification and characterization of an exotoxin from Staphylococcus aureus associated with toxic-shock syndrome. J Infect Dis 143:509–516. doi:10.1093/infdis/143.4.509.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    1. Gordon NC,
    2. Price JR,
    3. Cole K,
    4. Everitt R,
    5. Morgan M,
    6. Finney J,
    7. Kearns AM,
    8. Pichon B,
    9. Young B,
    10. Wilson DJ,
    11. Llewelyn MJ,
    12. Paul J,
    13. Peto TEA,
    14. Crook DW,
    15. Walker AS,
    16. Golubchik T
    . 2014. Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing. J Clin Microbiol 52:1182–1191. doi:10.1128/JCM.03117-13.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Reuter S,
    2. Ellington MJ,
    3. Cartwright EJP,
    4. Köser CU,
    5. Török ME,
    6. Gouliouris T,
    7. Harris SR,
    8. Brown NM,
    9. Holden MTG,
    10. Quail M,
    11. Parkhill J,
    12. Smith GP,
    13. Bentley SD,
    14. Peacock SJ
    . 2013. Rapid bacterial whole-genome sequencing to enhance diagnostic and public health microbiology. JAMA Intern Med 173:1397–1404. doi:10.1001/jamainternmed.2013.7734.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. World Health Organization
    . 2014. Antimicrobial resistance. Draft global action plan on antimicrobial resistance. Report by the secretariat. World Health Organization, Geneva, Switzerland. Accessed 24 March 2016. http://apps.who.int/gb/ebwha/pdf_files/EB136/B136_20-en.pdf. Accessed 24 March 2016.
  31. 31.↵
    1. Sangal V,
    2. Holt KE,
    3. Yuan J,
    4. Brown DJ,
    5. Filliol-Toutain I,
    6. Weill F,
    7. Kim D,
    8. da Silveira WD,
    9. Pickard D,
    10. Thomson NR,
    11. Parkhill J,
    12. Yu J
    . 2013. Global phylogeny of Shigella sonnei strains from limited single nucleotide polymorphisms (SNPs) and development of a rapid and cost-effective SNP-typing scheme for strain identification by high-resolution melting analysis. J Clin Microbiol 51:303–305. doi:10.1128/JCM.02238-12.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Gill SR,
    2. Fouts DE,
    3. Archer GL,
    4. Mongodin EF,
    5. Deboy RT,
    6. Ravel J,
    7. Paulsen IT,
    8. Kolonay JF,
    9. Brinkac L,
    10. Beanan M,
    11. Dodson RJ,
    12. Daugherty SC,
    13. Madupu R,
    14. Angiuoli SV,
    15. Durkin AS,
    16. Haft DH,
    17. Vamathevan J,
    18. Khouri H,
    19. Utterback T,
    20. Lee C,
    21. Dimitrov G,
    22. Jiang L,
    23. Qin H,
    24. Weidman J,
    25. Tran K,
    26. Kang K,
    27. Hance IR,
    28. Nelson KE,
    29. Fraser CM
    . 2005. Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain. J Bacteriol 187:2426–2438. doi:10.1128/JB.187.7.2426-2438.2005.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Lindqvist M,
    2. Isaksson B,
    3. Grub C,
    4. Jonassen TØ,
    5. Hällgren A
    . 2012. Detection and characterisation of SCCmec remnants in multiresistant methicillin-susceptible Staphylococcus aureus causing a clonal outbreak in a Swedish county. Eur J Clin Microbiol Infect Dis 31:141–147. doi:10.1007/s10096-011-1286-y.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Guinane CM,
    2. Ben Zakour NL,
    3. Tormo-Mas MA,
    4. Weinert LA,
    5. Lowder BV,
    6. Cartwright RA,
    7. Smyth DS,
    8. Smyth CJ,
    9. Lindsay JA,
    10. Gould KA,
    11. Witney A,
    12. Hinds J,
    13. Bollback JP,
    14. Rambaut A,
    15. Penades JR,
    16. Fitzgerald JR
    . 2010. Evolutionary genomics of Staphylococcus aureus reveals insights into the origin and molecular basis of ruminant host adaptation. Genome Biol Evol 2:454–466. doi:10.1093/gbe/evq031.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Lowder BV,
    2. Guinane CM,
    3. Ben Zakour NL,
    4. Weinert LA,
    5. Conway-Morris A,
    6. Cartwright RA,
    7. Simpson AJ,
    8. Rambaut A,
    9. Nubel U,
    10. Fitzgerald JR
    . 2009. Recent human-to-poultry host jump, adaptation, and pandemic spread of Staphylococcus aureus. Proc Natl Acad Sci USA 106:19545–19550. doi:10.1073/pnas.0909285106.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Mwangi MM,
    2. Wu SW,
    3. Zhou Y,
    4. Sieradzki K,
    5. de Lencastre H,
    6. Richardson P,
    7. Bruce D,
    8. Rubin E,
    9. Myers E,
    10. Siggia ED,
    11. Tomasz A
    . 2007. Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing. Proc Natl Acad Sci USA 104:9451–9456. doi:10.1073/pnas.0609839104.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Howden BP,
    2. Seemann T,
    3. Harrison PF,
    4. McEvoy CR,
    5. Stanton J-L,
    6. Rand CJ,
    7. Mason CW,
    8. Jensen SO,
    9. Firth N,
    10. Davies JK,
    11. Johnson PDR,
    12. Stinear TP
    . 2010. Complete genome sequence of Staphylococcus aureus strain JKD6008, an ST239 clone of methicillin-resistant Staphylococcus aureus with intermediate-level vancomycin resistance. J Bacteriol 192:5848–5849. doi:10.1128/JB.00951-10.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Chua K,
    2. Seemann T,
    3. Harrison PF,
    4. Davies JK,
    5. Coutts SJ,
    6. Chen H,
    7. Haring V,
    8. Moore R,
    9. Howden BP,
    10. Stinear TP
    . 2010. Complete genome sequence of Staphylococcus aureus strain JKD6159, a unique Australian clone of ST93-IV community methicillin-resistant Staphylococcus aureus. J Bacteriol 192:5556–5557. doi:10.1128/JB.00878-10.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. García-Álvarez L,
    2. Holden MT,
    3. Lindsay H,
    4. Webb CR,
    5. Brown DF,
    6. Curran MD,
    7. Walpole E,
    8. Brooks K,
    9. Pickard DJ,
    10. Teale C,
    11. Parkhill J,
    12. Bentley SD,
    13. Edwards GF,
    14. Girvan EK,
    15. Kearns AM,
    16. Pichon B,
    17. Hill RL,
    18. Larsen AR,
    19. Skov RL,
    20. Peacock SJ,
    21. Maskell DJ,
    22. Holmes MA
    . 2011. Meticillin-resistant Staphylococcus aureus with a novel mecA homologue in human and bovine populations in the UK and Denmark: a descriptive study. Lancet Infect Dis 11:595–603. doi:10.1016/S1473-3099(11)70126-8.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Holden MTG,
    2. Feil EJ,
    3. Lindsay JA,
    4. Peacock SJ,
    5. Day NPJ,
    6. Enright MC,
    7. Foster TJ,
    8. Moore CE,
    9. Hurst L,
    10. Atkin R,
    11. Barron A,
    12. Bason N,
    13. Bentley SD,
    14. Chillingworth C,
    15. Chillingworth T,
    16. Churcher C,
    17. Clark L,
    18. Corton C,
    19. Cronin A,
    20. Doggett J,
    21. Dowd L,
    22. Feltwell T,
    23. Hance Z,
    24. Harris B,
    25. Hauser H,
    26. Holroyd S,
    27. Jagels K,
    28. James KD,
    29. Lennard N,
    30. Line A,
    31. Mayes R,
    32. Moule S,
    33. Mungall K,
    34. Ormond D,
    35. Quail MA,
    36. Rabbinowitsch E,
    37. Rutherford K,
    38. Sanders M,
    39. Sharp S,
    40. Simmonds M,
    41. Stevens K,
    42. Whitehead S,
    43. Barrell BG,
    44. Spratt BG,
    45. Parkhill J
    . 2004. Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance. Proc Natl Acad Sci USA 101:9786–9791. doi:10.1073/pnas.0402521101.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Neoh HM,
    2. Cui L,
    3. Yuzawa H,
    4. Takeuchi F,
    5. Matsuo M,
    6. Hiramatsu K
    . 2008. Mutated response regulator graR is responsible for phenotypic conversion of Staphylococcus aureus from heterogeneous vancomycin-intermediate resistance to vancomycin-intermediate resistance. Antimicrob Agents Chemother 52:45–53. doi:10.1128/AAC.00534-07.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Kuroda M,
    2. Ohta T,
    3. Uchiyama I,
    4. Baba T,
    5. Yuzawa H,
    6. Kobayashi I,
    7. Cui L,
    8. Oguchi A,
    9. Aoki K,
    10. Nagai Y,
    11. Lian J,
    12. Ito T,
    13. Kanamori M,
    14. Matsumaru H,
    15. Maruyama A,
    16. Murakami H,
    17. Hosoyama A,
    18. Mizutani-Ui Y,
    19. Takahashi NK,
    20. Sawano T,
    21. Inoue R,
    22. Kaito C,
    23. Sekimizu K,
    24. Hirakawa H,
    25. Kuhara S,
    26. Goto S,
    27. Yabuzaki J,
    28. Kanehisa M,
    29. Yamashita A,
    30. Oshima K,
    31. Furuya K,
    32. Yoshino C,
    33. Shiba T,
    34. Hattori M,
    35. Ogasawara N,
    36. Hayashi H,
    37. Hiramatsu K
    . 2001. Whole genome sequencing of meticillin-resistant Staphylococcus aureus. Lancet 357:1225–1240. doi:10.1016/S0140-6736(00)04403-2.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    1. Baba T,
    2. Takeuchi F,
    3. Kuroda M,
    4. Yuzawa H,
    5. Aoki K,
    6. Oguchi A,
    7. Nagai Y,
    8. Iwama N,
    9. Asano K,
    10. Naimi T,
    11. Kuroda H,
    12. Cui L,
    13. Yamamoto K,
    14. Hiramatsu K
    . 2002. Genome and virulence determinants of high virulence community-acquired MRSA. Lancet 359:1819–1827. doi:10.1016/S0140-6736(02)08713-5.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Gillaspy A,
    2. Worrell V,
    3. Orvis J,
    4. Roe B,
    5. Dyer D,
    6. Iandolo J
    . 2006. The Staphylococcus aureus NCTC 8325 genome, p 381–412. In Fischetti V, Novick R, Ferretti J, Portnoy D, Rood J (ed), Gram-positive pathogens, vol 2. ASM Press, Washington, DC.
    OpenUrl
  45. 45.↵
    1. Baba T,
    2. Bae T,
    3. Schneewind O,
    4. Takeuchi F,
    5. Hiramatsu K
    . 2008. Genome sequence of Staphylococcus aureus strain Newman and comparative analysis of staphylococcal genomes: polymorphism and evolution of two major pathogenicity islands. J Bacteriol 190:300–310. doi:10.1128/JB.01000-07.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Herron-Olson L,
    2. Fitzgerald JR,
    3. Musser JM,
    4. Kapur V
    . 2007. Molecular correlates of host specialization in Staphylococcus aureus. PLoS One 2:e1120. doi:10.1371/journal.pone.0001120.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Schijffelen MJ,
    2. Boel CE,
    3. van Strijp JA,
    4. Fluit AC
    . 2010. Whole genome analysis of a livestock-associated methicillin-resistant Staphylococcus aureus ST398 isolate from a case of human endocarditis. BMC Genomics 11:376. doi:10.1186/1471-2164-11-376.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Holden MTG,
    2. Lindsay JA,
    3. Corton C,
    4. Quail MA,
    5. Cockfield JD,
    6. Pathak S,
    7. Batra R,
    8. Parkhill J,
    9. Bentley SD,
    10. Edgeworth JD
    . 2010. Genome sequence of a recently emerged, highly transmissible, multi-antibiotic- and antiseptic-resistant variant of methicillin-resistant Staphylococcus aureus, sequence type 239 (TW). J Bacteriol 192:888–892. doi:10.1128/JB.01255-09.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Highlander SK,
    2. Hulten KG,
    3. Qin X,
    4. Jiang H,
    5. Yerrapragada S,
    6. Mason EO, Jr..,
    7. Shang Y,
    8. Williams TM,
    9. Fortunov RM,
    10. Liu Y,
    11. Igboeli O,
    12. Petrosino J,
    13. Tirumalai M,
    14. Uzman A,
    15. Fox GE,
    16. Cardenas AM,
    17. Muzny DM,
    18. Hemphill L,
    19. Ding Y,
    20. Dugan S,
    21. Blyth PR,
    22. Buhay CJ,
    23. Dinh HH,
    24. Hawes AC,
    25. Holder M,
    26. Kovar CL,
    27. Lee SL,
    28. Liu W,
    29. Nazareth LV,
    30. Wang Q,
    31. Zhou J,
    32. Kaplan SL,
    33. Weinstock GM
    . 2007. Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus. BMC Microbiol 7:99. doi:10.1186/1471-2180-7-99.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Croucher NJ,
    2. Harris SR,
    3. Fraser C,
    4. Quail MA,
    5. Burton J,
    6. van der Linden M,
    7. McGee L,
    8. von Gottberg A,
    9. Song JH,
    10. Ko KS,
    11. Pichon B,
    12. Baker S,
    13. Parry CM,
    14. Lambertsen LM,
    15. Shahinas D,
    16. Pillai DR,
    17. Mitchell TJ,
    18. Dougan G,
    19. Tomasz A,
    20. Klugman KP,
    21. Parkhill J,
    22. Hanage WP,
    23. Bentley SD
    . 2011. Rapid pneumococcal evolution in response to clinical interventions. Science 331:430–434. doi:10.1126/science.1198545.
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    1. McKenna A,
    2. Hanna M,
    3. Banks E,
    4. Sivachenko A,
    5. Cibulskis K,
    6. Kernytsky A,
    7. Garimella K,
    8. Altshuler D,
    9. Gabriel S,
    10. Daly M,
    11. DePristo MA
    . 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303. doi:10.1101/gr.107524.110.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Price MN,
    2. Dehal PS,
    3. Arkin AP
    . 2009. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol 26:1641–1650. doi:10.1093/molbev/msp077.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. Stamatakis A
    . 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312–1313. doi:10.1093/bioinformatics/btu033.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Ruimy R,
    2. Maiga A,
    3. Armand-Lefevre L,
    4. Maiga I,
    5. Diallo A,
    6. Koumare AK,
    7. Ouattara K,
    8. Soumare S,
    9. Gaillard K,
    10. Lucet JC,
    11. Andremont A,
    12. Feil EJ
    . 2008. The carriage population of Staphylococcus aureus from Mali is composed of a combination of pandemic clones and the divergent Panton-Valentine leukocidin-positive genotype ST152. J Bacteriol 190:3962–3968. doi:10.1128/JB.01947-07.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Kelley DR,
    2. Schatz MC,
    3. Salzberg SL
    . 2010. Quake: quality-aware detection and correction of sequencing errors. Genome Biol 11:R116. doi:10.1186/gb-2010-11-11-r116.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Luo R,
    2. Liu B,
    3. Xie Y,
    4. Li Z,
    5. Huang W,
    6. Yuan J,
    7. He G,
    8. Chen Y,
    9. Pan Q,
    10. Liu Y,
    11. Tang J,
    12. Wu G,
    13. Zhang H,
    14. Shi Y,
    15. Liu Y,
    16. Yu C,
    17. Wang B,
    18. Lu Y,
    19. Han C,
    20. Cheung DW,
    21. Yiu SM,
    22. Peng S,
    23. Xiaoqian Z,
    24. Liu G,
    25. Liao X,
    26. Li Y,
    27. Yang H,
    28. Wang J,
    29. Lam TW,
    30. Wang J
    . 2012. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience 1:18. doi:10.1186/2047-217X-1-18.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Hyatt D,
    2. Chen G,
    3. Locascio PF,
    4. Land ML,
    5. Larimer FW,
    6. Hauser LJ
    . 2010. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi:10.1186/1471-2105-11-119.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    . 1990. Basic local alignment search tool. J Mol Biol 215:403–410. doi:10.1016/S0022-2836(05)80360-2.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Enright AJ,
    2. Van Dongen S,
    3. Ouzounis CA
    . 2002. An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res 30:1575–1584. doi:10.1093/nar/30.7.1575.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Jarraud S,
    2. Cozon G,
    3. Vandenesch F,
    4. Bes M,
    5. Etienne J,
    6. Lina G
    . 1999. Involvement of enterotoxins G and I in staphylococcal toxic shock syndrome and staphylococcal scarlet fever. J Clin Microbiol 37:2446–2449.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. McClure JA,
    2. Conly JM,
    3. Lau V,
    4. Elsayed S,
    5. Louie T,
    6. Hutchins W,
    7. Zhang K
    . 2006. Novel multiplex PCR assay for detection of the staphylococcal virulence marker Panton-Valentine leukocidin genes and simultaneous discrimination of methicillin-susceptible from -resistant staphylococci. J Clin Microbiol 44:1141–1144. doi:10.1128/JCM.44.3.1141-1144.2006.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Sung JM,
    2. Lloyd DH,
    3. Lindsay JA
    . 2008. Staphylococcus aureus host specificity: comparative genomics of human versus animal isolates by multi-strain microarray. Microbiology 154:1949–1959. doi:10.1099/mic.0.2007/015289-0.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    1. Jarraud S,
    2. Mougel C,
    3. Thioulouse J,
    4. Lina G,
    5. Meugnier H,
    6. Forey F,
    7. Nesme X,
    8. Etienne J,
    9. Vandenesch F
    . 2002. Relationships between Staphylococcus aureus genetic background, virulence factors, agr groups (alleles), and human disease. Infect Immun 70:631–641. doi:10.1128/IAI.70.2.631-641.2002.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Johnson WM,
    2. Tyler SD,
    3. Ewan EP,
    4. Ashton FE,
    5. Pollard DR,
    6. Rozee KR
    . 1991. Detection of genes for enterotoxins, exfoliative toxins, and toxic shock syndrome toxin 1 in Staphylococcus aureus by the polymerase chain reaction. J Clin Microbiol 29:426–430.
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    1. Diep B,
    2. Stone G,
    3. Basuino L,
    4. Graber C,
    5. Miller A,
    6. des Etages S,
    7. Jones A,
    8. Palazzolo-Ballance A,
    9. Perdreau-Remington F,
    10. Sensabaugh G,
    11. DeLeo F,
    12. Chambers H
    . 2008. 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 197:1523–1530. doi:10.1086/587907.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Robinson DA,
    2. Enright MC
    . 2004. Evolution of Staphylococcus aureus by large chromosomal replacements. J Bacteriol 186:1060–1064.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe
David M. Aanensen, Edward J. Feil, Matthew T. G. Holden, Janina Dordel, Corin A. Yeats, Artemij Fedosejev, Richard Goater, Santiago Castillo-Ramírez, Jukka Corander, Caroline Colijn, Monika A. Chlebowicz, Leo Schouls, Max Heck, Gerlinde Pluister, Raymond Ruimy, Gunnar Kahlmeter, Jenny Åhman, Erika Matuschek, Alexander W. Friedrich, Julian Parkhill, Stephen D. Bentley, Brian G. Spratt, Hajo Grundmann ESCMID Study Group on Molecular Epidemiological Markers (ESGEM), European Staphylococcal Reference Laboratory Working Group
mBio May 2016, 7 (3) e00444-16; DOI: 10.1128/mBio.00444-16

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe
David M. Aanensen, Edward J. Feil, Matthew T. G. Holden, Janina Dordel, Corin A. Yeats, Artemij Fedosejev, Richard Goater, Santiago Castillo-Ramírez, Jukka Corander, Caroline Colijn, Monika A. Chlebowicz, Leo Schouls, Max Heck, Gerlinde Pluister, Raymond Ruimy, Gunnar Kahlmeter, Jenny Åhman, Erika Matuschek, Alexander W. Friedrich, Julian Parkhill, Stephen D. Bentley, Brian G. Spratt, Hajo Grundmann ESCMID Study Group on Molecular Epidemiological Markers (ESGEM), European Staphylococcal Reference Laboratory Working Group
mBio May 2016, 7 (3) e00444-16; DOI: 10.1128/mBio.00444-16
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511