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
    • COVID-19 Special Collection
    • 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
    • COVID-19 Special Collection
    • 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
Bacteriology

Whole-Genome Sequencing Reveals the Contribution of Long-Term Carriers in Staphylococcus aureus Outbreak Investigation

N. C. Gordon, B. Pichon, T. Golubchik, D. J. Wilson, J. Paul, D. S. Blanc, K. Cole, J. Collins, N. Cortes, M. Cubbon, F. K. Gould, P. J. Jenks, M. Llewelyn, J. Q. Nash, J. M. Orendi, K. Paranthaman, J. R. Price, L. Senn, H. L. Thomas, S. Wyllie, D. W. Crook, T. E. A. Peto, A. S. Walker, A. M. Kearns
Sandra S. Richter, Editor
N. C. Gordon
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
B. Pichon
cAntimicrobial Resistance and Healthcare Associated Infections Reference Unit, Public Health England, Colindale, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
T. Golubchik
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
D. J. Wilson
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
J. Paul
dPublic Health England, Royal Sussex County Hospital, Brighton, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
D. S. Blanc
eLausanne University Hospital, Service of Preventative Medicine, Lausanne, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
K. Cole
dPublic Health England, Royal Sussex County Hospital, Brighton, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
J. Collins
fNewcastle upon Tyne Hospitals NHS Foundation Trust, Newcastle, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
N. Cortes
gPortsmouth Hospitals NHS Trust, Portsmouth, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
M. Cubbon
hBrighton and Sussex University Hospitals NHS Trust, Brighton, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
F. K. Gould
fNewcastle upon Tyne Hospitals NHS Foundation Trust, Newcastle, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
P. J. Jenks
iPlymouth Hospitals NHS Trust, Plymouth, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
M. Llewelyn
hBrighton and Sussex University Hospitals NHS Trust, Brighton, United Kingdom
jBrighton and Sussex Medical School, Falmer, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
J. Q. Nash
kEast Kent Hospitals NHS Foundation Trust, Canterbury, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
J. M. Orendi
lRoyal Stoke University Hospital, University Hospitals of North Midlands NHS Trust, Stoke-on-Trent, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
K. Paranthaman
mPublic Health England, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
J. R. Price
dPublic Health England, Royal Sussex County Hospital, Brighton, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
L. Senn
eLausanne University Hospital, Service of Preventative Medicine, Lausanne, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
H. L. Thomas
mPublic Health England, London, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
S. Wyllie
gPortsmouth Hospitals NHS Trust, Portsmouth, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
D. W. Crook
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
mPublic Health England, London, United Kingdom
nNational Institute for Health Research Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance at University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
T. E. A. Peto
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
nNational Institute for Health Research Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance at University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. S. Walker
aNational Institute for Health Research Oxford Biomedical Research Centre, John Radcliffe Hospital, Oxford, United Kingdom
bNuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
nNational Institute for Health Research Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance at University of Oxford, Oxford, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. M. Kearns
cAntimicrobial Resistance and Healthcare Associated Infections Reference Unit, Public Health England, Colindale, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sandra S. Richter
Cleveland Clinic
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/JCM.00363-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Whole-genome sequencing (WGS) makes it possible to determine the relatedness of bacterial isolates at a high resolution, thereby helping to characterize outbreaks. However, for Staphylococcus aureus, the accumulation of within-host diversity during carriage might limit the interpretation of sequencing data. In this study, we hypothesized the converse, namely, that within-host diversity can in fact be exploited to reveal the involvement of long-term carriers (LTCs) in outbreaks. We analyzed WGS data from 20 historical outbreaks and applied phylogenetic methods to assess genetic relatedness and to estimate the time to most recent common ancestor (TMRCA). The findings were compared with the routine investigation results and epidemiological evidence. Outbreaks with epidemiological evidence for an LTC source had a mean estimated TMRCA (adjusted for outbreak duration) of 243 days (95% highest posterior density interval [HPD], 143 to 343 days) compared with 55 days (95% HPD, 28 to 81 days) for outbreaks lacking epidemiological evidence for an LTC (P = 0.004). A threshold of 156 days predicted LTC involvement with a sensitivity of 0.875 and a specificity of 1. We also found 6/20 outbreaks included isolates with differing antimicrobial susceptibility profiles; however, these had only modestly increased pairwise diversity (mean 17.5 single nucleotide variants [SNVs] [95% confidence interval {CI}, 17.3 to 17.8]) compared with isolates with identical antibiograms (12.7 SNVs [95% CI, 12.5 to 12.8]) (P < 0.0001). Additionally, for 2 outbreaks, WGS identified 1 or more isolates that were genetically distinct despite having the outbreak pulsed-field gel electrophoresis (PFGE) pulsotype. The duration-adjusted TMRCA allowed the involvement of LTCs in outbreaks to be identified and could be used to decide whether screening for long-term carriage (e.g., in health care workers) is warranted. Requiring identical antibiograms to trigger investigation could miss important contributors to outbreaks.

INTRODUCTION

To manage Staphylococcus aureus outbreaks effectively, infection control practitioners need to determine the relatedness of isolates from suspected cases. Whole-genome sequencing (WGS) has shown superior resolution compared with standard typing techniques (spa, pulsed-field gel electrophoresis [PFGE]) when used for individual outbreaks (1–4) and can also provide additional information about resistance, pathogenicity, and population structure (5–8). However, it has been argued that the accumulation of within-host diversity during S. aureus carriage could result in erroneous inferences about transmission. This has been cited as a potential weakness in applying sequencing to S. aureus outbreaks and may lead to the misinterpretation of genuine transmission routes (1, 9, 10).

However, rather than within-host diversity being a limitation on sequencing-based outbreak investigation, it could in fact be exploited to determine whether a long-term carrier is implicated in maintaining an outbreak. This information could be used by infection control practitioners when considering whether or not to deploy extended screening (e.g., of health care workers).

In this study, we tested the hypothesis that WGS can be used to predict the presence of a long-term carrier as an outbreak source. First, we examined individuals with newly acquired S. aureus nasal carriage to ascertain whether diversity is present at acquisition or develops over time. Next, we analyzed 20 S. aureus outbreaks, which were previously investigated using standard typing techniques, to assess the added utility of WGS. Finally, we compared WGS with epidemiological data to determine whether the presence of a long-term carrier maintaining the outbreak could be inferred from the WGS data.

RESULTS

Comparison of within-host diversity in newly acquired and long-term carriage.Eight subjects were identified with ≥3 consecutive bimonthly negative nasal swabs, followed by ≥1 year of swabs consistently positive for S. aureus. All isolates were methicillin-susceptible S. aureus (MSSA), representing 7 spa types, 5 sequence types, and 4 clonal complexes. The median time from the first to last positive sample was 490 days (range, 358 to 727 days). In total, 135 isolates were successfully sequenced from 16 samples. One isolate (case 1219, early sample) failed quality checks and was excluded.

In 6/8 subjects, there was a significant increase in mean pairwise diversity (MPWD) between the first and last samples (P < 0.05) (Fig. 1). In one participant (case 1236), the increase was not significant (P = 0.52), and for another (case 1375), there was a decrease which was marginally significant (P = 0.07). Overall, MPWD increased from 0.88 single nucleotide variants (SNVs) (95% confidence interval [CI], 0.65 to 1.11) to 3.30 (95% CI, 2.92 to 3.68) between the first and last samples (P < 0.001). An analysis of the phylogenetic trees (see supplemental material) showed highly clonal early populations, and in 2 participants, only a single strain was observed. One individual (case 1219) had a more diverse early sample (MPWD, 4.57; 95% CI, 3.10 to 6.04) compared with those of the other participants. This subject's first positive swab was at month 12, and they had completed a course of amoxicillin-clavulanic acid 1 day before their final negative swab. Therefore, it is possible that this was a false negative due to antibiotic suppression, meaning that there may have been up to 4 months of carriage prior to the first positive swab, accounting for the increased diversity.

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

All pairwise differences between early (<2 months since acquisition) and late (>12 months since acquisition) nasal swab samples from 7 patients with previously negative nasal swabs. Included for comparison are samples from patients positive at entry to the study (time of acquisition unknown).

Two participants (cases 1218 and 1219) shared the same address and had isolates of the same spa type. Participant 1219 (donor) became positive 2 months before participant 1218 (recipient). On a direct comparison of both early populations, we found that the recipient had an entirely clonal initial population, identical to 4/8 of the donor's strains (see supplemental material).

For an additional 13 participants positive at the study entry, within-host diversity as measured by MPWD ranged from 0 SNVs (3 individuals) to 26 SNVs. This may be due to differences in the acquisition time to the time of the first sample, which is unknown for these individuals.

Outbreak characteristics.Twenty outbreaks were included in the study (Table 1). Fourteen (70%) were hospital associated, including 5 in neonatal units, 4 in general wards, 1 from a surgical unit, 2 from maternity units, and 2 involved multiple wards or hospital sites. Six (30%) were community associated, including 4 households, 1 nursing home, and 1 school. The reasons for instigating an outbreak investigation were an increase in methicillin-resistant S. aureus (MRSA) carriage (8 outbreaks), Panton-Valentine leukocidin (PVL)-producing skin/soft tissue infection (7 outbreaks), surgical site infections (3 outbreaks), MRSA bacteremia (1 outbreak), and staphylococcal scalded skin syndrome (1 outbreak). Three (15%) were due to MSSA and 17 (85%) to MRSA. The median number of outbreak cases was 7 (interquartile range [IQR], 5 to 9). The median duration was 72 days (IQR, 44 to 188 days).

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

Description of 20 outbreaks analyzed by whole-genome sequencing

Overall, isolates from 391 cases were sequenced. Nine (2.3%) were from health care workers (HCWs), the remainder being from patients or household members. Outbreak samples represented 9 clonal complexes, 11 sequence types, and 12 spa types.

Phylogenetic analysis of outbreaks.Phylogenetic trees for each outbreak are provided in the supplemental material. Two outbreaks had isolates which were equally or more distant than comparator isolates despite having the outbreak pulsotype: outbreak D (one isolate, 53 SNVs from the index case compared with 21) and outbreak S (two isolates, 49 and 46 SNVs from the index case compared with 46). These were therefore considered to be sporadic nonoutbreak isolates and were excluded from further analysis.

The overall MPWD across all outbreak sample pairs for the remaining 388 isolates was 13.8 SNVs (95% CI, 13.6 to 13.9) compared with 4,444 SNVs for nonoutbreak spa-matched pairs (95% CI, 2,492 to 6,395) and 30,192 SNVs for nonoutbreak isolates from the same units (95% CI, 29,781 to 30,603). All outbreak isolates were ≤30 SNVs from the index case; 381/388 (98%) were ≤10 SNVs from their nearest neighbor. The 7 more distant isolates came from outbreaks lasting more than 6 months (B, G, and S). All isolates were mapped to a standard reference genome; mapping to an alternative reference strain (performed for 6 outbreaks) yielded only 2 additional SNVs overall (see supplemental material), with no effect on topology.

Time to most recent common ancestor and long-term carriers.Twelve outbreaks (60%) had epidemiological evidence of a long-term carrier (LTC). Three included cases with recurrent staphylococcal disease, in 5 an LTC was suspected due to nonoverlapping ward stays, and in 4, at least one case had postoutbreak long-term carriage (Fig. 2). The pairwise distances between isolates from outbreaks with evidence for an LTC ranged from 0 to 46 SNVs compared with 0 to 10 SNVs for outbreaks with no evidence for an LTC (Table 2). The mean duration-adjusted time to most recent common ancestor (TMRCA) for outbreaks with a suspected or proven LTC was 243 days (95% highest posterior density interval [HPD], 143 to 343 days) compared with 55 days (95% HPD, 28 to 81 days) for outbreaks with no evidence for an LTC (P = 0.004) (Fig. 2). Excluding postoutbreak carriage, an analysis of the receiver operating characteristic curve gave an area under the curve (AUC) of 0.953 (95% CI, 0.851 to 1). Using the Youden index to select the optimal threshold gave a cutoff value of 156 days, with a sensitivity of 0.875 and a specificity of 1.

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

Duration-adjusted TMRCA for outbreaks with (i) no evidence of a long-term carrier (direct contacts between all cases); (ii) likely LTC (indirect ward contacts or preoutbreak LTC); or (iii) LTC unclear/possible (evidence of a postoutbreak LTC).

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

Long-term carrier category, duration-adjusted TMRCA, and SNV range for outbreaks investigated using WGS

Relationship between PFGE pulsotype/antibiogram and SNV distance.Five outbreaks contained isolates differing by ≥1 band from the index case on PFGE. MPWD between outbreak isolates with identical PFGE pulsotypes was 13.6 SNVs (95% CI, 13.4 to 13.7) compared with 17.3 (95% CI, 17.0 to 17.6) between isolates with differing pulsotypes (P < 0.0001).

In 6/20 outbreaks, antimicrobial susceptibility differed across isolates, confirmed by the presence/absence of mobile resistance determinants identified using BLAST (11); however, these clearly belonged to the outbreak on phylogenetic analysis. MPWD between isolates sharing an antibiogram was 12.7 SNVs (95% CI, 12.5 to 12.8) compared with 17.5 (95% CI, 17.3 to 17.8) for isolates with differing antibiograms (P < 0.0001), although a substantial number of isolate pairs with different antibiograms had 0 SNVs between their core genomes (Fig. 3).

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

Pairwise SNV differences for all pairs within an outbreak, where isolates had differing antibiograms (a) or differing PFGE pulsotypes (b).

For other factors potentially related to outbreak diversity, there was no evidence of an association between MPWD and outbreak duration, reason for investigation, epidemiological setting, or MRSA phenotype (P > 0.05).

DISCUSSION

We have tested the use of WGS for S. aureus outbreak investigation using 20 outbreaks. By comparing observed outbreak SNV distances with nonoutbreak and spa/multilocus sequence type (MLST)-specific diversities, we were able to distinguish outbreak from nonoutbreak strains.

Our observation of minimal diversity in recent acquisitions of nasal carriage is reassuring for the application of WGS data to outbreaks. For the donor-recipient pair, we observed a narrow transmission bottleneck, with a clonal founding population despite a diverse donor population. Although this is a single case, the findings are supported by the minimal diversity seen in the early samples for the majority of carriage study subjects, and further evidence for a narrow transmission bottleneck is provided by the relatively short SNV distances observed across the outbreaks. Taken together, these findings suggest that, in an acute short-term outbreak, there will be insufficient time for diversity to accumulate.

If WGS is to be used routinely for outbreak investigation, these findings provide evidence that single-colony sequencing is likely to identify clusters reliably in this context, allowing for the ease of interpretation and ensuring that WGS remains an affordable alternative to standard typing, as a requirement for sequencing multiple colonies per case, as implied by previous investigators (1, 10), would rapidly escalate costs and render WGS too expensive for routine use.

Previous carriage studies have found greater distances than seen here (9, 11); however, these did not account for the estimated time of acquisition. We postulate that the existence of a significant cloud of diversity (4, 12) may be a marker of long-term carriage. Therefore, in outbreaks, higher diversity may indicate the involvement of an LTC, with outbreak diversity reflecting the donor cloud.

In support of this, we observed a significant difference in duration-adjusted TMRCAs between outbreaks with and without evidence of an LTC. The longest TMRCAs were in hospital outbreaks with indirect links between cases (i.e., nonoverlapping ward stays). The likelihood of “missed” cases in these outbreaks was considered low due to the enhanced screening, and the most likely reason for the reoccurrence of the outbreak strain was thought by the investigating teams to be either a reintroduction from the community (outbreak G) or from a staff member with long-term carriage (outbreaks A, I, N, and S). Staff carriage was proven in one outbreak (by sampling and subsequent termination of the outbreak on their exclusion), but in the remaining outbreaks, HCWs were either not sampled or HCW sampling was anonymized and positive results could not be linked definitively with the suspected carrier.

The outbreaks included in this study necessarily reflect the circulating S. aureus clones in the United Kingdom and the concerns of local infection control teams. The sampling frame is therefore enriched for MRSA and PVL-positive outbreaks and those from neonatal units. Despite this, there is a wide representation of sequence types.

In spite of the enhanced surveillance during each outbreak, there inevitably are missing transmission links, due to missed sampling, suppression from antimicrobial therapy, or delays in identifying contacts. One reason for missed samples may be the use of antibiograms as an initial screening tool for identifying putative outbreak isolates, as most investigating teams only collected isolates with identical or highly similar antimicrobial susceptibility profiles. However, in the six outbreaks where isolates were included with differing antibiograms, the core genomes were remarkably conserved. This is presumably due to the ready loss/gain of mobile genetic elements (13) and shows that reliance on antibiograms may lead to samples being wrongly excluded.

The variability of mobile elements is also important for interpreting genetic distances. Recombination events such as the gain/loss of a mobile element will introduce a large number of SNVs even though this represents a single genetic event. Current analysis tools which can accommodate this are computationally complex and, for large data sets, require sizable computing resources. A simpler approach is to exclude the “mobile-ome” from phylogenetic analyses and compare only the core genome, and the results above demonstrate that this is an acceptable strategy. Similarly, mapping to alternative reference strains (performed for six outbreaks) had minimal effects on SNV analysis and phylogeny, removing the need for identification of clonal complex or assembly of index case prior to phylogenetic analysis. This streamlined approach brings WGS closer to routine use, as a readily deployable method with a minimal burden of computational time and bioinformatics expertise.

In conclusion, we have demonstrated how a WGS-based approach can be applied to S. aureus outbreak investigations. We have shown that current sampling strategies provide sufficient information to determine whether isolates belong to an outbreak, and that, rather than confounding the investigation, within-host diversity can be utilized to identify the possible involvement of a long-term carrier, potentially enhancing the infection control response. Combining this with directed multisampling of suspected LTCs (1) may be a cost-effective method of using WGS to ensure that, where HCWs are implicated, potentially career-altering decisions are made using the best possible evidence.

MATERIALS AND METHODS

Comparison of within-host diversity in newly acquired and long-term carriage.Eight participants were identified from a population study of S. aureus nasal carriage in adults attending general practices in Oxfordshire (15), in which participants had nasal swabs taken at two-month intervals, with positive samples stored as mixed glycerol stocks taken by sweeping across multiple colonies on the primary plates to preserve the diversity of carried strains (11). The eight participants were negative for nasal carriage at recruitment and subsequently had consistently negative swabs for ≥6 months before acquiring a strain which they carried for at least 1 year. The first and last positive samples for each individual were retrieved from the mixed glycerol stocks. Samples were plated on Columbia blood agar (CBA) and incubated overnight at 37°C. For each time point, 8 individual colonies (12 for case no. 1218) were selected and further subcultured on a CBA plate and again incubated overnight at 37°C.

We also retrieved sequencing data from 13 participants previously investigated, (9) for whom the approximate time of acquisition was unknown. Each of these had 8 to 12 individual colonies sequenced.

Collection of outbreak isolates and epidemiological data.Nineteen outbreaks were purposively sampled in collaboration with the Public Health England (PHE) staphylococcal reference laboratory, representing a range of sequence types and epidemiological settings and including both MRSA and MSSA. One further outbreak was investigated in conjunction with Lausanne University Hospital, Switzerland (14, 16). Epidemiological information was obtained from each infection control team (specimen date, site, ward location, and where applicable, admission/discharge dates and previous screening results).

For each outbreak, additional background isolates were also included for comparison. We sequenced all isolates submitted to PHE as part of the outbreak investigation, including those identified as “nonoutbreak” by routine typing, to estimate the expected genetic diversity of the outbreak strain and to ensure that the apparent outbreak strains were not part of an ongoing clonal expansion. We also included non-epidemiologically linked isolates matched for spa type and/or MLST to provide a comparison for expected within-spa distances and to provide an outgroup for phylogenetic analysis.

Isolates were retrieved from single-colony frozen stocks held at the PHE reference laboratory, Colindale, London, or at Lausanne University Hospital. We used only the first isolate from each case and included isolates both from clinical samples and screening swabs.

Extraction and sequencing.DNA was extracted and sequenced as previously described (6) from a single colony subcultured on CBA and incubated for 18 to 24 h. Sequencing was performed using the Illumina HiSeq or MiSeq platforms.

Genome assembly and construction of phylogenetic trees.For all outbreaks, reads were aligned using Stampy v1.0.17 to a standard reference genome (MRSA252; GenBank no. NC_002952 ) (17). Six outbreaks were also mapped to clonal complex-specific reference genomes obtained from in-house collections or GenBank. Single nucleotide variants were identified across all mapped nonrepetitive sites using SAMtools v 0.1.18 mpileup, with the extended base-alignment quality flag and masking of mobile genetic elements. A consensus of ≥75% and ≥5 reads, including one in each direction, were required to support an SNV, and calls were required to be homozygous under a diploid model. Maximum likelihood trees were estimated from the mapped whole genomes using PhyML (18).

Outbreak analysis and calculation of TMRCA.The index case was defined as the earliest microbiologically confirmed case in each cluster. Outbreak cases were defined as those sharing related PFGE pulsotypes (19) plus a definite epidemiological link to the index or secondary cases (>24-h stay in the same ward or household/classroom/similar community situation with prolonged contact, e.g., childcare). For each outbreak case, the genetic distance in SNVs was calculated from the index case and the nearest neighbor. If an isolate was more distant from the index case than the nearest spa/MLST-matched comparator, it was considered sporadic and excluded from further outbreak analysis.

We classified each outbreak according to the possibility of long-term carrier involvement (LTC carrying for ≥6 weeks) as follows: (i) LTC not suspected, direct contact between cases, or no history of preexisting staphylococcal disease; (ii) evidence for a pre- or perioutbreak LTC, either ≥1 case with a history of recurrent staphylococcal disease or nonoverlapping hospital stays (ward case identified after a case-free interval, indicating a possible health care worker carrier); and (iii) evidence of a postoutbreak LTC, ≥1 case with positive nasal swab >6 weeks after initial swab (indicating a propensity for long-term carriage).

To evaluate the relationship between outbreak diversity and the likelihood of a long-term carrier, we estimated the time to most recent common ancestor (TMRCA) using BEAST v1.8.1 (20). We applied a simple HKY substitution model with constant population size and a standardized substitution rate of 3.3 × 10−6 substitutions per genome per year (7) (see supplemental material). To control for differences in outbreak duration, outbreaks were censored at 6 months, and the (censored) outbreak duration was subtracted from the calculated TMRCA to obtain a duration-adjusted TMRCA.

We compared SNV distances between isolates of identical pulsotypes and those differing by one or more band. To determine whether there was an increase in genetic diversity associated with the acquisition of antimicrobial resistance, we also interrogated the predicted antibiograms as previously described (21).

Statistical analyses were performed using Stata v13.1. Mean pairwise differences were modeled using normal linear regression using robust standard errors to account for dependence within person/outbreak. The ability of TMRCA to differentiate between outbreaks with evidence for an LTC compared with outbreaks with no evidence for an LTC was evaluated using a receiver operating characteristic curve analysis.

Accession number(s).The sequences reported in this paper have been deposited in the NCBI Sequence Read Archive under bioproject number PRJNA380544 .

ACKNOWLEDGMENTS

We thank the staff of the PHE Staphylococcal reference service (especially Marjorie Ganner) for help in identifying and locating samples. We also thank the infection control teams at each site for permission to use isolates and epidemiological data, in particular, the team at East Kent Hospitals NHS Foundation Trust.

This study was supported by the Health Innovation Challenge Fund (a parallel funding partnership between the Wellcome Trust [WT098615/Z/12/Z] and the Department of Health [grant HICF-T5-358]) and by the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Healthcare Associated Infections and Antimicrobial Resistance at University of Oxford in partnership with Public Health England (PHE) (grant number HPRU-2012-10041). The investigation of the Swiss outbreak was supported by a grant from the Swiss Foundation for Research (no. 31003A_150029).

N.C.G. is an MRC doctoral research fellow. D.W.C. and T.E.A.P. are NIHR Senior Investigators. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, the Department of Health, Wellcome Trust, the Medical Research Council, or Public Health England.

FOOTNOTES

    • Received 1 March 2017.
    • Returned for modification 25 March 2017.
    • Accepted 22 April 2017.
    • Accepted manuscript posted online 3 May 2017.
  • Supplemental material for this article may be found at https://doi.org/10.1128/JCM.00363-17 .

  • Copyright © 2017 Gordon 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. Harris SR,
    2. Cartwright EJ,
    3. Torok ME,
    4. Holden MT,
    5. Brown NM,
    6. Ogilvy-Stuart AL,
    7. Ellington MJ,
    8. Quail MA,
    9. Bentley SD,
    10. Parkhill J,
    11. Peacock SJ
    . 2013. Whole-genome sequencing for analysis of an outbreak of meticillin-resistant Staphylococcus aureus: a descriptive study. Lancet Infect Dis13:130–136. doi:10.1016/S1473-3099(12)70268-2.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Koser CU,
    2. Holden MT,
    3. Ellington MJ,
    4. Cartwright EJ,
    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 Med366:2267–2275. doi:10.1056/NEJMoa1109910.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Eyre DW,
    2. Golubchik T,
    3. Gordon NC,
    4. Bowden R,
    5. Piazza P,
    6. Batty EM,
    7. Ip CL,
    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 TE,
    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 Open2:e001124. doi:10.1136/bmjopen-2012-001124.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Tong SY,
    2. Holden MT,
    3. Nickerson EK,
    4. Cooper BS,
    5. Koser CU,
    6. Cori A,
    7. Jombart T,
    8. Cauchemez S,
    9. Fraser C,
    10. Wuthiekanun V,
    11. Thaipadungpanit J,
    12. Hongsuwan M,
    13. Day NP,
    14. Limmathurotsakul D,
    15. Parkhill J,
    16. Peacock SJ
    . 2015. Genome sequencing defines phylogeny and spread of methicillin-resistant Staphylococcus aureus in a high transmission setting. Genome Res25:111–118. doi:10.1101/gr.174730.114.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Alam MT,
    2. Petit RA III,
    3. Crispell EK,
    4. Thornton TA,
    5. Conneely KN,
    6. Jiang Y,
    7. Satola SW,
    8. Read TD
    . 2014. Dissecting vancomycin-intermediate resistance in staphylococcus aureus using genome-wide association. Genome Biol Evol6:1174–1185. doi:10.1093/gbe/evu092.
    OpenUrlCrossRefPubMed
  6. 6.↵
    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 CL,
    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 A109:4550–4555. doi:10.1073/pnas.1113219109.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Harris SR,
    2. Feil EJ,
    3. Holden MT,
    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. Science327:469–474. doi:10.1126/science.1182395.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Holden MT,
    2. Hsu LY,
    3. Kurt K,
    4. Weinert LA,
    5. Mather AE,
    6. Harris SR,
    7. Strommenger B,
    8. Layer F,
    9. Witte W,
    10. de Lencastre H,
    11. Skov R,
    12. Westh H,
    13. Zemlickova H,
    14. Coombs G,
    15. Kearns AM,
    16. Hill RL,
    17. Edgeworth J,
    18. Gould I,
    19. Gant V,
    20. Cooke J,
    21. Edwards GF,
    22. McAdam PR,
    23. Templeton KE,
    24. McCann A,
    25. Zhou Z,
    26. Castillo-Ramirez S,
    27. Feil EJ,
    28. Hudson LO,
    29. Enright MC,
    30. Balloux F,
    31. Aanensen DM,
    32. Spratt BG,
    33. Fitzgerald JR,
    34. Parkhill J,
    35. Achtman M,
    36. Bentley SD,
    37. Nubel U
    . 2013. A genomic portrait of the emergence, evolution, and global spread of a methicillin-resistant Staphylococcus aureus pandemic. Genome Res23:653–664. doi:10.1101/gr.147710.112.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Golubchik T,
    2. Batty EM,
    3. Miller RR,
    4. Farr H,
    5. Young BC,
    6. Larner-Svensson H,
    7. Fung R,
    8. Godwin H,
    9. Knox K,
    10. Votintseva A,
    11. Everitt RG,
    12. Street T,
    13. Cule M,
    14. Ip CL,
    15. Didelot X,
    16. Peto TE,
    17. Harding RM,
    18. Wilson DJ,
    19. Crook DW,
    20. Bowden R
    . 2013. Within-host evolution of Staphylococcus aureus during asymptomatic carriage. PLoS One8:e61319. doi:10.1371/journal.pone.0061319.
    OpenUrlCrossRef
  10. 10.↵
    1. Worby CJ,
    2. Lipsitch M,
    3. Hanage WP
    . 2014. Within-host bacterial diversity hinders accurate reconstruction of transmission networks from genomic distance data. PLoS Comput Biol10:e1003549. doi:10.1371/journal.pcbi.1003549.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Votintseva AA,
    2. Miller RR,
    3. Fung R,
    4. Knox K,
    5. Godwin H,
    6. Peto TE,
    7. Crook DW,
    8. Bowden R,
    9. Walker AS
    . 2014. Multiple-strain colonization in nasal carriers of Staphylococcus aureus. J Clin Microbiol52:1192–1200. doi:10.1128/JCM.03254-13.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Paterson GK,
    2. Harrison EM,
    3. Murray GG,
    4. Welch JJ,
    5. Warland JH,
    6. Holden MT,
    7. Morgan FJ,
    8. Ba X,
    9. Koop G,
    10. Harris SR,
    11. Maskell DJ,
    12. Peacock SJ,
    13. Herrtage ME,
    14. Parkhill J,
    15. Holmes MA
    . 2015. Capturing the cloud of diversity reveals complexity and heterogeneity of MRSA carriage, infection and transmission. Nat Commun6:6560. doi:10.1038/ncomms7560.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. McCarthy AJ,
    2. Loeffler A,
    3. Witney AA,
    4. Gould KA,
    5. Lloyd DH,
    6. Lindsay JA
    . 2014. Extensive horizontal gene transfer during Staphylococcus aureus co-colonization in vivo. Genome Biol Evol6:2697–2708. doi:10.1093/gbe/evu214.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Vogel V,
    2. Falquet L,
    3. Calderon-Copete SP,
    4. Basset P,
    5. Blanc DS
    . 2012. Short term evolution of a highly transmissible methicillin-resistant Staphylococcus aureus clone (ST228) in a tertiary care hospital. PLoS One7:e38969. doi:10.1371/journal.pone.0038969.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Miller RR,
    2. Walker AS,
    3. Godwin H,
    4. Fung R,
    5. Votintseva A,
    6. Bowden R,
    7. Mant D,
    8. Peto TE,
    9. Crook DW,
    10. Knox K
    . 2014. Dynamics of acquisition and loss of carriage of Staphylococcus aureus strains in the community: the effect of clonal complex. J Infect68:426–439. doi:10.1016/j.jinf.2013.12.013.
    OpenUrlCrossRef
  16. 16.↵
    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. mBio7:e02039-15. doi:10.1128/mBio.02039-15.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Holden MT,
    2. Feil EJ,
    3. Lindsay JA,
    4. Peacock SJ,
    5. Day NP,
    6. Enright MC,
    7. Foster TJ,
    8. Moore CE,
    9. Hurst L,
    10. Atkin R,
    11. Barron A,
    12. Bason N,
    13. Bentley SD,
    14. Chillingworth C,
    15. Chillingworth T,
    16. Churcher C,
    17. Clark L,
    18. Corton C,
    19. Cronin A,
    20. Doggett J,
    21. Dowd L,
    22. Feltwell T,
    23. Hance Z,
    24. Harris B,
    25. Hauser H,
    26. Holroyd S,
    27. Jagels K,
    28. James KD,
    29. Lennard N,
    30. Line A,
    31. Mayes R,
    32. Moule S,
    33. Mungall K,
    34. Ormond D,
    35. Quail MA,
    36. Rabbinowitsch E,
    37. Rutherford K,
    38. Sanders M,
    39. Sharp S,
    40. Simmonds M,
    41. Stevens K,
    42. Whitehead S,
    43. Barrell BG,
    44. Spratt BG,
    45. Parkhill J
    . 2004. Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance. Proc Natl Acad Sci U S A101:9786–9791. doi:10.1073/pnas.0402521101.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Guindon S,
    2. Delsuc F,
    3. Dufayard JF,
    4. Gascuel O
    . 2009. Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol537:113–137. doi:10.1007/978-1-59745-251-9_6.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Tenover FC,
    2. Arbeit RD,
    3. Goering RV,
    4. Mickelsen PA,
    5. Murray BE,
    6. Persing DH,
    7. Swaminathan B
    . 1995. Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J Clin Microbiol33:2233–2239.
    OpenUrlFREE Full Text
  20. 20.↵
    1. Drummond AJ,
    2. Rambaut A
    . 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol7:214. doi:10.1186/1471-2148-7-214.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Gordon NC,
    2. Price JR,
    3. Cole K,
    4. Everitt R,
    5. Morgan M,
    6. Finney J,
    7. Kearns AM,
    8. Pichon B,
    9. Young B,
    10. Wilson DJ,
    11. Llewelyn MJ,
    12. Paul J,
    13. Peto TE,
    14. Crook DW,
    15. Walker AS,
    16. Golubchik T
    . 2014. Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing. J Clin Microbiol52:1182–1191. doi:10.1128/JCM.03117-13.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Whole-Genome Sequencing Reveals the Contribution of Long-Term Carriers in Staphylococcus aureus Outbreak Investigation
N. C. Gordon, B. Pichon, T. Golubchik, D. J. Wilson, J. Paul, D. S. Blanc, K. Cole, J. Collins, N. Cortes, M. Cubbon, F. K. Gould, P. J. Jenks, M. Llewelyn, J. Q. Nash, J. M. Orendi, K. Paranthaman, J. R. Price, L. Senn, H. L. Thomas, S. Wyllie, D. W. Crook, T. E. A. Peto, A. S. Walker, A. M. Kearns
Journal of Clinical Microbiology Jun 2017, 55 (7) 2188-2197; DOI: 10.1128/JCM.00363-17

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.
Whole-Genome Sequencing Reveals the Contribution of Long-Term Carriers in Staphylococcus aureus Outbreak Investigation
(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.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Whole-Genome Sequencing Reveals the Contribution of Long-Term Carriers in Staphylococcus aureus Outbreak Investigation
N. C. Gordon, B. Pichon, T. Golubchik, D. J. Wilson, J. Paul, D. S. Blanc, K. Cole, J. Collins, N. Cortes, M. Cubbon, F. K. Gould, P. J. Jenks, M. Llewelyn, J. Q. Nash, J. M. Orendi, K. Paranthaman, J. R. Price, L. Senn, H. L. Thomas, S. Wyllie, D. W. Crook, T. E. A. Peto, A. S. Walker, A. M. Kearns
Journal of Clinical Microbiology Jun 2017, 55 (7) 2188-2197; DOI: 10.1128/JCM.00363-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Carrier State
Disease Outbreaks
molecular typing
Staphylococcal Infections
Staphylococcus aureus
Whole Genome Sequencing
MRSA
Staphylococcus aureus
outbreaks
whole-genome sequencing

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

 

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

 

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

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