Genomic epidemiologic assessment implicates prolonged silent carriage, virulence factors and transmission between staff and patients in a NICU outbreak of MRSA
===============================================================================================================================================================

* Sharline Madera
* Nicole McNeil
* Paula Hayakawa Serpa
* Jack Kamm
* Christy Pak
* Carolyn Caughell
* Amy Nichols
* David Dynerman
* Lucy M. Li
* Estella Sanchez-Guerrero
* Maira Phelps
* Angela M. Detweiler
* Norma Neff
* Helen Reyes
* Steve Miller
* Deborah Yokoe
* Joseph L. DeRisi
* Lynn Ramirez-Avila
* Charles R. Langelier

## ABSTRACT

**Introduction** Methicillin-resistant *Staphylococcus aureus* (MRSA) is an important pathogen in neonatal intensive care units (NICU) that carries significant morbidity and mortality. Improving our understanding of MRSA transmission dynamics, especially among high risk patients, is an infection prevention priority.

**Methods** We investigated a cluster of clinical MRSA cases in the NICU using a combination of epidemiologic review and whole genome sequencing (WGS) of isolates from clinical and surveillance cultures obtained from patients and healthcare personnel (HCP).

**Results** Phylogenetic analysis identified two genetically distinct phylogenetic clades and revealed multiple silent transmission events between HCP and infants. The predominant outbreak strain harbored multiple virulence factors. Epidemiologic investigation and genomic analysis identified a HCP colonized with the dominant MRSA outbreak strain who cared for the majority of NICU patients who were infected or colonized with the same strain, including one NICU patient with severe infection seven months before the described outbreak. These results guided implementation of infection prevention interventions that prevented further transmission events.

**Conclusion** Silent transmission of MRSA between HCP and NICU patients likely contributed to a NICU outbreak involving a virulent MRSA strain. WGS enabled data-driven decision making to inform implementation of infection control policies that mitigated the outbreak. Prospective WGS coupled with epidemiologic analysis can be used to detect transmission events and prompt early implementation of control strategies.

**Summary** Whole genome sequencing identified a NICU outbreak of MRSA preceded by silent transmission between health care personnel and patients seven months prior to outbreak detection. This study suggests prospective genomic analysis may identify isolates with high virulence potential, offering an early opportunity for intervention.

Keywords
*   MRSA
*   Staphylococcus aureus
*   neonatal intensive care unit
*   hospital acquired infections
*   health care personnel
*   infection control
*   outbreak
*   antimicrobial resistance
*   whole genome sequencing

## INTRODUCTION

*Staphylococcus aureus* (*S. aureus)* remains a leading cause of hospital-acquired infections in the neonatal intensive care unit (NICU)[1]. Critically-ill and preterm neonates are particularly vulnerable to invasive infections with methicillin-resistant *S. aureus* (MRSA), which confers significant morbidity, mortality [2,3] and financial costs[3]. Despite this, the risk factors and mechanisms of transmission of MRSA in NICUs are incompletely understood.

While MRSA colonization has been identified as a risk factor for the development of invasive infection[4], the specific factors contributing to transmission of MRSA between NICU patients are not well defined. Factors previously identified include environmental reservoirs, as outbreak-associated MRSA strains can exhibit longer environmental persistence than sporadic-MRSA isolates[5]. Parental colonization has been described as a source of MRSA both via vertical[6] and horizontal[7,8] transmission; however, these transmission events are usually self-limited and have not been implicated in outbreaks. Exposure to colonized health-care personnel (HCP) or patients may play a role in precipitating outbreaks. In many hospitals, patients are tested for MRSA nares colonization upon admission, however HCP generally are not [9–11].

Here we report an outbreak investigation of MRSA in a NICU in which HCP were identified as a potential key link in transmission and colonization. Whole genome sequencing (WGS) further identified circulation of a highly virulent outbreak strain for at least 7 months prior to outbreak recognition in the NICU. We highlight an opportunity for prospective WGS in the surveillance setting to aid in the early identification of predominant circulating strains whose virulence may facilitate spread and invasive infection.

## METHODS

### Setting

The University of California San Francisco (UCSF) NICU contains 58 beds and is located within the 183-bed UCSF Benioff Children’s Hospital.

### Ethics Statement

This study was approved by the UCSF Institutional Review Board (IRB) under protocol 17-24056, which permitted WGS analysis of clinical isolates and review of clinical microbiology results and other electronic health record data.

### Microbial culture

Bacterial cultures were inoculated into BACTEC blood culture bottles (blood). MRSA surveillance cultures were performed using selective chromogenic agar (Thermo Scientific Spectra MRSA).

### Metagenomic sequencing

DNA was extracted from cultured isolates using the Zymo ZR Fungal/Bacterial DNA MiniPrep Kit, 100 ng of DNA from each sample was then sheared with NEB fragmentase and used to construct sequencing libraries using the NEBNext Ultra II Library Prep Kit (New England Biolabs). Adaptor ligated samples underwent 12 cycles of amplification with dual-indexing primers using the NEBNext Ultra II Q5 polymerase. Libraries were then quantified with Qubit and quality checked with a Bioanalyzer High Sensitivity DNA chip. Samples were pooled and sequenced on an Illumina NextSeq instrument using a NextSeq 500/500 Mid Output kit v2.5 (300 cycles).

### Bioinformatic analyses

Raw sequences were analyzed using the SNP Pipeline for Infectious Disease (SPID)[12]. SPID aligned samples against S. aureus reference genome USA300 TCH1516 using minimap2, followed by samtools to perform an mpileup. Julia code was then run to call a consensus allele at each position, and the SNP instances were computed between every pair of samples. Phylogenetic analysis was performed using randomized axelerated maximum likelihood (RAxML)[13]. Phylogenetic trees were further visualized with ETE Python API[14].

ARGannot was used to identify antimicrobial resistance genes from quality-filtered raw sequence data for each of the evaluated isolates[15]. Genes with greater than 99% gene coverage were included in the analyses. SCC mec elements were detected using SCC mec finder[16]. Multilocus sequence type (MLST) typing was performed using MLST 2.0[17]. Detection of *S. aureus* virulence genes was performed using VirulenceFinder[18].

## RESULTS

### Outbreak Description

After four infants in the NICU with invasive MRSA infections were identified in an eight-day period, an investigation was initiated. Retrospective review of microbiology data identified two additional positive MRSA cultures from infants in the same geographic zone of the NICU two months prior (Supplemental Table 1). NICU-wide surveillance screening of all infants and staff identified 16 other cases of MRSA colonization in asymptomatic infants or HCP (Supplemental Figure 1). Review of prior MRSA isolates identified an isolate from an infant with extensive MRSA infection 7 months prior (ID 10B, Supplemental Table 1). In total, 23 MRSA isolates from 18 patients and 5 HCP were evaluated, and baseline characteristics of each subject are listed in Supplemental Table 1. Antimicrobial susceptibility testing (AST) demonstrated similarities across some isolates, suggestive of an outbreak (Supplemental Figure 1).

### Genomic Analyses

All isolates underwent WGS. Assembly by alignment against the chromosomal sequence of the National Center for Biotechnology Information (NCBI) reference strain *S. aureus* USA300 TCH1516 confirmed an outbreak composed of two genetically distinct clades (Figure 1A). Pairwise single nucleotide polymorphism (SNP) was used as a measure of genetic relatedness between isolates. A cutoff of 23 SNPs was used to define isolate belonging to the same *S. aureus* phylogenetic clade[19]. Phylogenetic analysis revealed and outbreak composed of two distinct phylogenetic clades; clade 1 was composed of 10 infants and one HCP (11), and clade 2 composed of three infants and 1 HCP (15, Figure 1A). 8 of the 23 cases were deemed unrelated to either of the two clades based on SNP distance.

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/07/29/2021.07.28.21261307/F1/graphic-1.medium.gif)

[](http://medrxiv.org/content/early/2021/07/29/2021.07.28.21261307/F1/graphic-1)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/07/29/2021.07.28.21261307/F1/graphic-2.medium.gif)

[](http://medrxiv.org/content/early/2021/07/29/2021.07.28.21261307/F1/graphic-2)

Figure 1. Phylogenetic Analysis and Phenotypic Antibiogram of all Sequenced MRSA Cases.
Isolates are identified by case number followed by H for healthcare personnel, B for bacteremia, S for skin and soft tissue infection, and C for colonization. (A) Cases associated with clade 1 outbreak strain labeled red, and clade 2 outbreak strain labeled green. Reference *S. aureus* strains ECT-R2, N315, JH1, 04-02981, and USA 300 strains TCH1516, FPR3757, and ISMMS1, are shown. (B) Antimicrobial susceptibility pattern for all cases sequenced. Susceptibility to antibiotics are denoted as resistant (R) or intermediate (I). Sequence type (MLST) and Staphylococcal Cassette Chromosome (SCC) mec type are shown.

Clade 1 exhibited marked genomic similarity between the 11 isolates that make up the clade with a median SNP difference of 3, range 0-6. Similarly, four isolates comprising clade 2 shared striking genetic similarity with a mean SNP difference of 7, range 0-23. An isolate belonging to case ID 23EC was approximately 800 SNPs away from clade 2, and thus not likely part of a singular transmission event. Interestingly, clade 1 samples were more closely related to ST8 MRSA reference strains USA300 TCH 1516, ISMMS1, and FPR3757, with approximately 140 SNP differences; while clade 2 was more closely related to ST5 MRSA reference strain N315 with an approximately 700 SNP difference. MRSA reference strain N315 harbored a 17,000 SNP difference from clade 1, and ST8 MRSA reference strains carried roughly a 20,000 SNP difference from clade 2.

Confirmation of a highly related cluster of MRSA isolates through WGS prompted further epidemiologic investigation that revealed that eight out of the eleven infants associated outbreak clade 1 had been cared for by HCP 11 (Figure 2). HCP 11 also cared for patient 10B, who was identified 7 months prior to clade 1, raising the possibility that HCP 11 could be the source of the current outbreak. Almost all patients associated with the outbreak strains were born prematurely and exhibited low birthweight (Supplemental Table 1). None of the cases associated with outbreak clade 2 group exhibited active infection.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/07/29/2021.07.28.21261307/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/07/29/2021.07.28.21261307/F2)

Figure 2. Epidemiologic Tracing of MRSA Outbreak in the NICU.
Index case and 22 other cases associated with an MRSA outbreak. Healthcare personnel (HCP) are indicated by an H. Epidemiologic link to each of two HCP iis ndicated by yellow (11H) and purple (15H), respectively.

*In silico* MLST typing revealed that outbreak clade 1 was sequence type (ST)-8, and smaller outbreak clade 2 was ST-5 (Figure 1B). Staphylococcal cassette chromosome (SCC) *mec* typing revealed further similarities among outbreak clades, with clade 1 samples identified to harbor SCC*mec* IVa, and clade 2 cases identified to harbor SCC*mec* V (Figure 1B). Notably, MRSA isolate belonging to case ID 19C shared identical AST, MLST and SCCmec typing as clade 1 isolates but was approximately 270 SNPs away from clade 1, and therefore distinct from the predominant outbreak clade.

We next carried out antimicrobial-resistance, virulence and toxin gene assessment. Cases associated with both outbreak clade 1 and clade 2 exhibited similar antimicrobial-resistance gene profiles (Figure 3) and concordant phenotypic antimicrobial-resistance profile (Figure 1B). Interestingly, across non-outbreak associated isolates, isolates 17C and 18C, were found to lack the *mecA* gene, associated with resistance to beta-lactam antibiotics like methicillin and penicillin, despite being phenotypically resistant. Further, 18C also lacked the *blaZ* gene, associated with resistance to penicillin, despite being phenotypically resistant. Harboring of aminoglycoside resistance genes, *sat4A* or *spc*, was not associated with phenotypic antimicrobial-resistance.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/07/29/2021.07.28.21261307/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/07/29/2021.07.28.21261307/F3)

Figure 3. Antimicrobial resistance and toxin genes identified by WGS of MRSA isolates.
Genes associated with antibiotic resistance and toxins are shown for all cases sequenced. Gene presence is depicted in red and gene absence is depicted in blue. PCN denotes penicillin, MET methicillin, ERY erythromycin, CLIN clindamycin, CIP ciprofloxacin, DOX doxycycline, GEN gentamycin, and TMP trimethoprim. Genes associated with toxins are shown.

Presence of genes encoding virulence factors or toxins revealed similarities across clade 1 and clade 2 isolates. Interestingly, clade 1 isolates harbored *lukS/F*, which encodes the cytotoxin Panton-Valentine leukocidin (PVL), and the arginine catabolic mobile element (ACME), which are associated with severe invasive infection**[20]** or the promotion of persistence and skin-colonization**[21]**, respectively. Clade 1 was also found to harbor superantigen enterotoxin genes *seq* and *sek***[22]**. In further contrast to clade 1, clade 2 isolates contained the enterotoxin gene clade *sei, sem, sen, seo, seg*, and *seu* encoded by superantigen enterotoxin gene clade *egc***[23,24]**. All isolates were found to contain the gamma-hemolysin genes *hlgA/B/C* and the leukocidin genes *lukD/E*.

## DISCUSSION

Real-time deployment of WGS provided high resolution assessment of an outbreak involving dynamic exchange of genetically distinct MRSA strains between HCPs and patients in a NICU. Phylogenetic analysis uncovered two outbreak clades in the NICU, both involving HCP who had contact with most, if not all, clade associated isolates. Clade 1 was associated with invasive and severe infections, while clade 2 was associated with asymptomatic carriage.

Surprisingly, outbreak clade 1 included patient 10B, who was admitted to the NICU 7 months prior to the outbreak and was cared for by HCP 11, also found to be colonized by clade 1. This raises the possibility that HCP 11 may have played a role in the dissemination of the outbreak clade 1. Not all clade 1 isolates had an epidemiologic link to HCP 11, however, which raises the possibility of unrecognized inter-HCP transmission of the clade 1 strain or, alternatively, the presence of an environmental reservoir.

Assessment of *S. aureus* virulence factors in the strains from each clade suggested potential molecular determinants driving outbreak emergence. Clade 1 isolates associated with the pandemic clonal lineage USA300, classically regarded as a community-associated (CA-) MRSA strain. Clade 2 isolates most closely associated with strain NS315, a hospital-associated (HA-) MRSA strain. Both outbreak isolates harbored the lukS/F gene that encodes for the cytotoxin Panton-Valentine leukocidin (PVL), a protein first associated with epidemic CA-MRSA strains that are more likely to cause sepsis, necrotizing pneumonia, and necrotizing fasciitis[20]. However, unlike clade 2, the clade 1 strain harbored the genomic island referred to as ACME. First described in CA-MRSA strain USA300 in 2006, ACME likely originated from *S. epidermidis*[25]. This element is composed of 33 putative genes and two operons, *arc* and *opp*.

The arc operon encodes genes thought to be important in arginine catabolism and have been shown to confer increased survival in the acidic skin environment[26]. The AMCE *speG* gene, which encodes a spermidine acetyltransferase, has been shown to facilitate survival in the presence of lethal polyamines, spermidine and spermine that are produced by human skin[26]. Therefore, ACME-containing *S. aureus* are poised to execute epidemic, severe infections, likely explaining why clade 1 isolate was able to circulate in the NICU for over nine months, causing multiple episodes of severe infection along the way. Rapid identification of potentially high-virulence strains highlights an opportunity for preventative surveillance of this high-risk population and the HCP who provide their care.

Notably, isolates 16P and 20C also harbored the *lukS/F* gene and the ACME element. These isolates were not associated with a clade, and likely reflect singular acquisitions of highly virulent strains. Clade 1 and clade 2 isolates, on the other hand, had HCP to provide an epidemiologic link that could suggest repeated exposures to the highly virulent isolates that led to the pathologic infection and/or colonization of numerous neonates. This highlights the possibility that silent transmission of antimicrobial resistant organisms is taking place throughout the healthcare system, leaving our most vulnerable patients at greatest risk.

Genomic analysis also provided insight regarding antimicrobial resistance mechanisms of the evaluated isolates and revealed an unexpected discordance between antimicrobial resistance genes and phenotypic antimicrobial resistance. Resistance to methicillin and other β-lactam antibiotics in *S. aureus* is largely mediated by production of an altered penicillin binding protein (PBP)2a, which is encoded by *mecA* and is harbored in the large genetic element designated SCC*mec*[27]. Presence of *mecA* is considered hallmark in the identification of MRSA. Interestingly, infant isolates 17C and 18C did not harbor *mecA* (Figure 3), despite still exhibiting resistance to methicillin (Figure 1B). Borderline and low-level methicillin resistance (MIC 2-4 mg/L) in *mecA* deficient MRSA has been ascribed to overproduction of a penicillinase[28] or alteration in PBP expression[29]. Infant isolate 18C was further noted to lack the *blaZ* gene, but still exhibited phenotypic penicillin resistance. Nonetheless, this highlights the presence of alternative mechanisms of antibiotic resistance, which merit further investigation into alternative genetic determinants related to resistance.

Control of both outbreak clades was obtained through several measures including temporary closure of the NICU to new admissions, global contact isolation, positive case cohorting, reinforcement of hand hygiene and PPE practices, extensive environmental disinfection, and patient, parent, and HCP decolonization. All MRSA colonized patients and HCP underwent MRSA decolonization with nasal mupirocin (patients) or nasal povidone iodine (HCP) and topical chlorhexidine if they were eligible. Parents of MRSA colonized patients were offered decolonization with nasal povidone iodine and topical chlorhexidine. Ongoing prospective surveillance in the NICU has not identified further clade 1 or 2 associated MRSA isolates 6 months after the outbreak investigation, underscoring the importance and efficacy of decolonization, as has been noted in the case of parental reservoirs[30].

As evidenced in this work and others[31], conventional methods of identifying pathogens and outbreak investigation, which have largely relied on MLST[32], and prior to that, Pulse-field Gel Electrophoresis (PFGE)[10], lack the resolution necessary to identify the source and transmission of an outbreak. For instance, isolate belonging to ID 19C shared identical MLST typing as clade 1, as well as the same AST and SSC*mec* typing. The added resolution of WGS allowed us to accurately identify outbreak clades and infer direction of transmission, which proved pivotal in mitigating the outbreak.

Although the directionality of MRSA transmission cannot be confirmed, both outbreak clades identified were associated with HCP through genetic and epidemiologic linkage. Prior MRSA outbreak studies in the NICU have also identified the potential contribution of the HCP reservoir[10,33], but few, if any, have shown prolonged circulation in the hospital environment through silent transmission events that have led to significant patient morbidity, as this work demonstrates. The high risk of ongoing transmission events originating from undetected MRSA carriers suggests that increasing the frequency of surveillance efforts through genomic diagnostics may result in decreased transmission. This begs the consideration of prospective genomic surveillance that includes the HCP reservoir as a way to identify uniquely virulent MRSA strains with theoretical high epidemic potential (PVL and ACME carriage, etc). While economic modeling studies of prospective genomic MRSA surveillance suggests it’s cost-effectiveness[34], it may still pose various difficulties such as when to include HCP during prospective screening and who should be included, at what intervals should screening occur, and should clinical triggers be used to initiate HCP screening, among others[35]. More research into economic and clinical benefits of HCP inclusion in prospective genomic surveillance is necessary.

This study has limitations that should be considered. While this work highlights potential that silent transmission of MRSA may occur for months before detection of an outbreak, we were unable to sequence all MRSA isolates detected in the NICU during the seven-month gap between sample the index case and the outbreak, nor isolates prior to the index case. This limits our ability to definitively pinpoint the introduction and complex transmission dynamics of outbreak clade 1. Further, the source of introduction of clade 1 and clade 2 remain elusive and necessitates consideration of other potential reservoirs besides the affected HCPs, such as other HCPs, parents, and the environment. However, environmental samples were unavailable. Nevertheless, the high degree of genetic similarity of clade 1 isolates, evidenced by low SNP difference, suggest clade 1 isolates were part of a singular transmission event and thus highlight the dangers of silent transmission of highly virulent isolates can pose. Finally, we were not able to validate the discrepancy between AST and *mecA* exhibited by isolates belonging to ID 17C and 18C.

In summary, our results suggest that transmission between HCP and patients contributed to a NICU outbreak involving a uniquely virulent MRSA strain. WGS enabled data-driven implementation of infection prevention strategies that prevented transmission of the outbreak strains and identified a potential source. Prospective genomic epidemiology of hospital acquired infections may help identify occult transmission events that may precede outbreaks of MRSA and other pathogens. Our findings corroborate prior work demonstrating the utility of WGS in outbreak management[18,36,37]. Given the increasing availability and affordability of WGS, genomic epidemiology of CDC priority pathogens from staff and patients may be a useful measure for identifying early colonization and silent transmission events of virulent strains prior to outbreak emergence.

## Supporting information

Supplemental Figure 1 [[supplements/261307_file02.pdf]](pending:yes)

Supplemental Table 1 [[supplements/261307_file03.pdf]](pending:yes)

## Data Availability

Data is available upon request.NHLBI K23HL138461-01A1 (CL) NIH 5T32AI007641-19 (SM) 

## Footnotes

*   **Funding:** NHLBI K23HL138461-01A1 (CL)

*   NIH 5T32AI007641-19 (SM)

*   Received July 28, 2021.
*   Revision received July 28, 2021.
*   Accepted July 29, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## REFERENCES

1.  1.Lake JG, Weiner LM, Milstone AM, Saiman L, Magill SS, See I. Pathogen Distribution and Antimicrobial Resistance Among Pediatric Healthcare-Associated Infections Reported to the National Healthcare Safety Network, 2011–2014. Infect Control Hosp Epidemiol 2018; 39:1–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/ice.2017.236&link_type=DOI) 

2.  2.Shah J, Jefferies A, Yoon E, Lee S, Shah P, on behalf of the Canadian Neonatal Network. Risk Factors and Outcomes of Late-Onset Bacterial Sepsis in Preterm Neonates Born at < 32 Weeks’ Gestation. Am J Perinatol 2014; 32:675–682.
    
    

3.  3.Song X, Perencevich E, Campos J, Short BL, Singh N. Clinical and Economic Impact of Methicillin-Resistant Staphylococcus aureus Colonization or Infection on Neonates in Intensive Care Units. Infect Control Hosp Epidemiol 2010; 31:177–182.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/649797&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20001732&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 

4.  4.Zervou FN, Zacharioudakis IM, Ziakas PD, Mylonakis E. MRSA Colonization and Risk of Infection in the Neonatal and Pediatric ICU: A Meta-analysis. PEDIATRICS 2014; 133:e1015–e1023.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6MTE6IjEzMy80L2UxMDE1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMjkvMjAyMS4wNy4yOC4yMTI2MTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  5.Wagenvoort JHT, Sluijsmans W, Penders RJR. Better environmental survival of outbreak vs. sporadic MRSA isolates. J Hosp Infect 2000; 45:231–234.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/jhin.2000.0757&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10896803&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088308300009&link_type=ISI) 

6.  6.Jimenez-Truque N, Tedeschi S, Saye EJ, et al. Relationship Between Maternal and Neonatal Staphylococcus aureus Colonization. PEDIATRICS 2012; 129:e1252–e1259.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6MTE6IjEyOS81L2UxMjUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMjkvMjAyMS4wNy4yOC4yMTI2MTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

7.  7.Hollis RJ, Barr JL, Doebbeling BN, Pfaller MA, Wenzel RP. Familial Carriage of Methicillin-Resistant Staphylococcus aureus and Subsequent Infection in a Premature Neonate. Clin Infect Dis 1995; 21:328–332.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/clinids/21.2.328&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8562740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RM81700011&link_type=ISI) 

8.  8.Al-Tawfiq JA. Father-to-Infant Transmission of Community-Acquired Methicillin-Resistant Staphylococcus aureus in a Neonatal Intensive Care Unit. Infect Control Hosp Epidemiol 2006; 27:636–637.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/505097&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16755488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249036800019&link_type=ISI) 

9.  9.Boyce JM, Opal SM, Potter-Bynoe G, Medeiros AA. Spread of Methicillin-Resistant Staphylococcus aureus in a Hospital After Exposure to a Health Care Worker with Chronic Sinusitis. Clin Infect Dis 1993; 17:496–504.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/clinids/17.3.496&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8218696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993LW23700029&link_type=ISI) 

10. 10.Saiman L, Cronquist A, Wu F, et al. An Outbreak of Methicillin-Resistant Staphylococcus aureus in a Neonatal Intensive Care Unit. Infect Control Hosp Epidemiol 2003; 24:317–321.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/502217&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12785403&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182892700007&link_type=ISI) 

11. 11.Bertin ML, Vinski J, Schmitt S, et al. Outbreak of methicillin-resistant Staphylococcus aureus colonization and infection in a neonatal intensive care unit epidemiologically linked to a healthcare worker with chronic otitis. Infect Control Hosp Epidemiol 2006; 27:581–585.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/504933&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16755477&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249036800008&link_type=ISI) 

12. 12.Kamm J. SPID: SNP Pipeline for Infectious Disease. Available at: [https://github.com/czbiohub/Spid.jl](https://github.com/czbiohub/Spid.jl).
    
    

13. 13.Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014; 30:1312–1313.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24451623&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000336095100024&link_type=ISI) 

14. 14.Huerta-Cepas J, Serra F, Bork P. ETE 3: Reconstruction, Analysis, and Visualization of Phylogenomic Data. Mol Biol Evol 2016; 33:1635–1638.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msw046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26921390&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 

15. 15.Gupta SK, Padmanabhan BR, Diene SM, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother 2014; 58:212–220.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjU4LzEvMjEyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMjkvMjAyMS4wNy4yOC4yMTI2MTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

16. 16.Kaya H, Hasman H, Larsen J, et al. SCCmecFinder, a Web-Based Tool for Typing of Staphylococcal Cassette Chromosome mec in Staphylococcus aureus Using Whole-Genome Sequence Data. mSphere 2018; 3:e00612–17.
    
    

17. 17.Larsen MV, Cosentino S, Rasmussen S, et al. Multilocus Sequence Typing of Total-Genome-Sequenced Bacteria. J Clin Microbiol 2012; 50:1355–1361.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjUwLzQvMTM1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA3LzI5LzIwMjEuMDcuMjguMjEyNjEzMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

18. 18.Joensen KG, Scheutz F, Lund O, et al. Real-Time Whole-Genome Sequencing for Routine Typing, Surveillance, and Outbreak Detection of Verotoxigenic Escherichia coli. J Clin Microbiol 2014; 52:1501–1510.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjUyLzUvMTUwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA3LzI5LzIwMjEuMDcuMjguMjEyNjEzMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

19. 19.Uhlemann A-C, Dordel J, Knox JR, et al. 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 2014; 111:6738–6743.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTExLzE4LzY3MzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8yOS8yMDIxLjA3LjI4LjIxMjYxMzA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

20. 20.Gillet Y, Issartel B, Vanhems P, et al. Association between Staphylococcus aureus strains carrying gene for Panton-Valentine leukocidin and highly lethal necrotising pneumonia in young immunocompetent patients. Lancet Lond Engl 2002; 359:753–759.
    
    

21. 21.Planet PJ, LaRussa SJ, Dana A, et al. Emergence of the Epidemic Methicillin-Resistant Staphylococcus aureus Strain USA300 Coincides with Horizontal Transfer of the Arginine Catabolic Mobile Element and speG-mediated Adaptations for Survival on Skin. mBio 2013; 4:e00889–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.00889-13&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24345744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 

22. 22.Yarwood JM, McCormick JK, Paustian ML, Orwin PM, Kapur V, Schlievert PM. Characterization and Expression Analysis of Staphylococcus aureus Pathogenicity Island 3: IMPLICATIONS FOR THE EVOLUTION OF STAPHYLOCOCCAL PATHOGENICITY ISLANDS. J Biol Chem 2002; 277:13138–13147.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNzcvMTUvMTMxMzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8yOS8yMDIxLjA3LjI4LjIxMjYxMzA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

23. 23.Kuroda M, Ohta T, Uchiyama I, et al. Whole genome sequencing of meticillin-resistant Staphylococcus aureus. Lancet Lond Engl 2001; 357:1225–1240.
    
    

24. 24.Letertre C, Perelle S, Dilasser F, Fach P. Identification of a new putative enterotoxin SEU encoded by the egc cluster of Staphylococcus aureus. J Appl Microbiol 2003; 95:38–43.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2672.2003.01957.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12807452&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183584600005&link_type=ISI) 

25. 25.Diep BA, Gill SR, Chang RF, et al. Complete genome sequence of USA300, an epidemic clone of community-acquired meticillin-resistant Staphylococcus aureus. Lancet Lond Engl 2006; 367:731–739.
    
    

26. 26.ThURLow LR, Joshi GS, Clark JR, et al. Functional Modularity of the Arginine Catabolic Mobile Element Contributes to the Success of USA300 Methicillin-Resistant Staphylococcus aureus. Cell Host Microbe 2013; 13:100–107.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2012.11.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23332159&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314262900012&link_type=ISI) 

27. 27.Hiramatsu K. Molecular evolution of MRSA. Microbiol Immunol 1995; 39:531–543.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1348-0421.1995.tb02239.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7494490&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RQ85400001&link_type=ISI) 

28. 28.McDougal LK, Thornsberry C. The role of beta-lactamase in staphylococcal resistance to penicillinase-resistant penicillins and cephalosporins. J Clin Microbiol 1986; 23:832–839.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjg6IjIzLzUvODMyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMjkvMjAyMS4wNy4yOC4yMTI2MTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

29. 29.Tomasz A, Drugeon HB, de Lencastre HM, Jabes D, McDougall L, Bille J. New mechanism for methicillin resistance in Staphylococcus aureus: clinical isolates that lack the PBP 2a gene and contain normal penicillin-binding proteins with modified penicillin-binding capacity. Antimicrob Agents Chemother 1989; 33:1869–1874.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjEwOiIzMy8xMS8xODY5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMjkvMjAyMS4wNy4yOC4yMTI2MTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

30. 30.Milstone AM, Voskertchian A, Koontz DW, et al. Effect of Treating Parents Colonized With Staphylococcus aureus on Transmission to Neonates in the Intensive Care Unit: A Randomized Clinical Trial. JAMA 2020; 323:319.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2019.20785&link_type=DOI) 

31. 31.Hwang SM, Cho HW, Kim TY, et al. Whole-Genome Sequencing for Investigating a Health Care-Associated Outbreak of Carbapenem-Resistant Acinetobacter baumannii. Diagnostics 2021; 11:201.
    
    

32. 32.Yamaji R, Rubin J, Thys E, Friedman CR, Riley LW. Persistent Pandemic Lineages of Uropathogenic Escherichia coli in a College Community from 1999 to 2017. J Clin Microbiol 2018; 56. Available at: [https://journals.asm.org/doi/10.1128/JCM.01834-17](https://journals.asm.org/doi/10.1128/JCM.01834-17). Accessed 7 June 2021.
    
    

33. 33.Bertin ML, Vinski J, Schmitt S, et al. Outbreak of methicillin-resistant Staphylococcus aureus colonization and infection in a neonatal intensive care unit epidemiologically linked to a healthcare worker with chronic otitis. Infect Control Hosp Epidemiol 2006; 27:581–585.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/504933&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16755477&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249036800008&link_type=ISI) 

34. 34.Dymond A, Davies H, Mealing S, et al. Genomic Surveillance of Methicillin-resistant Staphylococcus aureus: A Mathematical Early Modeling Study of Cost-effectiveness. Clin Infect Dis 2020; 70:1613–1619.
    
    

35. 35.Grant PS, Charns LG, Rawot BW, Benedetti SG. Consideration to culture health care workers related to increased methicillin-resistant Staphylococcus aureus activity in a neonatal intensive care unit. Am J Infect Control 2008; 36:638–643.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18834741&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F29%2F2021.07.28.21261307.atom) 

36. 36.Snitkin ES, Zelazny AM, Thomas PJ, et al. Tracking a Hospital Outbreak of Carbapenem-Resistant Klebsiella pneumoniae with Whole-Genome Sequencing. Sci Transl Med 2012; 4:148ra116–148ra116.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InNjaXRyYW5zbWVkIjtzOjU6InJlc2lkIjtzOjE0OiI0LzE0OC8xNDhyYTExNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA3LzI5LzIwMjEuMDcuMjguMjEyNjEzMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

37. 37.Crawford E, Kamm J, Miller S, et al. Investigating Transfusion-related Sepsis Using Culture-Independent Metagenomic Sequencing. Clin Infect Dis 2020; 71:1179–1185.