Evaluating the impact of genomic epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) on hospital infection prevention and control decisions.

Genomic epidemiology enhances the ability to detect and refute methicillin-resistant Staphylococcus aureus (MRSA) outbreaks in healthcare settings, but its routine introduction requires further evidence of benefits for patients and resource utilization. We performed a 12 month prospective study at Cambridge University Hospitals NHS Foundation Trust in the UK to capture its impact on hospital infection prevention and control (IPC) decisions. MRSA-positive samples were identified via the hospital microbiology laboratory between November 2018 and November 2019. We included samples from in-patients, clinic out-patients, people reviewed in the Emergency Department and healthcare workers screened by Occupational Health. We sequenced the first MRSA isolate from 823 consecutive individuals, defined their pairwise genetic relatedness, and sought epidemiological links in the hospital and community. Genomic analysis of 823 MRSA isolates identified 72 genetic clusters of two or more isolates containing 339/823 (41 %) of the cases. Epidemiological links were identified between two or more cases for 190 (23 %) individuals in 34/72 clusters. Weekly genomic epidemiology updates were shared with the IPC team, culminating in 49 face-to-face meetings and 21 written communications. Seventeen clusters were identified that were consistent with hospital MRSA transmission, discussion of which led to additional IPC actions in 14 of these. Two outbreaks were also identified where transmission had occurred in the community prior to hospital presentation; these were escalated to relevant IPC teams. We identified 38 instances where two or more in-patients shared a ward location on overlapping dates but carried unrelated MRSA isolates (pseudo-outbreaks); research data led to de-escalation of investigations in six of these. Our findings provide further support for the routine use of genomic epidemiology to enhance and target IPC resources.

Antibiotic resistance determination using Enterococcus faecium whole-genome sequences: a diagnostic accuracy study using genotypic and phenotypic data.

BACKGROUND: DNA sequencing could become an alternative to in vitro antibiotic susceptibility testing (AST) methods for determining antibiotic resistance by detecting genetic determinants associated with decreased antibiotic susceptibility. Here, we aimed to assess and improve the accuracy of antibiotic resistance determination from Enterococcus faecium genomes for diagnosis and surveillance purposes. METHODS: In this retrospective diagnostic accuracy study, we first conducted a literature search in PubMed on Jan 14, 2021, to compile a catalogue of genes and mutations predictive of antibiotic resistance in E faecium. We then evaluated the diagnostic accuracy of this database to determine susceptibility to 12 different, clinically relevant antibiotics using a diverse population of 4382 E faecium isolates with available whole-genome sequences and in vitro culture-based AST phenotypes. Isolates were obtained from various sources in 11 countries worldwide between 2000 and 2018. We included isolates tested with broth microdilution, Vitek 2, and disc diffusion, and antibiotics with at least 50 susceptible and 50 resistant isolates. Phenotypic resistance was derived from raw minimum inhibitory concentrations and measured inhibition diameters, and harmonised primarily using the breakpoints set by the European Committee on Antimicrobial Susceptibility Testing. A bioinformatics pipeline was developed to process raw sequencing reads, identify antibiotic resistance genetic determinants, and report genotypic resistance. We used our curated database, as well as ResFinder, AMRFinderPlus, and LRE-Finder, to assess the accuracy of genotypic predictions against phenotypic resistance. FINDINGS: We curated a catalogue of 228 genetic markers involved in resistance to 12 antibiotics in E faecium. Very accurate genotypic predictions were obtained for ampicillin (sensitivity 99·7% [95% CI 99·5-99·9] and specificity 97·9% [95·8-99·0]), ciprofloxacin (98·0% [96·4-98·9] and 98·8% [95·9-99·7]), vancomycin (98·8% [98·3-99·2] and 98·8% [98·0-99·3]), and linezolid resistance (after re-testing false negatives: 100·0% [90·8-100·0] and 98·3% [97·8-98·7]). High sensitivity was obtained for tetracycline (99·5% [99·1-99·7]), teicoplanin (98·9% [98·4-99·3]), and high-level resistance to aminoglycosides (97·7% [96·6-98·4] for streptomycin and 96·8% [95·8-97·5] for gentamicin), although at lower specificity (60-90%). Sensitivity was expectedly low for daptomycin (73·6% [65·1-80·6]) and tigecycline (38·3% [27·1-51·0]), for which the genetic basis of resistance is not fully characterised. Compared with other antibiotic resistance databases and bioinformatic tools, our curated database was similarly accurate at detecting resistance to ciprofloxacin and linezolid and high-level resistance to streptomycin and gentamicin, but had better sensitivity for detecting resistance to ampicillin, tigecycline, daptomycin, and quinupristin-dalfopristin, and better specificity for ampicillin, vancomycin, teicoplanin, and tetracycline resistance. In a validation dataset of 382 isolates, similar or improved diagnostic accuracies were also achieved. INTERPRETATION: To our knowledge, this work represents the largest published evaluation to date of the accuracy of antibiotic susceptibility predictions from E faecium genomes. The results and resources will facilitate the adoption of whole-genome sequencing as a tool for the diagnosis and surveillance of antimicrobial resistance in E faecium. A complete characterisation of the genetic basis of resistance to last-line antibiotics, and the mechanisms mediating antibiotic resistance silencing, are needed to close the remaining sensitivity and specificity gaps in genotypic predictions. FUNDING: Wellcome Trust, UK Department of Health, British Society for Antimicrobial Chemotherapy, Academy of Medical Sciences and the Health Foundation, Medical Research Council Newton Fund, Vietnamese Ministry of Science and Technology, and European Society of Clinical Microbiology and Infectious Disease.

The clinical, genomic, and microbiological profile of invasive multi-drug resistant Escherichia coli in a major teaching hospital in the United Kingdom.

Escherichia coli is a ubiquitous component of the human gut microbiome, but is also a common pathogen, causing around 40, 000 bloodstream infections (BSI) in the United Kingdom (UK) annually. The number of E. coli BSI has increased over the last decade in the UK, and emerging antimicrobial resistance (AMR) profiles threaten treatment options. Here, we combined clinical, epidemiological, and whole genome sequencing data with high content imaging to characterise over 300 E. coli isolates associated with BSI in a large teaching hospital in the East of England. Overall, only a limited number of sequence types (ST) were responsible for the majority of organisms causing invasive disease. The most abundant (20 % of all isolates) was ST131, of which around 90 % comprised the pandemic O25b:H4 group. ST131-O25b:H4 isolates were frequently multi-drug resistant (MDR), with a high prevalence of extended spectrum ß-lactamases (ESBL) and fluoroquinolone resistance. There was no association between AMR phenotypes and the source of E. coli bacteraemia or whether the infection was healthcare-associated. Several clusters of ST131 were genetically similar, potentially suggesting a shared transmission network. However, there was no clear epidemiological associations between these cases, and they included organisms from both healthcare-associated and non-healthcare-associated origins. The majority of ST131 isolates exhibited strong binding with an anti-O25b antibody, raising the possibility of developing rapid diagnostics targeting this pathogen. In summary, our data suggest that a restricted set of MDR E. coli populations can be maintained and spread across both community and healthcare settings in this location, contributing disproportionately to invasive disease and AMR.

Dissemination of carbapenemase-producing Enterobacterales in Ireland from 2012 to 2017: a retrospective genomic surveillance study.

The spread of carbapenemase-producing Enterobacterales (CPE) is of major public health concern. The transmission dynamics of CPE in hospitals, particularly at the national level, are not well understood. Here, we describe a retrospective nationwide genomic surveillance study of CPE in Ireland between 2012 and 2017. We sequenced 746 national surveillance CPE samples obtained between 2012 and 2017. After clustering the sequences, we used thresholds based on pairwise SNPs, and reported within-host diversity along with epidemiological data to infer recent putative transmissions. All clusters in circulating clones, derived from high-resolution phylogenies, of a species (Klebsiella pneumoniae, Escherichia coli, Klebsiella oxytoca, Enterobacter cloacae, Enterobacter hormaechei and Citrobacter freundii) were individually examined for evidence of transmission. Antimicrobial resistance trends over time were also assessed. We identified 352 putative transmission events in six species including widespread and frequent transmissions in three species. We detected putative outbreaks in 4/6 species with three hospitals experiencing prolonged outbreaks. The bla OXA-48 gene was the main cause of carbapenem resistance in Ireland in almost all species. An expansion in the number of sequence types carrying bla OXA-48 was an additional cause of the increasing prevalence of carbapenemase-producing K. pneumoniae and E. coli.

Multi-Drug Resistant Staphylococcus aureus Carriage in Abattoir Workers in Busia, Kenya.

Abattoir workers have been identified as high-risk for livestock-associated Staphylococcus aureus carriage. This study investigated S. aureus carriage in abattoir workers in Western Kenya. Nasal swabs were collected once from participants between February-November 2012. S. aureus was isolated using bacterial culture and antibiotic susceptibility testing performed using the VITEK 2 instrument and disc diffusion methods. Isolates underwent whole genome sequencing and Multi Locus Sequence Types were derived from these data. S. aureus (n = 126) was isolated from 118/737 (16.0%) participants. Carriage was higher in HIV-positive (24/89, 27.0%) than HIV−negative participants (94/648, 14.5%; p = 0.003). There were 23 sequence types (STs) identified, and half of the isolates were ST152 (34.1%) or ST8 (15.1%). Many isolates carried the Panton-Valentine leucocidin toxin gene (42.9%). Only three isolates were methicillin resistant S. aureus (MRSA) (3/126, 2.4%) and the prevalence of MRSA carriage was 0.4% (3/737). All MRSA were ST88. Isolates from HIV-positive participants (37.0%) were more frequently resistant to sulfamethoxazole/trimethoprim compared to isolates from HIV-negative participants (6.1%; p < 0.001). Similarly, trimethoprim resistance genes were more frequently detected in isolates from HIV-positive (81.5%) compared to HIV-negative participants (60.6%; p = 0.044). S. aureus in abattoir workers were representative of major sequence types in Africa, with a high proportion being toxigenic isolates. HIV-positive individuals were more frequently colonized by antimicrobial resistant S. aureus which may be explained by prophylactic antimicrobial use.

Characteristics of Staphylococcus aureus Isolated from Patients in Busia County Referral Hospital, Kenya.

Staphylococcus aureus is an important pathogen associated with hospital, community, and livestock-acquired infections, with the ability to develop resistance to antibiotics. Nasal carriage by hospital inpatients is a risk for opportunistic infections. Antibiotic susceptibility patterns, virulence genes and genetic population structure of S. aureus nasal isolates, from inpatients at Busia County Referral Hospital (BCRH) were analyzed. A total of 263 inpatients were randomly sampled, from May to July 2015. The majority of inpatients (85.9%) were treated empirically with antimicrobials, including ceftriaxone (65.8%) and metronidazole (49.8%). Thirty S. aureus isolates were cultured from 29 inpatients with a prevalence of 11% (10.3% methicillin-susceptible S. aureus (MSSA), 0.8% methicillin resistant S. aureus (MRSA)). Phenotypic and genotypic resistance was highest to penicillin-G (96.8%), trimethoprim (73.3%), and tetracycline (13.3%) with 20% of isolates classified as multidrug resistant. Virulence genes, Panton-Valentine leukocidin (pvl), toxic shock syndrome toxin-1 (tsst-1), and sasX gene were detected in 16.7%, 23.3% and 3.3% of isolates. Phylogenetic analysis showed 4 predominant clonal complexes CC152, CC8, CC80, and CC508. This study has identified that inpatients of BCRH were carriers of S. aureus harbouring virulence genes and resistance to a range of antibiotics. This may indicate a public health risk to other patients and the community.

Large-Scale Evaluation of a Rapid Fully Automated Analysis Platform to Detect and Refute Outbreaks Based on MRSA Genome Comparisons.

Genomic epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) could transform outbreak investigations, but its clinical introduction is hampered by the lack of automated data analysis tools to rapidly and accurately define transmission based on sequence relatedness. We aimed to evaluate a fully automated bioinformatics system for MRSA genome analysis versus a bespoke researcher-led manual informatics pipeline. We analyzed 781 MRSA genomes from 777 consecutive patients identified over a 9-month period in a clinical microbiology laboratory in the United Kingdom. Outputs were bacterial species identification, detection of mec genes, assignment to sequence types (STs), identification of pairwise relatedness using a definition of ≤25 single nucleotide polymorphisms (SNPs) apart, and use of genetic relatedness to identify clusters. There was full concordance between the two analysis methods for species identification, detection of mec genes, and ST assignment. A total of 3,311 isolate pairs ≤25 SNPs apart were identified by at least one method. These had a median (range) SNP difference between the two methods of 1.2 SNPs (0 to 22 SNPs), with most isolate pairs (87%) varying by ≤2 SNPs. This similarity increased when the research pipeline was modified to use a clonal-complex-specific reference (median 0 SNP difference, 91% varying by ≤2 SNPs). Both pipelines clustered 338 isolates/334 patients into 66 unique clusters based on genetic relatedness. We conclude that the automated transmission detection tool worked at least as well as a researcher-led manual analysis and indicates how such tools could support the rapid use of MRSA genomic epidemiology in infection control practice. IMPORTANCE It has been clearly established that genome sequencing of MRSA improves the accuracy of health care-associated outbreak investigations, including the confirmation and exclusion of outbreaks and identification of patients involved. This could lead to more targeted infection control actions but its use in clinical practice is prevented by several barriers, one of which is the availability of genome analysis tools that do not depend on specialist knowledge to analyze or interpret the results. We evaluated a prototype of a fully automated bioinformatics system for MRSA genome analysis versus a bespoke researcher-led manual informatics pipeline, using genomes from 777 patients over a period of 9 months. The performance was at least equivalent to the researcher-led manual genomic analysis. This indicates the feasibility of automated analysis and represents one more step toward the routine use of pathogen sequencing in infection prevention and control practice.

Genomic epidemiology of SARS-CoV-2 in a UK university identifies dynamics of transmission.

Understanding SARS-CoV-2 transmission in higher education settings is important to limit spread between students, and into at-risk populations. In this study, we sequenced 482 SARS-CoV-2 isolates from the University of Cambridge from 5 October to 6 December 2020. We perform a detailed phylogenetic comparison with 972 isolates from the surrounding community, complemented with epidemiological and contact tracing data, to determine transmission dynamics. We observe limited viral introductions into the university; the majority of student cases were linked to a single genetic cluster, likely following social gatherings at a venue outside the university. We identify considerable onward transmission associated with student accommodation and courses; this was effectively contained using local infection control measures and following a national lockdown. Transmission clusters were largely segregated within the university or the community. Our study highlights key determinants of SARS-CoV-2 transmission and effective interventions in a higher education setting that will inform public health policy during pandemics.

Significant variability exists in the cytotoxicity of global methicillin-resistant Staphylococcus aureus lineages.

Staphylococcus aureus is a major human pathogen where the emergence of antibiotic resistant lineages, such as methicillin-resistant S. aureus (MRSA), is a major health concern. While some MRSA lineages are restricted to the healthcare setting, the epidemiology of MRSA is changing globally, with the rise of specific lineages causing disease in healthy people in the community. In the past two decades, community-associated MRSA (CA-MRSA) has emerged as a clinically important and virulent pathogen associated with serious skin and soft-tissue infections (SSTI). These infections are primarily cytotoxin driven, leading to the suggestion that hypervirulent lineages/multi-locus sequence types (STs) exist. To examine this, we compared the cytotoxicity of 475 MRSA isolates representing five major MRSA STs (ST22, ST93, ST8, ST239 and ST36) by employing a monocyte-macrophage THP-1 cell line as a surrogate for measuring gross cytotoxicity. We demonstrate that while certain MRSA STs contain highly toxic isolates, there is such variability within lineages to suggest that this aspect of virulence should not be inferred from the genotype of any given isolate. Furthermore, by interrogating the accessory gene regulator (Agr) sequences in this collection we identified several Agr mutations that were associated with reduced cytotoxicity. Interestingly, the majority of isolates that were attenuated in cytotoxin production contained no mutations in the agr locus, indicating a role of other undefined genes in S. aureus toxin regulation.

Defining nosocomial transmission of Escherichia coli and antimicrobial resistance genes: a genomic surveillance study.

BACKGROUND: Escherichia coli is a leading cause of bloodstream infections. Developing interventions to reduce E coli infections requires an understanding of the frequency of nosocomial transmission, but the available evidence is scarce. We aimed to detect and characterise transmission of E coli and associated plasmids in a hospital setting. METHODS: In this prospective observational cohort study, patients were admitted to two adult haematology wards at the Cambridge University Hospitals NHS Foundation Trust in England. Patients aged 16 years and older who were treated for haematological malignancies were included. Stool samples were collected from study participants on admission, once per week, and at discharge. We sequenced multiple E coli isolates (both extended spectrum ß-lactamase [ESBL]-producing and non-ESBL-producing) from each stool sample. A genetic threshold to infer E coli transmission was defined by maximum within-host single nucleotide polymorphism (SNP) diversity and the probability of drawing observed pairs of between-patient isolates at different SNP thresholds. Putative transmission clusters were identified when sequences were less than the genetic threshold. Epidemiological links for each transmission event were investigated. We sequenced all E coli positive blood samples from the two adult haematology wards. FINDINGS: We recruited 174 (51%) of 338 adult patients admitted to the wards between May 13 and Nov 13, 2015. We obtained and cultured 376 stool samples from 149 patients, of which 152 samples from 97 (65%) patients grew E coli. Whole-genome sequencing was done on 970 isolates. We identified extensive diversity in the bacterial population (90 sequence types) and mixed E coli sequence type carriage. 24 (26%) patients carried two sequence types, 12 (13%) carried three, and six (6%) patients carried four or more sequence types. Using a 17 SNP cutoff we identified ten clusters in 20 patients. The largest cluster contained seven patients, whereas four patients were included in multiple clusters. Strong epidemiological links were found between patients in seven clusters. 17 (11%) of 149 patients had stool samples positive for ESBL-producing E coli, the most common of which was associated with bla CTX-M-15 (12 [71%] of 17). Long-read sequencing revealed that bla CTX-M-15 was often integrated into the chromosome, with little evidence for plasmid transmission. Seven patients developed E coli bloodstream infection, four with identical strains to those in their stool; two of these had documented nosocomial acquisition. INTERPRETATION: We provide evidence of bacterial transmission and endogenous infection during routine care by integrating genomic and epidemiological data and by determining a genetic cutoff informed by within-host diversity in the studied population. Our findings challenge single colony-based investigations, and the widely accepted notion of plasmid spread. FUNDING: UK Department of Health, Wellcome Trust, UK National Institute for Health Research.

Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial.

BACKGROUND: A new variant of SARS-CoV-2, B.1.1.7, emerged as the dominant cause of COVID-19 disease in the UK from November, 2020. We report a post-hoc analysis of the efficacy of the adenoviral vector vaccine, ChAdOx1 nCoV-19 (AZD1222), against this variant. METHODS: Volunteers (aged ≥18 years) who were enrolled in phase 2/3 vaccine efficacy studies in the UK, and who were randomly assigned (1:1) to receive ChAdOx1 nCoV-19 or a meningococcal conjugate control (MenACWY) vaccine, provided upper airway swabs on a weekly basis and also if they developed symptoms of COVID-19 disease (a cough, a fever of 37·8°C or higher, shortness of breath, anosmia, or ageusia). Swabs were tested by nucleic acid amplification test (NAAT) for SARS-CoV-2 and positive samples were sequenced through the COVID-19 Genomics UK consortium. Neutralising antibody responses were measured using a live-virus microneutralisation assay against the B.1.1.7 lineage and a canonical non-B.1.1.7 lineage (Victoria). The efficacy analysis included symptomatic COVID-19 in seronegative participants with a NAAT positive swab more than 14 days after a second dose of vaccine. Participants were analysed according to vaccine received. Vaccine efficacy was calculated as 1 - relative risk (ChAdOx1 nCoV-19 vs MenACWY groups) derived from a robust Poisson regression model. This study is continuing and is registered with ClinicalTrials.gov, NCT04400838, and ISRCTN, 15281137. FINDINGS: Participants in efficacy cohorts were recruited between May 31 and Nov 13, 2020, and received booster doses between Aug 3 and Dec 30, 2020. Of 8534 participants in the primary efficacy cohort, 6636 (78%) were aged 18-55 years and 5065 (59%) were female. Between Oct 1, 2020, and Jan 14, 2021, 520 participants developed SARS-CoV-2 infection. 1466 NAAT positive nose and throat swabs were collected from these participants during the trial. Of these, 401 swabs from 311 participants were successfully sequenced. Laboratory virus neutralisation activity by vaccine-induced antibodies was lower against the B.1.1.7 variant than against the Victoria lineage (geometric mean ratio 8·9, 95% CI 7·2-11·0). Clinical vaccine efficacy against symptomatic NAAT positive infection was 70·4% (95% CI 43·6-84·5) for B.1.1.7 and 81·5% (67·9-89·4) for non-B.1.1.7 lineages. INTERPRETATION: ChAdOx1 nCoV-19 showed reduced neutralisation activity against the B.1.1.7 variant compared with a non-B.1.1.7 variant in vitro, but the vaccine showed efficacy against the B.1.1.7 variant of SARS-CoV-2. FUNDING: UK Research and Innovation, National Institute for Health Research (NIHR), Coalition for Epidemic Preparedness Innovations, NIHR Oxford Biomedical Research Centre, Thames Valley and South Midlands NIHR Clinical Research Network, and AstraZeneca.

A common protocol for the simultaneous processing of multiple clinically relevant bacterial species for whole genome sequencing.

Whole-genome sequencing is likely to become increasingly used by local clinical microbiology laboratories, where sequencing volume is low compared with national reference laboratories. Here, we describe a universal protocol for simultaneous DNA extraction and sequencing of numerous different bacterial species, allowing mixed species sequence runs to meet variable laboratory demand. We assembled test panels representing 20 clinically relevant bacterial species. The DNA extraction process used the QIAamp mini DNA kit, to which different combinations of reagents were added. Thereafter, a common protocol was used for library preparation and sequencing. The addition of lysostaphin, lysozyme or buffer ATL (a tissue lysis buffer) alone did not produce sufficient DNA for library preparation across the species tested. By contrast, lysozyme plus lysostaphin produced sufficient DNA across all 20 species. DNA from 15 of 20 species could be extracted from a 24-h culture plate, while the remainder required 48-72 h. The process demonstrated 100% reproducibility. Sequencing of the resulting DNA was used to recapitulate previous findings for species, outbreak detection, antimicrobial resistance gene detection and capsular type. This single protocol for simultaneous processing and sequencing of multiple bacterial species supports low volume and rapid turnaround time by local clinical microbiology laboratories.

Quantifying acquisition and transmission of Enterococcus faecium using genomic surveillance.

Nosocomial acquisition and transmission of vancomycin-resistant Enterococcus faecium (VREfm) is the driver for E. faecium carriage in hospitalized patients, which, in turn, is a risk factor for invasive infection in immunocompromised patients. In the present study, we provide a comprehensive picture of E. faecium transmission in an entire sampled patient population using a sequence-driven approach. We prospectively identified and followed 149 haematology patients admitted to a hospital in England for 6 months. Patient stools (n = 376) and environmental swabs (n = 922) were taken at intervals and cultured for E. faecium. We sequenced 1,560 isolates (1,001 stool, 559 environment) and focused our genomic analyses on 1,477 isolates (95%) in the hospital-adapted clade A1. Of 101 patients who provided two or more stool samples, 40 (40%) developed E. faecium carriage after admission based on culture, compared with 64 patients (63%) based on genomic analysis (73% VREfm). Half of 922 environmental swabs (447, 48%) were positive for VREfm. Network analysis showed that, of 111 patients positive for the A1 clade, 67 had strong epidemiological and genomic links with at least one other patient and/or their direct environment, supporting nosocomial transmission. Six patients (3.4%) developed an invasive E. faecium infection from their own gut-colonizing strain, which was preceded by nosocomial acquisition of the infecting isolate in half of these. Two informatics approaches (subtype categorization to define phylogenetic clusters and the development of an SNP cut-off for transmission) were central to our analyses, both of which will inform the future translation of E. faecium sequencing into routine outbreak detection and investigation. In conclusion, we showed that carriage and environmental contamination by the hospital-adapted E. faecium lineage were hyperendemic in our study population and that improved infection control measures will be needed to reduce hospital acquisition rates.

Definition of a genetic relatedness cutoff to exclude recent transmission of meticillin-resistant Staphylococcus aureus: a genomic epidemiology analysis.

BACKGROUND: Whole-genome sequencing (WGS) can be used in genomic epidemiology investigations to confirm or refute outbreaks of bacterial pathogens, and to support targeted and efficient infection control interventions. We aimed to define a genetic relatedness cutoff, quantified as a number of single-nucleotide polymorphisms (SNP), for meticillin-resistant Staphylococcus aureus (MRSA), above which recent (ie, within 6 months) patient-to-patient transmission could be ruled out. METHODS: We did a retrospective genomic and epidemiological analysis of MRSA data from two prospective observational cohort studies in the UK to establish SNP cutoffs for genetic relatedness, above which recent transmission was unlikely. We used three separate approaches to calculate these thresholds. First, we applied a linear mixed model to estimate the S aureus substitution rate and 95th percentile within-host diversity in a cohort in which multiple isolates were sequenced per individual. Second, we applied a simulated transmission model to this same genomic dataset. Finally, in a second cohort, we determined the genetic distance (ie, the number of SNPs) that would capture 95% of epidemiologically linked cases. We applied the three approaches to both whole-genome and core-genome sequences. FINDINGS: In the linear mixed model, the estimated substitution rate was roughly 5 whole-genome SNPs (wgSNPs) or 3 core-genome SNPs (cgSNPs) per genome per year, and the 95th percentile within-host diversity was 19 wgSNPs or 10 cgSNPs. The combined SNP cutoffs for detection of MRSA transmission within 6 months per this model were thus 24 wgSNPs or 13 cgSNPs. The simulated transmission model suggested that cutoffs of 17 wgSNPs or 12 cgSNPs would detect 95% of MRSA transmission events within the same timeframe. Finally, in the second cohort, cutoffs of 22 wgSNPs or 11 cgSNPs captured 95% of epidemiologically linked cases within 6 months. INTERPRETATION: On the basis of our results, we propose conservative cutoffs of 25 wgSNPs or 15 cgSNPS above which transmission of MRSA within the previous 6 months can be ruled out. These cutoffs could potentially be used as part of a genomic sequencing approach to the management of outbreaks of MRSA in conjunction with traditional epidemiological techniques. FUNDING: UK Department of Health, Wellcome Trust, UK National Institute for Health Research.

Defining metrics for whole-genome sequence analysis of MRSA in clinical practice.

Bacterial sequencing will become increasingly adopted in routine microbiology laboratories. Here, we report the findings of a technical evaluation of almost 800 clinical methicillin-resistant Staphylococcus aureus (MRSA) isolates, in which we sought to define key quality metrics to support MRSA sequencing in clinical practice. We evaluated the accuracy of mapping to a generic reference versus clonal complex (CC)-specific mapping, which is more computationally challenging. Focusing on isolates that were genetically related (<50 single nucleotide polymorphisms (SNPs)) and belonged to prevalent sequence types, concordance between these methods was 99.5 %. We use MRSA MPROS0386 to control for base calling accuracy by the sequencer, and used multiple repeat sequences of the control to define a permitted range of SNPs different to the mapping reference for this control (equating to 3 standard deviations from the mean). Repeat sequences of the control were also used to demonstrate that SNP calling was most accurate across differing coverage depths (above 35×, the lowest depth in our study) when the depth required to call a SNP as present was at least 4-8×. Using 786 MRSA sequences, we defined a robust measure for mec gene detection to reduce false-positives arising from contamination, which was no greater than 2 standard deviations below the average depth of coverage across the genome. Sequencing from bacteria harvested from clinical plates runs an increased risk of contamination with the same or different species, and we defined a cut-off of 30 heterozygous sites >50 bp apart to identify same-species contamination for MRSA. These metrics were combined into a quality-control (QC) flowchart to determine whether sequence runs and individual clinical isolates passed QC, which could be adapted by future automated analysis systems to enable rapid hands-off sequence analysis by clinical laboratories.

Evaluation of a fully automated bioinformatics tool to predict antibiotic resistance from MRSA genomes.

OBJECTIVES: The genetic prediction of phenotypic antibiotic resistance based on analysis of WGS data is becoming increasingly feasible, but a major barrier to its introduction into routine use is the lack of fully automated interpretation tools. Here, we report the findings of a large evaluation of the Next Gen Diagnostics (NGD) automated bioinformatics analysis tool to predict the phenotypic resistance of MRSA. METHODS: MRSA-positive patients were identified in a clinical microbiology laboratory in England between January and November 2018. One MRSA isolate per patient together with all blood culture isolates (total n = 778) were sequenced on the Illumina MiniSeq instrument in batches of 21 clinical MRSA isolates and three controls. RESULTS: The NGD system activated post-sequencing and processed the sequences to determine susceptible/resistant predictions for 11 antibiotics, taking around 11 minutes to analyse 24 isolates sequenced on a single sequencing run. NGD results were compared with phenotypic susceptibility testing performed by the clinical laboratory using the disc diffusion method and EUCAST breakpoints. Following retesting of discrepant results, concordance between phenotypic results and NGD genetic predictions was 99.69%. Further investigation of 22 isolate genomes associated with persistent discrepancies revealed a range of reasons in 12 cases, but no cause could be found for the remainder. Genetic predictions generated by the NGD tool were compared with predictions generated by an independent research-based informatics approach, which demonstrated an overall concordance between the two methods of 99.97%. CONCLUSIONS: We conclude that the NGD system provides rapid and accurate prediction of the antibiotic susceptibility of MRSA.

A One Health Study of the Genetic Relatedness of Klebsiella pneumoniae and Their Mobile Elements in the East of England.

BACKGROUND: Klebsiella pneumoniae is a human, animal, and environmental commensal and a leading cause of nosocomial infections, which are often caused by multiresistant strains. We evaluate putative sources of K. pneumoniae that are carried by and infect hospital patients. METHODS: We conducted a 6-month survey on 2 hematology wards at Addenbrooke's Hospital, Cambridge, United Kingdom, in 2015 to isolate K. pneumoniae from stool, blood, and the environment. We conducted cross-sectional surveys of K. pneumoniae from 29 livestock farms, 97 meat products, the hospital sewer, and 20 municipal wastewater treatment plants in the East of England between 2014 and 2015. Isolates were sequenced and their genomes compared. RESULTS: Klebsiella pneumoniae was isolated from stool of 17/149 (11%) patients and 18/922 swabs of their environment, together with 1 bloodstream infection during the study and 4 others over a 24-month period. Each patient carried 1 or more lineages that was unique to them, but 2 broad environmental contamination events and patient-environment transmission were identified. Klebsiella pneumoniae was isolated from cattle, poultry, hospital sewage, and 12/20 wastewater treatment plants. There was low genetic relatedness between isolates from patients/their hospital environment vs isolates from elsewhere. Identical genes encoding cephalosporin resistance were carried by isolates from humans/environment and elsewhere but were carried on different plasmids. CONCLUSION: We identified no patient-to-patient transmission and no evidence for livestock as a source of K. pneumoniae infecting humans. However, our findings reaffirm the importance of the hospital environment as a source of K. pneumoniae associated with serious human infection.

The Emergence of Successful Streptococcus pyogenes Lineages through Convergent Pathways of Capsule Loss and Recombination Directing High Toxin Expression.

Gene transfer and homologous recombination in Streptococcus pyogenes has the potential to trigger the emergence of pandemic lineages, as exemplified by lineages of emm1 and emm89 that emerged in the 1980s and 2000s, respectively. Although near-identical replacement gene transfer events in the nga (NADase) and slo (streptolysin O) loci conferring high expression of these toxins underpinned the success of these lineages, extension to other emm genotype lineages is unreported. The emergent emm89 lineage was characterized by five regions of homologous recombination additional to nga-slo, including complete loss of the hyaluronic acid capsule synthesis locus hasABC, a genetic trait replicated in two other leading emm types and recapitulated by other emm types by inactivating mutations. We hypothesized that other leading genotypes may have undergone similar recombination events. We analyzed a longitudinal data set of genomes from 344 clinical invasive disease isolates representative of locations across England, dating from 2001 to 2011, and an international collection of S. pyogenes genomes representing 54 different genotypes and found frequent evidence of recombination events at the nga-slo locus predicted to confer higher toxin genotype. We identified multiple associations between recombination at this locus and inactivating mutations within hasAB, suggesting convergent evolutionary pathways in successful genotypes. This included common genotypes emm28 and emm87. The combination of no or low capsule and high expression of nga and slo may underpin the success of many emergent S. pyogenes lineages of different genotypes, triggering new pandemics, and could change the way S. pyogenes causes disease.IMPORTANCEStreptococcus pyogenes is a genetically diverse pathogen, with over 200 different genotypes defined by emm typing, but only a minority of these genotypes are responsible for the majority of human infection in high-income countries. Two prevalent genotypes associated with disease rose to international dominance following recombination of a toxin locus that conferred increased expression. Here, we found that recombination of this locus and promoter has occurred in other diverse genotypes, events that may allow these genotypes to expand in the population. We identified an association between the loss of hyaluronic acid capsule synthesis and high toxin expression, which we propose may be associated with an adaptive advantage. As S. pyogenes pathogenesis depends both on capsule and toxin production, new variants with altered expression may result in abrupt changes in the molecular epidemiology of this pathogen in the human population over time.

The composition and functional protein subsystems of the human nasal microbiome in granulomatosis with polyangiitis: a pilot study.

BACKGROUND: Ear, nose and throat involvement in granulomatosis with polyangiitis (GPA) is frequently the initial disease manifestation. Previous investigations have observed a higher prevalence of Staphylococcus aureus in patients with GPA, and chronic nasal carriage has been linked with an increased risk of disease relapse. In this cross-sectional study, we investigated changes in the nasal microbiota including a detailed analysis of Staphylococcus spp. by shotgun metagenomics in patients with active and inactive granulomatosis with polyangiitis (GPA). Shotgun metagenomic sequence data were also used to identify protein-encoding genes within the SEED database, and the abundance of proteins then correlated with the presence of bacterial species on an annotated heatmap. RESULTS: The presence of S. aureus in the nose as assessed by culture was more frequently detected in patients with active GPA (66.7%) compared with inactive GPA (34.1%). Beta diversity analysis of nasal microbiota by bacterial 16S rRNA profiling revealed a different composition between GPA patients and healthy controls (P = 0.039). Beta diversity analysis of shotgun metagenomic sequence data for Staphylococcus spp. revealed a different composition between active GPA patients and healthy controls and disease controls (P = 0.0007 and P = 0.0023, respectively), and between healthy controls and inactive GPA patients and household controls (P = 0.0168 and P = 0.0168, respectively). Patients with active GPA had a higher abundance of S. aureus, mirroring the culture data, while healthy controls had a higher abundance of S. epidermidis. Staphylococcus pseudintermedius, generally assumed to be a pathogen of cats and dogs, showed an abundance of 13% among the Staphylococcus spp. in our cohort. During long-term follow-up of patients with inactive GPA at baseline, a higher S. aureus abundance was not associated with an increased relapse risk. Functional analyses identified ten SEED protein subsystems that differed between the groups. Most significant associations were related to chorismate synthesis and involved in the vitamin B12 pathway. CONCLUSION: Our data revealed a distinct dysbiosis of the nasal microbiota in GPA patients compared with disease and healthy controls. Metagenomic sequencing demonstrated that this dysbiosis in active GPA patients is manifested by increased abundance of S. aureus and a depletion of S. epidermidis, further demonstrating the antagonist relationships between these species. SEED functional protein subsystem analysis identified an association between the unique bacterial nasal microbiota clusters seen mainly in GPA patients and an elevated abundance of genes associated with chorismate synthesis and vitamin B12 pathways. Further studies are required to further elucidate the relationship between the biosynthesis genes and the associated bacterial species.