Skip to main content

Comparative genomics of quinolone-resistant Escherichia coli from broilers and humans in Norway

Abstract

Background

The usage of fluoroquinolones in Norwegian livestock production is very low, including in broiler production. Historically, quinolone-resistant Escherichia coli (QREC) isolated from Norwegian production animals rarely occur. However, with the introduction of a selective screening method for QREC in the Norwegian monitoring programme for antimicrobial resistance in the veterinary sector in 2014; 89.5% of broiler caecal samples and 70.7% of broiler meat samples were positive. This triggered the concern if there could be possible links between broiler and human reservoirs of QREC. We are addressing this by characterizing genomes of QREC from humans (healthy carriers and patients) and broiler isolates (meat and caecum).

Results

The most frequent mechanism for quinolone resistance in both broiler and human E. coli isolates were mutations in the chromosomally located gyrA and parC genes, although plasmid mediated quinolone resistance (PMQR) was also identified. There was some relatedness of the isolates within human and broiler groups, but little between these two groups. Further, some overlap was seen for isolates with the same sequence type isolated from broiler and humans, but overall, the SNP distance was high.

Conclusion

Based on data from this study, QREC from broiler makes a limited contribution to the incidence of QREC in humans in Norway.

Peer Review reports

Introduction

Fluoroquinolones are a group of broad-spectrum antimicrobials that have been used for treatment of a wide range of infections in both humans and animals. In humans there are multiple approved indications for the use of fluoroquinolones, among them bone and joint infections, gastrointestinal infections, urinary tract infections (UTI) and respiratory tract infections [1]. However, the human use of this group of antimicrobials is very low in Norway compared to other countries, and primarily restricted to patients with complicated urinary tract or intra-abdominal infections [2,3,4]. Fluoroquinolones have since 2019 been classified by the World Health Organization as critically important for human treatment [5], and have been restricted for use in animals by the EU Antimicrobial Advise ad hoc Expert Group (AMEG) since 2020 [6]. From 1993 to 2021 the proportion of sales of fluoroquinolones for food-producing terrestrial animals in Norway varied between 0.1% and 0.3%, while in humans the quinolone usage has decreased by more than 60% since 2012 [7]. This restrictive use of fluoroquinolones is reflected through low occurrence of quinolone-resistant Escherichia coli (QREC) in the national surveillance data on antimicrobial resistance, i.e. in the monitoring programme for antimicrobial resistance in bacteria from feed, food and animals (NORM-VET) and in the surveillance programme for antimicrobial resistance in human pathogens (NORM).

Quinolone resistance among E. coli from broilers investigated in NORM-VET has historically been very low, though with an increase from 3.4% in 2014 to 12.6% in 2020 [8,9,10,11]. These data are based on a non-selective method, where one indicator/random E. coli isolate is selected from one pooled caecal sample from ten broilers per flock, and then further susceptibility tested and interpreted using epidemiological cut-off values (ECOFFs) as defined by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) [12]. In E. coli isolates from human blood and UTIs in Norway in 2014, the occurrence of fluoroquinolone resistance based on clinical breakpoints was 12.6% and 8.7% of all isolates tested, respectively [8]. Since then, quinolone resistance seems to have been relatively stable between 10% and 13% in E. coli from blood and between 7.5% and 9% in E. coli from urine samples [13].

Fluoroquinolones act by inhibiting DNA synthesis of two essential type II topoisomerases, DNA gyrase and topoisomerase IV [14,15,16]. The mechanisms of quinolone resistance are either chromosomal point mutations in genes causing a reduced affinity to protein targets, reduced accumulation of drugs either by a decrease in uptake or an increase in efflux of the agent, or plasmid-mediated quinolone resistance (PMQR) genes encoding target protection [14]. Chromosomal point mutations are located within a small quinolone resistance-determining region (QRDR) of the genes gyrA, gyrB, parC and parE. The level of resistance increases with the number of mutations in these genes [14]. For E. coli, the most amino acid substitutions are at the Ser83 (S83) and Asp87 (D87) in GyrA and at Ser80 (S80) and Glu84 (E84) in ParC. The product of plasmid-located qnr genes, on the other hand, is a pentapeptide repeat family molecule that blocks the action of ciprofloxacin on DNA gyrase and topoisomerase IV [17, 18]. Several plasmid-located gene variants have been described; qnrA, qnrB, qnrC, qnrD, qnrE, qnrS, qnrVC, as well as aac(6’)-lb-cr, oqxAB and qepA [16, 19,20,21].

Some E. coli might cause extraintestinal disease such as UTIs, sepsis and respiratory tract infections and have been identified as extraintestinal pathogenic E. coli (ExPEC) [22]. A common method for characterizing E. coli is multi-locus sequence typing (MLST), molecular typing based on seven housekeeping genes, which is used for identifying and tracking specific sequence types (STs) that are associated with antimicrobial resistance and/or pathogenicity [22, 23].

Some E. coli STs have frequently been identified as human ExPEC including ST10, ST69, ST73, ST95, ST117, and ST131 [24, 25]. Moreover, some successful ExPEC STs, disseminated worldwide, have been detected as antimicrobial resistant lineages both from humans and from food/animal sources, like ST10, ST69, ST95, ST117, ST131, and ST405 [24, 26,27,28,29]. In a review from Manges [24], poultry and poultry meat was hypothesized as a reservoir for human ExPEC. This was further supported by a Swedish study [30, 31], which showed that identical or closely related AmpC or extended-spectrum beta-lactamase (ESBL)-producing E. coli assigned to ST10, ST38, ST69 and ST131 were found both on Swedish and imported chicken meat. The Swedish study indicates both that the introduction of the AmpC-/ESBL resistant E. coli strains in broiler was due to broiler import and that faecal contamination to meat occurring at slaughter represents a source for further spread to humans. Also, a recent Italian study described various human pandemic and emerging ExPEC STs such as ST10, ST23, ST69, ST95, ST117 and ST131 from poultry [28]. In that study, the quinolone resistance mechanisms detected in approximately 45–50% of the quinolone-resistant isolates were mainly due to chromosomal point mutations in gyrA, parC and/or parE, though qnrS1 and qnrB19 were also detected. Moreover, a study on fluoroquinolone-resistant avian pathogenic E. coli (APEC) from Korea identified chromosomal point mutations in gyrA, parC and parE as the main reason behind the quinolone resistance mechanism, while qnrS1 was the only PMQR gene detected in only a few isolates [29]. Some of these isolates were assigned to major lineages of ExPEC such as ST95 (n = 3) and ST69 (n = 2). One or more mutations in gyrA, parC or parE were identified in these isolates, in addition, two of the ST95 isolates carried the qnrS1 gene.

The aim of the present study was to explore whether QREC from broilers in Norway may contribute to the occurrence of QREC in humans. We collected, whole genome sequenced and compared isolates that were phenotypically resistant to quinolones. They were from both humans and broilers, 100 isolates from each. Broiler QREC isolates were selected in a manner that ensured representativeness among QREC in broilers in Norway. The human QREC isolates were from UTI or bacteraemia isolates, as well as isolated from stool samples from healthy carriers. We performed whole genome sequencing followed by MLST analysis, detection of genes associated with quinolone resistance and other antimicrobial resistance genes as well as phylogenetic analysis.

Results

Isolate overview

In total, 200 QREC genomes were sequenced and assembled. Ten of the human isolates were excluded, seven because they were phenotypically susceptible to ciprofloxacin and/or nalidixic acid, two due to low quality of sequences/reads, and one that turned out to be Klebsiella variicola (Mash Screen). The reads as well as the annotated assemblies were deposited in the European Nucleotide Archive (ENA) repository under study accession numbers PRJEB36302 and PRJEB33048 and at the Sequence Read Archive (SRA) under BioProject accession number PRJNA1117742 (Additional file 1). On average, the assemblies contained 111 contigs (min: 39, max: 491), and had an average N50 size of 173,710 bp (min: 57,055, max: 361,993). Details regarding the quality of the genome data are available in the MultiQC and Quast reports provided in Additional file 2 and 3, respectively.

All 190 isolates included in the analysis were either resistant to ciprofloxacin (R > 0.06 mg/L) and/or nalidixic acid (R > 8 mg/L) according to ECOFFs and inclusion criteria set for the study. The overall MIC distributions and the antimicrobial resistance profiles for both the broiler and human isolates are shown in Additional file 1. In total, 68 of the 90 human isolates (75.6%) were considered multidrug resistant (MDR) (i.e. resistant to three or more antimicrobial classes), of which 42 were isolated from bacteraemia and twelve from UTI cases. Among the 100 broiler isolates 23 (23%) were considered MDR, of which eleven were isolated from meat and twelve from caecum.

Quinolone resistance genes

Point mutations

Table 1 gives an overview of the detected chromosomal point mutations causing quinolone resistance in all the 190 QREC isolates. A detailed description of the mutations in each isolate is shown in Additional file 1.

Table 1 Overview of the chromosomal point mutations causing quinolone resistance by amino acid substitution in the Escherichia coli isolates (n = 190) from broilers and from humans

Overall, amino acid substitutions were more diverse in the human than in the broiler QREC isolates (16 different substitutions in all four genes and six different substitutions altogether in two genes, respectively). In broiler QREC isolates compared to human isolates, point mutations in one gene were more common (gyrA (n = 76) or parE (n = 1)) followed by two genes (gyrA + parE (n = 23)). The human isolates were more diverse when it came to point mutations: no mutations (n = 9), point mutations in gyrA only (n = 52), point mutations in two genes (gyrA + parC (n = 17) or gyrA + parE (n = 4)), or point mutations in three genes (n = 8). Interestingly, none of the broiler isolates had substitutions in gyrB nor parC. For both broiler and human isolates, the amino acid substitution S83L in GyrA were most frequent (n = 96 and n = 77, respectively). Thereafter, the most frequent amino acid substitutions D475E in ParE (n = 22) and D87N in GyrA (n = 11) among the broiler isolates, and D87N in GyrA (n = 26), S80I in ParC (n = 15) and S58I in ParC (n = 15) among the human isolates.

Acquired genes

In total, five plasmid-located quinolone resistance genes were detected; five in broiler and 13 in human QREC isolates (Table 2). Two of the genes (qnrB19 and qnrS1) were found in both human and poultry isolates, while the remaining three genes (aac(6’)-Ib-cr, qnrB1 and qnrB5) were found in human isolates, only. An overview of the total acquired genes detected in the QREC isolates are shown in Additional file 1.

Of the five broiler isolates with acquired quinolone resistance genes, four carried the qnrS1 gene together with the amino acid substitution S83L in GyrA, while the last one carried the qnrB19 gene together with the same amino acid substitution S83L in GyrA. The qnrS1 gene was also detected in five (35.7%) of the 13 human isolates, but only one had the amino acid substitution S83L in GyrA, and none had the same sequence type (ST) as the qnrS1 carrying broiler isolates. Four (28.6%) of the human isolates carried the qnrB19 gene, of which one had three amino acid substitutions at S83L and D87N in GyrA and S80I in ParC, and none had the same ST as that of the qnrB19 carrying broiler isolates.

For acquired genes that were exclusively found in human isolates, the aac(6’)-lb-cr gene was detected in three (21.4%) isolates. One isolate had the additional amino acid substitution S83L in GyrA and another one had two amino acid substitutions (S83L in GyrA and S58I in ParC), and one isolate contained qnrB1 but no additional amino acid substitutions. The last isolate contained the qnrB5 gene together with three amino acid substitutions at S83L and D87N in GyrA and S80I in ParC.

Table 2 Quinolone resistance mechanisms and sequence types among broiler (n = 5) and human (n = 13) isolates carrying acquired quinolone resistance genes

Multi-locus sequence types (MLST) and phylogenetic analysis

In total, the QREC isolates were assigned into 62 unique E. coli MLSTs, as illustrated in Fig. 1. The human QREC isolates were assigned into 38 STs, where 62 isolates grouped into ten different STs and 28 were singletons. Broiler QREC isolates where assigned to 32 STs, 81 isolates grouped into 13 different STs and 19 were singletons. The most prevalent E. coli STs were ST131 (n = 26, 22 human, 4 broiler), ST355 (n = 22, all broiler), ST10 (n = 16, 4 human, 12 broiler), ST69 (n = 13, all human), ST162 (n = 12, 4 human, 8 poultry), and ST453 (n = 9, 1 human, 8 poultry). In total, 30 different STs were exclusively found in human isolates, 24 were found in broiler isolates, only, while eight STs were identified in both species (STs 453, 155, 1642, 162, 57, 10, 117 and 131).

A maximum likelihood tree based on core gene SNPs was made (Fig. 2) to further illustrate the relationship between isolates. The pangenome analysis of the 190 genomes detected 19 853 unique genes. Of these, 3156 were defined as core genes, i.e. present in at least 95% of the genomes. As expected, isolates with the same ST predominantly clustered closely together, with the notable exception of ST10. The STs that contained isolates from both species did not display any within-cluster clade separation between isolates stemming from the two species.

Fig. 1
figure 1

eBURST analysis. Minimum Spanning Tree of all MLSTs using PHYLOViZ v 2.0 showing the clonal relationship of broiler (n = 100) and human (n = 90) quinolone-resistant Escherichia coli isolates. Size of the circle corresponds to the number of isolates. The colour of the circle borders corresponds to group founder (light green) and common node (light blue)

Fig. 2
figure 2

Maximum likelihood SNP tree of 100 broiler and 90 human quinolone-resistant Escherichia coli isolates. Branch supports (ultrafast bootstrap approximation) are denoted with black or white nodes. The coloured tips and shapes on the tree denote species of origin and material. The tip labels denote the sequence types from the MLST scheme hosted by EnteroBase. The colouring on the outer rings denotes the presence/absence of mutations leading to amino acid (AA) substitutions in chromosomal genes and the presence/absence of plasmid-mediated genes leading to quinolone resistance. The grey shades denote groups of STs that contains isolates from both hosts, and which contains five or more isolates in total, ST10, ST117, ST131, ST155, ST162, and ST453. The tree was generated with IQTree from SNPs in core genes from Roary aligned with MAFFT. The evolutionary model used was GTRFASCR9 - GTRGAMMA. The tree is midpoint rooted for better visualization

Core genome phylogeny on selected sequence types

Core genome cluster analyses was performed for the STs that contained either five isolates or more, or isolates from both poultry and humans; they were ST10, ST117, ST131, ST155, ST162 and ST453 (Fig. 3). Within this detailed analysis, human and broiler isolates produced distinctly separate clades with the same resistance mechanisms present for all STs, except for ST162. The average genome coverage, SNP range and SNP distance are described for each of the core genome alignments in Table 3.

Table 3 Phylogenetic analysis per sequence type including quinolone-resistant Escherichia coli isolates from both broiler and humans

Within ST162, one of the human QREC isolates clustered with two broiler isolates. This subcluster had an SNP-range of 16–33, a mean SNP distance of 24.7 and a median SNP distance of 25.

Fig. 3
figure 3

Maximum likelihood core gene midpoint rooted SNP trees for clades ST10, ST117, ST131, ST155, ST162 and ST453. These six clades comprise quinolone-resistant Escherichia coli isolates from both humans and broiler. Core genome SNPs were identified with ParSNP, and recombinant sites were removed with Gubbins. ST - sequence type. Genotypes explaining the quinolone resistance mechanism are represented in the heat map

Discussion

In this study, we explored to what extent the QREC population from broilers in Norway could have contributed to QREC in humans. To this end, a dataset comprising QREC isolates from broiler meat and caecum, and from healthy humans and from patients with bacteraemia or UTIs was sequenced and compared. In a large-scale phylogenetic comparison based on gene alignment, it was clear that there was no substantial intermingling between broiler and human isolates. However, some STs did contain isolates that could seem similar, and they were further explored with a phylogenetic tree based on whole genome alignments, applying a threshold of around 20 SNPs, which is a commonly used threshold for source attribution [32]. As is clear from several of the subtrees created, only two contained subgroups of isolates that came close to the criterion; that is ST162 and ST453. When exploring relatedness of isolate sequences there should not only be a low number of SNP differences, but the isolates in question should also be found inside a monophyletic group [32]. This is clearly not the case for ST453. On the other hand, ST162 does have one such group, where one of the human ST162 QREC isolates clustered closely together with two of the broiler isolates (one from caecum and one from meat, with 25 SNP differences to the closest broiler isolate). A limitation that is often left out when comparing sequences, though is how much of the genomes are actually included [33]. For the STs comprising both human and broiler isolates, the average genome coverage ranged from 75.2 to 89.3%. The ST containing a potential monophyletic group with low SNP distance had the highest genome coverage. This means that more than 10% of the genome was not similar enough to be included in the comparison. E. coli ST162, which had the highest average genome coverage of 89.3%, included four human and nine broiler QREC isolates. Our finding is supported by a British study where fluoroquinolone-resistant E. coli from cattle were compared with those causing bacteriuria in humans living in the same geographical area [34]. They identified possible farm-to-human sharing of ST744 and ST162 but concluded that they had limited impact on community bacteriuria.

However, a clear limitation to our study is that we have no data on level of exposure among humans to broilers or chicken meat. In addition, we did not perform further analyses on these strains collected from 2012 to 2016 regarding PMQR genes, plasmid characterization and more detailed phylogeny, since chromosomal point mutations appeared to be the major mechanism behind quinolone resistance in E. coli from both human and broiler.

The results in the present study did not point to broilers as a major source for QREC in humans. Both differences and similarities in STs between humans and broiler were however, detected. There was a higher diversity of STs in human QREC isolates as compared to those from broiler. Also, within some STs, such as ST355, ST453, ST131, and ST393, quinolone resistance mechanisms seemed to be more uniform than for other STs. QREC isolates assigned to ST355 and ST453 all, but one human ST453 isolate, derived from broiler. ST131 and ST393 QREC isolates derived solely from human isolates, except four broiler meat QREC isolates that were assigned in ST131. Similar findings of E. coli ST355 were also documented by Röderova et al. [35], who reported ST355 to be one of the major clones prevalent in Czech retail turkeys carrying qnrS1 and qnrS1 together with qnrB19. In contrast to our study, Ramadan et al. [36], found that E. coli isolates from chicken in Egypt were more diverse regarding STs compared to E. coli isolates from humans and beef. Our findings are not surprising since the broiler production in Norway is dependent on import of breeding animals originating from a common ancestor [37].

Some of the STs identified in this study, ST10, ST69, ST117 and ST131, belong to successful clones disseminated worldwide and have been detected in antimicrobial resistant lineages of ExPEC from humans and from food/animal sources [22, 24, 26,27,28,29, 38]. Human E. coli isolates including a sublineage of ST131, namely ST131-H30, are known to be multidrug resistant (MDR) [22, 25, 38]. This is somewhat consistent with our study where 15 of the human ST131 QREC isolates were considered MDR, of which two of the four ST131-H30 found were MDR. Five ST131 QREC isolates were found in healthy carriers where three were identified as sublineage ST131-H30 (data not shown). Manges [24] suggested that E. coli ST131 from broiler may be a possible reservoir for human ExPEC, however, in our study this does not seem to be the case as they are phylogenetically distant and the resistance patterns are different. In contrast to our study, a comparable Algerian study demonstrated clonal relationship between ciprofloxacin-resistant E. coli ST131 strains in samples from human and broiler [39]. Liu et al. [40] investigated broiler meat samples and clinical cultures for one year and their findings suggested that a ST131-H22 sublineage has been established in food of animal origin. This corresponds well to our study where the four broiler isolates belonged to ST131-H22 while most of our ST131 human isolates belonged to ST131-H41 (n = 17 of 22 ST131 from human) (data not shown). Interestingly, the ST131-H41 sublineage has been described in E. coli isolated from UTIs and related to extended-cephalosporin resistance [41, 42].

When considering the main mechanisms behind quinolone resistance in both broiler and human QREC isolates, the most frequently seen were mutations in the QRDR of the chromosomally located gyrA and parC genes, though PMQR was also identified. Our results are in concordance with the findings in Börjesson et al. [43] who investigated imported breeding animals to the Swedish broiler production chain and identified that the main mechanism for QREC was due to a single mutation in gyrA. Others have reported that chromosomal point mutations are the main mechanisms of QREC in the broiler production chain [28, 44,45,46]. Studies on other animals such as cats, dogs, pigs and cattle support that chromosomal point mutations in QRDR is the most commonly identified key mechanism [47,48,49,50]. However, in some of these studies PMQR genes were also identified, most frequently a qnr gene variant [47, 48, 50].

Only a few studies have characterized PMQR genes from E. coli from broilers in Europe. In our study, PMQR genes were detected in QREC from 5% of the broiler samples, qnrS1 (n = 4, all ST453) and qnrB19 (n = 1, all ST349), and always together with the point mutation GyrA S83L. Similar findings in broilers have been described in an Italian study of E. coli isolated from poultry flocks with colibacillosis [51]. In contrast to our study on broilers, a Czech study on turkeys detected PMQR genes, mainly qnrB and qnrS, in as much as 58% of the E. coli [52]. A recent study from Nigeria detected, similar to our study, no clonal relationship between commensal PMQR carrying E. coli isolates from poultry, poultry workers, and poultry farms and market environments [53].

Conclusion

Based on data from this study, QREC from broiler makes a limited contribution to the incidence of QREC in humans in Norway. There was some relatedness between QREC isolated from humans and broilers, but little was shared across the two groups. Still, some overlap was seen of QREC isolates with the same STs isolated from broiler and humans, but overall, the SNP distance was too high to suggest substantial transfer. Some of the broiler QREC isolates with the same STs as human isolates, belonged to phylogroups that were closely related to human ExPEC, namely B2 and D, but they were not closely related. Only one human QREC isolate, E. coli ST162 carrying GyrA S83L substitution, clustered closely to two broiler QREC isolates. Further, the major quinolone resistance mechanisms were chromosomal point mutations in gyrA and parC. Overall, only a few PMQR genes were detected in the included isolates. These results demonstrate that the human and broiler QREC isolates are not highly related, and that the high occurrence of quinolone resistance is not horizontally transferable. In four of the five broiler QREC isolates carrying PMQR genes, the qnrS1 gene was detected in E. coli ST453. More studies are needed using newer technology like long read sequencing to look at the location of these qnr-genes and possible plasmid epidemiology to see if the same plasmids can be found in both human and broiler isolates.

Materials and methods

Bacterial isolates

Broiler isolates

A selective method for screening for quinolone-resistant E. coli (QREC) was implemented on broiler samples in NORM-VET in 2014, resulting in QREC isolates available from 89.5% of caecal flock samples and 70.7% of chicken meat samples. This collection of isolates consisted of 188 caecal and 140 meat isolates [8] and was susceptibility tested under the auspices of the NORM-VET programme. The isolates were classified as QREC with MICs for nalidixic acid and/or ciprofloxacin above ECOFF (EUCAST accessed 04.05.2022, i.e. ciprofloxacin MIC > 0.06 mg/L and/or nalidixic acid MIC > 8 mg/L).

The phylogroup of all 328 isolates was found using the method by Clermont et al. [54] and as described in Mo et al. [55]. Based on phylogroups, MIC values for ciprofloxacin and nalidixic acid, and resistance profile beyond quinolones, the isolates were divided into 86 different groups. Isolates were selected from this list based on their proportional distribution in the dataset, and to have at least one isolate per group. Overall, 47 of the 86 groups had only one isolate, all of which were included. The remaining 53 isolates were randomly chosen from the remaining groups proportionally to the frequency of isolates per group. A total of 100 QREC isolates from chicken meat (n = 47) and broiler caecal flock samples (n = 53) were included for further analyses in the present study. The broiler isolates metadata is described in Additional file 1.

Human isolates

A total of 100 human QREC isolates were included. The collection consisted of isolates from UTI and bacteraemia cases (n = 77) in 2012–2014 and retrieved through NORM. Further, 23 isolates were included from a study performed at the Norwegian Institute for Public Health (NIPH) from 2014 to 2016 where healthy carriers were screened for the occurrence of QREC [56].

The human isolates were originally susceptibility tested by the disk diffusion method. For full comparability to the broiler isolates, human isolates were retested using the broth microdilution method applied in NORM-VET, and with the Sensititre™ TREK EUVSEC plate (Trek Diagnostic System Ltd., United Kingdom). The EUVSEC panel contains the following 14 antimicrobial agents: ampicillin, azithromycin, cefotaxime, ceftazidime, chloramphenicol, ciprofloxacin, colistin, gentamicin, meropenem, nalidixic acid, sulfamethoxazole, tetracycline, tigecycline, and trimethoprim. The fully susceptible E. coli ATCC 25,922 was included as quality control. ECOFFs (EUCAST accessed 04.05.2022) were used to categorize the isolates as susceptible or resistant. The human isolates metadata is described in Additional file 1.

Whole genome sequencing

Bacteria were inoculated in Luria Bertani (LB) broth (Merck, Darmstadt, Germany) and incubated at 37 °C for 19–21 h in a shaker at 180 rpm. Genomic DNA was extracted through either automatic DNA extraction with QIASymphony DSP DNA Mini Kit or by manual extraction with DNeasy Blood and Tissue kit (both QIAGEN, Hilden, Germany), and according to the producer’s protocols. The DNA concentration and quality were measured using a Qubit™ 4 Fluorometer with Qubit™ dsDNA BR (Broad-Range) Assay kit (Thermo Fisher Scientific) and NanoDrop™ One Spectrophotometer (Thermo Fisher Scientific).

Library preparation and sequencing of the isolates were done at the Norwegian Sequencing Centre [57] using Illumina™ Nextera XT library prep (Illumina, Inc., San Diego, California, USA) on an Illumina™ HiSeq 2500 (Illumina) with Rapid Run spiked with PhiX generating 250 bp paired-end reads with four reads per isolate.

Assembly and characterization

Analyses were performed using the Bifrost pipeline (https://doi.org/10.5281/zenodo.4043861). Specifically, quality control was done using FastQC v0.11.9 [58] (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and MultiQC v1.9 [59] to collate the FastQC data. Contamination screening was performed using Mash Screen v2.1 [60]. Genome assembly was performed using Trimmomatic v0.39 [61], followed by SPAdes assembly v3.14.0 [62] with parameters coverage cutoff set to “auto”, “--careful” settings and excluding contigs shorter than 500 nucleotides. The SPAdes assemblies were further run through Pilon v1.23 [63] for finishing, before genome annotation with Prokka v1.14.5 [64]. QUAST v5.0.2 [65] was used for evaluating the assemblies.

Multi-locus sequence typing (MLST) and detection of antimicrobial resistance genes were done using ARIBA [66]. The scheme hosted by Enterobase was used for MLST [23], while the MEGARes [67] and ResFinder (accessed 28 March 2019) [68, 69] databases were used for the detection of ARGs (covering both chromosomal point mutations and acquired, respectively).

Phylogenetic analysis

An eBURST analysis [70, 71] of the MLST results was performed to do a preliminary investigation of the clonal relationship of the QREC isolates using PHYLOViZ v 2.0 [72]. Further, the genomes were used as input for pangenome analyses using the core gene track in the phylogeny pipeline ALPPACA v0.4.1 [73]. In short, the core genes were identified in the assemblies by running Prokka v1.14.5 [64]. Panaroo v1.2.2 [74] was run in sensitive mode to build a graphical representation of the pangenome that were used to error correct the annotated genomes. SeqKit [75] was used to remove duplicated genomes and SNPs were filtered with SNP-sites v2.5.1 [76]. A maximum likelihood (ML) tree was generated using IQtree v1.6.12 [77, 78], and snp-dists v0.6.3 (https://github.com/tseemann/snp-dists) was used to calculate SNP distances.

Genomes belonging to ST groups with more than five isolates from both human and broilers were analysed separately using the core genome track of the ALPPACA v0.4.1 pipeline. In short, the core genome was identified using ParSNP v1.5.3 [79] and duplicated genomes were removed using SeqKit [75]. Gubbins v2.4.1 [80] was used to identify recombinant areas that were further masked with maskrc-svg v0.5 (https://github.com/kwongj/maskrc-svg) before filtering with SNP-sites v2.5.1 [76] and generating the ML trees in IQtree v1.6.12 [77, 78].

Data management

The results from the antimicrobial resistance genes and MLST analysis was collated in R v4.0.5 to produce a summary report. Figures and tables were generated using R v4.0.5 [81], the ggtree package v2.4.2 [82,83,84] and ggtreeExtra v1.0.4 [85].

Data availability

The dataset supporting the conclusions of this article is available in the European Nucleotide Archive (ENA) repository (https://www.ebi.ac.uk/ena/browser/home), under study accession numbers PRJEB36302 and PRJEB33048, and at Sequence Read Archive (SRA) at National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/sra) under BioProduct accession number PRJNA1117742.

The dataset supporting the conclusions of this article is included within the article (and its additional file(s)).

References

  1. Bolon MK. The newer fluoroquinolones. Med Clin North Am. 2011;95(4):793–817. viii.

    Article  CAS  PubMed  Google Scholar 

  2. NORM/NORM-VET. NORM/NORM-VET. 2022. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway. Tromsø / Oslo; 2023. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  3. Norwegian Directory of Health. Antibiotics in primary care 2024 [ https://www.helsedirektoratet.no/retningslinjer/antibiotika-i-primaerhelsetjenesten

  4. Norwegian Directory of Health. Antibiotics in the hospital 2024 [ https://www.helsedirektoratet.no/retningslinjer/antibiotika-i-sykehus

  5. WHO. Critically important antimicrobials for human medicine, 6th revision. Licence: CC BY-NC-SA 3.0 IGO. World Health Organization; 2019.

  6. Categorisation of antibiotics used in. Animals promotes responsible use to protect public and animal health. [press release]. European Medicines Agency; 2020.

  7. NORM/NORM-VET. NORM/NORM-VET. 2021. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway. Tromsø / Oslo; 2022. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  8. NORM/NORM-VET. NORM/NORM-VET. 2014. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway. Tromsø / Oslo Tromsø / Oslo; 2015. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  9. NORM/NORM-VET. NORM/NORM-VET. 2016. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway. Tromsø / Oslo Tromsø / Oslo; 2017. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  10. NORM/NORM-VET. NORM/NORM-VET, Oslo. Tromsø / Oslo.; 2019. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  11. NORM/NORM-VET. NORM/NORM-VET. 2020. Usage of Antimicrobial Agents and Occurrence of Antimicrobial Resistance in Norway. Tromsø / Oslo; 2021. Report No.: ISSN:1502–2307 (print) / 1890–9965 (electronic).

  12. EUCAST. The European Committee on Antimicrobial Susceptibility Testing [Available from: www.eucast.org.

  13. NORM-atlas. [Available from: www.norm-atlas.no.

  14. Fabrega A, Madurga S, Giralt E, Vila J. Mechanism of action of and resistance to quinolones. Microb Biotechnol. 2009;2(1):40–61.

    Article  CAS  PubMed  Google Scholar 

  15. Vila J, Fabrega A, Roca I, Hernandez A, Martinez JL. Efflux pumps as an important mechanism for quinolone resistance. Adv Enzymol Relat Areas Mol Biol. 2011;77:167–235.

    Article  CAS  PubMed  Google Scholar 

  16. Ruiz J. Transferable mechanisms of Quinolone Resistance from 1998 onward. Clin Microbiol Rev. 2019;32(4).

  17. Tran JH, Jacoby GA. Mechanism of plasmid-mediated quinolone resistance. Proc Natl Acad Sci USA. 2002;99(8):5638–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Tran JH, Jacoby GA, Hooper DC. Interaction of the plasmid-encoded quinolone resistance protein QnrA with Escherichia coli topoisomerase IV. Antimicrob Agents Chemother. 2005;49(7):3050–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Martínez-Martínez L, Pascual A, Jacoby GA. Quinolone resistance from a transferable plasmid. Lancet. 1998;351(9105):797–9.

    Article  PubMed  Google Scholar 

  20. Guan X, Xue X, Liu Y, Wang J, Wang Y, Wang J, et al. Plasmid-mediated quinolone resistance - current knowledge and future perspectives. J Int Med Res. 2013;41(1):20–30.

    Article  CAS  PubMed  Google Scholar 

  21. Jacoby GA, Strahilevitz J, Hooper DC. Plasmid-mediated quinolone resistance. Microbiol Spectr. 2014;2(5).

  22. Manges AR, Geum HM, Guo A, Edens TJ, Fibke CD, Pitout JDD. Global Extraintestinal Pathogenic Escherichia coli (ExPEC) lineages. Clin Microbiol Rev. 2019;32(3).

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Manges AR. Escherichia coli and urinary tract infections: the role of poultry-meat. Clin Microbiol Infect. 2016;22(2):122–9.

    Article  CAS  PubMed  Google Scholar 

  25. Doumith M, Day M, Ciesielczuk H, Hope R, Underwood A, Reynolds R, et al. Rapid identification of major Escherichia coli sequence types causing urinary tract and bloodstream infections. J Clin Microbiol. 2015;53(1):160–6.

    Article  CAS  PubMed  Google Scholar 

  26. Jamborova I, Dolejska M, Vojtech J, Guenther S, Uricariu R, Drozdowska J, et al. Plasmid-mediated resistance to cephalosporins and fluoroquinolones in various Escherichia coli sequence types isolated from rooks wintering in Europe. Appl Environ Microbiol. 2015;81(2):648–57.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Stoesser N, Sheppard AE, Pankhurst L, De Maio N, Moore CE, Sebra R, et al. Evolutionary history of the global emergence of the Escherichia coli epidemic clone ST131. mBio. 2016;7(2):e02162.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Massella E, Giacometti F, Bonilauri P, Reid CJ, Djordjevic SP, Merialdi G, et al. Antimicrobial Resistance Profile and ExPEC Virulence potential in Commensal Escherichia coli of multiple sources. Antibiotics. 2021;10(4):351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kim JH, Lee HJ, Jeong OM, Kim DW, Jeong JY, Kwon YK, et al. High prevalence and variable fitness of fluoroquinolone-resistant avian pathogenic Escherichia coli isolated from chickens in Korea. Avian Pathology: J WVPA. 2021;50(2):151–60.

    Article  CAS  Google Scholar 

  30. Borjesson S, Egervarn M, Lindblad M, Englund S. Frequent occurrence of extended-spectrum beta-lactamase- and transferable ampc beta-lactamase-producing Escherichia coli on domestic chicken meat in Sweden. Appl Environ Microbiol. 2013;79(7):2463–6.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Egervarn M, Borjesson S, Byfors S, Finn M, Kaipe C, Englund S, et al. Escherichia coli with extended-spectrum beta-lactamases or transferable AmpC beta-lactamases and Salmonella on meat imported into Sweden. Int J Food Microbiol. 2014;171:8–14.

    Article  PubMed  Google Scholar 

  32. Pightling AW, Pettengill JB, Luo Y, Baugher JD, Rand H, Strain E. Interpreting whole-genome sequence analyses of Foodborne Bacteria for Regulatory Applications and Outbreak investigations. Front Microbiol. 2018;9:1482.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kravik IH, Kaspersen H, Sjurseth SK, Jonsson M, David B, Aspholm M, et al. High sequence similarity between avian pathogenic E. Coli isolates from individual birds and within broiler chicken flocks during colibacillosis outbreaks. Vet Microbiol. 2022;267:109378.

    Article  CAS  PubMed  Google Scholar 

  34. Mounsey O, Schubert H, Findlay J, Morley K, Puddy EF, Gould VC, et al. Limited phylogenetic overlap between fluoroquinolone-resistant < i > Escherichia coli isolated on dairy farms and those causing bacteriuria in humans living in the same geographical region. J Antimicrob Chemoth. 2021;76(12):3144–50.

    Article  CAS  Google Scholar 

  35. Roderova M, Halova D, Papousek I, Dolejska M, Masarikova M, Hanulik V, et al. Characteristics of Quinolone Resistance in Escherichia coli isolates from humans, animals, and the Environment in the Czech Republic. Front Microbiol. 2016;7:2147.

    PubMed  Google Scholar 

  36. Ramadan H, Jackson CR, Frye JG, Hiott LM, Samir M, Awad A, et al. Antimicrobial Resistance, genetic diversity and Multilocus sequence typing of Escherichia coli from humans, Retail Chicken and Ground Beef in Egypt. Pathogens. 2020;9(5):357.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Mo SS, Norström M, Slettemeås JS, Løvland 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.

    Article  CAS  PubMed  Google Scholar 

  38. Manges AR. Escherichia coli causing bloodstream and other extraintestinal infections: tracking the next pandemic. Lancet Infect Dis. 2019;19(12):1269–70.

    Article  PubMed  Google Scholar 

  39. Agabou A, Lezzar N, Ouchenane Z, Khemissi S, Satta D, Sotto A, et al. Clonal relationship between human and avian ciprofloxacin-resistant Escherichia coli isolates in North-Eastern Algeria. Eur J Clin Microbiol Infect Dis. 2016;35(2):227–34.

    Article  CAS  PubMed  Google Scholar 

  40. Liu CM, Stegger M, Aziz M, Johnson TJ, Waits K, Nordstrom L et al. Escherichia coli ST131-H22 as a Foodborne Uropathogen. mBio. 2018;9(4).

  41. Choi HJ, Jeong SH, Shin KS, Kim YA, Kim YR, Kim HS, et al. Characteristics of Escherichia coli urine isolates and risk factors for Secondary Bloodstream Infections in patients with urinary tract infections. Microbiol Spectr. 2022;10(4):e0166022.

    Article  PubMed  Google Scholar 

  42. Matsumura Y, Johnson JR, Yamamoto M, Nagao M, Tanaka M, Takakura S, et al. CTX-M-27- and CTX-M-14-producing, ciprofloxacin-resistant Escherichia coli of the H30 subclonal group within ST131 drive a Japanese regional ESBL epidemic. J Antimicrob Chemother. 2015;70(6):1639–49.

    Article  CAS  PubMed  Google Scholar 

  43. Börjesson S, Guillard T, Landén A, Bengtsson B, Nilsson O. Introduction of quinolone resistant Escherichia coli to Swedish broiler population by imported breeding animals. Vet Microbiol. 2015.

  44. Kmet V, Kmetova M. High level of quinolone resistance in Escherichia coli from healthy chicken broilers. Folia Microbiol (Praha). 2010;55(1):79–82.

    Article  CAS  PubMed  Google Scholar 

  45. Jones-Dias D, Manageiro V, Francisco AP, Martins AP, Domingues G, Louro D, et al. Assessing the molecular basis of transferable quinolone resistance in Escherichia coli and Salmonella spp. from food-producing animals and food products. Vet Microbiol. 2013;167(3–4):523–31.

    Article  CAS  PubMed  Google Scholar 

  46. Seo KW, Lee YJ. Molecular characterization of fluoroquinolone-resistant Escherichia coli from broiler breeder farms. Poult Sci. 2021;100(8):101250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. de Jong A, Muggeo A, El Garch F, Moyaert H, de Champs C, Guillard T. Characterization of quinolone resistance mechanisms in Enterobacteriaceae isolated from companion animals in Europe (ComPath II study). Vet Microbiol. 2018;216:159–67.

    Article  PubMed  Google Scholar 

  48. Kaspersen H, Sekse C, Zeyl Fiskebeck E, Slettemeås JS, Simm R, Norström M et al. Dissemination of Quinolone-Resistant Escherichia coli in the Norwegian broiler and Pig Production Chains and possible persistence in the broiler production environment. Appl Environ Microbiol. 2020;86(7).

  49. Hordijk J, Veldman K, Dierikx C, van Essen-Zandbergen A, Wagenaar JA, Mevius D. Prevalence and characteristics of quinolone resistance in Escherichia coli in veal calves. Vet Microbiol. 2012;156(1–2):136–42.

    Article  CAS  PubMed  Google Scholar 

  50. Yang F, Zhang S, Shang X, Wang L, Li H, Wang X. Characteristics of quinolone-resistant Escherichia coli isolated from bovine mastitis in China. J Dairy Sci. 2018;101(7):6244–52.

    Article  CAS  PubMed  Google Scholar 

  51. Niero G, Bortolaia V, Vanni M, Intorre L, Guardabassi L, Piccirillo A. High diversity of genes and plasmids encoding resistance to third-generation cephalosporins and quinolones in clinical Escherichia coli from commercial poultry flocks in Italy. Vet Microbiol. 2018;216:93–8.

    Article  CAS  PubMed  Google Scholar 

  52. Hricova K, Roderova M, Pudova V, Hanulik V, Halova D, Julinkova P, et al. Quinolone-resistant Escherichia coli in Poultry Farming. Cent Eur J Public Health. 2017;25(2):163–7.

    Article  CAS  PubMed  Google Scholar 

  53. Aworh MK, Kwaga JKP, Hendriksen RS, Okolocha EC, Harrell E, Thakur S. Quinolone-resistant Escherichia coli at the interface between humans, poultry and their shared environment- a potential public health risk. One Health Outlook. 2023;5(1).

  54. Clermont O, Christenson JK, Denamur E, Gordon DM. The Clermont Escherichia coli phylo-typing method revisited: improvement of specificity and detection of new phylo-groups. Environ Microbiol Rep. 2013;5(1):58–65.

    Article  CAS  PubMed  Google Scholar 

  55. Mo SS, Slettemeås JS, Berg ES, Norström 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.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Ulstad CR, Solheim M, Berg S, Lindbaek M, Dahle UR, Wester AL. Carriage of ESBL/AmpC-producing or ciprofloxacin non-susceptible Escherichia coli and Klebsiella spp. in healthy people in Norway. Antimicrob Resist Infect Control. 2016;5:57.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Norwegian Sequencing Centre. [Available from: www.sequencing.uio.no.

  58. Babraham Bioinformatics. FastQC [ http://www.bioinformatics.babraham.ac.uk/projects/fastqc/

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Ondov BD, Starrett GJ, Sappington A, Kostic A, Koren S, Buck CB, et al. Mash screen: high-throughput sequence containment estimation for genome discovery. Genome Biol. 2019;20(1):232.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biology: J Comput Mol cell Biology. 2012;19(5):455–77.

    Article  CAS  Google Scholar 

  63. Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE. 2014;9(11):e112963.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.

    Article  CAS  PubMed  Google Scholar 

  65. Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Hunt M, Mather AE, Sanchez-Buso L, Page AJ, Parkhill J, Keane JA, et al. ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads. Microb Genomics. 2017;3(10):e000131.

    Article  Google Scholar 

  67. Lakin SM, Dean C, Noyes NR, Dettenwanger A, Ross AS, Doster E, et al. MEGARes: an antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res. 2017;45(D1):D574–80.

    Article  CAS  PubMed  Google Scholar 

  68. Zankari E, Allesoe R, Joensen KG, Cavaco LM, Lund O, Aarestrup FM. PointFinder: a novel web tool for WGS-based detection of antimicrobial resistance associated with chromosomal point mutations in bacterial pathogens. J Antimicrob Chemother. 2017;72(10):2764–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004;186(5):1518–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Francisco AP, Bugalho M, Ramirez M, Carriço JA. Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinformatics. 2009;10:152.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Nascimento M, Sousa A, Ramirez M, Francisco AP, Carriço JA, Vaz C. PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods. Bioinformatics. 2017;33(1):128–9.

    Article  CAS  PubMed  Google Scholar 

  73. Kaspersen H, Fiskebeck EZ. ALPPACA - A tooL for Prokaryotic Phylogeny and Clustering Analysis. J Open Source Softw. 2022;7(79):4677.

    Article  Google Scholar 

  74. Tonkin-Hill G, MacAlasdair N, Ruis C, Weimann A, Horesh G, Lees JA, et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol. 2020;21(1):180.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Shen W, Le S, Li Y, Hu F, SeqKit:. A cross-platform and Ultrafast Toolkit for FASTA/Q file manipulation. PLoS ONE. 2016;11(10):e0163962.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Page AJ, Taylor B, Delaney AJ, Soares J, Seemann T, Keane JA et al. SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. Microb Genomics. 2016;2(4).

  77. Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74.

    Article  CAS  PubMed  Google Scholar 

  78. Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

  80. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43(3):e15–e.

    Article  PubMed  Google Scholar 

  81. Team RC. R: A Language and Environment for Statistical Computing. In: Computing RFfS, editor. Vienna, Austria2021.

  82. Yu G. Using ggtree to visualize data on Tree-Like structures. Curr Protocols Bioinf. 2020;69(1):e96.

    Article  Google Scholar 

  83. Yu G, Lam TT-Y, Zhu H, Guan Y. Mol Biol Evol. 2018;35(12):3041–3. Two Methods for Mapping and Visualizing Associated Data on Phylogeny Using Ggtree.

  84. Yu G, Smith DK, Zhu H, Guan Y, Lam TTY. GGTREE: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8(1):28–36.

    Article  Google Scholar 

  85. Xu S, Dai Z, Guo P, Fu X, Liu S, Zhou L, et al. ggtreeExtra: Compact visualization of richly annotated phylogenetic data. Mol Biol Evol. 2021;38(9):4039–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are thankful for to the hospital laboratories in Norway included in the NORM surveillance programme and the NORM-VET monitoring programme for access to the isolates, and volunteers to the healthy carrier study. Thank you to Charlotte Berg Andersen, Linn Fritzner Hansen, Iselin Høgbakk Ophus, June Gottenberg Skaalvik, and Kine Stenerud for performing the phylotyping on all broiler QREC isolates collected in the NORM-VET 2014 monitoring programme as part of their bachelor thesis. Further we would like to express our gratitude to the excellent work in the lab at NVI by Solveig Sølverød Mo and Aina Steihaug Barstad for assisting in collecting and performing the antimicrobial susceptibility testing of the QREC through NORM-VET. Finally, we are grateful for all help provided by Håkon Kaspersen at the NVI to finalize tables and figures in R.

Funding

The project was funded by the Norwegian Research Funding for Agriculture and Food Industry, grant number 244140 and the Norwegian Veterinary Institute.

Author information

Authors and Affiliations

Authors

Contributions

Jannice Schau Slettemeås (JSS), Camilla Sekse (CS), Marianne Sunde (MS), Madelaine Norström (MN), Astrid Louise Wester (ALW), Umaer Naseer (UN), Gunnar Skov Simonsen (GSS), Charlotte Rosenberg Ulstad (CRU), Anne Margrete Urdahl (AMU), Karin Lagesen (KL): JSS, CS, MS, GSS, ALW, AMU, and KL conceived, designed, and executed the study. JSS wrote the first draft of the manuscript. CS, JSS and MN were responsible for the selection of broiler strains and acquisition of general information on the broiler strains. GSS was responsible for the selection of human clinical strains. CRU performed the microbiology laboratory analyses of the human isolates. JSS and CS performed the microbiology laboratory analyses of the broiler isolates. JSS, CS and KL performed bioinformatics analyses. The data and results were analysed by JSS, CS and KL. All authors critically revised the manuscript, and JSS finalized the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jannice Schau Slettemeås.

Ethics declarations

Ethics approval and consent to participate

The healthy carrier study included both QREC and ESBLs and was approved by the Norwegian Regional Committees for Medical and Health Research Ethics (ID 2014/419). Clinical isolates were collected in accordance with the NORM Resistance Register Regulation. For the human isolates informed consent was obtained from all subjects and/or their legal guardian(s) as stated in previous studies involving the original sampling and bacterial isolation by Ulstad et al. (2016) [56]. Broiler isolates were collected under the auspices of the monitoring programme NORM-VET (Norwegian monitoring programme for antimicrobial resistance in bacteria from feed, food and animals) where caeca samples from broiler were taken at slaughter and meat samples were collected at retail. The NORM-VET monitoring programme is following the Commission Implementing Decision (EU) 2020/1729.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Slettemeås, J.S., Sekse, C., Sunde, M. et al. Comparative genomics of quinolone-resistant Escherichia coli from broilers and humans in Norway. BMC Microbiol 24, 248 (2024). https://doi.org/10.1186/s12866-024-03412-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03412-3

Keywords