Whole genome analyses of CMY-2-producing Escherichia coli isolates from humans, animals and food in Germany
BMC Genomics volume 19, Article number: 601 (2018)
Resistance to 3rd-generation cephalosporins in Escherichia coli is mostly mediated by extended-spectrum beta-lactamases (ESBLs) or AmpC beta-lactamases. Besides overexpression of the species-specific chromosomal ampC gene, acquisition of plasmid-encoded ampC genes, e.g. blaCMY-2, has been described worldwide in E. coli from humans and animals. To investigate a possible transmission of blaCMY-2 along the food production chain, we conducted a next-generation sequencing (NGS)-based analysis of 164 CMY-2-producing E. coli isolates from humans, livestock animals and foodstuff from Germany.
The data of the 164 sequenced isolates revealed 59 different sequence types (STs); the most prevalent ones were ST38 (n = 19), ST131 (n = 16) and ST117 (n = 13). Two STs were present in all reservoirs: ST131 (human n = 8; food n = 2; animal n = 6) and ST38 (human n = 3; animal n = 9; food n = 7). All but one CMY-2-producing ST131 isolates belonged to the clade B (fimH22) that differed substantially from the worldwide dominant CTX-M-15-producing clonal lineage ST131-O25b clade C (fimH30). Plasmid replicon types IncI1 (n = 61) and IncK (n = 72) were identified for the majority of blaCMY-2-carrying plasmids. Plasmid sequence comparisons showed a remarkable sequence identity, especially for IncK plasmids. Associations of replicon types and distinct STs were shown for IncK and ST57, ST429 and ST38 as well as for IncI1 and ST58. Additional β-lactamase genes (blaTEM, blaCTX-M, blaOXA, blaSHV) were detected in 50% of the isolates, and twelve E. coli from chicken and retail chicken meat carried the colistin resistance gene mcr-1.
We found isolates of distinct E. coli clonal lineages (ST131 and ST38) in all three reservoirs. However, a direct clonal relationship of isolates from food animals and humans was only noticeable for a few cases. The CMY-2-producing E. coli-ST131 represents a clonal lineage different from the CTX-M-15-producing ST131-O25b cluster. Apart from the ST
-driven spread, plasmid-mediated spread, especially via IncI1 and IncK plasmids, likely plays an important role for emergence and transmission of blaCMY-2 between animals and humans.
The production of extended-spectrum β-lactamases (ESBLs) is the worldwide most important mechanism of resistance to 3rd-generation cephalosporins in Escherichia coli. . AmpC beta-lactamases are also able to hydrolyse 3rd-generation cephalosporins. Induction or overexpression of chromosomally-located, species-specific ampC genes and the acquisition of plasmid-encoded ampC genes (e.g. blaCMY-like, blaACC-like, blaDHA-like) have been described in E. coli . The most common plasmidic ampC gene reported in Enterobacteriaceae including E. coli is blaCMY-2 [3, 4]. It originates from the chromosomal ampC gene of Citrobacter freundii, and had been mobilized onto plasmids of different replicon types (IncK, IncI1, IncA/C and IncFIA-FIB) by the insertion sequence (IS) ISEcp1 that also provides the promotor for high-level expression of blaCMY-2 .
In recent years an increasing number of CMY-2-producing E. coli was noticed in the European livestock production. Especially in poultry, there is a prevalence of more than 30% among 3rd-generation cephalosporin-resistant E. coli, whereas only ~ 1% of the 3rd-generation cephalosporin-resistant E. coli from humans harboured the blaCMY-2 gene [3, 5, 6]. Recent studies reported the finding of blaCMY-2 on similar IncK or IncI1 plasmids in unrelated E. coli isolates from poultry and other livestock animals, meat products, humans and companion animals [7,8,9,10]. This finding points towards a zoonotic potential for the dissemination of this resistance determinant via the food production chain. In recent years, a big variety of clonal lineages of E. coli has been reported worldwide, some of them have proved to be dominant. Of interest are ESBL-producing E. coli of sequence type (ST)131 that are prevalent in humans but not in livestock animals and food [3, 11, 12]. Using sequence based analyses ST131 has been grouped into different clades, which are usually associated with specific fimH alleles: clade A (fimH41 ST131-O16), clade B (fimH22 ST131-O25b) and clade C (fimH30 including ST131-O25b fimH30-R/fimH30-Rx) [13, 14]. Especially strains of ST131-O25b with presence of ESBL-type CTX-M-15 and resistance to fluoroquinolones (fimH30-Rx) have been reported worldwide, and frequently represent a cause of infections, particularly urinary tract infections, in human patients [15, 16]. However, CMY-2 production has been hitherto rarely described for ST131 isolates [7, 17].
To investigate a possible transmission of blaCMY-2 along the food production chain, we conducted next-generation sequencing (NGS)-based analysis of CMY-2-producing E. coli isolates from humans, livestock animals and foodstuff from Germany.
Antibiotic susceptibilities, resistance and virulence genes
All 164 CMY-2-producing E. coli isolates of our study were resistant to ampicillin, cefotaxime, ceftazidime and cefoxitin but remained susceptible to imipenem and meropenem with one exception (isolate no. 10–16 with non-susceptibility to imipenem and meropenem). Additional resistance to ciprofloxacin was detected in 25% (41/164) of the isolates. The proportion of resistance to ciprofloxacin in CMY-2-producing E. coli from humans was 43.5% (20/46) remarkably higher than with 17.8% (21/118) in isolates from livestock animals and meat products.
Presence of resistance genes in the whole genome sequences of the 164 E. coli isolates was investigated by ResFinder (Additional file 1: Table S1). All genomes contained blaCMY-2. Additional β-lacatamase genes were present quite frequently: blaTEM-like (n = 72), blaCTX-M-1/9-group (n = 8), blaOXA-1-like (n = 5), blaSHV-like (n = 5) and blaOXA-10-like (n = 1). Several isolates carried plasmid-mediated quinolone resistance (PMQR) genes that contribute to ciprofloxacin-resistance: aac(6′)-Ib-cr (n = 7), qnrS1 (n = 4), qnrB1 (n = 1), and qnrB19 (n = 1). Furthermore, the plasmid-mediated colistin resistance gene mcr-1 was observed in twelve isolates from livestock animals (chickens n = 3; pigs n = 5; turkey n = 1) and food (n = 3). Subsequently performed susceptibility tests (broth microdilution test according to EUCAST criteria v. 7.1) confirmed resistance to colistin (MIC 4–8 mg/L) of these 12 isolates.
Additionally, the presence of shigatoxin gene stx-2 was observed in one isolate (isolate no 6–16, serotype O141:H49).
Our collection of 164 blaCMY-2-positive E. coli isolates showed a high diversity of sequence types (ST). For isolates from human patients (n = 46), isolates from livestock animals (n = 63), and isolates from food (n = 55) we determined 31, 29, and 20 different STs, respectively (Table 1). The most prevalent STs were ST38 (n = 19), ST131 (n = 14) and ST117 (n = 14). In isolates from human patients the most prevalent ST was ST131 (n = 6; 13%), its proportion among isolates from livestock animals and food was 12.2 and 3.8%, respectively. However, in isolates from livestock animals and food products ST38 was most prevalent (chicken meat n = 9; chickens n = 6; cattle n = 1). Several STs were present in E. coli from all three sources: ST38 (humans n = 3; chicken meat n = 9; chickens n = 6; cattle n = 1), ST131 (humans n = 6; chicken meat n = 2; chickens n = 5; turkey n = 1), ST117 (human n = 1; chicken meat n = 5; turkey meat n = 1; chickens n = 7) and ST10 (humans n = 4; chicken meat n = 1; turkey meat n = 1; pig meat n = 2; chickens n = 1; pigs n = 1). In contrast, other STs, e.g. ST429, were only observed in chicken meat products and chicken isolates (n = 3/3) (Table 1).
The application of a core genome multilocus sequence typing (cgMLST) scheme achieved a much higher resolution (comparison of 2547 alleles vs. seven MLST alleles) regarding the genetic relatedness of the CMY-2-producing E. coli isolates. As expected, isolates forming one ST by the classical multilocus sequence typing (MLST) scheme grouped together in the cgMLST scheme, but due to the higher discrimination capability of the cgMLST, subgroups within the respective cgMLST cluster were observed (Fig. 1). The highest diversity was observed for ST10 isolates; the allele differences between ST10s isolates varied between 326 and 982). In general, isolates from the same source (human, livestock animal or food) grouped more closely together, but also clusters of isolates from animals and meat products could be observed, e.g. ST117 and ST2040 (Fig. 1). Most isolates from human patients showed a noticeably higher allele distance to isolates from livestock animals and food. Exceptions were observed for ST131 and ST1463. ST131 isolates from the three sources grouped into mixed clades with an allele difference of 15–205. In contrast, ST38 isolates from human patients separated from the livestock/food isolates by a minimum of 513 alleles, whereas the difference between livestock and food isolates was 1–190 alleles. The most related isolates were identified in ST429 (livestock/food isolates, allele difference 0–46), ST117 and ST3778 (allele difference 0–79). ST117 and ST3778 differ in three nucleotide substitutions in the gyr allele (gyr4- > gyr14: 174C > T; 288C > T; 396 T > C), and therefore, form one group in the cgMLST scheme. All but one (isolate no. 473–14, human tracheal secretion isolate) isolate of these two STs were from livestock (chickens n = 10; pig n = 1) and food (chicken meat n = 5; turkey meat n = 1).
In an additional analysis we compared our ST429 isolates to all available ST429 E. coli isolates from EnteroBase (http://enterobase.warwick.ac.uk/). Our cgMLST analysis demonstrated that a cluster of E. coli isolates from chicken meat products and chicken from Denmark, Germany and France showed a high genetic relatedness (Additional file 2: Figure S1).
Isolates from selected cgMLST clusters (ST131, ST1196, ST429 and ST117/3778) were analysed by a single nucleotide polymorphisms (SNP)-based mapping approach. For each cgMLST cluster the best fitting reference for the mapping-based SNP analysis was identified using refRank. E. coli strain MDR_56 (GenBank accession no. NZ_CP019903.1) was determined as the best available reference for ST117/ST3778. The analysis of ST117/ST3778 isolates revealed three clusters with differences from 0 to 47 SNPs in each cluster (Additional file 3: Figure S2a). Notably, cluster A (ST3778) isolates no. 510–16 and 484–16 from different farms in the federal state of Brandenburg showed no SNPs (exclusion distance d = 0) and originated from chicken faeces (isolated in 2011) and an environmental sample from a pig farm (isolated 2011) within a distance of ca. 50 km. Another isolate (no. 485–16) from chicken faeces, isolated in 2011 in the federal state Saxony-Anhalt, represented a single allele difference and only four SNPs compared to isolates no. 510–16 and 484–16. Cluster B and C equally presented closely related ST117 isolates (3–20 SNPs in cluster B and 11–33 SNPs in cluster C observed) and originated from chicken (n = 4) and chicken meat (n = 1) and turkey meat (n = 1), isolated in 2012 in the federal states Lower Saxony and Hesse (Additional file 3: Figure S2a).
Within ST131, three closely related isolates from livestock in farms of different federal states (all isolated 2012) were observed: no. 464–16 (chicken, Lower-Saxony), no. 500–16 (turkey, Bavaria) and no. 495–16 (chicken, Bavaria) demonstrating 0 (no. 464–16 to no. 500–16) to 3 (no. 500–16 to no. 495–16) allele differences and 1 to 6 SNPs differences, respectively (Additional file 3: Figure S2b). Interestingly, two isolates from human patients (no. 399–15 and no. 1045–14; isolated in November 2013 and August 2014) differed by only 28 and 27 SNPs, respectively, from a chicken isolate (no. 468–16, isolated in August 2012). These three closely related isolates were collected in Western Germany (radius 200 km) within a period of three years indicating a trans-sectoral spread.
The seven identified isolates assigned to ST1196 represented closely related strains from livestock animals, varying from 0 to 125 allele differences in the cgMLST scheme and grouping into two clusters (Additional file 3: Figure S2c). Cluster A included three isolates from turkeys. These isolates were collected in three stables from two farms in the same federal state in 2014. Differences between the isolates were from 1 to 3 alleles and 0 to 3 SNPs (by mapping against for ST1196 isolates best determined reference NC_020518.1 E. coli str. K-12 substr. MDS42). Cluster B consisted of three isolates from pigs, isolated in the federal states North Rhine-Westphalia and Saxony-Anhalt in 2011, showing 0 allele and 0 to 2 SNPs differences. One additional isolate from chicken (no. 472–16) was sampled in 2011 in the same geographical region as the pig isolates (cluster B) but differed by 99 SNPs or by 88 alleles to cluster A isolates.
The two closely related ST1463 isolates from human and turkey meat, respectively, exhibited high similarity in the cgMLST approach (19 allele differences), and only 21 SNPs were detected after mapping (for ST1463 isolates best determined reference: GenBank accession no. NC_020518.1, data not shown). The isolates were identified in Bavaria and North Rhine-Westphalia in the years 2012 and 2014, respectively.
All E. coli-ST429 isolates were from livestock animals and food and showed only a small number of allele variants (Additional file 3: Figure S2d). Two isolates from chickens were collected in the same federal state within a 16 month period between 2011 and 2012 and showed 0 allele and 0 SNP differences. Isolates from chicken meat - which were recovered from meat products bought 2012 in Lower Saxony and Bavaria - presented only 6 and 8 alleles difference, respectively, to the other isolates of this ST from chicken. Furthermore, isolate no. 8–16 from a diseased chicken (isolated in 2012 in Bavaria) showed a slightly distant similarity (32 to 46 allele differences) to all other E. coli-ST429 isolates.
fimH analysis of ST131-E. coli isolates
Typing of the fimH gene of E. coli can be used to elucidate the population structure within a ST and is used to differentiate the subclades of ST131 . Our analysis showed that all but one CMY-2-producing ST131-E. coli (n = 14) belonged to the clade B, and were associated with fimH22. The exceptional isolate no. 177–10 (human origin) harboured fimH30 (clade C) instead and carried an additional ESBL gene blaCTX-M-15. Isolates belonging to clade A (fimH41) were not identified. Sequence data suggested a plasmid location of blaCMY-2 in isolate 177–10 but not on an IncI1 or IncK2 plasmid as found for most CMY-2 E. coli isolates in our study (see below). A screening of nearly 1000 available E. coli-ST131 sequences from the years 2008 to 2017 (sequences were obtained from the International Nucleotide Sequence Database Collaboration (INSDC) by the SS+ NCBI Bacteria Genome Browser and E. coli assemblies from EnteroBase (http://enterobase.warwick.ac.uk/species/index/ecoli)) revealed only 15 blaCMY-2-positive strains (Additional file 4: Figure S3). The majority (n = 9) of these 15 isolates clustered in the applied cgMLST scheme in clade B (fimH22), along with our study isolates; only one isolate harboured a fimH41 allele (clade A). The remaining five isolates carried fimH30 (clade C) but no additionally blaCTX-M-15 gene and IncK or IncI1 replicon sequences were present. In general, most of the compared ST131 isolates from EnteroBase contained fimH30 with blaCTX-M-15 (n = 377), fimH30 without blaCTX-M-15 (n = 287) or fimH30 with blaCTX-M-27 (n = 87). The fimH variants grouped exclusively according to their variants in the cgMLST schema. A total of 789 fimH30 isolates were identified, in addition to 66 fimH41, 16 fimH27 and 84 fimH22; most of the blaCMY-2-positive isolates cluster in the latter clade H22 (Additional file 4: Figure S3). Isolate no. 177–10 clustered together with fimH30 ESBL-positive strains in clade C.
In general, the most frequent replicon types identified by PlasmidFinder were FIB (n = 133), FII (n = 88), FIA (n = 33) and FIC (n = 55), I1 (n = 81) and B/O/K/Z (n = 78); further replicon types were A/C (n = 7), N (n = 10) and I2 (n = 12). Contig aligned to several completely assembled blaCMY-2-carrying plasmids from GenBank database enabled the identification of blaCMY-2 location.
In 72 E. coli isolates (human n = 14; livestock n = 25; food n = 37) blaCMY-2 has been linked to plasmids of the recently defined incompatibility group IncK2  (Additional file 1: Table S1). In 26 isolates the link between blaCMY-2 and the IncK2 replicon sequence was confirmed on the same de novo assembled contig. In the remaining 46 isolates multiple occurrences of identical insertion sequences (IS) in the sequence data prevented a complete de novo assembly. In this case, a mapping of contigs to previously published blaCMY-2-carrying plasmid sequences (pTMSA1088 and pDV45) was conducted to deduce the respective plasmid contigs. Only one plasmid backbone with minor (SNPs) and major (mobile genetic element (MGE) insertions and variable shufflon region) alterations was identified among all IncK2 plasmid sequences. This backbone was identical to the sequence of the annotated and fully assembled blaCMY-2-carrying plasmids pDV45 (85.9 kb) and pTMSA1088 (79.3 kb). These two plasmids differed by a ca. 6.5 kb sized fragment inserted between the genes yfbA and psiB and presented a highly variable shufflon region . We confirmed these plasmid sizes by S1 nuclease pulsed-field gel electrophoresis (PFGE) for selected isolates of this study. The Inc type was confirmed by PBRT .
All 72 IncK2 plasmids showed the previously described genetic environment of blaCMY-2, consisting of ISEcp1 upstream and blc and sugE downstream of blaCMY-2 (Fig. 2a) [18, 19]. However, in two isolates a sequence alteration was observed: Isolate no. 66–08 showed an insertion element (ISKpn26-like) integrated between ISEcp1 and blaCMY-2. Isolate no. 35–16 had a truncated ISEcp1 element due to an ISKpn26 integration. Transconjugants harbouring both blaCMY-2-carrying plasmids were resistant to cefotaxime, ceftazidime and cefoxitin confirming the functionality of this beta-lactamase gene despite sequence alterations upstream of blaCMY-2.
Our genome data showed high nucleotide sequence identity (99.8–100%) when the shufflon region was excluded for 13 isolates (human n = 2; poultry n = 5; poultry meat n = 6) compared with plasmid pTMSA1088 (Genbank accession no: KR905386.1); plasmid sizes ranged from 79.3 kb to 80.6 kb. Insertion of an IS2 element into different locations of the plasmid sequence was found in two of the 13 isolates. The majority of IncK2 plasmids (53/72 isolates) were identical to plasmid pDV45 (KR905384.1). This highly conserved plasmid sequence was found in isolates from human patients (n = 11), broiler chicken (n = 13), pigs (n = 3), broiler meat (n = 25) and turkey meat (n = 1). A nucleotide sequence identity of 98.6–100% was observed when the shufflon region and unique MGEs were excluded. In eleven of these isolates, the plasmids showed unique insertion sites of various MGEs (Additional file 6: Table S2). One isolate (no. 466–16, ST117, chicken) carried both, a pDV45-like plasmid and an IncI1 plasmid (type IncI1 pST12b, see below) carrying blaCMY-2.
In six isolates (all E. coli-ST429; chickens n = 3 and chicken meat n = 3) the IncK2 plasmid backbone was reconstructed, but showed an additional unique nucleotide sequence compared to plasmid pDV45 (Additional file 5: Figure S4). The plasmid sizes ranged after de novo assembly from 114.5–120 kb (confirmed by S1 nuclease PFGE for isolate no. 8–16) and varied due to inserted MGEs (Additional file 6: Table S2). The additional nucleotide sequence as compared to pDV45 included an undescribed TnAs3-like-element, interrupted by an IS1326 element, a mercury resistance operon, genes encoding plasmid stability proteins (stbAB) and a colicin operon.
The plasmid-wide maximum common-gene-approach by a gene-by-gene comparison was conducted for all IncK2 plasmids similar to pDV45 and pTMSA1088. Additionally available sequence data (from Genbank) of further IncK plasmids harbouring blaCMY-2 were included into this comparison. Results revealed a high level of similarity: The amount of differences varied between 0 and a maximum of 20 genes. Nevertheless the majority of investigated plasmid sequences exhibited variations in 0 up to 2 genes. Additional blaCMY-2-containing IncK2 plasmids sequences from previous studies conducted in Europe showed identical allele pattern to plasmid sequences from this study or differed by a maximum of 4 alleles in pairwise comparisons to plasmids from Germany (Fig. 3).
Further replicon typing by the PlasmidFinder software revealed in 81 isolates the presence of an IncI1 sequence (Additional file 1: Table S1). This replicon sequence and blaCMY-2 were found on one contig in 32 isolates. For another 29 isolates, this linkage was confirmed by alignments to reference sequences of IncI1 plasmids. Thus, 61 isolates (human n = 17; chicken n = 20; pig n = 2; turkey n = 1; chicken meat n = 16; turkey meat n = 4; pork n = 1) carried blaCMY-2 on IncI1 plasmids. The existing pMLST scheme for IncI1 plasmids was used to classify these plasmids into different backbone lineages: pST12 (n = 39), pST2 (n = 12), pST43 (n = 2), pST23, pST27, pST26, pST36, pST55, pST240, pST241, and pST242 (each n = 1). Plasmids of pST12 presented two types of the blaCMY-2 integration into the plasmid backbone: IncI1 pST12a (n = 25) and IncI1 pST12b (n = 12). The plasmid backbone of four isolates, identified as pST12 with blaCMY-2, could not be determined.
IncI1 pST12a plasmids (chicken/broiler n = 12; human n = 5; pig n = 1; turkey n = 1; chicken meat n = 6) revealed a high similarity (98.5–99.8% pairwise identity) to plasmid pCVM29188_101 (Genbank accession no: NC_011077.1; 101 kb) from a Salmonella enterica isolate of poultry origin (USA, 2003) and only differed from that by absence of a 2462 bp ISEc22 in the pilus region. The direct integration site of blaCMY-2 and the flanking genes differed in comparison to the genetic environment of blaCMY-2 found in IncI1 pST12b plasmids similar to p12-4374_96 (Genbank accession no: CP012929.1) (Fig. 2a). Plasmid pCVM29188_101 contained a 3915 bp sized fragment (ISEcp1- blaCMY-2 - blc - sugE) that was integrated into the yagA gene. In plasmid p12-4374_96 the blaCMY-2 fragment was 3123 bp in size and contained ISEcp1- blaCMY-2 and a 65 bp truncated version of blc. The integration site of this blaCMY-2 fragment is located in the reading frame for a hypothetical/uncharacterized protein. In addition, insertions of MGEs were observed in eight isolates of the type IncI1 pST12a (Additional file 6: Table S2).
IncI1 pST12b plasmids (human n = 1; chickens n = 5; chicken meat n = 3; turkey meat n = 3) exhibited a high level of sequence identity (93.7–99.9% pairwise identity) with a previously described 96 kb plasmid (p12-4374_96: Genbank accession no: CP012929.1) from a human clinical Salmonella strain from Canada in 2012 . Unique integration sites for MGEs were observed in IncI1 plasmids of five isolates (Additional file 6: Table S2).
The second most common IncI1-pST was pST2 (humans n = 3; pig n = 1; chicken n = 1; chicken meat n = 6; pork n = 1), and these plasmids showed high similarity (nucleotide sequence identity 96.7–99.9%) to plasmid pC-6 (Genbank accession no: KT186369.1) from a human urine sample (2009, Denmark). In nine of the twelve pST2 plasmids, unique integration sites of MGEs were observed (Additional file 6: Table S2). Interestingly, only in plasmids of E. coli-ST58 isolates (no. RL89, no. RL225, no. RL226 and no. 47–16) an insertion of ISEc12 into the reading frame of a hypothetical protein was found.
The plasmid backbones of IncI1 pST43 (n = 2), pST23, pST27, pST26, pST36, pST55, pST240, pST241, pST242 (each n = 1) were not reconstructed but the direct genetic environment of blaCMY-2 in these plasmids was identical to the ISEcp1-blaCMY-2-blc-sugE structure (Fig. 2a).
IncA/C plasmids and chromosomal integration of bla CMY-2
Replicon sequence for IncA/C group and blaCMY-2 gene were found on the same de novo assembled contig of four isolates (humans n = 2; turkey meat n = 1; pork n = 1). The contigs and aligned reads showed high resemblance (up to 98.9% nucleotide sequence identity at 90.9% coverage) to the 160 kb-IncA/C plasmid pSAN1–1736 from Salmonella enterica from bovine skin (GenBank accession no. CP014658.1). S1-PFGE of a transconjugant of isolate no.761–14 confirmed the presence of a blaCMY-2-carrying IncA/C plasmid with an estimated size of 170 kb. The blaCMY-2 gene was in all four isolates flanked by ISEcp1 upstream and blc and sugE downstream (Fig. 2a).
For ten (6.1%) of the 164 CMY-2-producing E. coli isolates (humans n = 7; turkeys n = 2; pig n = 1) the blaCMY-2 gene was found to be integrated into chromosomal genes. As expected blaCMY-2 transfer by broth mate conjugation was unsuccessful for these isolates. The integration of blaCMY-2 appeared for eight of ten isolates at a unique site and two isolates (no. 90–16, 216–16) at the same site in the chromosome (Fig. 2b). Four of these eight isolates (no. 501–16, 502–16, 509–16 and 739–12) showed a twice as high coverage of the ISEcp1-blaCMY-2 sequence compared to adjacent sequence regions, indicating a second integration site of blaCMY-2. Detailed analysis of the contigs and reads revealed for these four isolates the same additional blaCMY-2-containing genetic environment (Fig. 2c): IS4 was adjacent to the ISEcp1-blaCMY-2 sequence; downstream of blaCMY-2 a truncated blc gene (Δblc) was observed, which was followed by a truncated yggR gene (ΔyggR) and a Shigella flexneri 2a SRL pathogenicity island-like region (95% identity). This conserved structure was identical to the previously reported genetic environment TN44889 of blaCMY-2 (GenBank accession no. FM246884) and described to be chromosomally integrated by Fang et al. [18, 21]. However, a direct chromosomal integration of the additional blaCMY-2-containing genetic environment could not be confirmed in the four isolates due to the above mentioned mobile genetic elements being adjacent to this sequence.
In ten further isolates (humans n = 5; turkey n = 1; pigs n = 2; cattle n = 1; turkey meat n = 1) we found only the putative, chromosomally integrated blaCMY-2-environment described by Fang et al. ; plasmid replicons were not identified. Likewise no transfer of blaCMY-2 by broth mate conjugation appeared. In addition, these isolates showed a nucleotide sequence downstream of TN44889 identical to the draft genome of E. coli Co6114 (GenBank accession no. CP016034.1) (Fig. 2c) .
In the present study, we compared core genome sequences and plasmid sequences of 164 blaCMY-2-harbouring E. coli isolates from meat products (n = 63), livestock animals (n = 55), and human patients (n = 46) to assess national or international transmission of CMY-2-producing strains or blaCMY-2-carrying plasmids. Until today, only a few NGS-based transmission analyses of CMY-2-producing E. coli across sectors have been published despite knowing prevalence rates of more than 20% in poultry and poultry meat [7, 9, 10]. In human patients in Europe, ca. 1% of the 3rd generation cephalosporin-resistant E. coli produce CMY-2 [5, 23, 24]. However, in recent studies from Asia much higher rates and an increasing trend among 3rd-generation cephalosporin-resistant E. coli isolates have been reported [25, 26]. Moreover, studies from Denmark, Norway and the Netherlands suggest that blaCMY-2 transmission along the food chain is probable due to either horizontal transfer events or clonal spread [7, 9, 10, 27]. Interestingly, whole genome analyses of ESBL-producing E. coli-ST410 that were found in various sources (human patients, healthy humans, livestock animals, pets and wastewater samples) indicated a high stability of this clonal lineage and a transfer between different reservoirs [27, 28].
Distribution and proportion of clonal lineages – MLST and cgMLST
In this study MLST analysis was done first to provide an overview about the present CMY-2-E. coli population structure, and to enable the comparison with results of previous studies. The 164 CMY-2-producing isolates could be assigned to 59 different STs, and the isolates from human patients showed a higher ST variety than the isolates from livestock animals or meat products (Table 1). The more frequently observed STs, e.g. ST38 (n = 19), ST131 (n = 14), ST117 (n = 14), and ST69 (n = 8), were also found to be prevalent in studies on ESBL-producing E. coli from retail poultry meat, livestock animals and human patients indicating a successful spread of clonal lineages independent from distinct resistance genes [3, 29]. These STs were previously described in other European countries (Norway, Denmark, Sweden, Switzerland and Netherlands) demonstrating the high variability of blaCMY-2-harbouring E. coli lineages in all sources (human/livestock/food) [7,8,9, 30]. Interestingly, we found CMY-2-producing E. coli-ST131 in equal proportions in human patients (n = 6) and chickens (n = 6) and in some meat products (n = 2). In contrast, ESBL-producing E. coli-ST131 isolates have been found in previous studies only occasionally in samples of animal origin [7, 31]. On the other hand, CMY-2-producing E. coli-ST38 isolates were found exclusively in poultry and poultry meat, whereas this ST has been rarely described in ESBL-producing isolates from humans .
A putative clonal transfer along the food chain cannot be drawn from results of typing methods with a comparably low discrimination such as MLST, even in combination with additional data, such as resistance pattern, plasmid replicon type and beta-lactamase gene presence . Hence, we conducted a genome-wide gene-by-gene comparison. The clustering of isolates in our cgMLST scheme was in good concordance with their respective ST groups, but not always identical (e.g. ST10, Fig. 1). ST10 E. coli isolates were noticeable less genetically related than isolates of other STs. This high diversity within ST10 is known and can be explained by the chosen most common MLST scheme (“Warwick scheme”) and the corresponding housekeeping genes, which are not sufficiently discriminative for ST10 compared to the housekeeping genes from the two other available MLST schemes (“Pasteur and Michigan scheme”) [32,33,34,35]. This highlights the importance of using typing techniques with a higher discrimination than MLST to avoid wrong assumptions about genetic relatedness or transmission events. Nevertheless, we used the ST nomenclature to enable an international comparison and a differentiation of isolates that grouped together in our cgMLST scheme.
Within a ST, isolates from human patients mostly showed a larger distance to isolates from livestock or food, but all seemed to have a common ancestor. Sporadic clonal transmission events between different sources (human, livestock, food) in the past followed by parallel independent micro-evolution within the different sources might explain this differentiation between isolates from human and livestock/food isolates . In fact, only within ST1463 and ST131 a very high similarity of isolates from human and animal/food origin was detected in cgMLST and in subsequently performed SNP-based analysis (Additional file 3: Figure S2b). Since these isolates were either from the same geographical region (federal state) or from the same period of time (6–12 month distance), a clonal transmission seems likely. Recently, a research group from Norway reported nearly identical strains of CMY-2-producing E. coli-ST38 (SNP differences 1–13) from chickens and from human patients indicating a clonal transmission as well .
Furthermore, we found a close relationship of isolates from poultry and poultry meat products, especially for E. coli isolates of ST117, ST3778 and ST429 (Additional file 3: Figure S2a,d). The subsequently performed SNP analysis of isolates of ST117 and ST3778 revealed high similarities with 0 to 40 SNPs differences between the isolates. The identified numbers of discriminating SNPs are equal or marginally higher than SNP differences reported previously from isolates derived from clonal outbreaks [10, 36]. In previous studies, epidemiologically linked isolates exhibited SNPs ranging from ≤6, (isolates of a big German EHEC outbreak) ≤ 4 (epidemiologically linked cases of E. coli O157 from Scotland) or ≤ 23 (E. coli O157 pork-associated outbreak in Alberta, Canada) [36,37,38,39].
In a study by Mellmann et al. a threshold of < 11 alleles difference in a pairwise comparison of isolates was used to identify nosocomial transmissions of E. coli in an ad-hoc cgMLST analysis (based on 2325 alleles) . Applying this threshold to the analyzed ST429 isolates with blaCMY-2, the relationship of nearly all isolates would indicate clonal transmission. Taken into consideration that plasmid analyses of these isolates also showed highly similar plasmid sequences to the p486–16 plasmid sequence, a clonal transmission of an E. coli ST429 strain carrying a blaCMY-2-carrying IncK2 plasmid seems likely. This supports the hypothesis of clonal transmission of blaCMY-2 along the poultry production chain and across sectors. The study design did not enable a direct observation of individual animals from breeding to slaughter, whereas this argument remains speculative to a certain extent. However, veterinary studies in other countries support this hypothesis. In the Danish conventional broiler production, the spread of ESBL-producing E. coli clones and ESBL gene-carrying plasmids were observed in imported broiler parent flocks, despite cephalosporins have been never approved for use in poultry . Furthermore, a relatively stable colonization with CMY-2-producing E. coli was observed in a Norwegian study through the whole broiler production chain from grandparent animals to retail meat . Cross-contaminations during the successive fattening of several herds in the same stable due to insufficient cleaning and cross-contaminations in big slaughterhouses, where animals from many farms were processed, facilitate the spread of distinct clones .
Resistance and virulence
Studies on ESBL-producing E. coli of human origin reported rates of 60–70% fluoroquinolone resistance [11, 43]. Our data showed that 43.5% of the blaCMY-2-harbouring isolates from human origin were resistant to ciprofloxacin. Differences in the observed fluoroquinolone resistance rates could be explained by the occurrence of different proportions of E. coli-ST131 in ESBL-E. coli compared to CMY-2-E. coli . The prevalence of E. coli-ST131 is observed up to 70% among ESBL-positive E. coli associated with nosocomial or ambulant urinary tract infections . Further, CMY-2-producing E. coli-ST131 isolates belong mainly to the fluoroquinolone-susceptible ST131 clade B (fimH22; see also next chapter) and not to clade C (fimH30); the latter comprises mostly strains resistant to fluoroquinolones . However, rates of fluoroquinolone-resistant CMY-2-E. coli from livestock animals or meat products were comparable to the rates of resistance to nalidixic acid and ciprofloxacin among E. coli from broilers in Sweden, but noticeable lower than in human E. coli isolates .
The recently described colistin resistance gene mcr-1 was observed additionally in twelve CMY-2-positive E. coli strains from this study in chicken (n = 3), pig (n = 5), turkey (n = 1) and chicken meat products (n = 3). Mcr-1 in Germany is found predominantly in the poultry production chain . The higher number of findings of mcr-1 in isolates from pigs compared to isolates from chicken in the present study is most likely biased by the isolate selection process for this study.
In addition, in one isolate from a diseased pig (no. 6–16, ST10, serotype O141:H49) analysis exhibited the presence of stx-2. Shiga toxin-production has been rarely described in ESBL/AmpC-producing E. coli so far [46, 47].
CMY-2-producing E. coli-ST131
We identified a substantial amount of E. coli-ST131 isolates (human patients 6/46, 13.0%; poultry 6/45, 13.3%; food 2/63, 3.2%), collected from 2010 till 2016. Regarding ESBL-producing E. coli, ST131 is the most frequent ST in humans but occurs only sporadically in pets, livestock or wild animals (0–1.8% of 3rd generation cephalosporin-resistant E. coli) [3, 11, 12]. All but one CMY-2-producing E. coli-ST131 isolates from the present study harboured fimH22 and clustered in clade B. The single non-fimH22 isolate (no. 177–10) belonged to the successful ST131 subclone H30-Rx and carried the putative chromosomally encoded ESBL gene blaCTX-M-15. The comparison of 15 published whole genome sequences of E. coli-ST131 isolated between 2008 and 2017 with blaCMY-2, using the ad-hoc cgMLST approach, showed that nine isolates from other countries likewise contained fimH22 and differed in 3 to 113 target genes from CMY-2-producing E. coli-ST131 isolates from this study (Additional file 4: Figure S3). The affiliation of our CMY-2-producing E. coli-ST131 into clade B (fimH22) in contrast to the ESBL-producing E. coli-ST131 of clade C (ST131-O25b fimH30-R/fimH30-Rx) suggests a CMY-2-producing E. coli-ST131 population with an independent evolution in different reservoirs. ST131-O25b with ESBL-production often occurs in human patients and isolates might be adapted to colonization and infection in the human bladder . In contrast, CMY-2-producing E. coli-ST131 isolates were described rarely in humans [31, 49, 50] and recent studies from Germany, the Netherlands and Denmark reported only a few isolates in poultry [7, 17]. The much higher proportion of ST131 in our CMY-2-producing E. coli collection remains unclear; we did not recognize any obvious bias in the sampling strategy of the underlying studies for the isolate collection which might have influenced enrichment for ST131 clade B (fimH22) isolates in our sample.
The integration of the blaCMY-2 gene occurred in a narrower spectrum of genetic structures and plasmid backbones in relation to the comparably high variety of clonal lineages. blaCMY-2 has been frequently reported on IncI1, IncK and IncA/C plasmids in E. coli and sporadically on IncF, IncX, IncI2 or IncHI2 plasmids [21, 51,52,53]. However, in our collection we only identified IncI1, IncK and IncA/C plasmids carrying blaCMY-2. The most prevalent plasmid types were of incompatibility groups IncK2 (n = 76) and IncI1 (n = 61). Previous studies reported that dissemination of blaCMY-2 in E. coli is mainly driven by plasmids of these two Inc groups within Europe [7,8,9, 54]. In contrast, IncA/C plasmids seem to play a more prominent role in the spread of blaCMY-2 in North America [18, 55,56,57]. In a Dutch study all blaCMY-2-harbouring plasmids belonged to IncA/C indicating the importance of clonal spread in isolated populations such as broiler stocks .
Our analysis revealed that all IncK plasmids belonged to the recently described IncK2 type . Sequence identity among these plasmids was high and the most frequent plasmid backbone structures were similar (98.6–100% identity) to two previously published plasmids pTMSA1088 (n = 13) (KR905386.1) and pDV45 (n = 54) (KR905384.1) . Plasmids similar to pTMSA1088 were found in isolates of ST38 (n = 11), ST10 (n = 1) and ST117 (n = 1). This specific ST38/pTMSA1088-like combination has been reported before in broiler meat production of several countries (Netherlands, Denmark, Norway, Sweden) and raised the discussion if a highly conserved IncK2 plasmid carrying blaCMY-2 is genetically linked to ST38 [7,8,9]. In a direct gene-by-gene comparison of all pTMSA1088-like and pDV45-like plasmids from the present study with blaCMY-2 plasmids from the Netherlands, Denmark, Norway and Switzerland, we confirmed the expected high relatedness. The low number of variations observed in the genes of these plasmids indicates a horizontal transmission, and for the ST38/pTMSA1088-like combination a dissemination of the plasmids linked to a clonal transmission across the European broiler production is highly probable (Fig. 3) [7, 8, 10, 59]. Moreover, a modification of the plasmid backbone structure in comparison to plasmid pDV45 was observed in six ST429 isolates (Additional file 5: Figure S4). The occurrence of this backbone structure only in ST429 suggests a clonal transmission linked to the plasmid. Identical IS insertions in the plasmids and SNP-based comparisons of the core genomes (Additional file 3: Figure S2d) support this suggestion. The modified plasmid backbone structure contained additional genes encoding a plasmid segregation system (stbAB) which enables a stable vertical transmission during cell division as well as horizontal transmission (conjugation), and thus it promotes efficient blaCMY-2 plasmid propagation .
In contrast to the IncK2 plasmids, we detected more diverse IncI1 plasmid structures. IncI1 plasmids carrying blaCMY-2 were reported worldwide and differ considerably in their plasmid backbone structures [7, 18, 61, 62]. In the present study most IncI1 plasmids could be assigned to pST12 (39/61) and pST2 (12/61). The majorities of IncI1 plasmids of each respective pSTs were highly similar in their conserved backbone structures and differed only in a few SNPs or an MGE insertion. Among IncI1-pST12 plasmids two different types (IncI1 pST12a and IncI1 pST12b) were identified differing in the direct genetic surrounding of blaCMY-2 and an additional plasmid addiction system. These differences in the integration site of blaCMY-2 indicate two independent mobilization events into highly similar plasmid backbones.
In general, the distribution of blaCMY-2-carrying IncK2 and IncI1 plasmids was heterogeneous among isolates from human patients, livestock animals and food as well as in various E. coli-STs (Fig. 1) which supports the hypothesis that the horizontal transfer via plasmids plays the major role for transmission of blaCMY-2 between the reservoirs. Given the capability of these plasmids to transfer themselves within E. coli and into other Enterobacteriaceae genera, the ingested blaCMY-2-carrying strains from meat products might transfer their resistance gene to the human host-adapted strains .
bla CMY-2 in the chromosome
In 20 isolates, mainly from human patients (12/20), the localization of blaCMY-2 on the chromosome was either confirmed or suspected to be chromosomal. While for ten isolates blaCMY-2 was found to be located on the chromosome by sequence and read analysis, for further ten isolates the localization of blaCMY-2 on the chromosome was suspected due to the identity to a previously described blaCMY-2 structure: A recent publication from China reports the integration of this blaCMY-2-containing fragment TN44889 into the chromosome of E. coli strains (9/469) isolated from companion and livestock animals . It is known that blaCMY-2 can integrate into the chromosome of Proteus mirabilis or S. enterica but only in few studies isolates with blaCMY-2 integrated into the chromosome of E. coli were identified [64, 65]. The variety of confirmed and presumable chromosomal integration sites and the differences in the blaCMY-2genetic environments in general point towards independent events of chromosomal integration of this resistance gene. However, the high similarity (2 to 3 SNP-differences) of three ST1196 isolates from livestock animals from our collection (Additional file 3: Figure S2c) with an identical insertion site of ISApl1 suggested propagation by clonal transmission in this specific setting. The reason for the comparable high number of isolates from humans (n = 12) with chromosomal blaCMY-2 integration sites in this study remains unknown.
The study has several limitations to be considered. First, the sampling of the 164 CMY-2-producing E. coli isolates was conducted over a period of 5 years in the scope of different studies (human patients, livestock animals, food) of the national research project RESET. Regarding studies on livestock animals, we selected only one isolate from each farm or isolates from the same farm with different properties, e.g. different phylogenetic groups, to achieve a preferably heterogeneous sample of the distribution of blaCMY-2-carrying isolates in German farms. Therefore, the selection of the isolates, especially from broiler production, is not completely randomized but as diverse as possible. Second, the long sampling period of five years was necessary due to the comparably low prevalence of CMY-2-producing E. coli among 3rd-generation cephalosporin-resistant E. coli from human infections and colonizations. Third, putative links between isolates from animals, food and humans could not be further elucidated due to limited clinical data of patients and missing information on origin and transport ways of livestock animals and meat products.
The whole genome sequence analysis of 164 CMY-2-producing E. coli isolated from human patients, livestock animals and meat products in Germany revealed a high diversity of STs across all sources; most frequent types were ST38, ST131 and ST117. The CMY-2-producing ST131 isolates from human patients and livestock animals belonged to another sub-lineage (clade B; fimH22) than the worldwide prevalent multidrug-resistant and ESBL-producing ST131-O25b lineage (clade C; fimH30) known from humans. Our data suggest that clonal transmission of blaCMY-2 is a rare event while the horizontal transfer of temporally stable blaCMY-2-carrying IncK2 and IncI1 plasmids is more likely the dominant way of transmission between humans and animals. At present, CMY-2-carrying E. coli are rarely detected from human infections in Germany and Central Europe. However, the large reservoir of blaCMY-2 -carrying plasmids in livestock animals, especially poultry, poses a serious future risk of a more pronounced potential of CMY-2-producing E. coli isolates causing infections in humans and animals.
We included 164 CMY-2-producing E. coli from Germany collected 2008–2016. The majority of isolates (n = 149) was collected and isolated between 2011 and 2014 in the scope of the national research consortium “RESET” (www.reset-verbund.de) by different project partners: Robert Koch Institute (RKI), German Federal Institute for Risk Assessment (BfR), Friedrich-Loeffler-Institut (FLI) and Freie Universität Berlin (FU). No concrete sampling plan for the collection of the isolates was developed; rather all isolates were collected in the scope of different studies including isolates from healthy livestock animals, diseased livestock animals, food samples and human patients, and represent therefore a random collection of CMY-2-producing E. coli isolates from Germany (Additional file 1: Table S1). To enable a higher comparability of studies on different animal species, the RESET consortium harmonized protocols for sampling and basic phenotypical and genotypical analyses were developed in advance. FU provided 50 CMY-2-producing E. coli isolates from chicken (n = 9 in 2011 and n = 28 in 2012), turkey (n = 1 in 2012 and n = 3 in 2014), pig (n = 6 in 2011 and n = 2 in 2012) and cattle (n = 1 in 2012) for WGS analysis. These isolates were collected in longitudinal- and cross-sectional studies in different livestock farms during the years 2011–2014 [66,67,68,69,70]; only isolates from unique production sites were included. In addition, BfR provided 63 CMY-2-producing E. coli isolates that were identified from 2500 food samples of different origins and matrices, collected and investigated by different German state laboratories, in the scope of a cross-sectional study on ESBL-E. coli in foodstuff in Germany (unpublished data). In brief, the 63 isolates were from chicken meat samples (n = 53), turkey meat samples (n = 7) and pig meat samples (n = 3), and were isolated between 2011 and 2013. Furthermore, FLI identified five CMY-2-producing E. coli isolates from diseased livestock animal (chicken n = 3, pig n = 1 and cattle n = 1) that were included in the present study . These five isolates were detected between 2008 and 2015 within the GERM-Vet program, a monitoring program collecting data of resistance of pathogenic bacteria and including only isolates from clinically diseased, non antibiotically pretreated animals. Finally, CMY-2-producing E. coli from human clinical samples (n = 46) were collected by RKI between 2008 and 2016. These isolates were obtained from different studies [5, 11] or were sent from German laboratories to RKI for confirmatory diagnostics (Additional file 1: Table S1).
Antimicrobial susceptibility testing and resistance gene screening
Susceptibilities to ampicillin, cefotaxime, ceftazidime, cefoxitin, ciprofloxacin, imipenem and meropenem were tested for all E. coli isolates by disk diffusion (Oxoid Ltd., Basingstoke, United Kingdom) according to the manufacturer’s instructions, and interpretation was done according to EUCAST criteria (http://www.eucast.org/clinical_breakpoints; v7.1). The phylogenetic group of the E. coli isolates was determined by a PCR-based assay . The presence of β-lactamase genes (blaTEM, blaSHV, blaCTX-M-group-1-2-9, and blaCMY) was confirmed by PCR and sequencing as previously described [73, 74]; only CMY-2 positive isolates were included in the present study.
bla CMY-2 transferability and plasmid typing
A broth mate conjugation assay was performed for 60 randomly selected isolates using the sodium azide-resistant E. coli K12 J53 recipient strain . Transconjugants were selected on Lysogeny-Broth (LB) agar supplemented with 200 mg/L sodium azide and 100 mg/L ampicillin. Selected transconjugants were tested for the presence of blaCMY-2 by PCR as previously described, and the plasmid sizes were determined using S1 nuclease-restriction followed by pulsed-field gel electrophoresis (PFGE) [73, 76]. Plasmidic DNA was extracted by the Qiagen Plasmid mini Kit (Hilden, Germany) and followed by PCR-based replicon typing (PBRT-Kit, Diatheva, Italy).
Whole genome sequencing and sequence reconstruction
All 164 blaCMY-2-containing E. coli isolates were sequenced using an Illumina MiSeq platform (Illumina, San Diego, U.S.A). Isolate preparation, DNA extraction and sequencing were carried out by RKI, BfR and LGC Genomics GmbH (Berlin, Germany). At the RKI genomic DNA of 58 E. coli isolates was extracted using the DNA Blood & Tissue Kit according to the manufacturer’s instruction (Qiagen, Hilden, Germany), followed by library preparation, using Nextera XT library (Illumina, USA), and sequenced on an Illumina MiSeq platform using the MiSeq v3 reagent kit (Illumina, USA) with 2 × 300 bp paired-end reads. Furthermore, 62 E. coli isolates were sent to the company LGC Genomics that performed DNA extraction and sequencing, using the Illumina MiSeq platform MiSeq v3 reagent kit (Illumina, USA) with 2 × 300 bp paired-end reads. At the BfR, extraction of genomic DNA of 44 isolates was carried out using PureLink® genomic DNA Kit from Invitrogen™ (by ThermoFisher Scientific).
The following RKI internal procedure was used: The resulting raw reads were processed by Trimmomatic (v. 0.0.9; default parameters except maxinfo 15:0.5) and assembled by A5-miseq (v. 0.0.9 beta; default parameters) [77, 78]. The quality of the read data was assessed by using the PHRED scores. Raw sequence data had an average PHRED score above 30. The average read quality was increased by trimming to a PHRED score of above 35.
All contigs were submitted to the CGE Finder Series (Centre for Genomic Epidemiology, Technical University of Denmark (DTU), https://cge.cbs.dtu.dk/services/). Different analysis tools (MLST 1.8, pMLST 1.4, PlasmidFinder 1.3, ResFinder 2.1, VirulenceFinder 1.5) were applied to extract the multilocus sequence type (ST) according to Wirth et al., plasmid multilocus sequence type (pST), plasmid replicon types and information on genes mediating resistance to β-lactams and fluoroquinolones and distinct virulence genes (stx1/2), respectively [32, 79,80,81,82,83].
For deeper phylogenetic analysis (i) a gene-based comparison approach and (ii) a single nucleotide polymorphisms (SNP)-based mapping analysis were performed. For this purpose, an ad-hoc core genome multilocus sequence typing scheme (cgMLST) was created using SeqSphere+ (v. 4.0.1, Ridom GmbH, Münster, Germany) as described before [84, 85]. E. coli O157:H7 str. Sakai (GenBank accession no. NC_002695.1) was used as reference genome and further 43 E. coli genomes from the National Centre for Biotechnology Information (NCBI) were selected as query genomes. Open reading frames (ORFs) were predicted and extracted by using the cgMLST Target Definer v1.4 of SeqSphere+, resulting in 2547 shared genes among these isolates, which were defined as core genome (Additional file 7: Table S3). Obtained cgMLST target gene variants were used to visualize the phylogenetic distance by calculating a Neighbour-Joining tree (parameters were: pairwise ignoring missing values; % columns difference) based on the distance matrix of the core genome differences. The tree was visualized using iTOL (v. 3.5.4) (http://itol.embl.de/) .
Isolates which presented a high relationship in the cgMLST scheme were further investigated by single nucleotide polymorphisms (SNP)-based mapping analysis. A suitable reference was identified using refRank as described previously (https://gitlab.com/s.fuchs/refRank) . Paired-end reads were mapped to the identified best reference genome for respective selected isolates (ST131: NZ_CP019008.1; ST117/ST3778: NZ_CP019903.1; ST1196: NC_020518.1; ST429: NC_013654.1) by BWA-SW (v. 0.7.15-r1140; default parameters) . Variant calling was performed using VarScan (parameter: min-coverage 10; min-reads2 6; min-avg-qual 20; min-var-freq 0.8; min-var-for-hom 0.75; p-value 0.01; strand-filter 0) . SNPs were filtered using SNPfilter (https://gitlab.com/s.fuchs/snpfilter) . Maximum likelihood trees were calculated using RAxML with a GTR GAMMA nucleotide model (rapid hill-climbing, using 100 starting trees) . Phylogenetic trees were visualized by iTOL (v. 3.5.4).
Regarding the identification of the blaCMY-2-harbouring contigs, all contigs were aligned using Geneious (v. 10.0.5, Biomatters Ltd., Auckland, New Zealand) to the blaCMY-2 gene (GenBank accession no. X91840.1). Contigs containing blaCMY-2 were hereafter aligned to several fully sequenced blaCMY-2-carrying plasmids of different Incompatibility groups (Inc) from GenBank using Geneious to identify highest resemblance to a reference plasmid: IncK2: pDV45 (GenBank accession no. KR905384.1), pTMSA1088 (GenBank accession no. KR905386.1); IncI1: p12-4374_96 (GenBank accession no. CP012929.1), pC-6 (GenBank accession no. KT186369.1), pCVM29188_101 (GenBank accession no. NC_011077.1); IncA/C: pSAN1–1736 (GenBank accession no. CP014658.1).
All contigs of a respective isolate were aligned to the identified best fitting reference plasmid. Additionally, read alignment of trimmed reads (Trimmomatic v. 0.0.9; default parameters except maxinfo 15:0.5) to the selected reference plasmid using the Geneious (v. 10.0.5) mapper (medium-low sensitivity, no iterations) was performed. The results obtained were then compared by aligning with the respective reference plasmids and examined for MGEs, using the ISfinder (https://www-is.biotoul.fr/index.php) .
To compare IncK plasmids from this study with available IncK plasmids carrying blaCMY-2 from nucleotide sequence databases, a maximum common gene approach was conducted. For that purpose ORFs from the blaCMY-2-carrying IncK2 plasmid pDV45 were predicted and extracted by using the cgMLST Target Definer v1.4 of SeqSphere+, resulting in 89 genes, of which 67 were found in all the plasmids to be compared. Target gene distance was visualized by a minimum spanning tree (based on 67 targets, pairwise ignoring missing values).
Nucleotide sequence accession numbers
Sequence data were submitted to the European Nucleotide Archive (http://www.ebi.ac.uk/ena) under the study accession number PRJEB23663.
German Federal Institute for Risk Assessment
Centre for Genomic Epidemiology
core-genome multilocus sequence typing
- E. coli :
European Committee on Antimicrobial Susceptibility Testing
Freie Universität Berlin
International Nucleotide Sequence Database Collaboration
Mobile genetic element
Minimum inhibitory concentration
Multilocus sequence typing
National Centre for Biotechnology Information
Open reading frames
Plasmid- Multilocus sequence typing
PCR-based replicon typing
Plasmid sequence type
Robert Koch Institute
single nucleotide polymorphism
Whole genome sequencing
Pfeifer Y, Cullik A, Witte W. Resistance to cephalosporins and carbapenems in gram-negative bacterial pathogens. Int J Med Microbiol. 2010;300(6):371–9.
Jacoby GA. AmpC beta-lactamases. Clin Microbiol Rev. 2009;22(1):161–82.
Ewers C, Bethe A, Semmler T, Guenther S, Wieler LH. 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. 2012;18(7):646–55.
Denisuik AJ, Lagace-Wiens PR, Pitout JD, Mulvey MR, Simner PJ, Tailor F, et al. Molecular epidemiology of extended-spectrum beta-lactamase-, AmpC beta-lactamase- and carbapenemase-producing Escherichia coli and Klebsiella pneumoniae isolated from Canadian hospitals over a 5 year period: CANWARD 2007-11. J Antimicrob Chemother. 2013;68(Suppl 1):i57–65.
Valenza G, Nickel S, Pfeifer Y, Eller C, Krupa E, Lehner-Reindl V, et al. Extended-spectrum-beta-lactamase-producing Escherichia coli as intestinal colonizers in the German community. Antimicrob Agents Chemother. 2014;58(2):1228–30.
Jorgensen RL, Nielsen JB, Friis-Moller A, Fjeldsoe-Nielsen H, Schonning K. Prevalence and molecular characterization of clinical isolates of Escherichia coli expressing an AmpC phenotype. J Antimicrob Chemother. 2010;65(3):460–4.
Hansen KH, Bortolaia V, Nielsen CA, Nielsen JB, Schonning K, Agerso Y, et al. Host-specific patterns of genetic diversity among IncI1-Igamma and IncK plasmids encoding CMY-2 beta-lactamase in Escherichia coli isolates from humans, poultry meat, poultry, and dogs in Denmark. Appl Environ Microbiol. 2016;82(15):4705–14.
Seiffert SN, Carattoli A, Schwendener S, Collaud A, Endimiani A, Perreten V. Plasmids carrying bla CMY-2/4 in Escherichia coli from poultry, poultry meat, and humans belong to a novel IncK subgroup designated IncK2. Front Microbiol. 2017;8:407.
Berg ES, Wester AL, Ahrenfeldt J, Mo SS, Slettemeas JS, Steinbakk M, et al. Norwegian patients and retail chicken meat share cephalosporin-resistant Escherichia coli and IncK/bla CMY-2 resistance plasmids. Clin Microbiol Infect. 2017;23(6):407.e9-407.e15. https://doi.org/10.1016/j.cmi.2016.12.035. [Epub 2017 Jan 7]
de Been M, Lanza VF, de Toro M, Scharringa J, Dohmen W, Du Y, et al. Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages. PLoS Genet. 2014;10(12):e1004776.
Pietsch M, Eller C, Wendt C, Holfelder M, Falgenhauer L, Fruth A, et al. Molecular characterisation of extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli isolates from hospital and ambulatory patients in Germany. Vet Microbiol. 2017;200:130–7.
Nicolas-Chanoine MH, Bertrand X, Madec JY. Escherichia coli ST131, an intriguing clonal group. Clin Microbiol Rev. 2014;27(3):543–74.
Pitout JD, DeVinney R. Escherichia coli ST131: a multidrug-resistant clone primed for global domination. F1000Res. 2017;6(F1000 Faculty Rev):195. https://doi.org/10.12688/f1000research.10609.1.
Price LB, Johnson JR, Aziz M, Clabots C, Johnston B, Tchesnokova V, et al. The epidemic of extended-spectrum-beta-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx. MBio. 2013;4(6):e00377–13.
Rogers BA, Sidjabat HE, Paterson DL. Escherichia coli O25b-ST131: a pandemic, multiresistant, community-associated strain. J Antimicrob Chemother. 2011;66(1):1–14.
Livermore DM, Canton R, Gniadkowski M, Nordmann P, Rossolini GM, Arlet G, et al. CTX-M: changing the face of ESBLs in Europe. J Antimicrob Chemother. 2007;59(2):165–74.
Day MJ, Rodriguez I, van Essen-Zandbergen A, Dierikx C, Kadlec K, Schink AK, et al. Diversity of STs, plasmids and ESBL genes among Escherichia coli from humans, animals and food in Germany, the Netherlands and the UK. J Antimicrob Chemother. 2016;71(5):1178–82.
Verdet C, Gautier V, Chachaty E, Ronco E, Hidri N, Decre D, et al. Genetic context of plasmid-carried bla CMY-2-like genes in Enterobacteriaceae. Antimicrob Agents Chemother. 2009;53(9):4002–6.
Kang MS, Besser TE, Call DR. Variability in the region downstream of the bla CMY-2 beta-lactamase gene in Escherichia coli and Salmonella enterica plasmids. Antimicrob Agents Chemother. 2006;50(4):1590–3.
Labbe G, Edirmanasinghe R, Ziebell K, Nash JH, Bekal S, Parmley EJ, et al. Complete genome and plasmid sequences of three canadian isolates of Salmonella enterica subsp. enterica Serovar Heidelberg from human and food sources. Genome Announc. 2016;4(1)
Fang LX, Sun J, Li L, Deng H, Huang T, Yang QE, et al. Dissemination of the chromosomally encoded CMY-2 cephalosporinase gene in Escherichia coli isolated from animals. Int J Antimicrob Agents. 2015;46(2):209–13.
Both A, Huang J, Kaase M, Hezel J, Wertheimer D, Fenner I, et al. First report of Escherichia coli co-producing NDM-1 and OXA-232. Diagn Microbiol Infect Dis. 2016;86(4):437–8.
Mo SS, Norstrom M, Slettemeas JS, Lovland A, Urdahl AM, Sunde M. Emergence of AmpC-producing Escherichia coli in the broiler production chain in a country with a low antimicrobial usage profile. Vet Microbiol. 2014;171(3–4):315–20.
Hamprecht A, Rohde AM, Behnke M, Feihl S, Gastmeier P, Gebhardt F, et al. Colonization with third-generation cephalosporin-resistant Enterobacteriaceae on hospital admission: prevalence and risk factors. J Antimicrob Chemother. 2016;71(10):2957–63.
Jean SS, Hsueh PR, Group SA-P. Distribution of ESBLs, AmpC beta-lactamases and carbapenemases among Enterobacteriaceae isolates causing intra-abdominal and urinary tract infections in the Asia-Pacific region during 2008-14: results from the study for monitoring antimicrobial resistance trends (SMART). J Antimicrob Chemother. 2017;72(1):166–71.
Harris PNA, Ben Zakour NL, Roberts LW, Wailan AM, Zowawi HM, Tambyah PA, et al. Whole genome analysis of cephalosporin-resistant Escherichia coli from bloodstream infections in Australia, New Zealand and Singapore: high prevalence of CMY-2 producers and ST131 carrying blaCTX-M-15 and blaCTX-M-27. J Antimicrob Chemother. 2017;
Falgenhauer L, Imirzalioglu C, Ghosh H, Gwozdzinski K, Schmiedel J, Gentil K, et al. Circulation of clonal populations of fluoroquinolone-resistant CTX-M-15-producing Escherichia coli ST410 in humans and animals in Germany. Int J Antimicrob Agents. 2016;47(6):457–65.
Schaufler K, Semmler T, Wieler LH, Wohrmann M, Baddam R, Ahmed N, et al. Clonal spread and interspecies transmission of clinically relevant ESBL-producing Escherichia coli of ST410--another successful pandemic clone? FEMS Microbiol Ecol. 2016;92(1)
Manges AR. Escherichia coli and urinary tract infections: the role of poultry-meat. Clin Microbiol Infect. 2016;22(2):122–9.
Arvand M, Moser V, Pfeifer Y. Prevalence of extended-spectrum-beta-lactamase-producing Escherichia coli and spread of the epidemic clonal lineage ST131 in nursing homes in Hesse. Germany J Antimicrob Chemother. 2013;68(11):2686–8.
Naseer U, Haldorsen B, Simonsen GS, Sundsfjord A. Sporadic occurrence of CMY-2-producing multidrug-resistant Escherichia coli of ST-complexes 38 and 448, and ST131 in Norway. Clin Microbiol Infect. 2010;16(2):171–8.
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.
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:560.
Reid SD, Herbelin CJ, Bumbaugh AC, Selander RK, Whittam TS. Parallel evolution of virulence in pathogenic Escherichia coli. Nature. 2000;406(6791):64–7.
Clermont O, Gordon D, Denamur E. Guide to the various phylogenetic classification schemes for Escherichia coli and the correspondence among schemes. Microbiology. 2015;161(Pt 5):980–8.
Parsons BD, Zelyas N, Berenger BM, Chui L. Detection, characterization, and typing of Shiga toxin-producing Escherichia coli. Front Microbiol. 2016;7:478.
Grad YH, Lipsitch M, Feldgarden M, Arachchi HM, Cerqueira GC, Fitzgerald M, et al. Genomic epidemiology of the Escherichia coli O104:H4 outbreaks in Europe, 2011. Proc Natl Acad Sci U S A. 2012;109(8):3065–70.
Holmes A, Allison L, Ward M, Dallman TJ, Clark R, Fawkes A, et al. Utility of whole-genome sequencing of Escherichia coli O157 for outbreak detection and epidemiological surveillance. J Clin Microbiol. 2015;53(11):3565–73.
Sabat AJ, Budimir A, Nashev D, Sa-Leao R, van Dijl J, Laurent F, et al. Overview of molecular typing methods for outbreak detection and epidemiological surveillance. Euro Surveill. 2013;18(4):20380.
Mellmann A, Bletz S, Boking T, Kipp F, Becker K, Schultes A, et al. Real-time genome sequencing of resistant Bacteria provides precision infection control in an institutional setting. J Clin Microbiol. 2016;54(12):2874–81.
Agerso Y, Jensen JD, Hasman H, Pedersen K. Spread of extended spectrum cephalosporinase-producing Escherichia coli clones and plasmids from parent animals to broilers and to broiler meat in a production without use of cephalosporins. Foodborne Pathog Dis. 2014;11(9):740–6.
Projahn M, Daehre K, Roesler U, Friese A. Extended-spectrum-beta-lactamase- and plasmid-encoded cephamycinase-producing Enterobacteria in the broiler hatchery as a potential mode of pseudo-vertical transmission. Appl Environ Microbiol. 2017;83(1)
Cassier P, Lallechere S, Aho S, Astruc K, Neuwirth C, Piroth L, et al. Cephalosporin and fluoroquinolone combinations are highly associated with CTX-M beta-lactamase-producing Escherichia coli: a case-control study in a French teaching hospital. Clin Microbiol Infect. 2011;17(11):1746–51.
Borjesson S, Guillard T, Landen A, Bengtsson B, Nilsson O. Introduction of quinolone resistant Escherichia coli to Swedish broiler population by imported breeding animals. Vet Microbiol. 2016;194:74–8.
Irrgang A, Roschanski N, Tenhagen BA, Grobbel M, Skladnikiewicz-Ziemer T, Thomas K, et al. Prevalence of mcr-1 in E. coli from livestock and food in Germany, 2010–2015. PLoS One. 2016;11(7):e0159863.
Arvand M, Bettge-Weller G, Fruth A, Uphoff H, Pfeifer Y. Extended-spectrum beta-lactamase-producing Shiga toxin gene (stx1)-positive Escherichia coli O91:H14 carrying bla CTX-M-15 on an IncI1-ST31 plasmid isolated from a human patient in Germany. Int J Med Microbiol. 2015;305(3):404–7.
Donaldson SC, Straley BA, Hegde NV, Sawant AA, DebRoy C, Jayarao BM. Molecular epidemiology of ceftiofur-resistant Escherichia coli isolates from dairy calves. Appl Environ Microbiol. 2006;72(6):3940–8.
Vimont S, Boyd A, Bleibtreu A, Bens M, Goujon JM, Garry L, et al. The CTX-M-15-producing Escherichia coli clone O25b: H4-ST131 has high intestine colonization and urinary tract infection abilities. PLoS One. 2012;7(9):e46547.
Oteo J, Cercenado E, Cuevas O, Bautista V, Delgado-Iribarren A, Orden B, et al. AmpC beta-lactamases in Escherichia coli: emergence of CMY-2-producing virulent phylogroup D isolates belonging mainly to STs 57, 115, 354, 393, and 420, and phylogroup B2 isolates belonging to the international clone O25b-ST131. Diagn Microbiol Infect Dis. 2010;67(3):270–6.
Sidjabat HE, Seah KY, Coleman L, Sartor A, Derrington P, Heney C, et al. Expansive spread of IncI1 plasmids carrying bla CMY-2 amongst Escherichia coli. Int J Antimicrob Agents. 2014;44(3):203–8.
Karczmarczyk M, Wang J, Leonard N, Fanning S. Complete nucleotide sequence of a conjugative IncF plasmid from an Escherichia coli isolate of equine origin containing bla CMY-2 within a novel genetic context. FEMS Microbiol Lett. 2014;352(1):123–7.
Carattoli A. Resistance plasmid families in Enterobacteriaceae. Antimicrob Agents Chemother. 2009;53(6):2227–38.
Deng H, Si HB, Zeng SY, Sun J, Fang LX, Yang RS, et al. Prevalence of extended-spectrum cephalosporin-resistant Escherichia coli in a farrowing farm: ST1121 clone harboring IncHI2 plasmid contributes to the dissemination of bla CMY-2. Front Microbiol. 2015;6:1210.
Voets GM, Fluit AC, Scharringa J, Schapendonk C, van den Munckhof T, Leverstein-van Hall MA, et al. Identical plasmid AmpC beta-lactamase genes and plasmid types in E. coli isolates from patients and poultry meat in the Netherlands. Int J Food Microbiol. 2013;167(3):359–62.
Jahanbakhsh S, Smith MG, Kohan-Ghadr HR, Letellier A, Abraham S, Trott DJ, et al. Dynamics of extended-spectrum cephalosporin resistance in pathogenic Escherichia coli isolated from diseased pigs in Quebec. Canada Int J Antimicrob Agents. 2016;48(2):194–202.
Jamborova I, Dolejska M, Zurek L, Townsend AK, Clark AB, Ellis JC, et al. Plasmid-mediated resistance to cephalosporins and quinolones in Escherichia coli from American crows in the USA. Environ Microbiol. 2017;19(5):2025–36.
Bortolaia V, Hansen KH, Nielsen CA, Fritsche TR, Guardabassi L. High diversity of plasmids harbouring bla CMY-2 among clinical Escherichia coli isolates from humans and companion animals in the upper Midwestern USA. J Antimicrob Chemother. 2014;69(6):1492–6.
Dame-Korevaar A, Fischer EAJ, Stegeman A, Mevius D, van Essen-Zandbergen A, Velkers F, et al. Dynamics of CMY-2 producing E. coli in a broiler parent flock. Vet Microbiol. 2017;203:211–4.
Mo SS, Slettemeas JS, Berg ES, Norstrom M, Sunde M. Plasmid and host strain characteristics of Escherichia coli resistant to extended-spectrum cephalosporins in the Norwegian broiler production. PLoS One. 2016;11(4):e0154019.
Guynet C, de la Cruz F. Plasmid segregation without partition. Mob Genet Elements. 2011;1(3):236–41.
da Silva KC, Cunha MP, Cerdeira L, de Oliveira MG, de Oliveira MC, Gomes CR, et al. High-virulence CMY-2- and CTX-M-2-producing avian pathogenic Escherichia coli strains isolated from commercial turkeys. Diagn Microbiol Infect Dis. 2017;87(1):64–7.
Tagg KA, Iredell JR, Partridge SR. Complete sequencing of IncI1 sequence type 2 plasmid pJIE512b indicates mobilization of bla CMY-2 from an IncA/C plasmid. Antimicrob Agents Chemother. 2014;58(8):4949–52.
Huddleston JR. Horizontal gene transfer in the human gastrointestinal tract: potential spread of antibiotic resistance genes. Infect Drug Resist. 2014;7:167–76.
Shahada F, Sekizuka T, Kuroda M, Kusumoto M, Ohishi D, Matsumoto A, et al. Characterization of Salmonella enterica serovar typhimurium isolates harboring a chromosomally encoded CMY-2 beta-lactamase gene located on a multidrug resistance genomic island. Antimicrob Agents Chemother. 2011;55(9):4114–21.
Mac Aogain M, Rogers TR, Crowley B. Identification of emergent bla CMY-2 -carrying Proteus mirabilis lineages by whole-genome sequencing. New Microbes New Infect. 2016;9:58–62.
Hille K, Ruddat I, Schmid A, Hering J, Hartmann M, von Munchhausen C, et al. Cefotaxime-resistant E. coli in dairy and beef cattle farms-joint analyses of two cross-sectional investigations in Germany. Prev Vet Med. 2017;142:39–45.
Laube H, Friese A, von Salviati C, Guerra B, Kasbohrer A, Kreienbrock L, et al. Longitudinal monitoring of extended-spectrum-beta-lactamase/AmpC-producing Escherichia coli at German broiler chicken fattening farms. Appl Environ Microbiol. 2013;79(16):4815–20.
Von Salviati C, Friese A, Roschanski N, Laube H, Guerra B, Kasbohrer A, et al. Extended-spectrum beta-lactamases (ESBL)/AmpC beta-lactamases-producing Escherichia coli in German fattening pig farms: a longitudinal study. Berl Munch Tierarztl Wochenschr. 2014;127(9–10):412–9.
Hering J, Fromke C, von Munchhausen C, Hartmann M, Schneider B, Friese A, et al. Cefotaxime-resistant Escherichia coli in broiler farms-a cross-sectional investigation in Germany. Prev Vet Med. 2016;125:154–7.
Hering J, Hille K, Fromke C, von Munchhausen C, Hartmann M, Schneider B, et al. Prevalence and potential risk factors for the occurrence of cefotaxime resistant Escherichia coli in German fattening pig farms--a cross-sectional study. Prev Vet Med. 2014;116(1–2):129–37.
Wallmann J, Schroter K, Wieler LH, Kroker R. National antibiotic resistance monitoring in veterinary pathogens from sick food-producing animals: the German programme and results from the 2001 pilot study. Int J Antimicrob Agents. 2003;22(4):420–8.
Clermont O, Bonacorsi S, Bingen E. Rapid and simple determination of the Escherichia coli phylogenetic group. Appl Environ Microbiol. 2000;66(8):4555.
Eller C, Simon S, Miller T, Frick JS, Prager R, Rabsch W, et al. Presence of beta-lactamases in extended-spectrum-cephalosporin-resistant Salmonella enterica of 30 different serovars in Germany 2005-11. J Antimicrob Chemother. 2013;68(9):1978–81.
Grobner S, Linke D, Schutz W, Fladerer C, Madlung J, Autenrieth IB, et al. Emergence of carbapenem-non-susceptible extended-spectrum beta-lactamase-producing Klebsiella pneumoniae isolates at the university hospital of Tubingen. Germany J Med Microbiol. 2009;58(Pt 7):912–22.
Clowes RC, Rowley D. Some observations on linkage effects in genetic recombination in Escherichia coli K-12. J Gen Microbiol. 1954;11(2):250–60.
Barton BM, Harding GP, Zuccarelli AJ. A general method for detecting and sizing large plasmids. Anal Biochem. 1995;226(2):235–40.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Coil D, Jospin G, Darling AE. A5-miseq: an updated pipeline to assemble microbial genomes from Illumina MiSeq data. Bioinformatics. 2015;31(4):587–9.
Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, 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.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67(11):2640–4.
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.
Joensen KG, Scheutz F, Lund O, Hasman H, Kaas RS, Nielsen EM, et al. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol. 2014;52(5):1501–10.
Joensen KG, Tetzschner AM, Iguchi A, Aarestrup FM, Scheutz F. Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J Clin Microbiol. 2015;53(8):2410–26.
Junemann S, Sedlazeck FJ, Prior K, Albersmeier A, John U, Kalinowski J, et al. Updating benchtop sequencing performance comparison. Nat Biotechnol. 2013;31(4):294–6.
Ferdous M, Kooistra-Smid AM, Zhou K, Rossen JW, Friedrich AW. Virulence, antimicrobial resistance properties and phylogenetic background of non-H7 enteropathogenic Escherichia coli O157. Front Microbiol. 2016;7:1540.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.
Becker L, Fuchs S, Pfeifer Y, Semmler T, Eckmanns T, Korr G, et al. Whole genome sequence analysis of CTX-M-15 producing Klebsiella isolates allowed dissecting a polyclonal outbreak scenario. Front Microbiol. 2018;9:322.
Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26(5):589–95.
Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Suzuki S, Shibata N, Yamane K, Wachino J, Ito K, Arakawa Y. Change in the prevalence of extended-spectrum-beta-lactamase-producing Escherichia coli in Japan by clonal spread. J Antimicrob Chemother. 2009;63(1):72–9.
Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference Centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34:D32–6.
We thank Sibylle Müller-Bertling and Kirstin Ganske for excellent technical support. We further thank Dr. Birgit Strommenger for her support with SeqSphere+ analyses, Dr. Angelika Fruth for the laboratory confirmation of the shigatoxin gene stx-2 and Katja Hille for the excellent coordinative work for the national research consortium “RESET” (www.reset-verbund.de). We also thank the laboratory staff for collecting samples and performing primary diagnostics for AmpC-producing E. coli (see also “RESET Study Group”).
We also thank all members of the RESET Study Group: Anke Bühling (MVZ Gemeinschaftslabor Cottbus GbR, Cottbus, Germany), Burkhard Domurath (Werner Wicker Klinik, Bad Wildungen, Germany), Christina von Salviati-Claudius (Freie Universität Berlin, Institute for Animal Hygiene and Environmental Health, Berlin, Germany), Constanze Wendt (Labor Limbach, Heidelberg, Germany), Giuseppe Valenza (Bavarian Health and Food Safety Authority, Erlangen, Germany), Hans Günther Wahl (Medizinisches Labor Wahl, Lüdenscheid, Germany), Helene M. Lu (Freie Universität Berlin, Institute for Animal Hygiene and Environmental Health, Berlin, Germany), Henriette Laube (Freie Universität Berlin, Institute for Animal Hygiene and Environmental Health, Berlin, Germany), Johanna Hering (University of Veterinary Medicine Hannover, Institute for Biometry, Epidemiology and Information Processing, Hannover, Germany), Josef Hellkamp (Klinikum Oldenburg, Institut für Laboratoriumsdiagnostik und Mikrobiologie, Oldenburg, Germany), Mardjan Arvand (Hessisches Landesprüfungs- u. Untersuchungsamt im Gesundheitswesen, Dillenburg, Germany, present address: Robert Koch Institute, FG 14 Hospital Hygiene, Infection Prevention and Control, Berlin, Germany), Michael Kresken (Antiinfectives Intelligence GmbH, Rheinbach, Germany and University of Applied Sciences, Cologne, Germany), Rainer Podschun (Christian-Albrechts-Universität zu Kiel, Institut für Infektionsmedizin, Kiel, Germany), Sandra Schneider (Charité - Institut für Hygiene und Umweltmedizin, Berlin, Germany), Steffen Tobisch and Marc Witt (MVZ Laborverbund GmbH, Berlin, Germany), Tim Eckmanns (Robert Koch Institute, Berlin, Germany), Ulrike Hachmann (MVZ Labor Münster, Münster, Germany), Uwe Bührlen (Labor Prof. Gisela Enders MVZ GbR, Stuttgart, Germany).
This work was partly funded by a grant of the Federal Ministry of Education and Research, Germany, grant no. 01KI1013 to MP, AI, NR, and GBM. Sample collection, study design, analysis, data interpretation and manuscript writing was carried out by the authors.
Availability of data and materials
Sequence data have been submitted to the European Nucleotide Archive (ENA) under the accession number PRJEB23663.
Ethics approval and consent to participate
Ethical approval was not necessary since patient data were not included in this project. All isolates from animals, farm environment and food were collected in the scope of other studies. Since these isolates originated from faeces samples and environmental samples, no ethical approval was necessary.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics of 164 CMY-2 -producing Escherichia coli isolates from different sources, 2008–2013, Germany (selection for whole genome sequencing). (XLSX 26 kb)
Figure S1. Neighbour-joining tree of ST429 E. coli isolates based on an ad-hoc cgMLST including 2547 alleles. The tree was built with SeqSphere+. Included isolate sequences originated from this study and were obtained as contigs assembled by http://enterobase.warwick.ac.uk/ (ERR2091318, ERR209121, ERR2091324, ERR2091328, ERR2091334, ERR2091342, ERR2091342, ERR2091348, ERR2091349, ERR2091357, ERR2091358, ERR2091419, ERR2091421, DRR102690, SRR3050857, SRR3098809, SRR3987496, ERR1619552, ERR1622238, ERR1622239, ERR1622406, ERR1622406, ERR1622407, ERR1595423, ERR1543414, ERR277049, ERR1415546, ERR1163310, ERR435146, SRR2000414). All isolates were investigated for the presence of blaCMY-2 and IncK2 RNAI sequence. Positive isolates are marked by red border. All blaCMY-2 carrying isolates exhibited a p486–16-like IncK2 plasmid sequence. The allele distance between two blaCMY-2 -carrying isolates is shown. (PDF 580 kb)
Figure S2. SNP-based maximum-likelihood-trees of isolates of selected sequence types. a: Maximum-likelihood-tree of ST117 and ST3778 isolates. SNPfilter (d = 0, reference: NZ_CP019903.1 E. coli strain MDR_56) based tree, calculated with RAxML, GTR Gamma and rapid hill-climbing and 100 starting trees. b: Maximum-likelihood-tree of ST131 O25b:H4 fimH22 isolates. SNPfilter (d = 0, reference: NZ_CP019008.1 E. coli strain Ecol_AZ159) based tree, calculated with RAxML, GTR Gamma and rapid hill-climbing and 100 starting trees. c: Maximum-likelihood-tree of ST1196 isolates. SNPfilter (d = 0, reference: NC_020518.1 E. coli str. K-12 substr. MDS42) based tree, calculated with RAxML, GTR Gamma and rapid hill-climbing and 100 starting trees. d: Maximum-likelihood-tree of ST429 isolates. SNPfilter (d = 0, reference: NC_013654.1 E. coli strain SE15) based tree, calculated with RAxML, GTR Gamma and rapid hill-climbing and 100 starting trees. (PDF 174 kb)
Figure S3. Minimum spanning tree of E. coli ST131 isolates from this study and INSDC and Enterobase based on an ad-hoc cgMLST including 2547 alleles. All fimH alleles and blaCMY and blaCTX-M genes are color-coded. Isolates from this study are marked with a red ring. (PDF 971 kb)
Figure S4. Comparison of plasmid sequence p486–16 with other IncK2 plasmids created with EasyFig v.2.2.2 (http://mjsull.github.io/Easyfig/). Used reference plasmid sequence were pTMSA1088 (Genbank: KR905386.1), pDV45 (KR905384.1). (PDF 38 kb)
Table S2. Unique integration sites for MGEs in IncK and IncI plasmids of the present collection of 164 CMY-2 producing E. coli isolates of various origins. (XLSX 15 kb)
Table S3. List of chromosomal core genome genes used. (XLSX 136 kb)
About this article
Cite this article
Pietsch, M., Irrgang, A., Roschanski, N. et al. Whole genome analyses of CMY-2-producing Escherichia coli isolates from humans, animals and food in Germany. BMC Genomics 19, 601 (2018). https://doi.org/10.1186/s12864-018-4976-3