You are viewing the site in preview mode

Skip to main content

Tracing the transmission of carbapenem-resistant Enterobacterales at the patient: ward environmental nexus

Abstract

Introduction

Colonisation and infection with Carbapenem-resistant Enterobacterales (CRE) in healthcare settings poses significant risks, especially for vulnerable patients. Genomic analysis can be used to trace transmission routes, supporting antimicrobial stewardship and informing infection control strategies. Here we used genomic analysis to track the movement and transmission of CREs within clinical and environmental samples.

Methods

25 isolates were cultured from clinical patient samples or swabs, that tested positive for OXA-48-like variants using the NG-Test® CARBA-5 test and whole genome sequenced (WGS) using Oxford Nanopore Technologies (ONT). 158 swabs and 52 wastewater samples were collected from the ward environment. 60 isolates (matching clinical isolate genera; Klebsiella, Enterobacter, Citrobacter and Escherichia) were isolated from the environmental samples using selective agar. Metagenomic sequencing was undertaken on 36 environmental wastewater and swab samples.

Results

21/25 (84%) clinical isolates had > 1 blaOXA gene and 19/25 (76%) harboured > 1 blaNDM gene. Enterobacterales were most commonly isolated from environmental wastewater samples 27/52 (51.9%), then stick swabs 5/43 (11.6%) and sponge swabs 5/115 (4.3%). 11/60 (18%) environmental isolates harboured > 1 blaOXA gene and 1.9% (1/60) harboured blaNDM-1. blaOXA genes were found in 2/36 (5.5%) metagenomic environmental samples.

Conclusions

Potential for putative patient-patient and patient-ward transmission was shown. Metagenomic sampling needs optimization to improve sensitivity.

Introduction

Colonisation and subsequent healthcare associated infection (HCAI) with multi-drug resistant organisms (MDROs) is a concern for vulnerable patient groups, such as the elderly, or immunocompromised, within the hospital setting. Many of these HCAI incidents could be preventable with enhanced infection and control (IPC) measures [1]. The first line therapy for use against MDROs are the carbapenems, and infection with carbapenem-resistant Enterobacterales (CRE) are associated with high patient mortality rates [2]. CREs are difficult to treat as carbapenemase enzymes can hydrolyse almost all β-lactam antibiotics [3].

Oxacillinase-48-type carbapenemases (OXA-48) and New Delhi metallo-β-lactamase (NDM) are common CRE resistance mechanisms, which are now found globally, are highly mobile and no longer confined to the original bacterial species that they were characterised from [2, 4]. CRE resistance is most often transferred between bacterial isolates on mobile genetic elements, including plasmids, and a wide range of plasmid types have been seen in CRE organisms [4]. CREs are often associated with other genes that confer β-lactam resistances, such as blaSHV and blaCTX, as they can be found on the same plasmids, complicating detection [5].

CREs have been found in both community settings and hospital environments around the globe [6]. Whilst colonised patients do not need antibiotic therapy, they still pose a transmission risk, and so both colonisation and infection must be considered when undertaking IPC. CREs have been isolated from high-touch surfaces, such as door handles, medical equipment trollies, as well as bed sheets and rails, and also from hospital wastewater, including sink U-bends [2, 7,8,9,10]. The colonisation of hospital wastewater by CREs and other MDROs may be particularly problematic, as this water will pass out of the hospital and into the general wastewater system. Without adequate processing at wastewater treatment plants, and the risk of sewage being released into water systems, this may lead to the distribution of these MDROs and their AMR genes back into environmental and community settings [11, 12]. This is especially likely to be problematic in low- and middle-income settings, where the treatment of wastewater may not be adequate to remove these organisms [13].

Surveillance is crucial for the containment of pathogen and AMR outbreaks, especially in hospital settings with vulnerable patients and multiple-occupancy bed bays [14]. The current standard testing for CREs is culture and phenotypic antibiotic sensitivity testing (AST), and molecular methods such as PCR or rapid diagnostic tests (RDTs), such as lateral flow devices (including the NG-Test CARBA-5) [15].

As the COVID-19 pandemic showed with real-time tracking of variants, genomic analysis provides vital enhanced surveillance of transmission patterns [16]. Phenotypic and other molecular tests for resistance genes cannot identify the genetic relatedness of isolates, and thus cannot accurately track potential transmission. The use of whole genome sequencing (WGS) and metagenomics in the diagnosis and surveillance of MDROs can identify phylogeny, novel drug resistance mutations and inform design of targeted diagnostics [17].

blaOXA- and blaNDM-mediated CRE colonisations and infections were detected in several patients occupying single-bed rooms and multiple-occupancy bed bays on wards at a North London tertiary referral hospital between 2022 and 2023. Temporal and spatial associations between patients indicated the possibility of ongoing transmission events but routine phenotypic and molecular testing was unable to pinpoint transmission routes. In this study we aimed to use genomic analysis to build a picture of the movement and transmission of CRE species, AMR and plasmids within clinical pathogens and environmental species found in the ward environment.

Methods

Sample collection

Environmental sampling was requested by the hospital Trust and IPC lead for investigation as part of the extended standard of care. Further characterisation by molecular typing of the CREs that were part of the outbreak was included as an extension of the routine standard of diagnostic care pathway. Patient metadata was obtained through the electronic clinical infection database (elCID). AST profiles of the clinical bacterial pathogens isolated relevant to the CRE surveillance cases were undertaken following the EUCAST Clinical breakpoints (v12.0) for Gram Negative bacteria [18].

Twenty-eight clinical isolates were obtained from 20 patients with a positive NG-Test® CARBA-5 (NG Biotech Laboratories) immunochromatographic lateral flow test for blaOXA and/or blaNDM resistance on multiple wards between 7 February 2022 and 20 January 2023 [19]. Four were from infection sites and 24 were CRE screen samples. These CRE screens were rectal swabs, which were plated onto speciation agar, any suspected to be CREs were tested using the CARBA-5 lateral flow test, as per Health Services Laboratories ‘healthcare associated infection detection of carbapenemase producing organisms by culture’ standard operating procedure. Phenotypic pathogen species data collected from patient samples was obtained using the MALDI-TOF (Bruker). MALDI-TOF-MS of bacterial isolates was undertaken from pure isolates no older than 24 h from culture. Isolates were spotted in duplicate and identifications with corresponding Log scores ≥ 2.0 “high-confidence to the species level” were considered only, and reported in the results. Patient spatio-temporal metadata was collected, see supplementary materials S1. For ward specialties, see supplementary materials S2.

For environmental sampling, a site visit to ward 7D was undertaken prior to collection, to evaluate ward layout and staff and patient routes of travels and sampling locations. Samples were taken from every multiple-occupancy bed bay and single-bed room, covering affected patient areas and non-affected areas. From each bay and room every sink and drain was sampled (in both bed and bathroom areas) and shower drains sampled. Non-clinical rooms were also covered, such as the shared use pantry, staff toilets, nurses station, storage rooms, sluice and workstations on wheels. For the layout of ward 7D see supplementary materials S3.

Two hundred µL of wastewater samples were collected into sterile sample containers pre-dosed with 1 mL of a neutralising buffer comprising: 3% (w/v) Tween 80, 0.3% (w/v) Lecithin, 1.0% (w/v) Sodium thiosulfate, 1.5% (w/v) K2HPO4, KH2PO4 0.05% (w/v), 1% (w/v) Poly-[sodium-4-styrenesulfonate], 0.1% (v/v) Triton® × 100 (Sigma-Aldrich, UK) and prepared in Phosphate-buffered saline (PBS) solution (Oxoid, UK) as previously described [20]. Samples were refrigerated (2–8 °C) within 2 h of collection and processed within 24 h. Aliquots (0.5 mL, 0.1 mL) from the neat and serial 1/10 and 1/100 dilutions from the original wastewater samples were surface-plated onto selective agars: Colorex™ mSuperCARBA™ (EO Labs, UK), Brilliance CRE, E. coli Coliform, Pseudomonas CN (Oxoid, UK) and non-selective Columbia Blood Agar (Oxoid, UK). Plates were incubated aerobically at 37 °C for 72 h and inspected daily. A further 45 µL of the original samples were vacuum-filter concentrated as previously described [21] via 47 mm diameter nitrocellulose membranes (0.45um pore size) and the membrane transferred to selective agars used above prior to incubation.

Cotton-tipped stick swabs (SS352, Appleton Woods) and sponge swabs (TS/15-B, Technical Service Consultants Ltd.), both pre-moistened with a neutraliser buffer were used to collect samples from difficult-to-access areas, as described above. Stick swab samples were transferred to sterile universal tubes containing 9 mL of diluent buffer (saline), 1 mL of neutralising buffer and 3–5 glass beads. The swab contents were released by bead-washing (vortex mixing) for 30 s. Aliquots and serial dilutions of the resulting homogenised suspension was plated on selective and non-selective agars as above.

Suspect colonies were harvested for streak-purification onto non-selective agars and confirmed by MALDI-TOF mass spectrometry (MS) (Maldi-TOF Biotyper IVD system Bruker Daltronics). The remaining portion of the environmental samples were preserved in 500 µL 1 × DNA/RNA Shield (Zymo Research Corporation).

DNA extraction and quantification

Clinical and environmental isolates were grown on Columbia horse blood plates (Oxoid Limited), then the DNA extracted using the DNeasy Blood & Tissue Miniprep Kit (Qiagen), following manufacturer’s instructions [22]. Metagenomic DNA from a subset of environmental samples, due to limited resources, was extracted directly from the environmental swab and water samples using the ZymoBIOMICS™ DNA Miniprep Kit (Zymo Research Corporation), following manufacturer’s instructions [23]. A sample of ZymoBIOMICS™ Microbial Community Standard (Zymo Research Corporation) was included, following manufacturer’s instructions [24]. DNA was assessed for concentration using the Qubit™ dsDNA BR Assay Kit (Thermo Fisher). Molecular weight and DNA integrity was confirmed using the Genomic DNA ScreenTape and reagents on the TapeStation 4150 (Agilent Technologies Inc.).

Library preparation

Clinical and environmental isolate DNA libraries were prepared using the Rapid Barcoding Kit 96 (SQK-RBK110.96) with a DNA input of between 50 and 200 ng, following manufacturer’s instructions [25]. For the environmental swab samples, metagenomic DNA libraries were prepared using the ONT Rapid PCR Barcoding Kit (SQK-RPB004) with a DNA input of 1–5 ng and following the manufacturers’ instructions [26]. ZymoBIOMICS™ Microbial Community DNA Standards (Zymo Research Corporation) were included, following manufacturer’s instructions [27].

Sequencing and basecalling

Up to 24 barcoded clinical and environmental isolate samples, or 12 barcoded environmental metagenomic samples were run together on a flow cell version R9.4.1 (Oxford Nanopore Technologies) using a MinION device for 72 h, using the default parameters on the MinKNOW software (v23.04.6). Basecalling was performed either by the MinKNOW software alongside sequencing or using the Guppy basecalling software (v6.5.7) [41], using the flip-flop high accuracy algorithm, with a minimum Q score of 8 and minimum depth of 40x.

Data analysis

Fastq files were quality checked (QC) using FastQC (v0.21.1) and MultiQC (v1.15) and those with a read depth of  > 40x  were included in further analysis [28, 29]. 25 clinical isolates, 60 environmental isolates, 36 metagenomics environmental samples were included in this analysis. Barcodes were trimmed from the reads using Guppy. All clinical and environmental isolates, and environmental metagenomic samples were analysed for the presence of AMR genes using KmerResistance 2.2 [30, 31] and for species using KmerFinder 3.2 (v3.0.2) [30, 32, 33]. Both programmes match query sequences to k-mer databases, for resistance genes and species respectively. Plasmids were identified using PlasmidFinder 2.1 (v2.0.1) [30, 34]. MLST types were identified using MLST (v2.0) [30, 35,36,37,38,39,40].

Clinical and environmental isolate fastq files were aligned to the reference genome for their species (C. freundii: GCF_003812345.1, C. portucalensis: GCA_023374935.1, C. youngae: GCF_015139575.1, E. cloacae: GCF_905331265.2, E. bugandensis: GCF_020042625.1, E. hormaechei: GCF_024218835.1, E. coli: GCF_000005845.2, K. pneumoniae: GCF_000240185.1, K. michigenensis: GCF_015139575.1) using MiniMap2 (v2.26) [41], then sorted and indexed using Samtools (v1.17) [42]. Alignments were visualised using Artemis (v18.1.0) [43]. Depth and coverage were calculated using Samtools. Consensus fasta files were created using Samtools and then dendograms for each species (with three or more isolates) were created using Parsnp (utilising maximal unique matches) (v1.7.4) [44,45,46,47] and visualised using iTOL (v6) [48].

Speciation for the results and discussion were as per the genomic speciation. Sequence data were deposited under BioProject PRJEB76684 on the European Nucleotide Archive and outlined in supplementary materials S4, S5 and S6.

Results

MALDI-TOF vs WGS for isolate speciation

All 25 clinical and 48 environmental isolates were subjected to both MALDI-TOF and WGS analysis. Non-concordant C. freundii by MALDI-TOF were speciated as C. portucalensis (one clinical isolate and one environmental isolate) or C. youngae (three environmental isolates) in WGS. One MALDI-TOF call of C. braakii/freundii was called as C. youngae using WGS. One environmental isolate that MALDI-TOF identified as C. freundii was speciated as P. mirabilis, possibly a mixed culture in which the P. mirabilis may have been present on the original selective plate, but then grew on the non-selective plate used for DNA extraction and sequencing. Non-matching E. cloacae were speciated as E. asburiae (one environmental isolate) and E. hormaechei (nine, including all five of the clinical isolates and four environmental isolates). All MALDI-TOF and WGS speciation for E. coli and K. pneumoniae was concordant. All the MALDI-TOF-called K. oxytoca were speciated as either K. michigenensis (4/5 isolates) or K. grimontii (1/5 isolates) using WGS. See Table 1 and supplementary materials S5.

Table 1 Concordance of MALDI-TOF speciation vs WGS for both clinical and environmental isolates. All non-concordant samples were either of the C. freundii or E. cloacae complex

Environmental sampling

210 environmental samples were collected, including 43 stick swabs, 115 sponge swabs and 52 wastewater samples. In total, 195 bacteria were isolated from 36 (16.9%) environmental samples, 76 isolates were Enterobacterales: 6 from stick swabs, 5 from sponge swabs and 65 from non-potable water samples. There was a significant difference in the proportion of samples from which Enterobacterales species were isolated from, when the one-way ANOVA was applied (p =  ≤ 0.0001). When swab types were compared using Tukey’s multiple comparisons test, there was no significant difference between stick and sponge swabs (p = 0.4286), but a significant difference between stick swabs and wastewater (p = ≤ 0.0001) and sponge swabs and wastewater (p = ≤ 0.0001). Enterobacterales were isolated from 51.9% of water samples, compared with 11.6% of stick swabs and 4.3% of sponge swabs (see Table 2).

Table 2 Environmental sample types collected, and the number of each sample type from which Enterobacterales were isolated

Enterobacterales were isolated from 36 different environmental samples (mean = 2.2 isolates per site, SD = 1.4). See Table 3. Table 4 details the sample sites from which each species was isolated, Table 5 describes the ward areas in which they were isolated. For full data see supplementary materials S6.

Table 3 Number of each species isolated from each different type of environmental sample
Table 4 Environmental isolates present divided by sample site and type depicting an environmental reservoir density by location
Table 5 Environmental isolates present divided by room indicating the environmental reservoir densities by room-type

Table 6 details the environmental sample swab types, total number of bacterial species, and the total number of reads obtained from each of the 36 metagenomic sequencing sample sites that Enterobacterales were also present. There was no significant difference in the number of sequencing reads obtained, when comparing the swab type, or when comparing the swab location type. When comparing rooms, there were no significant differences, apart from Bay (beds 1–4) and Bay bed 20 (p = 0.0321). All tests were undertaken using One Way ANOVA and Tukey’s multiple comparison test. Speciation was based on genomic data, so only those that passed QC are included in the analysis.

Table 6 Environmental swab metagenomic samples, from which Enterobacterales were isolated, showing which type of sample was taken (sponge swab, stick swab or wastewater sample), number of bacterial species identified (including Gram negatives and Gram positives), and number of reads obtained from metagenomic sequencing

AMR genes

Clinical isolates

Of the clinical isolates, blaOXA-48 was identified in 2 of the total 28 (7.1%) isolates (one C. portucalensis and one K. pneumoniae). blaNDM-1 was identified in nine (32.1%) clinical isolates (one C. portucalensis, five E. hormaechei, one E. coli and two K. pneumoniae). blaOXA genes on the CARBA-5 panel were found in 84% of clinical isolates and blaNDM genes on the CARBA-5 panel were found in 76% of clinical isolates, see Table 7. See supplementary materials S7 for full list of genes found.

Table 7 Clinical isolates with blaOXA-48, at least one blaOXA (on the CARBA-5 panel), total number of isolates with any blaOXA gene, blaNDM-1, at least one blaNDM (on the CARBA-5 panel) and total number of isolates with any blaNDM gene

Environmental isolates

Of the environmental isolates, blaOXA-48 was identified in 11 of the total 60 (18%) isolates (See Table 8). blaNDM-1 was identified in one (1.7%) environmental isolate, a K. michigenensis. See supplementary materials S8 for full list of genes found.

Table 8 Environmental isolates with blaOXA-48, at least one blaOXA (on the CARBA-5 panel), total number of isolates with any blaOXA gene, blaNDM-1, at least one blaNDM (on the CARBA-5 panel) and total number of isolates with any blaNDM gene

Metagenomic environmental samples

blaOXA-48 was identified in 1/36 (2.8%) metagenomic environmental sample, and three other blaOXA genes on the CARBA-5 panel (blaOXA-204, blaOXA-370 and blaOXA-515) were found together in 1/36 (2.8%) other metagenomic environmental sample. No blaNDM genes were identified in any of the 36 metagenomic environmental samples.

Plasmids

Clinical isolates

The clinical isolates had a mean number of 5.6 plasmids (SD = 2.3) and a total of 20 different plasmids were identified. Table 9 describes the plasmid types and species they were found in. Col440II was the only plasmid found across all clinical isolate species. IncFIB(AP001918) was the most commonly identified plasmid.

Table 9 Plasmid types found in percentage of clinical isolates

Environmental isolates

The environmental isolates had a mean number of 4.3 plasmids (SD = 2.3) and a total of 42 different plasmids were identified. Supplementary materials S9 describes the plasmids and species they were found in. Col(IRGK) was the most commonly identified plasmid in 35/60 (59.3%) isolates, no plasmid was ubiquitous across all environmental isolate species.

Metagenomic environmental samples

Plasmids for Enterobacterales were found in 24/36 (66.6%) environmental metagenomic samples (see Supplementary materials S10). The mean number of plasmids was 3.1 (SD = 3.6).

Metagenomic species

When species were analysed in the metagenomic samples, reads mapping to the Enterobacterales order were found in 64%-75% of sites (see Table 10). When analysed, there was no significant difference in the number of reads that mapped to Citrobacter, Enterobacter, Escherichia or Klebsiella genera, using One Way ANOVA and Tukey’s multiple comparison test.

Table 10 Number of metagenomic samples sites for which reads mapped to Enterobacterales genus sequences

Transmission and patient metadata

One large cluster of C. freundii isolates that mapped to strain N16-03880 were all from the same room Bay (beds 29–32) but came from both the bay HWB DWT as well as the bathroom HWB DWT samples, suggesting cross-contamination of the two sinks (see Fig. 1A, C. freundii cluster #1). These, along with a toilet bowl swab from Bay (beds 15–18) were the only C. freundii found to harbour blaOXA-48 genes (see Fig. 1A). All the C. youngae environmental isolates mapped to strain CF10 (see Fig. 1B, C. youngae cluster #1). The C. portucalensis isolates all mapped to different strains. The clinical isolate (patient Z1) mapped mainly (85% query coverage) to strain PNUCL1, and the environmental isolates mapped to FDAARGOS_617 and SWHIN_111. All harboured blaOXA genes, and the patient isolate additionally harboured blaNDM genes (see Fig. 1C).

Fig. 1
figure 1figure 1

Dendograms for species sequenced in this study

The E. asburiae environmental isolates all mapped to different strains: A2563, 2497 and RHBSTW-00542 (see Fig. 1D). The E. cloacae environmental isolates all mapped to strain GGT036. Neither species were found to harbour any blaOXA or blaNDM genes (see Fig. 1E). All of the E. hormaechei clinical isolates mapped to strain SH19PTE2 and were found to harbour blaNDM genes, whereas the environmental isolates each mapped to a different strain: Y323, F2, AR_0365 and RHBSTW-00086, none of which harboured either blaOXA or blaNDM genes (see Fig. 1F).

Except patient Z2, all patient E. coli isolates mapped most closely to strain 035152. The 035152 isolates had a similar AMR profile, with all (9/9) harbouring blaOXA-181 and 7/9 harbouring blaOXA-484. Patient Z2, which mapped most closely to E. coli VREC0864, clustered with EI106, an environmental isolate from a water sample from the pantry sink’s drain waste trap, which also mapped most closely to VREC0864 (see Fig. 1G, E. coli cluster #1).

The environmental isolates that mapped to K. michigenensis mapped most closely to strains RHB20-CO2, K518 and RHBSTW-00409 (two isolates from the same HWB DWT in Bay (beds 15–18). Only one was found to harbour blaOXA and blaNDM genes (see Fig. 1H). Except Z5 (strain RHBSTW-00113), all of the K. pneumoniae clinical isolates mapped most closely to strain B5617 (see Fig. 1I, K. pneumoniae cluster #2). Patient isolates Z3, Z4, Z5, Z6, Z17, Z18 and Z19 (all strain B5617) showed a similar AMR profile, all harbouring blaOXA-181, blaOXA-232 and blaOXA-484, and blaNDM-5. There was a cluster of environmental isolates that mapped most closely to K. pneumoniae strain 3987, all isolated from hand wash basin drain waste traps, but in differing patient bed areas around the ward (see Fig. 1I, K. pneumoniae cluster #1).

9/20 (45%) patients were housed on ward 7D during their time as an inpatient and 9/9 (100%) patients had stayed in Bay (beds 15–18), Bay (beds 29–32) or in one of the single-occupancy bed rooms, each of which were noted to have had higher numbers of Enterobacterales isolated from the environmental samples. 4/9 (44%) of these patients had been housed in more than one of the bays, one patient had been on both Bay (beds 15–18), bay (beds 29–32) and Bay 19.

Fig. 1. Phylogram of CRE species isolated in this study and metadata including location (of either sample collection for environmental isolates, or beds inhabited for clinical isolates), strain and presence or absence of blaOXA and blaNDM genes. A) C. freundii, B) C. portucalensis, C) C. youngae, D) E. asburiae, E) E. cloacae, F) E. hormaechei, G) E. coli, H) K. michigenensis and I) K. pneumoniae. Branch lengths indicate number of substitutions divided by the length of the genome sequence. HWB = hand wash basin, DWT = drain waste trap, WoW = workstation on wheels. blaOXA = blaOXA gene on CARBA-5 panel, blaNDM = blaNDM gene on CARBA-5 panel, EI = Environmental isolate. The top branch of each phylogram is the reference genome (e.g. HS11286 for K. pneumoniae), against which the clinical and environmental isolates were mapped.

Discussion

There was evidence for the potential for patient-patient transmission for E. hormaechei, E. coli and K. pneumoniae. This conclusion is supported by the patient spatio-temporal data collected (supplementary materials S1). This study also identified environment-patient transmission. The E. coli isolate from patient Z2 clustered with EI106, an environmental isolate from a wastewater sample taken from the communal pantry sink’s DWT, although they did not share the same plasmids. Whilst the clinical isolate was found to harbour the blaNDM-1 gene, no blaOXA or blaNDM genes were found in the environmental isolate. This may have been a loss or gain of function; bacteria found within the environment are less likely to encounter antibiotics or their residues compared with those in a hospitalised patient, so that isolate may have lost the plasmid containing these resistance genes. Conversely, they may gain plasmids from commensal organisms within the host [49].

It was noted that Enterobacterales were commonly isolated from environmental samples in Bay (beds 15–18), Bay (beds 29–32) and the single-occupancy bays. All the patients in this study located on ward 7D had been housed in at least one of these bays and 44% had been housed in more than one of these bays, suggesting potential host-related reservoirs. One report suggests that patients are on average 73% more likely to acquire a HAI if the patient previously occupying their room was colonised or infected, [50]. This suggests enhanced location-specific IPC would be beneficial when colonised or infected patients have been identified, and that isolation of these patients may not be enough. Enterobacterales were also isolated from environmental samples in non-clinical areas such as the shared-use pantry and the staff office. Studies have shown that HCWs can be colonised when handling patients and infected materials, it’s also possible that colonised patients or HCWs may have caused reservoirs in the pantry due to transmission of organisms via the faecal-oral route [51].

It is well known that it can be difficult to standardise environmental sampling, and that the recovery rate can vary between sampling tools. In this study, multiple environmental surfaces were swabbed with both cotton and sponge swabs and then isolates of interest grown on selective agar to minimise selection bias and maximise the chance of isolating a CRE. In this study, there was a much greater recovery from wastewater. For recovery from dry surfaces, sponge swabs appeared to isolate more Enterobacteriales than cotton swabs, perhaps due to the larger surface area for absorption (surface area of swab site remained the same).

Even with comprehensive patient metadata, such as bed movements, it can be difficult to confidently infer transmission of a clonal isolate. Environmental isolates are especially complex, as the bacteria may have been present for long periods, for example by forming hard to remove biofilms in U-bends, which is in itself difficult to monitor and identify provenance [52]. For human pathogens, the ward environment may not be optimal for growth, so their doubling time may be slower. Genetic-relatedness cutoffs to determine phylogeny tend to be calculated depending on the sample number, type and environments, as well as the species, as some have faster molecular clocks than others, and so the number of SNPs difference does not necessarily reflect the closeness of two isolates. These data are often missing for isolates extracted from the environment [53].

Whilst all clinical isolates tested positive for OXA-48-like or NDM-like variants using the CARBA-5 test, only 84% of the clinical isolates were found to harbour blaOXA genes that the CARBA-5 panel tested for when whole genome sequenced. Studies have shown the specificity of the CARBA-5 test to vary from 96% (from blood cultures) to 100% from isolates and rectal swabs [54,55,56,57]. The clinical sites that the swabs were taken from in this study varied and the majority came from CRE screens, usually rectal, rather than sites of infection. It is possible that the CARBA-5 panel picked up blaOXA genes from other colonising species present in CRE screening samples, as these were not sterile sites. Most of the clinical and environmental isolates found to harbour the blaOXA-48 gene also had other blaOXA genes present, as well as other resistance genes, such as blaCTX and blaSHV. These are often found together on the same plasmids and are likely to be transferred between bacteria collectively [4]. Indeed, all of the C. portucalensis isolates identified in this study were found to harbour blaOXA genes, with the clinical isolate also containing blaNDM-1.

Whilst Col440II, which does not carry resistance or virulence genes, was the only plasmid found across all clinical isolate species, IncFIB(AP001918), the most commonly identified plasmid, is linked with resistance to several antimicrobial classes, including β-lactams, aminoglycosides, sulfonamides and tetracyclines, but was less commonly identified in environmental isolates [58, 59]. The prevalence of IncFIB(AP001918) in clinical samples suggests that genomic analysis of plasmids as well as isolates is important for enhanced IPC surveillance.

In this study, the initial MALDI-TOF speciation for C. freundii, E. cloacae and K. oxytoca isolates did not all match with the WGS speciation. Both methods were more consistent when comparing environmental rather than clinical isolates. All three species complexes contain multiple closely related species, which can make it difficult to fully resolve using conventional methods. C. portucalensis is a relatively newly described clinical pathogen but has the capacity to harbour and transmit AMR genes, thus identifying it to species level may be important in the future [60].

The use of genomic analysis in enhanced outbreak surveillance technologies provides greater detail on the potential transmission of MDROs and their associations with patients, HCWs and the ward environment. Sequencing platforms such as ONT shows promise, especially for use in low- and middle-income countries, as long read sequencing enables read lengths of thousands, rather than hundreds, of base pairs which is especially useful when resolving speciation in metagenomic samples.

Conclusions

Understanding the resistance genes, plasmids and sequencing types present in an environment can provide greater resolution than phenotypic and other molecular methods, helping to identify targeted IPC interventions in outbreak situations. As a result of the evidence from this study, highlighting the presence of CREs in wastewater, the hospital estates team has since replaced all of the sink U-bends on the sampled ward, as well as reviewed and revised ward IPC practices to reduce potential transmission risks. Due to the number of colonised patients found, this study also recommends the use of more widespread CRE screening for hospitalised patients, to enable interventions to reduce the risks from human and environmental reservoirs and therefore reduce risks to vulnerable patients.

Overview:

  • Putative patient-patient and patient-ward transmission was identified, utilising WGS and patient metadata

  • CREs were more commonly isolated from wastewater samples than either stick swabs or sponge swabs

  • All patients in this study tested positive for OXA-48-like carbapenemase variants using the NG-Test® CARBA-5 rapid diagnostic test, 84% of patient isolates harboured a blaOXA gene present on the CARBA-5 panel

  • blaOXA and blaNDM genes were identified in fewer environmental CRE isolates, compared with clinical isolates, suggesting either different populations, or a loss/gain of plasmids or genes

  • ONT sequencing can expedite clinical decisions whilst awaiting reference laboratory results, providing economic and patient care benefits.

Availability of data and materials

All datasets have been made freely available and can be found within the supplementary materials. Further patient metadata are available from the corresponding author upon reasonable request. Sequence data were deposited under BioProject PRJEB76684 on the European Nucleotide Archive.

Abbreviations

AMR:

Antimicrobial resistance

AST:

Antibiotic sensitivity test

CRE:

Carbapenem-resistant Enterobacterales

DWT:

Drain waste trap

HCAI:

Healthcare associated infection

HWB:

Hand wash basin

IPC:

Infection prevention and control

MDRO:

Multi-drug resistant organism

NDM:

New Delhi metallo-beta lactamase

ONT:

Oxford Nanopore Technologies

RDT:

Rapid diagnostic test

SD:

Standard deviation

WoW:

Workstation on wheels

References

  1. World Health Organization. Guidelines for the prevention and control of carbapenem-resistant enterobacteriaceae, Acinetobacter baumannii and Pseudomonas aeruginosa in health care facilities. 2017. http://apps.who.int/bookorders.

  2. Mills MC, Lee J. The threat of carbapenem-resistant bacteria in the environment: Evidence of widespread contamination of reservoirs at a global scale. Environ Pollut. 2019;255:113143.

    Article  CAS  PubMed  Google Scholar 

  3. Tängdén T, Giske CG. Global dissemination of extensively drug-resistant carbapenemase-producing enterobacteriaceae: clinical perspectives on detection, treatment and infection control. J Int Med. 2015;277:501–12.

    Article  Google Scholar 

  4. Kopotsa K, OseiSekyere J, Mbelle NM. Plasmid evolution in carbapenemase-producing Enterobacteriaceae: a review. Annals New York Acad Sci. 2019;1457:61–91.

    Article  CAS  Google Scholar 

  5. Nordmann P, Gniadkowski M, Giske CG, Poirel L, Woodford N, Miriagou V, et al. Identification and screening of carbapenemase-producing Enterobacteriaceae. Clin Microbiol Infect. 2012;18:432–8.

    Article  CAS  PubMed  Google Scholar 

  6. Gupta N, Limbago BM, Patel JB, Kallen AJ. Carbapenem-resistant enterobacteriaceae: epidemiology and prevention. Clin Infect Dis. 2011;53:60–7.

    Article  PubMed  Google Scholar 

  7. Apanga PA, Ahmed J, Tanner W, Starcevich K, VanDerslice JA, Rehman U, et al. Carbapenem-resistant Enterobacteriaceae in sink drains of 40 healthcare facilities in Sindh, Pakistan: a cross-sectional study. PLoS ONE. 2022;17(2):e0263297.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Islam CA, Islam MA, Hasan M, Hossain RI, Nabi MI, Rahman A, et al. Environmental spread of New Delhi metallo-lactamase-1-producing multidrug-resistant bacteria in Dhaka, Bangladesh environmental spread of New Delhi metallo-β-lactamase-1-producing multidrug-resistant bacteria in Dhaka. Appl Environ Microbiol. 2017;83:00793–817.

    Article  Google Scholar 

  9. Zurfluh K, Bagutti C, Brodmann P, Alt M, Schulze J, Fanning S, et al. Wastewater is a reservoir for clinically relevant carbapenemase- and 16s rRNA methylase-producing Enterobacteriaceae. Int J Antimicrob Agents. 2017;50(3):436–40.

    Article  CAS  PubMed  Google Scholar 

  10. Zenati K, Touati A, Bakour S, Sahli F, Rolain JM. Characterization of NDM-1- and OXA-23-producing Acinetobacter baumannii isolates from inanimate surfaces in a hospital environment in Algeria. J Hosp Infect. 2016;92(1):19–26.

    Article  CAS  PubMed  Google Scholar 

  11. Perry M, van Bunnik B, McNally L, Wee B, Munk P, Warr A, et al. Antimicrobial resistance in hospital wastewater in Scotland: a cross-sectional metagenomics study. The Lancet. 2019;394:S1.

    Article  Google Scholar 

  12. Department of Health and Social Care. Sewage in water: a growing public health problem. 2022. https://www.gov.uk/government/news/sewage-in-water-a-growing-public-health-problem. (Accessed 23 Jun 2024)

  13. Klümper U, Gionchetta G, Catão E, Bellanger X, Dielacher I, Elena AX, et al. Environmental microbiome diversity and stability is a barrier to antimicrobial resistance gene accumulation. Commun Biol. 2024;7(1):706.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Struelens MJ, Ludden C, Werner G, Sintchenko V, Jokelainen P, Ip M. Real-time genomic surveillance for enhanced control of infectious diseases and antimicrobial resistance. Front Sci. 2024;25:2.

    Google Scholar 

  15. Clinical and laboratory standards institute. Performance standards for antimicrobial susceptibility testing. In: CLSI supplement M100. 31st ed. 2021.

  16. Chen Z, Azman AS, Chen X, Zou J, Tian Y, Sun R, et al. Global landscape of SARS-CoV-2 genomic surveillance and data sharing. Nat Genet. 2022;54(4):499–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Boolchandani M, D’Souza AW, Dantas G. Sequencing-based methods and resources to study antimicrobial resistance. Nat Rev Genetics. 2019;20:356–70.

    CAS  PubMed  Google Scholar 

  18. European committee on antimicrobial susceptibility testing. Breakpoint tables for interpretation of MICs and zone diameters (v12). 2022.

  19. Hardy Diagnostics. NG-Test CARBA 5. 2021.

  20. Ali S, Wilson APR. Effect of poly-hexamethylene biguanide hydrochloride (PHMB) treated non-sterile medical gloves upon the transmission of Streptococcus pyogenes, carbapenem-resistant E coli, MRSA and Klebsiella pneumoniae from contact surfaces. BMC Infect Dis. 2017;17(1):1–8.

    Article  Google Scholar 

  21. Yetiş Ö, Ali S, Karia K, Bassett P, Wilson P. Enhanced monitoring of healthcare shower water in augmented and non-augmented care wards showing persistence of Pseudomonas aeruginosa despite remediation work. J Med Microbiol. 2023;72(5):001698.

    Article  Google Scholar 

  22. Qiagen. sample to insight DNeasy ® blood & tissue handbook. 2023.

  23. Zymo research. ZymoBIOMICSTM DNA Miniprep Kit. https://files.zymoresearch.com/protocols/_d4300t_d4300_d4304_zymobiomics_dna_miniprep_kit.pdf. Accessed 22 Aug 2023.

  24. Zymo research. ZymoBIOMICS® microbial community standard. https://files.zymoresearch.com/protocols/_d6300_zymobiomics_microbial_community_standard.pdf. Accessed 22 Aug 2023.

  25. Oxford nanopore technologies. Rapid sequencing gDNA - barcoding (SQK-RBK110.96). 2022. https://community.nanoporetech.com/docs/prepare/library_prep_protocols/rapid-barcoding-kit-96-sqk-rbk110-96/v/rbk_9126_v110_revm_24mar2021

  26. Oxford nanopore technologies. Rapid PCR barcoding kit (SQK-RPB004) protocol. 2020. p. 1–4. https://community.nanoporetech.com/protocols/rapid-pcr-barcoding/checklist_example.pdf. Accessed 28 Jun 2021.

  27. Zymo research. ZymoBIOMICSTM microbial community DNA standard. https://files.zymoresearch.com/protocols/_d6305_d6306_zymobiomics_microbial_community_dna_standard.pdf. Accessed 22 Aug 2023.

  28. Babraham bioinformatics. FastQC. 2019. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 28 Jun 2021

  29. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Clausen PTLC, Aarestrup FM, Lund O. Rapid and precise alignment of raw reads against redundant databases with KMA. BMC Bioinformatics. 2018;19(1):1–8.

    Article  Google Scholar 

  31. Clausen PTLC, Zankari E, Aarestrup FM, Lund O. Benchmarking of methods for identification of antimicrobial resistance genes in bacterial whole genome data. J Antimicrob Chemother. 2016;71(9):2484–8.

    Article  CAS  PubMed  Google Scholar 

  32. Hasman H, Saputra D, Sicheritz-Ponten T, Lund O, Svendsen CA, Frimodt-Moller N, et al. Rapid whole-genome sequencing for detection and characterization of microorganisms directly from clinical samples. J Clin Microbiol. 2014;52(1):139–46.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Larsen MV, Cosentino S, Lukjancenko O, Saputra D, Rasmussen S, Hasman H, et al. Benchmarking of methods for genomic taxonomy. J Clin Microbiol. 2014;52(5):1529–39.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Carattoli A, Zankari E, Garciá-Fernández A, Larsen MV, Lund O, Villa L, et al. In Silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, et al. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50(4):1355–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Bartual SG, Seifert H, Hippler C, Luzon MAD, Wisplinghoff H, Rodríguez-Valera F. Development of a multilocus sequence typing scheme for characterization of clinical isolates of Acinetobacter baumannii. J Clin Microbiol. 2005;43(9):4382–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Griffiths D, Fawley W, Kachrimanidou M, Bowden R, Crook DW, Fung R, et al. Multilocus sequence typing of Clostridium difficile. J Clin Microbiol. 2010;48(3):770–8.

    Article  CAS  PubMed  Google Scholar 

  38. Lemee L, Dhalluin A, Pestel-Caron M, Lemeland JF, Pons JL. Multilocus sequence typing analysis of human and animal Clostridium difficile isolates of various toxigenic types. J Clin Microbiol. 2004;42(6):2609–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, et al. Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006;60(5):1136–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Jaureguy F, Landraud L, Passet V, Diancourt L, Frapy E, Guigon G, et al. Phylogenetic and genomic diversity of human bacteremic Escherichia coli strains. BMC Genomics. 2008;9:1–14.

    Article  Google Scholar 

  41. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18):3094–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Carver T, Harris SR, Berriman M, Parkhill J, McQuillan JA. Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics. 2012;28(4):464–9.

    Article  CAS  PubMed  Google Scholar 

  44. Treangen TJ, Ondov BD, Koren S, Phillippy AM. The harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 2014;15(11):1–5.

    Article  Google Scholar 

  45. Price MN, Dehal PS, Arkin AP. FastTree 2 - Approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5(3):e9490.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Edgar RC. Muscle: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006;172(4):2665–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Letunic I, Bork P. Interactive tree of life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Dunai A, Spohn R, Farkas Z, Lazar V, Györkei Á, Apjok G, et al. Rapid decline of bacterial drug-resistance in an antibiotic-free environment through phenotypic reversion. Elife. 2019. https://doiorg.publicaciones.saludcastillayleon.es/10.7554/eLife.47088.001.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Rawlinson S, Ciric L, Cloutman-Green E. How to carry out microbiological sampling of healthcare environment surfaces? A review of current evidence. J Hospital. 2019;103:363–74.

    CAS  Google Scholar 

  51. Wohrley JD, Bartlett AH. The role of the environment and colonization in healthcare-associated infections. In: healthcare-associated infections in children. Springer International Publishing; 2019. p. 17–36.

  52. Yin W, Xu S, Wang Y, Zhang Y, Chou SH, Galperin MY, et al. Ways to control harmful biofilms: prevention, inhibition, and eradication. Crit Rev Microbio. 2021;47:57–78.

    Article  CAS  Google Scholar 

  53. Mustapha MM, Srinivasa VR, Griffith MP, Cho ST, Evans DR, Waggle K, et al. Genomic diversity of hospital-acquired infections revealed through prospective whole-genome sequencing-based surveillance. mSystems. 2022;7(3):e01384-e1421.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Yoon J, Kim CH, Yoon SY, Lim CS, Lee CK. Application of a multiplex immunochromatographic assay for rapid identification of carbapenemases in a clinical microbiology laboratory: performance and turn-around-time evaluation of NG-test Carba 5. BMC Microbiol. 2021;21(1):1–7.

    Article  CAS  Google Scholar 

  55. Saito K, Mizuno S, Nakano R, Tanouchi A, Mizuno T, Nakano A, et al. Evaluation of NG-Test CARBA 5 for the detection of carbapenemase-producing Gram-negative bacilli. J Med Microbiol. 2022;71(6):001557.

    Article  CAS  Google Scholar 

  56. Takissian J, Bonnin RA, Naas T, Dortet L. NG-test carba 5 for rapid detection of carbapenemase-producing Enterobacterales from positive blood cultures. Antimicrob Agents Chemother. 2019;63(5):10–128.

    Article  Google Scholar 

  57. Vasilakopoulou A, Karakosta P, Vourli S, Kalogeropoulou E, Pournaras S. Detection of kpc, ndm and vim-producing organisms directly from rectal swabs by a multiplex lateral flow immunoassay. Microorganisms. 2021;9(5):942.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Shankar C, Muthuirulandi Sethuvel DP, Neeravi AR, Venkatesan M, Devanga Ragupathi NK, Anandan S, et al. Identification of plasmids by PCR based replicon typing in bacteremic Klebsiella pneumoniae. Microb Pathog. 2020;1:148.

    Google Scholar 

  59. Zhou W, Lin R, Zhou Z, Ma J, Lin H, Zheng X, et al. Antimicrobial resistance and genomic characterization of Escherichia coli from pigs and chickens in Zhejiang. China Front Microbiol. 2022;24:13.

    Google Scholar 

  60. Cao X, Xie H, Huang D, Zhou W, Liu Y, Shen H, et al. Detection of a clinical carbapenem-resistant Citrobacter portucalensis strain and the dissemination of C. portucalensis in clinical settings. J Glob Antimicrob Resist. 2021;27:79–81.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

This research was funded by the Royal Free London NHS Foundation Trust.

Author information

Authors and Affiliations

Authors

Contributions

EW, LE, TM and DM designed the study. EW oversaw the clinical aspects of this project and coordinated the collection and processing of the clinical isolates. SQ, JH, VP, LC, CO were instrumental in the ward environmental analysis and collection of the environmental samples. AW, SA, OY, EC, IW, MM, SY, KK, ES undertook phenotypic analysis on the environmental isolates. LE undertook clinical and environmental isolate and metagenomic DNA extraction, WGS and metagenomic sequencing, with advice from SR. LE was responsible for data analysis and the writing of this manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Linzy Elton.

Ethics declarations

Ethical approval and consent to participate

This diagnostic service evaluation was undertaken at the request of the Division of Nursing in support of the Infection Prevention and Control team as part of enhanced standard of care practices for IPC outbreak investigation, in accordance with the Declaration of Helsinki. This project is under HRA approval for ELCID and registered with Health Services London. IRAS project ID: 283,831, National HRA REC reference: 20/HRA/4928, Local REC: RFL R&D ref: 134,895, HSL Project number: 1529.

Consent for publications

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Elton, L., Williams, A., Ali, S. et al. Tracing the transmission of carbapenem-resistant Enterobacterales at the patient: ward environmental nexus. Ann Clin Microbiol Antimicrob 23, 108 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12941-024-00762-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12941-024-00762-8

Keywords