Introduction

Increased antimicrobial resistance (AMR) constitutes a major public health concern in healthcare and community settings [1]. At the healthcare facility level, multidrug-resistant organisms (MDRO) can spread by physical contact between patients and healthcare workers or through contaminated surfaces and medical devices. Furthermore, the transfer of patients between hospitals and other healthcare institutions, such as nursing homes (NHs), within and between healthcare networks plays an important role in the spread of MDROs [2]. As data on MDRO carriage in Belgian healthcare institutions are scarce and transmission between institutions has not been studied on a large scale, our study explored ESBL-E, CPE and VRE prevalence and transmission in three tertiary care hospitals and 13 NHs in Flanders, Belgium, using a combination of conventional culture techniques and WGS.

Materials and methods

Study design

A cross-sectional MDRO prevalence survey, the i-4-1-Health project, was organized in the cross-border region of the southern part of the Netherlands and Flanders, Belgium. In this paper, the results obtained between October 2017 and February 2019 from the Flemish hospitals (n = 3) and NHs (n = 13) are presented.

Sampling and microbiological analysis

The microbiological method of the project is described by Kluytmans – van den Bergh et al. [3].

Whole-genome sequencing

Sequencing was performed on an Illumina MiSeq sequencer, using the MiSeq reagent kit v2 generating 250-bp paired-end reads. wgMLST was performed using the BioNumerics software v7.6.3. (bioMérieux, France). Clonal relatedness between isolates was determined based on the similarity of wgMLST allelic profiles. By combining both sequencing and epidemiological data of well-described outbreaks, similarity thresholds for clonal relatedness were determined. Sequencing data should be combined with the metadata to make conclusions for an outbreak, especially in the ‘grey zone’(Table S1). Using a BLAST-based approach requiring at least 95% identity with the reference sequence and at least 95% reference length coverage, genes known to confer resistance were identified from the assembled genomes. Reference sequences and mutations from the Center for Genomic Epidemiology’s ResFinder (database version 2019-08-21) and PointFinder (database version 2019-07-02) databases were used, respectively.

Results

Sample and patient information

In total, 2536 perianal or stoma swabs were obtained from patients in Flemish tertiary care hospitals (n = 1878) and NH residents (n = 658). After the exclusion of six swabs (0.2%) due to poor sampling quality and 125 swabs (6.7%) due to duplicate sampling from previously sampled hospitalized patients, 2405 swabs remained for analysis. Of those, 1748 (73%) were obtained in hospitals and 657 (27%) in NHs (Table 1).

Table 1 MDRO prevalence in perianal swabs from hospitalized patients (H) and nursing home residents (NH)

Prevalence of ESBL-E, CPE and VRE in hospitals and NHs

The percentage of screened patients or residents is based on the number of patients of residents present on a ward. The genotypic confirmed ESBL-E, CPE and VRE isolates are shown in Table 1. Of the MRDOs, ESBL-E were cultured most often (11.4% in hospitals and 14.0% in NHs) compared to a lower prevalence of VRE (0.9% in hospitals and 0.3% in NHs) and CPE (1.0% in hospitals and undetected in NHs). Significant differences in prevalence of MDRO between the individual hospitals (p < 0.001) and between NHs (p < 0.001) were detected, while no significant variation in prevalence between the hospital and NH sector was observed.

Phenotype-genotype correlation for MDRO detection

In all VRE (15 E. faecium and 3 E. avium strains with an MIC for vancomycin of > 4 mg/l), a resistance gene could be detected (14 vanA genes and 4 vanB genes). However, for ESBL-E and CPE, phenotypic resistance could not always be confirmed by the presence of a resistance gene or mutation, resulting in lower genetic confirmation rates, especially for CPE (Table 2). The presence of other β-lactamase genes and porins are usually responsible for an ESBL-like phenotype.

Table 2 Confirmation rates of phenotypical ESBL-E and CPE strains by the presence of resistance genes

ESBL-E were detected in 291 of the 2405 swabs, corresponding with 328 ESBL-E strains (predominantly E. coli and K. pneumoniae) as 37 swabs contained multiple ESBL-E strains. Over half (51.5%) of the ESBL production capability could be attributed to the presence of the blaCTX-M-15 gene.

To study local and regional clonal relatedness of ESBL-E. coli, separate spanning trees were created for each hospital and its surrounding NHs (Fig. 1A–C). NH 3 is located in a range of < 10 km; NHs 4, 5 and 11 in a range of 11–25 km; NHs 2, 6, 8, 10, 12 and 13 in a range of 26–50 km and NHs 1, 7 and 9 in a range of 51–75 km from the tertiary care institutions.

Fig. 1
figure 1

Three minimum spanning trees based on wgMLST analysis of 233 ESBL-E. coli isolates per geographical region around the three hospitals (hospital 1 (A), hospital 2 (B) and hospital 3 (C) and their surrounding NHs). Isolates are represented by circles connected by branches proportional to the allelic distance, colours of the circle represent the origin location, the shading represent clonal clusters and ST number represent sequence types

The first region enclosed 84 ESBL-E. coli isolates from hospital 1 and four NHs (Fig. 1A). Five clonal clusters could be visualized, of which two are large ones with over 10 ST131 ESBL-E. coli isolates. The three smaller clusters comprise two isolates each and are within hospital 1 (two clusters of ST2797 and ST38) and NH 7 (a cluster of ST405).

In the second region, the isolates obtained from the hospitalized patients predominate in five of the seven clusters (Fig. 1B). Four clonal clusters (clusters of ST131, ST69, ST23 and ST10) contain isolates obtained from hospital 2 or NH13. The large cluster of ST131 comprises isolates of hospital 2, NH4, NH6 and NH13.

In the third region around hospital 3, a high level of clonality and clustering among the collected ESBL-E.coli could be observed (Fig. 1C). Five of nine clonal clusters consist of isolates obtained within hospital 3 (cluster of ST393 and ST131), NH8 (cluster of ST131) and NH12 (cluster of ST127 and ST131). NH3 (seven isolates in total in the two clusters) and hospital 3 (ten isolates in total in the two clusters) have many related isolates and are located in a range of < 10 km. In all regions, there were clusters of isolates in the hospital and the attached nursing homes indicative of transmission between healthcare institutes.

A minimum spanning tree of all ESBL-K. pneumoniae isolates (n = 66) is visualized in Fig. 2. One cluster contains two isolates (ST405) from NH 8 (one) and NH 10 (one) which differ by 3 alleles and have a similarity of 99.9%. NH 8 and 10 are from the same region (within 11–25 km).

Fig. 2
figure 2

A minimum spanning tree based on wgMLST analysis of 66 ESBL-producing Klebsiella pneumoniae isolates. Isolates are represented by circles connected by branches proportional to the allelic distance, colours of the circle represent the location, the shading represent clonal clusters, branches and numbers represent allelic differences between isolates and ST number represent sequence types

Carbapenemase-producing Enterobacterales

Twenty-three (1.3%) CPE isolates were detected; all of them originated from hospitalized patients. These CPE isolates consisted of six different species of which K. pneumoniae (n = 9) and Citrobacter spp. (n = 7) were predominant. The presence of the blaOXA-48 gene explained the carbapenem resistance in all but one isolate, a Citrobacter freundii strain carrying a blaVIM-4 resistance gene. All the blaOXA-48 positive isolates contained the ca. 64kb IncL/M plasmid. Analysis of clonal relatedness of these CPE revealed a cluster of 4 ST15 K. pneumoniae CPE in ward 3 of hospital 1 (Fig. 3). In total, 10 blaOXA-48 producing bacteria of three different species, K. pneumoniae, Citrobacter spp. and E. cloacae complex were detected in this specific ward of hospital 1.

Fig. 3
figure 3

A minimum spanning tree based on wgMLST analysis of 9 OXA-48-positive Klebsiella pneumoniae in two hospitals. Isolates are represented by circles connected by branches proportional to the allelic distance, colours of the circle represent the origin location, the shading represent clonal clusters, branches and numbers represent allelic differences between isolates and ST numbers represent sequence types

Vancomycin-resistant Enterococci

Eighteen (0.7%) VRE isolates were detected in the hospitals and NHs. In total, 14 isolates were vanA-positive and four isolates contained vanB. Hospital 3 had three individual clusters of two isolates (ST80, ST203, ST117) on three different wards (Fig. 4). The cluster of five isolates with ST612 consisted of strains of hospital 2 from two different wards and NH 4. Hospital 2 and NH 4 are located within a range of 11–25 km.

Fig. 4
figure 4

A minimum spanning tree based on wgMLST analysis of 15 vancomycin resistant Enterococcus faecium isolates. Isolates are represented by circles connected by branches proportional to the allelic distance, colours of the circles represent the origin location, the shading represent clonal clusters, branches and numbers represent allelic differences between isolates and ST number represent sequence types

Discussion

The prevalence of ESBL-E was found to be 11.4% in participating hospitals and 14.0% in NHs. Notably, the NH data indicates a rising trend in ESBL-E carriage, with rates increasing from 6.2% in 2011 to 11.3% in 2015 [4]. This upward trajectory may also be observed in Belgian hospitals, as evidenced by a noteworthy increase in the occurrence of ESBL-E. coli isolates in clinical samples, escalating from 3.9% in 2005 to 10.1% in 2018 [5].

Among the detected ESBL genes, blaCTX-M-15 was the most prevalent, followed by other blaCTX-M-type genes. Rodriguez-Villalobos et al. reported an increased presence of blaCTX-M-type genes in clinical isolates in Belgium in 2008 (54%) compared to 2006 (23%) with 51.5% of the ESBL-E. coli isolates and 31% of the ESBL-K. pneumoniae isolates harbouring the blaCTX-M-15 gene in 2008 [6]. In this study, the percentage of E. coli containing the blaCTX-M-15 gene was also 51.5%. However, a remarkable increase could be noted in the K. pneumoniae isolates as blaCTX-M-15 was now present in 90.9% of ESBL-K. pneumoniae strains. The predominant ST type among the blaCTX-M-15 gene harbouring E. coli strains is ST131. This pandemic clone is pathogenic with an ability to cause frequent infections in the community, especially in the elderly, and with the remarkable feature of a low level of core genome recombinations [7], as also visible in the minimum spanning tree.

Unlike ESBL-E, CPE’s were not detected in the studied NHs, which is in line with previous Belgian studies conducted in NHs in 2011 [8] and 2015 [4]. In contrast to NHs, CPE could be detected in the hospitals with a mean prevalence of 1% and predominance of K. pneumoniae harbouring blaOXA-48 in line with previous findings [9]. The blaOXA-48 gene containing CPE isolates in our study (E. coli, K. pneumoniae, E. cloacae complex, Citrobacter spp., K. aerogenes and M. morgannii) all harboured the ca. 64kb IncL/M plasmid. The blaOXA-48 gene on the 64kb IncL/M plasmid is known to be widely distributed in Europe, especially in the Mediterranean area [10]. Therefore, focusing on the clonal spread of one bacterial species may be insufficient to control a plasmid-mediated outbreak [10], as illustrated by the different CPE species (K. pneumoniae, Citrobacter spp. and E. cloacae complex) isolated in ward 3 of hospital 1.

In total, 18 VRE isolates were detected with MIC values ranging from 12 to 256 mg/l. The detection of only two VRE in two different NHs is in accordance with the low prevalence described before [4], [11].

To conclude, our study shows, compared to previous studies, that the prevalence of MDRO in Flemish hospitals and NHs has increased over recent years and transmission of ESBL-E, CPE and/or VRE within hospital wards and in and between institutions could be observed. These results re-emphasize the importance of developing an inter-institutional collaboration for infection prevention in the healthcare network.