1. Introduction
Infections caused by extended-spectrum β-lactamase (ESBL)-producing Enterobacteriaceae have been linked to their dissemination through both food and the environment [1,2,3]. ESBL-producing bacteria, including E. coli are frequently resistant to multiple antimicrobials including 3rd and 4th generation cephalosporins in addition to quinolones and aminoglycosides [1,2,4]. These bacteria have been designated priority group 1 pathogens for the development of alternative antimicrobial therapies by the World Health Organization [5] as they can cause serious urinary tract [6,7] and bloodstream [8,9,10] infections, occasionally resulting in mortalities [11].
Genetic determinants of ESBL resistance are commonly associated with mobile genetic elements such as plasmids, transposons, and integrative conjugative elements (ICEs), which facilitate horizontal gene transfer among bacteria [4,12,13]. Plasmids have been frequently identified in E. coli [14], while ICEs appear to be far less common than in other pathogenic bacteria such as Staphylococcus aureus and Vibrio cholera [15,16].
Antimicrobials are used in livestock production for both prophylactic and therapeutic purposes, practices that can promote antimicrobial resistance (AMR) [17,18,19]. A number of reports have characterized AMR in isolated segments of the beef production system, often employing traditional methods of AMR surveillance [19,20,21]. Conventional AMR susceptibility tests have been instrumental in detecting and monitoring AMR status and emergence, but has limitations with regard to the extent that information on phylogeny, mobile genetic elements (MGEs), or virulence traits can contribute to epidemiological investigations [22,23].
Whole genome sequencing (WGS) can make significant contributions to AMR surveillance and in-depth analyses of ESBL-producing E. coli (ESBL-EC) of varying origin so as to identify potential “One Health” linkages [24,25]. This approach can help elucidate the genetic relationship among isolates [23] and to identify their origin within the food production chain [26,27]. Moreover, increasing accessibility and decreasing cost of high throughput WGS has made it the method of choice for routine surveillance and outbreak investigation [28].
An integrated One Health approach is needed to consider animal, food, environment, and human contributions to AMR ecology [29,30]. However, there is a dearth of information about the contemporaneous occurrence of AMR along the One Health continuum [31]. Considering the impact of ESBL infections and the prioritization of these bacterial agents, the present study specifically investigated presumptive ESBL-EC from a One Health continuum linked to the beef production system in southern Alberta, Canada. The objective was to characterize presumptive ESBL-EC with respect to phenotypic and genotypic AMR; phylogeny; and the prevalence of plasmids, ICE and virulence factors among isolates. We hypothesized that, although presumptive ESBL-EC from human and cattle sources along the One Health continuum may harbor some similar ARGs, MGEs, and virulence factors, the accessory genomes of isolates from beef cattle production environments and human sources would be genetically distinct.
2. Materials and Methods 2.1. Study Area, Sampled Sites, and E. coli Isolation
Detailed description of feedlot cattle feces, catch basin, surface streams as well as processing plant, waste water treatment plant, and clinical sources from which presumptive ESBL-EC were collected are published elsewhere [32,33,34]. Cattle feces (CFeces) and liquid samples from catch basins (CBasins) were collected from four feedlots (A, B, C, and D) monthly over 2 years (March 2014 – April 2016). Sampling of feces from cattle was conducted according to the protocol approved by the University of Calgary Animal Care Committee (AC14–0029). In feedlot C, liquid samples were also collected from a surrounding stream (SStreams). Sewage samples were collected from treatment plants located upstream (Calgary) and downstream (Medicine Hat) of the major feedlot cattle production region in Alberta, while human clinical isolates were randomly selected from anonymous individual patients for the same time period from the Division of Medical Microbiology, Calgary Laboratory Services (CLS) biorepository. This laboratory services about 1.5 million people in the City of Calgary and surrounding rural area. Beef processing plant (BProcessing) samples were obtained from ground beef and retail meat and from swabs from hides, beef trimmings, washed carcasses, and conveyers following the procedures of Gills and Jones [35].
Samples from CFeces and BProcessing were individually enriched overnight at 37 °C in E. coli broth supplemented with 2 µg/mL cefotaxime (MilliporeSigma, Burlington, MA, USA) to increase the likelihood of isolating presumptive ESBL from cattle and their environments, while enrichment was not necessary for human sources. Presumptive ESBL-EC from each feedlot sample were isolated by adding 0.5 g cattle feces to 4.5 mL of E. coli broth containing 2 µg/mL cefotaxime) for overnight enrichment at 37 °C (Cameron-Veas et al., 2015; Ibrahim et al., 2016). Each resultant culture was mixed, and a saturated swab was used to inoculate MacConkey plates supplemented with 1 µg/mL of ceftriaxone (MilliporeSigma) and incubated overnight. A maximum of three distinct red/magenta lactose-fermenting colonies were randomly isolated and subcultured onto tryptic soy agar (TSA)-ampicillin (32 µg/mL) (MilliporeSigma) for purification and subsequent archiving at −80 °C. Presumptive ESBL-EC from CBasins and SStreams were isolated via membrane filtration onto MacConky plates with cefatriaxone (1 µg/mL) using the US EPA Method 1603 [32,33]. The same procedure without enrichment or the use of selective plates was used to isolate generic E. coli. All isolates were archived in Brain Heart Infusion (BHI) broth containing 15% glycerol and stored at −80˚C. E. coli isolated from selective plates with or without enrichment were described as presumptive ESBL-EC and were the only isolates selected for whole genome sequencing.
2.2. Antimicrobial Susceptibility Tests
Briefly, resistance of presumptive ESBL-EC (n = 705) to ampicillin, amoxicillin/clavulanic acid, ceftiofor, ceftazadime, streptomycin, neomycin, ceftazidime/clavulanic acid, oxytetracycline, florfenicol, sulfisoxazole, and trimethoprim/sulfamethoxazole was evaluated using the Kirby–Bauer disk diffusion method according to the Clinical and Laboratory Standards Institute [36]. Zone diameters of inhibition were documented with BIOMIC V3 Microbiology System (Giles Scientific, Santa Barbara, CA, USA) and resistant, intermediate, and susceptible frequencies computed as described elsewhere [36]. On the basis of ESBL phenotype, AMR phenotype, isolate source, and sampling site, a representative subset of presumptive ESBL-EC isolates (n = 162) was selected for WGS from (i) beef production sources: CFeces n = 40; CBasin, n = 42; SStreams n = 21; BProcessing n = 4, and from (ii) humans sources: sewage MSewage n = 30 and clinical CHumans n = 25.
2.3. Whole Genome Sequencing and Assembly
From glycerol stocks, presumptive ESBL-EC were cultured on sheep blood agar overnight at 37 °C and a single colony was selected and subcultured on the same media. A single E. coli colony was selected and suspended (~2 × 109 cells/mL) in 10 mM Tris and 1 mM EDTA (TE, pH 8.0) at an OD600 of 2. From each suspension, 1 mL was dispensed into a microcentrifuge tube and centrifuged at 14,000× g for 2 min to form a cell pellet for DNA extraction. Genomic DNA were extracted using the DNeasy Blood and Tissue kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. DNA quality was assessed using Nanodrop 2000 (NanoDrop Technologies, Inc., Wilmington, DE, USA) and was considered acceptable with an A260/280 ~1.8–2.0 [37,38]. Genomic libraries were constructed using the Illumina Nextera XT DNA sample preparation kit with Nextera XT index kit for multiplexing according to the manufacturer’s instructions (Illumina Inc., San Diego, CA, USA). Libraries were then sequenced on the MiSeq Whole Genome sequencing platform (2 × 300 bp paired-end, Illumina Inc.) to obtain an average coverage depth of >100× and a read retention of >98%. The adaptor sequences were cropped with Trimmomatic 0.36 [39] using the following conditions: phred33, LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36. Sequences were assembled using SPAdes 3.10.1 [40] at default parameters with modifications: K-mer values (21,33,55,77,99,127), careful correction, length for average coverage calculation (5000), coverage cutoff ratio (0.33), repeat cutoff ratio (1.75), and contig length cutoff (1000). Sequences were annotated with PROKKA [41] at default settings with modifications: BLASTP e-value cutoff for annotation (1e-06), minimum contig size (200), locus tag prefix (EC), locus tag counter increment (1), GFF version (3), searching for ncRNAs with Infernal+Rfam (false) enabled, strain name (blank), species (coli), genus (Escherichia), kingdom (bacteria), genetic code (11), and plasmid (blank). Downstream WGS analyses were implemented using appropriate programs within a Galaxy pipeline (0.3.0) [42] for identification of resistance genes, plasmids, and virulence genes as detailed in Section 2.4.
2.4. Whole Genome Sequence Analyses
Pan genome analyses was conducted using Roary [43], followed by AMR mechanism identification using the National Center for Biotechnology Information (NCBI) and ResFinder databases [44]. ARGs were designated as present if sequence identity was ≥90%. Agreement between ARGs and AMR phenotypes were examined as described by McDermont et al. [24], where resistance or susceptibility to a given antimicrobial was compared with the presence or absence of a known corresponding resistance gene(s) and/or structural gene mutation(s). WGS sensitivity was computed as the number of isolates that harbored resistant determinants divided by the total number of isolates exhibiting clinical resistance, with the phenotypic resistance as the standard outcome. Specificity was calculated as the number of isolates that did not harbor genetic determinants divided by the total number of phenotypic susceptible isolates [24,45]. Separate binary logistic regression [46] was used to examine a subset of cattle versus human isolates (CFeces and CBasins vs. CHumans and MSewage) for differences in ARG prevalence. The few isolates obtained from the BProcessing environment did not allow for them to be included in more extensive data analyses.
Based on single nucleotide polymorphisms (SNPs) in the core genes, phylotyping was performed using SNVPhyl [47] to generate a maximum likelihood tree with the NCBI reference strain E. coli str. K-12 substr. MG1655. The maximum likelihood tree was generated for the entire 162 presumptive ESBL E. coli isolates as well as a subset of 108 isolates which harbored blaTEM, blaCTX-M, blaOXA, or blaSHV and were designated as true ESBLs. The reference strain E. coli MG1655 was used as it has been widely used by others [48]. The tree was subsequently visualized using iTOL [49]. Whole genome multi-locus sequence typing (wgMLST) and a minimum span tree was also generated using BioNumerics (BN v7.6) (BioNumerics, Applied Maths, Keistraat, Belgium) with default parameters for categorical data. Multi-locus sequence typing (MLST) was achieved with PubMLST E. coli#1 (Achtman) scheme [50], while virulence factors were identified from the Virulence Factors Data Base (VFDB) [51]. To investigate the prevalence of ICEs in E. coli, BLAST homology searches were conducted against 1032 ICE sequences obtained from ICEberg database 2.0 [15]. In silico E. coli subtyping was performed based on outer membrane lipopolysaccharide (O) and flagellar (H) surface antigens using an in-house ECTyper and EcOH [52,53] to designate isolates as either O or H groups.
A set of presumptive ESBL-ECs randomly selected from each source of the beef production system that identified as non-ESBL-EC by WGS (total n = 23 not carrying any blaTEM, blaCTX-M, blaOXA, and blaSHV) and those identified as true ESBL-EC (n = 24 carrying blaTEM, blaCTX-M, blaOXA, and blaSHV) were compared using a GView BLAST atlas [54] to elucidate possible similarities or differences in genomes of the two groups using E. coli MG1655 as reference. The draft genome sequences of the 162 E. coli isolates have been deposited to GenBank under Bio project PRJNA556083.
3. Results 3.1. Occurrence of Phenotypic Antimicrobial Resistance
Of the phenotypic resistances tested, oxytetracycline and ampicillin resistance was most common, while neomycin resistance was the least (Figure 1; Table S1). Although the CLSI minimum inhibitory concentration (MIC) ampicillin break-point was used for isolating presumptive ESBL-EC, not all E. coli were resistant to ampicillin based on subsequent disk diffusion susceptibility tests. Differences may have resulted from resistance gene silencing, loss of plasmids, or perhaps differences between MIC and disk diffusion assays. Among the isolates collected, 1.9% were not resistant to any of the antimicrobials, 5.6% were resistant to only one class, while 69.1% were resistant to five or more classes. Overall, a total of 92.6% of E. coli strains showed MDR in terms of resistance to three or more classes of antimicrobials. It is worthy to note that, despite enrichment, a limited number of presumptive ESBL-EC isolates were obtained from BProcessing environments (n = 4), an attestation to the sanitary procedures within the plant.
3.2. Antimicrobial Resistance Genes Prevalence
Tetracyclines: A total of 77 unique ARGs encoding diverse antimicrobial resistance were detected. The tetracycline determinants tet(A) and tet(B) were detected in both cattle and human sources, whereas tet(C) and tet(M) were detected in cattle sources, but not MSewage nor CHumans. In contrast, tet(32) was exclusive to isolates associated with human sources (Figure 2; Table S2).
Beta-lactams: β-lactam blaEC-8, blaEC-15, blaTEM-1, blaCTX-M-15, blaCTX-M-55, blaEC-18, and blaCMY-2 were present in all CBasins, CFeces, SStreams, CHumans, and MSewage isolates, whereas blaTEM-150 was detected only in cattle-sources (CBasins, CFeces, and SStreams). In contrast, blaCTX-M-1, blaCTX-M-14, and blaOXA-1 were detected only in CHumans and MSewage, and blaCTX-M-27 and blaCMY-42 were exclusive to human clinical isolates. Additionally, blaLAP-2 occurred in isolates from CBasins and CFeces, with blaSHV-11 and blaSHV-12 occurring exclusively in isolates from CFeces. Folate pathway inhibitors: For folate synthesis inhibitors, sul2 was ubiquitous but sul1 was not detected in BProcessing isolates, whereas sul3 was not found in either BProcessing or SStreams isolates. dfrA1 was present in MSewage, CBasins, and CFeces, while dfrA12 was unique to CBasins, CFeces, and SStreams. dfrA7 was unique to CBasins and dfrA18, and dfrA19 was unique to CHumans isolates. DfrA14, dfrA17, and dfrA27 were present in all CBasins, CFeces, SStreams, CHumans, and MSewage isolates. Aminoglycosides: Aminoglycoside resistance genes also revealed the ubiquity of aph(3’’)-Ib and aph(6)-Id. The genes aph(3’)-Ia and aadA1 were identified in CHumans, MSewage, CBasins, and CFeces isolates, whereas the aph(3’)-IIa and aph(6)-Ic were only present in CFeces and MSewage isolates. Aac(3)-IIa, aac(3)-IId, and Sat1(A7J11) occurred in CHumans and MSewage, whereas the latter two determinants were also identified in CFeces and CBasins isolates. Isolates from MSewage, CBasins, and CFeces possessed aadA16 and armA, and aac(3)-Via and aadA22, respectively. AadA2 were detected in CBasins and CFeces, while aac(6’)Ib-cr was only found in isolates originating from humans or sewage. Phenicols: Phenicol resistance determinants catB3 and catB4 were only found in CHumans and MSewage isolates, whereas cmlA1 and cmlA5 were unique to CBasin and CFeces isolates. In contrast, floR was found in isolates collected from all environments.
Quinolones: The plasmid-borne quinolone-resistant gene qnrS1 was also ubiquitous, while qnrB2 was only found in CHumans isolates (Figure 2; Table S2). Aside from plasmid-borne quinolone resistant qnr and aac(6’)Ib-cr genes, evaluation for mutations in the quinolone resistance-determining regions (QRDRs) of gyrA, gyrB, parC, and parE aligned with frequently reported mutations (gyrA: S83L, S83A, D87N, and D87Y; gyrB: E185D; parC: S80I and E84G; parE: S458T and I529L [55,56]), but additional mutations were found that may also have also conferred quinolone resistance (Table S4). All BProcessing isolates harbored only tet(C), sul2, blaEC-18, blaCMY-2, aph(3’’)-Ib, aph(6)-Id, and floR genes.
Overall: In cattle- and human-derived isolates, β-lactam (100% in both sources) and tetracycline genes (91.6% and 87.3%) were prevalent, with quinolone genes being the least prevalent (20.6% and 10.9%) (Figure S1 and S2; Table S2). When cattle (CFeces and CBasins) and human (CHumans and MSewage) isolates were compared, differences in prevalence were observed. The prevalence of 13 ARGs (sul1, sul2, blaEC, blaEC-8, blaEC-15, blaEC-18, blaCTX-M-15, blaCMY-2, aph(3’’)-Ib, aph(6)-Id, floR, blaTEM-1, and blaCTX-M-55) differed (p < 0.05; Table S5) between isolates originating from cattle vs. human sources and were generally more prevalent in human isolates, with the exception of blaTEM-1 and blaCTX-M-55.
Misc: Determinants for resistance to rifamycin (arr-2 and arr-3), lincosamides (lnu(F)), bleomycin (ble), streptothricin (sat1), and quaternary ammonium compounds (QAC), efflux small multidrug resistance transporter genes (qacG2, qacEΔ1 and qacL) were also detected. Rifamycin resistance was found in CFeces (7.5%) and MSewage (10%) isolates; bleomycin was found only in to MSewage (3.3%) isolates; and streptothricin in CBasins (2.9%), CHumans (4%), and MSewage (10%) isolates, while lincosamide determinants were present only in CFeces (1.9%) and MSewage (3.3) isolates. The QAC efflux genes were less frequent in CFeces (17.5%), CBasin (19.0%), and SStreams (28.6%) than in MSewage (50%) and CHuman (48%) isolates (Figure 2; Figure S1). Categorically mph(E), msr(E), erm(B), arr-3, and ble were only found in human-sourced isolates, while qacG2 were exclusively associated with cattle-sourced isolates.
3.3. Pan-Genome Analysis and Compairson of AMR Determinants as Identified by WGS to the Occurrence of Phenotypic Reisistance
E. coli genomes ranged from 4.6–5.7 Mbp (Table S3). Analyses of the 162 isolates identified a total of 24,338 genes, with 2763 core genes in all isolates, 327 soft core genes present in 95%–99% of isolates, 2789 shell genes present in 15%–95% of isolates, and 18,459 cloud genes present in less than 15% of isolates.
Sensitivity for AMR prediction by WGS compared to clinical resistance phenotype varied per antimicrobial class and antimicrobial agent, with ranges between 93.2% (ampicillin) and 100% (neomycin, trimethoprim/sulfamethoxazole, and enrofloxacin), while specificity ranged between 1.1% (enrofloxacin) to 100% (ampicillin) (Table 1; Table S6). WGS-positive predictive values also ranged between 7.7% (neomycin) and 100% (ampicillin), while negative predictive values ranged between 61.5% (ampicillin) and 100% (neomycin, trimethoprim/sulfamethoxazole, and enrofloxacin). WGS revealed that 90.7% of all the presumptive ESBL-EC were true ESBL-producers based on the presence of blaSHV, blaTEM, blaOXA, blaCTXM, or blaEC-13–15–18–19, with gene prevalence not differing by source (p = 0.89). In retrospect, both PCR and WGS identified comparable levels of true ESBL-EC 68.5% and 66.7% respectively, based on the presence of blaSHV, blaTEM, blaOXA, or blaCTXM genes.
3.4. Phylogenetics, MLST Cluster, and Serogroups Analysis
Evaluation of the genetic relatedness of the beef/cattle-associated and human-associated isolates using a maximum likelihood phylogenetic tree showed that 71.64% of the core genome of isolates from these sources was shared. Limited phylogenetic intermixing among isolates from segments along the One Health continuum was observed (Figure 3). Single nucleotide variant phylogenomics revealed that, in general, human-sourced ESBL-EC comprising CHumans and MSewage formed one major cluster within which a number of distinct clades and sub-clades occurred. The majority of CHumans and MSewage isolates were closely related, although about 7.5% of all beef/cattle-associated isolates (n = 107) were in this cluster. Interestingly, 92.5% of the isolates from CFeces, CBasins, and SStreams formed a separate cluster of mixed clades with only a few MSewage (n = 9) and CHumans (n = 2) isolates in this clade. Isolates from CBasins and CFeces often clustered together. Although the few BProcessing isolates formed a distinct sub-clade from CBasins and SStreams, they also clustered separately from isolates obtained from other points within the beef production chain. Apart from an isolate from CFeces, that was within 18 SNP of MSewage isolates (=2) (Figure S3), the lowest SNP pairwise differences between CHumans and MSewage vs. all cattle-sourced isolates were 552 and 718, respectively. Minimum span tree (MST) generated with Bionumerics and whole genome MLST (wgMLST) corroborated our phylogenetic analysis and showed distinct clusters of cattle- versus human-sourced isolates (Figure 4; Figure S4). Subsequently, a separate phylogenetic tree which featured only true ESBL isolates also showed limited intermixing between cattle and human isolates (Figure S5). BLAST atlas of isolates identified as non-ESBL versus true ESBL-EC by WGS did not reveal much diversity among non-ESBL and ESBL-EC sequences (Figure 5), although the genetic elements that appeared to be commonly absent in both E. coli groups were associated with prophage elements when compared to the reference genome (data not shown).
Using MLST, 58 sequence types were identified for all but five isolates, which did not match any known ST groups. CFeces isolates generally belonged to ST groups 224 (37.5%), 10 (20%), 23 (5%), and unknown ST (5%) (Figure 6; Table S7). CBasins isolates were mainly identified as members of ST10 (23.8%), 224 (9.5%), and 515 (9.5%). For SStreams isolates, ST69 (14.3%) and ST155 (14.3%) were the most common, followed by STs 244, 515, and 7618 (all at 9.5%). All four isolates from BProcessing belonged to ST20. The pandemic ST group 131 exclusively occurred in MSewage and CHumans, accounting for the largest portion of the isolates at 26.7% and 52%, respectively. In MSewage isolates, other common ST groups included 38, 44, 405, and 1193, while other common ST in CHumans isolates were ST 648 and 38.
In silico typing using the ECTyper and EcOH databases yielded 67/162 and 37/162 unknown and novel O-antigen loci, respectively. About half of the clinical isolates (52%; n = 13) were O25:ST131, while three MSewage isolates were of the wzx-Onovel31:ST131 subtype. The majority of cattle-sourced isolates were found to be O8, O89, O9, O15, and O99. Two BProcessing O group isolates remained unclassified, while the other two isolates both classified as O128:H2. All presumptive ESBL-EC identified by MLST as ST224 (n = 21) showed an unknown O-antigen with a H23 flagella antigen. 3.5. Prevalence of Plasmids and Integrative Conjugative Elements
Most of the presumptive ESBL-EC (98.8%) harbored at least one of the 38 different incompatibility (Inc) and colicinogenic (Col) plasmid types detected (Figure 7; Figure S6). The single most predominant plasmid type across sources was the IncFIB(AP001918) occurring in CFeces (27.5%), CBasin (40.5%), SStreams (19.0%), BProcessing (25%), MSewage (86.7%), and CHumans (64%). Other frequently detected plasmid types along the entire continuum included the IncFII (47.9%), ColRNAI (25.2%), IncFIA (21.5%), and p0111 (19.6%). Prevalence of plasmids Col(MG828), ColRNAI, Col156, IncFIA, IncFIB(AP001918), IncFII, and IncFIC(FII) was higher in human-sourced (CHumans and MSewage) than cattle-sourced isolates (CFeces plus CBasins), while p0111 was higher in cattle isolates (p ≤ 0.018; Tables S8 and S9). CFeces and MSewage isolates showed the highest plasmid diversity (24/38 each), followed by SStreams (16/38). Plasmids Col(BS512), Col(KPHS6), Col8282, Col(MGD2), IncFII(29), IncFII(pCoo), and IncI2 were unique to human-associated E. coli, while plasmids Col(IRGK), ColE10, IncFIA(HI1), IncFII(pHN7A8), IncHI2, IncHI2A, and TrfA were unique to cattle-associated isolates.
The ICEberg database used to describe ICE revealed that all 162 E. coli isolates harbored these elements. In fact, twenty-two (7.9%) of all the ICE reported in this study were found in all 162 genomes (Table 2; Table S9). Most of these elements (ICEVchBan8, ICEValHN437, ICEValHN396, ICEValA056–2, and ICEVflInd1) were members of the SXT/R391 ICE family originally found in Vibrio spp. Additionally, a proportion of these ubiquitous ICE have been previously reported in Pseudomonas spp. (ICE6441, ICE6440, PFGI-1, and ICEPaeLESB58–1), Salmonella spp. (ICESenTy2–1, ICESb2, and ICESb1), and Yersinia spp. (YAPI and ICEYe1). Notably, ICEEc2 was the only element that was originally described in E. coli. Many of the elements were of unknown function, with a few associated with toxin-antitoxin systems, pathogenicity, and mercuric or antimicrobial resistance.
3.6. Prevalence of Virulence Genes
A total of 207 different virulence determinants were identified (Table S11). Some of the most frequently encountered (70–100%) virulence genes in presumptive ESBL-EC were associated with the components of the type II (T2SS) and III secretion system (T3SS), ferric enterobactin esterase, ferric enterobactin transport ATP-binding protein, enterobactin biosynthesis, adhesion, fimbral adhesins, E. coli common pilus (ECP), and invasion-related outer membrane protein.
Other virulence genes which play various roles in E. coli pathogenicity, including the pap fimbral adhesins that were associated with the CFeces, MSewage, and CHuman isolates, with higher frequencies in human-sourced E. coli. Fimbrial adhesins f17 variants and cdt were identified in only CFeces and CBasins. The pic (8%), sat (52%), senB (28%), set (8%), sfa (4%–12%), and tpc (8%) genes were also found in CHumans, but not in cattle-sourced isolates (Table S11). Hemolysin genes (hlyA, hlyB, hlyC, and hlyD) were more often associated with CHumans (32%) isolates, but were also present in most other sources [CFeces (2.5%), CBasins (4.8%), SStreams (4.8%), MSewage (6.7%)]).
Overall, 100% of the isolates carried entB, entC, entF, entS, fepA, fepB, fepC, and fes and 99.4% carried entA, fepD, fepG, and ompA. Only five genes (espL1, espR1, espX1, espX4, and espX5) were more prevalent (p < 0.0001) in cattle (comprising CBasins and CFeces isolates) than human sources (comprising CHuman and MSewage isolates), while 46 virulence determinants including aslA, fimA, fyuA, gspL, irp12, iucABCD, kpsMD, shuATX, ybtAEPQSTUX, chuSTUVWY, espY1Y2Y3, papBCDFGHIJKX, and hlyABCD were dominant in human-sourced isolates (p ≤ 0.03; Table S12).
4. Discussion 4.1. Antimicrobial Resistance Determinants
The overall ubiquity of β-lactam and tetracycline genes among cattle- and human-sourced isolates across the One Health continuum was expected as both physicians and veterinarians use these classes of antimicrobials to treat humans and beef cattle, respectively [17,57,58]. Tetracycline resistance genes, especially the efflux pumps encoded tet(A), have been detected in high frequencies in generic E. coli and ESBL-EC from bovines and humans from United States and China [58,59,60] and from Canadian feedlot environments using metagenomics [61]. However, the ribosomal protection mechanisms conferred by tet(M) exclusive to cattle isolates has rarely been reported in fecal samples from beef cattle. Metagenomics studies of manure, soil, and wastewater from dairy and beef production systems [61] within the same region in Canada did not identify tet(M). Since the present study investigated only presumptive ESBL-EC, it is possible that tet(M) may be present in the selected ESBL-EC subpopulation of the feedlot environment, likely acquired through horizontal gene transfer (HGT) as proposed by Bryan et al. [58,62]. The predominance of sul2, aph(3’’)-Ib, aph(6)-Id, and floR has also been previously observed in cattle- and human-sources [63,64,65]. These resistances seem to arise with AMU and likely reflect the use of aminoglycosides, sulfonamides, quinolones, and/or phenicols for therapy in humans and cattle.
The exclusivity of bleomycin anticancer antimicrobial inactivator ble [66] in MSewage ESBL-EC likely reflect the sole use of this antimicrobial in humans and not cattle, as sewage would receive human excrement [67,68]. Although this appears to be the first report of the ble gene in municipal sewage in Canada, emerging carbapenemase NDM-1 has been associated with bleomycin resistance in clinical Enterobacteriaceae and attributed to selective pressure of bleomycin or bleomycin-like molecules [69].
Quaternary ammonia compounds are often used as disinfectants in processing plants and clinical settings. A number of qac subgroup members of small multidrug proteins have been previously identified in Enterobactericaeae including qacE, qacEΔ1, qacF, qacG, and qacH. These genes are usually associated with plasmid-mediated class 1 integrons, which carry a variety of ARGs [70,71]. The qac genes are frequently found in combination with genes coding for β-lactams, aminoglycoside, sulfonamide, chloramphenicol, and trimethoprim resistance [4,72]. The regular association of qac genes with human isolates in this study is probably due to frequent use of hydrogen peroxide, quaternary ammonium compounds, and sodium hypochlorite-based chemical disinfectants in clinical and household settings [73,74]. Considering that only four isolates from BProcessing were evaluated in this study, the likelihood of detecting qac genes was limited even though these disinfectants are often used in meat processing plants. Overall, the generally low occurrence of QACs in CFeces, CBasins, and SStreams in contrast to occurrence in human isolates may reflect the rare use of these sanitizers in feedlot environments.
Enrichment was required to isolate ESBL-EC from feedlot-related samples but not for samples from human sources. It may be deduced that ESBL-EC in human-sourced samples were levels of magnitude higher, possibly selected on the basis of clinical illness and/or resistant infections. Variants of blaCTX-M-1,-14,-27 were only associated with ESBL isolates from humans, a finding which is compatible with other reports of blaCTX-M-14,-15,-27 association with the global spread of blaCTX-M in human clinical isolates [1,75,76]. Although blaCTX-M-15 was more frequently detected in human- than cattle-related isolates, their occurrence in both cattle- and human-sources may possibly imply transmission between humans and animals [75,77] as ESBL-EC can be zoonotic [78]. In this study, blaSHV genes were rather rare in contrast to other Canadian investigations where blaSHV2 variants have been frequently reported [75,79]. An interesting finding of the current work is that both blaSHV-11 and 12 were detected in a single isolate from CFeces. E. coli blaSHV-12 was reported in one of the first clinical ESBL-EC cases in a dog with recurrent urinary tract infection in Spain in 2000 [80,81]. These findings implicate that it is possible that even minimal interactions may cause resistance gene transfer within a One Health continuum [81].
We also report what is believed to be the first detection of the recently described blaLAP-2 class A β-lactamase gene [82] from E. coli in Canada. The LAP-2 gene (originally described in 2007) encodes narrow spectrum β-lactamase resistance. All the isolates from CFeces (n = 2) and CBasins (n = 1) with LAP-2 also harbored qnrS1, blaTEM-1, and dfrA14, together with several other ARGs. Co-carriage of LAP-2 and plasmid mediated qnrS1 has been previously reported, and LAP has been found to be associated with the same gene cluster together with insertion elements and transposons in clinical E. coli [83] and Enterobacter cloacae [82,84]. The present study found that isolates carrying LAPwere resistant to 8–10 antimicrobial agents and carried at least five distinct plasmid types, implying a high potential for horizontal gene transfer of MDR.
4.2. Genotype-Phenotype Antimicrobial Resistance Concordance
There was a strong agreement (≥93.2%) between clinical resistance phenotypes and genetic determinants across the One Health continuum for all antimicrobials tested. Notably, sensitivity and specificity for quinolone (based on qnr and aac(6’)Ib-cr genes) were initially low (data not shown). Upon examination of QRDR mutations, sensitivity catapulted to 100%, with each quinolone resistant isolate harboring a resistant gene and/or mutation. In contrast, WGS specificity reduced to 1.1%, with a large proportion of quinolone susceptible isolates also possessing a gene and/or point mutations. These differences between WGS sensitivity and specificity based only on resistance genes versus a combination of genes and mutations may be attributed to the possibility that not all QRDR mutations confer clinical resistance to quinolones or have yet to be confirmed as causing quinolone resistance. A recent study by Varughese et al. [85] found that the commonly reported S83L QRDR mutations did not yield quinolone resistance; although Varughese et al. [85] also noted that coupled with other mutations, S83L conferred quinolone resistance in uropathogenic E. coli.
To assess and facilitate adoption of WGS as the gold standard for AMR surveillance, a number of studies investigating the sensitivity and specificity of WGS in E. coli [86,87] have reported high agreement between a variety of gene determinants and associated phenotypes. In this study, tetracycline, ampicillin, streptomycin, trimethoprim/sulfamethoxazole, and florfenicol resistance exhibited a high concordance for both sensitivity and specificity. An important aspect of WGS for gene detection is the utilization of already identified sequences as references. It is possible that discrepancies between WGS and phenotypic resistance may result from the presence and expression of novel genes which encode for associated resistance phenotypes. Indeed, previous studies have identified phenotypic resistance which were unaccounted for by known resistance genes (e.g., chloramphenicol, gentamicin, streptomycin, and cefoxitin [24]). It is also possible that WGS sensitivity was impacted by the classification of intermediate resistant and sensitive isolates as resistance phenotype negative. In this study, an appreciable proportion of intermediate ceftazidime resistant (37 out 40 intermediate) isolates were found to harbor genes associated with ceftazidime resistance. McDermott et al. [24] made a similar observation, where intermediate resistant isolates were found to harbor resistance genes. The absence of phenotypic resistance in the presence of resistance genes could also be explained by differences in gene expression among isolates [88].
4.3. Insights from E. coli Phylogenetics, MLST and Serogroups
Investigation of the genetic relatedness of presumptive ESBL-EC along the One Health continuum showed generally distinct clades, with evidence of genomes clustering by origin from cattle or human sources, as observed by others [25]. Genetic distinction among ESBL-EC originating from cattle versus human sources along the continuum via phylogenetics implies that isolates are separately adapted to their respective environments and that exchange of isolates between these niches maybe limited [25,89]. For example, ESBL-EC ST131 is reported to be highly human–host specific as it is frequently linked to health-care system exposures rather than environmental sources [90]. It could also be inferred that the segments of the beef production system investigated might not be directly implicated in recent human ESBL-EC infections in this region [25]. Ludden et al. [25] and Salinas et al. [91] also observed distinct clusters between E. coli isolates from livestock and humans. However, a study by Mulvey et al. [92] highlighted the possibility of linkage of MDR E. coli occurring in cattle and human infections in Canadian hospital on the basis of plasmid blaCMY A/C replicon fingerprint similarities >90%. In this study, comparative WGS discriminated between isolates originating from cattle and farm environments which were phylogenetically related, as were isolates from CHumans and MSewage. The frequent clustering of CFeces and CBasins ESBL-EC suggested that these isolates were closely related, suggesting that CBasins were catching flow from feedlot pens and reducing the dissemination of ESBL-EC to the broader environment. However, since a number of CFeces and SStreams were intermixed, it can be deduced that some ESBL-EC may be potentially released into the broader environment as water from catch basins is often utilized to irrigate surrounding crops.
Interestingly, different STs prevailed in different sources, even between CBasins and SStreams. Similar to other reports, ST244 were dominant [75] in cattle-sources, while hyperendemic ExPEC ST131 (implicated in bloodstream and urinary tract infection) [76] and, to a lesser degree, STs 38 and 648 were predominant in human-sources. E. coli ST648 exhibiting an ESBL phenotype have been reported globally in human patients [7,75]. In concordance with the frequent occurrence of ST131 human sources in this study, ESBL-EC isolated from 11 different Canadian medical centers showed a high occurrence of the clonal complex ST131 producing CTX-M-15, -14, and -27 from Brampton, Calgary, and Winnipeg, illustrating its propensity to cause nosocomial infections in Canada [93]. Additionally, the overall findings of distinct ST occurrence between sources agrees with other reports [25,75], where ST131 was only found in E. coli from dogs and not in farm animals [75]. The unique predominance of specific ST also possibly reflects the spread of clonal lineages or the outcome of selective pressures within the diverse segments of the continuum [94,95].
Serogroups in ESBL-EC from CFeces were mostly unknown O serogroups (50%), with a few O8, O89, and O99, while human-sourced ESBL-EC were mostly O25:H4-ST131. Except for a single instance of O26:H11 and wzx-/wzy-Onovel26:H28, our study did not identify any O157 nor any of the other top 6 non-O157 serogroups (O45, O103, O111, O26, O121, and O145) [96]. O157 and the top 6 non-O157 serogroups cause the majority of Shiga toxigenic E. coli (STEC)-associated foodborne infections [96,97]. The infrequent occurrence of these STEC serogroups in ESBL-EC along the continuum could be due to the unspecific nature of the ESBL-EC isolation procedure employed, unlike the specific procedures that are used to isolate STEC, such as direct PCR [98] and sequential immunomagnetic separation [99]. Notwithstanding, the O128 serogroup identified in BProcessing isolates from hides have been linked to serious human infections [97] and may serve as a risk of product contamination during processing. The majority of the O25:ST131 clones from human clinical isolates were CTX-M 15 producers, similar to findings of Aslantaş and Yilmaz [100], who characterized ESBL-EC from Canadian dogs. The specific differences between serogroups along the continuum aligns with MLST and phylogenetic differences in this study and further suggests little to no transmission between cattle and humans.
4.4. Mobile Elements: Plasmid and Integrative Conjugative Elements
It is thought that livestock and animal-derived foods are potential sources of ESBL-EC in humans because the same ESBL genes or plasmids have been detected in livestock and farm workers [78,101]. In this study, the greatest plasmid overlaps occurring among all sources of the One Health continuum were IncFIB(AP001918) and IncFII. Johnson et al. [102] noted that, regardless of E. coli source, IncFIB was predominant in avian, human, and poultry meat isolates in the US, although plasmid replicons and colicin-related genes differed among E. coli sources. Likewise, Rodriguez-Siek et al. [103] also found differences in the actual prevalence of plasmid, virulence, phylogeny, and serogroup traits between human and avian E. coli, a finding that is compatible with the differences observed between human and cattle isolates in the present study. In Canada, there is little information regarding the occurrence and frequency of plasmid-mediated transfer of ARGs and virulence determinants among E. coli originating from different segments of the One Health continuum. Nevertheless, plasmid-mediated horizontal transfer of multiple ARGs is important for resistance circulation among bacteria in various sources of the continuum [57,92], considering that ~99% of the isolates in our study bore at least one plasmid.
Surprisingly, ICE were found throughout all the genomes sourced from the One Health continuum, even though reports of cattle-related E. coli-sourced ICE are rare [15]. Apart from the characteristic set of core genes encoding functions essential for self-transmission and maintenance, ICE are known to often carry cargo genes that impart various fitness and adaptive advantages to hosts [104,105]. Consequently, besides other well-characterized MGEs including plasmids and integrons, ICE may play a crucial role in the evolution and adaptation of E. coli to various niches, including those environments subject to the selective pressure of antimicrobials. In this study, genes associated with MDR, metal resistance, and pathogenicity functions were noted for some ICE. The presence of diverse ICE families may also imply exchange of genetic material among different bacterial species. For example, the SXT/R391 family of ICE often in Vibrio spp. [106] were present in E. coli genomes. Further characterization of ARG-bearing MGEs and their roles in gene exchange via conjugation would shed more light on the extent and frequency at which these elements transfer AMR and virulence genes among E. coli.
4.5. Occurrence of Virulence Genes
Often, STEC are characterized for their stx virulence determinants due to their role in enterohemorrhagic colitis and hemolytic urea syndrome. The repertoire of virulence factors detected in this study included common genes that encode adhesins essential for biofilm formation, host cell invasion, tissue degradation, and host cell death. For example, the fimbria adhesin ecpABCDE that encodes the E. coli common pilus essential for early-stage biofilm development and host cell recognition [107] were common in most isolates (91.4%). Likewise, the type 1 fimbriae determinants fimBCDEFGHI involved in early stage biofilm formation on host mucosa and abiotic surfaces was also common. The detection of T3SS effector protein genes in >70% of isolates in this study demonstrate that all these isolates likely harbored the E. coli locus for enterocyte effacement (LEE), the most important pathogenicity island for enterohemorrhagic and enteropathogenic E. coli [108].
The exclusive occurrence of f17 fimbriae genes in pathogenic E. coli in CFeces and associated CBasins perhaps reflects point-source dissemination of virulence genes from CFeces to CBasins. The f17 fimbriae mediates binding to host intestinal microvilli and has been linked to diarrhea and septicemia outbreaks in calves and lambs [109]. Additionally, the higher occurrence of esp in cattle-sourced E. coli and the higher occurrence of pap, fim, hly, and iuc in human-sourced isolates probably reflects the importance of these factors in E. coli survival or pathogenicity in their respective hosts. All factors considered, the detection of unique virulence genes specific to niche/environmental hosts may support that virulent ESBL-EC occurring in humans did not likely originate from beef cattle. It may also be deduced that, while some virulence genes are ubiquitous due to their importance for fitness and adaptation to biotic and abiotic conditions, other virulence genes may be unique to isolates based on the source characteristic as well as isolate adaptation [26,110].
5. Conclusions This study used epidemiologically robust study design and sample collection to examine geographically and temporally related presumptive ESBL-EC isolates from beef cattle, water sources, and human clinical samples in southern Alberta. Comparison of the genotypic and phenotypic characteristics of these isolates along a One Health continuum revealed differences in prevalence of similar ARGs, MGEs, and virulence factors together with phylogenetic differences, implying that ESBL-EC originating from beef cattle may not play a significant role in ESBL-EC infections in humans in southern Alberta. The ubiquity of diverse plasmid and ICE families could be indicative of HGT, even within the same niche and warrants further investigation. This study also documents the first detection of the blaLAP-2 β-lactamase in cattle feces and catch basins and of the ble anticancer resistance genes in municipal sewage in Canada. From specificity and sensitivity outcomes, this study noted that WGS promises to be a robust method for AMR investigations, but warrants further validation and improvement in AMR phenotype prediction accuracy. Future investigation of ESBL-EC based on sound epidemiological project designs will continue to add to the body of work described here and may provide useful insights at to the significance of similarities and differences identified within the various segments of the One Health continuum.
Enlarge this image.Figure 1. Phenotypic antimicrobial resistance profiles of presumptive extended-spectrum β-lactamase (ESBL) E. coli obtained from different sources as described in materials and methods. OXYT-tetracycline; AMPI-ampicillin; CTIO-ceftiofur; STEP-streptomycin; SULF-sulfisoxazole; TMSZ-trimethoprim/sulfamethoxazole; FLOR-florfenicol; ENRO-enrofloxacin; CTZD-ceftazidime; AMCL-amoxicillin/clavulanic acid; NMYN-neomycin; MDR-multidrug resistance. MDR referred to resistance to two or more classes of antibiotics.
Enlarge this image.Figure 2. Heat map of resistance genes associated with various classes of antimicrobials as detected by the whole genome sequences of 162 presumptive ESBL-producing E. coli across a One Health continuum: Isolates originated from multiple segments of the One Health continuum including human clinical, municipal sewage, beef processing plant, cattle feces, catch basin, and surrounding streams. aac-(6')-Ib-cr classified under aminoglycoside also encodes quinolone resistance. The overall prevalence of each determinant is displayed as a numeric percentage at the bottom of each determinant.
Enlarge this image.Figure 3. Phylogenetic tree generated on the basis of single-nucleotide polymorphisms (SNPs) of the core genes of 162 E. coli isolates obtained from multiple segments of the beef production system and human-associated isolates with reference genome E. coli str. K-12 substr. MG1655 (GenBank accession # GI: 545778205/U00096.3). Inner, middle ring, and outer rings are representative of isolate sources, phenotypic antimicrobial, and multidrug resistance, respectively.
Enlarge this image.Figure 4. Minimum spanning tree (MST) based on whole-genome multi-locus sequence typing (wgMLST) profiles of 162 presumptive ESBL-producing E. coli genomes along the One Health continuum from cattle feces, catch basins, surface streams, beef processing plants, municipal sewage, and humans: The MST includes a total of 9580 wgMLST loci and was constructed with BioNumerics (version 7.6.2). Each circle corresponds to a unique wgMLST profile and is colored based on sample origin. The size of the circle is proportional to the number of isolates sharing the same wgMLST profile, while the branch lengths correspond to the number of allele differences between isolates. Boxes highlight three major clusters of isolates of cattle and human origin.
Enlarge this image.Figure 5. Blast atlas (A-circular) showing sequences similarity and diversity among true ESBL and presumptive non-ESBL E. coli using E. coli MG1655 as a reference: Isolates shown in this analysis are all cattle-derived (cattle feces, catch basins, surface water, and processing plant) and do not include any human-derived presumptive ESBL-EC because WGS did not identify any non-ESBL-EC from humans or municipal sewage sources based on the presence of blaSHV, blaTEM, blaOXA, or blaCTX-M. Inner ring comprises GC skew, GC content, and sequence of reference E. coli MG1655, while middle and outer rings comprise ESBL and non-ESBL E. coli, respectively.
Enlarge this image.Figure 6. Percentage prevalence of specific sequence types (ST10 to ST7618 and unknown STs) identified in multiple segments of the One Health continuum using the E. coli Achtman scheme.
Enlarge this image.Figure 7. Prevalence of plasmid replicons per source of presumptive ESBL-producing E. coli from cattle feces, catch basins, surface streams, municipal sewage, humans, and beef processing along a One Health continuum.
| Phenotype: Resistant | Phenotype: Susceptible | |||||||
|---|---|---|---|---|---|---|---|---|
| Gene Positive | Gene Negative | Gene Positive | Gene Negative | Sensitivity (%) | Specificity (%) | PPV (%) | NPV (%) | |
| Aminoglycoside | ||||||||
| Streptomycin | 111 | 4 | 16 | 31 | 96.5 | 66.0 | 87.4 | 88.6 |
| Neomycin | 8 | 0 | 96 | 58 | 100.0 | 37.7 | 7.7 | 100.0 |
| Beta-lactam/beta-lactam inhibitor | ||||||||
| Cephems | ||||||||
| Ceftazidime | 61 | 2 | 74 | 25 | 96.8 | 25.3 | 45.2 | 92.6 |
| Ceftiofor | 127 | 3 | 8 | 24 | 97.7 | 75.0 | 94.1 | 88.9 |
| Penicillin | ||||||||
| Ampicillin | 136 | 10 | 0 | 16 | 93.2 | 100.0 | 100.0 | 61.5 |
| Folate pathway inhibitors | ||||||||
| Trimethoprim/sulfamethoxazole | 93 | 0 | 29 | 40 | 100.0 | 58.0 | 76.2 | 100.0 |
| Sulfisoxazole | 112 | 5 | 6 | 39 | 95.7 | 86.7 | 94.9 | 88.6 |
| Phenicol | ||||||||
| Florfenicol | 67 | 3 | 16 | 76 | 95.7 | 82.6 | 80.7 | 96.2 |
| Quinolones | ||||||||
| Enrofloxacin | 70 | 0 | 91 | 1 | 100.0 | 1.1 | 43.5 | 100.0 |
| Tetracyclines | ||||||||
| Oxytetracycline | 144 | 5 | 2 | 11 | 95.7 | 84.6 | 98.6 | 68.8 |
Sensitivity (%) equals number of isolates that harbored resistant determinants divided by the total number of isolates exhibiting clinical resistance phenotypes, while specificity (%) was calculated as the number of isolates that did not harbor genetic determinants divided by the total number of phenotypically susceptible isolates. Whole genome sequence of positive predictive value (PPV%) and negative predictive value (NPV%).
| ICE name | Function | Family | Source Bacteria |
|---|---|---|---|
| ICEVchBan8 | toxin-antitoxin system | SXT/R391 | Vibrio cholerae MZO-3 |
| HAI2 | Unknown (-) | Unclassified | Erwinia carotovora subsp. atroseptica SCRI1043 |
| ICECroICC168–1 | Unknown (-) | Unclassified | Citrobacter rodentium ICC168 |
| ICE6441 | Unknown (-) | Unclassified | Pseudomonas aeruginosa FFUP PS CB5 |
| ICE6440 | Unknown (-) | Unclassified | Pseudomonas aeruginosa HSV3483 |
| AICEScatt35120 | Unknown (-) | Unclassified | Streptomyces cattleya DSM |
| ICEValHN437 | Unknown (-) | SXT/R391 | Vibrio alginolyticus HN437 |
| ICEValHN396 | Unknown (-) | SXT/R391 | Vibrio alginolyticus HN396 |
| ICEValA056–2 | Unknown (-) | SXT/R391 | Vibrio alginolyticus 103826 |
| ICEKpnHS11286–2 | Unknown (-) | Unclassified | Klebsiella pneumoniae HS11286 |
| Tn6098 | Unknown (-) | Unclassified | Lactococcus lactis subsp. lactis KF147 |
| YAPI | Pathogenicity | ICEYe1 | Yersinia pseudotuberculosis IP 31758 |
| YAPI | Pathogenicity | ICEYe2 | Yersinia pseudotuberculosis 32777 |
| ICESenTy2–1 | Unknown (-) | SPI-7 | Salmonella enterica subsp. enterica serovar Typhi str. Ty2 |
| PFGI-1 | Unknown (-) | PAPI-1 | Pseudomonas fluorescens Pf-5 |
| ICEDzeEch1591–1 | Unknown (-) | ICEKp1 | Dickeya zeae Ech1591 |
| ICESb2 | Unknown (-) | SPI-7 | Salmonella bongori 2022/77 |
| ICESb1 | Unknown (-) | SPI-7 | Salmonella bongori CEIM46082 |
| ICEYe1 | Unknown (-) | ICEYe1 | Yersinia enterocolitica Y69 |
| ICEEc2 | Unknown (-) | Unclassified | Escherichia coli BEN374 |
| ICEPaeLESB58–1 | Mercuric resistance | ICEclc | Pseudomonas aeruginosa LESB58 |
| ICEVflInd1 | Antibiotic resistance; toxin-antitoxin system | SXT/R391 | Vibrio fluvialis Ind1 |
Integrative conjugative elements (ICE) were identified by sequence homology via BLAST against the ICEberg database.
Supplementary Materials
The following are available online at https://www.mdpi.com/2076-2607/8/3/448/s1, Figure S1: Occurrence (%) of genes that encode resistance for antibiotics in whole genome sequences of presumptive ESBL E. coli per source of isolates. Figure S2: Prevalence of genetic determinants of antimicrobial resistance in human- versus cattle-sourced E. coli isolates. Figure S3: (A) Pairwise SNP differences between clinical human isolates (n = 25) and municipal sewage isolates (n = 30), processing plants (n = 4), surrounding streams (n = 21), catch basins (n = 42), and cattle feces (n = 40). (B) Pairwise SNP differences between municipal sewage isolates (n = 30) and clinical human isolates (n = 25), processing plants (n = 4), surrounding streams (n = 21), catch basins (n = 42), and cattle feces (n = 40). Figure S4: Whole genome MLST of presumptive ESBL-producing E. coli genomes along the One Health continuum from human clinical, municipal sewage, beef processing, cattle feces, catch basins, and surface streams generated using 9580 wgMLST loci, colored based on isolate origin. Figure S5: Phylogenetic tree generated based on single-nucleotide polymorphisms (SNPs) of the core genes of 108 E. coli isolates that proved to be true ESBL producers (blaTEM, blaCTX-M, blaOXA, and blaSHV) obtained from cattle feces, catch basin, surface streams, municipal sewage, and human clinical isolates from the One Health continuum, with reference genome E. coli str. K-12 substr. MG1655 (GenBank accession # GI: 545778205/U00096.3). Figure S6: Relative abundance of 38 plasmid types in ESBL-producing E. coli from the One Health continuum of the beef-production system. Table S1: Antimicrobial resistance phenotypes of presumptive ESBL-E. coli. Table S2: Resistance determinants identified in presumptive ESBL-E. coli. TableS3: Whole genome sequencing assembly data for 162 presumptive ESBL-E. coli isolates. Table S4: Identification of mutations in the quinolone resistance-determining regions of presumptive ESBL-E. coli. Table S5: Least square means comparing antimicrobial resistance gene prevalence in cattle versus human isolates. Table S6: Genotype and phenotype comparison of presumptive ESBL E. coli isolates from multiple sources of the One Health continuum. Table S7: Serogroup and ST type of in presumptive ESBL-E. coli. Table S8: Differences between plasmid least square means estimates of the effect of isolate from cattle versus human sources. Table S9: Plasmids identified in presumptive ESBL-E. coli. Table S10: Integrative conjugative elements identified in presumptive ESBL-E. coli. Table S11: Virulence determinants identified in presumptive ESBL-E. coli. Table S12: Differences between virulence determinants identified in presumptive ESBL-E. coli isolated from different sources across a One Health continuum.
Author Contributions
Study conceptualization: E.H.A., C.A.N., and T.A.M.; methodology: E.H.A., T.A.M., and C.N.-B., R.Z., and R.R.; investigation: E.H.A., N.G., C.B., S.R.C., S.J.H., L.T., and D.C.; formal analysis: E.H.A. and M.W.; resources: T.A.M., C.N.-B., and K.A.; Supervision: T.A.M., C.N.-B.; writing-original draft: E.H.A.; writing-review and editing: E.H.A., M.W., C.N.-B., R.Z., S.R.C., N.G., L.T., S.J.H., D.C., C.W.B., K.A., C.A.N., R.R., and T.A.M. All authors have read and agreed to the published version of the manuscript.
Funding
This study was supported by funding from the Beef Cattle Research Council (BCRC Project FOS 10.13) and Genomic Research and Development Initiative (GRDI)-AMR program of the Federal Government of Canada. EHA is supported by the University of Manitoba PhD Graduate Fellowship. The funding agencies had no role in study design, sample collection, analysis, interpretation, and writing of the manuscript.
Acknowledgments
We thank Ruth Barbieri, Wendi Smart, and Jillian Rumore for technical assistance. We also gratefully acknowledge Francis Zvomuya, Dept. of Soil Science, University of Manitoba for his assistance with data analysis. We wish to also thank the cooperating commercial feedlots for allowing us access to collect samples.
Conflicts of Interest
CWB is part owner and managing partner of Feedlot Health Management Services and Southern Alberta Veterinary Services. SJH is an employee at Feedlot Health Management Services, Okotoks, Alberta, Canada. Feedlot Health is a private company that provides expert consultation regarding management and production of calf grower calves and feedlot cattle, including developing veterinary protocols to support animal health. Feedlot Health also conducts in-house and contract research related to dairy calf grower and feedlot production.
Abbreviations
| AMR | Antimicrobial resistance |
| ARG | Antimicrobial resistance gene |
| BProcessing | Beef processing plant |
| CBasins | Catch basins |
| CFeces | Cattle feces |
| CHumans | Clinically ill humans |
| ECP | E. coli common pilus |
| ESBL | Extended spectrum beta-lactamase |
| ESBL-EC | Extended spectrum beta-lactamase E. coli |
| ExPEC | Extraintestinal Pathogenic Escherichia coli |
| HGT | Horizontal gene transfer |
| ICE | Integrative Conjugative Element |
| LEE | Locus for enterocyte effacement |
| MDR | Multidrug resistance |
| MGE | Mobile genetic element |
| MSewage | Municipal sewage treatment |
| QRDR | Quinolone resistance-determining region |
| SNPs | single nucleotide polymorphisms |
| SStreams | Surrounding streams |
| STEC | Shiga toxigenic E. coli |
| WGS | Whole genome sequence |
1. Nicolas-Chanoine, M.H.; Bertrand, X.; Madec, J.Y. Escherichia coli ST131, an intriguing clonal group. Clin. Microbiol. Rev. 2014, 27, 543-574.
2. Johnson, J.R.; Nicolas-Chanoine, M.H.; DebRoy, C.; Castanheira, M.; Robicsek, A.; Hansen, G.; Weissman, S.; Urban, C.; Platell, J.; Trott, D.; et al. Comparison of Escherichia coli ST131 pulsotypes, by epidemiologic traits, 1967-2009. Emerg. Infect. Dis. 2012, 18, 598-607.
3. Ojer-Usoz, E.; González, D.; Vitas, A.I. Clonal diversity of ESBL-Producing Escherichia coli isolated from environmental, human and food samples. Int. J. Environ. Res. Public Health 2017, 14.
4. Rehman, M.A.; Yin, X.; Lepp, D.; Laing, C.; Ziebell, K.; Talbot, G.; Topp, E.; Diarra, M.S. Genomic analysis of third generation cephalosporin resistant Escherichia coli from dairy cow manure. Vet. Sci. 2017, 4.
5. World Health Organization. Global Priority List of Antibiotic-Resistant Bacteria to Guide Research, Discovery, and Development of New Antibiotics, WHO Essential Medicines and Health Products Guidelines. 2017. Available online: https://www.who.int/medicines/publications/global-priority-list-antibiotic-resistant-bacteria/en/ (accessed on 27 February 2017).
6. Hilty, M.; Betsch, B.Y.; Bögli-Stuber, K.; Heiniger, N.; Stadler, M.; Küffer, M.; Kronenberg, A.; Rohrer, C.; Aebi, S.; Endimiani, A.; et al. Transmission dynamics of extended-spectrum β-lactamase-producing Enterobacteriaceae in the tertiary care hospital and the household setting. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 2012, 55, 967-975.
7. Ewers, C.; Bethe, A.; Stamm, I.; Grobbel, M.; Kopp, P.A.; Guerra, B.; Stubbe, M.; Doi, Y.; Zong, Z.; Kola, A.; et al. CTX-M-15-D-ST648 Escherichia coli from companion animals and horses: Another pandemic clone combining multiresistance and extraintestinal virulence? J. Antimicrob. Chemother. 2014, 69, 1224-1230.
8. Liu, M.; Li, M.; Wu, L.; Song, Q.; Zhao, D.; Chen, Z.; Kang, M.; Xie, Y. Extended-spectrum β-lactamase-producing E. coli septicemia among rectal carriers in the ICU. Medicine 2018, 97, e12445.
9. Gudiol, C.; Royo-Cebrecos, C.; Tebe, C.; Abdala, E.; Akova, M.; Álvarez, R.; Maestro-de la Calle, G.; Cano, A.; Cervera, C.; Clemente, W.T.; et al. Clinical efficacy of β-lactam/β-lactamase inhibitor combinations for the treatment of bloodstream infection due to extended-spectrum β-lactamase-producing Enterobacteriaceae in haematological patients with neutropaenia: A study protocol for a retrospective observational study (BICAR). BMJ Open 2017, 7, e013268.
10. Malande, O.O.; Nuttall, J.; Pillay, V.; Bamford, C.; Eley, B. A ten-year review of ESBL and non-ESBL Escherichia coli bloodstream infections among children at a tertiary referral hospital in South Africa. PLoS ONE 2019, 14, e0222675.
11. Repessé, X.; Artiguenave, M.; Paktoris-Papine, S.; Espinasse, F.; Dinh, A.; Charron, C.; El Sayed, F.; Geri, G.; Vieillard-Baron, A. Epidemiology of extended-spectrum beta-lactamase-producing Enterobacteriaceae in an intensive care unit with no single rooms. Ann. Intensive Care 2017, 7, 73.
12. Roche, D.; Fléchard, M.; Lallier, N.; Répérant, M.; Brée, A.; Pascal, G.; Schouler, C.; Germon, P. ICEEc2, a new integrative and conjugative element belonging to the pKLC102/PAGI-2 family, identified in Escherichia coli strain BEN374. J. Bacteriol. 2010, 192, 5026-5036.
13. Partridge, S.R.; Kwong, S.M.; Firth, N.; Jensen, S.O. Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev. 2018, 31.
14. Carattoli, A. Resistance plasmid families in Enterobacteriaceae. Antimicrob. Agents Chemother. 2009, 53, 2227-2238.
15. Liu, M.; Li, X.; Xie, Y.; Bi, D.; Sun, J.; Li, J.; Tai, C.; Deng, Z.; Ou, H.Y. ICEberg 2.0: An updated database of bacterial integrative and conjugative elements. Nucleic Acids Res. 2019, 47, D660-D665.
16. Bi, D.; Xu, Z.; Harrison, E.M.; Tai, C.; Wei, Y.; He, X.; Jia, S.; Deng, Z.; Rajakumar, K.; Ou, H.Y. ICEberg: A web-based resource for integrative and conjugative elements found in Bacteria. Nucleic Acids Res. 2012, 40, D621-D626.
17. Schmidt, J.W.; Griffin, D.; Kuehn, L.A.; Brichta-Harhay, D.M. Influence of therapeutic ceftiofur treatments of feedlot cattle on fecal and hide prevalences of commensal Escherichia coli resistant to expanded-spectrum cephalosporins, and molecular characterization of resistant isolates. Appl. Environ. Microbiol. 2013, 79, 2273-2283.
18. Aslam, M.; Diarra, M.S.; Service, C.; Rempel, H. Antimicrobial resistance genes in Escherichia coli isolates recovered from a commercial beef processing plant. J. Food Prot. 2009, 72, 1089-1093.
19. Kojima, A.; Ishii, Y.; Ishihara, K.; Esaki, H.; Asai, T.; Oda, C.; Tamura, Y.; Takahashi, T.; Yamaguchi, K. Extended-spectrum-β-lactamase-producing Escherichia coli strains isolated from farm animals from 1999 to 2002: Report from the Japanese Veterinary Antimicrobial Resistance Monitoring Program. Antimicrob. Agents Chemother. 2005, 49, 3533-3537.
20. Schmidt, J.W.; Agga, G.E.; Bosilevac, J.M.; Brichta-Harhay, D.M.; Shackelford, S.D.; Wang, R.; Wheeler, T.L.; Arthur, T.M. Occurrence of antimicrobial-resistant Escherichia coli and Salmonella enterica in the beef cattle production and processing continuum. Appl. Environ. Microbiol. 2015, 81, 713-725.
21. Zaheer, R.; Cook, S.R.; Klima, C.L.; Stanford, K.; Alexander, T.; Topp, E.; Read, R.; McAllister, T.A. Effect of Subtherapeutic vs Therapeutic Administration of Macrolides on Antimicrobial Resistance in Mannheimia haemolytica and Enterococci isolated from Beef Cattle. Front. Microbiol. 2013, 4, 133.
22. van Belkum, A.; Dunne, W.M. Next-generation antimicrobial susceptibility testing. J. Clin. Microbiol. 2013, 51, 2018-2024.
23. Quainoo, S.; Coolen, J.P.M.; van Hijum, S.A.F.T.; Huynen, M.A.; Melchers, W.J.G.; van Schaik, W.; Wertheim, H.F.L. Whole-Genome Sequencing of bacterial pathogens: The future of nosocomial outbreak analysis. Clin. Microbiol. Rev. 2017, 30, 1015-1063.
24. McDermott, P.F.; Tyson, G.H.; Kabera, C.; Chen, Y.; Li, C.; Folster, J.P.; Ayers, S.L.; Lam, C.; Tate, H.P.; Zhao, S. Whole-genome sequencing for detecting antimicrobial resistance in Nontyphoidal Salmonella. Antimicrob. Agents Chemother. 2016, 60, 5515-5520.
25. Ludden, C.; Raven, K.E.; Jamrozy, D.; Gouliouris, T.; Blane, B.; Coll, F.; de Goffau, M.; Naydenova, P.; Horner, C.; Hernandez-Garcia, J.; et al. One Health genomic surveillance of Escherichia coli demonstrates distinct lineages and mobile genetic elements in isolates from humans versus livestock. mBio 2019, 10.
26. 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.
27. Rusconi, B.; Sanjar, F.; Koenig, S.S.; Mammel, M.K.; Tarr, P.I.; Eppinger, M. Whole genome sequencing for genomics-guided investigations of Escherichia coli O157:H7 Outbreaks. Front. Microbiol. 2016, 7, 985.
28. Vincent, C.; Usongo, V.; Berry, C.; Tremblay, D.M.; Moineau, S.; Yousfi, K.; Doualla-Bell, F.; Fournier, E.; Nadon, C.; Goodridge, L.; et al. Comparison of advanced whole genome sequence-based methods to distinguish strains of Salmonella enterica serovar Heidelberg involved in foodborne outbreaks in Québec. Food Microbiol. 2018, 73, 99-110.
29. WHO. Integrated Surveillance of Antimicrobial Resistance in Foodborne Bacteria: Application of a One Health Approach: Guidance from the WHO Advisory Group on Integrated Surveillanec of Antimicrobial Resistance (AGISAR). 2017. Available online: https://www.who.int/foodsafety/publications/agisar_guidance2017/en/.
30. Europrean Commission. A European One Health Action Plan Against Antimicrobial Resistance (AMR). 2017. Available online: https://ec.europa.eu/health/amr/antimicrobial-resistance_en/.
31. Rousham, E.K.; Unicomb, L.; Islam, M.A. Human, animal and environmental contributors to antibiotic resistance in low-resource settings: Integrating behavioural, epidemiological and One Health approaches. Proc. Biol. Sci. 2018, 285.
32. Beukers, A.G.; Zaheer, R.; Cook, S.R.; Chaves, A.V.; Ward, M.P.; Tymensen, L.; Morley, P.S.; Hannon, S.; Booker, C.W.; Read, R.R.; et al. Comparison of antimicrobial resistance genes in feedlots and urban wastewater. Can. J. Vet. Res. 2018, 82, 24-38.
33. Tymensen, L.; Booker, C.W.; Hannon, S.J.; Cook, S.R.; Jokinen, C.C.; Zaheer, R.; Read, R.; Boerlin, P.; McAllister, T.A. Plasmid distribution among Escherichia coli from livestock and associated wastewater: Unraveling factors that shape the presence of genes conferring Third-Generation Cephalosporin resistance. Environ. Sci. Technol. 2019, 53, 11666-11674.
34. Zaheer, R.; Cook, S.R.; Barbieri, R.; Goji, N.; Cameron, A.; Petkau, A.; Polo, R.O.; Tymensen, L.; Stamm, C.; Song, J.; et al. Surveillance of Enterococcus spp. reveals distinct species and antimicrobial resistance diversity across a One-Health continuum. Sci. Rep. 2020, 10, 3937.
35. Gill, C.O.; Jones, T. Microbiological sampling of carcasses by excision or swabbing. J. Food Prot. 2000, 63, 167-173.
36. CLSI. Clinical and Laboratory Standards Istitute. Performance Standards for Antimicrobial Susceptibility Testing; Twenty-Fifth Informational Supplement; CLSI document M100S-26; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2016; Volume 26.
37. Desjardins, P.; Conklin, D. NanoDrop microvolume quantitation of nucleic acids. J. Vis. Exp. 2010, 2565.
38. Gallagher, S.R. Quantitation of DNA and RNA with absorption and fluorescence Spectroscopy. Curr. Protoc. Immunol. 2017, 116, A.3L.1-A.3L.14.
39. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114-2120.
40. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455-477.
41. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068-2069.
42. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Cech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018, 46, W537-W544.
43. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691-3693.
44. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640-2644.
45. Babiker, A.; Mustapha, M.M.; Pacey, M.P.; Shutt, K.A.; Ezeonwuka, C.D.; Ohm, S.L.; Cooper, V.S.; Marsh, J.W.; Doi, Y.; Harrison, L.H. Use of online tools for antimicrobial resistance prediction by whole-genome sequencing in methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococci (VRE). J. Glob. Antimicrob. Resist. 2019, 19, 136-143.
46. SAS-Institute. User's Guide Version 5.0 SAS Institute Cary; SAS/STAT: Cary, NC, USA, 2015.
47. Petkau, A.; Mabon, P.; Sieffert, C.; Knox, N.C.; Cabral, J.; Iskander, M.; Weedmark, K.; Zaheer, R.; Katz, L.S.; Nadon, C.; et al. SNVPhyl: A single nucleotide variant phylogenomics pipeline for microbial genomic epidemiology. Microb. Genom. 2017, 3, e000116.
48. Blattner, F.R.; Plunkett, G.; Bloch, C.A.; Perna, N.T.; Burland, V.; Riley, M.; Collado-Vides, J.; Glasner, J.D.; Rode, C.K.; Mayhew, G.F.; et al. The complete genome sequence of Escherichia coli K-12. Science 1997, 277, 1453-1462.
49. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242-W245.
50. Jolley, K.A.; Bray, J.E.; Maiden, M.C.J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018, 3, 124.
51. Chen, L.; Yang, J.; Yu, J.; Yao, Z.; Sun, L.; Shen, Y.; Jin, Q. VFDB: A reference database for bacterial virulence factors. Nucleic Acids Res. 2005, 33, D325-D328.
52. Browne, A.S.; Biggs, P.J.; Wilkinson, D.A.; Cookson, A.L.; Midwinter, A.C.; Bloomfield, S.J.; Hranac, C.R.; Rogers, L.E.; Marshall, J.C.; Benschop, J.; et al. Use of genomics to investigate historical importation of Shiga toxin-producing Escherichia coli serogroup O26 and nontoxigenic variants into New Zealand. Emerg. Infect. Dis. 2019, 25, 489-500.
53. Ingle, D.J.; Valcanis, M.; Kuzevski, A.; Tauschek, M.; Inouye, M.; Stinear, T.; Levine, M.M.; Robins-Browne, R.M.; Holt, K.E. In silico serotyping of E. coli from short read data identifies limited novel O-loci but extensive diversity of O: H serotype combinations within and between pathogenic lineages. Microb. Genom. 2016, 2, e000064.
54. Petkau, A.; Stuart-Edwards, M.; Stothard, P.; Van Domselaar, G. Interactive microbial genome visualization with GView. Bioinformatics 2010, 26, 3125-3126.
55. Mbelle, N.M.; Feldman, C.; Osei Sekyere, J.; Maningi, N.E.; Modipane, L.; Essack, S.Y. The resistome, mobilome, virulome and phylogenomics of multidrug-resistant Escherichia coli clinical isolates from Pretoria, South Africa. Sci. Rep. 2019, 9, 16457.
56. Hopkins, K.L.; Davies, R.H.; Threlfall, E.J. Mechanisms of quinolone resistance in Escherichia coli and Salmonella: Recent developments. Int. J. Antimicrob. Agents 2005, 25, 358-373.
57. Cameron, A.; McAllister, T.A. Antimicrobial usage and resistance in beef production. J. Anim. Sci. Biotechnol. 2016, 7, 68.
58. Bryan, A.; Shapir, N.; Sadowsky, M.J. Frequency and distribution of tetracycline resistance genes in genetically diverse, nonselected, and nonclinical Escherichia coli strains isolated from diverse human and animal sources. Appl. Environ. Microbiol. 2004, 70, 2503-2507.
59. Wilkerson, C.; Samadpour, M.; van Kirk, N.; Roberts, M.C. Antibiotic resistance and distribution of tetracycline resistance genes in Escherichia coli O157:H7 isolates from humans and bovines. Antimicrob. Agents Chemother. 2004, 48, 1066-1067.
60. Xu, Y.; Sun, H.; Bai, X.; Fu, S.; Fan, R.; Xiong, Y. Occurrence of multidrug-resistant and ESBL-producing atypical enteropathogenic Escherichia coli in China. Gut Pathog. 2018, 10, 8.
61. Noyes, N.R.; Yang, X.; Linke, L.M.; Magnuson, R.J.; Cook, S.R.; Zaheer, R.; Yang, H.; Woerner, D.R.; Geornaras, I.; McArt, J.A.; et al. Characterization of the resistome in manure, soil and wastewater from dairy and beef production systems. Sci. Rep. 2016, 6, 24645.
62. Poyart, C.; Celli, J.; Trieu-Cuot, P. Conjugative transposition of Tn916-related elements from Enterococcus faecalis to Escherichia coli and Pseudomonas fluorescens. Antimicrob. Agents Chemother. 1995, 39, 500-506.
63. Gozi, K.S.; Froes, J.R.; Deus Ajude, L.P.T.; da Silva, C.R.; Baptista, R.S.; Peiró, J.R.; Marinho, M.; Mendes, L.C.N.; Nogueira, M.C.L.; Casella, T. Dissemination of multidrug-resistant commensal. Front. Microbiol. 2019, 10, 1394.
64. Seenama, C.; Thamlikitkul, V.; Ratthawongjirakul, P. Multilocus sequence typing and blaESBL characterization of extended-spectrum beta-lactamase-producing Escherichia coli isolated from healthy humans and swine in Northern Thailand. Infect. Drug Resist. 2019, 12, 2201-2214.
65. Frye, J.G.; Jackson, C.R. Genetic mechanisms of antimicrobial resistance identified in Salmonella enterica, Escherichia coli, and Enteroccocus spp. isolated from U.S. food animals. Front. Microbiol. 2013, 4, 135.
66. Wright, G.D. Antibiotic resistance in the environment: A link to the clinic? Curr. Opin. Microbiol. 2010, 13, 589-594.
67. Reinthaler, F.F.; Galler, H.; Feierl, G.; Haas, D.; Leitner, E.; Mascher, F.; Melkes, A.; Posch, J.; Pertschy, B.; Winter, I.; et al. Resistance patterns of Escherichia coli isolated from sewage sludge in comparison with those isolated from human patients in 2000 and 2009. J. Water Health 2013, 11, 13-20.
68. Pärnänen, K.M.M.; Narciso-da-Rocha, C.; Kneis, D.; Berendonk, T.U.; Cacace, D.; Do, T.T.; Elpers, C.; Fatta-Kassinos, D.; Henriques, I.; Jaeger, T.; et al. Antibiotic resistance in European wastewater treatment plants mirrors the pattern of clinical antibiotic resistance prevalence. Sci. Adv. 2019, 5, eaau9124.
69. Dortet, L.; Nordmann, P.; Poirel, L. Association of the emerging carbapenemase NDM-1 with a bleomycin resistance protein in Enterobacteriaceae and Acinetobacter baumannii. Antimicrob. Agents Chemother. 2012, 56, 1693-1697.
70. Slipski, C.J.; Zhanel, G.G.; Bay, D.C. Biocide selective TolC-independent efflux pumps in Enterobacteriaceae. J. Membr. Biol. 2018, 251, 15-33.
71. Gaze, W.H.; Abdouslam, N.; Hawkey, P.M.; Wellington, E.M. Incidence of class 1 integrons in a quaternary ammonium compound-polluted environment. Antimicrob. Agents Chemother. 2005, 49, 1802-1807.
72. Jaglic, Z.; Cervinkova, D. Genetic basis of resistance to quaternary ammonium compounds--the qac genes and their role: A review. Vet. Med. 2012, 57.
73. Health Quality Ontario. Portable ultraviolet light surface-disinfecting devices for prevention of hospital-acquired infections: A health technology assessment. Ont. Health Technol. Assess. Ser. 2018, 18, 1-73.
74. Han, J.H.; Sullivan, N.; Leas, B.F.; Pegues, D.A.; Kaczmarek, J.L.; Umscheid, C.A. Cleaning hospital room surfaces to prevent health care-associated infections: A technical brief. Ann. Intern. Med. 2015, 163, 598-607.
75. Cormier, A.; Zhang, P.L.C.; Chalmers, G.; Weese, J.S.; Deckert, A.; Mulvey, M.; McAllister, T.; Boerlin, P. Diversity of CTX-M-positive Escherichia coli recovered from animals in Canada. Vet. Microbiol. 2019, 231, 71-75.
76. Pitout, J.D.; DeVinney, R. ST131: A multidrug-resistant clone primed for global domination. F1000Res 2017, 6.
77. Grönthal, T.; Österblad, M.; Eklund, M.; Jalava, J.; Nykäsenoja, S.; Pekkanen, K.; Rantala, M. Sharing more than friendship - transmission of NDM-5 ST167 and CTX-M-9 ST69 Escherichia coli between dogs and humans in a family, Finland, 2015. Euro Surveill. 2018, 23.
78. 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.
79. Pouget, J.G.; Coutinho, F.J.; Reid-Smith, R.J.; Boerlin, P. Characterization of bla(SHV) genes on plasmids from Escherichia coli and Salmonella enterica isolates from Canadian food animals (2006-2007). Appl. Environ. Microbiol. 2013, 79, 3864-3866.
80. Carattoli, A.; Lovari, S.; Franco, A.; Cordaro, G.; Di Matteo, P.; Battisti, A. Extended-spectrum beta-lactamases in Escherichia coli isolated from dogs and cats in Rome, Italy, from 2001 to 2003. Antimicrob. Agents Chemother. 2005, 49, 833-835.
81. Teshager, T.; Domínguez, L.; Moreno, M.A.; Saénz, Y.; Torres, C.; Cardeñosa, S. Isolation of an SHV-12 beta-lactamase-producing Escherichia coli strain from a dog with recurrent urinary tract infections. Antimicrob. Agents Chemother. 2000, 44, 3483-3484.
82. Huang, Z.; Mi, Z.; Wang, C. A novel beta-lactamase gene, LAP-2, produced by an Enterobacter cloacae clinical isolate in China. J. Hosp. Infect. 2008, 70, 95-96.
83. Le, V.; Nhu, N.T.; Cerdeno-Tarraga, A.; Campbell, J.I.; Tuyen, H.T.; Nhu, T.o.H.; Tam, P.T.; Schultsz, C.; Thwaites, G.; Thomson, N.R.; et al. Genetic characterization of three qnrS1-harbouring multidrug-resistance plasmids and qnrS1-containing transposons circulating in Ho Chi Minh City, Vietnam. J. Med. Microbiol. 2015, 64, 869-878.
84. Park, Y.J.; Yu, J.K.; Kim, S.I.; Lee, K.; Arakawa, Y. Accumulation of plasmid-mediated fluoroquinolone resistance genes, qepA and qnrS1, in Enterobacter aerogenes co-producing RmtB and class A beta-lactamase LAP-1. Ann. Clin. Lab. Sci. 2009, 39, 55-59.
85. Varughese, L.R.; Rajpoot, M.; Goyal, S.; Mehra, R.; Chhokar, V.; Beniwal, V. Analytical profiling of mutations in quinolone resistance determining region of gyrA gene among UPEC. PLoS ONE 2018, 13, e0190729.
86. Stoesser, N.; Batty, E.M.; Eyre, D.W.; Morgan, M.; Wyllie, D.H.; Del Ojo Elias, C.; Johnson, J.R.; Walker, A.S.; Peto, T.E.; Crook, D.W. Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data. J. Antimicrob. Chemother. 2013, 68, 2234-2244.
87. Tyson, G.H.; McDermott, P.F.; Li, C.; Chen, Y.; Tadesse, D.A.; Mukherjee, S.; Bodeis-Jones, S.; Kabera, C.; Gaines, S.A.; Loneragan, G.H.; et al. WGS accurately predicts antimicrobial resistance in Escherichia coli. J. Antimicrob. Chemother. 2015, 70, 2763-2769.
88. Palmer, A.C.; Chait, R.; Kishony, R. Nonoptimal gene expression creates latent potential for antibiotic resistance. Mol. Biol. Evol. 2018, 35, 2669-2684.
89. Gouliouris, T.; Raven, K.E.; Ludden, C.; Blane, B.; Corander, J.; Horner, C.S.; Hernandez-Garcia, J.; Wood, P.; Hadjirin, N.F.; Radakovic, M.; et al. Genomic Surveillance of Enterococcus faecium Reveals Limited Sharing of Strains and Resistance Genes between Livestock and Humans in the United Kingdom. mBio 2018, 9.
90. Manges, A.R. Escherichia coli causing bloodstream and other extraintestinal infections: Tracking the next pandemic. Lancet Infect. Dis. 2019, 19, 1269-1270.
91. Salinas, L.; Cárdenas, P.; Johnson, T.J.; Vasco, K.; Graham, J.; Trueba, G. Diverse commensal Escherichia coli clones and plasmids disseminate antimicrobial resistance genes in domestic animals and children in a Semirural Community in Ecuador. mSphere 2019, 4.
92. Mulvey, M.R.; Susky, E.; McCracken, M.; Morck, D.W.; Read, R.R. Similar cefoxitin-resistance plasmids circulating in Escherichia coli from human and animal sources. Vet. Microbiol. 2009, 134, 279-287.
93. Peirano, G.; Richardson, D.; Nigrin, J.; McGeer, A.; Loo, V.; Toye, B.; Alfa, M.; Pienaar, C.; Kibsey, P.; Pitout, J.D. High prevalence of ST131 isolates producing CTX-M-15 and CTX-M-14 among extended-spectrum-beta-lactamase-producing Escherichia coli isolates from Canada. Antimicrob. Agents Chemother. 2010, 54, 1327-1330.
94. Pietsch, M.; Irrgang, A.; Roschanski, N.; Brenner Michael, G.; Hamprecht, A.; Rieber, H.; Käsbohrer, A.; Schwarz, S.; Rösler, U.; Kreienbrock, L.; et al. Whole genome analyses of CMY-2-producing Escherichia coli isolates from humans, animals and food in Germany. BMC Genom. 2018, 19, 601.
95. Mathers, A.J.; Peirano, G.; Pitout, J.D. The role of epidemic resistance plasmids and international high-risk clones in the spread of multidrug-resistant Enterobacteriaceae. Clin. Microbiol. Rev. 2015, 28, 565-591.
96. Adator, E.H.; Cheng, M.; Holley, R.; McAllister, T.; Narvaez-Bravo, C. Ability of Shiga toxigenic Escherichia coli to survive within dry-surface biofilms and transfer to fresh lettuce. Int. J. Food Microbiol. 2018, 269, 52-59.
97. Wang, J.; Stanford, K.; McAllister, T.A.; Johnson, R.P.; Chen, J.; Hou, H.; Zhang, G.; Niu, Y.D. Biofilm formation, virulence gene profiles, and antimicrobial resistance of nine serogroups of non-O157 shiga toxin-producing Escherichia coli. Foodborne Pathog. Dis. 2016, 13, 316-324.
98. Cernicchiaro, N.; Cull, C.A.; Paddock, Z.D.; Shi, X.; Bai, J.; Nagaraja, T.G.; Renter, D.G. Prevalence of Shiga toxin-producing Escherichia coli and associated virulence genes in feces of commercial feedlot cattle. Foodborne Pathog. Dis. 2013, 10, 835-841.
99. Stanford, K.; Johnson, R.P.; Alexander, T.W.; McAllister, T.A.; Reuter, T. Influence of season and feedlot location on prevalence and virulence factors of seven serogroups of Escherichia coli in feces of western-Canadian slaughter cattle. PLoS ONE 2016, 11, e0159866.
100. Aslantaş, Ö.; Yilmaz, E. Prevalence and molecular characterization of extended-spectrum β-lactamase (ESBL) and plasmidic AmpC β-lactamase (pAmpC) producing Escherichia coli in dogs. J. Vet. Med. Sci. 2017, 79, 1024-1030.
101. Sharp, H.; Valentin, L.; Fischer, J.; Guerra, B.; Appel, B.; Käsbohrer, A. Estimation of the transfer of ESBL-producing Escherichia coli to humans in Germany. Berl. Munch. Tierarztl. Wochenschr. 2014, 127, 464-477.
102. Johnson, T.J.; Wannemuehler, Y.M.; Johnson, S.J.; Logue, C.M.; White, D.G.; Doetkott, C.; Nolan, L.K. Plasmid replicon typing of commensal and pathogenic Escherichia coli isolates. Appl. Environ. Microbiol. 2007, 73, 1976-1983.
103. Rodriguez-Siek, K.E.; Giddings, C.W.; Doetkott, C.; Johnson, T.J.; Fakhr, M.K.; Nolan, L.K. Comparison of Escherichia coli isolates implicated in human urinary tract infection and avian colibacillosis. Microbiology 2005, 151, 2097-2110.
104. Obi, C.C.; Vayla, S.; de Gannes, V.; Berres, M.E.; Walker, J.; Pavelec, D.; Hyman, J.; Hickey, W.J. The Integrative Conjugative Element clc (ICEclc) of Pseudomonas aeruginosa JB2. Front. Microbiol. 2018, 9, 1532.
105. Johnson, C.M.; Grossman, A.D. Integrative and Conjugative Elements (ICEs): What They Do and How They Work. Annu. Rev. Genet. 2015, 49, 577-601.
106. Wozniak, R.A.; Fouts, D.E.; Spagnoletti, M.; Colombo, M.M.; Ceccarelli, D.; Garriss, G.; Déry, C.; Burrus, V.; Waldor, M.K. Comparative ICE genomics: Insights into the evolution of the SXT/R391 family of ICEs. PLoS Genet. 2009, 5, e1000786.
107. Garnett, J.A.; Martínez-Santos, V.I.; Saldaña, Z.; Pape, T.; Hawthorne, W.; Chan, J.; Simpson, P.J.; Cota, E.; Puente, J.L.; Girón, J.A.; et al. Structural insights into the biogenesis and biofilm formation by the Escherichia coli common pilus. Proc. Natl. Acad. Sci. USA 2012, 109, 3950-3955.
108. Cameron, E.A.; Curtis, M.M.; Kumar, A.; Dunny, G.M.; Sperandio, V. Microbiota and pathogen proteases modulate type III secretion activity in enterohemorrhagic Escherichia coli. mBio 2018, 9.
109. Bihannic, M.; Ghanbarpour, R.; Auvray, F.; Cavalié, L.; Châtre, P.; Boury, M.; Brugère, H.; Madec, J.Y.; Oswald, E. Identification and detection of three new F17 fimbrial variants in Escherichia coli strains isolated from cattle. Vet. Res. 2014, 45, 76.
110. Puhar, A.; Sansonetti, P.J. Type III secretion system. Curr. Biol. 2014, 24, R784-R791.
Emelia H. Adator1, Matthew Walker2, Claudia Narvaez-Bravo1, Rahat Zaheer3, Noriko Goji4, Shaun R. Cook5, Lisa Tymensen5, Sherry J. Hannon6, Deirdre Church6, Calvin W. Booker6, Kingsley Amoako4, Celine A. Nadon2, Ron Read7 and Tim A. McAllister1,3,*
1Department of Food and Human Nutritional Sciences, University of Manitoba, Winnipeg, MB R3T 2N2, Canada
2National Microbiology Laboratory, Winnipeg, MB R3E 3R2, Canada
3Lethbridge Research and Development Centre, Lethbridge, AB T1J 4B1, Canada
4Canadian Food Inspection Agency, National Center for Animal Disease, Lethbridge Laboratory, Lethbridge, AB T1J 3Z4, Canada
5Alberta Agriculture and Forestry, Lethbridge, AB T1J 4V6, Canada
6Feedlot Health Management Services Ltd., Okotoks, AB T1S 2A2, Canada
7Cumming School of Medicine, University of Calgary, Calgary, AB T2N 4N1 Canada
*Author to whom correspondence should be addressed.
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Goji, Noriko; Cook, Shaun R; Tymensen, Lisa; Hannon, Sherry J; Church, Deirdre; Booker, Calvin W
; Amoako, Kingsley; Nadon, Celine A; Read, Ron; McAllister, Tim A