Original Research ARTICLE
Epidemiology of Extended-Spectrum β-Lactamase-Producing E. coli and Vancomycin-Resistant Enterococci in the Northern Dutch–German Cross-Border Region
- 1University of Groningen, University Medical Center Groningen, Department of Medical Microbiology, Groningen, Netherlands
- 2Laboratory of Medical Microbiology and Infectious Diseases, Isala Clinics, Zwolle, Netherlands
- 3Department of Microbiology, CERTE Medische Diagnostiek en Advies, Groningen, Netherlands
- 4Department of Intensive Care Medicine, Martini Hospital, Groningen, Netherlands
- 5Centre for Infectious Diseases Friesland, Izore, Leeuwarden, Netherlands
- 6Institute of Medical Microbiology, University Hospital Münster, Münster, Germany
- 7Institute of Hospital Hygiene, Klinikum Oldenburg, Oldenburg, Germany
Objectives: To reveal the prevalence and epidemiology of extended-spectrum β-lactamase (ESBL)- and/or plasmid AmpC (pAmpC)- and carbapenemase (CP) producing Enterobacteriaceae and vancomycin-resistant enterococci (VRE) across the Northern Dutch–German border region.
Methods: A point-prevalence study on ESBL/pAmpC/CP producing Enterobacteriaceae and VRE was carried out in hospitalized patients in the Northern Netherlands (n = 445, 2012–2013) and Germany (n = 242, 2012). Healthy individuals from the Dutch community (n = 400, 2010–2012) were also screened. In addition, a genome-wide gene-by-gene approach was applied to study the epidemiology of ESBL-Escherichia coli and VRE.
Results: A total of 34 isolates from 27 patients (6.1%) admitted to Dutch hospitals were ESBL/pAmpC positive and 29 ESBL-E. coli, three pAmpC-E. coli, one ESBL-Enterobacter cloacae, and one pAmpC-Proteus mirabilis were found. In the German hospital, 18 isolates (16 E. coli and 2 Klebsiella pneumoniae) from 17 patients (7.7%) were ESBL positive. In isolates from the hospitalized patients CTX-M-15 was the most frequently detected ESBL-gene. In the Dutch community, 11 individuals (2.75%) were ESBL/pAmpC positive: 10 ESBL-E. coli (CTX-M-1 being the most prevalent gene) and one pAmpC E. coli. Six Dutch (1.3%) and four German (3.9%) hospitalized patients were colonized with VRE. Genetic relatedness by core genome multi-locus sequence typing (cgMLST) was found between two ESBL-E. coli isolates from Dutch and German cross-border hospitals and between VRE isolates from different hospitals within the same region.
Conclusion: The prevalence of ESBL/pAmpC-Enterobacteriaceae was similar in hospitalized patients across the Dutch–German border region, whereas VRE prevalence was slightly higher on the German side. The overall prevalence of the studied pathogens was lower in the community than in hospitals in the Northern Netherlands. Cross-border transmission of ESBL-E. coli and VRE seems unlikely based on cgMLST analysis, however continuous monitoring is necessary to control their spread and stay informed about their epidemiology.
International travel and patient care are risk factors for dissemination of bacteria including multidrug-resistant microorganisms (MDRO), such as extended-spectrum β-lactamase (ESBL) and carbapenemase (CP)-producing Enterobacteriaceae (Woerther et al., 2013; Muller et al., 2015), and vancomycin-resistant enterococci (VRE). The prevalence of the latter has increased in the last years due to successful polyclonal subpopulations of hospital-associated (HA) Enterococcus faecium [previously designated clonal complex (CC) CC17] and which are also associated with amoxicillin resistance (ARE) (Arias and Murray, 2012). These populations are distinct from E. faecium isolates in the community and isolates from non-human sources (Galloway-Pena et al., 2012; Lebreton et al., 2013).
Netherlands and Germany as bordering countries with possible transfer of patients between them, created a cooperative network to prevent the spread of MDRO and to harmonize guidelines in healthcare settings1 (Muller et al., 2015). Surveillance studies to monitor the prevalence, resistance patterns, and molecular background of MDRO in hospitals and the community are essential to get insights into their epidemiology to implement infection prevention measures. Bacterial whole-genome sequencing (WGS) has been demonstrated to be very useful for epidemiological surveillance and detection of antimicrobial resistance (Kwong et al., 2015). The gene-by-gene approach uses a defined set of genes to extract an allele-based profile which makes it scalable and portable between laboratories (Maiden et al., 2013; Leopold et al., 2014). A core genome multi-locus sequence typing (cgMLST) scheme has been developed for E. faecium to distinguish between epidemiologically related and unrelated isolates (de Been et al., 2015). Although there is no cgMLST scheme nor threshold publicly approved yet for Escherichia coli, there are several tools available that allow to define an ad hoc cgMLST.
The aim of this study was to perform a point-prevalence study on ESBL/plasmid-mediated AmpC (pAmpC)/CP-producing Enterobacteriaceae and HA E. faecium (VRE and ARE) in hospitals in the Northern Dutch–German border region and to determine the predominant resistance genes. In addition, stool community samples from the Northern Netherlands were screened for the same resistant pathogens. A cgMLST was used to study hospital and cross-border dissemination of ESBL-E. coli and VRE.
Materials and Methods
A prospective point prevalence study was conducted in four of the largest hospitals (in total 3550 beds) in the Northern Netherlands between November 2012 and February 2013, covering a total population of approximately 2.85 million people. The Hospital Ethical Committee of the University Medical Center Groningen (UMCG) was informed and patients were approached to voluntarily participate in the study. Patients included in this study provided their written informed consent and a questionnaire concerning epidemiological and clinical data. The following high-risk wards for antibiotic-resistant microorganisms were selected: intensive care units (ICU), vascular surgery, internal medicine hematology/oncology, and dialysis wards (both for in- and outpatients). Gynecology and neurology (low-risk wards) were also included for comparison. From the largest German university hospital in the same (border) region, patients from four ICUs, a surgical ward, and a hematology/oncology ward were screened during October and November 2012 and included in the study. After consent agreement, all admitted patients from the studied wards were screened until completing a minimum of 100 samples per hospital.
The study in healthy people living in the Northern Netherlands was conducted retrospectively, using control patients included in a previous case–control study on microorganisms causing gastroenteritis. Control subjects were patients attending their general practitioner for a variety of medical questions, but no gastrointestinal problems, in the period between August 2010 and December 2012 (Bruijnesteijn van Coppenraet et al., 2015). No prevalence study was performed in the community in Germany.
A total of 445 rectal swabs (Copan ESwabTM) were taken from hospitalized patients (median age = 66 years, range 18–99 years) in the Northern Netherlands, 51.7% (n = 230) from men and 48.3% (n = 215) from women. A total of 328 (73.7%) patients were screened at high-risk wards and 117 (26.3%) patients were screened at low-risk wards (Table 1). In the German university hospital, 242 patients (median age = 64 years, range 0–94 years) were included, 64.5% (n = 156) men and 35.5% (n = 86) women. Of these 242 patients, 140 were screened only for ESBL, 22 only for VRE, and 80 for both. From the Dutch community study, 400 frozen faeces samples were included; 41% (n = 164) from men and 59% (n = 236) from women, 12% of the samples were from children. The median age of the healthy individuals was 47.5 years (range 0–84 years).
TABLE 1. Distribution of ESBL/pAmpC producing Enterobacteriaceae, and amoxicillin- and vancomycin-resistant E. faecium among the different wards in Dutch hospitals.
Microbiological Detection, Identification, and Susceptibility Testing
Dutch Hospitals and Retrospective Dutch Community Study
Rectal swabs (Dutch hospitalized patients) and approximately 50 μg of feces per sample (Dutch community patients) were enriched in selective broths: VRA broth containing BHI (brain heart infusion) with 20 mg/L amphoterin-B, 20 mg/L aztreonam, 20 mg/L colistin, and 16 mg/L amoxicillin and TSB-VC broth containing tryptic soy broth with 8 mg/L vancomycin and 0.25 mg/L cefotaxime. Both broths were incubated for 24 h at 35 ± 1°C. Subsequently, 10 μL of VRA broth was subcultured on VRE Brilliance agar (Oxoid®) and BMEG-2 agar (blood agar containing 64 mg/L meropenem, 2 mg/L gentamicin, 10 mg/L oxacillin, and 20 mg/L amphotericin-B) for identification of VRE and all ARE, respectively. Ten microliters of TSB-VC broth was subcultured onto ME/CF/CX comparted plates, containing iso-sensitest agar with 1 mg/L meropenem, 1 mg/L ceftazidime, or 1 mg/L cefotaxime, respectively, plus 20 mg/L vancomycin and 20 mg/L amphotericin-B (Mediaproducts, Groningen), for selection of ESBL/pAmpC/CP-producing bacteria. Plates were incubated for 24 h at 35 ± 1°C, except for VRE Brilliance agar plates that were incubated for 48 h.
Suspected colonies on VRE Brilliance, BMEG-2 and ME/CF/CX agar plates were streaked on blood agar (one isolate per morphotype). Species identification was done by matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) (Bruker Daltonik GmbH, Bremen). Confirmed Enterococcus spp. and Enterobacteriaceae spp., were tested for antibiotic susceptibility using VITEK®2 (bioMérieux) automatic system and EUCAST clinical breakpoints.
Rectal swabs were directly plated on chromID® ESBL agar (bioMérieux) for ESBL screening and enriched EnterococcoselTM Broth (Bile Esculin Azide Broth) (BD; Becton, Dickinson and Company) was used for VRE screening and subsequently cultured on chromID® VRE agar (bioMérieux).
Species identification and antibiotic susceptibility testing was done by MALDI-TOF MS (Bruker Daltonik GmbH, Bremen) and VITEK®2 (bioMérieux), respectively, following EUCAST criteria. Confirmation of ESBL was performed using disk diffusion (cefotaxime 30 μg, cefotaxime 30 μg plus clavulanic acid 10 μg, ceftazidime 30 μg, ceftazidime plus clavulanic acid 10 μg, cefepime 30 μg, cefepime 30 μg plus clavulanic acid 10 μg, and cefoxitin 30 μg) (Mast Diagnostics, Bootle, UK).
PCRs and Microarray
Enterococci isolates from Netherlands were screened by in-house PCR for IS16 (a marker for specific hospital-associated strains), vanA and vanB genes as described previously (Clark et al., 1993; Werner et al., 2011). The GenoType Enterococcus (Hain Lifescience GmbH) was used in enterococci isolates from Germany, which detects species and genotypes vanA, vanB, vanC1, and vanC2. ESBL and VRE positive isolates were sent to our hospital for further characterization.
Enterobacteriaceae isolates resistant to third generation cephalosporins and natural chromosomal AmpC producers intermediate or resistant to cefepime were selected for DNA extraction using the UltraClean Microbial DNA Isolation Kit (Mo Bio Laboratories, Inc.) and further characterized for the presence of ESBL/AmpC genes using a DNA-array (Check-MDR CT103, Check-points, Wageningen, Netherlands) (Garcia-Cobos et al., 2015).
Whole-Genome Sequencing of VRE and ESBL-E. coli
WGS was performed for all ESBL-E. coli and VRE isolates. For each isolate, several colonies (about 5 μL) of the culture were suspended in 300 μL microbead solution, which was subjected to DNA extraction with the Ultraclean Microbial DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA, United States). The DNA concentration and purity were measured using a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and the Qubit double-stranded DNA (dsDNA) HS and BR assay kits (Life Technologies, Carlsbad, CA, United States). One nanogram of bacterial DNA was used for library preparation. The DNA library was prepared using the Nextera XT library preparation kit with the Nextera XT v2 index kit (Illumina, San Diego, CA, United States). The library fragment length was aimed at fragments with a median size of 575 bases and was assessed with the Genomic DNA ScreenTape assay with the 2200 TapeStation system (Agilent Technologies, Waldbronn, Germany). Subsequently, the library was sequenced on a MiSeq sequencer, using the MiSeq reagent kit v2 generating 250-bp paired-end reads. Sequencing was aimed at a coverage of at least 60-fold. MiSeq data were processed with MiSeq control software v188.8.131.52 and MiSeq Reporter v2.4 (Illumina, San Diego, CA, United States). Reads were quality-trimmed using the CLC Genomics Workbench software version 9.0.1 (CLC bio, Aarhus, Denmark) using default settings except for the following modifications: “trim using quality scores was set to 0.02” and “discard reads below length was set to 15.” Subsequently, trimmed-reads were de novo assembled with an optimal word size of 29 and a minimum contig length of 500. Metrics on raw read and assembly level are provided in Supplementary Table S1.
cgMLST of VRE and ESBL-E. coli
A genome-wide gene-by-gene comparison approach was used to determine the genetic relatedness using SeqSphere+ version 3.4.0 (Ridom GmbH, Münster, Germany) (Leopold et al., 2014).
Genome assemblies from the VRE isolates were analyzed using the E. faecium cgMLST scheme previously published, considering a cluster alert distance of 20 different alleles (de Been et al., 2015).
An ad hoc cgMLST and whole-genome MLST (wgMLST) scheme was determined for E. coli isolates using the MLST+ target definer function with default parameters2 and E. coli K-12 as a reference (GenBank accession no. NC_010473.1). The filters applied to reference genome were: “minimum length filter” that discards genes shorter than 50 bases; “start codon filter” that discards all genes that contain no start codon at the beginning of the gene: “stop codon filter” that discards all genes that contain no stop codon, more than one stop codon or if the stop codon is not at the end of the gene: “homologous gene filter” that discards all genes that have fragments that occur in multiple copies in a genome (with identity ≥90% and more than 100 bases overlap); “gene overlap filter” that discards the shorter of two overlapping flanking genes if these genes overlap >4 bp. The remaining genes were then used in a pairwise comparison using BLAST (Leopold et al., 2014) with 45 query genomes (Supplementary Table S2a). All genes of the reference genome that were common in all query genomes with a sequence identity of ≥90 and 100 overlap, and with the default parameter stop codon percentage filter turned on, formed the final cgMLST scheme; this discards all genes that have internal stop codons in >20% of the query genomes. Additionally, 26 plasmid sequences (Supplementary Table S2b) were added to exclude such genes are part of the cgMLST typing scheme. The final cgMLST scheme consisted of 1.771 targets/genes, and 2329 accessory genes were additionally included for the wgMLST scheme (Supplementary Tables S3, S4). The minimum coverage of the genome assemblies was 20 times (Supplementary Table S1) and the percentage of good targets/genes included in the cgMLST were 97.6% for E. coli and 98.6% for E. faecium (Supplementary Tables S5, S6).
Furthermore, to determine the genetic relatedness, the genetic distance for the E. coli isolates was calculated as the proportion of allele differences: dividing the number of allele differences between two genomes by the total number of genes commonly shared by those two genomes (Kluytmans-van den Bergh et al., 2016). In this study, thresholds for genetic distance were described to discriminate between epidemiologically related and unrelated E. coli isolates as 0.0095 when using wgMLST and 0.0105 for cgMLST.
Escherichia coli STs were determined uploading genome assemblies to SeqSphere+ software following the scheme of Wirth et al. (2006). Sequence genomes with no conclusive results for the 7-gene MLST were uploaded to the Enterobase database3. Additionally, E. coli major phylogenetic groups (A, B1, B2, and D) were analyzed in silico by using MLST+ Target Definer function of SeqSphere+, including the chuA, yjaA, and TSPE4.C2 loci (Clermont et al., 2000).
Genome assemblies were also uploaded to the Center for Genomic Epidemiology to extract information on resistance genes (ResFinder) and virulence factors (VirulenceFinder), and species confirmation for VRE and ESBL-E. coli (KmerFinder), and serotype (SerotypeFinder) and plasmid replicons (PlasmidFinder) for ESBL-E. coli (Zankari et al., 2012; Carattoli et al., 2014; Hasman et al., 2014; Joensen et al., 2014, 2015; Larsen et al., 2014).
In the Dutch hospital prevalence study, associations between ESBL and ARE carriage and the following variables were analyzed: length of hospital stay, antibiotic use, and (low- or high-risk) ward. Information was gathered by the questionnaires. Statistical analyses were performed using SPSS for Windows, v. 20.0. Univariate analyses were performed using the Fisher’s exact or Chi-square methods for categorical variables. The Mann–Whitney U-test was used as a non-parametric tests in variables with no normal distribution. Results with a p-value of ≤0.05 were considered to be statistically significant. All p-values are two-tailed.
Thirty-four isolates from 27 of the 445 included patients admitted to hospitals in the Northern Netherlands (6.1%) were confirmed ESBL and/or pAmpC positive. A total of 85.2% (23/27), 14.8% (4/27), and 3.7% (1/27) of these patients were positive for ESBL, pAmpC, and both, respectively. Among the 34 isolates, 32 were E. coli, of which 29 were ESBL positive and three were pAmpC producers. Resistance genes detected in the E. coli isolates are shown in Table 2. CTX-M-15 (n = 8) and CTX-M-14 (n = 8) were the most prevalent ones. The other two isolates were an Enterobacter cloacae, containing a CTX-M-1-like gene and a pAmpC CMY-II producing Proteus mirabilis. At high-risk wards, 19 patients (5.8%) were found with ESBL/pAmpC isolates compared to eight patients (6.8%) at low-risk wards (p = 0.68; NS). No association was found between ESBL/pAmpC carriage and antibiotic use, length of hospital stay or ward (Table 1).
TABLE 2. Molecular characterization of the E. coli isolates from the community and hospital patients in the Netherlands and Germany.
In the German hospital, a total of 18 isolates from 17 patients (17/220; 7.7%) were ESBL positive. Sixteen isolates were E. coli and two were Klebsiella pneumoniae. Of these, 12 E. coli and one K. pneumoniae isolates were available for molecular testing. Six out of 12 (50%) E. coli isolates and the K. pneumoniae isolate had a CTX-M-15 gene (Table 2).
In the retrospective Dutch community study, 11 patients (11/400; 2.75%) were ESBL/pAmpC positive: 10 ESBL-E. coli (CTX-M-1 being the most prevalent gene) and one pAmpC E. coli (Table 2). Overall, no carbapenem resistance was observed neither in the community nor in the hospitals.
Escherichia coli MLST and Phylogenetic Groups
Among ESBL/pAmpC-E. coli isolates from Dutch hospitals, the most prevalent ST was ST131 (CC, ST131; n = 5, 15.6%), all of them belonging to phylogroup B2 (Table 2). In the Dutch community isolates, 10 different STs were found, most of them belonging to CC ST10 (n = 3, 27.3%) and one isolate to ST131 (phylogroup B2). In the German hospital, the most prevalent STs were ST38 (33.3%) and ST10 (33.3%) (Table 2).
Amoxicillin- and Vancomycin-Resistant E. faecium (ARE and VRE)
In the Dutch hospitals, 105 patients (105/445; 23.6%) were colonized with ARE, including six patients (6/445; 1.3%) with VRE. All ARE were positive for IS16 and all VRE were vanB positive. Colonization of ARE (and VRE) was associated with high-risk wards (p < 0.001), prolonged hospitalization (p < 0.001), and use of antibiotics (p = 0.05), especially penicillins and fluoroquinolones (p < 0.001) (Table 3).
TABLE 3. Variables associated with carriage of amoxicillin-resistant E. faecium (ARE) and extended-spectrum β-lactamase (ESBL)- and/or plasmid AmpC (pAmpC)- producing Enterobacteriaecae.
In the border German university hospital four (4/102; 3.9%) VRE isolates were isolated. Three of them were vanA positive and one was vanB positive.
In the retrospective Dutch community study, six ARE (6/400; 1.5%) were found, three of them were IS16 positive. Only one vanA-VRE (1/400; 0.25%) was found, this strain was amoxicillin susceptible and IS16 negative.
cgMLST and wgMLST Comparison of ESBL-E. coli Isolates from the Community and Hospitals
Genome assemblies of 55 ESBL-E. coli [Dutch community (n = 11), Dutch hospitals (n = 32), and German hospital (n = 12)] of this study were analyzed by a gene-by-gene approach and the allelic distance from the cgMLST and wgMLST were visualized in a minimum spanning tree (Figure 1 and Supplementary Figure S1, respectively).
FIGURE 1. Minimum spanning tree of ESBL-E. coli isolates from hospitals and the community. Distance based on a cgMLST of 1771 genes using the parameters “pairwise ignoring missing values” during calculation. Each circle represents a genotype, colors indicate geographical origin and community or hospital. Orange: hospital—Netherlands; blue: hospital—Germany; green: community—Netherlands. Number of different alleles are indicated on the edges between connected isolates (nodes). The cut-off values for defining a group was 35 alleles. Isolates are presented by their ID and ST.
Six groups of isolates with a lower number of different alleles (≤35) by cgMLST were further analyzed. Supplementary Table S7 summarizes the origin of the isolates in every group and the core and whole-genome genetic distance. Those groups formed by isolates with an epidemiological link (isolated from the same patient; group 1, 4, 5a, 6a, and 7), showed a core and whole-genome genetic distance lower than 0.0030 and 0.0046, respectively. In addition, isolates of group 5b, although with unknown epidemiological link, had a core genetic distance of 0.0063 and a whole-genome genetic distance of 0.0076. Both isolates were positive for CTX-M-14, however, no plasmid replicons were found in one of them (isolate 38_Esco_HA-DE) (Supplementary Table S7).
Among those groups including isolates with non (or unknown) epidemiological link, the core genome genetic distance was between 0.0122 and 0.0199 and the whole-genome genetic distance was between 0.0104 and 0.0208 (groups 2, 3, 6b, and 6c; Figure 1). Resistance and virulence profiles of the isolates are shown in Supplementary Table S8.
cgMLST Comparison of VRE Isolates from the Community and Hospitals
A minimum spanning tree was created for the 11 VRE isolates [Dutch community (n = 1), Dutch hospitals (n = 6), and German hospital (n = 4)]. Two clusters of isolates from different patients were observed (Figure 2). One cluster of four vanB-VRE isolates from the Dutch hospital belonged to cluster type (CT) 110 (ST17); two isolates were from the same ward in hospital A and the other two isolates were isolated from different wards in hospital B. The other cluster of two vanA-VRE isolates were isolated from different wards from the German hospital (CT 20, ST203). The resistance and virulence genotypes of VRE isolates are shown in Supplementary Table S8.
FIGURE 2. Minimum spanning tree of VREfcm, cgMLST based on 1423 genes using the parameters “pairwise ignoring missing values” during distance calculation. Each circle represents a genotype and colors indicate cluster types (CT). Number of different alleles are indicated on the edges between connected isolates (nodes). Isolates are presented by their ID, ST, and CT.
Nucleotide Sequence Accession Number
Sequence data obtained in this study has been deposited at the National Center for Biotechnology Information under BioProject no. PRJNA352198.
This study shows the molecular epidemiology of ESBL/pAmpC-E. coli and HA E. faecium in hospitals in the Northern Dutch–German border region and the community in the Northern Netherlands. Dutch hospitals showed a prevalence for ESBL/pAmpC – Enterobacteriaceae, VRE, and ARE of 6.1, 1.3, and 23.6%, respectively, whereas the prevalence in the community was 2.75, 0.25, and 1.5%, respectively. The German hospital had an ESBL/pAmpC prevalence of 7.7 and 3.9% for VRE.
A previous study reported a prevalence of ESBL-producing bacteria of 4.9% in the Netherlands (Overdevest et al., 2014), comparable to the 6.1% prevalence observed in Dutch hospitals in this study. A prevalence of 5.6% ESBL-producing E. coli isolates in hospitalized and ambulatory patients in Germany has been reported recently (Pietsch et al., 2015), which is slightly lower than the 7.7% observed in the present study.
Furthermore, we observed an ESBL-E. coli prevalence of 2.5% in the Northern Netherlands community, which is low compared to previous studies in other regions, in which the prevalence in the community ranged from 4.7% (2009) to 10.1% (2011) (Reuland et al., 2013; van Hoek et al., 2015). This difference may have several reasons. First, ESBL prevalence may vary between regions and over time, and natural eradication of resistant Enterobacteriaceae might occur over time in the community (Bar-Yoseph et al., 2016). Additionally, samples included in this study were only chosen from patients without any gastrointestinal complaints, a factor which otherwise has been described to be associated with high ESBL prevalence (Reuland et al., 2013).
The majority of the resistance genes found in our community isolates were CTX-M-1 which is broadly disseminated among animals in Europe, especially in cattle and pigs, followed by the CTX-M-15 gene, commonly associated with human origin (Reuland et al., 2013; Pietsch et al., 2015). The latter was the most frequent gene among the Dutch and German hospital isolates, in concordance with previous studies (Ewers et al., 2012; Reuland et al., 2013; Pietsch et al., 2015).
The pAmpC prevalence in E. coli in our study was 0.3%, comparable to the prevalence of 0.6% what was reported in the study of van Hoek et al. (2015) (0.6% pAmpC Enterobacteriaceae) and somewhat lower to findings of Reuland et al. (2015) (1.3% pAmpC-E. coli). The most common pAmpC gene found in hospital and community isolates were CMY-II, which is together with DHA frequently found in human isolates (Reuland et al., 2015).
ESBL-producing E. coli belonging to CC ST131-phylogroup B2 are usually associated with more virulent strains (Overdevest et al., 2015). These were frequently found in the Dutch hospitals included in the present study but only sporadically in the community samples. This CC ST131-phylogroup B2 was also prevalent in a study carried out in hospitals in the Rotterdam region (van der Bij et al., 2011). CC ST10 was predominant among the ESBL-producing E. coli in the community, the same CC was also described to be prevalent in another Dutch study in community patients (Reuland et al., 2013).
We observed an overall ARE and VRE prevalence in hospitalized patients of 23.6 and 1.3%, respectively. Similar observations were made in a study performed in Dutch hospitals in 2008 reporting ARE carriage rates of 10–16% upon admission and 15–39% on acquisition in hematology and gastroenterology/nephrology wards (de Regt et al., 2008). The clinical significance of enterococcal infections and active VRE screening has been a matter of discussion. However, in immunocompromised patients, high morbidity and mortality rates have been reported in infections caused by enterococci (Zhou et al., 2013). In this study, ARE/VRE carriage was associated with prolonged hospitalization and antibiotic use, which is in line with previous literature (Torell et al., 1999). We found a high carriage rate of ARE in high-risk wards (30.2%). Notably, these patients may be at risk for a subsequent infection. Since 2011, VRE started to become a problem in multiple hospitals in the Netherlands: a total of 14 hospitals were affected with outbreaks of VRE in October 2012 (van den Brink, 2013). However, in this study a prevalence of VRE (vanB) carriage of only 1.3% was found. This is probably due to extensive infection prevention measures and successful outbreak management control. The prevalence of 1.3% is similar to what has been previously published in the Netherlands, with prevalence rates ranging from 1.4 to 2% in the 1990s (Endtz et al., 1997; van den Braak et al., 2000). The VRE prevalence in the German hospital was slightly higher (3.9%), though it is known that Germany has a higher VRE prevalence compared to the Netherlands4.
In our Dutch community one vanA-VRE was found, that was amoxicillin susceptible and IS16 negative, indicative for a non-hospital origin (Galloway-Pena et al., 2012; Lebreton et al., 2013). Endtz et al. (1997) reported a higher number of VRE in the community (2%), however, this study did not include information about amoxicillin resistance nor IS16 which makes it difficult to determine if they had a hospital or non-hospital origin (Galloway-Pena et al., 2012; Lebreton et al., 2013).
The cgMLST analysis in our study showed heterogeneity among E. coli species, and isolates were genetically distributed independently of their origin. The hospital and community ESBL-E. coli isolates included in this study did not show any genetic relatedness except for the ones isolated from the same patient and for two isolates (group 5b) from patients in different hospitals across the Dutch–German border, in a distance of approximately 200 km and with no known epidemiological link. The patient from the Dutch hospital was a dialyses outpatient (isolation date December 2012) whereas the patient from the German hospital was admitted to ICU (isolation date November 2012). Interestingly, both isolates harbored the same ESBL gene and virulence factors.
Genetic relatedness was found between four VRE isolates (CT110) from patients from two different Dutch hospitals (Figure 2), which indicates transmission between wards, but also between hospitals in a close geographical region similar to findings of a previous population-based study of VRE using WGS that also showed intra- and inter-regional spread of closely related VRE isolates (Pinholt et al., 2015). Although no genetic relatedness was found between VRE isolates of Dutch and German hospitals, the numbers of VRE isolates were too low to draw definite conclusions. It is known that several VRE CT co-circulate in Germany and the Netherlands (data not shown). However, only some laboratories have implemented the use of cgMLST in their routine to analyze VRE outbreaks and more epidemiological studies are needed to investigate cross-border transmission of VRE.
To our knowledge, there are no similar studies that compare and investigate the molecular epidemiology of ESBL-E. coli and VRE in hospitals and the community by WGS. Recently, the same approach has been used to study the clonality of ESBL-producing Enterobacteriaceae from environmental and stool samples from farmers suggesting possible cross-transmission between the farmers and the environment. This was only based on number of allele differences (Fischer et al., 2016) which makes it difficult to interpret results without considering the total number of genes included in the cgMLST scheme. In our study, we determined the genetic relatedness between ESBL-E. coli using cgMLST or wgMLST comparison and genetic distance calculation. These results were in concordance with the genetic distance thresholds of 0.0095 (wgMLST) and 0.0105 (cgMLST) previously established for E. coli based on known existing epidemiological links by analyzing more than 2000 ESBL-Enterobacteriaceae isolates from Dutch hospitals (Kluytmans-van den Bergh et al., 2016).
In another study, a cgMLST approach for several MDR bacteria was prospectively used for taking relevant infection control decisions in a hospital setting (Mellmann et al., 2016). A threshold of >10 differing alleles was defined to exclude nosocomial transmission of MDR E. coli (Mellmann et al., 2016). If we would have applied this threshold we would have missed the genetic relatedness between isolates belonging to group 5b, presenting 11 different alleles (Figure 1 and Supplementary Table S7). This highlights that thresholds based on number of allele differences are only applicable to specific collections within a study, whereas the genetic distance calculation seems to give a more objective result, independently of the analyzed population.
We acknowledge this study has some limitations. No community study in the German cross-border region, neither ARE monitoring in the German hospital were performed. Laboratory methods for isolation of ESBL Enterobacteriaceae and VRE differed between Dutch and German hospitals since no enrichment broth was used in Germany, however, selective media agar was used in both regions. Since this study was anonymous, some epidemiological data were not available which makes it more difficult to draw conclusions regarding genetic relatedness among isolates between patients.
In conclusion, the results of this study suggest that ESBL/pAmpC-E. coli circulate in the hospital and the community, although a higher prevalence of ESBL/pAmpC-E. coli was observed in hospitals compared to the community in the Northern Netherlands. Hospitals in the Northern Dutch–German region showed a similar prevalence of ESBL/pAmpC-Enterobacteriaceae. VRE prevalence was still low in the hospital as well as in the community in the Northern Netherlands. The German hospital showed a slightly higher VRE prevalence compared to hospitals in the Northern Netherlands. Nosocomial but no cross-border transmission of VRE was observed in this study. Epidemiologically related ESBL-E. coli and VRE were uncommon across the Northern Dutch–German border in the studied population, as only two ESBL-E. coli isolates from a Dutch and a German hospital were genetically similar. Cooperation between bordering countries and continuous monitoring using high discriminatory typing methods are still necessary to keep the epidemiology of resistant pathogens updated thereby helping to control their spread.
These results were partially presented at the ECCMID conference 2016, Amsterdam.
The “METc (medical ethical committee) UMCG” believes that this research is not a research involving humans as is ment in the Law on Medical Scientific Research involving Human Beings (WMO). Therefore, the METc UMCG has decided that no WMO approval is needed.
GR, DB, LM, NH, TS, LB, JW, JvZ, and RK have all contributed to the collection of all samples, they have all reviewed the article carefully. XZ, SG-C, GK, JA, DB, JR, and AF have contributed to the study design and reviewed the article. XZ and SG-C wrote the manuscript together and analyzed the data. JR and AF reviewed the manuscript several times, before sending to the other authors for review.
This study was supported by the Interreg IVa-funded projects EurSafety Health-net (III-1-02=73) and SafeGuard (III-2-03=025), part of a Dutch–German cross-border network supported by the European Commission, the German Federal States of Nordrhein-Westfalen and Niedersachsen, and the Dutch Provinces of Overijssel, Gelderland, and Limburg.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Dr. C. R. C. Doorenbos, M. Zigterman, W. Postma, P. Rurenga, N. Welles, A. Woudstra, S. de Vries, and A. J. Stellingwerf for their participation in this study.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.01914/full#supplementary-material
FIGURE S1 | Minimum spanning tree of ESBL-E. coli isolates from hospitals and the community. Distance based on a wgMLST of 4100 genes (cgMLST of 1771 genes and 2329 accessory genes) using the parameters “pairwise ignoring missing values” during calculation. Each circle represents a genotype, colors indicate geographical origin and community or hospital. Orange: hospital— Netherlands; blue: hospital—Germany; green: community—Netherlands. Number of different alleles are indicated on the edges between connected isolates (nodes). The same groups considered by cgMLST analysis are highlighted. Isolates are presented by their ID and ST.
TABLE S1 | Metrics on raw read data and assemblies.
TABLE S2a | Finished E. coli query genomes used in this study to develop and ad hoc cgMLST scheme (n = 45). One representative isolate of every ST from every collection community NL (n = 10), Dutch hospitals (n = 20) and German hospital (n = 6) of the present study and 9 E. coli genomes from Dutch patients and farmers previously published (de Been et al., 2014).
TABLE S2b | Finished plasmid genomes for exclusion of genes with BLAST matches >90% and >100 bp length found within the query sequences used in this study to develop a cgMLST scheme.
TABLE S3 | Genes included in the cgMLST scheme of E. coli.
TABLE S4 | Accessory genes included in the wgMLST scheme of E. coli.
TABLE S5 |E. coli cgMLST allele types for distance calculation and percentage of good targets/genes.
TABLE S6 |E. faecium cgMLST allele types for distance calculation and percentage of good targets/genes.
TABLE S7 | Genetic distance for pairwise comparisons of grouped ESBL- E. coli isolates.
TABLE S8 | Results of ResFinder, VirulenceFinder, PlasmidFinder, and SerotypeFinder for E. coli and E. faecium isolates.
- ^Antimicrobial Resistance Interactive Database (EARS-Net) (http://ecdc.Europa.eu/en/)
Bar-Yoseph, H., Hussein, K., Braun, E., and Paul, M. (2016). Natural history and decolonization strategies for ESBL/carbapenem-resistant Enterobacteriaceae carriage: systematic review and meta-analysis. J. Antimicrob. Chemother. 71, 2729–2739. doi: 10.1093/jac/dkw221
Bruijnesteijn van Coppenraet, L. E., Dullaert-de Boer, M., Ruijs, G. J., van der Reijden, W. A., van der Zanden, A. G., Weel, J. F., et al. (2015). Case-control comparison of bacterial and protozoan microorganisms associated with gastroenteritis: application of molecular detection. Clin. Microbiol. Infect. 21, 592.e9–e19. doi: 10.1016/j.cmi.2015.02.007
Carattoli, A., Zankari, E., Garcia-Fernandez, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 58, 3895–3903. doi: 10.1128/AAC.02412-14
Clark, N. C., Cooksey, R. C., Hill, B. C., Swenson, J. M., and Tenover, F. C. (1993). Characterization of glycopeptide-resistant enterococci from U.S. hospitals. Antimicrob. Agents Chemother. 37, 2311–2317. doi: 10.1128/AAC.37.11.2311
Clermont, O., Bonacorsi, S., and Bingen, E. (2000). Rapid and simple determination of the Escherichia coli phylogenetic group. Appl. Environ. Microbiol. 66, 4555–4558. doi: 10.1128/AEM.66.10.4555-4558.2000
de Been, M., Lanza, V. F., de Toro, M., Scharringa, J., Dohmen, W., Du, Y., et al. (2014). Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages. PLOS Genet. 10:e1004776. doi: 10.1371/journal.pgen.1004776
de Been, M., Pinholt, M., Top, J., Bletz, S., Mellmann, A., van Schaik, W., et al. (2015). A core genome MLST scheme for high-resolution typing of Enterococcus faecium. J. Clin. Microbiol. 53, 3788–3797. doi: 10.1128/JCM.01946-15
de Regt, M. J., van der Wagen, L. E., Top, J., Blok, H. E., Hopmans, T. E., Dekker, A. W., et al. (2008). High acquisition and environmental contamination rates of CC17 ampicillin-resistant Enterococcus faecium in a dutch hospital. J. Antimicrob. Chemother. 62, 1401–1406. doi: 10.1093/jac/dkn390
Endtz, H. P., van den Braak, N., van Belkum, A., Kluytmans, J. A., Koeleman, J. G., Spanjaard, L., et al. (1997). Fecal carriage of vancomycin-resistant enterococci in hospitalized patients and those living in the community in the Netherlands. J. Clin. Microbiol. 35, 3026–3031.
Ewers, C., Bethe, A., Semmler, T., Guenther, S., and Wieler, L. H. (2012). Extended-spectrum beta-lactamase-producing and AmpC-producing Escherichia coli from livestock and companion animals, and their putative impact on public health: a global perspective. Clin. Microbiol. Infect. 18, 646–655. doi: 10.1111/j.1469-0691.2012.03850.x
Fischer, J., Hille, K., Ruddat, I., Mellmann, A., Kock, R., and Kreienbrock, L. (2016). Simultaneous occurrence of MRSA and ESBL-producing Enterobacteriaceae on pig farms and in nasal and stool samples from farmers. Vet. Microbiol. 200, 107–113. doi: 10.1016/j.vetmic.2016.05.021
Galloway-Pena, J., Roh, J. H., Latorre, M., Qin, X., and Murray, B. E. (2012). Genomic and SNP analyses demonstrate a distant separation of the hospital and community-associated clades of Enterococcus faecium. PLOS ONE 7:e30187. doi: 10.1371/journal.pone.0030187
Garcia-Cobos, S., Kock, R., Mellmann, A., Frenzel, J., Friedrich, A. W., and Rossen, J. W. (2015). Molecular typing of Enterobacteriaceae from pig holdings in north-western Germany reveals extended- spectrum and AmpC beta-lactamases producing but no carbapenem resistant ones. PLOS ONE 10:e0134533. doi: 10.1371/journal.pone.0134533
Hasman, H., Saputra, D., Sicheritz-Ponten, T., Lund, O., Svendsen, C. A., Frimodt-Moller, N., et al. (2014). Rapid whole-genome sequencing for detection and characterization of microorganisms directly from clinical samples. J. Clin. Microbiol. 52, 139–146. doi: 10.1128/JCM.02452-13
Joensen, K. G., Scheutz, F., Lund, O., Hasman, H., Kaas, R. S., Nielsen, E. M., et al. (2014). Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J. Clin. Microbiol. 52, 1501–1510. doi: 10.1128/JCM.03617-13
Joensen, K. G., Tetzschner, A. M., Iguchi, A., Aarestrup, F. M., and Scheutz, F. (2015). Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J. Clin. Microbiol. 53, 2410–2426. doi: 10.1128/JCM.00008-15
Kluytmans-van den Bergh, M. F., Rossen, J. W., Bruijning-Verhagen, P. C., Bonten, M. J., Friedrich, A. W., Vandenbroucke-Grauls, C. M., et al. (2016). Whole genome multilocus sequence typing of extended-spectrum beta-lactamase-producing Enterobacteriaceae. J. Clin. Microbiol. 54, 2919–2927. doi: 10.1128/JCM.01648-16
Larsen, M. V., Cosentino, S., Lukjancenko, O., Saputra, D., Rasmussen, S., Hasman, H., et al. (2014). Benchmarking of methods for genomic taxonomy. J. Clin. Microbiol. 52, 1529–1539. doi: 10.1128/JCM.02981-13
Lebreton, F., van Schaik, W., McGuire, A. M., Godfrey, P., Griggs, A., Mazumdar, V., et al. (2013). Emergence of epidemic multidrug-resistant Enterococcus faecium from animal and commensal strains. mBio 4:e00534-13. doi: 10.1128/mBio.00534-13
Leopold, S. R., Goering, R. V., Witten, A., Harmsen, D., and Mellmann, A. (2014). Bacterial whole-genome sequencing revisited: portable, scalable, and standardized analysis for typing and detection of virulence and antibiotic resistance genes. J. Clin. Microbiol. 52, 2365–2370. doi: 10.1128/JCM.00262-14
Maiden, M. C., Jansen van Rensburg, M. J., Bray, J. E., Earle, S. G., Ford, S. A., Jolley, K. A., et al. (2013). MLST revisited: the gene-by-gene approach to bacterial genomics. Nat. Rev. Microbiol. 11, 728–736. doi: 10.1038/nrmicro3093
Mellmann, A., Bletz, S., Boking, T., Kipp, F., Becker, K., Schultes, A., et al. (2016). Real-time genome sequencing of resistant bacteria provides precision infection control in an institutional setting. J. Clin. Microbiol. 54, 2874–2881. doi: 10.1128/JCM.00790-16
Muller, J., Voss, A., Kock, R., Sinha, B., Rossen, J. W., Kaase, M., et al. (2015). Cross-border comparison of the dutch and German guidelines on multidrug-resistant gram-negative microorganisms. Antimicrob. Res. Infect. Control 4:7. doi: 10.1186/s13756-015-0047-6
Overdevest, I. T., Bergmans, A. M., Verweij, J. J., Vissers, J., Bax, N., Snelders, E., et al. (2015). Prevalence of phylogroups and O25/ST131 in susceptible and extended-spectrum beta-lactamase-producing Escherichia coli isolates, the Netherlands. Clin. Microbiol. Infec. 21, 570.e1–e4. doi: 10.1016/j.cmi.2015.02.020
Overdevest, I. T., Heck, M., van der Zwaluw, K., Huijsdens, X., van Santen, M., Rijnsburger, M., et al. (2014). Extended-spectrum beta-lactamase producing Klebsiella spp. in chicken meat and humans: a comparison of typing methods. Clin. Microbiol. Infec. 20, 251–255. doi: 10.1111/1469-0691.12277
Pietsch, M., Eller, C., Wendt, C., Holfelder, M., Falgenhauer, L., Fruth, A., et al. (2015). Molecular characterisation of extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli isolates from hospital and ambulatory patients in Germany. Vet. Microbiol. 200, 130–137. doi: 10.1016/j.vetmic.2015.11.028
Pinholt, M., Larner-Svensson, H., Littauer, P., Moser, C. E., Pedersen, M., Lemming, L. E., et al. (2015). Multiple hospital outbreaks of vanA Enterococcus faecium in Denmark, 2012-13, investigated by WGS, MLST and PFGE. J. Antimicrob. Chemother. 70, 2474–2482. doi: 10.1093/jac/dkv142
Reuland, E. A., Halaby, T., Hays, J. P., de Jongh, D. M., Snetselaar, H. D., van Keulen, M., et al. (2015). Plasmid-mediated AmpC: prevalence in community-acquired isolates in Amsterdam, the Netherlands, and risk factors for carriage. PLOS ONE 10:e0113033. doi: 10.1371/journal.pone.0113033
Reuland, E. A., Overdevest, I. T., Al Naiemi, N., Kalpoe, J. S., Rijnsburger, M. C., Raadsen, S. A., et al. (2013). High prevalence of ESBL-producing Enterobacteriaceae carriage in dutch community patients with gastrointestinal complaints. Clin. Microbiol. Infect. 19, 542–549. doi: 10.1111/j.1469-0691.2012.03947.x
Torell, E., Cars, O., Olsson-Liljequist, B., Hoffman, B. M., Lindback, J., and Burman, L. G. (1999). Near absence of vancomycin-resistant enterococci but high carriage rates of quinolone-resistant ampicillin-resistant enterococci among hospitalized patients and nonhospitalized individuals in Sweden. J. Clin. Microbiol. 37, 3509–3513.
van den Braak, N., Ott, A., van Belkum, A., Kluytmans, J. A., Koeleman, J. G., Spanjaard, L., et al. (2000). Prevalence and determinants of fecal colonization with vancomycin-resistant Enterococcus in hospitalized patients in the Netherlands. Infect. Control Hosp. Epidemiol. 21, 520–524. doi: 10.1086/501797
van der Bij, A. K., Peirano, G., Goessens, W. H., van der Vorm, E. R., van Westreenen, M., and Pitout, J. D. (2011). Clinical and molecular characteristics of extended-spectrum-beta-lactamase-producing Escherichia coli causing bacteremia in the rotterdam area, Netherlands. Antimicrob. Agents Chemother. 55, 3576–3578. doi: 10.1128/AAC.00074-11
van Hoek, A. H., Schouls, L., van Santen, M. G., Florijn, A., de Greeff, S. C., and van Duijkeren, E. (2015). Molecular characteristics of extended-spectrum cephalosporin-resistant Enterobacteriaceae from humans in the community. PLOS ONE 10:e0129085. doi: 10.1371/journal.pone.0129085
Werner, G., Fleige, C., Geringer, U., van Schaik, W., Klare, I., and Witte, W. (2011). IS element IS16 as a molecular screening tool to identify hospital-associated strains of Enterococcus faecium. BMC Infect. Dis. 11:80. doi: 10.1186/1471-2334-11-80
Wirth, T., Falush, D., Lan, R., Colles, F., Mensa, P., Wieler, L. H., et al. (2006). Sex and virulence in Escherichia coli: an evolutionary perspective. Mol. Microbiol. 60, 1136–1151.. doi: 10.1111/j.1365-2958.2006.05172.x
Woerther, P. L., Burdet, C., Chachaty, E., and Andremont, A. (2013). Trends in human fecal carriage of extended-spectrum beta-lactamases in the community: toward the globalization of CTX-M. Clin. Microbiol. Rev. 26, 744–758. doi: 10.1128/CMR.00023-13
Zankari, E., Hasman, H., Cosentino, S., Vestergaard, M., Rasmussen, S., Lund, O., et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 67, 2640–2644. doi: 10.1093/jac/dks261
Zhou, X., Arends, J. P., Span, L. F., and Friedrich, A. W. (2013). Algorithm for pre-emptive glycopeptide treatment in patients with haematologic malignancies and an Enterococcus faecium bloodstream infection. Antimicrob. Resist. Infect. Control 2:24. doi: 10.1186/2047-2994-2-24
Keywords: WGS, cgMLST, VRE, ESBL, hospital, community, prevalence, cross-border research
Citation: Zhou X, García-Cobos S, Ruijs GJHM, Kampinga GA, Arends JP, Borst DM, Möller LV, Holman ND, Schuurs TA, Bruijnesteijn van Coppenraet LE, Weel JF, van Zeijl JH, Köck R, Rossen JWA and Friedrich AW (2017) Epidemiology of Extended-Spectrum β-Lactamase-Producing E. coli and Vancomycin-Resistant Enterococci in the Northern Dutch–German Cross-Border Region. Front. Microbiol. 8:1914. doi: 10.3389/fmicb.2017.01914
Received: 24 April 2017; Accepted: 20 September 2017;
Published: 05 October 2017.
Edited by:John R. Battista, Louisiana State University, United States
Reviewed by:Zhuofei Xu, Huazhong Agricultural University, China
Jose M. Munita, University of Texas Health Science Center at Houston, United States
Copyright © 2017 Zhou, García-Cobos, Ruijs, Kampinga, Arends, Borst, Möller, Holman, Schuurs, Bruijnesteijn van Coppenraet, Weel, van Zeijl, Köck, Rossen and Friedrich. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work.