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
    • Current Issue
    • Accepted Manuscripts
    • Archive
    • Minireviews
  • 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 JCM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • 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
Journal of Clinical Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • Archive
    • Minireviews
  • 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 JCM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Epidemiology

A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak

Mitchell J. Sullivan, Deena R. Altman, Kieran I. Chacko, Brianne Ciferri, Elizabeth Webster, Theodore R. Pak, Gintaras Deikus, Martha Lewis-Sandari, Zenab Khan, Colleen Beckford, Angela Rendo, Flora Samaroo, Robert Sebra, Ramona Karam-Howlin, Tanis Dingle, Camille Hamula, Ali Bashir, Eric Schadt, Gopi Patel, Frances Wallach, Andrew Kasarskis, Kathleen Gibbs, Harm van Bakel
John P. Dekker, Editor
Mitchell J. Sullivan
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Deena R. Altman
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USADepartment of Medicine, Division of Infectious Diseases, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Deena R. Altman
Kieran I. Chacko
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brianne Ciferri
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elizabeth Webster
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Theodore R. Pak
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Theodore R. Pak
Gintaras Deikus
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martha Lewis-Sandari
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zenab Khan
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Colleen Beckford
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Angela Rendo
Department of Pathology, Clinical Microbiology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Flora Samaroo
Department of Pathology, Clinical Microbiology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert Sebra
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ramona Karam-Howlin
Department of Medicine, Division of Infectious Diseases, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tanis Dingle
Department of Pathology, Clinical Microbiology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Camille Hamula
Department of Pathology, Clinical Microbiology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ali Bashir
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric Schadt
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gopi Patel
Department of Medicine, Division of Infectious Diseases, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Frances Wallach
Department of Medicine, Division of Infectious Diseases, Northwell Long Island Jewish Medical Center, New Hyde Park, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew Kasarskis
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kathleen Gibbs
Division of Neonatology and Department of Pediatrics, The Children’s Hospital of Philadelphia, The University of Pennsylvania, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Harm van Bakel
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York City, New York, USAIcahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York City, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Harm van Bakel
John P. Dekker
National Institute of Allergy and Infectious Diseases
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/JCM.01261-19
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices. In this study, we established a long-read technology-based WGS screening program of all first-episode methicillin-resistant Staphylococcus aureus (MRSA) blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads enabled detailed characterization of an outbreak lasting several months of a CC5/ST105/USA100 clone among 18 infants in a neonatal intensive care unit (NICU). Available hospital-wide genome surveillance data traced the origins of the outbreak to three patients admitted to adult wards during a 4-month period preceding the NICU outbreak. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its progression, which was characterized by multiple subtransmissions and likely precipitated by equipment sharing between adults and infants. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance, and persistence. This included a DNA recognition domain recombination in the hsdS gene of a type I restriction modification system that altered DNA methylation. Transcriptome sequencing (RNA-Seq) profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall, our findings demonstrate the utility of long-read sequencing for hospital surveillance and for characterizing accessory genomic elements that may impact MRSA virulence and persistence.

INTRODUCTION

Health care-associated infections (HAI) with methicillin-resistant Staphylococcus aureus (MRSA) are common, impair patient outcomes, and increase health care costs (1, 2). MRSA is highly clonal, and much of our understanding of its dissemination has relied on lower-resolution molecular strain typing methods such as pulsed-field gel electrophoresis (PFGE), S. aureus protein A (spa) typing, and multilocus sequence typing (MLST) (3) and typically includes characterization of accessory genome elements that define certain lineages and are implicated in their virulence. Examples of the latter include the arginine catabolic mobile element (ACME), S. aureus pathogenicity island 5 (SaPI5), and the Panton-Valentine leukocidin (PVL)-carrying ϕSa2 prophage in the community-associated (CA) CC8/USA300 lineage (4, 5). Molecular typing facilitates rapid screening but has limited resolution to identify transmissions within clonal lineages. Moreover, genetic changes can lead to alteration or loss of typing elements (6–9). As such, whole-genome sequencing (WGS) has emerged as the gold standard for studying lineage evolution and nosocomial outbreaks (10, 11). Transmission analysis with WGS has been performed largely retrospectively to date (12–15), although prospective screening with resulting interventions has also been described (10, 16).

In addition to lineage and outbreak analysis, WGS has furthered our understanding of S. aureus pathogenicity by delineating virulence and drug resistance determinants (17, 18), including those related to adaptation to the hospital environment (17, 19). Many of these elements are found in nonconserved “accessory” genome elements that include endogenous prophages, mobile genetic elements (MGEs), and plasmids (20, 21). The repetitive nature of many of these elements means that they are often fragmented and/or incompletely represented in most WGS studies to date due to limitations of commonly used short-read sequencing technologies, curbing insights into their evolution (21). Recent advances in throughput of long-read sequencing technologies now enable routine assembly of complete genomes (22, 23) and analysis of core and accessory genome elements (13, 18), including DNA methylation patterns (24), but these technologies have not yet been widely used for prospective MRSA surveillance.

Here, we describe the results of a complete genome-based screening program of MRSA blood isolates. During a 16-month period, we obtained finished-quality genomes for first blood isolates from all bacteremic patients. In addition to providing detailed contemporary insights into prevailing lineages and genome characteristics, we characterized widespread variation across accessory genome elements, impacting loci encoding virulence and resistance factors, including those commonly used as molecular strain typing markers. We also found multiple transmission events not recognized based on epidemiological information. During an outbreak event in the neonatal intensive care unit (NICU), our program was able to provide actionable information that discriminated outbreak-related transmissions, identified individual subtransmission events, and traced the NICU outbreak origin to adult hospital wards. Finally, comparative genome and gene expression analyses of the outbreak clone to hospital background strains identified genetic and epigenetic changes, including acquisition of accessory genome elements, which may have contributed to the persistence of the outbreak clone.

MATERIALS AND METHODS

Ethics statement.This study was reviewed and approved by the Institutional Review Board of the Icahn School of Medicine at Mount Sinai and the MSH Pediatric Quality Improvement Committee at Mount Sinai Hospital.

Case review.An investigation of the characteristics of the patients included review of existing medical records for relevant clinical data. Unique ventilator identification numbers and the real-time location system (RTLS) enabled mapping of ventilator locations over time.

Bacterial isolate identification and susceptibility testing.Isolates were grown and identified as part of standard clinical testing procedures in the Mount Sinai Hospital Clinical Microbiology Laboratory (CML, and stored in tryptic soy broth (TSB) with 15% glycerol at −80°C. Species confirmation was performed with matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) (Bruker Biotyper; Bruker Daltonics). Vitek 2 (bioMérieux) automated broth microdilution antibiotic susceptibility profiles were obtained for each isolate according to Clinical and Laboratory Standards Institute (CLSI) 2015 guidelines and reported according to CLSI guidelines (25). Susceptibility to mupirocin was determined by Etest (bioMérieux), and susceptibility to chlorhexidine was tested with discs (Hardy) impregnated with 5 μl of a 20% chlorhexidine gluconate solution (Sigma-Aldrich).

DNA preparation and sequencing.For each isolate, single colonies were selected and grown separately on tryptic soy agar (TSA) plates with 5% sheep blood (blood agar) (Thermo Fisher Scientific) under nonselective conditions. After growth overnight, cells underwent high-molecular-weight DNA extraction using the Qiagen DNeasy Blood & Tissue kit (Qiagen) according to the manufacturer’s instructions, with modified lysis conditions: bacterial cells were pelleted and resuspended in 300 μl enzymatic lysis buffer (20 mM Tris-Cl [pH 8.0], 2 mM sodium EDTA), 3 μl of 100 mg/ml RNase A (AM2286; Ambion), and 10 μl of 100 mg/ml lysozyme (L1667-1G; Sigma) for 30 min at 37°C, followed by the addition of 25 μl proteinase K (20 mg/ml) (Qiagen) and 300 μl Buffer AL and further incubation for 1 h at 56°C, after which two rounds of bead beating (BioSpec) were performed of 1 min each using 0.1-mm silica beads (MP Bio) (13).

Quality control, DNA quantification, library preparation, and sequencing were performed as described previously (13). Briefly, DNA was gently sheared using Covaris G-tube spin columns into ∼20,000-bp fragments and end-repaired before ligating SMRTbell adapters (Pacific Biosciences). The resulting library was treated with an exonuclease cocktail to remove unligated DNA fragments, followed by two additional purification steps with AMPure XP beads (Beckman Coulter) and Blue Pippin (Sage Science) size selection to deplete SMRTbells of <7,000 bp. Libraries were then sequenced using P5 enzyme chemistry on the Pacific Biosciences RS-II platform to >200× genome-wide coverage.

Complete genome assembly and finishing.PacBio SMRT sequencing data were assembled using a custom genome assembly and finishing pipeline (26). Briefly, sequencing data were first assembled with HGAP3 version 2.2.0 (22). Contigs with less than 10× coverage and small contigs that were completely encompassed in larger contigs were removed. Remaining contigs were circularized and reoriented to the origin of replication (ori) using Circlator (27) and aligned to the nonredundant nucleotide collection using BLAST+ (28) to identify plasmid sequences. In cases where chromosomes or plasmids did not assemble into complete circularized contigs, manual curation was performed using Contiguity (29). Genes were annotated using PROKKA (30) and visualized using ChromoZoom (31) and the Integrated Genome Browser (IGB) (32). Interproscan (33) was used to annotate protein domains and Gene Ontology (GO) categories for annotated genes.

Resolution of large genomic inversions.To resolve inversion events catalyzed by two prophage elements (Staphylococcus ϕSa1 and Staphylococcus aureus ϕSa5) with large (>40 kbp) nearly identical regions present in some of the assembled genomes, we developed a phasing approach that took advantage of unique variants present in each element (see Fig. S1 in the supplemental material). Raw (i.e., uncorrected) PacBio reads were first mapped to one of the repeat copies using BWA-MEM (34). Variants were then called with Freebayes (35), and high-quality single nucleotide variants with two distinct alleles of approximately equal reads coverage were identified. Analogous to procedures used in haplotype phasing, we then determined which variant alleles were colocated in the same repeat element: if at least three-quarters of the raw reads containing a particular allele also encompassed distinct allele(s) of a neighboring variant(s), the alleles were considered linked. In all cases, this resulted in two distinct paths through the repeated prophage elements that were each linked to a unique sequence flanking each repeat. We then used this information to correct assembly errors and identify bona fide inversion events between isolate genomes. Final verification of corrected assembly was performed by examining the phasing of the raw reads with HaploFlow (36).

Phylogenetic reconstruction and molecular typing.Phylogenetic analyses were based on whole-genome alignments with parsnp (37), using the filter for recombination. The VCF file of all variants identified by parsnp was then used to determine pairwise single nucleotide variant (SNV) distances between the core genomes of all strains. For visualization of the whole-genome alignments, isolate genomes were aligned using sibelia (38) and processed by ChromatiBlock (https://github.com/mjsull/chromatiblock). Indels, structural variants, and SNVs in pairwise complete genome alignments were identified using NucDiff (39).

The multilocus sequence type was determined from whole-genome sequences using the RESTful interface to the PubMLST S. aureus database (40). Staphylococcal protein A (spa) typing was performed using a custom script (https://github.com/mjsull/spa_typing). SCCmec typing was done using SCCmecFinder (41). Changes to ACME and SaPI5 were determined using BLASTn and Easyfig. Presence or absence of genes in each locus was determined using BLASTx (42), and a gene was considered to be present if 90% of the reference sequence was aligned with at least 90% identity. Prophage regions were detected using PHASTER. Each region was then aligned to a manually curated database of S. aureus phage integrases using BLASTx to identify their integrase group.

Annotation of antibiotic resistance determinants.Antibiotic resistance genes and variants were annotated by comparing to a manually curated database of 39 known S. aureus resistance determinants for 17 antibiotics compiled from the literature. BLAST (42) was used to identify the presence of genes in each isolate genome, with a sequence identity cutoff of ≥90% and an E value cutoff of ≤1e−10. Resistance variants were identified by BLAST alignment to the reference sequence of the antibiotic resistance determinant. Only exact matches to variants identified in literature were considered.

Identification of NICU outbreak subgroups.Changes between each outbreak isolate and the p133 reference isolate were identified using GWviz (https://github.com/mjsull/GWviz), which uses nucdiff (39) to identify all genomic variants between pairs of strains and then uses PROKKA gene annotations to determine the effect of the change on coding regions. nucdiff in turn uses nucmer to find alignments between two genomes and then identifies large structural rearrangements by looking at the organization of nucmer alignments and smaller changes such as SNVs or indels by finding differences between the aligned regions. To confirm indel changes between isolates in the outbreak, small indels were filtered based on underlying read data. Briefly, raw PacBio reads were aligned back to each outbreak genome assembly using BWA-MEM (34). GWviz was then used to determine the number and proportion of raw reads supporting variants in each strain. Variants were selected for further delineation of outbreak subgroups if they were present in two or more isolate genomes and supported by at least ten raw reads in each genome, of which at least 75% confirmed the variant.

A graph of SNV distances between isolates was obtained from a multiple alignment of all outbreak isolates. The minimum spanning tree was then constructed using the minimum spanning tree functionality in the Python library networkX (https://networkx.github.io/).

Identification of genetic variants unique to the NICU outbreak clone.To determine SNVs unique to the outbreak isolate the marginal ancestral states of the ST105 isolates were determined using RAxML (43) from a multiple alignment of all ST105s generated using parsnp. We identified all SNVs that had accumulated from the most recent common ancestor of the outbreak strain and the closest related nonoutbreak ST105, and the most recent common ancestor (MRCA) of all outbreak strains. SNVs causing nonsynonymous mutations or changes to the promoter region of a gene (defined as <500 bp upstream of the start site) were plotted. Orthology was assigned using BLASTkoala (44).

Core and accessory gene content in ST105 outbreak and nonoutbreak strains was determined using ROARY. Genes found in more than two outbreak strains and less than 33% of the other ST105 genomes were then plotted along with select methylation data. Phylogenetic reconstruction of ST105 was performed using parsnp, and the resulting tree and gene presence information was visualized using m.viridis.py (https://github.com/mjsull/m.viridis) which uses the python ETE toolkit (45).

DNA methylation profiling.SMRT raw reads were mapped to the assembled genomes and processed using smrtanalysis v5.0 (Pacific Biosciences, Menlo Park, CA). Interpulse durations (IPDs) were measured and processed as previously described (24, 46) to detect modified N6-methyladenine (m6A) nucleotides.

RNA preparation and sequencing.For RNA extraction, overnight cultures in tryptic soy broth (TSB) were diluted (optical density at 600 nm [OD600] of 0.05), grown to late-log phase (OD600 of ∼0.80) in TSB, and stabilized in RNALater (Thermo Fisher). Total RNA was isolated and purified using the RNeasy Mini kit (Qiagen) according to the manufacturer’s instructions, except that two cycles of 2-min bead beating with 1 ml of 0.1-mm silica beads in a mini bead beater (BioSpec) were used to disrupt cell walls. Isolated RNA was treated with 1 μl (1 unit) of Baseline Zero DNase (Epicentre) at 37°C for 30 min, followed by rRNA depletion using the Epicenter Ribo-Zero Magnetic Gold kit (Illumina), according to the manufacturer’s instructions.

RNA quality and quantity were assessed using the Agilent Bioanalyzer and Qubit RNA Broad Range assay kit (Thermo Fisher), respectively. Barcoded directional RNA sequencing libraries were prepared using the TruSeq Stranded Total RNA Sample Preparation kit (Illumina). Libraries were pooled and sequenced on the Illumina HiSeq platform in a 100-bp single-end read run format with six samples per lane.

Differential gene expression analysis.Raw reads were first trimmed by removing Illumina adapter sequences from 3′ ends using cutadapt (47) with a minimum match of 32 bp and allowing for 15% error rate. Trimmed reads were mapped to the reference genome using Bowtie2 (48), and htseq-count (49) was used to produce strand-specific transcript count summaries. Read counts were then combined into a numeric matrix and used as input for differential gene expression analysis with the Bioconductor EdgeR package (50). Normalization factors were computed on the data matrix using the weighted trimmed mean of M values (TMM) method (51). Data were fitted to a design matrix containing all sample groups, and pairwise comparisons were performed between the groups of interest. P values were corrected for multiple testing using the Benjamin-Hochberg (BH) method and used to select genes with significant expression differences (q < 0.05).

Data availability.All genome data and assemblies are available in GenBank under accession numbers CP030375 to CP030714, QNXD00000000, QNXE00000000, QNXF00000000, QNXG00000000, QNXH00000000, QNXI00000000, QNXJ00000000, QNXK00000000, QNXL00000000, QNXM00000000, and VZDK00000000 (see Table S1).

RESULTS

Complete genome surveillance reveals accessory genome variation among clonal MRSA lineages.To characterize the genetic diversity of MRSA blood infections at The Mount Sinai Hospital (MSH) in New York City, NY, we sequenced the first positive isolate from all 132 MSH inpatients diagnosed with MRSA bacteremia between fall 2014 and winter 2015. Single-molecule real-time (SMRT) long-read length RS-II WGS was used to obtain finished-quality chromosomes for 122 of 132 isolates (92%), along with 145 unique plasmids across isolates (see Table S1 in the supplemental material). The remaining isolates had one or more chromosomal contigs that could not be closed with available long-read sequencing data. We reconstructed a phylogeny from a multigenome alignment (Fig. 1A; see also Fig. S2A), which identified two major clades corresponding to S. aureus clonal complexes 8 (CC8; 45.5% of isolates) and 5 (CC5; 50% of isolates) based on the prevailing multilocus sequence types (STs) in each clade (ST8 and ST105/ST5, respectively). The CC8 isolates further partitioned among the endemic community-associated (CA) USA300 (80%) and the hospital-associated (HA) USA500 (20%) lineages (Fig. 1B), while CC5 isolates mainly consisted of USA100 (75.8%) and USA800 (15.2%) HA lineages (Fig. 1C). Overall, the phylogeny was consistent with the major MRSA lineages found in New York City, NY, and the United States (52).

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

Phylogeny of MRSA bacteremia surveillance isolates. (A) Maximum likelihood phylogenetic tree based on SNV distances in core genome alignments of 132 primary MRSA bacteremia isolates. CC8 and CC5 clades are shaded in red and blue, respectively. Multilocus sequence types (MLST) for each branch are shown as colored blocks, with a key at the bottom left. (B) Enlarged version of the CC8 clade from panel A. The isolate identifier is indicated next to each branch, together with blocks denoting the spa type, SCCmec type, and the presence (blue) or absence (yellow) of intact ACME, lukFS, and SaPI5 loci. The ACME type is indicated in each box. The lukFS locus is represented by two blocks indicating the presence of lukF and lukS, respectively. (C) Same as panel B, but for the CC5 clade. *, spa type II isolate with an inserted element in the locus. Four transmission events between patients are highlighted in red and labeled T1 to T4. Scale bars indicate the number of substitutions per site in the phylogeny.

We further examined larger (>500 bp) structural variation that may be missed by short-read-based WGS approaches (13, 14, 53). The multigenome alignment indicated that between 80.8% and 88.9% of the sequence in each genome was contained in core syntenic blocks shared among all 132 genomes see (Fig. S2A). Another 9.5% to 16.8% was contained in accessory blocks found in at least two but not all genomes. Many of these accessory genome elements were lineage specific and associated with prophage regions and plasmids (Fig. S2B). Finally, 0.8% to 4.5% of the sequence was not found in syntenic blocks and included unique elements gained by individual isolates. The extent of core and accessory genome variability impacted loci that are commonly used for molecular strain typing. Divergence from the dominant spa type was apparent in 8 (13.3%) of CC8 and 9 (13.6%) of CC5 lineage isolates. MLST loci were more stable in comparison with changes in 1.5% and 7.6% of isolates in each lineage, respectively. Notably, there were also widespread changes at ACME, PVL, and SaPI5 (Fig. 1B) in USA300 isolates, which are signature elements of this CA lineage (4, 5); 33.3% (16 of 48) either carried inactivating mutations or had partially or completely lost one or more elements (Fig. 1B). The multiple independent events of ACME, PVL, and SaPI5 loss throughout the USA300 clade may reflect its ongoing adaptation to hospital environments, as these elements are typically absent in HA lineages. Interestingly, we found one case of a PVL-positive USA100 isolate (Fig. 1C) that may have resulted from homologous recombination between a ϕSa2 and ϕSa2 PVL prophage (see Fig. S3). Thus, complete genomes of MRSA blood isolates demonstrate the mobility of the accessory genome in ways that impact commonly used S. aureus lineage definitions.

Identification of transmission events among adults and an outbreak in the NICU.We next compared isolate genomes to identify transmissions between patients. To establish similarity thresholds for complete genomes obtained from long-read SMRT sequencing data, considering both intrahost diversity and genetic drift, we first examined baseline single nucleotide variant (SNV) distances within each lineage. Median pairwise genome differences ranged from 101 SNVs for USA800 to 284 SNVs for USA100 (see Fig. S4A). We also examined the extent of divergence among 30 bacteremia isolate pairs collected within a span of 1 month to 1.4 years from individual patients. Pairwise distances for within-patient isolates were substantially lower than the median for each lineage (Fig. S4B to E), consistent with persistent carriage of the same clone (54, 55), with no more than 10 SNVs separating isolate pairs. Small (<5 bp) indels were more common than SNVs and mostly associated with homopolymer regions that can be problematic to resolve with third-generation sequencing technologies, indicating that they likely reflected sequencing errors. Notably, several patients showed variation between isolates collected within a span of several days (Fig. S4B to E), indicative of intrahost genetic diversity. As such, we considered intrahost diversity and genetic drift in aggregate and set a conservative distance of ≤7 SNVs to define transmission events in our genome phylogeny. At this threshold, we identified one USA300 and three USA100 transmissions involving six adults and three infants (Fig. 1B and C, labeled T1 to T4). Complete pairwise genome alignments for each event confirmed the absence of structural variants. In the USA300 transmission case (T1), the presumed index patient p5 was bacteremic with the same clone on two occasions ∼3 months apart (see Fig. S5). The recipient (p33) was admitted to the same ward and overlapped with p5 for 7 days prior to the time of bacteremia. The USA100 isolates in transmission T2 were collected ∼4 months apart, and although the patients had overlapping stays, they did not share a ward or other clear epidemiological links (Fig. S5). In transmission T3, patients shared a ward for several days (Fig. S5).

The final transmission involving 3 infants (T4) was part of a larger outbreak in the NICU, where positive clinical MRSA cultures from three infants within 5 weeks had prompted an investigation and consultation with the New York State Department of Health (NYSDOH). During 4 months, an additional 41 clinical and surveillance cultures from 20 infants tested positive for MRSA, bringing the total to 46 isolates from 22 infants. Three further isolates were obtained from incubators and an intravenous (i.v.) box, from a total of 123 environmental swabs (2.4%). Positive nasal surveillance cultures were also obtained from 2 of 130 (1.5%) health care workers (HCWs) who had provided direct care to newly MRSA-colonized infants. The NYSDOH performed PFGE on 22 isolates, of which 14 patients and 3 environmental isolates had nearly indistinguishable band patterns (data not shown). This included p90 and p110 in transmission T4 (Fig. 1C) (p125 was not tested). The USA100 (ST105) outbreak clone was resistant to fluoroquinolones, clindamycin, gentamicin, and mupirocin and susceptible to vancomycin, trimethoprim-sulfamethoxazole, and doxycycline (see Fig. S6 and Table S1). This pattern was uncommon (18.2%) among USA100 isolates in our study and was therefore used as an initial screening criteria for cases. None of the HCW isolates matched the MLST or antibiogram of the outbreak clone, and both staff members were successfully decolonized with nasal mupirocin and chlorhexidine gluconate (CHG) baths.

Complete genome surveillance resolves outbreak origin and progression.During the outbreak, we expanded our genomic screening program to include the first isolate of suspected outbreak cases. From day 354 onwards, we obtained 23 additional complete genomes (Table S1). Of these, 19 genomes from 16 infants and three environmental isolates matched the ST105 outbreak strain type, bringing the total to 22 outbreak genomes from 16 infants and the environment. The infection prevention team and the NYSDOH were informed of isolate genomes meeting our transmission threshold within 10 to 14 days of a positive test, which helped delineate the final case set and determine when the outbreak ended. Further analyses were performed retrospectively to reconstruct the chain of events that initiated and sustained the outbreak. To this end, we first reconstructed a phylogenetic tree based on core genome alignments of all ST105 isolates in our study, which grouped all 22 isolates with matching antibiograms and/or PFGE patterns in one well-defined clade (Fig. 2A). Surprisingly, this clade also contained 3 MRSA isolates obtained from adult bacteremia patients in other hospital wards prior to the first NICU case. The outbreak clade genomes were ≤15 SNVs apart, and the clade as a whole differed from other ST105 isolates by ≥41 SNVs. We therefore considered the 3 adult isolates to be part of a larger clonal outbreak that spanned 7 months. The availability of complete genome sequences provided additional genomic variants that contributed to strain diversity within the outbreak clade, including indels, structural variants, and a large megabase-size inversion (Fig. 2A). Based on these variant patterns, we distinguished 4 distinct subgroups. A minimum spanning tree based only on core genome SNVs that is more commonly used in outbreak investigations using short-read sequencing data largely recapitulated the same grouping (Fig. 2B) but with reduced resolution between subgroups A and B.

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

NICU outbreak subgroups and association with adult bacteremia patients. (A) Maximum likelihood phylogenetic tree based on SNV distances in core genome alignments of 31 ST105 primary bacteremia isolates (black) and 25 outbreak isolates (red). The core genome makes up 76.1% to 82.6% of each genome. The scale bar indicates the number of substitutions per site. The patient (p) or environmental (e) isolate identifier is shown next to each branch (a/b suffixes indicate multiple isolates from the same patient). Variants present in two or more NICU outbreak isolates, derived from full-length pairwise alignments to the p133 genome, are shown as colored boxes. Variants are colored according to outbreak subgroups inferred from common variant patterns, as indicated on the right. For each variant the genomic location, affected genes, and type of mutation is shown above the matrix. A 2-Mbp inversion in the adult isolates and a 2,411-bp region containing two substitutions and a deletion in subgroup Bare highlighted in the location bar in orange and purple, respectively. (B) Minimum spanning tree of the 25 outbreak isolates based on SNVs identified in the complete genome alignment of all ST105 isolates. The 15 labeled nodes represent individual isolates. The larger central node corresponds to ten isolates with identical core genomes, which includes the p133 reference. Nodes are colored according to the outbreak subgroups shown in panel A. Numbers at edges represent core genome SNV distances.

We then used the available epidemiology and genomic data to reconstruct an outbreak timeline (Fig. 3A). The three initial adult cases had overlapping stays and shared wards, and their isolates clustered together in subgroup A. Several of the earliest clinical isolates from infants p141, p150, and p151 that coincided with the spread of the outbreak to the NICU were not available for genomic analysis (marked X in Fig. 3A). The missing isolate from p141 was susceptible to gentamicin and differed from the PFGE pattern of the outbreak clone by five bands. The other two missing isolates from p150 and 151 matched the outbreak clone antibiogram and were therefore considered to be part of the outbreak. Subsequent cases were identified by positive surveillance cultures on days 357 to 386, and their isolates clustered in subgroup C. All but one of the infants in this subgroup stayed in NICU room 2 before or at the time of culture positivity. The three positive environmental isolates were also obtained from this room, suggesting that a local bioburden led to a high volume of colonized infants in a short time. Construction in the NICU and a resulting disruption of infection prevention practices was believed to play a role in the initial transmissions of MRSA.

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

Timeline of the NICU outbreak. (A) Overview of outbreak patient stays and isolates collected during the NICU outbreak. Rows correspond to patients with admission periods shown as horizontal bars. Solid fill patterns denote NICU stays and striped patterns indicate stays in other MSH wards. Fill colors correspond to NICU rooms (solid) or hospital wards (striped). Clinical or surveillance isolates collected during each stay are indicated by symbols, with a key shown below. Patient identifiers and isolate symbols are colored by outbreak subgroup. Timeline scale and key interventions are shown at the top. SRV, start of biweekly surveillance cultures; TC, terminal cleaning; SIM, in situ simulation. (B) Same as panel A, but with ventilator movements between patients and locations overlaid as lines. Ventilators are numbered and shown in distinct colors. Solid lines correspond to periods that a ventilator was in use by an outbreak patient. Dashed lines indicate when a ventilator was present in the NICU but not used by an outbreak patient. Dotted lines indicate when a ventilator was not in use by an outbreak patient and not present in the NICU. Background colors are muted to facilitate tracking of ventilator movements.

The increase in new cases on surveillance prompted a terminal clean (TC) of the NICU on day 395. During this time, all infants were temporarily transferred to two different locations. Infant p148 who was colonized with the outbreak clone was placed across the hall from p141 in the pediatric intensive care unit (PICU). A positive surveillance culture in the same subgroup (B) as p141 was obtained for p148 shortly afterwards (Fig. 3A), suggesting that a transmission had occurred during the TC. New positive surveillance cultures were subsequently found for three additional infants (p142, p143, and p146). Each had been admitted after the TC and stayed in room 3 before or at the time of culture positivity. Their isolates comprised subgroup D, suggesting that the outbreak clone spread to this location from the closely related subgroup C linked to room 2 (Fig. 2B and 3A). Thus, each outbreak subgroup (A to D) was associated with a specific area (adult wards, PICU, and NICU rooms 2 and 3, respectively), indicating that location sharing was a dominant factor in the spread of the outbreak clone.

The continued transmissions after the first TC prompted in situ simulation and a second TC (Fig. 3A). The simulation efforts reinforced the importance of compliance to infection prevention strategies, patient cohorting, enhanced environmental disinfection, and limiting patient census to decrease bioburden (56). Only one new case (p124) was detected after the second TC. Infant p124 was located the PICU at the time of detection, and based on the genomic profile (subgroup C) and earlier positive isolates, the transmission was believed to have occurred prior to the final TC and in situ simulation. As such, the workflow improvements were effective in halting the outbreak. The weekly surveillance cultures ended after three consecutive weeks of negative cultures (day 452). The last colonized patient was discharged 2 months later, and we did not detect the outbreak clone in our hospital-wide genomic screening program in the subsequent 2 years. While the majority of cases were positive by surveillance, there was morbidity related to the outbreak; five infants developed clinical infections, with three bacteremias, one pneumonia, and one surgical site infection. There were no deaths related to the outbreak.

Ventilator sharing implicated in the origin and progression of the NICU outbreak.Location and HCW sharing could not account for the link between adult and pediatric cases, which were housed in different buildings and cared for by different HCWs. We focused on a potential role of ventilators in the outbreak based on the observation that (i) all NICU outbreak cases were on invasive or noninvasive ventilator support prior to culture positivity, (ii) the three adult patients were ventilated for at least part of their hospitalizations, and (iii) prior to identification of the NICU outbreak, ventilators were shared between adult and pediatric wards. Ventilator exchange between units was discontinued after the first NICU cases were identified.

The ventilators that were present in the NICU at the time environmental surveillance was performed tested negative for MRSA, but we could not rule out earlier contamination or contributions of other ventilators. Retrospective analysis of equipment usage logs and tracking data provided by the hospital’s real-time location system (RTLS) identified six units that were shared between outbreak cases (Fig. 3B, numbered 1 to 6). Ventilator 1 was briefly used by adult p64 and then transferred to several locations before it was moved into the NICU and later used by infant p150. The first NICU isolate that matched the outbreak clone by antibiogram was isolated from this patient soon after (Fig. 3A and B). Ventilator 4 was used by adult p91 several weeks before this patient developed bacteremia, except for a 2-day period when it was used by infant p151, shortly before the first NICU outbreak case (Fig. 3B). Infant p151 was in the neighboring PICU at this time and remained there until a positive surveillance isolate was obtained. Finally, ventilator 2 was used by adult p64 in two separate hospital visits, but was only moved to the NICU after the outbreak had already spread there.

Within the NICU, the sequential use of ventilator 6 by patients p90 and p133, the timing of their respective culture positivity, and the similarity of their isolate genomes all supported a role for this ventilator in the transmission to p133. Likewise, ventilators 2 and/or 5 may have been a factor in the spread from room 2 (subgroup C) to room 3 (subgroup D), especially considering that both rooms were cleaned just prior to the transmission (Fig. 3A). Ventilator 5 may also have been a transmission vector from p150 to p110. Ventilator 3 was used by p141 and later by p149; however, it is unclear if it played a role in the outbreak, as the first two isolates obtained from p141 after ventilator 3 exposure did not match the outbreak. Altogether, the epidemiological and genomic data suggest that ventilators not only played a role in spreading the outbreak from adult wards to the NICU but were also a factor in subsequent subtransmissions within the NICU.

Mutations in the outbreak clone alter expression of virulence and persistence factors.Given the extended duration of the outbreak, we next sought to identify genomic features that could have contributed to its persistence. A comparison of complete genomes found 42 nonsynonymous or deleterious SNVs and indels in the outbreak clone that were not present in any of the ST105 hospital background strains, affecting 35 genes or their promoter regions (Fig. 4A). The products of these genes were primarily involved in nucleotide, amino acid, and energy metabolism as well as environmental signal processing and drug resistance. Several genes encoding cell wall proteins were also affected, including gatD, which is involved in amidation of peptidoglycan (57). Pan-genome analysis with Roary (58) further revealed 71 genes exclusive to the outbreak strain or infrequently (<33%) present in other MLST105 isolates (Fig. 4B). Most of these genes were associated with three prophage regions and a 43.5-kbp plasmid. The additional genes in prophage A encoded only phage replication or hypothetical proteins. Among the genes in prophage B was an extra copy of clpB, which promotes stress tolerance, intracellular replication, and biofilm formation (59). Prophage C included an extra copy of the sep gene encoding an enterotoxin P-like protein associated with an increased risk of MRSA bacteremia in colonized patients (60). The 43.5-kbp plasmid contained the mupirocin (mupA) and gentamicin (aacA-aphD) resistance genes (see Fig. S6B) that explained the distinct susceptibility profile of the outbreak clone. High-level mupirocin resistance (HLR) conferred by mupA has been linked to transmissions in previous studies (61, 62). Pan-genome analysis also revealed a unique variant of the hsdS gene in the outbreak strain, which encodes the specificity subunit of a type I restriction modification (RM) system. S. aureus typically contains two type I RM systems that vary in sequence specificity based on the configuration of two target recognition domains (TRDs), which have been categorized according to the DNA motifs they recognize (63). Closer examination revealed that a recombination event in one of the two hsdS gene copies present in USA100 (see Fig. S7) changed its TRDs from the typical CC5 “B-D” configuration (recognizing the “AGG-5-GAT” motif present at 738 sites, overlapping 595 genes and 120 promoter regions) to an “A-D” configuration (recognizing the CCAY-5-GAT motif present at 304 sites, overlapping 287 genes and 15 promoter regions), resulting in altered genome-wide m6A DNA methylation profiles compared to those of other ST105 isolates (Fig. 4B).

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

Differentiating features of the NICU outbreak clone compared to the USA100 hospital background. (A) Map of nonsynonymous SNVs in genes and promoter regions that are unique to the outbreak clone. Gene identifiers or names are shown next to their genomic location. The SNV type is indicated by colors with a key shown at the top right. KEGG pathways with two or more genes are indicated on the right (green boxes) and corresponding gene descriptions on the far right. (B) Pan-genome analysis of MLST105 isolates showing all genes present in the outbreak clone and absent from at least half of the nonoutbreak isolates collected during our study. A maximum likelihood phylogenetic tree based on SNV distances in core genome alignments is shown on the left with patient (p) or environmental (e) isolate identifiers. Changes in the m6A methylation profile due to the hsdS recombination in the outbreak strain are highlighted in green/blue. Gene presence (yellow) or absence (red) is indicated in a matrix organized by genomic location (top). Gene names and descriptions are shown at the top and bottom of the matrix, respectively. See key on bottom left for more details. (C) Hierarchical clustering of 35 genes with significant expression differences (false-discovery rate [FDR] q < 0.05) between three control and three outbreak strains. Columns correspond to control or outbreak isolates, with labels at the top. Gene names and descriptions are shown on the right. Color shades and intensity represent the difference in normalized log2 counts per million (CPM) relative to the average gene expression level, with a color key shown below.

We reasoned that the (epi)genetic changes in the outbreak clone could alter gene expression patterns and provide further insights into the effects of these changes. We therefore compared the gene expression profiles of three representative outbreak isolates (i.e., cases) to the three most similar nonoutbreak ST105 strains (i.e., controls) during late-log-phase growth. The control strains shared the 43.5-kbp plasmid and most of the prophage elements with the outbreak strain and demonstrated similar growth characteristics (see Fig. S8). Differential gene expression analysis showed altered expression of 35 genes (Fig. 4C). Two of these genes were mutated in the outbreak clone: a SNP in the promoter region of sdhC and a duplication of clpB. Methylation changes were found in six genes (17.1%), which was lower than the rate of 27.3% across all genes. Thus, most expression changes appear to be indirect results of (epi)genetic changes. Multiple upregulated genes in the outbreak clone encoded proteins involved in stress and heat shock responses. This included clpB, which was increased in copy number in the outbreak versus control strains, but also dnaK and clpC, which have been linked to biofilm formation in S. aureus and adherence to eukaryotic cells (64, 65). Expression of the gene encoding staphylococcal superantigen-like protein 5 (SSL5) was also increased. SSL5 is known to inhibit leukocyte activation by chemokines and anaphylatoxins (66). Among the downregulated genes, the agrABC genes of the accessory gene regulator (agr) locus stood out. agr is the major virulence regulator in S. aureus (67), and decreased agr function in clinical isolates is associated with attenuated virulence and increased biofilm and surface protein expression (68).

DISCUSSION

In this study, we implemented a complete genome screening program at a large quaternary urban medical center, with the aim of tracking circulating clones, to identify transmission events and to understand the genomic epidemiology of endemic strains impacting human health. The availability of complete genomes allowed us to precisely map all genetic changes between strains, highlighting the presence of substantial structural variation in lineages that are considered highly clonal. The extent of variation due to recombinations in prophages, mobilization of genetic elements, and large genomic inversions also impacted classical spa, MLST, and signature virulence and resistance elements used in S. aureus molecular typing schemes. As such, the stability of these elements should be considered when using such schemes for lineage analysis. Complete reconstruction of outbreak genomes provided comprehensive variation data to map subtransmission events during a NICU MRSA outbreak. Finally, the combination of genetic and gene expression differences between the NICU outbreak clone and USA100 hospital background revealed genomic features that may have contributed to its persistence.

Complete genome analysis of the outbreak clone revealed a pattern of genetic changes that matched patient locations, suggesting that transmission bottlenecks and local environmental contamination led to a unique genetic signature at each site. Some isolates and isolate subgroups were separated by >10 variants, which is relatively high considering a reported core genome mutation rate of 2.7 to 3.3 mutations per Mb per year (14, 54). This suggests that the outbreak may have originated from a genetically heterogeneous source, such as a patient with a history of persistent MRSA colonization that accumulated intrahost variants. It is also possible that the combination of selection pressures and transmission bottlenecks contributed to the diversification of the outbreak clone. Considering all available data, we think the most likely scenario is that the NICU outbreak originated from patient p64 and then spread to other adult patients through direct or indirect contact in shared wards. Ventilator 1, used by adult p64 and infant p150, was the most likely vector for entry into the NICU. Ventilator 4 may have provided a potential second entry route via p151, with subsequent transmissions to p141 and p148 (p151 and p141 had an overlapping stay in the PICU). Such a secondary introduction may explain why the p141 and p148 isolates were more distantly related to all other NICU isolates. We were not able to confirm this scenario, as the isolates from p151 were no longer available. All subsequent cases could be explained by location relative to other MRSA colonized patients or sharing of MRSA-exposed ventilators.

The outbreak strain genome differed from the hospital background by multiple mutations of core genes as well as accessory gene gain and loss. Hundreds of genes were impacted by DNA methylation changes in the gene body or promoter regions, but such genes were depleted rather than enriched among differentially expressed genes. As such, the impact of the methylation changes on the outbreak clone (if any) was unclear. Nonetheless, a common theme among the genetic and expression changes was the relevance of genes involved in biofilm formation, persistence, and quorum sensing. Although the collective impact of the mutations will require further investigation, we speculate that these changes may have contributed an increased persistence of the outbreak clone in the environment.

Complete genome data from our hospital-wide screening program provided key information for outbreak management and investigation that could not have been obtained by molecular typing. First, it provided conclusive differentiation of outbreak from nonoutbreak isolates within 10 to 14 days while the outbreak was ongoing, which helped delineate the final case set, identify transmission events, and determine when the outbreak ended. Second, a retrospective analysis of all genetic differences between outbreak cases allowed us to identify subtransmissions and better understand the chain of events that led to each subtransmission. Third, a retrospective analysis integrating all hospital-wide genomic surveillance data indicated that the NICU outbreak had originated much earlier in unrelated adult wards in a different building and, together with electronic location tracking data, helped identify ventilators as likely transmission vectors.

There are some limitations to our study. Our genomic survey was limited to first positive single-patient bacteremias, and additional transmissions may have been missed by excluding nonblood isolates. Moreover, by sequencing single colony isolates, we likely did not fully capture intrahost heterogeneity. Although such heterogeneity may be less common among bacteremias, we did encounter variation within some patients, which was considered when establishing our transmission thresholds. Finally, while we believe that we reconstructed the most likely transmission routes and vectors, we could not definitively link ventilators to the outbreak because we did not detect MRSA contamination. This may be explained by the fact that some ventilators were no longer available for testing when environmental surveillance was initiated, while others were only tested after the transmission events they were implicated in. Nonetheless, it is possible that other factors such as spread by HCWs and/or other vectors contributed as well.

In conclusion, we find that the application of routine genome sequencing in the clinical space provides significant benefits for infection prevention and control. In addition to providing contemporary data on the genomic characteristics of circulating lineages, timely directed intervention and containment of identified transmission events can help prevent further outbreak progression. Although our screening program was limited in scope to bacteremias, early detection of a transmission event between the adult and NICU ward could conceivably have allowed staff to intervene earlier. Accumulating a larger repository of complete and unique genome references and variants associated with successful spreading strains may be key to future outbreak detection and prevention programs by providing high-resolution feature sets for prospective and retrospective data mining purposes.

ACKNOWLEDGMENTS

This research was supported in part by R01 AI119145 (to H.V.B.), the Icahn Institute for Genomics and Data Science (to A.K. and E.S.), the NIAID-supported NRSA Institutional Research Training Grant for Global Health Research (T32 AI07647 to D.R.A), the CTSA/NCATS KL2 Program (KL2TR001435, Icahn School of Medicine at Mount Sinai, to D.R.A), the New York State Department of Health Empire Clinical Research Investigator Program (awarded to Judith A. Aberg, Icahn School of Medicine at Mount Sinai, to D.R.A.), and F30 AI122673 (to T.R.P.). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. This research was also supported in part through the computational resources and staff expertise provided by the Department of Scientific Computing at the Icahn School of Medicine at Mount Sinai.

We thank Karen Southwick, Eleanor Adams, Lin Ying, and John Kornblum from the New York State Department of Health (NYSDOH) for their consultation during the outbreak and providing PFGE results for a subset of the outbreak isolates.

M.J.S, D.R.A., and H.V.B. contributed to the methodology; M.J.S, D.R.A., K.I.C., B.C., E.W., T.R.P., G.D., M.L.-S., Z.K., C.B., A.R., F.S., K.G., and H.V.B. collected the data; M.J.S, D.R.A., K.I.C., and B.C. curated the data; M.J.S, D.R.A., and K.I.C. performed analyses; M.J.S, D.R.A., K.I.C., and H.V.B. constructed the figures; M.J.S, D.R.A., and H.V.B. wrote the draft; all authors reviewed and edited the manuscript; M.J.S, D.R.A., H.V.B., and K.G. designed the study; H.V.B. and K.G. supervised the study and were the project administrators; H.V.B., D.R.A., T.R.P., A.K. and E.S. acquired the funding.

FOOTNOTES

    • Received 31 July 2019.
    • Returned for modification 20 August 2019.
    • Accepted 23 September 2019.
    • Accepted manuscript posted online 2 October 2019.
  • Supplemental material for this article may be found at https://doi.org/10.1128/JCM.01261-19.

  • Copyright © 2019 Sullivan 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. Klevens RM,
    2. Morrison MA,
    3. Nadle J,
    4. Petit S,
    5. Gershman K,
    6. Ray S,
    7. Harrison LH,
    8. Lynfield R,
    9. Dumyati G,
    10. Townes JM,
    11. Craig AS,
    12. Zell ER,
    13. Fosheim GE,
    14. McDougal LK,
    15. Carey RB,
    16. Fridkin SK
    , Active Bacterial Core surveillance (ABCs) MRSA Investigators. 2007. Invasive methicillin-resistant Staphylococcus aureus infections in the United States. JAMA 298:1763–1771. doi:10.1001/jama.298.15.1763.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Grundmann H,
    2. Aires-de-Sousa M,
    3. Boyce J,
    4. Tiemersma E
    . 2006. Emergence and resurgence of meticillin-resistant Staphylococcus aureus as a public-health threat. Lancet 368:874–885. doi:10.1016/S0140-6736(06)68853-3.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. DeLeo FR,
    2. Chambers HF
    . 2009. Reemergence of antibiotic-resistant Staphylococcus aureus in the genomics era. J Clin Invest 119:2464–2474. doi:10.1172/JCI38226.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    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 methicillin-resistant Staphylococcus aureus. Lancet 367:731–739. doi:10.1016/S0140-6736(06)68231-7.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. DeLeo FR,
    2. Otto M,
    3. Kreiswirth BN,
    4. Chambers HF
    . 2010. Community-associated methicillin-resistant Staphylococcus aureus. Lancet 375:1557–1568. doi:10.1016/S0140-6736(09)61999-1.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Glaser P,
    2. Martins-Simões P,
    3. Villain A,
    4. Barbier M,
    5. Tristan A,
    6. Bouchier C,
    7. Ma L,
    8. Bes M,
    9. Laurent F,
    10. Guillemot D,
    11. Wirth T,
    12. Vandenesch F
    . 2016. Demography and intercontinental spread of the USA300 community-acquired methicillin-resistant Staphylococcus aureus lineage. mBio 7:e02183. doi:10.1128/mBio.02183-15.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Montgomery CP,
    2. Boyle-Vavra S,
    3. Daum RS
    . 2009. The arginine catabolic mobile element is not associated with enhanced virulence in experimental invasive disease caused by the community-associated methicillin-resistant Staphylococcus aureus USA300 genetic background. Infect Immun 77:2650–2656. doi:10.1128/IAI.00256-09.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Uhlemann A-C,
    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 U S A 111:6738–6743. doi:10.1073/pnas.1401006111.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Planet PJ,
    2. Diaz L,
    3. Kolokotronis S-O,
    4. Narechania A,
    5. Reyes J,
    6. Xing G,
    7. Rincon S,
    8. Smith H,
    9. Panesso D,
    10. Ryan C,
    11. Smith DP,
    12. Guzman M,
    13. Zurita J,
    14. Sebra R,
    15. Deikus G,
    16. Nolan RL,
    17. Tenover FC,
    18. Weinstock GM,
    19. Robinson DA,
    20. Arias CA
    . 2015. Parallel epidemics of community-associated methicillin-resistant Staphylococcus aureus USA300 infection in North and South America. J Infect Dis 212:1874–1882. doi:10.1093/infdis/jiv320.
    OpenUrlCrossRefPubMed
  10. 10.↵
    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
  11. 11.↵
    1. Price J,
    2. Gordon NC,
    3. Crook D,
    4. Llewelyn M,
    5. Paul J
    . 2013. The usefulness of whole genome sequencing in the management of Staphylococcus aureus infections. Clin Microbiol Infect 19:784–789. doi:10.1111/1469-0691.12109.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Azarian T,
    2. Cook RL,
    3. Johnson JA,
    4. Guzman N,
    5. McCarter YS,
    6. Gomez N,
    7. Rathore MH,
    8. Morris JG,
    9. Salemi M
    . 2015. Whole-genome sequencing for outbreak investigations of methicillin-resistant Staphylococcus aureus in the neonatal intensive care unit: time for routine practice? Infect Control Hosp Epidemiol 36:777–785. doi:10.1017/ice.2015.73.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Altman DR,
    2. Sebra R,
    3. Hand J,
    4. Attie O,
    5. Deikus G,
    6. Carpini KWD,
    7. Patel G,
    8. Rana M,
    9. Arvelakis A,
    10. Grewal P,
    11. Dutta J,
    12. Rose H,
    13. Shopsin B,
    14. Daefler S,
    15. Schadt E,
    16. Kasarskis A,
    17. van Bakel H,
    18. Bashir A,
    19. Huprikar S
    , Others. 2014. Transmission of methicillin-resistant Staphylococcus aureus via deceased donor liver transplantation confirmed by whole genome sequencing. Am J Transplant 14:2640–2644. doi:10.1111/ajt.12897.
    OpenUrlCrossRefPubMed
  14. 14.↵
    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
  15. 15.↵
    1. Snitkin ES,
    2. Zelazny AM,
    3. Thomas PJ,
    4. Stock F, NISC Comparative Sequencing Program Group,
    5. Henderson DK,
    6. Palmore TN,
    7. Segre JA
    . 2012. Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing. Sci Transl Med 4:148ra116. doi:10.1126/scitranslmed.3004129.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Eyre DW,
    2. Golubchik T,
    3. Gordon NC,
    4. Bowden R,
    5. Piazza P,
    6. Batty EM,
    7. Ip CLC,
    8. Wilson DJ,
    9. Didelot X,
    10. O'Connor L,
    11. Lay R,
    12. Buck D,
    13. Kearns AM,
    14. Shaw A,
    15. Paul J,
    16. Wilcox MH,
    17. Donnelly PJ,
    18. Peto TEA,
    19. Walker AS,
    20. Crook DW
    . 2012. A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance. BMJ Open 2:e001124. doi:10.1136/bmjopen-2012-001124.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    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 U S A 104:9451–9456. doi:10.1073/pnas.0609839104.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Benson MA,
    2. Ohneck EA,
    3. Ryan C,
    4. Alonzo F, 3rd,
    5. Smith H,
    6. Narechania A,
    7. Kolokotronis S-O,
    8. Satola SW,
    9. Uhlemann A-C,
    10. Sebra R,
    11. Deikus G,
    12. Shopsin B,
    13. Planet PJ,
    14. Torres VJ
    . 2014. Evolution of hypervirulence by a MRSA clone through acquisition of a transposable element. Mol Microbiol 93:664–681. doi:10.1111/mmi.12682.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Senn L,
    2. Clerc O,
    3. Zanetti G,
    4. Basset P,
    5. Prod'hom G,
    6. Gordon NC,
    7. Sheppard AE,
    8. Crook DW,
    9. James R,
    10. Thorpe HA,
    11. Feil EJ,
    12. Blanc DS
    . 2016. The stealthy superbug: the role of asymptomatic enteric carriage in maintaining a long-term hospital outbreak of ST228 methicillin-resistant Staphylococcus aureus. mBio 7:e02039. doi:10.1128/mBio.02039-15.
    OpenUrlCrossRef
  20. 20.↵
    1. Lindsay JA,
    2. Holden MT
    . 2004. Staphylococcus aureus: superbug, super genome? Trends Microbiol 12:378–385. doi:10.1016/j.tim.2004.06.004.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Sela U,
    2. Euler CW,
    3. Correa da Rosa J,
    4. Fischetti VA
    . 2018. Strains of bacterial species induce a greatly varied acute adaptive immune response: the contribution of the accessory genome. PLoS Pathog 14:e1006726. doi:10.1371/journal.ppat.1006726.
    OpenUrlCrossRef
  22. 22.↵
    1. Chin CS,
    2. Alexander DH,
    3. Marks P,
    4. Klammer AA,
    5. Drake J,
    6. Heiner C,
    7. Clum A,
    8. Copeland A,
    9. Huddleston J,
    10. Eichler EE,
    11. Turner SW,
    12. Korlach J
    . 2013. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 10:563–569. doi:10.1038/nmeth.2474.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Madoui M-A,
    2. Engelen S,
    3. Cruaud C,
    4. Belser C,
    5. Bertrand L,
    6. Alberti A,
    7. Lemainque A,
    8. Wincker P,
    9. Aury J-M
    . 2015. Genome assembly using Nanopore-guided long and error-free DNA reads. BMC Genomics 16:327. doi:10.1186/s12864-015-1519-z.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Fang G,
    2. Munera D,
    3. Friedman DI,
    4. Mandlik A,
    5. Chao MC,
    6. Banerjee O,
    7. Feng Z,
    8. Losic B,
    9. Mahajan MC,
    10. Jabado OJ,
    11. Deikus G,
    12. Clark TA,
    13. Luong K,
    14. Murray IA,
    15. Davis BM,
    16. Keren-Paz A,
    17. Chess A,
    18. Roberts RJ,
    19. Korlach J,
    20. Turner SW,
    21. Kumar V,
    22. Waldor MK,
    23. Schadt EE
    . 2012. Genome-wide mapping of methylated adenine residues in pathogenic Escherichia coli using single-molecule real-time sequencing. Nat Biotechnol 30:1232–1239. doi:10.1038/nbt.2432.
    OpenUrlCrossRefPubMed
  25. 25.↵
    CLSI. 2015. Performance standards for antimicrobial susceptibility testing; 25th informational supplement. CLSI Document M100-S25. Clinical and Laboratory Standards Institute, Wayne PA.
  26. 26.↵
    1. Altman DR,
    2. Sullivan MJ,
    3. Chacko KI,
    4. Balasubramanian D,
    5. Pak TR,
    6. Sause WE,
    7. Kumar K,
    8. Sebra R,
    9. Deikus G,
    10. Attie O,
    11. Rose H,
    12. Lewis M,
    13. Fulmer Y,
    14. Bashir A,
    15. Kasarskis A,
    16. Schadt EE,
    17. Richardson AR,
    18. Torres VJ,
    19. Shopsin B,
    20. van Bakel H
    . 2018. Genome plasticity of agr-defective Staphylococcus aureus during clinical infection. Infect Immun 86:e00331-18. doi:10.1128/IAI.00331-18.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Hunt M,
    2. Silva ND,
    3. Otto TD,
    4. Parkhill J,
    5. Keane JA,
    6. Harris SR
    . 2015. Circlator: automated circularization of genome assemblies using long sequencing reads. Genome Biol 16:294. doi:10.1186/s13059-015-0849-0.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Camacho C,
    2. Coulouris G,
    3. Avagyan V,
    4. Ma N,
    5. Papadopoulos J,
    6. Bealer K,
    7. Madden TL
    . 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi:10.1186/1471-2105-10-421.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Sullivan MJ,
    2. Ben Zakour NL,
    3. Forde BM,
    4. Stanton-Cook M,
    5. Beatson SA
    . 2015. Contiguity: contig adjacency graph construction and visualisation. PeerJ PrePrints 3:e1037v1. doi:10.7287/peerj.preprints.1037v1.
    OpenUrlCrossRef
  30. 30.↵
    1. Seemann T
    . 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069. doi:10.1093/bioinformatics/btu153.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Pak TR,
    2. Roth FP
    . 2013. ChromoZoom: a flexible, fluid, web-based genome browser. Bioinformatics 29:384–386. doi:10.1093/bioinformatics/bts695.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Nicol JW,
    2. Helt GA,
    3. Blanchard SG, Jr,
    4. Raja A,
    5. Loraine AE
    . 2009. The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics 25:2730–2731. doi:10.1093/bioinformatics/btp472.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Jones P,
    2. Binns D,
    3. Chang H-Y,
    4. Fraser M,
    5. Li W,
    6. McAnulla C,
    7. McWilliam H,
    8. Maslen J,
    9. Mitchell A,
    10. Nuka G,
    11. Pesseat S,
    12. Quinn AF,
    13. Sangrador-Vegas A,
    14. Scheremetjew M,
    15. Yong S-Y,
    16. Lopez R,
    17. Hunter S
    . 2014. InterProScan 5: genome-scale protein function classification. Bioinformatics 30:1236–1240. doi:10.1093/bioinformatics/btu031.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Li H
    . 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 1303.3997 [q-bioGN] https://arxiv.org/abs/1303.3997.
  35. 35.↵
    1. Garrison E,
    2. Marth G
    . 2012. Haplotype-based variant detection from short-read sequencing. arXiv 1207.3907 [q-bioGN] https://arxiv.org/abs/1207.3907.
  36. 36.↵
    1. Bachmann NL,
    2. Sullivan MJ,
    3. Jelocnik M,
    4. Myers GSA,
    5. Timms P,
    6. Polkinghorne A
    . 2015. Culture-independent genome sequencing of clinical samples reveals an unexpected heterogeneity of infections by Chlamydia pecorum. J Clin Microbiol 53:1573–1581. doi:10.1128/JCM.03534-14.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Treangen TJ,
    2. Ondov BD,
    3. Koren S,
    4. Phillippy AM
    . 2014. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol 15:524. doi:10.1186/s13059-014-0524-x.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Minkin I,
    2. Patel A,
    3. Kolmogorov M,
    4. Vyahhi N,
    5. Pham S
    . 2013. Sibelia: a scalable and comprehensive synteny block generation tool for closely related microbial genomes, p 215–229. Algorithms in bioinformatics. Springer, Berlin, Germany.
  39. 39.↵
    1. Khelik K,
    2. Lagesen K,
    3. Sandve GK,
    4. Rognes T,
    5. Nederbragt AJ
    . 2017. NucDiff: in-depth characterization and annotation of differences between two sets of DNA sequences. BMC Bioinformatics 18:338. doi:10.1186/s12859-017-1748-z.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Jolley KA,
    2. Bray JE,
    3. Maiden MCJ
    . 2017. A RESTful application programming interface for the PubMLST molecular typing and genome databases. Database (Oxford) 2017:bax060. doi:10.1093/database/bax060.
    OpenUrlCrossRef
  41. 41.↵
    1. Kaya H,
    2. Hasman H,
    3. Larsen J,
    4. Stegger M,
    5. Johannesen TB,
    6. Allesøe RL,
    7. Lemvigh CK,
    8. Aarestrup FM,
    9. Lund O,
    10. Larsen AR
    . 2018. SCCmecFinder, a web-based tool for typing of staphylococcal cassette chromosome mec in Staphylococcus aureus using whole-genome sequence data. mSphere 3:e00612-17. doi:10.1128/mSphere.00612-17.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    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
  43. 43.↵
    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
  44. 44.↵
    1. Kanehisa M,
    2. Sato Y,
    3. Morishima K
    . 2016. BlastKOALA and GhostKOALA: KEGG Tools for functional characterization of genome and metagenome sequences. J Mol Biol 428:726–731. doi:10.1016/j.jmb.2015.11.006.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Huerta-Cepas J,
    2. Serra F,
    3. Bork P
    2016. ETE 3: reconstruction, analysis, and visualization of phylogenomic data. Mol Biol Evol 33:1635–1638. doi:10.1093/molbev/msw046.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Flusberg BA,
    2. Webster DR,
    3. Lee JH,
    4. Travers KJ,
    5. Olivares EC,
    6. Clark TA,
    7. Korlach J,
    8. Turner SW
    . 2010. Direct detection of DNA methylation during single-molecule, real-time sequencing. Nat Methods 7:461–465. doi:10.1038/nmeth.1459.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Martin M
    . 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 17:10–12. doi:10.14806/ej.17.1.200.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357–359. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    . 2015. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 31:166–169. doi:10.1093/bioinformatics/btu638.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Robinson MD,
    2. McCarthy DJ,
    3. Smyth GK
    . 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139–140. doi:10.1093/bioinformatics/btp616.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Robinson MD,
    2. Oshlack A
    . 2010. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol 11:R25. doi:10.1186/gb-2010-11-3-r25.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Pardos de la Gandara M,
    2. Curry M,
    3. Berger J,
    4. Burstein D,
    5. Della-Latta P,
    6. Kopetz V,
    7. Quale J,
    8. Spitzer E,
    9. Tan R,
    10. Urban C,
    11. Wang G,
    12. Whittier S,
    13. de Lencastre H,
    14. Tomasz A
    . 2016. MRSA causing infections in hospitals in greater metropolitan New York: major shift in the dominant clonal type between 1996 and 2014. PLoS One 11:e0156924. doi:10.1371/journal.pone.0156924.
    OpenUrlCrossRef
  53. 53.↵
    1. Copin R,
    2. Shopsin B,
    3. Torres VJ
    . 2018. After the deluge: mining Staphylococcus aureus genomic data for clinical associations and host-pathogen interactions. Curr Opin Microbiol 41:43–50. doi:10.1016/j.mib.2017.11.014.
    OpenUrlCrossRef
  54. 54.↵
    1. Young BC,
    2. Golubchik T,
    3. Batty EM,
    4. Fung R,
    5. Larner-Svensson H,
    6. Votintseva AA,
    7. Miller RR,
    8. Godwin H,
    9. Knox K,
    10. Everitt RG,
    11. Iqbal Z,
    12. Rimmer AJ,
    13. Cule M,
    14. Ip CLC,
    15. Didelot X,
    16. Harding RM,
    17. Donnelly P,
    18. Peto TE,
    19. Crook DW,
    20. Bowden R,
    21. Wilson DJ
    . 2012. Evolutionary dynamics of Staphylococcus aureus during progression from carriage to disease. Proc Natl Acad Sci U S A 109:4550–4555. doi:10.1073/pnas.1113219109.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Von Eiff C,
    2. Becker K,
    3. Machka K,
    4. Stammer H,
    5. Peters G
    . 2001. Nasal carriage as a source of Staphylococcus aureus bacteremia. N Engl J Med 344:11–16. doi:10.1056/NEJM200101043440102.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Gibbs K,
    2. DeMaria S,
    3. McKinsey S,
    4. Fede A,
    5. Harrington A,
    6. Hutchison D,
    7. Torchen C,
    8. Levine A,
    9. Goldberg A
    . 2018. A novel in situ simulation intervention used to mitigate an outbreak of methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit. J Pediatr 194:22.e5–27.e5. doi:10.1016/j.jpeds.2017.10.040.
    OpenUrlCrossRef
  57. 57.↵
    1. Figueiredo TA,
    2. Sobral RG,
    3. Ludovice AM,
    4. de Almeida JMF,
    5. Bui NK,
    6. Vollmer W,
    7. de Lencastre H,
    8. Tomasz A
    . 2012. Identification of genetic determinants and enzymes involved with the amidation of glutamic acid residues in the peptidoglycan of Staphylococcus aureus. PLoS Pathog 8:e1002508. doi:10.1371/journal.ppat.1002508.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Page AJ,
    2. Cummins CA,
    3. Hunt M,
    4. Wong VK,
    5. Reuter S,
    6. Holden MTG,
    7. Fookes M,
    8. Falush D,
    9. Keane JA,
    10. Parkhill J
    . 2015. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31:3691–3693. doi:10.1093/bioinformatics/btv421.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Frees D,
    2. Chastanet A,
    3. Qazi S,
    4. Sørensen K,
    5. Hill P,
    6. Msadek T,
    7. Ingmer H
    . 2004. Clp ATPases are required for stress tolerance, intracellular replication and biofilm formation in Staphylococcus aureus. Mol Microbiol 54:1445–1462. doi:10.1111/j.1365-2958.2004.04368.x.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Calderwood MS,
    2. Desjardins CA,
    3. Sakoulas G,
    4. Nicol R,
    5. Dubois A,
    6. Delaney ML,
    7. Kleinman K,
    8. Cosimi LA,
    9. Feldgarden M,
    10. Onderdonk AB,
    11. Birren BW,
    12. Platt R,
    13. Huang SS
    , CDC Prevention Epicenters Program. 2014. Staphylococcal enterotoxin P predicts bacteremia in hospitalized patients colonized with methicillin-resistant Staphylococcus aureus. J Infect Dis 209:571–577. doi:10.1093/infdis/jit501.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Hodgson JE,
    2. Curnock SP,
    3. Dyke KG,
    4. Morris R,
    5. Sylvester DR,
    6. Gross MS
    . 1994. Molecular characterization of the gene encoding high-level mupirocin resistance in Staphylococcus aureus J2870. Antimicrob Agents Chemother 38:1205–1208. doi:10.1128/aac.38.5.1205.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Udo EE,
    2. Jacob LE,
    3. Mathew B
    . 2001. Genetic analysis of methicillin-resistant Staphylococcus aureus expressing high- and low-level mupirocin resistance. J Med Microbiol 50:909–915. doi:10.1099/0022-1317-50-10-909.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    1. Cooper LP,
    2. Roberts GA,
    3. White JH,
    4. Luyten YA,
    5. Bower EKM,
    6. Morgan RD,
    7. Roberts RJ,
    8. Lindsay JA,
    9. Dryden D
    . 2017. DNA target recognition domains in the type I restriction and modification systems of Staphylococcus aureus. Nucleic Acids Res 45:3395–3406. doi:10.1093/nar/gkx067.
    OpenUrlCrossRef
  64. 64.↵
    1. Singh VK,
    2. Syring M,
    3. Singh A,
    4. Singhal K,
    5. Dalecki A,
    6. Johansson T
    . 2012. An insight into the significance of the DnaK heat shock system in Staphylococcus aureus. Int J Med Microbiol 302:242–252. doi:10.1016/j.ijmm.2012.05.001.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Chatterjee I,
    2. Becker P,
    3. Grundmeier M,
    4. Bischoff M,
    5. Somerville GA,
    6. Peters G,
    7. Sinha B,
    8. Harraghy N,
    9. Proctor RA,
    10. Herrmann M
    . 2005. Staphylococcus aureus ClpC is required for stress resistance, aconitase activity, growth recovery, and death. J Bacteriol 187:4488–4496. doi:10.1128/JB.187.13.4488-4496.2005.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Bestebroer J,
    2. van Kessel KPM,
    3. Azouagh H,
    4. Walenkamp AM,
    5. Boer IGJ,
    6. Romijn RA,
    7. van Strijp JAG,
    8. de Haas C
    . 2009. Staphylococcal SSL5 inhibits leukocyte activation by chemokines and anaphylatoxins. Blood 113:328–337. doi:10.1182/blood-2008-04-153882.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Novick RP
    . 2003. Autoinduction and signal transduction in the regulation of staphylococcal virulence. Mol Microbiol 48:1429–1449. doi:10.1046/j.1365-2958.2003.03526.x.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Shopsin B,
    2. Copin R
    . 2018. Staphylococcus aureus adaptation during infection, p 431–459. In Fong IW, Shlaes D, Drlica K (ed), Antimicrobial resistance in the 21st century. Springer International Publishing, Cham, Switzerland.
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak
Mitchell J. Sullivan, Deena R. Altman, Kieran I. Chacko, Brianne Ciferri, Elizabeth Webster, Theodore R. Pak, Gintaras Deikus, Martha Lewis-Sandari, Zenab Khan, Colleen Beckford, Angela Rendo, Flora Samaroo, Robert Sebra, Ramona Karam-Howlin, Tanis Dingle, Camille Hamula, Ali Bashir, Eric Schadt, Gopi Patel, Frances Wallach, Andrew Kasarskis, Kathleen Gibbs, Harm van Bakel
Journal of Clinical Microbiology Nov 2019, 57 (12) e01261-19; DOI: 10.1128/JCM.01261-19

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 Journal of Clinical Microbiology 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.
A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak
(Your Name) has forwarded a page to you from Journal of Clinical Microbiology
(Your Name) thought you would be interested in this article in Journal of Clinical Microbiology.
Share
A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak
Mitchell J. Sullivan, Deena R. Altman, Kieran I. Chacko, Brianne Ciferri, Elizabeth Webster, Theodore R. Pak, Gintaras Deikus, Martha Lewis-Sandari, Zenab Khan, Colleen Beckford, Angela Rendo, Flora Samaroo, Robert Sebra, Ramona Karam-Howlin, Tanis Dingle, Camille Hamula, Ali Bashir, Eric Schadt, Gopi Patel, Frances Wallach, Andrew Kasarskis, Kathleen Gibbs, Harm van Bakel
Journal of Clinical Microbiology Nov 2019, 57 (12) e01261-19; DOI: 10.1128/JCM.01261-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

MRSA
NICU outbreak
genome analysis

Related Articles

Cited By...

About

  • About JCM
  • Editor in Chief
  • Board of Editors
  • Editor Conflicts of Interest
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Resources for Clinical Microbiologists
  • Ethics
  • Contact Us

Follow #JClinMicro

@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

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

Print ISSN: 0095-1137; Online ISSN: 1098-660X