Multi-site implementation of whole genome sequencing for hospital infection control: A prospective genomic epidemiological analysis

Genomics for antimicrobial resistance-progress and future directions

Antimicrobial resistance (AMR) is a critical global public health threat, with bacterial pathogens of primary concern. Pathogen genomics has revolutionized the study of bacterial pathogens and provided deep insights into the mechanisms and dissemination of AMR, with the precision of whole-genome sequencing informing better control strategies. However, generating actionable data from genomic surveillance and diagnostic efforts requires integration at the public health and clinical interface that goes beyond academic efforts to identify resistance mechanisms, undertake post hoc analyses of outbreaks, and share data after research publications. In addition to timely genomics data, consideration also needs to be given to epidemiological sampling frames, analysis, and reporting mechanisms that meet International Organization for Standardization (ISO) standards and generation of reports that are interpretable and actionable for public health and clinical "end-users." Importantly, ensuring all countries have equitable access to data and technology is critical, through timely data sharing following the FAIR principles (findable, accessible, interoperable, and re-usable). In this review, we describe (i) advances in genomic approaches for AMR research and surveillance to understand emergence, evolution, and transmission of AMR and the key requirements to enable this work and (ii) discuss emerging and future applications of genomics at the clinical and public health interface, including barriers to implementation. Harnessing advances in genomics-enhanced AMR research and embedding robust and reproducible workflows within clinical and public health practice promises to maximize the impact of pathogen genomics for AMR globally in the coming decade.

How does date-rounding affect phylodynamic inference for public health?

Phylodynamic analyses infer epidemiological parameters from pathogen genome sequences for enhanced genomic surveillance in public health. Pathogen genome sequences and their associated sampling dates are the essential data in every analysis. However, sampling dates are usually associated with hospitalisation or testing and can sometimes be used to identify individual patients, posing a threat to patient confidentiality. To lower this risk, sampling dates are often given with reduced date-resolution to the month or year, which can potentially bias inference. Here, we introduce a practical guideline on when date-rounding biases the inference of epidemiologically important parameters across a diverse range of empirical and simulated datasets. We show that the direction of bias varies for different parameters, datasets, and tree priors, while compounding with lower date-resolution and higher substitution rates. We also find that bias decreases for datasets with longer sampling intervals, implying that our guideline is most applicable to emerging datasets. We conclude by discussing future solutions that prioritise patient confidentiality and propose a method for safer sharing of sampling dates that translates them them uniformly by a random number.

Bacteriocin production facilitates nosocomial emergence of vancomycin-resistant Enterococcus faecium

Gastrointestinal colonization by the nosocomial pathogen vancomycin-resistant Enterococcus faecium (VREfm) can lead to bloodstream infections with high mortality rates. Shifts in VREfm lineages found within healthcare settings occur, but reasons underlying these changes are not understood. Here we sequenced 710 VREfm clinical isolates collected between 2017 and 2022 from a large tertiary care centre. Genomic analyses revealed a polyclonal VREfm population, although 46% of isolates formed genetically related clusters, suggesting a high transmission rate. Comparing these data to a global collection of 15,631 publicly available VREfm genomes collected between 2002 and 2022 identified replacement of the sequence type (ST) 17 VREfm lineage by emergent ST80 and ST117 lineages at the local and global level. Comparative genomic and functional analyses revealed that emergent lineages encoded bacteriocin T8, which conferred a competitive advantage over bacteriocin T8-negative strains in vitro and upon colonization of the mouse gut. Bacteriocin T8 carriage was also strongly associated with strain emergence in the global genome collection. These data suggest that bacteriocin T8-mediated competition may have contributed to VREfm lineage replacement.

Rapid, reference-free identification of bacterial pathogen transmission using optimized split k-mer analysis

Infections caused by multidrug-resistant organisms (MDROs) are difficult to treat and often life threatening and place a burden on the healthcare system. Minimizing the transmission of MDROs in hospitals is a global priority with genomics proving to be a powerful tool for identifying the transmission of MDROs. To optimize the utility of genomics for prospective infection control surveillance, results must be available in real time, reproducible and simple to communicate to clinicians. Traditional reference-based approaches suffer from several limitations for prospective genomic surveillance. Whilst reference-free or pairwise genome comparisons avoid some of these limitations, they can be computationally intensive and time consuming. Split k-mer analysis (SKA) offers a viable alternative facilitating rapid reference-free pairwise comparisons of genomic data, but the optimum SKA parameters for the detection of transmission have not been determined. Additionally, the accuracy of SKA-based inferences has not been measured, nor whether modified quality control parameters are required. Here, we explore the performance of 60 SKA parameter combinations across 50 simulations to quantify the false negative and positive SNP proportions for Escherichia coli, Enterococcus faecium, Klebsiella pneumoniae and Staphylococcus aureus. Using the optimum parameter combination, we explore concordance between SKA, multilocus sequence typing (MLST), core genome MLST (cgMLST) and Snippy in a real-world dataset. Lastly, we investigate whether simulated plasmid gain or loss could impact SNP detection with SKA. This work identifies that the use of SKA with sequencing reads, a k-mer length of 19 and a minor allele frequency filter of 0.01 is optimal for MDRO transmission detection. Whilst SNP detection with SKA (when used with sequencing reads) undercalls SNPs compared to Snippy, it is significantly faster, especially with larger datasets. SKA has excellent concordance with MLST and cgMLST and is not impacted by simulated plasmid movement. We propose that the use of SKA for the detection of bacterial pathogen transmission is superior to traditional methodologies, capable of providing results in a much shorter timeframe.

Pathogen genomics in healthcare: overcoming barriers to proactive surveillance

Pathogen genomic surveillance in healthcare has the potential to enhance patient safety by detecting outbreaks earlier, thereby reducing morbidity and mortality. Despite benefits, there are barriers to adoption, including cost, expertise, and lack of standardized methodologies and incentives. This commentary advocates for 1) investment from healthcare payors, public health, and regulatory bodies and 2) additional research on genomic surveillance for improving patient outcomes and reducing infections. Effective implementation will require strategic investment and cross-sector collaboration.

Evaluation of nationwide analysis surveillance for methicillin-resistant Staphylococcus aureus within Genomic Medicine Sweden

Background. National epidemiological investigations of microbial infections greatly benefit from the increased information gained by whole-genome sequencing (WGS) in combination with standardized approaches for data sharing and analysis.Aim. To evaluate the quality and accuracy of WGS data generated by different laboratories but analysed by joint pipelines to reach a national surveillance approach.Methods. A national methicillin-resistant Staphylococcus aureus (MRSA) collection of 20 strains was distributed to nine participating laboratories that performed in-house procedures for WGS. Raw data were shared and analysed by three pipelines: 1928 Diagnostics, JASEN (GMS pipeline) and CLC-Genomics Workbench. The outcomes were compared according to quality, correct strain identification and genetic distances.Results. One isolate contained intraspecies contamination and was excluded from further analysis. The mean sequencing depth varied between sites and technologies. However, all analysis methods identified 12 strains that belonged to one of five outbreak clusters. The cut-off definition was set to <10 allele differences for core genome multilocus sequence typing (cgMLST) and <20 genetic differences for SNP analysis in a pairwise comparison.Conclusions. MRSA isolates, which are whole genome sequenced by different laboratories and analysed using the same bioinformatic pipelines, yielded comparable results for outbreak clustering for both cgMLST and SNP, using the 1928 analysis pipeline. In this study, JASEN was best suited to analyse Illumina data and CLC to analyse within respective technology. In the future, real-time sharing of data and harmonized analysis within the Genomic Medicine Sweden consortium will further facilitate investigations of outbreaks and transmission routes.

How long do bacteria, fungi, protozoa, and viruses retain their replication capacity on inanimate surfaces? A systematic review examining environmental resilience versus healthcare-associated infection risk by "fomite-borne risk assessment"

SUMMARYIn healthcare settings, contaminated surfaces play an important role in the transmission of nosocomial pathogens potentially resulting in healthcare-associated infections (HAI). Pathogens can be transmitted directly from frequent hand-touch surfaces close to patients or indirectly by staff and visitors. HAI risk depends on exposure, extent of contamination, infectious dose (ID), virulence, hygiene practices, and patient vulnerability. This review attempts to close a gap in previous reviews on persistence/tenacity by only including articles (n = 171) providing quantitative data on re-cultivable pathogens from fomites for a better translation into clinical settings. We have therefore introduced the new term "replication capacity" (RC). The RC is affected by the degree of contamination, surface material, temperature, relative humidity, protein load, organic soil, UV-light (sunlight) exposure, and pH value. In general, investigations into surface RC are mainly performed in vitro using reference strains with high inocula. In vitro data from studies on 14 Gram-positive, 26 Gram-negative bacteria, 18 fungi, 4 protozoa, and 37 viruses. It should be regarded as a worst-case scenario indicating the upper bounds of risks when using such data for clinical decision-making. Information on RC after surface contamination could be seen as an opportunity to choose the most appropriate infection prevention and control (IPC) strategies. To help with decision-making, pathogens characterized by an increased nosocomial risk for transmission from inanimate surfaces ("fomite-borne") are presented and discussed in this systematic review. Thus, the review offers a theoretical basis to support local risk assessments and IPC recommendations.

More than just the gene: investigating expression using a non-native plasmid and host and its impact on resistance conferred by β-lactamase OXA-58 isolated from a hospital wastewater microbiome

With the escalation of hospital-acquired infections by multidrug resistant bacteria, understanding antibiotic resistance is of paramount importance. This study focuses on the β-lactamase gene, blaOXA-58, an important resistance determinant identified in a patient-facing hospital wastewater system. This study aimed to characterize the behaviour of the OXA-58 enzyme when expressed using a non-native plasmid and expression host. blaOXA-58 was cloned using a pET28a(+)/Escherichia coli BL21(DE3) expression system. Nitrocefin hydrolysis and antimicrobial susceptibility of OXA-58-producing cells were assessed against penicillin G, ampicillin, meropenem, and amoxicillin. blaOXA-58 conferred resistance to amoxicillin, penicillin G, and ampicillin, but not to meropenem. This was unexpected given OXA-58's annotation as a carbapenemase. The presence of meropenem also reduced nitrocefin hydrolysis, suggesting it acts as a competitive inhibitor of the OXA-58 enzyme. This study elucidates the phenotypic resistance conferred by an antimicrobial resistance gene (ARG) obtained from a clinically relevant setting and reveals that successful functional expression of ARGs is multifaceted. This study challenges the reliability of predicting antimicrobial resistance based solely on gene sequence alone, and serves as a reminder of the intricate interplay between genetics and structural factors in understanding resistance profiles across different host environments.

One decade of point-prevalence surveys for carriage of extended-spectrum beta-lactamase-producing enterobacterales: whole genome sequencing based prevalence and genetic characterization in a large Dutch teaching hospital from 2013 to 2022

OBJECTIVES

To determine the prevalence, trends, and potential nosocomial transmission events of the hidden reservoir of rectal carriage of extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-E).

METHODS

From 2013 to 2022, yearly point prevalence surveys were conducted in a large Dutch teaching hospital. On the day of the survey, all admitted patients were screened for ESBL-E rectal carriage using peri-anal swabs and a consistent and sensitive selective culturing method. All Enterobacterales phenotypically suspected of ESBL production were analysed using whole genome sequencing for ESBL gene detection and clonal relatedness analysis.

RESULTS

On average, the ESBL-E prevalence was 4.6% (188/4,119 patients), ranging from 2.1 to 6.6% per year. The ESBL-prevalence decreased on average 5.5% per year. After time trend correction, the prevalence in 2016 and 2020 was lower compared to the other year. Among the ESBL-E, Escherichia coli (80%) and CTX-M genes (85%) predominated. Potential nosocomial transmission events could be found in 5.9% (11/188) of the ESBL-E carriers.

CONCLUSIONS

The ESBL-E rectal carriage prevalence among hospitalized patients was 4.6% with a downward trend from 2013 to 2022. The decrease in ESBL-E prevalence in 2020 could have been due to the COVID-19 pandemic and subsequent countrywide measures as no nosocomial transmission events were detected in 2020. However, the persistently low ESBL-E prevalences in 2021 and 2022 suggest that the decline in ESBL-E prevalence goes beyond the COVID-19 pandemic, indicating that overall ESBL-E carriage rates are declining over time. Continuous monitoring of ESBL-E prevalence and transmission rates can aid infection control policy to keep antibiotic resistance rates in hospitals low.

Bacteriocin production facilitates nosocomial emergence of vancomycin-resistant Enterococcus faecium

Vancomycin-resistant Enterococcus faecium (VREfm) is a prevalent healthcare-acquired pathogen. Gastrointestinal colonization can lead to difficult-to-treat bloodstream infections with high mortality rates. Prior studies have investigated VREfm population structure within healthcare centers. However, little is known about how and why hospital-adapted VREfm populations change over time. We sequenced 710 healthcare-associated VREfm clinical isolates from 2017-2022 from a large tertiary care center as part of the Enhanced Detection System for Healthcare-Associated Transmission (EDS-HAT) program. Although the VREfm population in our center was polyclonal, 46% of isolates formed genetically related clusters, suggesting a high transmission rate. We compared our collection to 15,631 publicly available VREfm genomes spanning 20 years. Our findings describe a drastic shift in lineage replacement within nosocomial VREfm populations at both the local and global level. Functional and genomic analysis revealed, antimicrobial peptide, bacteriocin T8 may be a driving feature of strain emergence and persistence in the hospital setting.

Genomic sequencing surveillance of patients colonized with vancomycin-resistant Enterococcus (VRE) improves detection of hospital-associated transmission

Background

Vancomycin-resistant enterococcal (VRE) infections pose significant challenges in healthcare. Transmission dynamics of VRE are complex, often involving patient colonization and subsequent transmission through various healthcare-associated vectors. We utilized a whole genome sequencing (WGS) surveillance program at our institution to better understand the contribution of clinical and colonizing isolates to VRE transmission.

Methods

We performed whole genome sequencing on 352 VRE clinical isolates collected over 34 months and 891 rectal screening isolates collected over a 9-month nested period, and used single nucleotide polymorphisms to assess relatedness. We then performed a geo-temporal transmission analysis considering both clinical and rectal screening isolates compared with clinical isolates alone, and calculated 30-day outcomes of patients.

Results

VRE rectal carriage constituted 87.3% of VRE acquisition, with an average monthly acquisition rate of 7.6 per 1000 patient days. We identified 185 genetically related clusters containing 2-42 isolates and encompassing 69.6% of all isolates in the dataset. The inclusion of rectal swab isolates increased the detection of clinical isolate clusters (from 53% to 67%, P<0.01). Geo-temporal analysis identified hotspot locations of VRE transmission. Patients with clinical VRE isolates that were closely related to previously sampled rectal swab isolates experienced 30-day ICU admission (17.5%), hospital readmission (9.2%), and death (13.3%).

Conclusions

Our findings describe the high burden of VRE transmission at our hospital and shed light on the importance of using WGS surveillance of both clinical and rectal screening isolates to better understand the transmission of this pathogen. This study highlights the potential utility of incorporating WGS surveillance of VRE into routine hospital practice for improving infection prevention and patient safety.

WITHDRAWN: Evaluation of Fourier Transform Infrared spectroscopy (IR Biotyper) as a complement to Whole genome sequencing (WGS) to characterise Enterobacter cloacae , Citrobacter freundii and Klebsiella pneumoniae isolates recovered from hospital sinks

The authors have withdrawn their manuscript due to becoming aware of methodology issues related to the curation of the training set used to determine cut-off values for Biotyper cluster assignation and lack of replicate measurements on different days for the isolates analysed. It is therefore unclear whether the conclusions of the manuscript are founded and no further work is possible to correct these issues as the instrument is no longer available to the authors. If you have any questions, please contact the corresponding author.

Consideration of within-patient diversity highlights transmission pathways and antimicrobial resistance gene variability in vancomycin-resistant Enterococcus faecium

BACKGROUND

WGS is increasingly being applied to healthcare-associated vancomycin-resistant Enterococcus faecium (VREfm) outbreaks. Within-patient diversity could complicate transmission resolution if single colonies are sequenced from identified cases.

OBJECTIVES

Determine the impact of within-patient diversity on transmission resolution of VREfm.

MATERIALS AND METHODS

Fourteen colonies were collected from VREfm positive rectal screens, single colonies were collected from clinical samples and Illumina WGS was performed. Two isolates were selected for Oxford Nanopore sequencing and hybrid genome assembly to generate lineage-specific reference genomes. Mapping to closely related references was used to identify genetic variations and closely related genomes. A transmission network was inferred for the entire genome set using Phyloscanner.

RESULTS AND DISCUSSION

In total, 229 isolates from 11 patients were sequenced. Carriage of two or three sequence types was detected in 27% of patients. Presence of antimicrobial resistance genes and plasmids was variable within genomes from the same patient and sequence type. We identified two dominant sequence types (ST80 and ST1424), with two putative transmission clusters of two patients within ST80, and a single cluster of six patients within ST1424. We found transmission resolution was impaired using fewer than 14 colonies.

CONCLUSIONS

Patients can carry multiple sequence types of VREfm, and even within related lineages the presence of mobile genetic elements and antimicrobial resistance genes can vary. VREfm within-patient diversity could be considered in future to aid accurate resolution of transmission networks.

Whole Genome Sequencing: Applications in Clinical Bacteriology

The success in determining the whole genome sequence of a bacterial pathogen was first achieved in 1995 by determining the complete nucleotide sequence of Haemophilus influenzae Rd using the chain-termination method established by Sanger et al. in 1977 and automated by Hood et al. in 1987. However, this technology was laborious, costly, and time-consuming. Since 2004, high-throughput next-generation sequencing technologies have been developed, which are highly efficient, require less time, and are cost-effective for whole genome sequencing (WGS) of all organisms, including bacterial pathogens. In recent years, the data obtained using WGS technologies coupled with bioinformatics analyses of the sequenced genomes have been projected to revolutionize clinical bacteriology. WGS technologies have been used in the identification of bacterial species, strains, and genotypes from cultured organisms and directly from clinical specimens. WGS has also helped in determining resistance to antibiotics by the detection of antimicrobial resistance genes and point mutations. Furthermore, WGS data have helped in the epidemiological tracking and surveillance of pathogenic bacteria in healthcare settings as well as in communities. This review focuses on the applications of WGS in clinical bacteriology.

Economic and health impact modelling of a whole genome sequencing-led intervention strategy for bacterial healthcare-associated infections for England and for the USA

Bacterial healthcare-associated infections (HAIs) are a substantial source of global morbidity and mortality. The estimated cost associated with HAIs ranges from $35 to $45 billion in the USA alone. The costs and accessibility of whole genome sequencing (WGS) of bacteria and the lack of sufficiently accurate, high-resolution, scalable and accessible analysis for strain identification are being addressed. Thus, it is timely to determine the economic viability and impact of routine diagnostic bacterial genomics. The aim of this study was to model the economic impact of a WGS surveillance system that proactively detects and directs interventions for nosocomial infections and outbreaks compared to the current standard of care, without WGS. Using a synthesis of published models, inputs from national statistics, and peer-reviewed articles, the economic impacts of conducting a WGS-led surveillance system addressing the 11 most common nosocomial pathogen groups in England and the USA were modelled. This was followed by a series of sensitivity analyses. England was used to establish the baseline model because of the greater availability of underpinning data, and this was then modified using USA-specific parameters where available. The model for the NHS in England shows bacterial HAIs currently cost the NHS around £3 billion. WGS-based surveillance delivery is predicted to cost £61.1 million associated with the prevention of 74 408 HAIs and 1257 deaths. The net cost saving was £478.3 million, of which £65.8 million were from directly incurred savings (antibiotics, consumables, etc.) and £412.5 million from opportunity cost savings due to re-allocation of hospital beds and healthcare professionals. The USA model indicates that the bacterial HAI care baseline costs are around $18.3 billion. WGS surveillance costs $169.2 million, and resulted in a net saving of ca.$3.2 billion, while preventing 169 260 HAIs and 4862 deaths. From a 'return on investment' perspective, the model predicts a return to the hospitals of £7.83 per £1 invested in diagnostic WGS in the UK, and US$18.74 per $1 in the USA. Sensitivity analyses show that substantial savings are retained when inputs to the model are varied within a wide range of upper and lower limits. Modelling a proactive WGS system addressing HAI pathogens shows significant improvement in morbidity and mortality while simultaneously achieving substantial savings to healthcare facilities that more than offset the cost of implementing diagnostic genomics surveillance.