Genomic surveillance of antimicrobial resistant bacterial colonisation and infection in intensive care patients

Genomics for antimicrobial resistance surveillance to support infection prevention and control in health-care facilities

Integration of genomic technologies into routine antimicrobial resistance (AMR) surveillance in health-care facilities has the potential to generate rapid, actionable information for patient management and inform infection prevention and control measures in near real time. However, substantial challenges limit the implementation of genomics for AMR surveillance in clinical settings. Through a workshop series and online consultation, international experts from across the AMR and pathogen genomics fields convened to review the evidence base underpinning the use of genomics for AMR surveillance in a range of settings. Here, we summarise the identified challenges and potential benefits of genomic AMR surveillance in health-care settings, and outline the recommendations of the working group to realise this potential. These recommendations include the definition of viable and cost-effective use cases for genomic AMR surveillance, strengthening training competencies (particularly in bioinformatics), and building capacity at local, national, and regional levels using hub and spoke models.

Steering and controlling evolution - from bioengineering to fighting pathogens

Control interventions steer the evolution of molecules, viruses, microorganisms or other cells towards a desired outcome. Applications range from engineering biomolecules and synthetic organisms to drug, therapy and vaccine design against pathogens and cancer. In all these instances, a control system alters the eco-evolutionary trajectory of a target system, inducing new functions or suppressing escape evolution. Here, we synthesize the objectives, mechanisms and dynamics of eco-evolutionary control in different biological systems. We discuss how the control system learns and processes information about the target system by sensing or measuring, through adaptive evolution or computational prediction of future trajectories. This information flow distinguishes pre-emptive control strategies by humans from feedback control in biotic systems. We establish a cost-benefit calculus to gauge and optimize control protocols, highlighting the fundamental link between predictability of evolution and efficacy of pre-emptive control.

Community carriage of ESBL-producing Escherichia coli and Klebsiella pneumoniae: a cross-sectional study of risk factors and comparative genomics of carriage and clinical isolates

The global prevalence of infections caused by extended-spectrum β-lactamase-producing Enterobacterales (ESBL-E) is increasing, and for Escherichia coli, observations indicate that this is partly driven by community-onset cases. The ESBL-E population structure in the community is scarcely described, and data on risk factors for carriage are conflicting. Here, we report the prevalence and population structure of fecal ESBL-producing E. coli and Klebsiella pneumoniae (ESBL-Ec/Kp) in a general adult population, examine risk factors, and compare carriage isolates with contemporary clinical isolates. Fecal samples obtained from 4,999 participants (54% women) ≥40 years in the seventh survey of the population-based Tromsø Study, Norway (2015, 2016), were screened for ESBL-Ec/Kp. In addition, we included 118 ESBL-Ec clinical isolates from the Norwegian surveillance program in 2014. All isolates were whole-genome sequenced. Risk factors associated with carriage were analyzed using multivariable logistic regression. ESBL-Ec gastrointestinal carriage prevalence was 3.3% [95% confidence interval (CI) 2.8%-3.9%, no sex difference] and 0.08% (0.02%-0.20%) for ESBL-Kp. For ESBL-Ec, travel to Asia was the only independent risk factor (adjusted odds ratio 3.46, 95% CI 2.18-5.49). E. coli ST131 was most prevalent in both collections. However, the ST131 proportion was significantly lower in carriage (24%) versus clinical isolates (58%, P < 0.001). Carriage isolates were genetically more diverse with a higher proportion of phylogroup A (26%) than clinical isolates (5%, P < 0.001), indicating that ESBL gene acquisition occurs in a variety of E. coli lineages colonizing the gut. STs commonly related to extraintestinal infections were more frequent in clinical isolates also carrying a higher prevalence of antimicrobial resistance, which could indicate clone-associated pathogenicity.IMPORTANCEESBL-Ec and ESBL-Kp are major pathogens in the global burden of antimicrobial resistance. However, there is a gap in knowledge concerning the bacterial population structure of human ESBL-Ec/Kp carriage isolates in the community. We have examined ESBL-Ec/Kp isolates from a population-based study and compared these to contemporary clinical isolates. The large genetic diversity of carriage isolates indicates frequent ESBL gene acquisition, while those causing invasive infections are more clone dependent and associated with a higher prevalence of antibiotic resistance. The knowledge of factors associated with ESBL carriage helps to identify patients at risk to combat the spread of resistant bacteria within the healthcare system. Particularly, previous travel to Asia stands out as a major risk factor for carriage and should be considered in selecting empirical antibiotic treatment in critically ill patients.

Extended-Spectrum β-Lactamase Genes Traverse the Escherichia coli Populations of Intensive Care Unit Patients, Staff, and Environment

Over a 3-month period, we monitored the population of extended-spectrum β-lactam-resistant Escherichia coli (ESBL-EC) associated with the patients, staff, and environment of an intensive care unit (ICU) in Guangzhou, China. Thirty-four clinical isolates were obtained from the same hospital 12 months later. A total of 165 isolates were characterized and whole-genome sequenced, with 24 isolates subjected to long-read sequencing. The diverse population included representatives of 59 different sequence types (STs). ICU patient and environmental isolates were largely distinct from staff isolates and clinical isolates. We observed five instances of highly similar isolates (0 to 13 single nucleotide polymorphisms [SNPs]) being obtained from different patients or bed unit environments. ESBL resistance in this collection was largely conferred by blaCTX-M genes, which were found in 96.4% of all isolates. The contexts of blaCTX-M genes were diverse, situated in multiple chromosomal positions and in various plasmids. We identified blaCTX-M-bearing plasmid lineages that were present in multiple STs across the surveillance, staff, and clinical collections. Closer examination of ISEcp1-blaCTX-M transposition units shed light on the dynamics of their transmission, with evidence for the acquisition of chromosomal copies of blaCTX-M genes from specific plasmid lineages and for the movement of blaCTX-M-55 from a ST1193 chromosome to a small mobilizable plasmid. A carbapenem-resistant ST167 strain isolated from a patient that had been treated with meropenem and piperacillin-tazobactam contained seven copies of blaCMY-146, which appears to have been amplified by IS1. Our data revealed limited persistence and movement of ESBL-EC strains in the ICU environment, but we observed circulating plasmid lineages playing an essential and ongoing role in shaping the cephalosporin-resistance landscape in the population examined. IMPORTANCE ESBL resistance significantly impacts clinical management of E. coli infections in hospitals globally. It is important to understand the structures of ESBL-EC populations carried by hospital patients and staff, their capacity to persist in hospital environments, and the dynamics of mobile genes that drive the spread of ESBL resistance. In our 3-month study, ESBL-EC strains found in the ICU environment were strongly associated with patient carriage but distinct from strains found in staff. However, plasmid lineages carrying blaCTX-M genes were found across the ICU populations and in a collection of clinical isolates obtained 1 year later. By examining their content and contexts, we have traced the recent histories of chromosomal and plasmid-borne ISEcp1-blaCTX-M transposition units in the ICU population. This information allowed us to implicate specific plasmid lineages in the acquisition of chromosomal blaCTX-M genes, even when the plasmids were no longer present, and to detect recent transposition of blaCTX-M-55 from a chromosome to a mobilizable plasmid. Similar high-resolution approaches to the study of mobile genetic elements will be essential if the transmission routes associated with the spread of ESBL resistance are to be understood and subjected to interventions.

Assembling the perfect bacterial genome using Oxford Nanopore and Illumina sequencing

A perfect bacterial genome assembly is one where the assembled sequence is an exact match for the organism's genome-each replicon sequence is complete and contains no errors. While this has been difficult to achieve in the past, improvements in long-read sequencing, assemblers, and polishers have brought perfect assemblies within reach. Here, we describe our recommended approach for assembling a bacterial genome to perfection using a combination of Oxford Nanopore Technologies long reads and Illumina short reads: Trycycler long-read assembly, Medaka long-read polishing, Polypolish short-read polishing, followed by other short-read polishing tools and manual curation. We also discuss potential pitfalls one might encounter when assembling challenging genomes, and we provide an online tutorial with sample data (github.com/rrwick/perfect-bacterial-genome-tutorial).

Associated factors for bacterial colonization in patients admitted to the intensive care unit of the Clinical Hospital of Infectious Diseases

Introduction

This study aimed to identify isolates from colonization and assess the risk factors for bacterial colonization and the risk of death in patients admitted to the intensive care unit (ICU) of the Constanţa County Infectious Diseases Hospital between September 2017 and September 2019.

Methods

This was a retrospective case-control study in a single center that included all patients admitted to the ICU in Constanţa, Romania, who underwent bacteriological screening upon admission and 7 days after admission, between September 2017 and September 2019. In total, 253 patients were included in this study. The nasal exudate, pharyngeal exudate, and rectal swab samples were screened.

Results

In this study, 253 patients were screened bacteriologically, of which 53 had bacterial colonization and 200 did not. Among the bacterial strains, Klebsiella spp. (43.39%) was the most frequently isolated. The predominant resistance mechanism detected in the bacterial isolates was extended-spectrum β-lactamase (ESBL). Multivariate analysis identified a Carmeli score of 3 as an independent risk factor for acquiring bacterial colonization in the ICU. The mortality rate of patients with bacterial colonization was 11.32% and 6% for the patients without colonization (p>0.05).

Conclusions

Our study revealed an increased prevalence of Enterobacterales colonization in the ICU. Risk factors for acquiring bacterial colonization differed depending on the type of bacterial colonization, such as ESBL, carbapenemases, methicillin-resistant Staphylococcus aureus (MRSA), and vancomycin-resistant enterococci (VRE). An independent risk factor for acquiring bacterial colonization was the Carmeli score of 3.

Nanopore-only assemblies for genomic surveillance of the global priority drug-resistant pathogen, Klebsiella pneumoniae

Oxford Nanopore Technologies (ONT) sequencing has rich potential for genomic epidemiology and public health investigations of bacterial pathogens, particularly in low-resource settings and at the point of care, due to its portability and affordability. However, low base-call accuracy has limited the reliability of ONT data for critical tasks such as antimicrobial resistance (AMR) and virulence gene detection and typing, serotype prediction, and cluster identification. Thus, Illumina sequencing remains the standard for genomic surveillance despite higher capital and running costs. We tested the accuracy of ONT-only assemblies for common applied bacterial genomics tasks (genotyping and cluster detection, implemented via Kleborate, Kaptive and Pathogenwatch), using data from 54 unique Klebsiella pneumoniae isolates. ONT reads generated via MinION with R9.4.1 flowcells were basecalled using three alternative models [Fast, High-accuracy (HAC) and Super-accuracy (SUP), available within ONT's Guppy software], assembled with Flye and polished using Medaka. Accuracy of typing using ONT-only assemblies was compared with that of Illumina-only and hybrid ONT+Illumina assemblies, constructed from the same isolates as reference standards. The most resource-intensive ONT-assembly approach (SUP basecalling, with or without Medaka polishing) performed best, yielding reliable capsule (K) type calls for all strains (100 % exact or best matching locus), reliable multi-locus sequence type (MLST) assignment (98.3 % exact match or single-locus variants), and good detection of acquired AMR genes and mutations (88-100 % correct identification across the various drug classes). Distance-based trees generated from SUP+Medaka assemblies accurately reflected overall genetic relationships between isolates. The definition of outbreak clusters from ONT-only assemblies was problematic due to inflation of SNP counts by high base-call errors. However, ONT data could be reliably used to 'rule out' isolates of distinct lineages from suspected transmission clusters. HAC basecalling + Medaka polishing performed similarly to SUP basecalling without polishing. Therefore, we recommend investing compute resources into basecalling (SUP model), wherever compute resources and time allow, and note that polishing is also worthwhile for improved performance. Overall, our results show that MLST, K type and AMR determinants can be reliably identified with ONT-only R9.4.1 flowcell data. However, cluster detection remains challenging with this technology.

Research needs for optimising wastewater-based epidemiology monitoring for public health protection

Wastewater-based epidemiology (WBE) is an unobtrusive method used to observe patterns in illicit drug use, poliovirus, and severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). The pandemic and need for surveillance measures have led to the rapid acceleration of WBE research and development globally. With the infrastructure available to monitor SARS-CoV-2 from wastewater in 58 countries globally, there is potential to expand targets and applications for public health protection, such as other viral pathogens, antimicrobial resistance (AMR), pharmaceutical consumption, or exposure to chemical pollutants. Some applications have been explored in academic research but are not used to inform public health decision-making. We reflect on the current knowledge of WBE for these applications and identify barriers and opportunities for expanding beyond SARS-CoV-2. This paper critically reviews the applications of WBE for public health and identifies the important research gaps for WBE to be a useful tool in public health. It considers possible uses for pathogenic viruses, AMR, and chemicals. It summarises the current evidence on the following: (1) the presence of markers in stool and urine; (2) environmental factors influencing persistence of markers in wastewater; (3) methods for sample collection and storage; (4) prospective methods for detection and quantification; (5) reducing uncertainties; and (6) further considerations for public health use.

ESBL plasmids in Klebsiella pneumoniae: diversity, transmission and contribution to infection burden in the hospital setting

Background

Resistance to third-generation cephalosporins, often mediated by extended-spectrum beta-lactamases (ESBLs), is a considerable issue in hospital-associated infections as few drugs remain for treatment. ESBL genes are often located on large plasmids that transfer horizontally between strains and species of Enterobacteriaceae and frequently confer resistance to additional drug classes. Whilst plasmid transmission is recognised to occur in the hospital setting, the frequency and impact of plasmid transmission on infection burden, compared to ESBL + strain transmission, is not well understood.

Methods

We sequenced the genomes of clinical and carriage isolates of Klebsiella pneumoniae species complex from a year-long hospital surveillance study to investigate ESBL burden and plasmid transmission in an Australian hospital. Long-term persistence of a key transmitted ESBL + plasmid was investigated via sequencing of ceftriaxone-resistant isolates during 4 years of follow-up, beginning 3 years after the initial study.

Results

We found 25 distinct ESBL plasmids. We identified one plasmid, which we called Plasmid A, that carried bla_CTX-M-15 in an IncF backbone similar to pKPN-307. Plasmid A was transmitted at least four times into different Klebsiella species/lineages and was responsible for half of all ESBL episodes during the initial 1-year study period. Three of the Plasmid A-positive strains persisted locally 3–6 years later, and Plasmid A was detected in two additional strain backgrounds. Overall Plasmid A accounted for 21% of ESBL + infections in the follow-up period.

Conclusions

Here, we systematically surveyed ESBL strain and plasmid transmission over 1 year in a single hospital network. Whilst ESBL plasmid transmission events were rare in this setting, they had a significant and sustained impact on the burden of ceftriaxone-resistant and multidrug-resistant infections. If onward transmission of Plasmid A-carrying strains could have been prevented, this may have reduced the number of opportunities for Plasmid A to transmit and create novel ESBL + strains, as well as reducing overall ESBL infection burden.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01103-0.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Molecular epidemiology and antimicrobial resistance phenotype of paediatric bloodstream infections caused by Gram-negative bacteria

Background

Gram-negative organisms are common causes of bloodstream infection (BSI) during the neonatal period and early childhood. Whilst several large studies have characterised these isolates in adults, equivalent data (particularly incorporating whole genome sequencing) is lacking in the paediatric population.

Methods

We perform an epidemiological and sequencing based analysis of Gram-negative bloodstream infections (327 isolates (296 successfully sequenced) from 287 patients) in children <18 years old between 2008 and 2018 in Oxfordshire, UK.

Results

Here we show that the burden of infection lies predominantly in neonates and that most infections are caused by Escherichia coli, Klebsiella spp. and Enterobacter hormaechei. There is no evidence in our setting that the proportion of antimicrobial resistant isolates is increasing in the paediatric population although we identify some evidence of sub-breakpoint increases in gentamicin resistance. The population structure of E. coli BSI isolates in neonates and children mirrors that in adults with a predominance of STs 131/95/73/69 and the same proportions of O-antigen serotypes. In most cases in our setting there is no evidence of transmission/point-source acquisition and we demonstrate the utility of whole genome sequencing to refute a previously suspected outbreak.

Conclusions

Our findings support continued use of current empirical treatment guidelines and suggest that O-antigen targeted vaccines may have a role in reducing the incidence of neonatal sepsis.

A vanA vancomycin-resistant Enterococcus faecium ST80 outbreak resulting from a single importation event

BACKGROUND

A marked genotype shift among vancomycin-resistant Enterococcus faecium (VREfm) from vanB to vanA in Australia between 2011 and 2015 is a well-known phenomenon. It is hypothesized that this was caused by multiple independent clones emerging simultaneously in different settings and/or regions.

OBJECTIVES

To gain insights into the circumstances surrounding the shift from vanB to vanA VREfm in one Australian hospital.

METHODS

The genomes of 69 vanA VREfm isolates from St George Hospital collected between 2009 and 2018 were studied. An expansion of ST80 vanA VREfm was noted following a single introduction. ST80 isolates were thus further characterized using hybrid sequencing and contextualized through comparisons with other published Australian ST80 isolates. Phylogenies were constructed with plasmid sequences compared with the index isolate.

RESULTS

The 2011 expansion of ST80 vanA VREfm isolates in our institution originated from the 2009 index isolate, from a patient transferred from overseas. Phylogenetic analysis with other Australian ST80 vanA VREfm isolates showed that the 2011 expansion event was unique, with limited spread to adjacent local health districts. Plasmid analysis showed multiple variants, which can also be traced back to the 2009 isolate, consistent with ongoing plasmid adaptation over time.

CONCLUSIONS

These findings confirm an expansion event following a VREfm introduction event leading to a sustained clonal and plasmid outbreak over several years. Moreover, it demonstrates the complexity of countrywide replacement events. This study also highlights the use of hybrid sequencing in establishing an epidemiological relationship to the index isolate that was initially inapparent.

Recovery of small plasmid sequences via Oxford Nanopore sequencing

Oxford Nanopore Technologies (ONT) sequencing platforms currently offer two approaches to whole-genome native-DNA library preparation: ligation and rapid. In this study, we compared these two approaches for bacterial whole-genome sequencing, with a specific aim of assessing their ability to recover small plasmid sequences. To do so, we sequenced DNA from seven plasmid-rich bacterial isolates in three different ways: ONT ligation, ONT rapid and Illumina. Using the Illumina read depths to approximate true plasmid abundance, we found that small plasmids (<20 kbp) were underrepresented in ONT ligation read sets (by a mean factor of ~4) but were not underrepresented in ONT rapid read sets. This effect correlated with plasmid size, with the smallest plasmids being the most underrepresented in ONT ligation read sets. We also found lower rates of chimaeric reads in the rapid read sets relative to ligation read sets. These results show that when small plasmid recovery is important, ONT rapid library preparations are preferable to ligation-based protocols.