Next Article in Journal
Characterization of Environmental and Cultivable Antibiotic-Resistant Microbial Communities Associated with Wastewater Treatment
Next Article in Special Issue
Antibiotic-Resistant Enterobacteriaceae in Wastewater of Abattoirs
Previous Article in Journal
Prognostic Value of Procalcitonin and C-Reactive Protein in 1608 Critically Ill Patients with Severe Influenza Pneumonia
Previous Article in Special Issue
Impact of On-Farm Interventions against CTX-Resistant Escherichia coli on the Contamination of Carcasses before and during an Experimental Slaughter
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 10
Issue 4
10.3390/antibiotics10040351
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources

by
Elisa Massella
1,
Federica Giacometti
1,*,
Paolo Bonilauri
2,
Cameron J. Reid
3,
Steven P. Djordjevic
3,
Giuseppe Merialdi
4,
Cristina Bacci
5,
Laura Fiorentini
6,
Paola Massi
6,
Lia Bardasi
4,
Silva Rubini
7,
Federica Savini
1,
Andrea Serraino
1 and
Silvia Piva
1
1
Department of Veterinary Medical Sciences, Faculty of Veterinary Medicine, University of Bologna, Ozzano Emilia, 40064 Bologna, Italy
2
Experimental Zooprophylactic Institute of Lombardy and Emilia Romagna, 42124 Reggio Emilia, Italy
3
The iThree Institute, University of Technology Sydney, Ultimo, NSW 2007, Australia
4
Experimental Zooprophylactic Institute of Lombardy and Emilia Romagna, 40127 Bologna, Italy
5
Department of Veterinary Medical Sciences, Faculty of Veterinary Medicine, University of Parma, 43126 Parma, Italy
6
Experimental Zooprophylactic Institute of Lombardy and Emilia Romagna, 47122 Forlì, Italy
7
Experimental Zooprophylactic Institute of Lombardy and Emilia Romagna, 44124 Ferrara, Italy
*
Author to whom correspondence should be addressed.
Antibiotics 2021, 10(4), 351; https://doi.org/10.3390/antibiotics10040351
Submission received: 8 March 2021 / Revised: 18 March 2021 / Accepted: 19 March 2021 / Published: 26 March 2021
(This article belongs to the Special Issue Antimicrobial Resistance: From Farm to Fork)
Browse Figures
Versions Notes

Abstract

:
We recently described the genetic antimicrobial resistance and virulence profile of a collection of 279 commensal E. coli of food-producing animal (FPA), pet, wildlife and human origin. Phenotypic antimicrobial resistance (AMR) and the role of commensal E. coli as reservoir of extra-intestinal pathogenic Escherichia coli (ExPEC) virulence-associated genes (VAGs) or as potential ExPEC pathogens were evaluated. The most common phenotypic resistance was to tetracycline (76/279, 27.24%), sulfamethoxazole/trimethoprim (73/279, 26.16%), streptomycin and sulfisoxazole (71/279, 25.45% both) among the overall collection. Poultry and rabbit were the sources mostly associated to AMR, with a significant resistance rate (p > 0.01) to quinolones, streptomycin, sulphonamides, tetracycline and, only for poultry, to ampicillin and chloramphenicol. Finally, rabbit was the source mostly associated to colistin resistance. Different pandemic (ST69/69*, ST95, ST131) and emerging (ST10/ST10*, ST23, ST58, ST117, ST405, ST648) ExPEC sequence types (STs) were identified among the collection, especially in poultry source. Both ST groups carried high number of ExPEC VAGs (pandemic ExPEC STs, mean = 8.92; emerging ExPEC STs, mean = 6.43) and showed phenotypic resistance to different antimicrobials (pandemic ExPEC STs, mean = 2.23; emerging ExPEC STs, mean = 2.43), suggesting their role as potential ExPEC pathogens. Variable phenotypic resistance and ExPEC VAG distribution was also observed in uncommon ExPEC lineages, suggesting commensal flora as a potential reservoir of virulence (mean = 3.80) and antimicrobial resistance (mean = 1.69) determinants.

1. Introduction

Escherichia coli represents a commensal colonizer of human and animal gastrointestinal microbiota [1] and it is the most frequently isolated Gram-negative pathogen impacting human health [2]. E. coli is a worrisome public health threat due to its outstanding variability in pathotypes (enteropathogenic E. coli—EPEC, enterohaemorrhagic E. coli—EHEC, enterotoxigenic E. coli—ETEC, enteroaggregative E. coli—EAEC, enteroinvasive E. coli—EIEC, diffusely adherent E. coli—DAEC, uropathogenic E. coli—UPEC, meningitis associated E. coli—MNEC, sepsis associated E. coli—SEPEC), multiple infection sites (intestinal or extraintestinal), clinical symptomatology [1] and concerning antimicrobial resistance (AMR) profile, especially to carbapenemase and extended-spectrum beta-lactams [3,4]. The different E. coli pathotypes have been associated to specific phylogroups, each one showing distinct phylogenetic relatedness [5,6].
In particular, ExPEC represents one of the most common causes of bloodstream infections and community/hospital associated urinary tract infections worldwide [7,8]). They are also responsible for other extraintestinal diseases, being an important cause of neonatal meningitis [1]. Multiresistant ExPEC strains constitute ongoing healthcare concern and are associated with an increase in infection severity, treatment failure, hospitalisations and mortality, with growing costs for health care [2]. ExPEC mostly belonged to phylogroup B2 and, to a lesser extent, to phylogroup D [9,10].
It has been hypothesised that ExPEC are opportunistic pathogens. They may occupy a niche in human and animal intestinal microflora, showing their virulence potential while colonising extraintestinal sites [11]. Discrimination with molecular epidemiological approaches between potential ExPEC and commensals could be challenging to afford, because of the share of large genomic fractions and different VAGs, involved in common fitness function [12,13]. Therefore, ExPEC could be hidden among commensal flora, which could also represent a reservoir of virulence genes for these pathogens.
Commensal E. coli is also recognised as an AMR barometer among Gram-negative bacteria, due to its ubiquity and genomic plasticity. Indeed, E. coli represents the prevailing organism able to grow in aerobic conditions in the gastrointestinal microbiota of warm-blooded animals [14] and is also an environmentally adapted bacterium [15]. Its genomic plasticity allows constant and efficient exchange of genomic fractions, including genes conferring resistance to antimicrobials, with other enteric bacteria and the environment. Hence, AMR data acquired from E. coli indicator is considered representative of the overall bacterial population [16].
AMR epidemiological surveillance is an important tool, signaling changes in current bacterial AMR trends [17]. The data gathered from epidemiological surveillance allow the implementation of preventive and control strategies, including suitable antimicrobial stewardship programs, therapeutic guidelines and infection control policies [18].
E. coli indicator is routinely used in European AMR monitoring to oversee AMR in FPAs and related food since 2014 [16]. Attention has been focused on specific FPA categories (poultry, swine, turkey and bovine), because of their high population size and important meat demand (http://www.fao.org/faostat/en/, last accessed 30 June 2020). Additionally, overall antimicrobial consumption in the European Union (EU) is mostly associated to livestock, hence exposed to high antimicrobial selective pressure [19].
FPAs could play a primary role in AMR diffusion through direct/indirect contact with humans, related food and environmental (land and water) manure contamination [20,21,22,23,24]. Furthermore, they could represent an important AMR reservoir, where foodborne pathogens may obtain new antimicrobial resistance genes and develop novel resistance patterns [25].
In the last few decades AMR sources other than livestock have proved to be particularly concerning.
AMR bacteria and genes have been identified in human related environments (livestock, companion animals, animal and non-animal origin food) as well as remote ecosystems (wildlife), where antimicrobial selective pressure is supposed to be absent [23,26,27,28,29,30]. These data suggest the complexity of AMR dynamics, influenced by the thriving antimicrobial resistance gene (ARG) bacterial trade and the interconnection between different ecosystems.
AMR arises from selective pressure induced by antimicrobial treatments. Resistant bacteria and genes can be gathered, maintained and transmitted through horizontal gene transfer (HGT) within any environmental niche, determining the diffusion of novel AMR profiles in the overall bacterial population [31]. The complex interconnection between different ecosystems, sharing common habitats and water sources, is responsible for the wide AMR geographical distribution. Water and soil faecal contamination plays a primary role in AMR spread, establishing a link between various environments [32]. Further investigations are needed to achieve a better comprehension of AMR dynamics, focusing on the identification of potential AMR sources and transmission paths that are posing a risk to human health.
In our recent paper [33] we described the population structure, ARG and VAG carriage in a collection of 279 E. coli indicator of animal (livestock, pets, wild animals), food (vegetable and animal origin products) and human origin collected in Italy. E. coli strains were grouped in 12 sources (dairy, beef, wild boar, rabbit, poultry, swine, companion animal, vegetable, fishery, mollusc, wild animal and human), according to their origin. Concerning genetic AMR profile have been identified, including to highest priority critically important antimicrobials (HP-CIAs). Different typical ExPEC VAGs and pandemic and emerging ExPEC STs were observed among the overall collection.
Considering these previous findings, the current study aims to investigate (i) phenotypic antimicrobial resistance profile of the collection and concordance with genetic AMR profile previously identified, in order to establish the potential AMR risk associated with animals (livestock, companion animals, wildlife), food and human, (ii) ExPEC virulence potential of commensal E. coli, association with concerning AMR profile and their role of ExPEC VAG reservoir, (iii) phylogroup distribution, considering the revisited phylotyping method proposed by Clermont et al., (2013) [6], for further epidemiological evaluations.

2. Results

2.1. Occurrence of E. coli Strains in Animal, Food and Human Samples

E. coli was identified in 169/415 (40.72%) samples analysed. In particular, E. coli was detected in all companion animal (12/12) and swine (6/6) samples and in 25/27 (92.59%) human faeces analysed (Table 1). High E. coli occurrence was observed in poultry (25/33, 75.76%), rabbit (10/14, 71.43%), wild boar (22/31, 70.97%) and beef (21/34, 61.76%) samples. Fishery and vegetables were the sources with the lowest occurrence of E. coli, identified in only 24/94 (25.53%) and 24/164 (14.63%) samples, respectively. It was not possible to establish E. coli occurrence in mollusc, dairy and wild animal sources, since related strains were provided from the Experimental Institute for Zooprophylaxis in Lombardy and Emilia Romagna.

2.2. Antimicrobial Resistance Phenotypes

Among the 279 E. coli strains, 107 (38.35%) showed resistance to at least one antimicrobial and 79 (28.32%) were MDR (Table 2). The most common phenotypic resistances were to tetracycline (76/279, 27.24%), sulfamethoxazole/trimethoprim (73/279, 26.16%), streptomycin and sulfisoxazole (71/279, 25.45% both), ampicillin (66/279, 23.66%), followed by nalidixic acid (48/279, 17.20%), enrofloxacin (44/279, 15.77%), chloramphenicol (23/279, 8.24%) and gentamicin (11/279, 3.94%). Lastly, 3GC (ceftiofur and ceftazidime) and colistin resistances were observed in 6/279 (2.15% each) strains. The most common antimicrobials implicated in MDR were tetracycline, sulphonamides, streptomycin and ampicillin.
Four out of 279 strains (1.43%) were designated as ESBL producers. Notably, all ESBL producers were MDR, including to other HP-CIAs (nalidixic acid, enrofloxacin, colistin).
Considering the different sources investigated, the mean number of resistance was: rabbit, 6; poultry, 4; dairy, companion animal, swine and human, 2 each; beef, fishery and wild animal, 1 each. The lowest resistance (mean ≤ 1) was observed in mollusc, vegetable and wild boar sources, with 4, 2 and 1 resistant strains, respectively. Notably, rabbit and poultry were the sources displaying the most extensive AMR, with 23/23 (100%) and 23/25 (92%) strains resistant to at least one antimicrobial agent, respectively. Most of the strains (22/23, 95.7% in rabbit; 16/25, 64% in poultry) were MDR from 3 up to 8 different antimicrobial classes. Furthermore, rabbit was the niche mainly associated to colistin resistance (3/6).
3-GC resistance was mainly reported in dairy (2) and beef (2) strains. Interestingly, the dairy source carried the highest number of ESBL producers (2/4).
Wild animal, vegetable, fishery and companion animal niches generally displayed low resistance when compared to the other sources. MDR strains were detected with a lesser extent (companion animal, 3/12, 25%; wild animal, 6/25, 24%; fishery, 3/24, 12.5%; vegetable, 2/24, 8.33%), though resistance to HP-CIAs (Qs, 3-GCs) was observed in 6 strains (5 wild animals, 1 companion animal).
Notably, one ESBL producer was associated to wildlife (wild animal and wild boar sources).
Wild boars and molluscs revealed the lowest antimicrobial resistance among all the collection. Five strains showed resistance to at least one antimicrobial agent (1/22, 4.54%, wild boar; 4/25, 16%, mollusc), of which 4 (wild boar, 1; mollusc, 3) were MDR. Resistance to HP-CIAs was detected in 1 wild board and in 1 mollusc FQ resistant strains. The detailed phenotypic antimicrobial profile of the overall collection is shown in Figure 1.
The phenotypic AMR pattern identified among the collection was generally concordant with the genetic AMR profile previously identified (Table 3). Discrepancies only occurred with nalidixic acid/enrofloxacin (6 strains), chloramphenicol (1 strain) and colistin (1 strain) antimicrobials.

2.3. Phylogroups

Seven different phylogroups (A, B1, B2, C, D, E and F) were identified in the collection, as described in Table 4. Phylogroup B1 was the most common (127/279, 45.52%), followed by C (43/279, 15.41%), A (38/279, 13.62%), E (25/279, 8.96%), B2 (20/279, 7.17%), F (13/279, 4.66%) and D (9/279, 3.23%).
Four strains were designated as unknown since it was not possible to assign them to any phylogroup.
One specific phylotype was usually representative in each source despite variability in phylogroup distribution and abundance. Phylogroup B1 was mainly associated to beef (22/24, 91.67%), rabbit (20/23, 86.96%), wild animal (18/25, 72%), mollusc (15/25, 60%), vegetable (14/24, 58.3%), companion animal (6/12, 50%) and dairy (9/25, 36%) sources. Phylogroup A was most common in fishery (10/24, 41.67%), meanwhile phylogroup C was mostly represented in poultry (8/25, 32%). Swine strains were equally associated to phylogroup A (8/25, 32%), B1 (8/25, 32%) and C (8/25, 32%). Interestingly, “pathogen” phylogroups were mainly associated to human (F, 7; E, 3; B2/D, 2 both) and wild boar (E, 11; B2, 6) strains.
Notably, strains showing resistance to the highest number of antimicrobials (≥5) were mainly associated to “non pathogen” phylogroups (B1, 33; A–C, 6 both) and with a lesser extent to the “pathogen” ones (F, 6; E, 3; D, 2).
Comparison between molecular and in silico phylotyping data showed important concordance. Only 4 strains, properly assigned to a phylotype with ARIBA [34], were designated as unknown with the quadruplex PCR method (Figure 2).

2.4. Expec Virulence Potential

Thirty-one typical ExPEC VAGs were identified among the overall collection (Table 6). All strains carried between 1 and 17 VAGs (Figure 2), assembled in variable virulence profiles (Figure 1). Among the collection, VAG occurrence was as follow: fimH (270/279; 96.77%), iss (162/279; 58.06%), traT (125/279; 44.8%), sitA (100/279; 35.84%), fyuA (55/279; 19.71%), irp2/iucD (53/279; 18.90%), iroN/malX (50/279; 17.92%), ompT (45/279; 16.13%), iutA (41/279; 14.70%), cvaC (37/279; 13.26%), iha/kpsMT (23/279; 8.24%), ireA (21/279; 7.52%), vat (16/279; 5.73%), usp (15/279; 5.38%), papC/tsh/ibeA/pic (11/279; 3.94%), neuC (10/279; 3.58%), papG (7/279; 2.51%), sat (6/279; 2.15%), cdtB (5/279; 1.79%), bmaE/gimB (4/279; 1.43%), hlyE (3/279; 1.08%), sfaS/cnf1 (1/279; 0.36%).
Thirteen strains belonged to pandemic ExPEC STs, namely 10 ST69/69* (dairy, 1; wild boar, 2; poultry, 1; mollusc, 3; human, 3), 2 ST95 (human) and 1 ST131 (poultry). As expected, pandemic ExPEC STs showed the highest number of VAGs (mean = 8.92) and AMR phenotypes (mean = 2.23).
Interestingly, 7/13 (53.85%) pandemic ExPEC STs were associated to AMR, of which 6 were multidrug-resistant. Notably, ST131 showed the highest number of resistance to 9 antimicrobials. HP-CIA resistance was observed in 3 FQ resistant ST69 strains of poultry, mollusc and human origin. All pandemic ExPEC STs belonged to typical ExPEC phylogroups (B2, 3; D, 5; E, 5).
Thirty-seven emerging ExPEC STs were detected in different sources, namely 19 ST10/ST10* (dairy, 3; beef, 1; vegetable, 1; companion animal, 1; swine, 7; poultry, 2; mollusc, 2; human, 2), 6 ST23 (poultry, 3; mollusc, 1; wild animal, 2), 4 ST58 (dairy, 1; beef, 2; companion animal, 1), 6 ST117 (1 wild boar; 5, poultry), 1 ST405 (human) and 1 ST648 (wild animal). Most of the strains were associated to commensal phylotypes (A, 4; B1, 4; C, 21), with the exception of ST117 (5, F; 1F), ST405 (E) and ST648 (F) belonging to typical ExPEC phylogroups. VAG carriage and AMR phenotype mean were 6.43 and 2.43 respectively. MDR was observed in 14/37 (37.84%) strains. Notably 10 Q and 1 colistin resistant strains were identified.
The remaining 229 strains carried the lowest AMR phenotypes (mean = 1.69) and VAG number (mean = 3.80) among the 3 groups identified. Variable phylogroup distribution (A = 34; B1 = 123; B2 = 17; C = 22; D = 4; E = 18; F = 7) was observed. Fifty-three out of 59 MDR strains detected belonged to commensal phylogroups. Uncommon ExPEC STs showed the highest HP-CIA resistance, with 35 Q, 4 colistin resistant strains and 4 ESBL producers.
Pandemic, emerging and uncommon ExPEC STs showed variable distribution of genes coding for the 5 functional categories involved in ExPEC pathogenesis (Table 5).
Interestingly, eleven strains carried genes coding for all 5 functional categories. Threestrains belonged to typical ExPEC STs (ST131, 1; ST95, 2) and phylotypes (B2, 9; D, 1; E, 1), meanwhile 8 strains were associated to uncommon ExPEC STs (ST141, ST680, ST1170, ST80, ST706*, ST420*, ST136, ST372) but typical ExPEC phylogroups (B2, 6; D, 1; F, 1). MDR was observed in ST131 (1), ST95 (2), ST706* (1) strains.
Complete phenotypic AMR and genetic virulence profile, STs and phylogroup distribution among the collection are reported in Supplementary Materials (Tables S2–S4).

2.5. Statistical Results

Statistical analysis identified poultry and rabbit as the most important AMR sources among the collection. In these niches, the occurrence of resistance to at least one molecule and of MDR was statistically significant (p-value < 0.01) if compared to the other sources. Poultry and rabbit strains exhibited significantly higher (p-value < 0.01) resistant rates to Qs, streptomycin, sulfisoxazole, sulfamethoxazole/trimethoprim and tetracyline than the remaining collection. Furthermore, significant ampicillin and chloramphenicol resistance was observed in poultry (p-value < 0.01).
Considering the other sources, no significant resistance was detected, except for tetracycline resistance in swine strains (p < 0.01).
Important differences in phylogroup richness and eveness were identified among the sources (Table 4), with the highest diversities in human (H′ = 1.814; D′ = 0.850) and poultry (H′ = 1.691; D′ = 0.817) and the lowest in rabbit (H′ = 0.470; D′ = 0.245) and beef (H′ = 0.345; D′ = 0.163). Association between specific phylogroups and phenotypic antimicrobial patterns was investigated. Statistical analysis identified a significant association for ampicillin resistance—E/F phylotype and sulfisoxazole resistance—F phylotype (p < 0.05 both).

3. Discussion

AMR is a complex and ever-changing phenomenon, whose dynamics are still not completely understood. Epidemiological surveillance is an essential weapon against AMR threat, allowing to define key point sources where AMR could develop, evolve and spread. AMR monitoring in different environments represents a fundamental cornerstone for epidemiological evaluation and preventive/control measure implementation [18]. Indeed, multiple sources (animals, agriculture, human, food) have been investigated as a potential reservoir of AMR in the last decades [23,26,28,29,30,35].
Our previous study aimed to provide an overview of AMR genetic profile currently circulating in E. coli indicator in Italy. E. coli from different origins have been considered, including environments playing a well documented (FPAs, food, pets, human) or emerging (wildlife) role in AMR dynamics. Greater emphasis was given to livestock due to their important epidemiological role in AMR spread through human contact [36,37], environmental manure contamination [21,38] and related foodstuff [23,39,40,41,42]. FPAs may play an important role in the development of new antimicrobial resistance patterns, increasing the AMR gene pool available for foodborne pathogens [25].
Particular attention was given to HP-CIAs (including extended-spectrum cephalosporins, quinolones and colistin), recognised as the last treatment option for serious human infections with a possible FPA/food origin [7].
Antimicrobial susceptibility testing results were interpreted using ECOFF in order to perform an epidemiological evaluation about phenotypic AMR circulating in E. coli indicator of different origins.
In the present paper, wild type (WT) strains, without an acquired phenotypically detectable resistance mechanism, were defined as “susceptible”. Non-WT strains, related to acquired or mutational resistance to antimicrobials, were addressed as “resistant”.
Unfortunately, ECOFFs were not available for all the molecules tested, hence clinical breakpoints were adopted. Clinical breakpoints represented predictors of clinical success of antimicrobial treatments, but could not afford evaluation related to emergence and evolution in bacterial resistance profiles [43,44]. Nevertheless, they could provide worthwhile information about consolidated resistance patterns in the bacterial population. Changing in consolidated clinical resistant profiles could indirectly suggest the emergence of new resistance patterns among bacteria.

3.1. Tetracycline, Sulfonamide, Streptomycin and Ampicillin Resistance Is the Most Common among the Overall Collection

The most common phenotypic resistance was to tetracycline, sulfonamides (sulfisoxazole and sulfamethoxazole/trimethoprim), streptomycin and ampicillin among FPAs/food, companion animal, wildlife and humans.
The wide diffusion of resistance to tetracycline, sulfonamides and ampicillin is probably attributable to antimicrobial usage in different human and animal sectors in Europe. Sulfonamides and tetracycline represent the most common treatment in FPAs (swine, poultry, beef, rabbit and dairy), as well as ampicillin in FPAs (swine, poultry, beef, rabbit and dairy), companion animals and human [45]. In particular, the zootechnical system is associated to the highest antimicrobial consumption in Europe [19], turning livestock into an important reservoir of AMR bacteria and genes. Indeed, the high frequency of resistance to sulfonamide-ampicillin-tetracycline in swine, poultry, beef, rabbit and dairy strains was generally congruent with AMR profiles previously reported in Europe [16,46,47,48,49,50,51,52]. The same resistance pattern (sulfonamide-ampicillin-tetracycline resistance) was also observed in the strictly connected human environment, linked in turn with companion animals. High ampicillin resistance rate probably derived from the wide use of these molecules in human and pets [45], whereas sulfonamide and tetracycline resistance presumably originated from agricultural settings. Sulfonamide-ampicillin-tetracycline resistance has also been previously observed in human and pets strains [53,54].
Interestingly, streptomycin resistance was the third-highest phenotypic resistance (together with sulfisoxazole) (24.45%) among all the collection and it was identified in at least one strain of each source. These data are generally consistent with previous studies in animal and human [32,47,55,56,57,58,59]. Exceptions are represented by pig, poultry and beef, usually showing low level of aminoglycoside resistance. Streptomycin resistance is inconsistent with the seldom use of aminoglycosides in both human and animal therapeutic treatment [45,60], which account for low gentamycin resistance among the collection (3.94%) instead. So far, the evidence leads to the hypothesis of streptomycin resistance wide diffusion being referable to ARGs co-selection mechanisms, mediated by multiresistant plasmids [61] and transposons [62,63]. Additionally, streptomycin is one of the oldest antimicrobials discovered in 1943 [64] and its past use could explain the active resistance circulation detected nowadays.

3.2. The Potential Role of Aquaculture, Vegetable and Wildlife as AMR Source Sentinels

Generally, aquaculture, vegetable and wildlife (wild animal and wild boar) are sources slightly connected to AMR. In aquaculture (mollusc and fishery) antimicrobial treatments are rare [65], meanwhile vegetable and wildlife are not directly exposed to antimicrobials. The consistency of AMR profile identified in these niches with those observed in the remaining collection (livestock, human, companion animals) suggested an important AMR environmental pollution from settings highly associated to antimicrobial use, in line with findings reported by Giacometti et al., (2021) [66]. Water and environment faecal contamination appear the most important path allowing ARGs dissemination in ecosystems where antimicrobial use is rare or absent [32]. Irrigational water/manure [57,67] and sewage/runoff from land [68,69,70] have already been implicated in AMR emergence in aquaculture and vegetable products, respectively. Likewise, remoteness of area, zootechnical activity and human proximity represented the most important parameters influencing AMR occurrence in wildlife [71,72].
Interestingly, resistance to sulfonamide-ampicillin-tetracycline in aquaculture [56,73], vegetable [56,57,74] and wildlife [75,76] has already been reported in different countries suggesting the wide extension of AMR diffusion problem in outlying environments.
Considering all these assumptions, aquaculture, wildlife and vegetable could represent important AMR source-sentinel, giving useful information about specific AMR spread and the degree of AMR environmental contamination. However, the low E. coli indicator occurrence in fishery, mollusc and vegetable sources could suggest wildlife as the most promising and concrete AMR spread indicator. Additionally, low E. coli occurrence indicates fishery, mollusc and vegetable as low risk food-sources associated to transmission of potential pathogenic E. coli to human, when compared to swine, poultry and rabbit niches, where E. coli was frequently identified.

3.3. Phenotypic Q Resistance Is the Most Common Among HP-Cias

HP-CIA resistance represents a worrisome event, associated with significant morbility/mortality and treatment alternatives reduction, and has required considerable control during food production [7]. Qs are included in the HP-CIA list and represent the last treatment options in serious Salmonella spp. and E. coli infections [7]. Notably, Q resistance was mainly observed in FPAs, with a lesser extent in human and wildlife strains, and not detected in companion animals and vegetable.
Considering FPAs, the most common HP-CIA resistance was to Qs (21.05%, nalidixic acid; 18.71%, enrofloxacin), identified in all the livestock categories but in fishery source. The highest FQ resistance was observed in rabbit (nalidixic acid, 15/25; enrofloxacin, 13/25) and poultry (nalidixic acid, 13/25; enrofloxacin, 11/25) strains, presumably because of the common use of these antimicrobials in related breeding systems [65]. Resistance rates in these niches were significantly higher (p < 0.01) than those detected in the other FPA categories, where Q detection was notably lower. Our data are in accordance with previous studies reporting generally high Q resistance in poultry and rabbit in contrast to other livestock [47,48,49,52,77,78,79].
As for livestock, the highest HP-CIA resistance observed in human and wildlife was to Qs (24% and 15%, respectively). Notably, our data showed considerable higher Q resistance in human, as well as in wildlife, than those reported in previous studies [53,75,76,80,81,82,83].
These findings are particularly concerning, considering the increasing trend of Q resistance in human clinical E. coli in Europe, with Italy as the country showing the highest resistance [8].

3.4. Poultry and Rabbit Are the Sources Mainly Associated to AMR

AMR was mainly associated to livestock, with poultry and rabbit as the sources showing significant higher resistant rate (p < 0.01 or p < 0.05) to most molecules (Qs, streptomycin, sulfisoxazole, sulfametoxazole/trimethoprim, ampicillin and chloramphenicol in poultry), when compared to other niches.
These data are not surprising considering poultry intensive breeding systems are characterised by notably higher population density compared to the other FPAs categories [84]. Population size is strictly interconnected with antimicrobial usage, often characterised by collective treatments, and therefore influencing antimicrobial selective pressure [85]. Indeed, poultry have been frequently associated with generally higher AMR than other livestock in European countries [86]. These considerations could be extended also to rabbit farming.
Poultry AMR profile was generally in accordance with previous studies reporting tetracycline, Qs (nalidixic acid and enrofloxacin), ampicillin, sulfamethoxazole, trimethoprim resistance as frequently detected, meanwhile those to 3GCs and colistin were overall low. Contrary to our results, resistance to chloramphenicol was variable, meanwhile aminoglycoside resistance was generally rare [48,49,50,78,79]. Our findings in rabbit were in accordance with the seldom AMR data available in rabbit breeding system, reporting resistance to tetracyclines, ampicillin, aminoglycosides (streptomycin, gentamicin), Qs (ciprofloxacin, nalidixic acid) and sulfonamides (sulfadiazine, sulfamethoxazole/trimethoprim) as common phenotypic AMR profile [47,87,88,89]. Interestingly, rabbit was the source mostly associated to colistin resistance. Recent studies described colistin resistance in rabbit farms [88], also in Italy, where colistin treatment is reported as common [90]. In our knowledge, this is the first report identifying phenotypic colistin resistance associated to rabbit meat, highlighting the possible involvement of rabbit breeding system in colistin resistance diffusion through the food-chain.
Interestingly, poultry and rabbit resistant strains were mostly isolated from food products of animal origin (poultry meat, 13; rabbit meat, 10) and with a lesser extent from animals (poultry, faeces, 2; rabbit, intestine, 5). The low number of samples analysed could not afford a proper evaluation of AMR epidemiological risk associated to animal and related foodstuff. However, our data could suggest an important involvement of the food chain in AMR transmission, presumably from zootechnical settings up to consumers. This hypothesis is supported by different studies, suggesting poultry as a feasible origin of AMR clinical E. coli identified in human [91,92,93,94]. The scarce European data on AMR in rabbit breeding and foodstuff [47,87,88,89] encouraged further investigations to better elucidate the role of rabbit in AMR transmission to human.

3.5. Phenotypic Pattern Is Generally Concordant with Genotypic AMR Profile

The phenotypic AMR pattern identified among the collection was generally concordant with the genetic AMR profile previously identified (Table 5). Discrepancies only occurred with Qs (6 strains), chloramphenicol (1 strain) and colistin (1 strain) antimicrobials. Qs resistant strains carrying qnr genes generally provided low level of resistance [95], but in our collection qnrS1 gene (3) did not present with the expected phenotype. Their importance is mostly associated to the selection of specific chromosomal mutations, favoring the emergence of strain with higher FQ resistance [96].
FQ inefficacy is frequently associated to mutations in the quinolone resistance-determining region (QRDR: gyrA/gyrB and parC/parE), coding for the drug targeting enzymes DNA gyrase (GyrA/GyrB) and topoisomerase IV (ParC/ParE) subunits.
Generally, single or combined mutations in gyrA, parC, parE were associated with the phenotypic resistance profile. Interestingly, in some cases, substitutions Ala-56 → Thr in the ParC (ATC → ACC in parC) and Ile-355 → Thr in the ParE (ATC → ACC in parE) proteins, previously associated to resistant E. coli [97,98], were identified in nalidixic acid and enrofloxacin susceptible strains.
Chloramphenicol resistance was usually explained by the identification of genes coding for chloramphenicol acetyltransferases (cat) or specific chloramphenicol exporters (cml) [99]. Interestingly, cat and/or cml genes were not identified in one chloramphenicol resistant strain. Carriage of mdfA and acrA-acrB-TolC genes, coding for the aspecific multidrug transporter MdfA [100] and AcrAB-TolC [101], could explain the observed resistant phenotype.
In the present study, colistin resistance was associated to plasmid mediated mcr1 gene, coding for a phosphoethanolamine transferase. This enzyme is responsible for a cationic modification of the LPS of the bacterial outer cell membrane, target of polymyxin antimicrobials [102].
One colistin susceptible strain carried Val-161 → Gly mutation in PmrB, associated to both polymyxin resistance and susceptibility [103,104,105]. pmrB, together with pmrA gene, coded for the PmrAB two-component system activation, typically involved in bacterial survival against cellular mediated immune response [106]. Mutations in these genes could potentially lead to resistance, changing lipopolysaccharide charge and reducing polymyxin attachment to the external surface of Gram-negative bacteria [107].

3.6. Important Phylogroup Variability Occurred among the Different Sources

Clermont revisited phylotyping method [6,108] identified up to 7 different phylotypes among the collection, with important variability among FPAs, pets, human and wildlife.
As expected, B1 (127/279; 45.52%), C (43/279; 15.41%) and A (38/279; 13.62%) were the most common phylogroups among all the collection and also in most sources. These phylotypes are considered “generalists” of multiple hosts and are commonly identified in the commensal population [109].
Notably, typical ExPEC phylotypes were mainly observed in wild boar (17/23, 73.91%; 11, E; 6, B2) and human (14/25, 56%; 7, F; 3, E; 2, B2; 2, D) sources. Our findings are generally in accordance with previous studies describing wild boar strains mostly belonging to B2 and D (including derivatives) [110,111,112] and D (including derivatives) as the second most common phylogroup in human commensal isolates worldwide [6].
Important variability in phylogroup distribution has been observed among sources (table), especially in human (H′ = 1.814; D = 0.850) and some FPA categories (poultry H′ = 1.691, D = 0,817; dairy: H′ = 1.485, D = 0.763; fishery: H′ = 1.457, D = 0.764) showing the highest phylogroup richness.
A proper comparison between phylogroup distribution in our collection and those reported in previous studies is hard to perform, due to the phylotyping method used and the multiple factors influencing phylogroup appearance.
Phylogroup determination is mostly performed according to Clermont scheme, targeting two genes (chuA, yiaA) and a DNA fragment (TSpE4.C2) [5]. This method allows to identify only the four main phylogroup A, B1, B2, D. The lack of proper C, E, F identification prevents detailed epidemiological evaluation of E. coli population genetics in different environments. Furthermore, important information regarding AMR and virulence potential associated to these minor phylogroups are lacking.
Differences in phylogroup distribution are mainly ascribable to geographical (location and climate) and host (diet, gut morphology, body mass) factors [113], explaining variable animals and human phylogroup identification in different studies [109,110,113,114,115,116,117,118,119,120].
Generally, human commensal E. coli are worldwide mostly associated to phylogroup A and D [6,114], meanwhile FPAs mostly belonged to phylogroup A and B1 [109,110]. Differences in wild animal phylotype distribution have been observed, with D and B1 as the most common in birds, A in wild rabbits and D and B2 in wild boars [110].

3.7. Clermont Quadruplex PCR Is a Valid Alternative to In Silico Phylotyping Technique

Molecular and in silico phylotyping data were compared, considering that C, E, F phylotypes, identified with the revisited Clermont scheme [108], would be included in the four main phylogroups (A, B1, B2, D) in the original triplex phylogroup assignment method [5].
Important concordance between results obtained with the two phylotyping techniques was observed. Only exceptions were represented by 4 strains designated as unknown with the quadruplex PCR method and properly assigned to a phylotype with ARIBA.
These findings suggest Clermont quadruplex PCR as a suitable alternative for phylogroup determination to in silico phylotyping. Indeed, quadruplex PCR implementation had contained costs, is time sparing and does not need specific and expensive technology for its implementation.

3.8. Commensal E. coli Conceal Potential Multiresistant Pandemic and Emerging Expec Pathogens

ExPEC are responsible for a majority of human extraintestinal infections worldwide. Actual ExPEC pathogenic potential is interconnected with VAG number [121], as with the expression of specific VFs involved in ExPEC pathogenesis. Indeed, studies in animal models of extraintestinal infection suggest VF profile as an in vivo virulence predictor [122,123].
Successful ExPEC pathogens must perform a range of functions during pathogenesis, namely adhesion/colonisation, host defence evasion, multiplication, tissue damage and diffusion [1]. Specific VFs are essential in each phase of pathogenesis and are generally divided in 5 functional categories [11]: (i) adhesins, important adherent factors promoting host cell contact, hence adhesion and colonisation [124]; (ii) invasins, mediating cell invasion into the host tissues [124]; (iii) iron acquisition systems, allowing iron uptake in low iron conditions (i.e., in host fluids and tissue) [125]; (iv) toxins, responsible of tissue lesion and promoter of local bacteria diffusion, cytotoxicity and insensitivity to neutrophils [126]; (v) protectins, structural components of the bacterial outer membrane involved in host defense evasion (including resistance to innate immunity and serum survival increase) [127].
Notably, 50 pandemic and emerging ExPEC lineages were identified in our collection. Both groups showed high number of typical ExPEC VAGs (pandemic ExPEC mean = 8.92; emerging ExPEC mean = 6.43), including those coding for adhesins (fimH and papG), invasins (ibeA and gimB), iron acquisition systems (ireA, iroN, fyuA, irp2, iucD, iutA, sitA), serum survival protectins (iss and traT) and toxins (usp, vat, pic, sat, hlyE).
ST131, ST95 and ST69/69* are members of the predominant clonal ExPEC group, with ST131 as the most common ExPEC lineage isolated worldwide [128].
B2 ST131 (1) and B2 ST95 (2) strains carried an outstanding number of VAGs (17 and 12 respectively), coding for VFs of all the functional categories involved in ExPEC pathogenesis.
Their genetic virulence pattern showed similarities to clinical ST131 and ST95 prototypic virulence profile [129], suggesting the virulence potential of these strains.
Notably, ST95 strains carried papG_II, a specific papG allele coding for the P pilus tip adhesin, responsible for UPEC adhesion in the urinary tract [130,131]. papG_II adhesin recognises glycolipid receptors located in the human kidney [132] and is frequently associated with pyelonephritis [133]. Also toxin genes (usp and vat), typically identified in UPEC [35,134,135,136], were observed in both ST131 and ST95 strains.
On the other side, ST69/69* strains showed significant variability in VAGs carriage (between 2 and 12) with rather different genetic traits compared to those commonly described in clinical strains [129]. For example, toxin genes (cdtB, hlyD, cnf1) were not identified, though they are commonly reported in this lineage [137].
The high number of strains carrying typical ExPEC VAGs is particularly worrisome. It has been hypothesised that ExPEC represent facultative pathogens. They may occupy a niche in human and animal intestinal microflora, showing their virulence potential while colonising extraintestinal sites [11]. Discrimination with molecular epidemiological approaches between potential ExPEC and commensals could be challenging to afford, because of the share of large genomic fractions and different VAGs, involved in common fitness function [11].
ExPEC could be hidden among commensal flora, arising and showing their virulence with favourable conditions occurrence. Further studies are needed to better investigate commensal E. coli virulence traits and to elucidate their role in ExPEC infections. Identification of commensal E. coli virulence profile could be useful to determine sources potentially associated to ExPEC transmission.
Additionally, the highest phenotypic AMR was observed in both pandemic (mean = 2.23; MDR = 46.15%) and emerging (mean = 2.43; MDR = 37.84%) STs. A possible reason could be the potential co-carriage of virulence and antimicrobial determinants on the same genetic platform (plasmids, transposons, integrons), co-mobilized under antimicrobial selective pressure [138]. A typical example is represented by F plasmid replicons that are considered the E. coli major virulence-associated plasmids, as well as multiple ARG carriers [139,140,141]. Interestingly, HP-CIA resistance was mainly described in uncommon ExPEC lineages, though this AMR profile is more likely identified in ExPEC STs [142].

3.9. E. coli Indicator Are Expec VAG Reservoir

Although uncommon ExPEC lineages group was associated to the lowest number of VAGs (mean = 3.80), the virulence profile showed extreme variability in VAG number (from 1 up to 15 VAGs) and presence of VFs related to different functional categories (from 1 up to 5). In our findings, uncommon ExPEC lineages seem unlikely associated to virulence potential. However, ExPEC genetic traits in the commensal population may be acquired from potential ExPEC pathogens, coexisting in the gastrointestinal tract. Therefore, commensal E. coli could represent a critical VAG reservoir, increasing virulent armament of potential ExPEC pathogen or allowing the acquisition of virulence traits in traditionally avirulent bacteria. A notable exception is represented by 8 strains, carrying genes coding for all functional category members. Phylogroup association (6 B2, 1 D, 1 F) and virulence profile were similar to those observed in B2 ST131 and B2 ST95 strains, suggesting that in some cases unconventional ExPEC lineages may display ExPEC potential.

3.10. Pandemic and Emerging Expec STs Mainly Belonged to Poultry Source

Variable ExPEC profile has been observed in various sources, with poultry as the niche mostly associated to ExPEC potential. Indeed 11/25 (44%) poultry (ST117, 5; ST23, 3; ST131, 1; ST69, 1; ST10, 1) strains belonged to emerging or pandemic ExPEC STs and carried high number of VAGs (poultry: mean = 11.73).
Poultry has been suspected to represent an ExPEC animal reservoir for humans, due to the close genetic relationship between human EXPEC and avian pathogenic E. coli (APEC). Indeed, ExPEC and APEC may share genome content, phylogeny and genetic virulence profile [143]. Additionally, experimental studies suggest APEC as potential pathogen in mammals, as well as human-derived ExPEC showing their virulence in avian animal models [144]. Concerningly, different studies identified potential ExPEC STs in poultry meat, suggesting its role as a vector of potential ExPEC pathogen to human [145,146].
Faeces contamination of carcasses during slaughtering procedures seems to be the most feasible path of ExPEC diffusion through the food chain [147,148] as APEC are assumed to coexist with commensal microflora in gastrointestinal tract of asyntomatic animals [149].

4. Materials and Methods

4.1. Samples Collection

This study analysed a collection of 279 commensal E. coli, arranged in 12 sources: poultry (n = 25), swine (n = 25), vegetable (n = 24), fishery (n = 24), mollusc (n = 25), wild animal (n = 25) and human (n = 25).
For bacterial strain gathering, a total of 433 samples of food, animal and human origin (beef, 33; wild boar, 31; vegetable, 164; fishery, 94; companion animal, 12; swine, 6; poultry, 33; rabbit, 14; human, 27) were assembled and processed in the period between 2010 and 2018. After analysis, 169 E. coli were collected (beef, 21; wild boar, 22; vegetable, 24; fishery, 24; companion animal, 12; swine, 6; poultry, 25; rabbit, 10; human, 25).
Additionally, 110 presumptive E. coli isolates were provided by the Department of Veterinary Medical Science—University of Bologna—Service of Food Safety (dairy, 25) and Experimental Institute for Zooprophylaxis in Lombardy and Emilia Romagna (beef, 3; rabbit, 13; swine, 19; mollusc, 25; wild animal, 25) and included in the study.
The strains retain their original name and the respective genome sequences are available in the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/, last accessed 30 Jane 2020). Accession numbers are listed in Supplementary Materials Tables S2–S4.

4.2. Genetic Features of the Collection

In our previous study an overall description of E. coli population structure, genetic virulence and antimicrobial resistance profile and phylogroups was described and herafter summarised. One hundred-eight out of 279 (38.71%) strains carried between 1 and 18 ARGs, with consistent variability in genetic resistance profile among sources. The most common ARGs were tetracycline resistance gene tetA (57/279, 20.43%), sulfonamide resistance gene sul2 (45/279, 16.12%), penicillin resistance gene blaTEM-1b (43/279, 15.41%) and streptomycin resistance genes strA/B (42/279; 15.10%) among the overall collection. In that study, genes and chromosomal point mutations conferring resistance to HP-CIAs were widely detected, including to quinolones (Qs) (qnrS1, qnrB19, gyrA/parC/parE), colistin (mcr1, pmrB) and 3d generation cephalosporins (3GC) (blaCMY-2, ampC). Q resistant determinants were mainly observed in rabbit (17/57; 29.82%) and poultry (13/57; 22.81%) sources, meanwhile colistin ARGs were detected in rabbit (3), swine (3) and dairy (1) strains. 3GC resistance genes were identified in 2 beef and 1 human strains, respectively. Extended-spectrum beta-lactamase (ESBL) genes (blaCTX-M1,15,55) were rarely observed and identified in animal (2 dairy and 1 wild animal strains) and human (1) samples. One hundred and eleven different VAGs were detected among the collection. All strains carried between 1 and 37 VAGs. Virulence profile often included typical ExPEC VAGs. Notably, different pandemic (ST69, ST95, ST131) and emerging (ST10, ST23, ST58, ST117, ST405, ST648) ExPEC lineages were observed, especially in poultry meat strains (ST131, ST69, ST10, ST23, and ST117). Considering Clermont triplex PCR phylotyping method [5], the most common phylogroup identified among the collection was B1 (130; 46.6%), followed by A (81, 29%) D (47; 16.8%) and B2 (21; 7.5%).

4.3. Bacterial Isolation and Molecular Identification

A comprehensive description of the isolation and the molecular identification protocols are available in Massella et al., (2020) [33].
Briefly, lactose fermenting colonies were selected on MacConkey’s (Oxoid, Basington, UK) and Levine’s (Oxoid, Basington, UK) agar plates and were incubated for 18–24 h at 37 ± 1 °C. In food sample processing, an additional enrichment step in EC-Broth (Oxoid, Basington, UK) preceded strain isolation on agar plates. Gram stain and standard biochemical test (indole probe) were used for presumptive E. coli identification. Genomic DNA was extracted using a commercial kit (DNeasy Blood and Tissue Kit, Qiagen, Hilden, Germany), following the manufacturer’s instruction. Species identification was performed using the multiplex PCR designated by Horakova et al., (2008) [150].

4.4. Antimicrobial Resistance

Antimicrobial resistance phenotyping was performed with the Kirby-Bauer disk agar diffusion method in accordance with the Clinical and Laboratory Standards Institute guidelines (CLSI) [151]. The antimicrobial panel was chosen considering the antimicrobial genetic profile of the collection, the importance of antimicrobial classes in the treatment of human infections and the intrinsic resistance of E. coli [152]. Particular attention was focused on HP-CIAs, whose resistance is highly suspected to be linked to food-producing sector [7].
The following antimicrobials were tested: nalidixic acid (30 µg), ampicillin (10 µg) ceftazidime (10 µg), ceftiofur (30 µg), chloramphenicol (30 µg), enrofloxacin (10 µg), gentamicin (10 µg), streptomycin (10 µg), sulfisoxazole (300 µg), tetracycline (30 µg) and trimethoprim-sulfamethoxazole (25 µg).
Results were interpreted referring to the epidemiological cut-off values (ECOFFs) for E. coli proposed by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) (http://www.eucast.org; last accessed 30 June 2020). The EUCAST clinical breakpoints for Enterobacteriaceae [153] or the CLSI clinical breakpoints for Enterobacteriaceae [151,154,155] were considered when cut-off values were not available (Supplementary Materials Table S1).
Colistin MIC was determined by the broth microdilution method using customized Sensititre™ 96-well microtitre plates (Trek Diagnostic Systems, East Grinstead, UK). Following manufacturer instructions, 10 µl of bacterial suspension (0.5 McFarland) were placed in 11 mL Mueller-Hinton Broth cation adjusted (Oxoid, Basington, UK). Fifty microliters of the final suspension were put into all wells of the same strip within 30 min after its preparation. Plates were incubated for 18–24 h at 35 °C. Strains were considered susceptible/resistant considering EUCAST ECOFF (2 mg/L) for E. coli (http://www.eucast.org; last accessed 30 June 2020).
The combination disk test [152] was implemented for the evaluation of ESBL producing E. coli. Briefly, the strains were tested for cephalosporins (ceftazidime and cefuroxime) alone and in combination with clavulanic acid, performing disk agar diffusion method. When an increase ≥5 mm in zone diameter was observed in the presence of clavulanic acid compared with the cephalosporin alone, the strain was considered an ESBL producer.
E. coli ATCC 25,922 was used as control strain for antimicrobial testing and as negative control for the evaluation of ESBL profile.
Strains were considered multidrug resistant (MDR) when showing resistance to three or more antimicrobial classes [156].
Lastly, we evaluated concordance between E. coli phenotypic AMR pattern and genotypic AMR profile, identified in the previous publication [33].

4.5. Phylogrouping

Phylogroup evaluation was performed using the quadruplex PCR and primers described by Clermont et al., (2013) [108], able to identify 7 different phylogroups (A, B1, B2, C, D, E and F) and cryptic clade I. Strains generically designated as cryptic clade members in the quadruplex PCR were then screened to establish the specific cryptic clade (II-III-IV-V), according to Clermont et al., (2011) method [157].
Finally, results were compared to phylogroup distribution obtained with in silico phylotyping technique reported in the previous publication [33].

4.6. Expec Virulence Potential

Evaluation of ExPEC virulence potential considered different epidemiological traits, such as ST and virulence gene profile (both identified in the previous study) and phylogroup belonging (according to revisited Clermont scheme) [6]. The collection was divided into 3 main groups according to ST epidemiological role in ExPEC infections [158]: the first was characterised by pandemic ExPEC STs; the second consisted of emerging ExPEC STs; the third included the remaining lineages not specifically associated to ExPEC pathogens. Attention was focused on B2 and D strains, as these phylogroups have been historically associated to ExPEC pathogens [108].
Genes coding for typical ExPEC virulence factors (VFs) of 5 functional categories (adhesins, invasins, iron acquisition systems, toxins and protectins) [11] were selected from the original virulence profile. Additional genes were also included, such as ones coding for vacuolating autotransporter toxin (vat), genetic island associated with newborn meningitis (gimB) [125] and uropathogenic specific protein (usp) [159] (Table 6).
The resulting ExPEC genetic profile, associated to a specific phylogroup and ST belonging, was then considered in order to establish ExPEC pathogenic potential of the collection. Association of virulence and important AMR traits (MDR, HP-CIA resistance, ESBL profile) was also investigated.

4.7. Statistical Analysis

Differences in the occurrence of resistance to antimicrobial agents in E. coli strains among different sources were tested by Pearson’s chi-square. The same statistical test was used to assert differences in phylogroup distribution among E. coli of different niches and association between specific resistance and phylogroup. A p value < 0.05 was considered statistically significant.
Diversities in phylogroup assignment in terms of number and phylotype among E. coli of different origin were calculated using Shannon–Weaver diversity index (H′), H′ = −∑si = 1Pi ln(Pi), where Pi is the percentage of strains belonging to the i-th phylotype of the total number of strains surveyed [163], and Simpson’s diversity index (D′), D = −∑si = 1Ni (Ni − 1)/N (N − 1), where Ni is the number of strains in the i-th phylotype and N is the total number of strains [163].

5. Conclusions

In conclusion, our study provides an important overview of phenotypic AMR, ExPEC virulence potential (according to the revisited Clermont scheme) and phylogroup distribution in commensal E. coli of animal, food and human origin.
AMR pattern in human, companion animal and most FPA categories reflected general phenotypic resistance trends and antimicrobial stewardship in Europe. Identification of human and animal (livestock and companion animal) AMR profiles in niches with a rare (fishery, mollusc) or absent (vegetable, wild animal, wild boar) direct exposure to antimicrobials suggested widespread environmental AMR pollution. Some sources (wild animal and wild boar) may represent important AMR source sentinel. AMR was mainly associated to FPAs, already suspected to play a major role in AMR diffusion.
The important virulence profile identified among the collection proposed commensal E. coli as ExPEC VAG source, from which potential pathogen may acquire new virulence traits. Identification of different pandemic/emerging ExPEC STs and similarities in virulence profile between commensal and clinical human ExPEC lineages could suggest ExPEC virulence potential of commensal bacteria.
Rabbit and poultry were the most concerning sources, associated to the highest AMR among all the collection and suggested as potential AMR reservoirs. Additionally, different pandemic/emerging ExPEC STs and important virulence profiles were observed in poultry strains, as already described in the literature. In poultry and rabbit sources, identification of concerning AMR and virulence profile in an important number of meat origin strains suggested the food chain as a feasible transmission path of potential multiresistant pathogens to human.
Phylogrouping revealed a complex phylotype distribution, attributed to different host and geographical factors. General concordance was observed between phenotypic and genetic AMR profile, as between molecular and in silico and phylogrouping, suggesting the revisited Clermont method as a reliable and cheaper phylotyping technique alternative.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/antibiotics10040351/s1, Table S1: AMR breakpoint and cut-off considered in the present study; Table S2: phenotypic AMR profile of the collection; Table S3: genetic AMR profile of the collection; Table S4: genetic virulence profile of the collection.

Author Contributions

Conceptualization, E.M., A.S., S.P., F.G.; methodology, E.M., S.P., F.G.; validation: E.M.; formal analysis, E.M.; investigation: E.M., S.P., F.G., A.S.; resources, A.S., S.P., F.G.; data curation, E.M., S.P., F.G., A.S.; writing—original draft preparation, E.M.; review and editing, E.M., A.S., S.P., F.G., C.J.R., S.P.D., P.B., G.M., C.B., L.F., P.M., L.B., S.R., F.S.; visualization: E.M.; supervision, A.S., S.P., F.G.; project administration, A.S., S.P., F.G.; funding acquisition, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This manuscript has not received any funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the University of Bologna Bioethics Committee (protocol code 0252770, 2 November 2020).

Informed Consent Statement

Informed consent was obtained from all participants involved in the study.

Data Availability Statement

All relevant data is contained within the article. All datasets generated for this study are included in the article/Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaper, J.B.; Nataro, J.P.; Mobley, H.L.T. Pathogenic Escherichia coli. Nat. Rev. Microbiol. Genet. 2004, 2, 123–140. [Google Scholar] [CrossRef]
  2. Poolman, J.T.; Wacker, M. Extraintestinal Pathogenic Escherichia coli, a Common Human Pathogen: Challenges for Vaccine Development and Progress in the Field. J. Infect. Dis. 2016, 213, 6–13. [Google Scholar] [CrossRef] [Green Version]
  3. McNulty, C.A.M.; Lecky, D.M.; Xu-McCrae, L.; Nakiboneka-Ssenabulya, D.; Chung, K.-T.; Nichols, T.; Thomas, H.L.; Thomas, M.; Alvarez-Buylla, A.; Turner, K.; et al. CTX-M ESBL-producing Enterobacteriaceae: Estimated prevalence in adults in England in 2014. J. Antimicrob. Chemother. 2018, 73, 1368–1388. [Google Scholar] [CrossRef]
  4. Rodríguez-Baño, J.; López-Cerero, L.; Navarro, M.D.; De Alba, P.D.; Pascual, A. Faecal carriage of extended-spectrum β-lactamase-producing Escherichia coli: Prevalence, risk factors and molecular epidemiology. J. Antimicrob. Chemother. 2008, 62, 1142–1149. [Google Scholar] [CrossRef] [Green Version]
  5. Clermont, O.; Bonacorsi, S.; Bingen, E. Rapid and Simple Determination of the Escherichia coli Phylogenetic Group. Appl. Environ. Microbiol. 2000, 66, 4555–4558. [Google Scholar] [CrossRef] [Green Version]
  6. Clermont, O.; Christenson, J.K.; Denamur, E.; Gordon, D.M. The Clermont Escherichia coli phylo-typing method revisited: Improvement of specificity and detection of new phylo-groups. Environ. Microbiol. Rep. 2013, 5, 58–65. [Google Scholar] [CrossRef]
  7. WHO. Critically Important Antimicrobials for Human Medicine, 6th Revision; WHO: Geneva, Switzerland, 2019. [Google Scholar]
  8. ECDC. Surveillance Report. Surveillance of Antimicrobial Resistance in Europe 2017; ECDC: Solna, Sweden, 2018. [Google Scholar]
  9. Johnson, J.R.; Russo, T.A. Extraintestinal pathogenic Escherichia coli: “The other bad E. coli”. J. Lab. Clin. Med. 2002, 139, 155–162. [Google Scholar] [CrossRef]
  10. Johnson, J.R.; Kuskowski, M.A.; O’Bryan, T.T.; Maslow, J.N. Epidemiological Correlates of Virulence Genotype and Phylogenetic Background among Escherichia coli Blood Isolates from Adults with Diverse-Source Bacteremia. J. Infect. Dis. 2002, 185, 1439–1447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Köhler, C.-D.; Dobrindt, U. What defines extraintestinal pathogenic Escherichia coli? Int. J. Med. Microbiol. 2011, 301, 642–647. [Google Scholar] [CrossRef] [PubMed]
  12. Clermont, O.; Olier, M.; Hoede, C.; Diancourt, L.; Brisse, S.; Keroudean, M.; Glodt, J.; Picard, B.; Oswald, E.; Denamur, E. Animal and human pathogenic Escherichia coli strains share common genetic backgrounds. Infect. Genet. Evol. 2011, 11, 654–662. [Google Scholar] [CrossRef] [PubMed]
  13. Johnson, J.R.; Owens, K.; Gajewski, A.; Clabots, C. Escherichia coli Colonization Patterns among Human Household Members and Pets, with Attention to Acute Urinary Tract Infection. J. Infect. Dis. 2008, 197, 218–224. [Google Scholar] [CrossRef] [Green Version]
  14. Tenaillon, O.; Skurnik, D.; Picard, B.; Denamur, E. The population genetics of commensal Escherichia coli. Nat. Rev. Microbiol. Genet. 2010, 8, 207–217. [Google Scholar] [CrossRef]
  15. Jang, J.; Hur, H.-G.; Sadowsky, M.; Byappanahalli, M.; Yan, T.; Ishii, S. Environmental Escherichia coli: Ecology and public health implications—A review. J. Appl. Microbiol. 2017, 123, 570–581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. EFSA. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2017. EFSA J. 2019, 17, e05598. [Google Scholar]
  17. Fuhrmeister, A.S.; Jones, R.N. The Importance of Antimicrobial Resistance Monitoring Worldwide and the Origins of SENTRY An-timicrobial Surveillance Program. In Open Forum Infectious Diseases; Oxford University Press: Oxford, UK, 2019. [Google Scholar]
  18. Tacconelli, E.; Sifakis, F.; Harbarth, S.; Schrijver, R.; van Mourik, M.; Voss, A.; Sharland, M.; Rajendran, N.B.; Rodríguez-Baño, J.; Bielicki, J.; et al. Surveillance for control of antimicrobial resistance. Lancet Infect. Dis. 2018, 18, e99–e106. [Google Scholar] [CrossRef] [Green Version]
  19. ECDC/EFSA/EMA. Second joint report on the integrated analysis of the consumption of antimicrobial agents and occur-rence of antimicrobial resistance in bacteria from humans and food-producing animals. EFSA J. 2017, 15, 4872. [Google Scholar]
  20. Seiffert, S.N.; Tinguely, R.; Lupo, A.; Neuwirth, C.; Perreten, V.; Endimiani, A. High prevalence of extend-ed-spectrum-cephalosporin-resistant Enterobacteriaceae in poultry meat in Switzerland: Emergence of CMY-2-and VEB-6-possessing Proteus mirabilis. Antimicrob. Agents Chemother. 2013, 57, 6406–6408. [Google Scholar] [CrossRef] [Green Version]
  21. Blaak, H.; De Kruijf, P.; Hamidjaja, R.A.; Van Hoek, A.H.; de Roda Husman, A.M.; Schets, F.M. Prevalence and characteristics of ESBL-producing E. coli in Dutch recreational waters influenced by wastewater treatment plants. Vet. Microbiol. 2014, 171, 448–459. [Google Scholar] [CrossRef]
  22. Dahms, C.; Hübner, N.-O.; Kossow, A.; Mellmann, A.; Dittmann, K.; Kramer, A. Occurrence of ESBL-Producing Escherichia coli in Livestock and Farm Workers in Mecklenburg-Western Pomerania, Germany. PLoS ONE 2015, 10, e0143326. [Google Scholar] [CrossRef]
  23. Economou, V.; Gousia, P. Agriculture and food animals as a source of antimicrobial-resistant bacteria. Infect. Drug Resist. 2015, 8, 49–61. [Google Scholar] [CrossRef] [Green Version]
  24. Fletcher, S. Understanding the contribution of environmental factors in the spread of antimicrobial resistance. Environ. Health Prev. Med. 2015, 20, 243–252. [Google Scholar] [CrossRef] [Green Version]
  25. Verraes, C.; Van Boxstael, S.; Van Meervenne, E.; Van Coillie, E.; Butaye, P.; Catry, B.; De Schaetzen, M.-A.; Van Huffel, X.; Imberechts, H.; Dierick, K.; et al. Antimicrobial Resistance in the Food Chain: A Review. Int. J. Environ. Res. Public Health 2013, 10, 2643–2669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Anaissie, E.J.; Penzak, S.R.; Dignani, M.C. The hospital water supply as a source of nosocomial infections: A plea for action. Arch. Intern. Med. 2002, 162, 1483–1492. [Google Scholar] [CrossRef] [Green Version]
  27. Ewers, C.; Bethe, A.; Semmler, T.; Guenther, S.; Wieler, L.H. 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, 646–655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Carroll, D.; Wang, J.; Fanning, S.; McMahon, B.J. Antimicrobial Resistance in Wildlife: Implications for Public Health. Zoonoses Public Health 2015, 62, 534–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Van Schaik, W. The human gut resistome. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2015, 370, 1670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Madec, J.-Y.; Haenni, M. Antimicrobial resistance plasmid reservoir in food and food-producing animals. Plasmid 2018, 99, 72–81. [Google Scholar] [CrossRef]
  31. Von Wintersdorff, C.J.H.; Penders, J.; Van Niekerk, J.M.; Mills, N.D.; Majumder, S.; Van Alphen, L.B.; Savelkoul, P.H.M.; Wolffs, P.F.G. Dissemination of Antimicrobial Resistance in Microbial Ecosystems through Horizontal Gene Transfer. Front. Microbiol. 2016, 7, 173. [Google Scholar] [CrossRef] [Green Version]
  32. Guenther, S.; Ewers, C.; Wieler, L.H. Extended-Spectrum Beta-Lactamases Producing E. coli in Wildlife, yet Another Form of Environmental Pollution? Front. Microbiol. 2011, 2, 246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Massella, E.; Reid, C.J.; Cummins, M.L.; Anantanawat, K.; Zingali, T.; Serraino, A.; Piva, S.; Giacometti, F.; Djordjevic, S.P. Snapshot Study of Whole Genome Sequences of Escherichia coli from Healthy Companion Animals, Livestock, Wildlife, Humans and Food in Italy. Antibiotics 2020, 9, 782. [Google Scholar] [CrossRef] [PubMed]
  34. Hunt, M.; Mather, A.E.; Sánchez-Busó, L.; Page, A.J.; Parkhill, J.; Keane, J.A.; Harris, S.R. ARIBA: Rapid antimicrobial re-sistance genotyping directly from sequencing reads. Microb. Genom. 2017, 3, e000131. [Google Scholar] [CrossRef] [PubMed]
  35. Ewers, C.; Li, G.; Wilking, H.; Kiessling, S.; Alt, K.; Antáo, E.M.; Laturnus, C.; Diehl, I.; Glodde, S.; Homeier, T.; et al. Avian pathogenic, uropathogenic, and newborn meningitis-causing Escherichia coli: How closely related are they? Int. J. Med. Microbiol. 2007, 297, 163–176. [Google Scholar] [CrossRef]
  36. Dierikx, C.; van der Goot, J.; Fabri, T.; van Essen-Zandbergen, A.; Smith, H.; Mevius, D. Extended-spectrum-β-lactamase-and AmpC-β-lactamase-producing Escherichia coli in Dutch broilers and broiler farmers. J. Antimicrob. Chemother. 2013, 68, 60–67. [Google Scholar] [CrossRef]
  37. Pletinckx, L.J.; Verhegghe, M.; Crombe, F.; Dewulf, J.; De Bleecker, Y.; Rasschaert, G.; Butaye, P.; Goddeeris, B.M.; De Man, I. Evidence of possible methicillin-resistant Staphylococcus aureus ST398 spread between pigs and other animals and people re-siding on the same farm. Prev. Vet. Med. 2013, 109, 293–303. [Google Scholar] [CrossRef]
  38. Van Hoek, A.H.; Veenman, C.; Van Overbeek, W.M.; Lynch, G.; De Roda Husman, A.M.; Blaak, H. Prevalence and charac-terization of ESBL- and AmpC-producing Enterobacteriaceae on retail vegetables. Int. J. Food Microbiol. 2015, 204, 1–8. [Google Scholar] [CrossRef]
  39. Davis, G.S.; Waits, K.; Nordstrom, L.; Weaver, B.; Aziz, M.; Gauld, L.; Grande, H.; Bigler, R.; Horwinski, J.; Porter, S.B.; et al. Intermingled Klebsiella pneumoniae Populations Between Retail Meats and Human Urinary Tract Infections. Clin. Infect. Dis. 2015, 61, 892–899. [Google Scholar] [CrossRef] [Green Version]
  40. Florez-Cuadrado, D.; Moreno, M.A.; Ugarte-Ruiz, M.; Dominguez, L. Antimicrobial Resistance in the Food Chain in the Eu-ropean Union. Adv. Food Nutr. Res. 2018, 86, 115–136. [Google Scholar]
  41. Hudson, J.A.; Frewer, L.J.; Jones, G.; Brereton, P.A.; Whittingham, M.J.; Stewart, G. The agri-food chain and antimicrobial resistance: A review. Trends Food Sci. Technol. 2017, 69, 131–147. [Google Scholar] [CrossRef] [Green Version]
  42. Manges, A. Escherichia coli and urinary tract infections: The role of poultry-meat. Clin. Microbiol. Infect. 2016, 22, 122–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bywater, R.; Silley, P.; Simjee, S. Antimicrobial breakpoints—Definitions and conflicting requirements. Vet. Microbiol. 2006, 118, 158–159. [Google Scholar] [CrossRef] [PubMed]
  44. Simjee, S.; Silley, P.; Werling, H.O.; Bywater, R. Potential confusion regarding the term ’resistance’ in epidemiological surveys. J. Antimicrob. Chemother. 2007, 61, 228–229. [Google Scholar] [CrossRef] [PubMed]
  45. EMA. European Surveillance of Veterinary Antimicrobial Consumption (ESVAC). 2019. Available online: https://www.ema.europa.eu/en/veterinary-regulatory/overview/antimicrobial-resistance/european-surveillance-veterinary-antimicrobial-consumption-esvac (accessed on 2 October 2019).
  46. De Verdier, K.; Nyman, A.; Greko, C.; Bengtsson, B. Antimicrobial resistance and virulence factors in Escherichia coli from swedish dairy calves. Acta Vet. Scand. 2012, 54, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Dotto, G.; Giacomelli, M.; Grilli, G.; Ferrazzi, V.; Carattoli, A.; Fortini, D.; Piccirillo, A. High Prevalence of oqxAB in Escherichia coli Isolates from Domestic and Wild Lagomorphs in Italy. Microb. Drug Resist. 2014, 20, 118–123. [Google Scholar] [CrossRef]
  48. EFSA. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2016. EFSA J. 2018, 16, 5182. [Google Scholar]
  49. EFSA. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2014. EFSA J. 2016, 14, 4380. [Google Scholar]
  50. EFSA. The European Union Summary Report on Antimicrobial Resistance in zoonotic and indicator bacteria from humans, animals and food in 2017/2018. EFSA J. 2020, 18, 6007. [Google Scholar]
  51. Skockova, A.; Bogdanovicoa, K.; Kolackova, I.; Karpiskova, R. Antimicrobial-resistant and extended-spectrum be-ta-lactamase-producing Escherichia coli in raw cow’s milk. J. Food Prot. 2015, 78, 72–77. [Google Scholar] [CrossRef]
  52. EFSA. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2015. EFSA J. 2017, 15, 4694. [Google Scholar]
  53. Bailey, J.K.; Pinyon, J.L.; Anantham, S.; Hall, R.M. Commensal Escherichia coli of healthy humans: A reservoir for antibi-otic-resistance determinants. J. Med. Microbiol. 2010, 59, 1331–1339. [Google Scholar] [CrossRef]
  54. Carvalho, A.; Barbosa, A.; Arais, L.; Ribeiro, P.; Carneiro, V.; Cerqueira, A. Resistance patterns, ESBL genes, and genetic re-latedness of Escherichia coli from dogs and owners. Braz. J. Microbiol. 2016, 47, 150–158. [Google Scholar] [CrossRef] [Green Version]
  55. Costa, D.; Poeta, P.; Sáenz, Y.; Coelho, A.C.; Matos, M.; Vinué, L.; Rodrigues, J.; Torres, C. Prevalence of antimicrobial resistance and resistance genes in faecal Escherichia coli isolates recovered from healthy pets. Vet. Microbiol. 2008, 127, 97–105. [Google Scholar] [CrossRef]
  56. Grevskott, D.H.; Svanevik, C.S.; Sunde, M.; Wester, A.L.; Lunestad, B.T. Marine Bivalve Mollusks As Possible Indicators of Multidrug-Resistant Escherichia coli and Other Species of the Enterobacteriaceae Family. Front. Microbiol. 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  57. Holvoet, K.; Sampers, I.; Callens, B.; Dewulf, J.; Uyttendaele, M. Moderate Prevalence of Antimicrobial Resistance in Escherichia coli Isolates from Lettuce, Irrigation Water, and Soil. Appl. Environ. Microbiol. 2013, 79, 6677–6683. [Google Scholar] [CrossRef] [Green Version]
  58. Nobili, G.; Franconieri, I.; Basanisi, M.; La Bella, G.; Tozzoli, R.; Caprioli, A.; La Salandra, G. Short communication: Isolation of Shiga toxin-producing Escherichia coli in raw milk and mozzarella cheese in southern Italy. J. Dairy Sci. 2016, 99, 7877–7880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Skurnik, D.; Bonnet, D.; Bernède-Bauduin, C.; Michel, R.; Guette, C.; Becker, J.M.; Balaire, C.; Chau, F.; Mohler, J.; Jarlier, V. Characteristics of human intestinal Escherichia coli with changing environments. Environ. Microbiol. 2008, 10, 2132–2137. [Google Scholar] [CrossRef] [PubMed]
  60. ECDC. Antimicrobial Consumption. Annual Epidemiological Report for 2017; ECDC: Solna, Sweden, 2018. [Google Scholar]
  61. Rozwandowicz, M.; Brouwer, M.S.M.; Fischer, J.; Wagenaar, J.A.; Gonzalez-Zorn, B.; Guerra, B.; Mevius, D.J.; Hordijk, J. Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae. J. Antimicrob. Chemother. 2018, 73, 1121–1137. [Google Scholar] [CrossRef] [Green Version]
  62. Pezzella, C.; Ricci, A.; Digiannatale, E.; Luzzi, I.; Carattoli, A. Tetracycline and Streptomycin Resistance Genes, Transposons, and Plasmids in Salmonella enterica Isolates from Animals in Italy. Antimicrob. Agents Chemother. 2004, 48, 903–908. [Google Scholar] [CrossRef] [Green Version]
  63. Reid, C.J.; Chowdhury, P.R.; Djordjevic, S.P. Tn6026 and Tn6029 are found in complex resistance regions mobilised by diverse plasmids and chromosomal islands in multiple antibiotic resistant Enterobacteriaceae. Plasmid 2015, 80, 127–137. [Google Scholar] [CrossRef]
  64. Silver, L.L. Challenges of Antibacterial Discovery. Clin. Microbiol. Rev. 2011, 24, 71–109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. FVE. Antimicrobial Use in Food-Producing Animals. Replies to EFSA/EMA Questions on the Use of Antimicrobials in Food-Producing Animals in EU and Possible Measures to Reduce Antimicrobial Use. 2016. Available online: https://www.ema.europa.eu/en/documents/report/annex-replies-efsa/ema-questions-use-antimicrobials-food-producing-animals-eu-possible-measures-reduce-antimicrobial_en.pdf (accessed on 30 June 2020).
  66. Federica, G.; Pezzi, A.; Galletti, G.; Tamba, M.; Merialdi, G.; Piva, S.; Serraino, A.; Rubini, S. Antimicrobial resistance patterns in Salmonella enterica subsp. enterica and Escherichia coli isolated from bivalve molluscs and marine environment. Food Control. 2021, 121, 107590. [Google Scholar]
  67. Walczak, J.J.; Xu, S. Manure as a Source of Antibiotic-Resistant Escherichia coli and Enterococci: A Case Study of a Wisconsin, USA Family Dairy Farm. Water Air Soil Pollut. 2011, 219, 579–589. [Google Scholar] [CrossRef]
  68. Alves, M.S.; Pereira, A.; Araújo, S.M.; Castro, B.B.; Correia, A.C.M.; Henriques, I. Seawater is a reservoir of multi-resistant Escherichia coli, including strains hosting plasmid-mediated quinolones resistance and extended-spectrum beta-lactamases genes. Front. Microbiol. 2014, 5, 426. [Google Scholar] [CrossRef]
  69. Baliere, C.; Rince, A.; Blanco, J.; Dahbi, G.; Harel, J.; Vogeleer, P.; Giard, J.C.; Mariani-Kurkdjian, P.; Gourmelon, M. Prevalence and Characterization of Shiga Toxin-Producing and Enteropathogenic Escherichia coli in Shellfish-Harvesting Areas and Their Watersheds. Front. Microbiol. 2015, 6, 1356. [Google Scholar] [CrossRef] [PubMed]
  70. Martinez, J.L. Environmental pollution by antibiotics and by antibiotic resistance determinants. Environ. Pollut. 2009, 157, 2893–2902. [Google Scholar] [CrossRef] [PubMed]
  71. Allen, H.K.; Donato, J.; Wang, H.H.; Cloud-Hansen, K.A.; Davies, J.; Handelsman, J. Call of the wild: Antibiotic resistance genes in natural environments. Nat. Rev. Microbiol. Genet. 2010, 8, 251–259. [Google Scholar] [CrossRef]
  72. Skurnik, D.; Ruimy, R.; Andremont, A.; Amorin, C.; Rouquet, P.; Picard, B.; Denamur, E. Effect of human vicinity on antimi-crobial resistance and integrons in animal faecal Escherichia coli. J. Antimicrob Chemother. 2006, 57, 1215–1219. [Google Scholar] [CrossRef]
  73. Sousa, M.; Torres, C.; Barros, J.; Somalo, S.; Igrejas, G.; Poeta, P. Gilthead Seabream (Sparus aurata) as Carriers of SHV-12 and TEM-52 Extended-Spectrum Beta-Lactamases-Containing Escherichia coli Isolates. Foodborne Pathog. Dis. 2011, 8, 1139–1141. [Google Scholar] [CrossRef]
  74. Campos, J.; Mourao, J.; Pestana, N.; Peixe, L.; Novais, C.; Antunes, P. Microbiological quality of ready-to-eat salads: An un-derestimated vehicle of bacteria and clinically relevant antibiotic resistance genes. Int J. Food Microbiol. 2013, 166, 464–470. [Google Scholar] [CrossRef]
  75. Literak, I.; Dolejska, M.; Radimersky, T.; Klimes, J.; Friedman, M.; Aarestrup, F.M.; Hasman, H.; Cizek, A. Antimicrobial-resistant faecal Escherichia coli in wild mammals in central Europe: Multiresistant Escherichia coli producing extended-spectrum beta-lactamases in wild boars. J. Appl. Microbiol. 2010, 108, 1702–1711. [Google Scholar] [CrossRef]
  76. Giacopello, C.; Foti, M.; Mascetti, A.; Grosso, F.; Ricciardi, D.; Fisichella, V.; Piccolo, F.L. Antimicrobial resistance patterns of Enterobacteriaceae in European wild bird species admitted in a wildlife rescue centre. Vet. Ital. 2016, 52, 139–144. [Google Scholar]
  77. De Jong, A.; Stephan, B.; Silley, P. Fluoroquinolone resistance of Escherichia coli and Salmonella from healthy livestock and poultry in the EU. J. Appl. Microbiol. 2011, 112, 239–245. [Google Scholar] [CrossRef]
  78. Hricova, K.; Roderova, M.; Pudova, V.; Hanulik, V.; Halova, D.; Julinkova, P.; Dolejska, M.; Papousek, I.; Bardon, J. Quino-lone-resistant Escherichia coli in Poultry Farming. Cent. Eur. J. Public Health 2017, 25, 163–167. [Google Scholar] [CrossRef]
  79. Literak, I.; Reitschmied, T.; Bujnakova, D.; Dolejska, M.; Cizek, A.; Bardon, J.; Pokludova, L.; Alexa, P.; Halova, D.; Jamborova, I. Broilers as a source of quinolone-resistant and extraintestinal pathogenic Escherichia coli in the Czech Republic. Microb. Drug Resist. 2013, 19, 57–63. [Google Scholar] [CrossRef] [PubMed]
  80. Costa, D.; Poeta, P.; Sáenz, Y.; Vinué, L.; Coelho, A.C.; Matos, M.; Rojo-Bezares, B.; Rodrigues, J.; Torres, C. Mechanisms of Antibiotic Resistance in Escherichia coli Isolates Recovered from Wild Animals. Microb. Drug Resist. 2008, 14, 71–77. [Google Scholar] [CrossRef] [PubMed]
  81. Guenther, S.; Grobbel, M.; Lübke-Becker, A.; Goedecke, A.; Friedrich, N.D.; Wieler, L.H.; Ewers, C. Antimicrobial resistance profiles of Escherichia coli from common European wild bird species. Veter. Microbiol. 2010, 144, 219–225. [Google Scholar] [CrossRef] [PubMed]
  82. Johnson, T.J.; Logue, C.M.; Johnson, J.R.; Kuskowski, M.A.; Sherwood, J.S.; Barnes, H.J.; Debroy, C.; Wannemuehler, Y.M.; Obata-Yasuoka, M.; Spanjaard, L.; et al. Associations between multidrug resistance, plasmid content, and virulence potential among extraintestinal pathogenic and commensal Escherichia coli from humans and poultry. Foodborne Pathog. Dis. 2012, 9, 37–46. [Google Scholar] [CrossRef] [Green Version]
  83. Tadesse, D.A.; Zhao, S.; Tong, E.; Ayers, S.; Singh, A.; Bartholomew, M.J.; Mcdermott, P.F. Antimicrobial drug resistance in Escherichia coli from humans and food animals, United States, 1950–2002. J. Emerg. Infect. Dis. 2012, 18, 741. [Google Scholar] [CrossRef]
  84. Robinson, T.P.; Thornton, P.K.; Franceschini, G.; Kruska, R.L.; Chiozza, F.; Notenbaert, A.M.O.; Cecchi, G.; Herrero, M.; Epprecht, M.; Fritz, S.; et al. Global Livestock Production Systems; FAO: Rome, Italy; ILRI: Nairobi, Kenya, 2011. [Google Scholar]
  85. Jensen, V.; Jacobsen, E.; Bager, F. Veterinary antimicrobial-usage statistics based on standardized measures of dosage. Prev. Vet. Med. 2004, 64, 201–215. [Google Scholar] [CrossRef]
  86. Chantziaras, I.; Boyen, F.; Callens, B.; Dewulf, J. Correlation between veterinary antimicrobial use and antimicrobial resistance in food-producing animals: A report on seven countries. J. Antimicrob. Chemother. 2014, 69, 827–834. [Google Scholar] [CrossRef] [Green Version]
  87. Blanco, J.E.; Blanco, M.; Blanco, J.; Rioja, L.; Ducha, J. Serotypes, toxins and antibiotic resistance of Escherichia coli strains isolated from diarrhoeic and healthy rabbits in Spain. Vet. Microbiol. 1994, 38, 193–201. [Google Scholar] [CrossRef]
  88. Freitas-Silva, J.; Inácio, Â.S.; Mourão, J.; Antunes, P.; Mendes, Â.; De Carvalho, A.P.; Vasconcelos, V.; Peixe, L.; Da Costa, P.M. Occurrence of mcr-1 in Escherichia coli from rabbits of intensive farming. Veter. Microbiol. 2018, 227, 78–81. [Google Scholar] [CrossRef] [Green Version]
  89. Medina, A.; De La Fuente, R.; Badiola, I.; De Rozas, A.P.; Orden, J.A.; Ruiz-Santa Quiteria, J.A. Antimicrobial Resistance in Commensal Escherichia coli Isolates from Rabbits. Rev. Cient. 2015, 25, 38–42. [Google Scholar]
  90. Agnoletti, F.; Brunetta, R.; Bano, L.; Drigo, I.; Mazzolini, E. Longitudinal study on antimicrobial consumption and resistance in rabbit farming. Int. J. Antimicrob. Agents. 2018, 51, 197–205. [Google Scholar] [CrossRef]
  91. Johnson, J.R.; Kuskowski, M.A.; Menard, M.; Gajewski, A.; Xercavins, M.; Garau, J. Similarity between Human and Chicken Escherichia coli Isolates in Relation to Ciprofloxacin Resistance Status. J. Infect. Dis. 2006, 194, 71–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Johnson, J.R.; Kuskowski, M.A.; Smith, K.; O’bryan, T.T.; Tatini, S. Antimicrobial-resistant and extraintestinal pathogenic Escherichia coli in retail foods. J. Infect. Dis. 2005, 191, 1040–1049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Johnson, J.R.; Porter, S.B.; Johnston, B.; Thuras, P.; Clock, S.; Crupain, M.; Rangan, U. Extraintestinal Pathogenic and Anti-microbial-Resistant Escherichia coli, Including Sequence Type 131 (ST131), from Retail Chicken Breasts in the United States in 2013. Appl. Environ. Microbiol. 2017, 83, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Johnson, J.R.; Sannes, M.R.; Croy, C.; Johnston, B.; Clabots, C.; Kuskowski, M.A.; Bender, J.; Smith, K.E.; Winokur, P.L.; Be-longia, E.A. Antimicrobial drug-resistant Escherichia coli from humans and poultry products, Minnesota and Wisconsin, 2002–2004. Emerg. Infect. Dis. 2007, 13, 838–846. [Google Scholar] [CrossRef] [PubMed]
  95. Hopkins, K.L.; Wootton, L.; Day, M.R.; Threlfall, E.J. Plasmid-mediated quinolone resistance determinant qnrS1 found in Salmonella enterica strains isolated in the UK. J. Antimicrob. Chemother. 2007, 59, 1071–1075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Hooper, D.C.; Jacoby, G.A. Mechanisms of drug resistance: Quinolone resistance. Ann. N. Y. Acad. Sci. 2015, 1354, 12–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Sáenz, Y.; Zarazaga, M.; Briñas, L.; Ruiz-Larrea, F.; Torres, C. Mutations in gyrA and parC genes in nalidixic acid-resistant Escherichia coli strains from food products, humans and animals. J. Antimicrob. Chemother. 2003, 51, 1001–1005. [Google Scholar] [CrossRef]
  98. Paltansing, S.; Kraakman, M.E.M.; Ras, J.M.C.; Wessels, E.; Bernards, A.T. Characterization of fluoroquinolone and cephalosporin resistance mechanisms in Enterobacteriaceae isolated in a Dutch teaching hospital reveals the presence of an Escherichia coli ST131 clone with a specific mutation in parE. J. Antimicrob. Chemother. 2012, 68, 40–45. [Google Scholar] [CrossRef] [Green Version]
  99. Schwarz, S.; Kehrenberg, C.; Doublet, B.; Cloeckaert, A. Molecular basis of bacterial resistance to chloramphenicol and florfenicol. FEMS Microbiol. Rev. 2004, 28, 519–542. [Google Scholar] [CrossRef] [Green Version]
  100. Poole, K. Efflux-Mediated Antimicrobial Resistance. J. Antimicrob. Chemother. 2005, 56, 20–51. [Google Scholar] [CrossRef] [Green Version]
  101. Piddock, L.J.V. Multidrug-resistance efflux pumps? Not just for resistance. Nat. Rev. Microbiol. Genet. 2006, 4, 629–636. [Google Scholar] [CrossRef] [PubMed]
  102. Poirel, L.; Jayol, A.; Nordmann, P. Polymyxins: Antibacterial Activity, Susceptibility Testing, and Resistance Mechanisms Encoded by Plasmids or Chromosomes. Clin. Microbiol. Rev. 2017, 30, 557–596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Quesada, A.; Porrero, M.C.; Téllez, S.; Palomo, G.; García, M.; Domínguez, L. Polymorphism of genes encoding PmrAB in colistin-resistant strains of Escherichia coli and Salmonella enterica isolated from poultry and swine. J. Antimicrob. Chemother. 2015, 70, 71–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Zankari, E.; Allesøe, R.; Joensen, K.G.; Cavaco, L.M.; Lund, O.; Aarestrup, F.M. 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, 2764–2768. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Rafique, M.; Potter, R.F.; Ferreiro, A.; Wallace, M.A.; Rahim, A.; Malik, A.A.; Siddique, N.; Abbas, M.A.; D’Souza, A.W.; Burnham, C.-A.D.; et al. Genomic Characterization of Antibiotic Resistant Escherichia coli Isolated From Domestic Chickens in Pakistan. Front. Microbiol. 2020, 10, 3052. [Google Scholar] [CrossRef] [PubMed]
  106. Gunn, J.S. The Salmonella PmrAB regulon: Lipopolysaccharide modifications, antimicrobial peptide resistance and more. Trends Microbiol. 2008, 16, 284–290. [Google Scholar] [CrossRef]
  107. Kempf, I.; Fleury, M.A.; Drider, D.; Bruneau, M.; Sanders, P.; Chauvin, C.; Madec, J.Y.; Jouy, E. What do we know about re-sistance to colistin in Enterobacteriaceae in avian and pig production in Europe? Int. J. Antimicrob. Agents 2013, 42, 379–383. [Google Scholar] [CrossRef]
  108. Carlos, C.; Pires, M.M.; Stoppe, N.C.; Hachich, E.M.; Sato, M.I.; Gomes, T.A.T.; Amaral, L.A.; Ottoboni, L.M.M. Escherichia coli phylogenetic group determination and its application in the identification of the major animal source of fecal contamination. BMC Microbiol. 2010, 10, 161. [Google Scholar] [CrossRef] [Green Version]
  109. Escobar-Páramo, P.; Le Menac’H, A.; Le Gall, T.; Amorin, C.; Gouriou, S.; Picard, B.; Skurnik, D.; Denamur, E. Identification of forces shaping the commensal Escherichia coli genetic structure by comparing animal and human isolates. Environ. Microbiol. 2006, 8, 1975–1984. [Google Scholar] [CrossRef]
  110. Smati, M.; Clermont, O.; Bleibtreu, A.; Fourreau, F.; David, A.; Daubié, A.; Hignard, C.; Loison, O.; Picard, B.; Denamur, E. Quantitative analysis of commensal Escherichia coli populations reveals host-specific enterotypes at the intra-species level. Microbiologyopen 2015, 4, 604–615. [Google Scholar] [CrossRef]
  111. Alonso, C.A.; González-Barrio, D.; Ruiz-Fons, F.; Ruiz-Ripa, L.; Torres, C. High frequency of B2 phylogroup among non-clonally related fecal Escherichia coli isolates from wild boars, including the lineage ST131. FEMS Microbiol. Ecol. 2017, 93, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Gordon, D.M.; Cowling, A. The distribution and genetic structure of Escherichia coli in Australian vertebrates: Host and ge-ographic effects. Microbiology 2003, 149, 3575–3586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Stoppe, N.D.C.; Silva, J.S.; Carlos, C.; Sato, M.I.Z.; Saraiva, A.M.; Ottoboni, L.M.M.; Torres, T.T. Worldwide Phylogenetic Group Patterns of Escherichia coli from Commensal Human and Wastewater Treatment Plant Isolates. Front. Microbiol. 2017, 8, 2512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Jakobsen, L.; Kurbasic, A.; Skjøt-Rasmussen, L.; Ejrnæs, K.; Porsbo, L.J.; Pedersen, K.; Jensen, L.B.; Emborg, H.-D.; Agersø, Y.; Olsen, K.E.; et al. Escherichia coli Isolates from Broiler Chicken Meat, Broiler Chickens, Pork, and Pigs Share Phylogroups and Antimicrobial Resistance with Community-Dwelling Humans and Patients with Urinary Tract Infection. Foodborne Pathog. Dis. 2010, 7, 537–547. [Google Scholar] [CrossRef]
  115. Ferreira, J.C.; Filho, R.A.C.P.; Kuaye, A.P.Y.; Andrade, L.N.; Chang, Y.-F.; Darini, A.L.C. Virulence potential of commensal multidrug resistant Escherichia coli isolated from poultry in Brazil. Infect. Genet. Evol. 2018, 65, 251–256. [Google Scholar] [CrossRef] [Green Version]
  116. Smati, M.; Clermont, O.; Le Gal, F.; Schichmanoff, O.; Jauréguy, F.; Eddi, A.; Denamur, E.; Picard, B.; for the Coliville Group. Real-Time PCR for Quantitative Analysis of Human Commensal Escherichia coli Populations Reveals a High Frequency of Subdominant Phylogroups. Appl. Environ. Microbiol. 2013, 79, 5005–5012. [Google Scholar] [CrossRef] [Green Version]
  117. Schmidt, V.M.; Pinchbeck, G.L.; Nuttall, T.; McEwan, N.; Dawson, S.; Williams, N.J. Antimicrobial resistance risk factors and characterisation of faecal E. coli isolated from healthy Labrador retrievers in the United Kingdom. Prev. Vet. Med. 2015, 119, 31–40. [Google Scholar] [CrossRef] [Green Version]
  118. Lescat, M.; Clermont, O.; Woerther, P.L.; Glodt, J.; Dion, S.; Skurnik, D.; Djossou, F.; Dupont, C.; Perroz, G.; Picard, B.; et al. Commensal Escherichia coli strains in Guiana reveal a high genetic diversity with host-dependant population structure. Environ. Microbiol. Rep. 2013, 5, 49–57. [Google Scholar] [CrossRef] [PubMed]
  119. Bok, E.; Mazurek, J.; Stosik, M.; Wojciech, M.; Baldy-Chudzik, K. Prevalence of virulence determinants and antimicrobial re-sistance among commensal Escherichia coli derived from dairy and beef cattle. Int. J. Environ. Res. Public Health 2015, 12, 970–985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  120. Hamelin, K.; Bruant, G.; El-Shaarawi, A.; Hill, S.; Edge, T.A.; Fairbrother, J.; Harel, J.; Maynard, C.; Masson, L.; Brousseau, R. Occurrence of Virulence and Antimicrobial Resistance Genes in Escherichia coli Isolates from Different Aquatic Ecosystems within the St. Clair River and Detroit River Areas. Appl. Environ. Microbiol. 2006, 73, 477–484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Picard, B.; Garcia, J.S.; Gouriou, S.; Duriez, P.; Brahimi, N.; Bingen, E.; Elion, J.; Denamur, E. The link between phylogeny and virulence in Escherichia coli extraintestinal infection. Infect. Immun. 1999, 67, 546–553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Johnson, J.R.; Kuskowski, M. Clonal origin, virulence factors, and virulence. Infect. Immun. 2000, 68, 424–425. [Google Scholar] [CrossRef] [Green Version]
  123. Sarowska, J.; Futoma-Koloch, B.; Jama-Kmiecik, A.; Frej-Madrzak, M.; Ksiazczyk, M.; Bugla-Ploskonska, G.; Choroszy-Krol, I. Virulence factors, prevalence and potential transmission of extraintestinal pathogenic Escherichia coli isolated from different sources: Recent reports. Gut Pathog. 2019, 11, 10. [Google Scholar] [CrossRef] [Green Version]
  124. Schierack, P.; Walk, N.; Ewers, C.; Wilking, H.; Steinrück, H.; Filter, M.; Wieler, L.H. ExPEC-typical virulence-associated genes correlate with successful colonization by intestinal E. coli in a small piglet group. Environ Microbiol. 2008, 10, 1742–1751. [Google Scholar] [CrossRef]
  125. Los, F.C.O.; Randis, T.M.; Aroian, R.V.; Ratner, A.J. Role of Pore-Forming Toxins in Bacterial Infectious Diseases. Microbiol. Mol. Biol. Rev. 2013, 77, 173–207. [Google Scholar] [CrossRef] [Green Version]
  126. Sedeek, D.M.; Rady, M.M.; Fedawy, H.S.; Rabie, N.S. Molecular Epidemiology and Sequencing of Avian Pathogenic Escherichia coli APEC in Egypt. Adv. Anim. Vet. Sci. 2020, 8, 499–505. [Google Scholar] [CrossRef]
  127. Nicolas-Chanoine, M.-H.; Bertrand, X.; Madec, J.-Y. Escherichia coli ST131, an Intriguing Clonal Group. Clin. Microbiol. Rev. 2014, 27, 543–574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Riley, L.W. Pandemic lineages of extraintestinal pathogenic Escherichia coli. Clin. Microbiol. Infect. 2014, 20, 380–390. [Google Scholar] [CrossRef] [Green Version]
  129. Dodson, K.W.; Pinkner, J.S.; Rose, T.; Magnusson, G.; Hultgren, S.J.; Waksman, G. Structural Basis of the Interaction of the Pyelonephritic E. coli Adhesin to Its Human Kidney Receptor. Cell 2001, 105, 733–743. [Google Scholar] [CrossRef] [Green Version]
  130. Lund, B.; Lindberg, F.; Marklund, B.I.; Normark, S. The PapG protein is the α-d-galactopyranosyl-(1→4)-β-d-galactopyranose-binding adhesin of uropathogenic Escherichia coli. Proc. Natl. Acad. Sci. USA 1987, 84, 5898–5902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  131. Lane, M.; Mobley, H. Role of P-fimbrial-mediated adherence in pyelonephritis and persistence of uropathogenic Escherichia coli (UPEC) in the mammalian kidney. Kidney Int. 2007, 72, 19–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  132. Katouli, M. Population structure of gut Escherichia coli and its role in development of extra-intestinal infections. Iran. J. Microbiol. 2010, 2, 59–72. [Google Scholar]
  133. Restieri, C.; Garriss, G.; Locas, M.-C.; Dozois, C.M. Autotransporter-Encoding Sequences Are Phylogenetically Distributed among Escherichia coli Clinical Isolates and Reference Strains. Appl. Environ. Microbiol. 2007, 73, 1553–1562. [Google Scholar] [CrossRef] [Green Version]
  134. Rijavec, M.; Müller-Premru, M.; Zakotnik, B.; Žgur-Bertok, D. Virulence factors and biofilm production among Escherichia coli strains causing bacteraemia of urinary tract origin. J. Med. Microbiol. 2008, 57, 1329–1334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Yamamoto, S.; Nakano, M.; Terai, A.; Yuri, K.; Nakata, K.; Nair, G.B.; Kurazono, H.; Ogawa, O. The presence of the virulence island containing the usp gene in uropathogenic Escherichia coli is associated with urinary tract infection in an experimental mouse model. J. Urol. 2001, 165, 1347–1351. [Google Scholar] [CrossRef]
  136. Johnson, J.R.; Manges, A.R.; O’Bryan, T.T.; Riley, L.W. A disseminated multidrug-resistant clonal group of uropathogenic Escherichia coli in pyelonephritis. Lancet 2002, 359, 2249–2251. [Google Scholar] [CrossRef]
  137. Villa, L.; García-Fernández, A.; Fortini, D.; Carattoli, A. Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. J. Antimicrob. Chemother. 2010, 65, 2518–2529. [Google Scholar] [CrossRef] [Green Version]
  138. Johnson, T.J.; Nolan, L.K. Pathogenomics of the virulence plasmids of Escherichia coli. Microbiol. Mol. Biol. Rev. 2009, 73, 750–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Da Silva, G.J.; Mendonça, N. Association between antimicrobial resistance and virulence in Escherichia coli. Virulence 2012, 3, 18–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  140. Moura, A.; Oliveira, C.; Henriques, I.; Smalla, K.; Correia, A. Broad diversity of conjugative plasmids in integron-carrying bacteria from wastewater environments. FEMS Microbiol. Lett. 2012, 330, 157–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  141. Pitout, J.D.D. Extraintestinal pathogenic Escherichia coli: An update on antimicrobial resistance, laboratory diagnosis and treatment. Expert Rev. Anti Infect. Ther. 2012, 10, 1165–1176. [Google Scholar] [CrossRef]
  142. Bauchart, P.; Germon, P.; Bree, A.; Oswald, E.; Hacker, J.; Dobrindt, U. Pathogenomic comparison of human extraintestinal and avian pathogenic Escherichia coli—Search for factors involved in host specificity or zoonotic potential. Microb. Pathog. 2010, 49, 105–115. [Google Scholar] [CrossRef]
  143. Li, G.; Cai, W.; Hussein, A.; Wannemuehler, Y.M.; Logue, C.M.; Nolan, L.K. Proteome response of an extraintestinal pathogenic Escherichia coli strain with zoonotic potential to human and chicken sera. J. Proteom. 2012, 75, 4853–4862. [Google Scholar] [CrossRef]
  144. Vincent, C.; Boerlin, P.; Daignault, D.; Dozois, C.M.; Dutil, L.; Galanakis, C.; Reid-Smith, R.J.; Tellier, P.-P.; Tellis, P.A.; Ziebell, K.; et al. Food Reservoir for Escherichia coli Causing Urinary Tract Infections. Emerg. Infect. Dis. 2010, 16, 88–95. [Google Scholar] [CrossRef]
  145. Bergeron, C.R.; Prussing, C.; Boerlin, P.; Daignault, D.; Dutil, L.; Reid-Smith, R.J.; Zhanel, G.G.; Manges, A.R. Chicken as res-ervoir for extraintestinal pathogenic Escherichia coli in humans, Canada. Emerg. Infect. Dis. 2012, 18, 415–421. [Google Scholar] [CrossRef]
  146. Mor-Mur, M.; Yuste, J. Emerging Bacterial Pathogens in Meat and Poultry: An Overview. Food Bioprocess Technol. 2010, 3, 24–35. [Google Scholar] [CrossRef]
  147. Lyhs, U.; Ikonen, I.; Pohjanvirta, T.; Raninen, K.; Perko-Mäkelä, P.; Pelkonen, S. Extraintestinal pathogenic Escherichia coli in poultry meat products on the Finnish retail market. Acta Vet. Scand. 2012, 54, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Mitchell, S.M.; Subbiah, M.; Jeffrey, L.; Ullman, J.L.; Frear, C.; Call, D.R. Evaluation of 27 different biochars for potential sequestration of antibiotic residues in food animal production environments. J. Environ. Chem. Eng. 2015, 3, 162–169. [Google Scholar] [CrossRef]
  149. Horakova, K.; Mlejnkova, H.; Mlejnek, P. Specific detection of Escherichia coli isolated from water samples using polymerase chain reaction targeting four genes: Cytochrome bd complex, lactose permease, β-d-glucuronidase, and β-d-galactosidase. J. Appl. Microbiol. 2008, 105, 970–976. [Google Scholar] [CrossRef]
  150. CLSI. M100: Performance Standards for Antimicrobial Susceptibility Testing, 28th ed.; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2018. [Google Scholar]
  151. EUCAST. EUCAST Guideline for the Detection of Resistance Mechanisms and Specific Resistances of Clinical and/or Epidemiological Importance, v2.0.; EUCAST: Växjö, Sweden, 2017. [Google Scholar]
  152. EUCAST. Breakpoint Tables for Interpretation of MICs and Zone Diameters, Version 9.0; EUCAST: Växjö, Sweden, 2019. [Google Scholar]
  153. CLSI. Performance Standards for Antimicrobial Disk and Dilution Susceptibility Tests for Bacteria Isolated from Animals; Informational Supplement: M31-S1, No. 17; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2004; Volume 24. [Google Scholar]
  154. CLSI. M100: Performance Standards for Antimicrobial Susceptibility Testing, 27th ed.; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2017. [Google Scholar]
  155. Magiorakos, A.-P.; Srinivasan, A.; Carey, R.; Carmeli, Y.; Falagas, M.; Giske, C.; Harbarth, S.; Hindler, J.; Kahlmeter, G.; Olsson-Liljequist, B.; et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: An international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 2012, 18, 268–281. [Google Scholar] [CrossRef] [Green Version]
  156. Clermont, O.; Gordon, D.M.; Brisse, S.; Walk, S.T.; Denamur, E. Characterization of the cryptic Escherichia lineages: Rapid identification and prevalence. Environ. Microbiol. 2011, 13, 2468–2477. [Google Scholar] [CrossRef] [Green Version]
  157. Manges, A.R.; Geum, H.M.; Guo, A.; Edens, T.J.; Fibke, C.D.; Pitout, J.D.D. Global Extraintestinal Pathogenic Escherichia coli (ExPEC) Lineages. Clin. Microbiol. Rev. 2019, 32. [Google Scholar] [CrossRef] [PubMed]
  158. Lindstedt, B.-A.; Finton, M.D.; Porcellato, D.; Brandal, L.T. High frequency of hybrid Escherichia coli strains with combined Intestinal Pathogenic Escherichia coli (IPEC) and Extraintestinal Pathogenic Escherichia coli (ExPEC) virulence factors isolated from human faecal samples. BMC Infect. Dis. 2018, 18, 544. [Google Scholar] [CrossRef] [PubMed]
  159. Kudinha, T.; Kong, F.; Johnson, J.R.; Andrew, S.D.; Anderson, P.; Gilbert, G.L. Multiplex PCR-based reverse line blot assay for simultaneous detection of 22 virulence genes in uropathogenic Escherichia coli. Appl Environ. Microbiol. 2012, 78, 1198–1202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  160. Tóth, I.; Nougayrède, J.P.; Dobrindt, U.; Ledger, T.N.; Boury, M.; Morabito, S.; Fujiwara, T.; Sugai, M.; Hacker, J.; Oswald, E. Cytolethal distending toxin type I and type IV genes are framed with lambdoid prophage genes in extraintestinal pathogenic Escherichia coli. Infect Immun. 2009, 77, 492–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  161. Ostblom, A.; Adlerberth, I.; Wold, A.E.; Nowrouzian, F.L. Pathogenicity island markers, virulence determinants malX and usp, and the capacity of Escherichia coli to persist in infants’ commensal microbiotas. Appl. Environ. Microbiol. 2011, 77, 2303–2308. [Google Scholar] [CrossRef] [Green Version]
  162. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; Urbana, University of Illinois Press: Champaign, IL, USA, 1949. [Google Scholar]
  163. Simpson, E.H. Measurement of Diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
Antibiotics 10 00351 g001 550
Figure 1. Heat map depicting phenotypic AMR and carriage of typical ExPEC VAGs in the strain collection. The dendogram on the left represents clustering of E. coli isolates according to their phenotypic AMR profile. Presence of phenotypic AMR and carriage of typical ExPEC VAGs are shown in red, with a green (for phenotypic AMR) or blue (for ExPEC VAGs) square indicating their absence. Additional strain information is provided in column 1- 3 and includes: source (column beside the dendogram), ST (first column after the heatmap) and phylogroup (second column after the heatmap). Sources, STs and phylogroups are colour-coded as described in the legend.
Figure 1. Heat map depicting phenotypic AMR and carriage of typical ExPEC VAGs in the strain collection. The dendogram on the left represents clustering of E. coli isolates according to their phenotypic AMR profile. Presence of phenotypic AMR and carriage of typical ExPEC VAGs are shown in red, with a green (for phenotypic AMR) or blue (for ExPEC VAGs) square indicating their absence. Additional strain information is provided in column 1- 3 and includes: source (column beside the dendogram), ST (first column after the heatmap) and phylogroup (second column after the heatmap). Sources, STs and phylogroups are colour-coded as described in the legend.
Antibiotics 10 00351 g001
Antibiotics 10 00351 g002 550
Figure 2. A mid-point rooted, maximum-likelihood phylogenetic tree of 279 commensal E. coli included in the study. The phylogenetic tree is rearranged from the original one, whose comprehensive description is available in the previous publication [33]. Branches are coloured by clade and subclade (red= clade 1, subclade 1; orange= clade 1, subclade 2; black= clade2, subclade 1; blue= clade 2, subclade 2; sea green= clade 2, subclade 3; light green= clade 2, subclade 4). Source (inner ring), phylogroup according to Clermont quadruplex PCR (ring 2), phylogroup according to in silico method (ring 3), emergent/pandemic/uncommon ST (ring 4), VAG number (ring 5), MDR profile (ring 6), HP-CIA resistance profile (outer ring) are annotated according to the legend. A comprehensive description of the phylogenetic tree is available in the previous publication [33].
Figure 2. A mid-point rooted, maximum-likelihood phylogenetic tree of 279 commensal E. coli included in the study. The phylogenetic tree is rearranged from the original one, whose comprehensive description is available in the previous publication [33]. Branches are coloured by clade and subclade (red= clade 1, subclade 1; orange= clade 1, subclade 2; black= clade2, subclade 1; blue= clade 2, subclade 2; sea green= clade 2, subclade 3; light green= clade 2, subclade 4). Source (inner ring), phylogroup according to Clermont quadruplex PCR (ring 2), phylogroup according to in silico method (ring 3), emergent/pandemic/uncommon ST (ring 4), VAG number (ring 5), MDR profile (ring 6), HP-CIA resistance profile (outer ring) are annotated according to the legend. A comprehensive description of the phylogenetic tree is available in the previous publication [33].
Antibiotics 10 00351 g002
Table
Table 1. Number of processed samples and number of isolated strains in the different sources considered in the study.
Table 1. Number of processed samples and number of isolated strains in the different sources considered in the study.
SourceNumber of SamplesNumber of Isolated E. coliE. coli Occurrence (%)
Beef332163.64
Wild Boar312270.97
Vegetable1642414.63
Fishery942425.53
Companion Animal1212100.00
Swine66100.00
Poultry332575.76
Rabbit141071.43
Human272592.59
Table
Table 2. AMR resistance rate to the molecules tested and ESBL profile among the sources investigated.
Table 2. AMR resistance rate to the molecules tested and ESBL profile among the sources investigated.
HP-CIA
SourcesnR ≥ 1MDRCNSCSTSXTTENAENREFTCAZCOLESBL
Dairy256 (24%)6 (24%)1 (4%)6 (24%)3 (12%)6 (24%)6 (24%)6 (24%)3 (12%)3 (12%)2 (8%)2 (8%)1 (4%)2 (8%)
Beef247 (29.17%)4 (16.66%)2 (8.33%)5 (20.83%)2 (8.33%)4 (16.67%)4 (16.67%)3 (12.5%)1 (4.17%)1 (4.17%)2 (8.33%)2 (8.33%)0 (0%)0 (0%)
Wild boar221 (4.55%)1 (4.55%)0 (0%)1 (4.55%)0 (0%)1 (4.55%)1 (4.55%)1 (4.55%)1 (4.55%)1 (4.55%)0 (0%)0 (0%)0 (0%)0 (0%)
Vegetable242 (8.33%)2 (8.33%)0 (0%)2 (8.33%)0 (0%)1 (4.17%)1 (4.17%)1 (4.55%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)
Fishery245 (20.83%)3 (12.5%)1 (4.17%)3 (12.5%)0 (0%)2 (8.33%)3 (12.5%)3 (12.5%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)
Companion animal124 (33.33%)3 (25%)0 (0%)1 (8.33%)1 (8.33%)3 (25%)3 (25%)4 (33.33%)1 (8.33%)1 (8.33%)0 (0%)0 (0%)0 (0%)0 (0%)
Swine2516 (64%)6 (24%)0 (0%)8 (32%)2 (8%)5 (20%)5 (20%)10 (40%)3 (12%)3 (12%)0 (0%)0 (0%)2 (8%)0 (0%)
Poultry2523 (92%)16 (64%)1 (4%)13 (52%)8 (32%)14 (56%)15 (60%)13 (52%)13 (52%)11 (44%)0 (0%)0 (0%)0 (0%)0 (0%)
Rabbit2323 (100%)22 (95.65%)6 (26.09%)22 (95.65%)3 (13.04%)22 (95.65%)22 (95.65%)22 (95.65%)15 (65.22%)13 (56.52%)0 (0%)0 (0%)3 (13.04%)0 (0%)
Mollusc254 (16%)3 (12%)0 (0%)2 (8%)0 (0%)1 (4%)1 (4%)2 (8%)1 (4%)1 (4%)0 (0%)0 (0%)0 (0%)0 (0%)
Human259 (36%)7 (28%)0 (0%)5 (20%)3 (12%)6 (24%)6 (24%)6 (24%)6 (24%)6 (24%)1 (4%)1 (4%)0 (0%)1 (4%)
Wild animal257 (28%)6 (24%)0 (0%)3 (12%)1 (4%)6 (24%)6 (24%)5 (20%)4 (16%)1 (4%)1 (4%)1 (4%)0 (0%)1 (4%)
Total279107 (38.35%)79 (28.32%)11 (3.94%)71 (25.45%)23 (8.24%)71 (25.45)73 (26.16%)76 (27.24%)48 (17.20%)44 (15.77%)6 (2.15%)6 (2.15%)6 (2.15%)4 (1.43%)
R ≥ 1: resistant to at least one antimicrobial; MDR: multiresistant isolate; CN: gentamicin; S: streptomycin; C: chloramphenicol; ST: sulfisoxazole; SXT: trimethoprim/sulfamethoxazole; TE: tetracycline; NA: nalidixic acid; ENR: enrofloxacin; EFT: ceftiofur; CAZ: ceftazidime; COL: colistin; ESBL: extended spectrum beta lactamase profile; HP-CIA: highest priority critically important antimicrobial.
Table
Table 3. Concordance between phenotypic and genetic AMR profile identified among the collection.
Table 3. Concordance between phenotypic and genetic AMR profile identified among the collection.
Phenotypic ResistancenGenetic Determinants
Beta-lactams
amipicillin60blaSHV73, blaTEM-1b (2)
blaTEM-1a (4)
blaTEM-1b (38)
blaTEM-1c (8)
blaTEM-1d (4)
blaTEM-214 (1)
blaTEM-220 (3)
ampicillin, ceftiofur, ceftazidime6ampC# (2)
blaCMY-2, blaCTX-M55, blaTEM-1b (1) (1)
blaCTX-M1 (1)
blaCTX-M1, blaTEM-1b (1)
blaCTX-M15, blaTEM-1b (1)
Chloramphenicol23catA1 (8)
catA1,cmlA1 (2)
catA2 (1)
cmlA1 (11)
mdfA, acrAB-TolC (1)
Aminoglycosides
gentamicin1aac(3)_IIa (1)
streptomycin61aadA1 (16)
aadA2 (1)
aadA2b (1)
strA (2)
strB(1)
aadA1, aadA2 (2)
strA, strB (20)
aadA1, aadA2b (6)
strA, strB, aadA5 (2)
aadA1, aadA2b, strB (1)
aadA1, strA, strB (7)
aadA1, aadA2b, strA, strB (2)
gentamicin, streptomycin10aac(3)_IIa, aadA1, aadA2b, strA, strB (1)
aac(3)_IIa, aadA1, strA, strB (1)
ant2_Ia, aadA1 (1)
aac(3)-IId, aadA5, strA, strB (2)
aac(3)_IIa, aac3_IV, aadA1, aadA2b, strA, strB (1)
aac(3)_IV, strA, strB (1)
aac(3)_IV, aadA1, strA, strB (2)
aac)(3)_IIa, strA, strB (1)
Sulphonamides
sulfisoxazole, sulfametoxazole/trimethoprim71sul1 (2)
sul1, dfrA1 (13)
sul1, dfrA12 (1)
sul2 (17)
sul2, dfrA1 (3)
sul2, dfrA14 (6)
sul2, dfrA17 (1)
sul1, sul2 (1)
sul1, sul2, dfrA1 (7)
sul1, sul2, dfrA17 (3)
sul1, sul2, dfrA7 (1)
sul1, sul2, dfrA1, dfrA14 (1)
sul1, sul2, dfrA1, dfrA7 (1)
sul3 (2)
sul2, sul3, dfrA14 (2)
sul3, dfrA1 (3)
sul3, dfrA12 (2)
sul3, dfrA14 (1)
sul1, sul2, sul3, dfrA1 (2)
sul1, sul3, dfrA1 (2)
Sulfamethoxazole/trimethoprim2dfrA1 (1)
dfrA12 (1)
Tetracycline76tetA (54)
tetB (19)
tetA, tetM (1)
tetA, tetB (2)
Colistin4mcr1 (4)
1 *pmrB# (1)
ESBL profile4blaCTX-M1 (2)
blaCTX-M15 (1)
blaCTX-M55 (1)
(Fluoro)quinolones
nalidixic acid4gyrA# (D87G) (3)
gyrA# (A84P), parC # (S57T) (1)
nalidixic acid, enrofloxacin45gyrA# (S83L) (21)
gyrA# (S83L-D87N), parC # (S80I) (13)
gyrA# (S83L-D87N), parC # (S80I-E84G) (2)
gyrA# (S83L-D87N), parC # (S80I-E84A) (1)
gyrA# (S83L-D87N), parC # (S80I-E84G-A56T) (1)
gyrA# (S83L-D87N), parC # (S80I), pare # (S458A) (3)
gyrA# (S83L-D87N), parC # (S80I), pare # (L416F) (1)
gyrA# (S83L-D87Y), parC # (S80I) (1)
gyrA# (S83L-D87Y), parC # (S80I), parE# (S458A) (1)
qnrB19, pare# (I355T) (1)
pare# (I355T) (5)
4 *parC# (A56T)(1)
qnrS1 (3)
#: chromosomal mutation; *: presence of genetic virulence determinant not associated to the expected phenotypic resistance.
Table
Table 4. Representation of phylogroup distribution, Shannon Index (H′) and Simposon index (D) among the sources investigated.
Table 4. Representation of phylogroup distribution, Shannon Index (H′) and Simposon index (D) among the sources investigated.
SourcenAB1B2CDEFunknownH′D
Dairy258 (32%)9 (36%)2 (8%)4 (16%)1 (4%)1 (4%)0 (0%)0 (0%)1.4850.763
Beef240 (0%)22 (91.67%)0 (0%)1 (4.17%)1 (4.17%)0 (0%)0 (0%)0 (0%)0.3450.163
Wild boar220(0%)3 (13.64%)6 (27.27%)0(0%)0 (0%)11 (50%)0 (0%)2 (9.09%)0.9730.68
Vegetable241(4.17%)14 (58.33%)1 (4.17%)6 (25%)0 (0%)1 (4.17%)0 (0%)1 (4.17%)1.1910.616
Fishery2410 (41.67%)4 (16.67%)2 (8.33%)5 (20.83%)0 (0%)3 (12.50%)0 (0%)0 (0%)1.4570.764
Companion animal120 (0%)6 (50%)2 (16.67%)2 (16.67%)0 (0%)2 (16.67%)0 (0%)0 (0%)1.0890.644
Swine258 (32%)8 (32%)0 (0%)8 (32%)0 (0%)1 (4%)0 (0%)0 (0%)1.2230.720
Poultry252 (8%)6 (24%)1 (4%)8 (32%)0 (0%)2 (8%)5 (20%)1 (4%)1.6910.817
Rabbit230 (0%)20 (86.96%)1 (4.35%)0 (0%)2 (8.70%)0 (0%)0 (0%)0 (0%)0.4700.245
Mollusc252 (8%)15 (60%)0 (0%)4 (16%)3 (12%)1 (4%)0 (0%)0 (0%)1.1850.617
Human256 (24%)2 (8%)2 (8%)3 (12%)2 (8%)3 (12%)7 (28%)0 (0%)1.8140.850
Wild animal251 (4%)18 (72%)3 (12%)2 (8%)0 (0%)0 (0%)1 (4%)0 (0%)0.9510.477
Total27938 (13.62%)127 (45.52%)20 (7.17%)43 (15.41%)9 (3.23%)25 (8.96%)13 (4.66%)4 (1.43%)
Table
Table 5. Association between pandemic (ST69/69*, ST95, ST131), emerging (ST10/ST10*, ST23, ST58, ST117, ST405, ST648), uncommon (the remaining ones) ExPEC lineages and functional category profile identified among the collection.
Table 5. Association between pandemic (ST69/69*, ST95, ST131), emerging (ST10/ST10*, ST23, ST58, ST117, ST405, ST648), uncommon (the remaining ones) ExPEC lineages and functional category profile identified among the collection.
LineagesnFunctional Category Profile
Pandemic ExPEC13adhesin, iron acquisiton system (2)
adhesin, iron acquisiton system, protectin (8)
adhesin, toxin, iron acquisition system, protectin, invasin (3)
Emerging ExPEC37adhesin (2)
iron acquisition system (1)
adhesin, protectin (11)
adhesin, iron acquisition system, protectin (18)
adhesin, iron acquisition system, toxin, protectin (5)
Uncommon ExPEC229adhesin (68)
protectin (4)
adhesin, invasin (1)
adhesin, toxin (3)
adhesin, protectin (68)
adhesin, iron acquisition system (9)
protectin, iron acquisition system (1)
adhesin, toxin, protectin (2)
adhesin, iron acquisition system, protectin (50)
adhesin, invasin, iron acquisition system, toxin (1)
adhesin, invasin, iron acquisition system, protectin (2)
adhesin, toxin, iron acquisition, system, protectin (12)
adhesin, toxin, iron acquisition system, protectin, invasin (8)
Table
Table 6. Representation of virulence genes identified in the collection and related virulence factors/functional categories. Each VAG is associated to specific ExPEC pathovars, where they have been commonly identified according to Kaper et al., (2004) [1], Sarowska et al., (2019) [124], Kudinha et al., (2012) [160] (bmaE), Tóth et al., (2009) [161] (cdtB), Ostblom et al., (2011) [162] (malX) and Schierack et al., (2008) [125] (gimB).
Table 6. Representation of virulence genes identified in the collection and related virulence factors/functional categories. Each VAG is associated to specific ExPEC pathovars, where they have been commonly identified according to Kaper et al., (2004) [1], Sarowska et al., (2019) [124], Kudinha et al., (2012) [160] (bmaE), Tóth et al., (2009) [161] (cdtB), Ostblom et al., (2011) [162] (malX) and Schierack et al., (2008) [125] (gimB).
Functional CategoryVirulence FactorGenePathotype
AdhesinM-agglutinin subunitbmaEUPEC
type 1 fimbrial adhesinfimHUPEC, NMEC, SEPEC, APEC
iron-regulated-gene-homologue adhesinihaUPEC
pilus associated with pyelonephritispapCUPEC, SEPEC, APEC
pilus associated with pyelonephritispapG_II; papG_IIIUPEC, SEPEC, APEC
S fimbrial adhesinsfaSUPEC, NMEC
bifunctional enterobactin receptor/adhesinihaUPEC, NMEC
temperature sensitive hemagglutinintshUPEC, NMEC, SEPEC, APEC
Inasininvasion of brain endotheliumibeANMEC, SEPEC, APEC
genetic island associated with newborn meningitisgimBNMEC
Iron acquisiton systemiron-responsive elementireAUPEC
catecholate siderophore receptoriroNUPEC, NMEC, SEPEC APEC
ferric yersinia uptakefyuAUPEC, NMEC
iron repressible proteinirp2NMEC
l-lysine 6 monooxigenaseiucDUPEC, APEC
ferric aerobactin receptor precursoriutAUPEC, NMEC
periplasmic iron binding proteinsitAUPEC, APEC
Toxincytolethal distending toxin BcdtBUPEC, NMEC, SEPEC APEC
cytotoxic necrotising factorcnf1UPEC, MNEC, SEPEC
haemolysin EhlyEUPEC
serine protease autotransporterspicUPEC
serine protease autotransporterssatUPEC
uropathogenic specific proteinuspUPEC
vacuolating autotransporter toxinvatUPEC, APEC
Protectinincreased serum survivalissNMEC, SEPEC, APEC
structural component of colicin V operoncvaCNMEC, SEPEC, APEC
group II capsule antigenskpsMT_IINMEC, SEPEC
outer membrane proteinompTUPEC, NMEC
transfer proteintraTNMEC, SEPEC
K1 capsular polysaccharideneuCNMEC, UPEC
Otherpathogenicity-associated island markermalXUPEC, NMEC, SEPEC, APEC
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Massella, E.; Giacometti, F.; Bonilauri, P.; Reid, C.J.; Djordjevic, S.P.; Merialdi, G.; Bacci, C.; Fiorentini, L.; Massi, P.; Bardasi, L.; et al. Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources. Antibiotics 2021, 10, 351. https://doi.org/10.3390/antibiotics10040351

AMA Style

Massella E, Giacometti F, Bonilauri P, Reid CJ, Djordjevic SP, Merialdi G, Bacci C, Fiorentini L, Massi P, Bardasi L, et al. Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources. Antibiotics. 2021; 10(4):351. https://doi.org/10.3390/antibiotics10040351

Chicago/Turabian Style

Massella, Elisa, Federica Giacometti, Paolo Bonilauri, Cameron J. Reid, Steven P. Djordjevic, Giuseppe Merialdi, Cristina Bacci, Laura Fiorentini, Paola Massi, Lia Bardasi, and et al. 2021. "Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources" Antibiotics 10, no. 4: 351. https://doi.org/10.3390/antibiotics10040351

APA Style

Massella, E., Giacometti, F., Bonilauri, P., Reid, C. J., Djordjevic, S. P., Merialdi, G., Bacci, C., Fiorentini, L., Massi, P., Bardasi, L., Rubini, S., Savini, F., Serraino, A., & Piva, S. (2021). Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources. Antibiotics, 10(4), 351. https://doi.org/10.3390/antibiotics10040351

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop